Biosynthesis of riccionidins and marchantins is regulated by R2R3-MYB transcription factors in Marchantia polymorpha

The landscape of transcription factor promoter activity during vegetative development in Marchantia

Transcription factors (TFs) are essential for the regulation of gene expression and cell fate determination. Characterizing the transcriptional activity of TF genes in space and time is a critical step toward understanding complex biological systems. The vegetative gametophyte meristems of bryophytes share some characteristics with the shoot apical meristems of flowering plants. However, the identity and expression profiles of TFs associated with gametophyte organization are largely unknown. With only ∼450 putative TF genes, Marchantia (Marchantia polymorpha) is an outstanding model system for plant systems biology. We have generated a near-complete collection of promoter elements derived from Marchantia TF genes. We experimentally tested reporter fusions for all the TF promoters in the collection and systematically analyzed expression patterns in Marchantia gemmae. This allowed us to build a map of expression domains in early vegetative development and identify a set of TF-derived promoters that are active in the stem-cell zone. The cell markers provide additional tools and insight into the dynamic regulation of the gametophytic meristem and its evolution. In addition, we provide an online database of expression patterns for all promoters in the collection. We expect that these promoter elements will be useful for cell-type-specific expression, synthetic biology applications, and functional genomics.

DELLA proteins recruit the Mediator complex subunit MED15 to coactivate transcription in land plants

DELLA proteins are negative regulators of the gibberellin response pathway in angiosperms, acting as central hubs that interact with hundreds of transcription factors (TFs) and regulators to modulate their activities. While the mechanism of TF sequestration by DELLAs to prevent DNA binding to downstream targets has been extensively documented, the mechanism that allows them to act as coactivators remains to be understood. Here, we demonstrate that DELLAs directly recruit the Mediator complex to specific loci in Arabidopsis, facilitating transcription. This recruitment involves DELLA amino-terminal domain and the conserved MED15 KIX domain. Accordingly, partial loss of MED15 function mainly disrupted processes known to rely on DELLA coactivation capacity, including cytokinin-dependent regulation of meristem function and skotomorphogenic response, gibberellin metabolism feedback, and flavonol production. We have also found that the single DELLA protein in the liverwort Marchantia polymorpha is capable of recruiting MpMED15 subunits, contributing to transcriptional coactivation. The conservation of Mediator-dependent transcriptional coactivation by DELLA between Arabidopsis and Marchantia implies that this mechanism is intrinsic to the emergence of DELLA in the last common ancestor of land plants.

Unveiling the molecular mechanisms of arsenic tolerance and resilience in the primitive bryophyte Marchantia polymorpha L

This study addresses the pressing issue of high arsenic (As) contaminations, which poses a severe threat to various life forms in our ecosystem. Despite this prevailing concern, all organisms have developed some techniques to mitigate the toxic effects of As. Certain plants, such as bryophytes, the earliest land plants, exhibit remarkable tolerance to wide range of harsh environmental conditions, due to their inherent competence. In this study, bryophytes collected from West Bengal, India, across varying contamination levels were investigated for their As tolerance capabilities. Assessment of As accumulation potential and antioxidant defense efficiency, including SOD, CAT, APX, GPX etc. revealed Marchantia polymorpha as the most tolerant species. It exhibited highest As accumulation, antioxidative proficiency, and minimal damage. Transcriptomic analysis of M. polymorpha exposed to 40 μM As(III) for 24 and 48 h identified several early responsive differentially expressing genes (DEGs) associated with As tolerance. These includes GSTs, GRXs, Hsp20s, SULTR1;2, ABCC2 etc., indicating a mechanism involving vacuolar sequestration. Interestingly, one As(III) efflux-transporter ACR3, an extrusion pump, known to combat As toxicity was found to be differentially expressed compared to control. The SEM-EDX analysis, further elucidated the operation of As extrusion mechanism, which contributes added As resilience in M. polymorpha. Yeast complementation assay using Δacr3 yeast cells, showed increased tolerance towards As(III), compared to the mutant cells, indicating As tolerant phenotype. Overall, these findings significantly enhance our understanding of As tolerance mechanisms in bryophytes. This can pave the way for the development of genetically engineered plants with heightened As tolerance and the creation of improved plant varieties.

MpTGA, together with MpNPR, regulates sexual reproduction and independently affects oil body formation in Marchantia polymorpha

In angiosperms, basic leucine-zipper (bZIP) TGACG-motif-binding (TGA) transcription factors (TFs) regulate developmental and stress-related processes, the latter often involving NON EXPRESSOR OF PATHOGENESIS-RELATED GENES (NPR) coregulator interactions. To gain insight into their functions in an early diverging land-plant lineage, the single MpTGA and sole MpNPR genes were investigated in the liverwort Marchantia polymorpha. We generated Marchantia MpTGA and MpNPR knockout and overexpression mutants and conducted morphological, transcriptomic and expression studies. Furthermore, we investigated MpTGA interactions with wild-type and mutagenized MpNPR and expanded our analyses including TGA TFs from two streptophyte algae. Mptga mutants fail to induce the switch from vegetative to reproductive development and lack gametangiophore formation. MpTGA and MpNPR proteins interact and Mpnpr mutant analysis reveals a novel coregulatory NPR role in sexual reproduction. Additionally, MpTGA acts independently of MpNPR as a repressor of oil body (OB) formation and can thereby affect herbivory. The single MpTGA TF exerts a dual role in sexual reproduction and OB formation in Marchantia. Common activities of MpTGA/MpNPR in sexual development suggest that coregulatory interactions were established after emergence of land-plant-specific NPR genes and contributed to the diversification of TGA TF functions during land-plant evolution.

The endomembrane system: how does it contribute to plant secondary metabolism?

New organelle acquisition through neofunctionalization of the endomembrane system (ES) with respect to plant secondary metabolism is a key evolutionary strategy for plant adaptation, which is overlooked due to the complexity of angiosperms. Bryophytes produce a broad range of plant secondary metabolites (PSMs), and their simple cellular structures, including unique organelles, such as oil bodies (OBs), highlight them as suitable model to investigate the contribution of the ES to PSMs. In this opinion, we review latest findings on the contribution of the ES to PSM biosynthesis, with a specific focus on OBs, and propose that the ES provides organelles and trafficking routes for PSM biosynthesis, transportation, and storage. Therefore, future research on ES-derived organelles and trafficking routes will provide essential knowledge for synthetic applications.

The maturation and aging trajectory of Marchantia polymorpha at single-cell resolution

Bryophytes represent a sister to the rest of land plants. Despite their evolutionary importance and relatively simple body plan, a comprehensive understanding of the cell types and transcriptional states that underpin the temporal development of bryophytes has not been achieved. Using time-resolved single-cell RNA sequencing, we define the cellular taxonomy of Marchantia polymorpha across asexual reproduction phases. We identify two maturation and aging trajectories of the main plant body of M. polymorpha at single-cell resolution: the gradual maturation of tissues and organs along the tip-to-base axis of the midvein and the progressive decline of meristem activities in the tip along the chronological axis. Specifically, we observe that the latter aging axis is temporally correlated with the formation of clonal propagules, suggesting an ancient strategy to optimize allocation of resources to producing offspring. Our work thus provides insights into the cellular heterogeneity that underpins the temporal development and aging of bryophytes.

A Polyphenol Oxidase Catalyzes Aurone Synthesis in Marchantia polymorpha

Aurones constitute one of the major classes of flavonoids, with a characteristic furanone structure that acts as the C-ring of flavonoids. Members of various enzyme families are involved in aurone biosynthesis in different higher plants, suggesting that during evolution plants acquired the ability to biosynthesize aurones independently and convergently. Bryophytes also produce aurones, but the biosynthetic pathways and enzymes involved have not been determined. The present study describes the identification and characterization of a polyphenol oxidase (PPO) that acts as an aureusidin synthase (MpAS1) in the model liverwort, Marchantia polymorpha. Crude enzyme assays using an M. polymorpha line overexpressing MpMYB14 with high accumulation of aureusidin showed that aureusidin was biosynthesized from naringenin chalcone and converted to riccionidin A. This activity was inhibited by N-phenylthiourea, an inhibitor specific to enzymes of the PPO family. Of the six PPOs highly induced in the line overexpressing MpMyb14, one, MpAS1, was found to biosynthesize aureusidin from naringenin chalcone when expressed in Saccharomyces cerevisiae. MpAS1 also recognized eriodictyol chalcone, isoliquiritigenin and butein, showing the highest activity for eriodictyol chalcone. Members of the PPO family in M. polymorpha evolved independently from PPOs in higher plants, indicating that aureusidin synthases evolved in parallel in land plants.

Interaction of PKR with STCS1: an indispensable step in the biosynthesis of lunularic acid in Marchantia polymorpha

Unlike bibenzyls derived from the vascular plants, lunularic acid (LA), a key precursor for macrocyclic bisbibenzyl synthesis in nonvascular liverworts, exhibits the absence of one hydroxy group within the A ring. It was hypothesized that both polyketide reductase (PKR) and stilbenecarboxylate synthase 1 (STCS1) were involved in the LA biosynthesis, but the underlined mechanisms have not been clarified. This study used bioinformatics analysis with molecular, biochemical and physiological approaches to characterize STCS1s and PKRs involved in the biosynthesis of LA. The results indicated that MpSTCS1s from Marchantia polymorpha catalyzed both C2→C7 aldol-type and C6→C1 Claisen-type cyclization using dihydro-p-coumaroyl-coenzyme A (CoA) and malonyl-CoA as substrates to yield a C6-C2-C6 skeleton of dihydro-resveratrol following decarboxylation and the C6-C3-C6 type of phloretin in vitro. The protein-protein interaction of PKRs with STCS1 (PPI-PS) was revealed and proved essential for LA accumulation when transiently co-expressed in Nicotiana benthamiana. Moreover, replacement of the active domain of STCS1 with an 18-amino-acid fragment from the chalcone synthase led to the PPI-PS greatly decreasing and diminishing the formation of LA. The replacement also increased the chalcone formation in STCS1s. Our results highlight a previously unrecognized PPI in planta that is indispensable for the formation of LA.

The origin of a land flora

The origin of a land flora fundamentally shifted the course of evolution of life on earth, facilitating terrestrialization of other eukaryotic lineages and altering the planet's geology, from changing atmospheric and hydrological cycles to transforming continental erosion processes. Despite algal lineages inhabiting the terrestrial environment for a considerable preceding period, they failed to evolve complex multicellularity necessary to conquer the land. About 470 million years ago, one lineage of charophycean alga evolved complex multicellularity via developmental innovations in both haploid and diploid generations and became land plants (embryophytes), which rapidly diversified to dominate most terrestrial habitats. Genome sequences have provided unprecedented insights into the genetic and genomic bases for embryophyte origins, with some embryophyte-specific genes being associated with the evolution of key developmental or physiological attributes, such as meristems, rhizoids and the ability to form mycorrhizal associations. However, based on the fossil record, the evolution of the defining feature of embryophytes, the embryo, and consequently the sporangium that provided a reproductive advantage, may have been most critical in their rise to dominance. The long timeframe and singularity of a land flora were perhaps due to the stepwise assembly of a large constellation of genetic innovations required to conquer the terrestrial environment.

MarpolBase Expression: A Web-Based, Comprehensive Platform for Visualization and Analysis of Transcriptomes in the Liverwort Marchantia polymorpha

The liverwort Marchantia polymorpha is equipped with a wide range of molecular and genetic tools and resources that have led to its wide use to explore the evo-devo aspects of land plants. Although its diverse transcriptome data are rapidly accumulating, there is no extensive yet user-friendly tool to exploit such a compilation of data and to summarize results with the latest annotations. Here, we have developed a web-based suite of tools, MarpolBase Expression (MBEX, https://marchantia.info/mbex/), where users can visualize gene expression profiles, identify differentially expressed genes, perform co-expression and functional enrichment analyses and summarize their comprehensive output in various portable formats. Using oil body biogenesis as an example, we demonstrated that the results generated by MBEX were consistent with the published experimental evidence and also revealed a novel transcriptional network in this process. MBEX should facilitate the exploration and discovery of the genetic and functional networks behind various biological processes in M. polymorpha and promote our understanding of the evolution of land plants.

Evolution and function of red pigmentation in land plants

BACKGROUND

Land plants commonly produce red pigmentation as a response to environmental stressors, both abiotic and biotic. The type of pigment produced varies among different land plant lineages. In the majority of species they are flavonoids, a large branch of the phenylpropanoid pathway. Flavonoids that can confer red colours include 3-hydroxyanthocyanins, 3-deoxyanthocyanins, sphagnorubins and auronidins, which are the predominant red pigments in flowering plants, ferns, mosses and liverworts, respectively. However, some flowering plants have lost the capacity for anthocyanin biosynthesis and produce nitrogen-containing betalain pigments instead. Some terrestrial algal species also produce red pigmentation as an abiotic stress response, and these include both carotenoid and phenolic pigments.

SCOPE

In this review, we examine: which environmental triggers induce red pigmentation in non-reproductive tissues; theories on the functions of stress-induced pigmentation; the evolution of the biosynthetic pathways; and structure-function aspects of different pigment types. We also compare data on stress-induced pigmentation in land plants with those for terrestrial algae, and discuss possible explanations for the lack of red pigmentation in the hornwort lineage of land plants.

CONCLUSIONS

The evidence suggests that pigment biosynthetic pathways have evolved numerous times in land plants to provide compounds that have red colour to screen damaging photosynthetically active radiation but that also have secondary functions that provide specific benefits to the particular land plant lineage.

The renaissance and enlightenment of Marchantia as a model system

The liverwort Marchantia polymorpha has been utilized as a model for biological studies since the 18th century. In the past few decades, there has been a Renaissance in its utilization in genomic and genetic approaches to investigating physiological, developmental, and evolutionary aspects of land plant biology. The reasons for its adoption are similar to those of other genetic models, e.g. simple cultivation, ready access via its worldwide distribution, ease of crossing, facile genetics, and more recently, efficient transformation, genome editing, and genomic resources. The haploid gametophyte dominant life cycle of M. polymorpha is conducive to forward genetic approaches. The lack of ancient whole-genome duplications within liverworts facilitates reverse genetic approaches, and possibly related to this genomic stability, liverworts possess sex chromosomes that evolved in the ancestral liverwort. As a representative of one of the three bryophyte lineages, its phylogenetic position allows comparative approaches to provide insights into ancestral land plants. Given the karyotype and genome stability within liverworts, the resources developed for M. polymorpha have facilitated the development of related species as models for biological processes lacking in M. polymorpha.

Moss transcription factors regulating development and defense responses to stress

Transcription factors control gene expression, leading to regulation of biological processes that determine plant development and adaptation to the environment. Land colonization by plants occurred 450-470 million years ago and was accompanied by an increase in the complexity of transcriptional regulation associated to transcription factor gene expansions. AP2/ERF, bHLH, MYB, NAC, GRAS, and WRKY transcription factor families increased in land plants compared with algae. In angiosperms, they play crucial roles in regulating plant growth and responses to environmental stressors. However, less information is available in bryophytes and only in a few cases is the functional role of moss transcription factors in stress mechanisms known. In this review, we discuss current knowledge of the transcription factor families involved in development and defense responses to stress in mosses and other bryophytes. By exploring and analysing the Physcomitrium patens public database and published transcriptional profiles, we show that a high number of AP2/ERF, bHLH, MYB, NAC, GRAS, and WRKY genes are differentially expressed in response to abiotic stresses and during biotic interactions. Expression profiles together with a comprehensive analysis provide insights into relevant transcription factors involved in moss defenses, and hint at distinct and conserved biological roles between bryophytes and angiosperms.

Bryophyte ultraviolet-omics: from genes to the environment

Ultraviolet (UV) radiation has contributed to the evolution of organisms since the origins of life. Bryophytes also have evolutionary importance as the first clearly identified lineage of land plants (embryophytes) colonizing the terrestrial environment, thus facing high UV and water scarcity, among other new challenges. Here we review bryophyte UV-omics, the discipline relating bryophytes and UV, with an integrative perspective from genes to the environment. We consider species and habitats investigated, methodology, response variables, protection mechanisms, environmental interactions, UV biomonitoring, molecular and evolutionary aspects, and applications. Bryophyte UV-omics shows convergences and divergences with the UV-omics of other photosynthetic organisms, from algae to tracheophytes. All these organisms converge in that UV damage may be limited under realistic UV levels, due to structural protection and/or physiological acclimation capacity. Nevertheless, bryophytes diverge because they have a unique combination of vegetative and reproductive characteristics to cope with high UV and other concomitant adverse processes, such as desiccation. This interaction has both evolutionary and ecological implications. In addition, UV effects on bryophytes depend on the species and the evolutionary lineage considered, with mosses more UV-tolerant than liverworts. Thus, bryophytes do not constitute a homogeneous functional type with respect to their UV tolerance.

Stress, senescence, and specialized metabolites in bryophytes

Life on land exposes plants to varied abiotic and biotic environmental stresses. These environmental drivers contributed to a large expansion of metabolic capabilities during land plant evolution and species diversification. In this review we summarize knowledge on how the specialized metabolite pathways of bryophytes may contribute to stress tolerance capabilities. Bryophytes are the non-tracheophyte land plant group (comprising the hornworts, liverworts, and mosses) and rapidly diversified following the colonization of land. Mosses and liverworts have as wide a distribution as flowering plants with regard to available environments, able to grow in polar regions through to hot desert landscapes. Yet in contrast to flowering plants, for which the biosynthetic pathways, transcriptional regulation, and compound function of stress tolerance-related metabolite pathways have been extensively characterized, it is only recently that similar data have become available for bryophytes. The bryophyte data are compared with those available for angiosperms, including examining how the differing plant forms of bryophytes and angiosperms may influence specialized metabolite diversity and function. The involvement of stress-induced specialized metabolites in senescence and nutrient response pathways is also discussed.

Liverwort oil bodies: diversity, biochemistry, and molecular cell biology of the earliest secretory structure of land plants

Liverworts are known for their large chemical diversity. Much of this diversity is synthesized and enclosed within oil bodies (OBs), a synapomorphy of the lineage. OBs contain the enzymes to biosynthesize and store large quantities of sesquiterpenoids and other compounds while limiting their cytotoxicity. Recent important biochemical and molecular discoveries related to OB formation, diversity, and biochemistry allow comparison with other secretory structures of land plants from an evo-devo perspective. This review addresses and discusses the most recent advances in OB origin, development, and function towards understanding the importance of these organelles in liverwort physiology and adaptation to changing environments. Our mapping of OB types and chemical compounds to the current liverwort phylogeny suggests that OBs were present in the most recent common ancestor of liverworts, supporting that OBs evolved as the first secretory structures in land plants. Yet, we require better sampling to define the macroevolutionary pattern governing the ancestral type of OB. We conclude that current efforts to find molecular mechanisms responsible for the morphological and chemical diversity of secretory structures will help understand the evolution of each major group of land plants, and open new avenues in biochemical research on bioactive compounds in bryophytes and vascular plants.

Biomass and Active Compounds Accumulation of the Medicinal Orchid Pleione bulbocodioides in Response to Light Intensity and Irrigation Frequency

Pseudobulbs of Pleione species are widely used as traditional medicine in Asian countries, but the mechanism of active compound accumulation remains unclear. In the present study, we investigated the accumulation of biomass and three active compounds (dactylorhin A, militarine and batatasin III) of Pleione bulbocodioides in response to different light intensities and irrigation frequencies. We found that single high light (65 % of full sunlight) or drought stress (14-day irrigation interval) increased active compounds accumulation but the combined effect of these two treatments decreased the total content of these three active compounds. This decrease was due to the plants under combined stress having a significantly lower photosynthetic rate, leaf area and longevity, leading to a dramatic decrease in pseudobulb biomass. Among all treatments, the highest total content of active compounds was observed in plants subjected to the high light level with a high water level (3-day irrigation interval), and plants under medium light intensity (30 % of full sunlight) also had considerable content of active compounds accumulation. To balance the quality and quantity of Pleione pseudobulbs during artificial cultivation, 30∼65 % of full sunlight with the avoidance of drought stress is recommended. Our results suggest the accumulation of the three active compounds is significantly influenced by light intensity and irrigation frequency, which may contribute to the artificial cultivation and quality control of medicinal Pleione.

KANADI promotes thallus differentiation and FR-induced gametangiophore formation in the liverwort Marchantia

In angiosperms, KANADI transcription factors have roles in the sporophyte generation regulating tissue polarity, organogenesis and shade avoidance responses, but are not required during the gametophyte generation. Whether these roles are conserved in the gametophyte-dominant bryophyte lineages is unknown, which we examined by characterising the sole KANADI ortholog, MpKAN, in the liverwort Marchantia polymorpha. In contrast to angiosperm orthologs, MpKAN functions in the gametophyte generation in Marchantia, where it regulates apical branching and tissue differentiation, but does not influence tissue polarity in either generation. MpKAN can partially rescue the sporophyte polarity defects of kanadi mutants in Arabidopsis, indicating that MpKAN has conserved biochemical activity to its angiosperm counterparts. Mpkan loss-of-function plants display defects in far-red (FR) light responses. Mpkan plants have reduced FR-induced growth tropisms, have a delayed transition to sexual reproduction and fail to correctly form gametangiophores. Our results indicate that MpKAN is a modulator of FR responses, which may reflect a conserved role for KANADI across land plants. Under FR, MpKAN negatively regulates MpDELLA expression, suggesting that MpKAN and MpDELLA act in a pathway regulating FR responses, placing MpKAN in a gene regulatory network exhibiting similarities with those of angiosperms.

Marchantia polymorpha model reveals conserved infection mechanisms in the vascular wilt fungal pathogen Fusarium oxysporum

Root-infecting vascular fungi cause wilt diseases and provoke devastating losses in hundreds of crops. It is currently unknown how these pathogens evolved and whether they can also infect nonvascular plants, which diverged from vascular plants over 450 million years ago. We established a pathosystem between the nonvascular plant Marchantia polymorpha (Mp) and the root-infecting vascular wilt fungus Fusarium oxysporum (Fo). On angiosperms, Fo exhibits exquisite adaptation to the plant xylem niche as well as host-specific pathogenicity, both of which are conferred by effectors encoded on lineage-specific chromosomes. Fo isolates displaying contrasting lifestyles on angiosperms - pathogenic vs endophytic - are able to infect Mp and cause tissue maceration and host cell killing. Using isogenic fungal mutants we define a set of conserved fungal pathogenicity factors, including mitogen activated protein kinases, transcriptional regulators and cell wall remodelling enzymes, that are required for infection of both vascular and nonvascular plants. Markedly, two host-specific effectors and a morphogenetic regulator, which contribute to vascular colonisation and virulence on tomato plants are dispensable on Mp. Collectively, these findings suggest that vascular wilt fungi employ conserved infection strategies on nonvascular and vascular plant lineages but also have specific mechanisms to access the vascular niche of angiosperms.

The liverwort Marchantia polymorpha, a model for all ages

The liverwort Marchantia polymorpha has been known to man for millennia due to its inclusion Greek herbals. Perhaps due to its familiarity and association with growth in, often, man-made disturbed habitats, it was readily used to address fundamental biological questions of the day, including elucidation of land plant life cycles in the late 18th century, the formulation of cell theory early in the 19th century and the discovery of the alternation of generations in land plants in the mid-19th century. Subsequently, Marchantia was used as model in botany classes. With the arrival of the molecular era, its organellar genomes, the chloroplast and mitochondrial, were some of the first to be sequenced from any plant. In the past two decades, molecular genetic tools have been applied such that genes may be manipulated seemingly at will. Here, are past, present, and some views to the future of Marchantia as a model.

Evolution and functional diversification of R2R3-MYB transcription factors in plants

R2R3-MYB genes (R2R3-MYBs) form one of the largest transcription factor gene families in the plant kingdom, with substantial structural and functional diversity. However, the evolutionary processes leading to this amazing functional diversity have not yet been clearly established. Recently developed genomic and classical molecular technologies have provided detailed insights into the evolutionary relationships and functions of plant R2R3-MYBs. Here, we review recent genome-level and functional analyses of plant R2R3-MYBs, with an emphasis on their evolution and functional diversification. In land plants, this gene family underwent a large expansion by whole genome duplications and small-scale duplications. Along with this population explosion, a series of functionally conserved or lineage-specific subfamilies/groups arose with roles in three major plant-specific biological processes: development and cell differentiation, specialized metabolism, and biotic and abiotic stresses. The rapid expansion and functional diversification of plant R2R3-MYBs are highly consistent with the increasing complexity of angiosperms. In particular, recently derived R2R3-MYBs with three highly homologous intron patterns (a, b, and c) are disproportionately related to specialized metabolism and have become the predominant subfamilies in land plant genomes. The evolution of plant R2R3-MYBs is an active area of research, and further studies are expected to improve our understanding of the evolution and functional diversification of this gene family.

Molecular mechanisms involved in functional macroevolution of plant transcription factors

Transcription factors (TFs) are key components of the transcriptional regulation machinery. In plants, they accompanied the evolution from unicellular aquatic algae to complex flowering plants that dominate the land environment. The adaptations of the body plan and physiological responses required changes in the biological functions of TFs. Some ancestral gene regulatory networks are highly conserved, while others evolved more recently and only exist in particular lineages. The recent emergence of novel model organisms provided the opportunity for comparative studies, producing new insights to infer these evolutionary trajectories. In this review, we comprehensively revisit the recent literature on TFs of nonseed plants and algae, focusing on the molecular mechanisms driving their functional evolution. We discuss the particular contribution of changes in DNA-binding specificity, protein-protein interactions and cis-regulatory elements to gene regulatory networks. Current advances have shown that these evolutionary processes were shaped by changes in TF expression pattern, not through great innovation in TF protein sequences. We propose that the role of TFs associated with environmental and developmental regulation was unevenly conserved during land plant evolution.

Developmental Stage Determines the Accumulation Pattern of UV-Absorbing Compounds in the Model Liverwort Marchantia polymorpha subsp. ruderalis under Controlled Conditions

The liverwort Marchantia polymorpha subsp. ruderalis is an emerging model plant, and some data are available on its responses to ultraviolet (UV) radiation. However, it is unknown if the developmental stage of the thalli modulates the effects of UV radiation on the contents of potentially protecting phenolic compounds. To fill this gap, liverwort samples were exposed or non-exposed to UV radiation for 38 days under controlled conditions, using three developmental stages: gemmae (G), one-month thalli (T1), and two-month thalli (T2). Then, the bulk level of methanol-soluble UV-absorbing compounds and the contents of six flavones (apigenin and luteolin derivatives) were measured. The UV responsiveness decreased with thallus age: G and T1 plants were the most UV-responsive and showed a strong increase in all the variables, with G plants more responsive than T1 plants. In UV-exposed T2 plants, only apigenin derivatives increased and more modestly, probably due to a lower acclimation capacity. Nevertheless, the thalli became progressively tougher due to a decreasing water content, representing a possible structural protection against UV. In UV-exposed plants, the temporal patterns of the accumulation of phenolic compounds were compound-specific. Most compounds decreased with thallus age, but di-glucuronide derivatives showed a bell-shaped pattern, with T1 plants showing the highest contents. A Principal Components Analysis (PCA) ordination of the different samples summarized the results found. The patterns described above should be taken into account to select thalli of an adequate developmental stage for experiments investigating the induction of phenolic compounds by UV radiation.

The liverwort oil body is formed by redirection of the secretory pathway

Eukaryotic cells acquired novel organelles during evolution through mechanisms that remain largely obscure. The existence of the unique oil body compartment is a synapomorphy of liverworts that represents lineage-specific acquisition of this organelle during evolution, although its origin, biogenesis, and physiological function are yet unknown. We find that two paralogous syntaxin-1 homologs in the liverwort Marchantia polymorpha are distinctly targeted to forming cell plates and the oil body, suggesting that these structures share some developmental similarity. Oil body formation is regulated by an ERF/AP2-type transcription factor and loss of the oil body increases M. polymorpha herbivory. These findings highlight a common strategy for the acquisition of organelles with distinct functions in plants, via periodical redirection of the secretory pathway depending on cellular phase transition.

Transcriptional and Morpho-Physiological Responses of Marchantia polymorpha upon Phosphate Starvation

Phosphate (Pi) is a pivotal nutrient that constraints plant development and productivity in natural ecosystems. Land colonization by plants, more than 470 million years ago, evolved adaptive mechanisms to conquer Pi-scarce environments. However, little is known about the molecular basis underlying such adaptations at early branches of plant phylogeny. To shed light on how early divergent plants respond to Pi limitation, we analyzed the morpho-physiological and transcriptional dynamics of Marchantia polymorpha upon Pi starvation. Our phylogenomic analysis highlights some gene networks present since the Chlorophytes and others established in the Streptophytes (e.g., PHR1–SPX1 and STOP1–ALMT1, respectively). At the morpho-physiological level, the response is characterized by the induction of phosphatase activity, media acidification, accumulation of auronidins, reduction of internal Pi concentration, and developmental modifications of rhizoids. The transcriptional response involves the induction of MpPHR1, Pi transporters, lipid turnover enzymes, and MpMYB14, which is an essential transcription factor for auronidins biosynthesis. MpSTOP2 up-regulation correlates with expression changes in genes related to organic acid biosynthesis and transport, suggesting a preference for citrate exudation. An analysis of MpPHR1 binding sequences (P1BS) shows an enrichment of this cis regulatory element in differentially expressed genes. Our study unravels the strategies, at diverse levels of organization, exerted by M. polymorpha to cope with low Pi availability.

Plant Evolution: Divergent Plants, Divergent Functions for C1HDZ Orthologs

Ortholog identification inferred by phylogenetic analyses does not always correlate with functional conservation. The recent functional characterization of the C1HDZ transcription factor in the early-diverging land plant Marchantia polymorpha reveals its role in biotic stress responses, contrary to its orthologs in flowering plants.

Systematic Tools for Reprogramming Plant Gene Expression in a Simple Model, Marchantia polymorpha

We present the OpenPlant toolkit, a set of interlinked resources and techniques to develop Marchantia as testbed for bioengineering in plants. Marchantia is a liverwort, a simple plant with an open form of development that allows direct visualization of gene expression and dynamics of cellular growth in living tissues. We describe new techniques for simple and efficient axenic propagation and maintenance of Marchantia lines with no requirement for glasshouse facilities. Marchantia plants spontaneously produce clonal propagules within a few weeks of regeneration, and lines can be amplified million-fold in a single generation by induction of the sexual phase of growth, crossing, and harvesting of progeny spores. The plant has a simple morphology and genome with reduced gene redundancy, and the dominant phase of its life cycle is haploid, making genetic analysis easier. We have built robust Loop assembly vector systems for nuclear and chloroplast transformation and genome editing. These have provided the basis for building and testing a modular library of standardized DNA elements with highly desirable properties. We have screened transcriptomic data to identify a range of candidate genes, extracted putative promoter sequences, and tested them in vivo to identify new constitutive promoter elements. The resources have been combined into a toolkit for plant bioengineering that is accessible for laboratories without access to traditional facilities for plant biology research. The toolkit is being made available under the terms of the OpenMTA and will facilitate the establishment of common standards and the use of this simple plant as testbed for synthetic biology.

The Evolution of Flavonoid Biosynthesis: A Bryophyte Perspective

The flavonoid pathway is one of the best characterized specialized metabolite pathways of plants. In angiosperms, the flavonoids have varied roles in assisting with tolerance to abiotic stress and are also key for signaling to pollinators and seed dispersal agents. The pathway is thought to be specific to land plants and to have arisen during the period of land colonization around 550-470 million years ago. In this review we consider current knowledge of the flavonoid pathway in the bryophytes, consisting of the liverworts, hornworts, and mosses. The pathway is less characterized for bryophytes than angiosperms, and the first genetic and molecular studies on bryophytes are finding both commonalities and significant differences in flavonoid biosynthesis and pathway regulation between angiosperms and bryophytes. This includes biosynthetic pathway branches specific to each plant group and the apparent complete absence of flavonoids from the hornworts.

Auronidins are a previously unreported class of flavonoid pigments that challenges when anthocyanin biosynthesis evolved in plants

Anthocyanins are key pigments of plants, providing color to flowers, fruit, and foliage and helping to counter the harmful effects of environmental stresses. It is generally assumed that anthocyanin biosynthesis arose during the evolutionary transition of plants from aquatic to land environments. Liverworts, which may be the closest living relatives to the first land plants, have been reported to produce red cell wall-bound riccionidin pigments in response to stresses such as UV-B light, drought, and nutrient deprivation, and these have been proposed to correspond to the first anthocyanidins present in early land plant ancestors. Taking advantage of the liverwort model species Marchantia polymorpha, we show that the red pigments of Marchantia are formed by a phenylpropanoid biosynthetic branch distinct from that leading to anthocyanins. They constitute a previously unreported flavonoid class, for which we propose the name "auronidin," with similar colors as anthocyanin but different chemistry, including strong fluorescence. Auronidins might contribute to the remarkable ability of liverworts to survive in extreme environments on land, and their discovery calls into question the possible pigment status of the first land plants.

Mining MYB transcription factors from the genomes of orchids (Phalaenopsis and Dendrobium) and characterization of an orchid R2R3-MYB gene involved in water-soluble polysaccharide biosynthesis

Members of the MYB superfamily act as regulators in a wide range of biological processes in plants. Despite this, the MYB superfamily from the Orchidaceae has not been identified, and MYB genes related to bioactive water-soluble polysaccharide (WSP) biosynthesis are relatively unknown. In this study, we identified 159 and 165 MYB genes from two orchids, Phalaenopsis equestris and Dendrobium officinale, respectively. The MYB proteins were classified into four MYB classes in both orchids: MYB-related (MYBR), R2R3-MYB, 3R-MYB and atypical MYB proteins. The MYBR proteins in both orchids were classified into five subfamilies and 12 genes were strongly up-regulated in response to cold stress in D. officinale. The R2R3-MYB proteins were both divided into 31 clades in P. equestris and D. officinale. Among these clades, nine contained MYB TFs related to secondary cell wall biosynthesis or testa mucilage biosynthesis in Arabidopsis thaliana. In D. officinale, 10 candidate genes showed an expression pattern corresponding to changes in the WSP content. Overexpression of one of these candidate genes (DoMYB75) in A. thaliana increased seed WSP content by about 14%. This study provides information about MYB genes in two orchids that will further help to understand the transcriptional regulation of WSP biosynthesis in these orchids as well as other plant species.

Mining MYB transcription factors from the genomes of orchids (Phalaenopsis and Dendrobium) and characterization of an orchid R2R3-MYB gene involved in water-soluble polysaccharide biosynthesis

Members of the MYB superfamily act as regulators in a wide range of biological processes in plants. Despite this, the MYB superfamily from the Orchidaceae has not been identified, and MYB genes related to bioactive water-soluble polysaccharide (WSP) biosynthesis are relatively unknown. In this study, we identified 159 and 165 MYB genes from two orchids, Phalaenopsis equestris and Dendrobium officinale, respectively. The MYB proteins were classified into four MYB classes in both orchids: MYB-related (MYBR), R2R3-MYB, 3R-MYB and atypical MYB proteins. The MYBR proteins in both orchids were classified into five subfamilies and 12 genes were strongly up-regulated in response to cold stress in D. officinale. The R2R3-MYB proteins were both divided into 31 clades in P. equestris and D. officinale. Among these clades, nine contained MYB TFs related to secondary cell wall biosynthesis or testa mucilage biosynthesis in Arabidopsis thaliana. In D. officinale, 10 candidate genes showed an expression pattern corresponding to changes in the WSP content. Overexpression of one of these candidate genes (DoMYB75) in A. thaliana increased seed WSP content by about 14%. This study provides information about MYB genes in two orchids that will further help to understand the transcriptional regulation of WSP biosynthesis in these orchids as well as other plant species.

A gain-of-function T-DNA insertion mutant of Marchantia polymorpha hyper-accumulates flavonoid riccionidin A

Marchantia polymorpha is a model liverwort for which many molecular biological techniques are now available. We previously developed the S-AgarTrap method for easy genetic transformation of M. polymorpha using spores. In this study, we report production of a T-DNA insertion mutant library (approx. 10,000 lines) for M. polymorpha using the S-AgarTrap method. We further isolate and characterize a gain-of-function mutant that hyper-accumulates the flavonoid riccionidin A. The present study demonstrates that the S-AgarTrap-mediated production of a T-DNA insertion mutant library is a powerful tool for molecular biology in M. polymorpha.

Conserved Biochemical Defenses Underpin Host Responses to Oomycete Infection in an Early-Divergent Land Plant Lineage

The expansion of plants onto land necessitated the evolution of robust defense strategies to protect against a wide array of microbial invaders. Whereas host responses to microbial colonization are extensively explored in evolutionarily young land plant lineages such as angiosperms, we know relatively little about plant-pathogen interactions in early-diverging land plants thought to better represent the ancestral state. Here, we define the transcriptional and proteomic response of the early-divergent liverwort Marchantia polymorpha to infection with the oomycete pathogen Phytophthora palmivora. We uncover a robust molecular response to oomycete colonization in Marchantia that consists of conserved land plant gene families. Direct macroevolutionary comparisons of host infection responses in Marchantia and the model angiosperm Nicotiana benthamiana further reveal a shared set of orthologous microbe-responsive genes that include members of the phenylpropanoid metabolic pathway. In addition, we identify a role for the Marchantia R2R3-MYB transcription factor MpMyb14 in activating phenylpropanoid (flavonoid) biosynthesis during oomycete infection. Mpmyb14 mutants infected with P. palmivora fail to activate phenylpropanoid biosynthesis gene expression and display enhanced disease susceptibility compared to wild-type plants. Conversely, the ectopic induction of MpMyb14 led to the accumulation of anthocyanin-like pigments and dramatically enhanced liverwort resistance to P. palmivora infection. Collectively, our results demonstrate that the Marchantia response to oomycete infection displays evolutionarily conserved features indicative of an ancestral pathogen deterrence strategy centered on phenylpropanoid-mediated biochemical defenses.

Physiological function of photoreceptor UVR8 in UV-B tolerance in the liverwort Marchantia polymorpha

The physiological importance of MpUVR8 in UV-B resistance and translocation in a UV-B-dependent manner from the cytosol into the nucleus is characterized in Marchantia polymorpha. UV RESISTANCE LOCUS 8 (UVR8) is an ultraviolet-B (UV-B) light receptor functioning for UV-B sensing and tolerance in Arabidopsis thaliana and other species. It is unclear whether UVR8 physiologically functions in UV-B-induced defense responses in Marchantia polymorpha, which belongs to the earliest diverging group of embryophyte lineages. Here, we demonstrate that UVR8 has a physiological function in UV-B tolerance and that there is a UVR8-dependent pathway involved. In addition, a UVR8-independent pathway is revealed. We examine the tissue-specific expression pattern of M. polymorpha UVR8 (MpUVR8), showing that it is highly expressed in the apical notch in thalli and gametangiophores, as well as in antheridial and archegonial heads. Furthermore, Mpuvr8^KO plant transformants, in which the MpUVR8 locus was disrupted, were produced and analyzed to understand the physiological and molecular function of MpUVR8. Analysis using these plants indicates the important roles of MpUVR8 and MpUVR8-regulated genes, and of MpUVR8-independent pathways in UV-B tolerance. Subcellular localization of Citrine-fused MpUVR8 in M. polymorpha cells was also investigated. It was found to translocate from the cytosol into the nucleus in response to UV-B irradiation. Our findings indicate strong conservation of the physiological function of UVR8 and the molecular mechanisms for UVR8-dependent signal transduction through regulation of gene expression in embryophytes.

Genome-wide analysis of MpBHLH12, a IIIf basic helix-loop-helix transcription factor of Marchantia polymorpha

The evolution of plants on land required adaptation to UV radiation and dry environments, and involved the appearance and/or rewiring of genetic connections, known as gene regulatory networks (GRNs), which consist of one or more transcription factors (TFs). The liverwort, Marchantia polymorpha, is a basal land plant, with a recently sequenced genome. The number of genes encoding basic helix-loop-helix (bHLH) family members is considerably higher in M. polymorpha than in charophyte green algae, suggesting the contribution of bHLH proteins to the evolution of GRNs associated with the adaptation of plants to land. Although an understanding of the evolutionary aspects of GRNs is fundamental for elucidating the mechanisms of environmental adaptation, the evolution of GRNs that led to land adaptation in plants remains poorly understood. In this study, we isolated a single gene encoding a IIIf bHLH TF from M. polymorpha, MpBHLH12. Transgenic M. polymorpha constitutively overexpressing MpBHLH12 showed smaller and fewer gemma cups than wild type, suggesting that MpBHLH12 is involved in the regulation of morphological development. Transcriptomic analysis of MpBHLH12 overexpressor (MpBHLH12ox) lines revealed an overlap with the GRN of MpMYB14, which regulates the biosynthesis of anthocyanins and phenolic compounds. However, MpBHLH12ox did not show anthocyanin accumulation. Results of the transient reporter assay suggest that MpBHLH12 could function in repression rather than activation. Our findings suggest that although the IIIf bHLH MpBHLH12 shows highest amino acid similarity with IIIf bHLH clade and is involved in developmental process and partly biosynthesis of phenolic compounds in M. polymorpha like Arabidopsis IIIf bHLH, the GRN involving MpBHLH12 would be distinct one from those of the IIIf bHLH TFs of seed plants.

The Origin and Evolution of Plant Flavonoid Metabolism

During their evolution, plants have acquired the ability to produce a huge variety of compounds. Unlike the specialized metabolites that accumulate in limited numbers of species, flavonoids are widely distributed in the plant kingdom. Therefore, a detailed analysis of flavonoid metabolism in genomics and metabolomics is an ideal way to investigate how plants have developed their unique metabolic pathways during the process of evolution. More comprehensive and precise metabolite profiling integrated with genomic information are helpful to emerge unexpected gene functions and/or pathways. The distribution of flavonoids and their biosynthetic genes in the plant kingdom suggests that flavonoid biosynthetic pathways evolved through a series of steps. The enzymes that form the flavonoid scaffold structures probably first appeared by recruitment of enzymes from primary metabolic pathways, and later, enzymes that belong to superfamilies such as 2-oxoglutarate-dependent dioxygenase, cytochrome P450, and short-chain dehydrogenase/reductase modified and varied the structures. It is widely accepted that the first two enzymes in flavonoid biosynthesis, chalcone synthase, and chalcone isomerase, were derived from common ancestors with enzymes in lipid metabolism. Later enzymes acquired their function by gene duplication and the subsequent acquisition of new functions. In this review, we describe the recent progress in metabolomics technologies for flavonoids and the evolution of flavonoid skeleton biosynthetic enzymes to understand the complicate evolutionary traits of flavonoid metabolism in plant kingdom.