A non-canonical BZR/BES transcription factor regulates the development of haploid reproductive organs in Marchantia polymorpha.

Gametogenesis, which is essential to the sexual reproductive system, has drastically changed during plant evolution. Bryophytes, lycophytes and ferns develop reproductive organs called gametangia-antheridia and archegonia for sperm and egg production, respectively. However, the molecular mechanism of early gametangium development remains unclear. Here we identified a 'non-canonical' type of BZR/BES transcription factor, MpBZR3, as a regulator of gametangium development in a model bryophyte, Marchantia polymorpha. Interestingly, overexpression of MpBZR3 induced ectopic gametangia. Genetic analysis revealed that MpBZR3 promotes the early phase of antheridium development in male plants. By contrast, MpBZR3 is required for the late phase of archegonium development in female plants. We demonstrate that MpBZR3 is necessary for the successful development of both antheridia and archegonia but functions in a different manner between the two sexes. Together, the functional specialization of this 'non-canonical' type of BZR/BES member may have contributed to the evolution of reproductive systems.

The cAMP signaling module regulates sperm motility in the liverwort Marchantia polymorpha.

Adenosine 3',5'-cyclic monophosphate (cAMP) is a universal signaling molecule that acts as a second messenger in various organisms. It is well established that cAMP plays essential roles across the tree of life, although the function of cAMP in land plants has long been debated. We previously identified the enzyme with both adenylyl cyclase (AC) and cAMP phosphodiesterase (PDE) activity as the cAMP-synthesis/hydrolysis enzyme COMBINED AC with PDE (CAPE) in the liverwort Marchantia polymorpha. CAPE is conserved in streptophytes that reproduce with motile sperm; however, the precise function of CAPE is not yet known. In this study, we demonstrate that the loss of function of CAPE in M. polymorpha led to male infertility due to impaired sperm flagellar motility. We also found that two genes encoding the regulatory subunits of cAMP-dependent protein kinase (PKA-R) were also involved in sperm motility. Based on these findings, it is evident that CAPE and PKA-Rs act as a cAMP signaling module that regulates sperm motility in M. polymorpha. Therefore, our results have shed light on the function of cAMP signaling and sperm motility regulators in land plants. This study suggests that cAMP signaling plays a common role in plant and animal sperm motility.

Dual regulation of cytochrome P450 gene expression by two distinct small RNAs, a novel tasiRNA and miRNA, in Marchantia polymorpha.

The miR390-derived TAS3 trans-acting short-interfering RNAs (tasiRNAs) module represents a conserved RNA silencing pathway in the plant kingdom; however, its characterization in the bryophyte Marchantia polymorpha is limited. This study elucidated that MpDCL4 processes MpTAS3 double-stranded RNA (dsRNA) to generate tasiRNAs, primarily from the 5'- and 3'-ends of dsRNA. Notably, we discovered a novel tasiRNA, tasi78A, can negatively regulate a cytochrome P450 gene, MpCYP78A101. Additionally, tasi78A was abundant in MpAGO1, and transient expression assays underscored the role of tasi78A in repressing MpCYP78A101. A microRNA, miR11700, also regulates MpCYP78A101 expression. This coordinate regulation suggests a role in modulating auxin signaling at apical notches of gemma, influencing the growth and sexual organ development of M. polymorpha and emphasizing the significance of RNA silencing in MpCYP78A101 regulation. However, phylogenetic analysis identified another paralog of the CYP78 family, Mp1g14150, which may have a redundant role with MpCYP78A101, explaining the absence of noticeable morphological changes in loss-of-function plants. Taken together, our findings provide new insights into the combined regulatory roles of miR390/MpTAS3/miR11700 in controlling MpCYP78A101 and expand our knowledge about the biogenesis and regulation of tasiRNAs in M. polymorpha.

Light-induced stomatal opening requires phosphorylation of the C-terminal autoinhibitory domain of plasma membrane H+-ATPase.

Plasma membrane H+-ATPase provides the driving force for light-induced stomatal opening. However, the mechanisms underlying the regulation of its activity remain unclear. Here, we show that the phosphorylation of two Thr residues in the C-terminal autoinhibitory domain is crucial for H+-ATPase activation and stomatal opening in Arabidopsis thaliana. Using phosphoproteome analysis, we show that blue light induces the phosphorylation of Thr-881 within the C-terminal region I, in addition to penultimate Thr-948 in AUTOINHIBITED H+-ATPASE 1 (AHA1). Based on site-directed mutagenesis experiments, phosphorylation of both Thr residues is essential for H+ pumping and stomatal opening in response to blue light. Thr-948 phosphorylation is a prerequisite for Thr-881 phosphorylation by blue light. Additionally, red light-driven guard cell photosynthesis induces Thr-881 phosphorylation, possibly contributing to red light-dependent stomatal opening. Our findings provide mechanistic insights into H+-ATPase activation that exploits the ion transport across the plasma membrane and light signalling network in guard cells.

Convergent evolution of water-conducting cells in Marchantia recruited the ZHOUPI gene promoting cell wall reinforcement and programmed cell death.

A key adaptation of plants to life on land is the formation of water-conducting cells (WCCs) for efficient long-distance water transport. Based on morphological analyses it is thought that WCCs have evolved independently on multiple occasions. For example, WCCs have been lost in all but a few lineages of bryophytes but, strikingly, within the liverworts a derived group, the complex thalloids, has evolved a novel externalized water-conducting tissue composed of reinforced, hollow cells termed pegged rhizoids. Here, we show that pegged rhizoid differentiation in Marchantia polymorpha is controlled by orthologs of the ZHOUPI and ICE bHLH transcription factors required for endosperm cell death in Arabidopsis seeds. By contrast, pegged rhizoid development was not affected by disruption of MpNAC5, the Marchantia ortholog of the VND genes that control WCC formation in flowering plants. We characterize the rapid, genetically controlled programmed cell death process that pegged rhizoids undergo to terminate cellular differentiation and identify a corresponding upregulation of conserved putative plant cell death effector genes. Lastly, we show that ectopic expression of MpZOU1 increases production of pegged rhizoids and enhances drought tolerance. Our results support that pegged rhizoids evolved independently of other WCCs. We suggest that elements of the genetic control of developmental cell death are conserved throughout land plants and that the ZHOUPI/ICE regulatory module has been independently recruited to promote cell wall modification and programmed cell death in liverwort rhizoids and in the endosperm of flowering plant seed.

Transcription of the Antisense Long Non-Coding RNA, SUPPRESSOR OF FEMINIZATION, Represses Expression of the Female-Promoting Gene FEMALE GAMETOPHYTE MYB in the Liverwort Marchantia polymorpha.

Sexual differentiation is a fundamental process in the life cycles of land plants, ensuring successful sexual reproduction and thereby contributing to species diversity and survival. In the dioicous liverwort Marchantia polymorpha, this process is governed by an autosomal sex-differentiation locus comprising FEMALE GAMETOPHYTE MYB (FGMYB), a female-promoting gene, and SUPPRESSOR OF FEMINIZATION (SUF), an antisense strand-encoded long non-coding RNA (lncRNA). SUF is specifically transcribed in male plants and suppresses the expression of FGMYB, leading to male differentiation. However, the molecular mechanisms underlying this process remain elusive. Here, we show that SUF acts through its transcription to suppress FGMYB expression. Transgene complementation analysis using CRISPR/Cas9D10A-based large-deletion mutants identified a genomic region sufficient for the sex differentiation switch function in the FGMYB-SUF locus. Inserting a transcriptional terminator sequence into the SUF-transcribed region resulted in the loss of SUF function and allowed expression of FGMYB in genetically male plants, leading to conversion of the sex phenotype from male to female. Partial deletions of SUF had no obvious impact on its function. Replacement of the FGMYB sequence with that of an unrelated gene did not affect the ability of SUF transcription to suppress sense-strand expression. Taken together, our findings suggest that the process of SUF transcription, rather than the resulting transcripts, is required for controlling sex differentiation in M. polymorpha.

Thermospermine Is an Evolutionarily Ancestral Phytohormone Required for Organ Development and Stress Responses in Marchantia Polymorpha.

Thermospermine suppresses auxin-inducible xylem differentiation, whereas its structural isomer, spermine, is involved in stress responses in angiosperms. The thermospermine synthase, ACAULIS5 (ACL5), is conserved from algae to land plants, but its physiological functions remain elusive in non-vascular plants. Here, we focused on MpACL5, a gene in the liverwort Marchantia polymorpha, that rescued the dwarf phenotype of the acl5 mutant in Arabidopsis. In the Mpacl5 mutants generated by genome editing, severe growth retardation was observed in the vegetative organ, thallus, and the sexual reproductive organ, gametangiophore. The mutant gametangiophores exhibited remarkable morphological defects such as short stalks, fasciation and indeterminate growth. Two gametangiophores fused together, and new gametangiophores were often initiated from the old ones. Furthermore, Mpacl5 showed altered responses to heat and salt stresses. Given the absence of spermine in bryophytes, these results suggest that thermospermine has a dual primordial function in organ development and stress responses in M. polymorpha. The stress response function may have eventually been assigned to spermine during land plant evolution.

Conserved CKI1-mediated signaling is required for female germline specification in Marchantia polymorpha.

In land plants, gametes derive from a small number of dedicated haploid cells.1 In angiosperms, one central cell and one egg cell are differentiated in the embryo sac as female gametes for double fertilization, while in non-flowering plants, only one egg cell is generated in the female sexual organ, called the archegonium.2,3 The central cell specification of Arabidopsis thaliana is controlled by the histidine kinase CYTOKININ-INDEPENDENT 1 (CKI1), which is a two-component signaling (TCS) activator sharing downstream regulatory components with the cytokinin signaling pathway.4,5,6,7 Our phylogenetic analysis suggested that CKI1 orthologs broadly exist in land plants. However, the role of CKI1 in non-flowering plants remains unclear. Here, we found that the sole CKI1 ortholog in the liverwort Marchantia polymorpha, MpCKI1, which functions through conserved downstream TCS components, regulates the female germline specification for egg cell development in the archegonium. In M. polymorpha, the archegonium develops three-dimensionally from a single cell accumulating MpBONOBO (MpBNB), a master regulator for germline initiation and differentiation.8 We visualized female germline specification by capturing the distribution pattern of MpBNB in discrete stages of early archegonium development, and found that MpBNB accumulation is restricted to female germline cells. MpCKI1 is required for the proper MpBNB accumulation in the female germline, and is critical for the asymmetric cell divisions that specify the female germline cells. These results suggest that CKI1-mediated TCS originated during early land plant evolution and participates in female germ cell specification in deeply diverged plant lineages.

The phosphorylated pathway of serine biosynthesis affects sperm, embryo, and sporophyte development, and metabolism in Marchantia polymorpha.

Serine metabolism is involved in various biological processes. Here we investigate primary functions of the phosphorylated pathway of serine biosynthesis in a non-vascular plant Marchantia polymorpha by analyzing knockout mutants of MpPGDH encoding 3-phosphoglycerate dehydrogenase in this pathway. Growth phenotypes indicate that serine from the phosphorylated pathway in the dark is crucial for thallus growth. Sperm development requires serine from the phosphorylated pathway, while egg formation does not. Functional MpPGDH in the maternal genome is necessary for embryo and sporophyte development. Under high CO2 where the glycolate pathway of serine biosynthesis is inhibited, suppressed thallus growth of the mutants is not fully recovered by exogenously-supplemented serine, suggesting the importance of serine homeostasis involving the phosphorylated and glycolate pathways. Metabolomic phenotypes indicate that the phosphorylated pathway mainly influences the tricarboxylic acid cycle, the amino acid and nucleotide metabolism, and lipid metabolism. These results indicate the importance of the phosphorylated pathway of serine biosynthesis in the dark, in the development of sperm, embryo, and sporophyte, and metabolism in M. polymorpha.

RAF-like protein kinases mediate a deeply conserved, rapid auxin response.

The plant-signaling molecule auxin triggers fast and slow cellular responses across land plants and algae. The nuclear auxin pathway mediates gene expression and controls growth and development in land plants, but this pathway is absent from algal sister groups. Several components of rapid responses have been identified in Arabidopsis, but it is unknown if these are part of a conserved mechanism. We recently identified a fast, proteome-wide phosphorylation response to auxin. Here, we show that this response occurs across 5 land plant and algal species and converges on a core group of shared targets. We found conserved rapid physiological responses to auxin in the same species and identified rapidly accelerated fibrosarcoma (RAF)-like protein kinases as central mediators of auxin-triggered phosphorylation across species. Genetic analysis connects this kinase to both auxin-triggered protein phosphorylation and rapid cellular response, thus identifying an ancient mechanism for fast auxin responses in the green lineage.

Harnessing Deep Learning to Analyze Cryptic Morphological Variability of Marchantia polymorpha.

Characterizing phenotypes is a fundamental aspect of biological sciences, although it can be challenging due to various factors. For instance, the liverwort Marchantia polymorpha is a model system for plant biology and exhibits morphological variability, making it difficult to identify and quantify distinct phenotypic features using objective measures. To address this issue, we utilized a deep-learning-based image classifier that can handle plant images directly without manual extraction of phenotypic features and analyzed pictures of M. polymorpha. This dioicous plant species exhibits morphological differences between male and female wild accessions at an early stage of gemmaling growth, although it remains elusive whether the differences are attributable to sex chromosomes. To isolate the effects of sex chromosomes from autosomal polymorphisms, we established a male and female set of recombinant inbred lines (RILs) from a set of male and female wild accessions. We then trained deep learning models to classify the sexes of the RILs and the wild accessions. Our results showed that the trained classifiers accurately classified male and female gemmalings of wild accessions in the first week of growth, confirming the intuition of researchers in a reproducible and objective manner. In contrast, the RILs were less distinguishable, indicating that the differences between the parental wild accessions arose from autosomal variations. Furthermore, we validated our trained models by an 'eXplainable AI' technique that highlights image regions relevant to the classification. Our findings demonstrate that the classifier-based approach provides a powerful tool for analyzing plant species that lack standardized phenotyping metrics.

A bHLH heterodimer regulates germ cell differentiation in land plant gametophytes.

Land plants are a monophyletic group of photosynthetic eukaryotes that diverged from streptophyte algae about 470 million years ago. During both the alternating haploid and diploid stages of the life cycle, land plants form multicellular bodies.1,2,3,4 The haploid multicellular body (gametophyte) produces progenitor cells that give rise to gametes and the reproductive organs.5,6,7,8 In the liverwort Marchantia polymorpha, differentiation of the initial cells of gamete-producing organs (gametangia) from the gametophyte is regulated by MpBONOBO (MpBNB), a member of the basic helix-loop-helix (bHLH) transcription factor subfamily VIIIa. In Arabidopsis thaliana, specification of generative cells in developing male gametophytes (pollen) requires redundant action of BNB1 and BNB2.9 Subfamily XI bHLHs, such as LOTUS JAPONICUS ROOTHAIRLESS LIKE1 (LRL1)/DEFECTIVE REGION OF POLLEN1 (DROP1) and LRL2/DROP2 in A. thaliana and the single LRL/DROP protein MpLRL in M. polymorpha, are the evolutionarily conserved regulators of rooting system development.10 Although the role of LRL1/DROP1 and LRL2/DROP2 in gametogenesis remains unclear, their loss leads to the formation of abnormal pollen devoid of sperm cells.11 Here, we show that BNBs and LRL/DROPs co-localize to gametophytic cell nuclei and form heterodimers. LRL1/DROP1 and LRL2/DROP2 act redundantly to regulate BNB expression for generative cell specification in A. thaliana after asymmetric division of the haploid microspore. MpLRL is required for differentiation of MpBNB-expressing gametangium initial cells in M. polymorpha gametophytes. Our findings suggest that broadly expressed LRL/DROP stabilizes BNB expression, leading to the formation of an evolutionarily conserved bHLH heterodimer, which regulates germ cell differentiation in the haploid gametophyte of land plants.

LysM-mediated signaling in Marchantia polymorpha highlights the conservation of pattern-triggered immunity in land plants.

Pattern-recognition receptor (PRR)-triggered immunity (PTI) wards off a wide range of pathogenic microbes, playing a pivotal role in angiosperms. The model liverwort Marchantia polymorpha triggers defense-related gene expression upon sensing components of bacterial and fungal extracts, suggesting the existence of PTI in this plant model. However, the molecular components of the putative PTI in M. polymorpha and the significance of PTI in bryophytes have not yet been described. We here show that M. polymorpha has four lysin motif (LysM)-domain-containing receptor homologs, two of which, LysM-receptor-like kinase (LYK) MpLYK1 and LYK-related (LYR) MpLYR, are responsible for sensing chitin and peptidoglycan fragments, triggering a series of characteristic immune responses. Comprehensive phosphoproteomic analysis of M. polymorpha in response to chitin treatment identified regulatory proteins that potentially shape LysM-mediated PTI. The identified proteins included homologs of well-described PTI components in angiosperms as well as proteins whose roles in PTI are not yet determined, including the blue-light receptor phototropin MpPHOT. We revealed that MpPHOT is required for negative feedback of defense-related gene expression during PTI. Taken together, this study outlines the basic framework of LysM-mediated PTI in M. polymorpha and highlights conserved elements and new aspects of pattern-triggered immunity in land plants.

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.

Three-Dimensional Morphological Analysis Revealed the Cell Patterning Bases for the Sexual Dimorphism Development in the Liverwort Marchantia polymorpha.

In land plants, sexual dimorphism can develop in both diploid sporophytes and haploid gametophytes. While developmental processes of sexual dimorphism have been extensively studied in the sporophytic reproductive organs of model flowering plants such as stamens and carpels of Arabidopsis thaliana, those occurring in gametophyte generation are less well characterized due to the lack of amenable model systems. In this study, we performed three-dimensional morphological analyses of gametophytic sexual branch differentiation in the liverwort Marchantia polymorpha, using high-depth confocal imaging and a computational cell segmentation technique. Our analysis revealed that the specification of germline precursors initiates in a very early stage of sexual branch development, where incipient branch primordia are barely recognizable in the apical notch region. Moreover, spatial distribution patterns of germline precursors differ between males and females from the initial stage of primordium development in a manner dependent on the master sexual differentiation regulator MpFGMYB. At later stages, distribution patterns of germline precursors predict the sex-specific gametangia arrangement and receptacle morphologies seen in mature sexual branches. Taken together, our data suggest a tightly coupled progression of germline segregation and sexual dimorphism development in M. polymorpha.

Stomatal regulators are co-opted for seta development in the astomatous liverwort Marchantia polymorpha.

The evolution of special types of cells requires the acquisition of new gene regulatory networks controlled by transcription factors (TFs). In stomatous plants, a TF module formed by subfamilies Ia and IIIb basic helix-loop-helix TFs (Ia-IIIb bHLH) regulates stomatal formation; however, how this module evolved during land plant diversification remains unclear. Here we show that, in the astomatous liverwort Marchantia polymorpha, a Ia-IIIb bHLH module regulates the development of a unique sporophyte tissue, the seta, which is found in mosses and liverworts. The sole Ia bHLH gene, MpSETA, and a IIIb bHLH gene, MpICE2, regulate the cell division and/or differentiation of seta lineage cells. MpSETA can partially replace the stomatal function of Ia bHLH TFs in Arabidopsis thaliana, suggesting that a common regulatory mechanism underlies setal and stomatal formation. Our findings reveal the co-option of a Ia-IIIb bHLH TF module for regulating cell fate determination and/or cell division of distinct types of cells during land plant evolution.

Auxin signaling is essential for organogenesis but not for cell survival in the liverwort Marchantia polymorpha.

Auxin plays pleiotropic roles in plant development via gene regulation upon its perception by the receptors TRANSPORT INHIBITOR RESPONSE 1/AUXIN SIGNALING F-BOX (TIR1/AFBs). This auxin-regulated transcriptional control mechanism originated in the common ancestor of land plants. Although the complete loss of TIR1/AFBs causes embryonic lethality in Arabidopsis thaliana, it is unclear whether the requirement for TIR1-mediated auxin perception in cell viability can be generalized. The model liverwort Marchantia polymorpha has a minimal auxin signaling system with only a single TIR1/AFB, MpTIR1. Here we show by genetic, biochemical, and transcriptomic analyses that MpTIR1 functions as an evolutionarily conserved auxin receptor. Null mutants and conditionally knocked-out mutants of MpTIR1 were viable but incapable of forming any organs and grew as cell masses. Principal component analysis performed using transcriptomes at various developmental stages indicated that MpTIR1 is involved in the developmental transition from spores to organized thalli, during which apical notches containing stem cells are established. In Mptir1 cells, stem cell- and differentiation-related genes were up- and downregulated, respectively. Our findings suggest that, in M. polymorpha, auxin signaling is dispensable for cell division but is essential for three-dimensional patterning of the plant body by establishing pluripotent stem cells for organogenesis, a derived trait of land plants.

R2R3-MYB transcription factor GEMMA CUP-ASSOCIATED MYB1 mediates the cytokinin signal to achieve proper organ development in Marchantia polymorpha.

Cytokinin, a plant hormone, plays essential roles in organ growth and development. The type-B response regulator-mediated cytokinin signaling is repressed by type-A response regulators and is conserved in the liverwort Marchantia polymorpha. Its signal coordinates the development of diverse organs on the thallus body, such as the gemma cup, rhizoid, and air pores. Here we report that the type-B response regulator MpRRB upregulates the expression of the R2R3-MYB transcription factor GEMMA CUP-ASSOCIATED MYB1 (MpGCAM1) in M. polymorpha. Whereas both Mpgcam1 and Mprrb knockout mutants exhibited defects in gemma cup formation, the Mpgcam1 Mprra double mutant, in which cytokinin signaling is activated due to the lack of type-A response regulator, also formed no gemma cups. This suggests that MpGCAM1 functions downstream of cytokinin signaling. Inducible overexpression of MpGCAM1 produced undifferentiated cell clumps on the thalli of both wild-type and Mprrb. However, smaller thalli were formed in Mprrb compared to the wild-type after the cessation of overexpression. These results suggest that cytokinin signaling promotes gemma cup formation and cellular reprogramming through MpGCAM1, while cytokinin signals also participate in activating cell division during thallus development.

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.