KATANIN-mediated microtubule severing is required for MTOC organisation and function in Marchantia polymorpha

Microtubule organising centres (MTOCs) are sites of localised microtubule nucleation in eukaryotic cells. Regulation of microtubule dynamics often involves KATANIN (KTN): a microtubule severing enzyme that cuts microtubules to generate new negative ends, leading to catastrophic depolymerisation. In Arabidopsis thaliana, KTN is required for the organisation of microtubules in the cell cortex, preprophase band, mitotic spindle and phragmoplast. However, as angiosperms lack MTOCs, the role of KTN in MTOC formation has yet to be studied in plants. Two unique MTOCs - the polar organisers - form on opposing sides of the preprophase nucleus in liverworts. Here, we show that KTN-mediated microtubule depolymerisation regulates the number and organisation of polar organisers formed in Marchantia polymorpha. Mpktn mutants that lacked KTN function had supernumerary disorganised polar organisers compared with wild type. This was in addition to defects in the microtubule organisation in the cell cortex, preprophase band, mitotic spindle and phragmoplast. These data are consistent with the hypothesis that KTN-mediated microtubule dynamics are required for the de novo formation of MTOCs, a previously unreported function in plants.

Elucidation of arsenic detoxification mechanism in Marchantia polymorpha: The role of ACR3

The arsenic-specific ACR3 transporter plays pivotal roles in As detoxification in yeast and a group of ancient tracheophytes, the ferns. Despite putative ACR3 genes being present in the genomes of bryophytes, whether they have the same relevance also in this lineage is currently unknown. In this study, we characterized the MpACR3 gene from the bryophyte Marchantia polymorpha L. through a multiplicity of functional approaches ranging from phylogenetic reconstruction, expression analysis, loss- and gain-of-function as well as genetic complementation with an MpACR3 gene tagged with a fluorescent protein. Genetic complementation demonstrates that MpACR3 plays a pivotal role in As tolerance in M. polymorpha, with loss-of-function Mpacr3 mutants being hypersensitive and MpACR3 overexpressors more tolerant to As. Additionally, MpACR3 activity regulates intracellular As concentration, affects its speciation and controls the levels of intracellular oxidative stress. The MpACR3::3xCitrine appears to localize at the plasma membrane and possibly in other endomembrane systems. Taken together, these results demonstrate the pivotal function of ACR3 detoxification in both sister lineages of land plants, indicating that it was present in the common ancestor to all embryophytes. We propose that Mpacr3 mutants could be used in developing countries as low-cost and low-technology visual bioindicators to detect As pollution in water.

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.

Myosin XI, a model of its conserved role in plant cell tip growth

In eukaryotic cells, organelle and vesicle transport, positioning, and interactions play crucial roles in cytoplasmic organization and function. These processes are governed by intracellular trafficking mechanisms. At the core of that trafficking, the cytoskeleton and directional transport by motor proteins stand out as its key regulators. Plant cell tip growth is a well-studied example of cytoplasm organization by polarization. This polarization, essential for the cell's function, is driven by the cytoskeleton and its associated motors. This review will focus on myosin XI, a molecular motor critical for vesicle trafficking and polarized plant cell growth. We will center our discussion on recent data from the moss Physcomitrium patens and the liverwort Marchantia polymorpha. The biochemical properties and structure of myosin XI in various plant species are discussed, highlighting functional conservation across species. We further explore this conservation of myosin XI function in the process of vesicle transport in tip-growing cells. Existing evidence indicates that myosin XI actively organizes actin filaments in tip-growing cells by a mechanism based on vesicle clustering at their tips. A hypothetical model is presented to explain the essential function of myosin XI in polarized plant cell growth based on vesicle clustering at the tip. The review also provides insight into the in vivo localization and dynamics of myosin XI, emphasizing its role in cytosolic calcium regulation, which influences the polymerization of F-actin. Lastly, we touch upon the need for additional research to elucidate the regulation of myosin function.

Efficient Production of Flavonoid Glucuronides in Escherichia coli Using Flavonoid O-Glucuronosyltransferases Characterized from Marchantia polymorpha

As a model liverwort, Marchantia polymorpha contains various flavone glucuronides with cardiovascular-promoting effects and anti-inflammatory properties. However, the related glucuronosyltransferases have not yet been reported. In this study, two bifunctional UDP-glucuronic acid/UDP-glucose:flavonoid glucuronosyltransferases/glucosyltransferases, MpUGT742A1 and MpUGT736B1, were identified from M. polymorpha. Extensive enzymatic assays found that MpUGT742A1 and MpUGT736B1 exhibited efficient glucuronidation activity for flavones, flavonols, and flavanones and showed promiscuous regioselectivity at positions 3, 6, 7, 3', and 4'. These enzymes catalyzed the production of a variety of flavonoid glucuronides with medicinal value, including apigenin-7-O-glucuronide and scutellarein-7-O-glucuronide. With the use of MpUGT736B1, apigenin-4'-O-glucuronide and apigenin-7,4'-di-O-glucuronide were prepared by scaled-up enzymatic catalysis and structurally identified by NMR spectroscopy. MpUGT742A1 also displayed glucosyltransferase activity on the 7-OH position of the flavanones using UDP-glucose as the sugar donor. Furthermore, we constructed four recombinant strains by combining the pathway for increasing the UDP-glucuronic acid supply with the two novel UGTs MpUGT742A1 and MpUGT736B1. When apigenin was used as a substrate, the extracellular apigenin-4'-O-glucuronide and apigenin-7,4'-di-O-glucuronide production obtained from the Escherichia coli strain BB2 reached 598 and 81 mg/L, respectively. Our study provides new candidate genes and strategies for the biosynthesis of flavonoid glucuronides.

The D-mannose/L-galactose pathway plays a predominant role in ascorbate biosynthesis in the liverwort Marchantia polymorpha but is not regulated by light and oxidative stress

Ascorbate plays an indispensable role in plants, functioning as both an antioxidant and a cellular redox buffer. It is widely acknowledged that the ascorbate biosynthesis in the photosynthetic tissues of land plants is governed by light-mediated regulation of the D-mannose/L-galactose (D-Man/L-Gal) pathway. At the core of this light-dependent regulation lies the VTC2 gene, encoding the rate-limiting enzyme GDP-L-Gal phosphorylase. The VTC2 expression is regulated by signals via the photosynthetic electron transport system. In this study, we directed our attention to the liverwort Marchantia polymorpha, representing one of the basal land plants, enabling us to conduct an in-depth analysis of its ascorbate biosynthesis. The M. polymorpha genome harbors a solitary gene for each enzyme involved in the D-Man/L-Gal pathway, including VTC2, along with three lactonase orthologs, which may be involved in the alternative ascorbate biosynthesis pathway. Through supplementation experiments with potential precursors, we observed that only L-Gal exhibited effectiveness in ascorbate biosynthesis. Furthermore, the generation of VTC2-deficient mutants through genome editing unveiled the inability of thallus regeneration in the absence of L-Gal supplementation, thereby revealing the importance of the D-Man/L-Gal pathway in ascorbate biosynthesis within M.  polymorpha. Interestingly, gene expression analyses unveiled a distinct characteristic of M. polymorpha, where none of the genes associated with the D-Man/L-Gal pathway, including VTC2, showed upregulation in response to light, unlike other known land plants. This study sheds light on the exceptional nature of M. polymorpha as a land plant that has evolved distinctive mechanisms concerning ascorbate biosynthesis and its regulation.

Plastid terminal oxidase (PTOX) protects photosystem I and not photosystem II against photoinhibition in Arabidopsis thaliana and Marchantia polymorpha

The plastid terminal oxidase PTOX controls the oxidation level of the plastoquinone pool in the thylakoid membrane and acts as a safety valve upon abiotic stress, but detailed characterization of its role in protecting the photosynthetic apparatus is limited. Here we used PTOX mutants in two model plants Arabidopsis thaliana and Marchantia polymorpha. In Arabidopsis, lack of PTOX leads to a severe defect in pigmentation, a so-called variegated phenotype, when plants are grown at standard light intensities. We created a green Arabidopsis PTOX mutant expressing the bacterial carotenoid desaturase CRTI and a double mutant in Marchantia lacking both PTOX isoforms, the plant-type and the alga-type PTOX. In both species, lack of PTOX affected the redox state of the plastoquinone pool. Exposure of plants to high light intensity showed in the absence of PTOX higher susceptibility of photosystem I to light-induced damage while photosystem II was more stable compared with the wild type demonstrating that PTOX plays both, a pro-oxidant and an anti-oxidant role in vivo. Our results shed new light on the function of PTOX in the protection of photosystem I and II.

MiRNAs differentially expressed in vegetative and reproductive organs of Marchantia polymorpha - insights into their expression pattern, gene structures and function

MicroRNAs regulate gene expression affecting a variety of plant developmental processes. The evolutionary position of Marchantia polymorpha makes it a significant model to understand miRNA-mediated gene regulatory pathways in plants. Previous studies focused on conserved miRNA-target mRNA modules showed their critical role in Marchantia development. Here, we demonstrate that the differential expression of conserved miRNAs among land plants and their targets in selected organs of Marchantia additionally underlines their role in regulating fundamental developmental processes. The main aim of this study was to characterize selected liverwort-specific miRNAs, as there is a limited knowledge on their biogenesis, accumulation, targets, and function in Marchantia. We demonstrate their differential accumulation in vegetative and generative organs. We reveal that all liverwort-specific miRNAs examined are encoded by independent transcriptional units. MpmiR11737a, MpmiR11887 and MpmiR11796, annotated as being encoded within protein-encoding genes, have their own independent transcription start sites. The analysis of selected liverwort-specific miRNAs and their pri-miRNAs often reveal correlation in their levels, suggesting transcriptional regulation. However, MpmiR11796 shows a reverse correlation to its pri-miRNA level, suggesting post-transcriptional regulation. Moreover, we identify novel targets for selected liverwort-specific miRNAs and demonstrate an inverse correlation between their expression and miRNA accumulation. In the case of one miRNA precursor, we provide evidence that it encodes two functional miRNAs with two independent targets. Overall, our research sheds light on liverwort-specific miRNA gene structure, provides new data on their biogenesis and expression regulation. Furthermore, identifying their targets, we hypothesize the potential role of these miRNAs in early land plant development and functioning.

Control of stem cell behavior by CLE-JINGASA signaling in the shoot apical meristem in Marchantia polymorpha

Land plants undergo indeterminate growth by the activity of meristems in both gametophyte (haploid) and sporophyte (diploid) generations. In the sporophyte of the flowering plant Arabidopsis thaliana, the apical meristems are located at the shoot and root tips in which a number of regulatory gene homologs are shared for their development, implying deep evolutionary origins. However, little is known about their functional conservation with gametophytic meristems in distantly related land plants such as bryophytes, even though genomic studies have revealed that the subfamily-level diversity of regulatory genes is mostly conserved throughout land plants. Here, we show that a NAM/ATAF/CUC (NAC) domain transcription factor, JINGASA (MpJIN), acts downstream of CLAVATA3 (CLV3)/ESR-related (CLE) peptide signaling and controls stem cell behavior in the gametophytic shoot apical meristem of the liverwort Marchantia polymorpha. In the meristem, strong MpJIN expression was associated with the periclinal cell division at the periphery of the stem cell zone (SCZ), whereas faint MpJIN expression was found at the center of the SCZ. Time course observation indicates that the MpJIN-negative cells are lost from the SCZ and respecified de novo at two separate positions during the dichotomous branching event. Consistently, the induction of MpJIN results in ectopic periclinal cell division in the SCZ and meristem termination. Based on the comparative expression data, we speculate that the function of JIN/FEZ subfamily genes was shared among the shoot apical meristems in the gametophyte and sporophyte generations in early land plants but was lost in certain lineages, including the flowering plant A. thaliana.

Conserved and non-conserved RNA-target modules in plants: lessons for a better understanding of Marchantia development

A wide variety of functional regulatory non-coding RNAs (ncRNAs) have been identified as essential regulators of plant growth and development. Depending on their category, ncRNAs are not only involved in modulating target gene expression at the transcriptional and post-transcriptional levels but also are involved in processes like RNA splicing and RNA-directed DNA methylation. To fulfill their molecular roles properly, ncRNAs must be precisely processed by multiprotein complexes. In the case of small RNAs, DICER-LIKE (DCL) proteins play critical roles in the production of mature molecules. Land plant genomes contain at least four distinct classes of DCL family proteins (DCL1-DCL4), of which DCL1, DCL3 and DCL4 are also present in the genomes of bryophytes, indicating the early divergence of these genes. The liverwort Marchantia polymorpha has become an attractive model species for investigating the evolutionary history of regulatory ncRNAs and proteins that are responsible for ncRNA biogenesis. Recent studies on Marchantia have started to uncover the similarities and differences in ncRNA production and function between the basal lineage of bryophytes and other land plants. In this review, we summarize findings on the essential role of regulatory ncRNAs in Marchantia development. We provide a comprehensive overview of conserved ncRNA-target modules among M. polymorpha, the moss Physcomitrium patens and the dicot Arabidopsis thaliana, as well as Marchantia-specific modules. Based on functional studies and data from the literature, we propose new connections between regulatory pathways involved in Marchantia's vegetative and reproductive development and emphasize the need for further functional studies to understand the molecular mechanisms that control ncRNA-directed developmental processes.

Biosynthesis of gibberellin-related compounds modulates far-red light responses in the liverwort Marchantia polymorpha

Gibberellins (GAs) are key phytohormones that regulate growth, development, and environmental responses in angiosperms. From an evolutionary perspective, all major steps of GA biosynthesis are conserved among vascular plants, while GA biosynthesis intermediates such as ent-kaurenoic acid (KA) are also produced by bryophytes. Here, we show that in the liverwort Marchantia polymorpha, KA and GA12 are synthesized by evolutionarily conserved enzymes, which are required for developmental responses to far-red light (FR). Under FR-enriched conditions, mutants of various biosynthesis enzymes consistently exhibited altered thallus growth allometry, delayed initiation of gametogenesis, and abnormal morphology of gamete-bearing structures (gametangiophores). By chemical treatments and liquid chromatography-tandem mass spectrometry (LC-MS/MS) analyses, we confirmed that these phenotypes were caused by the deficiency of some GA-related compounds derived from KA, but not bioactive GAs from vascular plants. Transcriptome analysis showed that FR enrichment induced the up-regulation of genes related to stress responses and secondary metabolism in M. polymorpha, which was largely dependent on the biosynthesis of GA-related compounds. Due to the lack of canonical GA receptors in bryophytes, we hypothesize that GA-related compounds are commonly synthesized in land plants but were co-opted independently to regulate responses to light quality change in different plant lineages during the past 450 million years of evolution.

An enhancer trap system to track developmental dynamics in Marchantia polymorpha

A combination of streamlined genetics, experimental tractability and relative morphological simplicity compared to vascular plants makes the liverwort Marchantia polymorpha an ideal model system for studying many aspects of plant biology. Here we describe a transformation vector combining a constitutive fluorescent membrane marker with a nuclear marker that is regulated by nearby enhancer elements and use this to produce a library of enhancer trap lines for Marchantia. Screening gemmae from these lines allowed the identification and characterization of novel marker lines, including markers for rhizoids and oil cells. The library allowed the identification of a margin tissue running around the thallus edge, highlighted during thallus development. The expression of this marker is correlated with auxin levels. We generated multiple markers for the meristematic apical notch region, which have different spatial expression patterns, reappear at different times during meristem regeneration following apical notch excision and have varying responses to auxin supplementation or inhibition. This reveals that there are proximodistal substructures within the apical notch that could not be observed otherwise. We employed our markers to study Marchantia sporeling development, observing meristem emergence as defining the protonema-to-prothallus stage transition, and subsequent production of margin tissue during the prothallus stage. Exogenous auxin treatment stalls meristem emergence at the protonema stage but does not inhibit cell division, resulting in callus-like sporelings with many rhizoids, whereas pharmacologically inhibiting auxin synthesis and transport does not prevent meristem emergence. This enhancer trap system presents a useful resource for the community and will contribute to future Marchantia research.

Insights into cryptochrome modulation of ABA signaling to mediate dormancy regulation in Marchantia polymorpha

The acquisition of dormancy capabilities has enabled plants to survive in adverse terrestrial environmental conditions. Dormancy accumulation and release is coupled with light signaling, which is well studied in Arabidopsis, but it is unclear in the distant nonvascular relative. We study the characteristics and function on dormancy regulation of a blue light receptor cryptochrome in Marchantia polymorpha (MpCRY). Here, we identified MpCRY via bioinformatics and mutant complement analysis. The biochemical characteristics were assessed by multiple protein-binding assays. The function of MpCRY in gemma dormancy was clarified by overexpression and mutation of MpCRY, and its mechanism was analyzed via RNA sequencing and quantitative PCR analyses associated with hormone treatment. We found that the unique MpCRY protein in M. polymorpha undergoes both blue light-promoted interaction with itself (self-interaction) and blue light-dependent phosphorylation. MpCRY has the specific characteristics of blue light-induced nuclear localization and degradation. We further demonstrated that MpCRY transcriptionally represses abscisic acid (ABA) signaling-related gene expression to suppress gemma dormancy, which is dependent on blue light signaling. Our findings indicate that MpCRY possesses specific biochemical and molecular characteristics, and modulates ABA signaling under blue light conditions to regulate gemma dormancy in M. polymorpha.

Transcriptome changes in chlorsulfuron-treated plants are caused by acetolactate synthase inhibition and not induction of a herbicide detoxification system in Marchantia polymorpha

A sensing mechanism in mammals perceives xenobiotics and induces the transcription of genes encoding proteins that detoxify these molecules. However, it is unclear if plants sense xenobiotics, and activate an analogous signalling system leading to their detoxification. Using the liverwort Marchantia polymorpha, we tested the hypothesis that there is a sensing system in plants that perceives herbicides resulting in the increased transcription of genes encoding proteins that detoxify these herbicides. Consistent with the hypothesis, we show that chlorsulfuron-treatment induces changes in the M. polymorpha transcriptome. However, these transcriptome changes do not occur in chlorsulfuron (CS)-treated target site resistant mutants, where the gene encoding the target carries a mutation that confers resistance to chlorsulfuron. Instead, we show that inactivation of the chlorsulfuron target, acetolactate synthase (ALS) (also known as acetohydroxyacid synthase (AHAS)), is required for the transcriptome response. These data demonstrate that the transcriptome changes in chlorsulfuron-treated plants are caused by disrupted amino acid synthesis and metabolism resulting from acetolactate synthase inhibition, and indicate that the transcriptome changes are not caused by a herbicide sensing mechanism.

Abscisic acid-mediated sugar responses are essential for vegetative desiccation tolerance in the liverwort Marchantia polymorpha

Low-molecular-weight sugars serve as protectants for cellular membranes and macromolecules under the condition of dehydration caused by environmental stress such as desiccation and freezing. These sugars also affect plant growth and development by provoking internal signaling pathways. We previously showed that both sugars and the stress hormone abscisic acid (ABA) enhance desiccation tolerance of gemma, a dormant propagule of the liverwort Marchantia polymorpha. To determine the role of ABA in sugar responses in liverworts, we generated genome-editing lines of M. polymorpha ABA DEFICIENT 1 (MpABA1) encoding zeaxanthin epoxidase, which catalyzes the initial reaction toward ABA biosynthesis. The generated Mpaba1 lines that accumulated only a trace amount of endogenous ABA showed reduced desiccation tolerance and reduced sugar responses. RNA-seq analysis of sucrose-treated gemmalings of M. polymorpha revealed that expression of a large part of sucrose-induced genes was reduced in Mpaba1 compared to the wild-type. Furthermore, Mpaba1 accumulated smaller amounts of low-molecular-weight sugars in tissues upon sucrose treatment than the wild-type, with reduced expression of genes for sucrose synthesis, sugar transporters, and starch-catabolizing enzymes. These results indicate that endogenous ABA plays a role in the regulation of the positive feedback loop for sugar-induced sugar accumulation in liverworts, enabling the tissue to have desiccation tolerance.

CIPK-B is essential for salt stress signalling in Marchantia polymorpha

Calcium signalling is central to many plant processes, with families of calcium decoder proteins having expanded across the green lineage and redundancy existing between decoders. The liverwort Marchantia polymorpha has fast become a new model plant, but the calcium decoders that exist in this species remain unclear. We performed phylogenetic analyses to identify the calcineurin B-like (CBL) and CBL-interacting protein kinase (CIPK) network of M. polymorpha. We analysed CBL-CIPK expression during salt stress, and determined protein-protein interactions using yeast two-hybrid and bimolecular fluorescence complementation. We also created genetic knockouts using CRISPR/Cas9. We confirm that M. polymorpha has two CIPKs and three CBLs. Both CIPKs and one CBL show pronounced salt-responsive transcriptional changes. All M. polymorpha CBL-CIPKs interact with each other in planta. Knocking out CIPK-B causes increased sensitivity to salt, suggesting that this CIPK is involved in salt signalling. We have identified CBL-CIPKs that form part of a salt tolerance pathway in M. polymorpha. Phylogeny and interaction studies imply that these CBL-CIPKs form an evolutionarily conserved salt overly sensitive pathway. Hence, salt responses may be some of the early functions of CBL-CIPK networks and increased abiotic stress tolerance required for land plant emergence.

All roads lead to Rome: alternative biosynthetic routes in plant specialised metabolism

This article is a Commentary on T‐T. Zhu et al. (2023), 237: 515–531.

A potential role of salicylic acid in the evolutionary behavior of Trichoderma as a plant pathogen: from Marchantia polymorpha to Arabidopsis thaliana

Recognition of the interaction of Trichoderma during the evolution of land plants plays a potential key role in the development of the salicylic acid defense pathway and the establishment of a mutualistic relationship. Marchantia polymorpha is a common liverwort considered in recent years as a model plant for evolutionary studies on plant-microorganism interactions. Despite the lack of research, remarkable results have been reported regarding the understanding of metabolic and evolutionary processes of beneficial and/or harmful interactions, owing to a better understanding of the origin and evolution of different plant defense pathways. In this study, we have carried out work on the direct and indirect interactions (exudates and volatiles) of M. polymorpha with different species of the fungal genus Trichoderma. These interactions showed different outcomes, including resistance or even growth promotion and disease. We have analyzed the level of tissue colonization and defense-related gene expression. Furthermore, we have used the pteridophyte Dryopteris affinis and the angiosperm Arabidopsis thaliana, as subsequent steps in plant evolution, together with the plant pathogen Rhizoctonia solani as a control of plant pathogenicity. Trichoderma virens, T. brevicompactum and T. hamatum are pathogens of M. polymorpha, while exudates of T. asperellum are harmful to the plant. The analysis of the expression of several defense genes in M. polymorpha and A. thaliana showed that there is a correlation of the transcriptional activation of SA-related genes with resistance or susceptibility of M. polymorpha to Trichoderma. Moreover, exogenous SA provides resistance to the virulent Trichoderma species. This beneficial fungus may have had an evolutionary period of interaction with plants in which it behaved as a plant pathogen until plants developed a defense system to limit its colonization through a defense response mediated by SA.

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.

Phylogenetic distribution and expression pattern analyses identified a divergent basal body assembly protein involved in land plant spermatogenesis

Land plant spermatozoids commonly possess characteristic structures such as the spline, which consists of a microtubule array, the multilayered structure (MLS) in which the uppermost layer is a continuum of the spline, and multiple flagella. However, the molecular mechanisms underpinning spermatogenesis remain to be elucidated. We successfully identified candidate genes involved in spermatogenesis, deeply divergent BLD10s, by computational analyses combining multiple methods and omics data. We then examined the functions of BLD10s in the liverwort Marchantia polymorpha and the moss Physcomitrium patens. MpBLD10 and PpBLD10 are required for normal basal body (BB) and flagella formation. Mpbld10 mutants exhibited defects in remodeling of the cytoplasm and nucleus during spermatozoid formation, and thus MpBLD10 should be involved in chromatin reorganization and elimination of the cytoplasm during spermiogenesis. We identified orthologs of MpBLD10 and PpBLD10 in diverse Streptophyta and found that MpBLD10 and PpBLD10 are orthologous to BLD10/CEP135 family proteins, which function in BB assembly. However, BLD10s evolved especially quickly in land plants and MpBLD10 might have acquired additional functions in spermatozoid formation through rapid molecular evolution.

Diversification of heat shock transcription factors expanded thermal stress responses during early plant evolution

The copy numbers of many plant transcription factor (TF) genes substantially increased during terrestrialization. This allowed TFs to acquire new specificities and thus create gene regulatory networks (GRNs) with new biological functions to help plants adapt to terrestrial environments. Through characterizing heat shock factor (HSF) genes MpHSFA1 and MpHSFB1 in the liverwort Marchantia polymorpha, we explored how heat-responsive GRNs widened their functions in M. polymorpha and Arabidopsis thaliana. An interspecies comparison of heat-induced transcriptomes and the evolutionary rates of HSFs demonstrated the emergence and subsequent rapid evolution of HSFB prior to terrestrialization. Transcriptome and metabolome analyses of M. polymorpha HSF-null mutants revealed that MpHSFA1 controls canonical heat responses such as thermotolerance and metabolic changes. MpHSFB1 also plays essential roles in heat responses, as well as regulating developmental processes including meristem branching and antheridiophore formation. Analysis of cis-regulatory elements revealed development- and stress-related TFs that function directly or indirectly downstream of HSFB. Male gametophytes of M. polymorpha showed higher levels of thermotolerance than female gametophytes, which could be explained by different expression levels of MpHSFA1U and MpHSFA1V on sex chromosome. We propose that the diversification of HSFs is linked to the expansion of HS responses, which enabled coordinated multicellular reactions in land plants.

Protein Kinase MpYAK1 Is Involved in Meristematic Cell Proliferation, Reproductive Phase Change and Nutrient Signaling in the Liverwort Marchantia polymorpha

Plant growth and development are regulated by environmental factors, including nutrient availability and light conditions, via endogenous genetic signaling pathways. Phosphorylation-dependent protein modification plays a major role in the regulation of cell proliferation in stress conditions, and several protein kinases have been shown to function in response to nutritional status, including dual-specificity tyrosine phosphorylation-regulated kinases (DYRKs). Although DYRKs are widely conserved in eukaryotes, the physiological functions of DYRKs in land plants are still to be elucidated. In the liverwort Marchantia polymorpha, a model bryophyte, four putative genes encoding DYRK homologous proteins, each of which belongs to the subfamily yet another kinase 1 (Yak1), plant-specific DYRK, DYRK2, or pre-mRNA processing protein 4 kinase, were identified. MpYAK1-defective male and female mutant lines generated by the clustered regularly interspaced short palindromic repeat (CRISPR)/CRISPR-associated nuclease 9 (Cas9) system showed smaller sizes of thalli than did the wild-type plants and repressed cell divisions in the apical notch regions. The Mpyak1 mutants developed rhizoids from gemmae in the gemma cup before release. The Mpyak1 lines developed sexual organs even in non-inductive short-day photoperiod conditions supplemented with far-red light. In nitrogen (N)-deficient conditions, rhizoid elongation was inhibited in the Mpyak1 mutants. In conditions of aeration with 0.08% CO2 (v/v) and N depletion, Mpyak1 mutants accumulated higher levels of sucrose and lower levels of starch compared to the wild type. Transcriptomic analyses revealed that the expression of peroxidase genes was differentially affected by MpYAK1. These results suggest that MpYAK1 is involved in the maintenance of plant growth and developmental responses to light conditions and nutrient signaling.

The wound-activated ERF15 transcription factor drives Marchantia polymorpha regeneration by activating an oxylipin biosynthesis feedback loop

The regenerative potential in response to wounding varies widely among species. Within the plant lineage, the liverwort Marchantia polymorpha displays an extraordinary regeneration capacity. However, its molecular pathways controlling the initial regeneration response are unknown. Here, we demonstrate that the MpERF15 transcription factor gene is instantly activated after wounding and is essential for gemmaling regeneration following tissue incision. MpERF15 operates both upstream and downstream of the MpCOI1 oxylipin receptor by controlling the expression of oxylipin biosynthesis genes. The resulting rise in the oxylipin dinor-12-oxo-phytodienoic acid (dn-OPDA) levels results in an increase in gemma cell number and apical notch organogenesis, generating highly disorganized and compact thalli. Our data pinpoint MpERF15 as a key factor activating an oxylipin biosynthesis amplification loop after wounding, which eventually results in reactivation of cell division and regeneration. We suggest that the genetic networks controlling oxylipin biosynthesis in response to wounding might have been reshuffled over evolution.

Circadian and diel regulation of photosynthesis in the bryophyte Marchantia polymorpha

Circadian rhythms are 24-h biological cycles that align metabolism, physiology, and development with daily environmental fluctuations. Photosynthetic processes are governed by the circadian clock in both flowering plants and some cyanobacteria, but it is unclear how extensively this is conserved throughout the green lineage. We investigated the contribution of circadian regulation to aspects of photosynthesis in Marchantia polymorpha, a liverwort that diverged from flowering plants early in the evolution of land plants. First, we identified in M. polymorpha the circadian regulation of photosynthetic biochemistry, measured using two approaches (delayed fluorescence, pulse amplitude modulation fluorescence). Second, we identified that light-dark cycles synchronize the phase of 24 h cycles of photosynthesis in M. polymorpha, whereas the phases of different thalli desynchronize under free-running conditions. This might also be due to the masking of the underlying circadian rhythms of photosynthesis by light-dark cycles. Finally, we used a pharmacological approach to identify that chloroplast translation might be necessary for clock control of light-harvesting in M. polymorpha. We infer that the circadian regulation of photosynthesis is well-conserved amongst terrestrial plants.

miniTurbo-based interactomics of two plasma membrane-localized SNARE proteins in Marchantia polymorpha

Marchantia polymorpha is a model liverwort and its overall low genetic redundancy is advantageous for dissecting complex pathways. Proximity-dependent in vivo biotin-labelling methods have emerged as powerful interactomics tools in recent years. However, interactomics studies applying proximity labelling are currently limited to angiosperm species in plants. Here, we established and evaluated a miniTurbo-based interactomics method in M. polymorpha using MpSYP12A and MpSYP13B, two plasma membrane-localized SNARE proteins, as baits. We show that our method yields a manifold of potential interactors of MpSYP12A and MpSYP13B compared to a coimmunoprecipitation approach. Our method could capture specific candidates for each SNARE. We conclude that a miniTurbo-based method is a feasible tool for interactomics in M. polymorpha and potentially applicable to other model bryophytes. Our interactome dataset on MpSYP12A and MpSYP13B will be a useful resource to elucidate the evolution of SNARE functions.

Distinct Functions of the Atypical Terminal Hydrophilic Domain of the HKT Transporter in the Liverwort Marchantia polymorpha

K+/Na+ homeostasis is important for land plants, particularly under salt stress. In this study, the structure and ion transport properties of the high-affinity K+ transporter (HKT) of the liverwort Marchantia polymorpha were investigated. Only one HKT gene, MpHKT1, was identified in the genome of M. polymorpha. Phylogenetic analysis of HKT proteins revealed that non-seed plants possess HKTs grouped into a clade independent of the other two clades including HKTs of angiosperms. A distinct long hydrophilic domain was found in the C-terminus of MpHKT1. Complementary DNA (cDNA) of truncated MpHKT1 (t-MpHKT1) encoding the MpHKT_Δ596-812 protein was used to examine the functions of the C-terminal domain. Both MpHKT1 transporters fused with enhanced green fluorescent protein at the N-terminus were localized to the plasma membrane when expressed in rice protoplasts. Two-electrode voltage clamp experiments using Xenopus laevis oocytes indicated that MpHKT1 mediated the transport of monovalent alkali cations with higher selectivity for Na+ and K+, but truncation of the C-terminal domain significantly reduced the transport activity with a decrease in the Na+ permeability. Overexpression of MpHKT1 or t-MpHKT1 in M. polymorpha conferred accumulation of higher Na+ levels and showed higher Na+ uptake rates, compared to those of wild-type plants; however, phenotypes with t-MpHKT1 were consistently weaker than those with MpHKT1. Together, these findings suggest that the hydrophilic C-terminal domain plays a unique role in the regulation of transport activity and ion selectivity of MpHKT1.

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.

Diminished Auxin Signaling Triggers Cellular Reprogramming by Inducing a Regeneration Factor in the Liverwort Marchantia polymorpha

Regeneration in land plants is accompanied by the establishment of new stem cells, which often involves reactivation of the cell division potential in differentiated cells. The phytohormone auxin plays pivotal roles in this process. In bryophytes, regeneration is enhanced by the removal of the apex and repressed by exogenously applied auxin, which has long been proposed as a form of apical dominance. However, the molecular basis behind these observations remains unexplored. Here, we demonstrate that in the liverwort Marchantia polymorpha, the level of endogenous auxin is transiently decreased in the cut surface of decapitated explants, and identify by transcriptome analysis a key transcription factor gene, LOW-AUXIN RESPONSIVE (MpLAXR), which is induced upon auxin reduction. Loss of MpLAXR function resulted in delayed cell cycle reactivation, and transient expression of MpLAXR was sufficient to overcome the inhibition of regeneration by exogenously applied auxin. Furthermore, ectopic expression of MpLAXR caused cell proliferation in normally quiescent tissues. Together, these data indicate that decapitation causes a reduction of auxin level at the cut surface, where, in response, MpLAXR is up-regulated to trigger cellular reprogramming. MpLAXR is an ortholog of Arabidopsis ENHANCER OF SHOOT REGENERATION 1/DORNRÖSCHEN, which has dual functions as a shoot regeneration factor and a regulator of axillary meristem initiation, the latter of which requires a low auxin level. Thus, our findings provide insights into stem cell regulation as well as apical dominance establishment in land plants.

Marchantia polymorpha L. Flavonoids Protect Liver From CCl4-Induced Injury by Antioxidant and Gene-Regulatory Effects

Marchantia polymorpha L. (MPL), a common type of liverwort, has been used as herbal medicine to improve liver function in China for many years. Although modern studies revealed that MPL contains various polyphenols, terpenoids, and bis[bibenzyls], its biological effects on liver function have never been systemically studied in any animal model. In this study, flavonoids were extracted from MPL and the components in the MPL flavonoids as well as the antioxidant capacity of MPL flavonoids were analysed. A rat model of liver injury was induced by intraperitoneal injection of 10% carbon tetrachloride (CCl4). MPL flavonoids were administered daily at a dose of 50, 100, and 200 mg/kg body weight to the rats for 2 weeks prior to injection of CC14. Treatment with MPL flavonoids, especially at a dose of 200 mg/kg, attenuated CCl4-induced increases in alanine aminotransferase, aspartate aminotransferase, alkaline phosphatase, γ-glutamyl transpeptidase, nitric oxide, malondialdehyde, tumour necrosis factor-α, interleukin-1β, and interleukin-6, as well as reductions in superoxide dismutase and glutathione peroxidase. Microarray analyses showed that co-treatment with MPL flavonoids and CCl4 up-regulated many antioxidant and anti-apoptotic genes, but down-regulated several pro-inflammatory genes, compared to treatment with CCl4 alone. PCR and western blot assays further identified that MPL flavonoids increased GPX1, TMX1, TXN, and XIAP expression, but decreased IL-1 and IL1RAP expression and inhibited Jak/stat3 signalling. In conclusion, MPL flavonoids exerted hepatoprotective effects via antioxidant and gene regulatory mechanisms. (Altern Ther Health Med.

An evolutionarily ancient fatty acid desaturase is required for the synthesis of hexadecatrienoic acid, which is the main source of the bioactive jasmonate in Marchantia polymorpha

Jasmonates are fatty acid-derived hormones that regulate multiple aspects of plant development, growth and stress responses. Bioactive jasmonates, defined as the ligands of the conserved COI1 receptor, differ between vascular plants and bryophytes (jasmonoyl-l-isoleucine (JA-Ile) and dinor-12-oxo-10,15(Z)-phytodienoic acid (dn-OPDA), respectively). The biosynthetic pathways of JA-Ile in the model vascular plant Arabidopsis thaliana have been elucidated. However, the details of dn-OPDA biosynthesis in bryophytes are still unclear. Here, we identify an orthologue of Arabidopsis fatty-acid-desaturase 5 (AtFAD5) in the model liverwort Marchantia polymorpha and show that FAD5 function is ancient and conserved between species separated by more than 450 million years (Myr) of independent evolution. Similar to AtFAD5, MpFAD5 is required for the synthesis of 7Z-hexadecenoic acid. Consequently, in Mpfad5 mutants, the hexadecanoid pathway is blocked, dn-OPDA concentrations are almost completely depleted and normal chloroplast development is impaired. Our results demonstrate that the main source of wounding-induced dn-OPDA in Marchantia is the hexadecanoid pathway and the contribution of the octadecanoid pathway (i.e. from OPDA) is minimal. Remarkably, despite extremely low concentrations of dn-OPDA, MpCOI1-mediated responses to wounding and insect feeding can still be activated in Mpfad5, suggesting that dn-OPDA may not be the only bioactive jasmonate and COI1 ligand in Marchantia.

Functions of COP1/SPA E3 Ubiquitin Ligase Mediated by MpCRY in the Liverwort Marchantia polymorpha under Blue Light

COP1/SPA1 complex in Arabidopsis inhibits photomorphogenesis through the ubiquitination of multiple photo-responsive transcription factors in darkness, but such inhibiting function of COP1/SPA1 complex would be suppressed by cryptochromes in blue light. Extensive studies have been conducted on these mechanisms in Arabidopsis whereas little attention has been focused on whether another branch of land plants bryophyte utilizes this blue-light regulatory pathway. To study this problem, we conducted a study in the liverwort Marchantia polymorpha and obtained a MpSPA knock-out mutant, in which Mpspa exhibits the phenotype of an increased percentage of individuals with asymmetrical thallus growth, similar to MpCRY knock-out mutant. We also verified interactions of MpSPA with MpCRY (in a blue light-independent way) and with MpCOP1. Concomitantly, both MpSPA and MpCOP1 could interact with MpHY5, and MpSPA can promote MpCOP1 to ubiquitinate MpHY5 but MpCRY does not regulate the ubiquitination of MpHY5 by MpCOP1/MpSPA complex. These data suggest that COP1/SPA ubiquitinating HY5 is conserved in Marchantia polymorpha, but dissimilar to CRY in Arabidopsis, MpCRY is not an inhibitor of this process under blue light.

DE-ETIOLATED1 has a role in the circadian clock of the liverwort Marchantia polymorpha

Previous studies of plant circadian clock evolution have often relied on clock models and genes defined in Arabidopsis. These studies identified homologues with seemingly conserved function, as well as frequent gene loss. In the present study, we aimed to identify candidate clock genes in the liverwort Marchantia polymorpha using a more unbiased approach. To identify genes with circadian rhythm we sequenced the transcriptomes of gemmalings in a time series in constant light conditions. Subsequently, we performed functional studies using loss-of-function mutants and gene expression reporters. Among the genes displaying circadian rhythm, a homologue to the transcriptional co-repressor Arabidopsis DE-ETIOLATED1 showed high amplitude and morning phase. Because AtDET1 is arrhythmic and associated with the morning gene function of AtCCA1/LHY, that lack a homologue in liverworts, we functionally studied DET1 in M. polymorpha. We found that the circadian rhythm of MpDET1 expression is disrupted in loss-of-function mutants of core clock genes and putative evening-complex genes. MpDET1 knock-down in turn results in altered circadian rhythm of nyctinastic thallus movement and clock gene expression. We could not detect any effect of MpDET1 knock-down on circadian response to light, suggesting that MpDET1 has a yet unknown function in the M. polymorpha circadian clock.

Distinct Functions of Ethylene and ACC in the Basal Land Plant Marchantia polymorpha

Ethylene is a gaseous phytohormone involved in various physiological processes, including fruit ripening, senescence, root hair development and stress responses. Recent genomics studies have suggested that most homologous genes of ethylene biosynthesis and signaling are conserved from algae to angiosperms, whereas the function and biosynthesis of ethylene remain unknown in basal plants. Here, we examined the physiological effects of ethylene, an ethylene precursor, 1-aminocyclopropane-1-carboxylic acid (ACC) and an inhibitor of ethylene perception, silver thiosulfate (STS), in a basal land plant, Marchantia polymorpha. M. polymorpha plants biosynthesized ethylene, and treatment with high concentrations of ACC slightly promoted ethylene production. ACC remarkably suppressed the growth of thalli (vegetative organs) and rhizoids (root-hair-like cells), whereas exogenous ethylene slightly promoted thallus growth. STS suppressed thallus growth and induced ectopic rhizoid formation on the dorsal surface of thalli. Thus, ACC and ethylene have different effects on the vegetative growth of M. polymorpha. We generated single and double mutants of ACC synthase-like (ACSL) genes, MpACSL1 and MpACSL2. The mutants did not show obvious defects in thallus growth, ACC content and ethylene production, indicating that MpACSL genes are not essential for the vegetative growth and biosynthesis of ACC and ethylene. Gene expression analysis suggested the involvement of MpACSL1 and MpACSL2 in stress responses. Collectively, our results imply ethylene-independent function of ACC and the absence of ACC-mediated ethylene biosynthesis in M. polymorpha.

Cloning and functional characterization of three flavonoid O-glucosyltransferase genes from the liverworts Marchantia emarginata and Marchantia paleacea

Flavonoid glycosides are important plant secondary metabolites with broad pharmacological activities. Flavonoid glycosides are generated from aglycones, in reactions catalyzed by typical uridine diphosphate-dependent glycosyltransferases (UGTs). Liverworts produce various types of flavonoid glycosides; however, only two UGTs have been characterized from liverworts to date. Here, we isolated three genes encoding UGTs (MeUGT1, MeUGT2, and MpalUGT1) from the liverwort species Marchantia emarginata and Marchantia paleacea through transcriptome sequencing. Recombinant MeUGT1, MeUGT2, and MpalUGT1 proteins heterologously produced in Escherichia coli exhibited catalytic activity towards multiple flavonoids. MeUGT1 and MpalUGT1 catalyzed the glycosylation of flavonols into the corresponding 3-O-glucosides with UDP-glucose as the sugar donor, while MeUGT2 exhibited a wider substrate specificity that included flavonols, flavones, and flavanones. When MeUGT2 was expressed in E. coli, the yield of flavonol 3-O-glucosides reached to 40-60% with feeding of the substrates kaempferol or quercetin under optimal conditions. Furthermore, heterologous expression of MeUGT1 in Arabidopsis thaliana increased the flavonol glycoside contents in the plants. Therefore, the UGTs characterized in this study could provide new data that will be useful for examining flavonoid biosynthesis in liverworts.

Major components of the KARRIKIN INSENSITIVE2-dependent signaling pathway are conserved in the liverwort Marchantia polymorpha

KARRIKIN INSENSITIVE2 (KAI2) was first identified as a receptor of karrikins, smoke-derived germination stimulants. KAI2 is also considered a receptor of an unidentified endogenous molecule called the KAI2 ligand. Upon KAI2 activation, signals are transmitted through the degradation of D53/SMXL proteins via MAX2-dependent ubiquitination. Although components in the KAI2-dependent signaling pathway, namely MpKAI2A and MpKAI2B, MpMAX2, and MpSMXL, exist in the genome of the liverwort Marchantia polymorpha, their functions remain unknown. Here, we show that early thallus growth is retarded and gemma dormancy in the dark is suppressed in Mpkai2a and Mpmax2 loss-of-function mutants. These defects are counteracted in Mpkai2a Mpsmxl and Mpmax2 Mpsmxl double mutants indicating that MpKAI2A, MpMAX2, and MpSMXL act in the same genetic pathway. Introduction of MpSMXLd53, in which a domain required for degradation is mutated, into wild-type plants mimicks Mpkai2a and Mpmax2 plants. In addition, the detection of citrine fluorescence in Nicotiana benthamiana cells transiently expressing a SMXL-Citrine fusion protein requires treatment with MG132, a proteasome inhibitor. These findings imply that MpSMXL is subjected to degradation, and that the degradation of MpSMXL is crucial for MpKAI2A-dependent signaling in M. polymorpha. Therefore, we claim that the basic mechanisms in the KAI2-dependent signaling pathway are conserved in M. polymorpha.

Fungal-Type Terpene Synthases in Marchantia polymorpha Are Involved in Sesquiterpene Biosynthesis in Oil Body Cells

The liverwort Marchantia polymorpha possesses oil bodies in idioblastic oil body cells scattered in its thallus. Oil bodies are subcellular organelles in which specific sesquiterpenes and bisbibenzyls are accumulated. Therefore, a specialized system for the biosynthesis and accumulation of these defense compounds specifically in oil bodies has been implied. A recent study on M. polymorpha genome sequencing revealed 10 genes that shared high similarities with fungal-type terpene synthases (TPSs). Eight of these fungal-type TPS-like genes in M. polymorpha (MpFTPSL1-6, -9 and -10) are located within a 376-kb stretch on chromosome 6 and share similarities of over 94% at the nucleotide level. Therefore, these genes have likely originated from recent gene duplication events. The expression of a subset of MpFTPSLs was induced under non-axenic growth on vermiculite, which increased the amounts of sesquiterpenes and number of oil bodies. The tdTomato fluorescent protein-based in-fusion reporter assay with MpFTPSL2 promoter revealed fluorescent signals specifically in oil body cells of the thallus, indicating that MpFTPSL2 functions in oil body cells. Recombinant MpFTPSL2 expression in Escherichia coli led to sesquiterpene synthesis from farnesyl pyrophosphate. Moreover, suppression of a subset of MpFTPSLs through RNA interference reduced sesquiterpene accumulation in thalli grown on vermiculite. Taken together, these results suggest that at least a subset of MpFTPSLs is involved in sesquiterpene synthesis in oil body cells.

Mass spectrometry-based quantification and spatial localization of small organic acid exudates in plant roots under phosphorus deficiency and aluminum toxicity

Low-molecular-weight organic acid (OA) extrusion by plant roots is critical for plant nutrition, tolerance to cations toxicity, and plant-microbe interactions. Therefore, methodologies for the rapid and precise quantification of OAs are necessary to be incorporated in the analysis of roots and their exudates. The spatial location of root exudates is also important to understand the molecular mechanisms directing OA production and release into the rhizosphere. Here, we report the development of two complementary methodologies for OA determination, which were employed to evaluate the effect of inorganic ortho-phosphate (Pi) deficiency and aluminum toxicity on OA excretion by Arabidopsis roots. OA exudation by roots is considered a core response to different types of abiotic stress and for the interaction of roots with soil microbes, and for decades has been a target trait to produce plant varieties with increased capacities of Pi uptake and Al tolerance. Using targeted ultra-performance liquid chromatography coupled with high-resolution tandem mass spectrometry (UPLC-HRMS/MS), we achieved the quantification of six OAs in root exudates at sub-micromolar detection limits with an analysis time of less than 5 min per sample. We also employed targeted (MS/MS) matrix-assisted laser desorption/ionization (MALDI) mass spectrometry imaging (MSI) to detect the spatial location of citric and malic acid with high specificity in roots and exudates. Using these methods, we studied OA exudation in response to Al toxicity and Pi deficiency in Arabidopsis seedlings overexpressing genes involved in OA excretion. Finally, we show the transferability of the MALDI-MSI method by analyzing OA excretion in Marchantia polymorpha gemmalings subjected to Pi deficiency.

Expression dynamics of dehydration tolerance in the tropical plant Marchantia inflexa

Tolerance to prolonged water deficit occurs along a continuum in plants, with dehydration tolerance (DhT) and desiccation tolerance (DT) representing some of the most extreme adaptations to water scarcity. Although DhT and DT presumably vary among individuals of a single species, this variability remains largely unstudied. Here, we characterized expression dynamics throughout a dehydration-rehydration time-course in six diverse genotypes of the dioecious liverwort Marchantia inflexa. We identified classical signatures of stress response in M. inflexa, including major changes in transcripts related to metabolism, expression of LEA and ELIP genes, and evidence of cell wall remodeling. However, we detected very little temporal synchronization of these responses across different genotypes of M. inflexa, which may be related to genotypic variation among samples, constitutive expression of dehydration-associated transcripts, the sequestration of mRNAs in ribonucleoprotein partials prior to drying, or the lower tolerance of M. inflexa relative to most bryophytes studied to date. Our characterization of intraspecific variation in expression dynamics suggests that differences in the timing of transcriptional adjustments contribute to variation among genotypes, and that developmental differences impact the relative tolerance of meristematic and differentiated tissues. This work highlights the complexity and variability of water stress tolerance, and underscores the need for comparative studies that seek to characterize variation in DT and DhT.

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.

A START domain-containing protein is involved in the incorporation of ER-derived fatty acids into chloroplast glycolipids in Marchantia polymorpha

The appropriate regulation of thylakoid lipid synthesis is essential for the function of chloroplasts. In plant cells, membrane lipids synthesized in the ER are utilized as a precursor for the synthesis of chloroplast glycolipids. This pathway is thought to be mediated by the transport of glycerolipids synthesized in the ER into chloroplasts. However, we have little knowledge about the proteins involved in the lipid transfer between these organelles in plant cells. Here we show a protein, STAR2, containing the START (Steroidogenic acute regulatory protein-related lipid transfer) domain known to function as a lipid transporter, is involved in the incorporation of ER-derived fatty acids into chloroplast glycolipids in Marchantia polymorpha. We found that STAR2 localizes on the chloroplast envelope membrane as a punctuate structure and is required for the increase of C₂₀ fatty acids, which are synthesized in the ER, in chloroplast glycolipids in response to phosphate deprivation. Our results indicate that STAR2 of M. polymorpha is likely to be involved in the lipid transfer from ER to chloroplast, presumably as a lipid transporter.

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.

Ethylene-independent functions of the ethylene precursor ACC in Marchantia polymorpha

The plant hormone ethylene has many roles in growth and development¹. In seed plants, the ethylene precursor 1-aminocyclopropane-1-carboxylic acid (ACC) is converted into ethylene by ACC oxidase (ACO), and treatment with ACC induces ethylene responses². However, non-seed plants lack ACO homologues^3-8, which led us to examine the relationship between ACC and ethylene in the liverwort Marchantia polymorpha. Here, we demonstrate that ACC and ethylene can induce divergent growth responses in Marchantia. Ethylene increases plant and gemma size, induces more gemma cups and promotes gemmae dormancy. As predicted, Mpctr1-knockout mutants display constitutive ethylene responses, whereas Mpein3-knockout mutants exhibit ethylene insensitivity. Compared with the wild type, Mpctr1 gemmae have more and larger epidermal cells, whereas Mpein3 gemmae have fewer and smaller epidermal cells, suggesting that ethylene promotes cell division and growth in developing gemmae. By contrast, ACC treatment inhibits gemma growth and development by suppressing cell division, even in the Mpein3-knockout alleles. Knockout mutants of one or both ACC SYNTHASE (ACS) gene homologues produce negligible levels of ACC, have more and larger gemma cups, and have more-expanded thallus branches. Mpacs2 and Mpacs1 Mpacs2 gemmae also display a high frequency of abnormal apical notches (meristems) that are not observed in ethylene mutants. These findings reveal that ethylene and ACC have distinct functions, and suggest that ACC is a signalling molecule in Marchantia. ACC may be an evolutionarily conserved signal that predates its efficient conversion to ethylene in higher plants.

Marchantia TCP transcription factor activity correlates with three-dimensional chromatin structure

Information in the genome is not only encoded within sequence or epigenetic modifications, but is also found in how it folds in three-dimensional space. The formation of self-interacting genomic regions, named topologically associated domains (TADs), is known as a key feature of genome organization beyond the nucleosomal level. However, our understanding of the formation and function of TADs in plants is extremely limited. Here we show that the genome of Marchantia polymorpha, a member of a basal land plant lineage, exhibits TADs with epigenetic features similar to those of higher plants. By analysing various epigenetic marks across Marchantia TADs, we find that these regions generally represent interstitial heterochromatin and their borders are enriched with Marchantia transcription factor TCP1. We also identify a type of TAD that we name 'TCP1-rich TAD', in which genomic regions are highly accessible and are densely bound by TCP1 proteins. Transcription of TCP1 target genes differs on the basis gene location, and those in TCP1-rich TADs clearly show a lower expression level. In tcp1 mutant lines, neither TCP1-bound TAD borders nor TCP1-rich TADs display drastically altered chromatin organization patterns, suggesting that, in Marchantia, TCP1 is dispensable for TAD formation. However, we find that in tcp1 mutants, genes residing in TCP1-rich TADs have a greater extent of expression fold change as opposed to genes that do not belong to these TADs. Our results suggest that, besides standing as spatial chromatin-packing modules, plant TADs function as nuclear microcompartments associated with transcription factor activities.

The cold-induced switch in direction of chloroplast relocation occurs independently of changes in endogenous phototropin levels

When exposed to fluctuating light intensity, chloroplasts move towards weak light (accumulation response), and away from strong light (avoidance response). In addition, cold treatment (5°C) induces the avoidance response even under weak-light conditions (cold-avoidance response). These three responses are mediated by the phototropin (phot), which is a blue-light photoreceptor and has also been reported to act as a thermosensory protein that perceives temperature variation. Our previous report indicated that cold-induced changes in phot biochemical activity initiate the cold-avoidance response. In this study, we further explored the induction mechanism of the cold-avoidance response in the liverwort Marchantia polymorpha and examined the relationship between changes in the amount of phot and the induction of the cold-avoidance response. The switch between the accumulation and avoidance responses occurs at a so-called ‘transitional’ light intensity. Our physiological experiments revealed that a cold-mediated decrease in the transitional light intensity leads to the induction of the cold-avoidance response. While artificial overexpression of phot decreased the transitional light intensity as much as cold treatment did, the amount of endogenous phot was not increased by cold treatment in wild-type M. polymorpha. Taken together, these findings show that the cold-avoidance response is initiated by a cold-mediated reduction of the transitional light intensity, independent of the amount of endogenous phot. This study provides a clue to understanding the mechanism underlying the switch in direction of chloroplast relocation in response to light and temperature.

Design principles of a minimal auxin response system

Auxin controls numerous growth processes in land plants through a gene expression system that modulates ARF transcription factor activity^1-3. Gene duplications in families encoding auxin response components have generated tremendous complexity in most land plants, and neofunctionalization enabled various unique response outputs during development^1,3,4. However, it is unclear what fundamental biochemical principles underlie this complex response system. By studying the minimal system in Marchantia polymorpha, we derive an intuitive and simple model where a single auxin-dependent A-ARF activates gene expression. It is antagonized by an auxin-independent B-ARF that represses common target genes. The expression patterns of both ARF proteins define developmental zones where auxin response is permitted, quantitatively tuned or prevented. This fundamental design probably represents the ancestral system and formed the basis for inflated, complex systems.

Regulation of the Poly(A) Status of Mitochondrial mRNA by Poly(A)-Specific Ribonuclease Is Conserved among Land Plants

Regulation of the stability and the quality of mitochondrial RNA is essential for the maintenance of mitochondrial and cellular functions in eukaryotes. We have previously reported that the eukaryotic poly(A)-specific ribonuclease (PARN) and the prokaryotic poly(A) polymerase encoded by AHG2 and AGS1, respectively, coordinately regulate the poly(A) status and the stability of mitochondrial mRNA in Arabidopsis. Mitochondrial function of PARN has not been reported in any other eukaryotes. To know how much this PARN-based mitochondrial mRNA regulation is conserved among plants, we studied the AHG2 and AGS1 counterparts of the liverwort, Marchantia polymorpha, a member of basal land plant lineage. We found that M. polymorpha has one ortholog each for AHG2 and AGS1, named MpAHG2 and MpAGS1, respectively. Their Citrine-fused proteins were detected in mitochondria of the liverwort. Molecular genetic analysis showed that MpAHG2 is essential and functionally interacts with MpAGS1 as observed in Arabidopsis. A recombinant MpAHG2 protein had a deadenylase activity in vitro. Overexpression of MpAGS1 and the reduced expression of MpAHG2 caused an accumulation of polyadenylated Mpcox1 mRNA. Furthermore, MpAHG2 suppressed Arabidopsis ahg2-1 mutant phenotype. These results suggest that the PARN-based mitochondrial mRNA regulatory system is conserved in land plants.

Regulation of Photosynthetic Carbohydrate Metabolism by a Raf-Like Kinase in the Liverwort Marchantia polymorpha

To optimize growth and development, plants monitor photosynthetic activities and appropriately regulate various cellular processes. However, signaling mechanisms that coordinate plant growth with photosynthesis remain poorly understood. To identify factors that are involved in signaling related to photosynthetic stimuli, we performed a phosphoproteomic analysis with Marchantia polymorpha, an extant bryophyte species in the basal lineage of land plants. Among proteins whose phosphorylation status changed differentially between dark-treated plants and those after light irradiation but failed to do so in the presence of a photosynthesis inhibitor, we identified a B4-group Raf-like kinase, named PHOTOSYNTHESIS-RELATED RAF (MpPRAF). Biochemical analyses confirmed photosynthesis-activity-dependent changes in the phosphorylation status of MpPRAF. Mutations in the MpPRAF gene resulted in growth retardation. Measurement of carbohydrates demonstrated both hyper-accumulation of starch and reduction of sucrose in Mppraf mutants. Neither inhibition of starch synthesis nor exogenous supply of sucrose alleviated the growth defect, suggesting serious impairment of Mppraf mutants in both the synthesis of sucrose and the repression of its catabolism. As a result of the compromised photosynthate metabolism, photosynthetic electron transport was downregulated in Mppraf mutants. A mutated MpPRAF with a common amino acid substitution for inactivating kinase activity was unable to rescue the Mppraf mutant defects. Our results provide evidence that MpPRAF is a photosynthesis signaling kinase that regulates sucrose metabolism.

MpTCP1 controls cell proliferation and redox processes in Marchantia polymorpha

TCP transcription factors are key regulators of angiosperm cell proliferation processes. It is unknown whether their regulatory growth capacities are conserved across land plants, which we examined in liverworts, one of the earliest diverging land plant lineages. We generated knockout mutants for MpTCP1, the single TCP-P clade gene in Marchantia polymorpha, and characterized its function by conducting cell proliferation and morphological analyses as well as messenger RNA expression, transcriptome, chemical, and DNA binding studies. Mptcp1^ge lines show a reduced vegetative thallus growth and extra tissue formation in female reproductive structures. Additionally, mutant plants reveal increased hydrogen peroxide (H₂ O₂ ) levels and an enhanced pigmentation in the thallus caused by formation of secondary metabolites, such as aminochromes. MpTCP1 proteins interact redox dependently with DNA and regulate the expression of a comprehensive redox network, comprising enzymes involved in H₂ O₂ metabolism. MpTCP1 regulates Marchantia growth in a context-dependent manner. Redox sensitivity of the DNA binding capacity of MpTCP1 proteins provides a mechanism to respond to altered redox conditions. Our data suggest that MpTCP1 activity could thereby have contributed to diversification of land plant morphologies and to adaptations to abiotic and biotic challenges, as experienced by liverworts during early land plant colonization.

Jasmonate-Related MYC Transcription Factors Are Functionally Conserved in Marchantia polymorpha

The lipid-derived phytohormone jasmonoyl-isoleucine regulates plant immunity, growth and development in vascular plants by activating genome-wide transcriptional reprogramming. In Arabidopsis (Arabidopsis thaliana), this process is largely orchestrated by the master regulator MYC2 and related transcription factors (TFs). However, the TFs activating this pathway in basal plant lineages are currently unknown. We report the functional conservation of MYC-related TFs between the eudicot Arabidopsis and the liverwort Marchantia polymorpha, a plant belonging to an early diverging lineage of land plants. Phylogenetic analysis suggests that MYC function first appeared in charophycean algae and therefore predates the evolutionary appearance of any other jasmonate pathway component. M. polymorpha possesses two functionally interchangeable MYC genes, one in females and one in males. Similar to AtMYC2, MpMYCs showed nuclear localization, interaction with JASMONATE-ZIM-DOMAIN PROTEIN repressors, and regulation by light. Phenotypic and molecular characterization of loss- and gain-of-function mutants demonstrated that MpMYCs are necessary and sufficient for activating the jasmonate pathway in M. polymorpha, but unlike their Arabidopsis orthologs, do not regulate fertility. Therefore, despite 450 million years of independent evolution, MYCs are functionally conserved between bryophytes and eudicots. Genetic conservation in an early diverging lineage suggests that MYC function existed in the common ancestor of land plants and evolved from a preexisting MYC function in charophycean algae.

Cytokinin signaling coordinates development of diverse organs in Marchantia polymorpha

Cytokinins play an essential role in plant growth and development. A recent study showed that the cytokinin signaling pathway was conserved in the liverwort Marchantia polymorpha, and that it controlled gemma cup and rhizoid formation during thallus development. Here we show that the type-B response regulator, MpRRB, is mainly localized in the nucleus. Moreover, observations of thalli revealed that the distribution of air pores and the shape of the thallus margin are impaired in cytokinin-deficient lines and those defective in cytokinin signaling. This suggests that cytokinins regulate cell division and/or differentiation of precursor cells derived from the apical cell, thereby coordinating development of various organs produced on the thallus.