Efficient CRISPR/Cas9-based genome editing and its application to conditional genetic analysis in Marchantia polymorpha

Regulation of ROP GTPase cycling between active and inactive states is essential for vegetative organogenesis in Marchantia polymorpha

Rho/Rac of plant (ROP) GTPases are plant-specific proteins that function as molecular switches, activated by guanine nucleotide exchange factors (GEFs) and inactivated by GTPase-activating proteins (GAPs). The bryophyte Marchantia polymorpha contains single copies of ROP (MpROP), GEFs [ROPGEF and SPIKE (SPK)] and GAPs [ROPGAP and ROP ENHANCER (REN)]. MpROP regulates the development of various tissues and organs, such as rhizoids, gemmae and air chambers. The ROPGEF KARAPPO (MpKAR) is essential for gemma initiation, but the functions of other ROP regulatory factors are less understood. This study focused on two GAPs: MpROPGAP and MpREN. Mpren single mutants showed defects in thallus growth, rhizoid tip growth, gemma development, and air-chamber formation, whereas Mpropgap mutants showed no visible abnormalities. However, Mpropgap Mpren double mutants had more severe phenotypes than the Mpren single mutants, suggesting backup roles of MpROPGAP in processes involving MpREN. Overexpression of MpROPGAP and MpREN resulted in similar gametophyte defects, highlighting the importance of MpROP activation/inactivation cycling (or balancing). Thus, MpREN predominantly, and MpROPGAP as a backup, regulate gametophyte development, likely by controlling MpROP activation in M. polymorpha.

Molecular Insights into MpAGO1 and Its Regulatory miRNA, miR11707, in the High-Temperature Acclimation of Marchantia polymorpha

As a model plant for bryophytes, Marchantia polymorpha offers insights into the role of RNA silencing in aiding early land plants navigate the challenges posed by high-temperature environments. Genomic analysis revealed unique ARGONAUTE1 ortholog gene (MpAGO1) in M. polymorpha, which is regulated by two species-specific microRNAs (miRNAs), miR11707.1 and miR11707.2. Comparative studies of small RNA profiles from M. polymorpha cellular and MpAGO1 immunoprecipitation (MpAGO1-IP) profiles at various temperatures, along with analyses of Arabidopsis AGO1 (AtAGO1), revealed that MpAGO1 has a low selectivity for a diverse range of small RNA species than AtAGO1. Protein structural comparisons revealed no discernible differences in the guide strand small RNA recognition middle domain, MID domain, of MpAGO1 and AtAGO1, suggesting the complexity of miRNA species specificity and necessitating further exploration. Small RNA profiling and size exclusion chromatography have pinpointed a subset of M. polymorpha miRNAs, notably miR11707, that remain in free form within the cell at 22°C but are loaded into MpAGO1 at 28°C to engage in RNA silencing. Investigations into the mir11707 gene editing (mir11707ge) mutants provided evidence of the regulation of miR11707 in MpAGO1. Notably, while MpAGO1 mRNA expression decreases at 28°C, the stability of the MpAGO1 protein and its associated miRNAs is essential for enhancing the RNA-inducing silencing complex (RISC) activity, revealing the importance of RNA silencing in enabling M. polymorpha to survive thermal stress. This study advances our understanding of RNA silencing in bryophytes and provides groundbreaking insights into the evolutionary resilience of land plants to climatic adversities.

An efficient sulfadiazine selection scheme for stable transformation in the model liverwort Marchantia polymorpha

Plant macroevolutionary studies leverage the phylogenetic position of non-flowering model systems like the liverwort Marchantia polymorpha to investigate the origin and evolution of key plant processes. To date, most molecular genetic studies in Marchantia rely on hygromycin and/or chlorsulfuron herbicide resistance markers for the selection of stable transformants. Here, we used a sulfonamide-resistant dihydropteroate synthase (DHPS) gene to enable sulfadiazine-based transformation selection in M. polymorpha. We demonstrate the reliability of sulfadiazine selection on its own and in combination with existing hygromycin and chlorsulfuron selection schemes through transgene stacking experiments. The utility of this system is further demonstrated through confocal microscopy of a triple transgenic line carrying fluorescent proteins labelling the plasma membrane, cortical microtubules, and the nucleus. Collectively, our findings and resources broaden the capacity to genetically manipulate the increasingly popular model liverwort M. polymorpha.

Streamlined regulation of chloroplast development in the liverwort Marchantia polymorpha

Chloroplasts develop from undifferentiated plastids in response to light. In angiosperms, after the perception of light, the Elongated Hypocotyl 5 (HY5) transcription factor initiates photomorphogenesis, and two families of transcription factors known as GOLDEN2-LIKE (GLK) and GATA are considered master regulators of chloroplast development. In addition, the MIR171-targeted SCARECROW-LIKE GRAS transcription factors also impact chlorophyll biosynthesis. The extent to which these proteins carry out conserved roles in non-seed plants is not known. Using the model liverwort Marchantia polymorpha, we show that GLK controls chloroplast biogenesis, and HY5 shows a small conditional effect on chlorophyll content. Chromatin immunoprecipitation sequencing (ChIP-seq) revealed that MpGLK has a broader set of targets than has been reported in angiosperms. We also identified a functional GLK homolog in green algae. In summary, our data support the hypothesis that GLK carries out a conserved role relating to chloroplast biogenesis in land plants and green algae.

MYB-related transcription factors control chloroplast biogenesis

Chloroplast biogenesis is dependent on master regulators from the GOLDEN2-LIKE (GLK) family of transcription factors. However, glk mutants contain residual chlorophyll, indicating that other proteins must be involved. Here, we identify MYB-related transcription factors as regulators of chloroplast biogenesis in the liverwort Marchantia polymorpha and angiosperm Arabidopsis thaliana. In both species, double-mutant alleles in MYB-related genes show very limited chloroplast development, and photosynthesis gene expression is perturbed to a greater extent than in GLK mutants. Genes encoding enzymes of chlorophyll biosynthesis are controlled by MYB-related and GLK proteins, whereas those allowing CO2 fixation, photorespiration, and photosystem assembly and repair require MYB-related proteins. Regulation between the MYB-related and GLK transcription factors appears more extensive in A. thaliana than in M. polymorpha. Thus, MYB-related and GLK genes have overlapping as well as distinct targets. We conclude that MYB-related and GLK transcription factors orchestrate chloroplast development in land plants.

The PLETHORA Homolog in Marchantia polymorpha is Essential for Meristem Maintenance, Developmental Progression, and Redox Homeostasis

To adapt to a terrestrial habitat, the ancestors of land plants must have made several morphological and physiological modifications, such as a meristem allowing for three-dimensional growth, rhizoids for water and nutrient uptake, air pore complexes or stomata that permit air exchange, and a defense system to cope with oxidative stress that occurs frequently in a terrestrial habitat. To understand how the meristem was determined during land plant evolution, we characterized the function of the closest PLETHORA homolog in the liverwort Marchantia polymorpha, which we named MpPLT. Through a transgenic approach, we showed that MpPLT is expressed not only in the stem cells at the apical notch but also in the proliferation zone of the meristem, as well as in cells that form the air-pore complex and rhizoids. Using the CRISPR method we then created mutants for MpPLT and found that the mutants are not only defective in meristem maintenance but also compromised in air-pore complex and rhizoid development. Strikingly, at later developmental stages, numerous gemma-like structures were formed in Mpplt mutants, suggesting developmental arrest. Further experiments indicated that MpPLT promotes plant growth by regulating MpWOX, which shared a similar expression pattern to MpPLT, and genes involved in auxin and cytokinin signaling pathways. Through transcriptome analyses, we found that MpPLT also has a role in redox homeostasis and that this role is essential for plant growth. Taken together, these results suggest that MpPLT has a crucial role in liverwort growth and development and hence may have played a crucial role in early land plant evolution.

The BNB-GLID module regulates germline fate determination in Marchantia polymorpha

Germline fate determination is a critical event in sexual reproduction. Unlike animals, plants specify the germline by reprogramming somatic cells at the late stages of their development. However, the genetic basis of germline fate determination and how it evolved during the land plant evolution are still poorly understood. Here, we report that the plant homeodomain finger protein GERMLINE IDENTITY DETERMINANT (GLID) is a key regulator of the germline specification in liverwort, Marchantia polymorpha. Loss of the MpGLID function causes failure of germline initiation, leading to the absence of sperm and egg cells. Remarkably, the overexpression of MpGLID in M. polymorpha induces the ectopic formation of cells with male germline cell features exclusively in male thalli. We further show that MpBONOBO (BNB), with an evolutionarily conserved function, can induce the formation of male germ cell-like cells through the activation of MpGLID by directly binding to its promoter. The Arabidopsis (Arabidopsis thaliana) MpGLID ortholog, MALE STERILITY1 (AtMS1), fails to replace the germline specification function of MpGLID in M. polymorpha, demonstrating that a derived function of MpGLID orthologs has been restricted to tapetum development in flowering plants. Collectively, our findings suggest the presence of the BNB-GLID module in complex ancestral land plants that has been retained in bryophytes, but rewired in flowering plants for male germline fate determination.

An ancient role for CYP73 monooxygenases in phenylpropanoid biosynthesis and embryophyte development

The phenylpropanoid pathway is one of the plant metabolic pathways most prominently linked to the transition to terrestrial life, but its evolution and early functions remain elusive. Here, we show that activity of the t-cinnamic acid 4-hydroxylase (C4H), the first plant-specific step in the pathway, emerged concomitantly with the CYP73 gene family in a common ancestor of embryophytes. Through structural studies, we identify conserved CYP73 residues, including a crucial arginine, that have supported C4H activity since the early stages of its evolution. We further demonstrate that impairing C4H function via CYP73 gene inactivation or inhibitor treatment in three bryophyte species-the moss Physcomitrium patens, the liverwort Marchantia polymorpha and the hornwort Anthoceros agrestis-consistently resulted in a shortage of phenylpropanoids and abnormal plant development. The latter could be rescued in the moss by exogenous supply of p-coumaric acid, the product of C4H. Our findings establish the emergence of the CYP73 gene family as a foundational event in the development of the plant phenylpropanoid pathway, and underscore the deep-rooted function of the C4H enzyme in embryophyte biology.

Contrasting and conserved roles of NPR pathways in diverged land plant lineages

The NPR proteins function as salicylic acid (SA) receptors in Arabidopsis thaliana. AtNPR1 plays a central role in SA-induced transcriptional reprogramming whereby positively regulates SA-mediated defense. NPRs are found in the genomes of nearly all land plants. However, we know little about the molecular functions and physiological roles of NPRs in most plant species. We conducted phylogenetic and alignment analyses of NPRs from 68 species covering the significant lineages of land plants. To investigate NPR functions in bryophyte lineages, we generated and characterized NPR loss-of-function mutants in the liverwort Marchantia polymorpha. Brassicaceae NPR1-like proteins have characteristically gained or lost functional residues identified in AtNPRs, pointing to the possibility of a unique evolutionary trajectory for the Brassicaceae NPR1-like proteins. We find that the only NPR in M. polymorpha, MpNPR, is not the master regulator of SA-induced transcriptional reprogramming and negatively regulates bacterial resistance in this species. The Mpnpr transcriptome suggested roles of MpNPR in heat and far-red light responses. We identify both Mpnpr and Atnpr1-1 display enhanced thermomorphogenesis. Interspecies complementation analysis indicated that the molecular properties of AtNPR1 and MpNPR are partially conserved. We further show that MpNPR has SA-binding activity. NPRs and NPR-associated pathways have evolved distinctively in diverged land plant lineages to cope with different terrestrial environments.

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, which 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.

Hyperspectral imaging of liverwort Marchantia polymorpha identifies MpWRKY10 as a key regulator defining Foliar pigmentation patterns

Foliar pigmentation patterns vary among plant species and growth conditions. In this study, we utilize hyperspectral imaging to assess foliar pigmentation in the bryophyte Marchantia polymorpha under nutrient stress and identify associated genetic factors. Using singular value decomposition (SVD) for feature selection, we quantitate color variations induced by deficiencies in phosphate, nitrate, magnesium, calcium, and iron. Pseudo-colored thallus images show that disrupting MpWRKY10 causes irregular pigmentation with auronidin accumulation. Transcriptomic profiling shows that MpWRKY10 regulates phenylpropanoid pathway enzymes and R2R3-MYB transcription factors during phosphate deficiency, with MpMYB14 upregulation preceding pigment accumulation. MpWRKY10 is downregulated in older, pigmented thalli under phosphate deficiency but maintained in young thalli, where it suppresses pigmentation genes. This downregulation is absent in pigmented thalli due to aging. Comparative transcriptome analysis suggests similar WRKY and MYB roles in nutrient response and pigmentation in red-leaf lettuce, alluding to conserved genetic factors controlling foliar pigmentation patterns under nutrient deficiency.

MpANT regulates meristem development in Marchantia polymorpha

Meristems are crucial for organ formation, but our knowledge of their molecular evolution is limited. Here, we show that AINTEGUMENTA (MpANT) in the euANT branch of the APETALA2-like transcription factor family is essential for meristem development in the nonvascular plant Marchantia polymorpha. MpANT is expressed in the thallus meristem. Mpant mutants show defects to maintain meristem identity and undergo meristem duplication, while MpANT overexpressers show ectopic thallus growth. MpANT directly upregulates MpGRAS9 in the SHORT-ROOT (SHR) branch of the GRAS family. In the vascular plant Arabidopsis thaliana, the euANT-branch genes PLETHORAs (AtPLTs) and AtANT are involved in the formation and maintenance of root/shoot apical meristems and lateral organ primordia, and AtPLTs directly target SHR-branch genes. In addition, euANTs bind through a similar DNA-binding motif to many conserved homologous genes in M. polymorpha and A. thaliana. Overall, the euANT pathway has an evolutionarily conserved role in meristem development.

MpMLO1 controls sperm discharge in liverwort

In bryophytes, sexual reproduction necessitates the release of motile sperm cells from a gametophyte into the environment. Since 1856, this process, particularly in liverworts, has been known to depend on water. However, the molecular mechanism underlying this phenomenon has remained elusive. Here we identify the plasma membrane protein MpMLO1 in Marchantia polymorpha, a model liverwort, as critical for sperm discharge from antheridia. The MpMLO1-expressing tip cells among the sperm-wrapping jacket cells undergo programmed cell death upon antheridium maturation to facilitate sperm discharge after the application of water and even hypertonic solutions. The absence of MpMLO1 leads to reduced cytoplasmic Ca2+ levels in tip cells, preventing cell death and consequently sperm discharge. Our findings reveal that MpMLO1-mediated programmed cell death in antheridial tip cells, regulated by cytosolic Ca2+ dynamics, is essential for sperm release, elucidating a key mechanism in bryophyte sexual reproduction and providing insights into terrestrial plant evolution.

A RabGAP negatively regulates plant autophagy and immune trafficking

Plants rely on autophagy and membrane trafficking to tolerate stress, combat infections, and maintain cellular homeostasis. However, the molecular interplay between autophagy and membrane trafficking is poorly understood. Using an AI-assisted approach, we identified Rab3GAP-like (Rab3GAPL) as a key membrane trafficking node that controls plant autophagy negatively. Rab3GAPL suppresses autophagy by binding to ATG8, the core autophagy adaptor, and deactivating Rab8a, a small GTPase essential for autophagosome formation and defense-related secretion. Rab3GAPL reduces autophagic flux in three model plant species, suggesting that its negative regulatory role in autophagy is conserved in land plants. Beyond autophagy regulation, Rab3GAPL modulates focal immunity against the oomycete pathogen Phytophthora infestans by preventing defense-related secretion. Altogether, our results suggest that Rab3GAPL acts as a molecular rheostat to coordinate autophagic flux and defense-related secretion by restraining Rab8a-mediated trafficking. This unprecedented interplay between a RabGAP-Rab pair and ATG8 sheds new light on the intricate membrane transport mechanisms underlying plant autophagy and immunity.

Rapid Propagation of Ca2+ Waves and Electrical Signals in the Liverwort Marchantia polymorpha

In response to both biotic and abiotic stresses, vascular plants transmit long-distance Ca2+ and electrical signals from localized stress sites to distant tissues through their vasculature. Various models have been proposed for the mechanisms underlying the long-distance signaling, primarily centered around the presence of vascular bundles. We here demonstrate that the non-vascular liverwort Marchantia polymorpha possesses a mechanism for propagating Ca2+ waves and electrical signals in response to wounding. The propagation velocity of these signals was approximately 1-2 mm s-1, equivalent to that observed in vascular plants. Both Ca2+ waves and electrical signals were inhibited by La3+ as well as tetraethylammonium chloride, suggesting the crucial importance of both Ca2+ channel(s) and K+ channel(s) in wound-induced membrane depolarization as well as the subsequent long-distance signal propagation. Simultaneous recordings of Ca2+ and electrical signals indicated a tight coupling between the dynamics of these two signaling modalities. Furthermore, molecular genetic studies revealed that a GLUTAMATE RECEPTOR-LIKE (GLR) channel plays a central role in the propagation of both Ca2+ waves and electrical signals. Conversely, none of the three two-pore channels were implicated in either signal propagation. These findings shed light on the evolutionary conservation of rapid long-distance Ca2+ wave and electrical signal propagation involving GLRs in land plants, even in the absence of vascular tissue.

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.

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.

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.

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.

Divergent evolution of the alcohol-forming pathway of wax biosynthesis among bryophytes

The plant cuticle is a hydrophobic barrier, which seals the epidermal surface of most aboveground organs. While the cuticle biosynthesis of angiosperms has been intensively studied, knowledge about its existence and composition in nonvascular plants is scarce. Here, we identified and characterized homologs of Arabidopsis thaliana fatty acyl-CoA reductase (FAR) ECERIFERUM 4 (AtCER4) and bifunctional wax ester synthase/acyl-CoA:diacylglycerol acyltransferase 1 (AtWSD1) in the liverwort Marchantia polymorpha (MpFAR2 and MpWSD1) and the moss Physcomitrium patens (PpFAR2A, PpFAR2B, and PpWSD1). Although bryophyte harbor similar compound classes as described for angiosperm cuticles, their biosynthesis may not be fully conserved between the bryophytes M. polymorpha and P. patens or between these bryophytes and angiosperms. While PpFAR2A and PpFAR2B contribute to the production of primary alcohols in P. patens, loss of MpFAR2 function does not affect the wax profile of M. polymorpha. By contrast, MpWSD1 acts as the major wax ester-producing enzyme in M. polymorpha, whereas mutations of PpWSD1 do not affect the wax ester levels of P. patens. Our results suggest that the biosynthetic enzymes involved in primary alcohol and wax ester formation in land plants have either evolved multiple times independently or undergone pronounced radiation followed by the formation of lineage-specific toolkits.

A subclass II bHLH transcription factor in Marchantia polymorpha gives insight into the ancestral land plant trait of spore formation

Sporopollenin is often said to be one of the toughest biopolymers known to man. The shift in dormancy cell wall deposition from around the diploid zygotes of charophycean algae to sporopollenin around the haploid spores of land plants essentially imparted onto land plants the gift of passive motility, a key acquisition that contributed to their vast and successful colonization across terrestrial habitats.1,2 A putative transcription factor controlling the land plant mode of sporopollenin deposition is the subclass II bHLHs, which are conserved and novel to land plants, with mutants of genes in angiosperms and mosses divulging roles relating to tapetum degeneration and spore development.3,4,5,6,7 We demonstrate that a subclass II bHLH gene, MpbHLH37, regulates sporopollenin biosynthesis and deposition in the model liverwort Marchantia polymorpha. Mpbhlh37 sporophytes show a striking loss of secondary wall deposits of the capsule wall, the elaters, and the spore exine, all while maintaining spore viability, identifying MpbHLH37 as a master regulator of secondary wall deposits of the sporophyte. Localization of MpbHLH37 to the capsule wall and elaters of the sporophyte directly designates these tissue types as a bona fide tapetum in liverworts, giving support to the notion that the presence of a tapetum is an ancestral land plant trait. Finally, as early land plant spore walls exhibit evidence of tapetal deposition,8,9,10,11,12 a tapetal capsule wall could have provided these plants with a developmental mechanism for sporopollenin deposition.

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.

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.

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.

Analysis of Plant-Specific ANTH Domain-Containing Protein in Marchantia polymorpha

Membrane trafficking is a fundamental mechanism for protein and lipid transport in eukaryotic cells and exhibits marked diversity among eukaryotic lineages with distinctive body plans and lifestyles. Diversification of the membrane trafficking system is associated with the expansion and secondary loss of key machinery components, including RAB GTPases, soluble N-ethylmaleimide-sensitive factor attachment protein receptors (SNAREs) and adaptor proteins, during plant evolution. The number of AP180 N-terminal homology (ANTH) proteins, an adaptor family that regulates vesicle formation and cargo sorting during clathrin-mediated endocytosis, increases during plant evolution. In the genome of Arabidopsis thaliana, 18 genes for ANTH proteins have been identified, a higher number than that in yeast and animals, suggesting a distinctive diversification of ANTH proteins. Conversely, the liverwort Marchantia polymorpha possesses a simpler repertoire; only two genes encoding canonical ANTH proteins have been identified in its genome. Intriguingly, a non-canonical ANTH protein is encoded in the genome of M. polymorpha, which also harbors a putative kinase domain. Similar proteins have been detected in sporadic lineages of plants, suggesting their ancient origin and multiple secondary losses during evolution. We named this unique ANTH group phosphatidylinositol-binding clathrin assembly protein-K (PICALM-K) and characterized it in M. polymorpha using genetic, cell biology-based and artificial intelligence (AI)-based approaches. Our results indicate a flagella-related function of MpPICALM-K in spermatozoids, which is distinct from that of canonical ANTH proteins. Therefore, ANTH proteins have undergone significant functional diversification during evolution, and PICALM-K represents a plant-unique ANTH protein that is delivered by neofunctionalization through exon shuffling.

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.

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.

Circadian clock does not play an essential role in daylength measurement for growth-phase transition in Marchantia polymorpha

Daylength is perceived as a seasonal cue to induce growth-phase transition at a proper time of a year. The core of the mechanism of daylength measurement in angiosperms lies in the circadian clock-controlled expression of regulators of growth-phase transition. However, the roles of the circadian clock in daylength measurement in basal land plants remain largely unknown. In this study, we investigated the contribution of circadian clock to daylength measurement in a basal land plant, the liverwort Marchantia polymorpha. In M. polymorpha, transition from vegetative to reproductive phase under long-day conditions results in differentiation of sexual branches called gametangiophores which harbor gametangia. First, we showed that a widely used wild-type accession Takaragaike-1 is an obligate long-day plant with a critical daylength of about 10 hours and requires multiple long days. Then, we compared the timing of gametangiophore formation between wild type and circadian clock mutants in long-day and short-day conditions. Mutations in two clock genes, MpTIMING OF CAB EXPRESSION 1 and MpPSEUDO-RESPONSE REGULATOR, had no significant effects on the timing of gametangiophore formation. In addition, when M. polymorpha plants were treated with a chemical which lengthens circadian period, there was no significant effect on the timing of gametangiophore formation, either. We next observed the timing of gametangiophore formation under various non-24-h light/dark cycles to examine the effect of phase alteration in circadian rhythms. The results suggest that daylength measurement in M. polymorpha is based on the relative amount of light and darkness within a cycle rather than the intrinsic rhythms generated by circadian clock. Our findings suggest that M. polymorpha has a daylength measurement system which is different from that of angiosperms centered on the circadian clock function.

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.

The Polycomb repressive complex 2 deposits H3K27me3 and represses transposable elements in a broad range of eukaryotes

The mobility of transposable elements (TEs) contributes to evolution of genomes. Their uncontrolled activity causes genomic instability; therefore, expression of TEs is silenced by host genomes. TEs are marked with DNA and H3K9 methylation, which are associated with silencing in flowering plants, animals, and fungi. However, in distantly related groups of eukaryotes, TEs are marked by H3K27me3 deposited by the Polycomb repressive complex 2 (PRC2), an epigenetic mark associated with gene silencing in flowering plants and animals. The direct silencing of TEs by PRC2 has so far only been shown in one species of ciliates. To test if PRC2 silences TEs in a broader range of eukaryotes, we generated mutants with reduced PRC2 activity and analyzed the role of PRC2 in extant species along the lineage of Archaeplastida and in the diatom P. tricornutum. In this diatom and the red alga C. merolae, a greater proportion of TEs than genes were repressed by PRC2, whereas a greater proportion of genes than TEs were repressed by PRC2 in bryophytes. In flowering plants, TEs contained potential cis-elements recognized by transcription factors and associated with neighbor genes as transcriptional units repressed by PRC2. Thus, silencing of TEs by PRC2 is observed not only in Archaeplastida but also in diatoms and ciliates, suggesting that PRC2 deposited H3K27me3 to silence TEs in the last common ancestor of eukaryotes. We hypothesize that during the evolution of Archaeplastida, TE fragments marked with H3K27me3 were selected to shape transcriptional regulation, controlling networks of genes regulated by PRC2.

Phosphatidic acid phosphohydrolase modulates glycerolipid synthesis in Marchantia polymorpha and is crucial for growth under both nutrient-replete and -deficient conditions

MAIN CONCLUSION

The phosphatidic acid phosphohydrolase of Marchantia polymorpha modulates plastid glycolipid synthesis through the ER pathway and is essential for normal plant development regardless of nutrient availability. Membrane lipid remodeling is one of the strategies plant cells use to secure inorganic phosphate (Pi) for plant growth, but many aspects of the molecular mechanism and its regulation remain unclear. Here we analyzed membrane lipid remodeling using a non-vascular plant, Marchantia polymorpha. The lipid composition and fatty acid profile during Pi starvation in M. polymorpha revealed a decrease in phospholipids and an increase in both galactolipids and betaine lipids. In Arabidopsis thaliana, phosphatidic acid phosphohydrolase (PAH) is involved in phospholipid degradation and is crucial for tolerance to both Pi and nitrogen starvation. We produced two M. polymorpha PAH (MpPAH) knockout mutants (Mppah-1 and Mppah-2) and found that, unlike Arabidopsis mutants, Mppah impaired plant growth with shorter rhizoids compared with wild-type plants even under nutrient-replete conditions. Mutation of MpPAH did not significantly affect the mole percent of each glycerolipid among total membrane glycerolipids from whole plants under both Pi-replete and Pi-deficient conditions. However, the fatty acid composition of monogalactosyldiacylglycerol indicated that the amount of plastid glycolipids produced through the endoplasmic reticulum pathway was suppressed in Mppah mutants. Phospholipids accumulated in the mutants under N starvation. These results reveal that MpPAH modulates plastid glycolipid synthesis through the endoplasmic reticulum pathway more so than what has been observed for Arabidopsis PAH; moreover, unlike Arabidopsis, MpPAH is crucial for M. polymorpha growth regardless of nutrient availability.

Activation of Strigolactone Biosynthesis by the DWARF14-LIKE/KARRIKIN-INSENSITIVE2 Pathway in Mycorrhizal Angiosperms, but Not in Arabidopsis, a Non-mycorrhizal Plant

Strigolactones (SLs) are a class of plant hormones that regulate many aspects of plant growth and development. SLs also improve symbiosis with arbuscular mycorrhizal fungi (AMF) in the rhizosphere. Recent studies have shown that the DWARF14-LIKE (D14L)/KARRIKIN-INSENSITIVE2 (KAI2) family, paralogs of the SL receptor D14, are required for AMF colonization in several flowering plants, including rice. In this study, we found that (-)-GR5, a 2'S-configured enantiomer of a synthetic SL analog (+)-GR5, significantly activated SL biosynthesis in rice roots via D14L. This result is consistent with a recent report, showing that the D14L pathway positively regulates SL biosynthesis in rice. In fact, the SL levels tended to be lower in the roots of the d14l mutant under both inorganic nutrient-deficient and -sufficient conditions. We also show that the increase in SL levels by (-)-GR5 was observed in other mycorrhizal plant species. In contrast, the KAI2 pathway did not upregulate the SL level and the expression of SL biosynthetic genes in Arabidopsis, a non-mycorrhizal plant. We also examined whether the KAI2 pathway enhances SL biosynthesis in the liverwort Marchantia paleacea, where SL functions as a rhizosphere signaling molecule for AMF. However, the SL level and SL biosynthetic genes were not positively regulated by the KAI2 pathway. These results imply that the activation of SL biosynthesis by the D14L/KAI2 pathway has been evolutionarily acquired after the divergence of bryophytes to efficiently promote symbiosis with AMF, although we cannot exclude the possibility that liverworts have specifically lost this regulatory system.

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.

DIENELACTONE HYDROLASE LIKE PROTEIN1 negatively regulates the KAI2-ligand pathway in Marchantia polymorpha

Karrikins are smoke-derived butenolides that induce seed germination and photomorphogenesis in a wide range of plants.1,2,3 KARRIKIN INSENSITIVE2 (KAI2), a paralog of a strigolactone receptor, perceives karrikins or their metabolized products in Arabidopsis thaliana.4,5,6,7 Furthermore, KAI2 is thought to perceive an unidentified plant hormone, called KAI2 ligand (KL).8,9 KL signal is transduced via the interaction between KAI2, MORE AXILLARY GROWTH2 (MAX2), and SUPPRESSOR of MORE AXILLARY GROWTH2 1 LIKE family proteins (SMXLs), followed by the degradation of SMXLs.4,7,10,11,12,13,14 This signaling pathway is conserved both in A. thaliana and the bryophyte Marchantia polymorpha.14 Although the KL signaling pathway is well characterized, the KL metabolism pathways remain poorly understood. Here, we show that DIENELACTONE HYDROLASE LIKE PROTEIN1 (DLP1) is a negative regulator of the KL pathway in M. polymorpha. The KL signal induces DLP1 expression. DLP1 overexpression lines phenocopied the Mpkai2a and Mpmax2 mutants, while dlp1 mutants phenocopied the Mpsmxl mutants. Mutations in the KL signaling genes largely suppressed these phenotypes, indicating that DLP1 acts upstream of the KL signaling pathway, although DLP1 also has KL pathway-independent functions. DLP1 exhibited enzymatic activity toward a potential substrate, suggesting the possibility that DLP1 works through KL inactivation. Investigation of DLP1 homologs in A. thaliana revealed that they do not play a major role in the KL pathway, suggesting different mechanisms for the KL signal regulation. Our findings provide new insights into the regulation of the KL signal in M. polymorpha and the evolution of the KL pathway in land plants.

Bioinformatic and literature assessment of toxicity and allergenicity of a CRISPR-Cas9 engineered gene drive to control Anopheles gambiae the mosquito vector of human malaria

BACKGROUND

Population suppression gene drive is currently being evaluated, including via environmental risk assessment (ERA), for malaria vector control. One such gene drive involves the dsxFCRISPRh transgene encoding (i) hCas9 endonuclease, (ii) T1 guide RNA (gRNA) targeting the doublesex locus, and (iii) DsRed fluorescent marker protein, in genetically-modified mosquitoes (GMMs). Problem formulation, the first stage of ERA, for environmental releases of dsxFCRISPRh previously identified nine potential harms to the environment or health that could occur, should expressed products of the transgene cause allergenicity or toxicity.

METHODS

Amino acid sequences of hCas9 and DsRed were interrogated against those of toxins or allergens from NCBI, UniProt, COMPARE and AllergenOnline bioinformatic databases and the gRNA was compared with microRNAs from the miRBase database for potential impacts on gene expression associated with toxicity or allergenicity. PubMed was also searched for any evidence of toxicity or allergenicity of Cas9 or DsRed, or of the donor organisms from which these products were originally derived.

RESULTS

While Cas9 nuclease activity can be toxic to some cell types in vitro and hCas9 was found to share homology with the prokaryotic toxin VapC, there was no evidence from previous studies of a risk of toxicity to humans and other animals from hCas9. Although hCas9 did contain an 8-mer epitope found in the latex allergen Hev b 9, the full amino acid sequence of hCas9 was not homologous to any known allergens. Combined with a lack of evidence in the literature of Cas9 allergenicity, this indicated negligible risk to humans of allergenicity from hCas9. No matches were found between the gRNA and microRNAs from either Anopheles or humans. Moreover, potential exposure to dsxFCRISPRh transgenic proteins from environmental releases was assessed as negligible.

CONCLUSIONS

Bioinformatic and literature assessments found no convincing evidence to suggest that transgenic products expressed from dsxFCRISPRh were allergens or toxins, indicating that environmental releases of this population suppression gene drive for malaria vector control should not result in any increased allergenicity or toxicity in humans or animals. These results should also inform evaluations of other GMMs being developed for vector control and in vivo clinical applications of CRISPR-Cas9.

RHO GTPase of plants regulates polarized cell growth and cell division orientation during morphogenesis

Cell polarity-broadly defined as the asymmetric distribution of cellular activities and subcellular components within a cell-determines the geometry of cell growth and division during development. RHO GTPase proteins regulate the establishment of cell polarity and are conserved among eukaryotes. RHO of plant (ROP) proteins are a subgroup of RHO GTPases that are required for cellular morphogenesis in plants. However, how ROP proteins modulate the geometry of cell growth and division during the morphogenesis of plant tissues and organs is not well understood. To investigate how ROP proteins function during tissue development and organogenesis, we characterized the function of the single-copy ROP gene of the liverwort Marchantia polymorpha (MpROP). M. polymorpha develops morphologically complex three-dimensional tissues and organs exemplified by air chambers and gemmae, respectively. Mprop loss-of-function mutants form defective air chambers and gemmae, indicating ROP function is required for tissue development and organogenesis. During air chamber and gemma development in wild type, the MpROP protein is enriched to sites of polarized growth at the cell surface and accumulates at the expanding cell plate of dividing cells. Consistent with these observations, polarized cell growth is lost and cell divisions are misoriented in Mprop mutants. We propose that ROP regulates both polarized cell growth and cell division orientation in a coordinated manner to orchestrate tissue development and organogenesis in land plants.

Liverwort bHLH transcription factors and the origin of stomata in plants

Stomata are distributed in nearly all major groups of land plants, with the only exception being liverworts. Instead of having stomata on sporophytes, many complex thalloid liverworts possess air pores in their gametophytes. At present, whether stomata in land plants are derived from a common origin remains under debate.^1,2,3 In Arabidopsis thaliana, a core regulatory module for stomatal development comprises members of the bHLH transcription factor (TF) family, including AtSPCH, AtMUTE, and AtFAMA of subfamily Ia and AtSCRM1/2 of subfamily IIIb. Specifically, AtSPCH, AtMUTE, and AtFAMA each successively form heterodimers with AtSCRM1/2, which in turn regulate the entry, division, and differentiation of stomatal lineages.^4,5,6,7 In the moss Physcomitrium patens, two SMF (SPCH, MUTE and FAMA) orthologs have been characterized, one of which is functionally conserved in regulating stomatal development.^8,9 We here provide experimental evidence that orthologous bHLH TFs in the liverwort Marchantia polymorpha affect air pore spacing as well as the development of the epidermis and gametangiophores. We found that the bHLH Ia and IIIb heterodimeric module is highly conserved in plants. Genetic complementation experiments showed that liverwort SCRM and SMF genes weakly restored a stomata phenotype in atscrm1, atmute, and atfama mutant backgrounds in A. thaliana. In addition, homologs of stomatal development regulators FLP and MYB88 also exist in liverworts and weakly rescued the stomatal phenotype of atflp/myb88 double mutant. These results provide evidence not only for a common origin of all stomata in extant plants but also for relatively simple stomata in the ancestral plant.

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.

PIN-FORMED is required for shoot phototropism/gravitropism and facilitates meristem formation in Marchantia polymorpha

PIN-FORMED auxin efflux transporters, a subclass of which is plasma membrane-localised, mediate a variety of land-plant developmental processes via their polar localisation and subsequent directional auxin transport. We provide the first characterisation of PIN proteins in liverworts using Marchantia polymorpha as a model system. Marchantia polymorpha possesses a single PIN-FORMED gene, whose protein product is predicted to be plasma membrane-localised, MpPIN1. To characterise MpPIN1, we created loss-of-function alleles and produced complementation lines in both M. polymorpha and Arabidopsis. In M. polymorpha, gene expression and protein localisation were tracked using an MpPIN1 transgene encoding a translationally fused fluorescent protein. Overexpression of MpPIN1 can partially complement loss of an orthologous gene, PIN-FORMED1, in Arabidopsis. In M. polymorpha, MpPIN1 influences development in numerous ways throughout its life cycle. Most notably, MpPIN1 is required to establish gemmaling dorsiventral polarity and for orthotropic growth of gametangiophore stalks, where MpPIN1 is basally polarised. PIN activity is largely conserved within land plants, with PIN-mediated auxin flow providing a flexible mechanism to organise growth. Specifically, PIN is fundamentally linked to orthotropism and to the establishment of de novo meristems, the latter potentially involving the formation of both auxin biosynthesis maxima and auxin-signalling minima.

Chemical Control of CRISPR Gene Editing via Conditional Diacylation Crosslinking of Guide RNAs

Conditional control of RNA structure and function has emerged as an effective toolkit. Here, a strategy based on a one-step introduction of diacylation linkers and azide groups on the 2'-OH of RNA is advance. Selected from eight phosphine reagents, it is found that 2-(diphenylphosphino)ethylamine has excellent performance in reducing azides via a Staudinger reduction to obtain the original RNA. It is demonstrated that the enzymatic activities of Cas13 and Cas9 can be regulated by chemically modified guide RNAs, and further achieved ligand-induced gene editing in living cells by a controllable CRISPR/Cas9 system.

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.

What can hornworts teach us?

The hornworts are a small group of land plants, consisting of only 11 families and approximately 220 species. Despite their small size as a group, their phylogenetic position and unique biology are of great importance. Hornworts, together with mosses and liverworts, form the monophyletic group of bryophytes that is sister to all other land plants (Tracheophytes). It is only recently that hornworts became amenable to experimental investigation with the establishment of Anthoceros agrestis as a model system. In this perspective, we summarize the recent advances in the development of A. agrestis as an experimental system and compare it with other plant model systems. We also discuss how A. agrestis can help to further research in comparative developmental studies across land plants and to solve key questions of plant biology associated with the colonization of the terrestrial environment. Finally, we explore the significance of A. agrestis in crop improvement and synthetic biology applications in general.

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.

Meristem dormancy in Marchantia polymorpha is regulated by a liverwort-specific miRNA and a clade III SPL gene

The shape of modular organisms depends on the branching architecture, which in plants is determined by the fates of generative centers called meristems. The branches of the liverwort Marchantia polymorpha are derived from two adjacent meristems that develop at thallus apices. These meristems may be active and develop branches or may be dormant and do not form branches. The relative number and position of active and dormant meristems define the overall shape and form of the thallus. We show that the clade III SQUAMOSA PROMOTER BINDING PROTEIN-LIKE (SPL) transcription factor MpSPL1 is required for meristem dormancy. The activity of MpSPL1 is regulated by the liverwort-specific Mpo-MR13 miRNA, which, in turn, is regulated by PIF-mediated signaling. An unrelated PIF-regulated miRNA, MIR156, represses a different SPL gene (belonging to clade IV) that inhibits branching during the shade avoidance response in Arabidopsis thaliana. This suggests that a conserved light signaling mechanism modulates branching architecture in liverworts and angiosperms and therefore is likely operated in the last common ancestor. However, PIF-mediated signaling represses the expression of different miRNA genes with different SPL targets during dichotomous, apical branching in liverworts and during lateral, subapical branching in angiosperms. We speculate that the mechanism that acts downstream of light and regulates meristem dormancy evolved independently in liverworts and angiosperms.

Control of vegetative reproduction in Marchantiapolymorpha by the KAI2-ligand signaling pathway

In vegetative reproduction of Marchantia polymorpha (M. polymorpha), propagules, called gemmae, are formed in gemma cups. Despite its significance for survival, control of gemma and gemma cup formation by environmental cues is not well understood. We show here that the number of gemmae formed in a gemma cup is a genetic trait. Gemma formation starts from the central region of the floor of the gemma cup, proceeds to the periphery, and terminates when the appropriate number of gemmae is initiated. The MpKARRIKIN INSENSITIVE2 (MpKAI2)-dependent signaling pathway promotes gemma cup formation and gemma initiation. The number of gemmae in a cup is controlled by modulating the ON/OFF switch of the KAI2-dependent signaling. Termination of the signaling results in the accumulation of MpSMXL, a suppressor protein. In the Mpsmxl mutants, gemma initiation continues, leading to the formation of a highly increased number of gemmae in a cup. Consistent with its function, the MpKAI2-dependent signaling pathway is active in gemma cups where gemmae initiate, as well as in the notch region of the mature gemma and midrib of the ventral side of the thallus. In this work, we also show that GEMMA CUP-ASSOCIATED MYB1 works downstream of this signaling pathway to promote gemma cup formation and gemma initiation. We also found that the availability of potassium affects gemma cup formation independently from the KAI2-dependent signaling pathway in M. polymorpha. We propose that the KAI2-dependent signaling pathway functions to optimize vegetative reproduction by adapting to the environment 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.

Comprehensive analysis of peptide-coding genes and initial characterization of an LRR-only microprotein in Marchantia polymorpha

In the past two decades, many plant peptides have been found to play crucial roles in various biological events by mediating cell-to-cell communications. However, a large number of small open reading frames (sORFs) or short genes capable of encoding peptides remain uncharacterized. In this study, we examined several candidate genes for peptides conserved between two model plants: Arabidopsis thaliana and Marchantia polymorpha. We examined their expression pattern in M. polymorpha and subcellular localization using a transient assay with Nicotiana benthamiana. We found that one candidate, MpSGF10B, was expressed in meristems, gemma cups, and male reproductive organs called antheridiophores. MpSGF10B has an N-terminal signal peptide followed by two leucine-rich repeat (LRR) domains and was secreted to the extracellular region in N. benthamiana and M. polymorpha. Compared with the wild type, two independent Mpsgf10b mutants had a slightly increased number of antheridiophores. It was revealed in gene ontology enrichment analysis that MpSGF10B was significantly co-expressed with genes related to cell cycle and development. These results suggest that MpSGF10B may be involved in the reproductive development of M. polymorpha. Our research should shed light on the unknown role of LRR-only proteins in land plants.

Application of a novel strong promoter from Chinese fir (Cunninghamia lanceolate) in the CRISPR/Cas mediated genome editing of its protoplasts and transgenesis of rice and poplar

Novel constitutive promoters are essential for plant biotechnology. Although in angiosperms, a number of promoters were applied in monocots or dicots genetic engineering, only a few promoters were used in gymnosperm. Here we identified two strong promoters (Cula11 and Cula08) from Chinese fir (C. lanceolate) by screening the transcriptomic data and preliminary promoter activity assays in tobacco. By using the newly established Chinese fir protoplast transient expression technology that enables in vivo molecular biology studies in its homologous system, we compared the activities of Cula11 and Cula08 with that of the commonly used promoters in genetic engineering of monocots or dicots, such as CaM35S, CmYLCV, and ZmUbi, and our results revealed that Cula11 and Cula08 promoters have stronger activities in Chinese fir protoplasts. Furthermore, the vector containing Cas gene driven by Cula11 promoter and sgRNA driven by the newly isolated CulaU6b polyIII promoters were introduced into Chinese fir protoplasts, and CRISPR/Cas mediated gene knock-out event was successfully achieved. More importantly, compared with the commonly used promoters in the genetic engineering in angiosperms, Cula11 promoter has much stronger activity than CaM35S promoter in transgenic poplar, and ZmUbi promoter in transgenic rice, respectively, indicating its potential application in poplar and rice genetic engineering. Overall, the novel putative constitutive gene promoters reported here will have great potential application in gymnosperm and angiosperm biotechnology, and the transient gene expression system established here will serve as a useful tool for the molecular and genetic analyses of Chinese fir genes.

CRISPR/Cas9 Technology for Potato Functional Genomics and Breeding

Cultivated potato (Solanum tuberosum L.) is one of the most important staple food crops worldwide. Its tetraploid and highly heterozygous nature poses a great challenge to its basic research and trait improvement through traditional mutagenesis and/or crossbreeding. The establishment of the clustered regularly interspaced short palindromic repeats (CRISPR) and CRISPR-associated protein 9 (Cas9) as a gene editing tool has allowed the alteration of specific gene sequences and their concomitant gene function, providing powerful technology for potato gene functional analysis and improvement of elite cultivars. This technology relies on a short RNA molecule called single guide RNA (sgRNA) that directs the Cas9 nuclease to induce a site-specific double-stranded break (DSB). Further, repair of the DSB by the error-prone non-homologous end joining (NHEJ) mechanism leads to the introduction of targeted mutations, which can be used to produce the loss of function of specific gene(s). In this chapter, we describe experimental procedures to apply the CRISPR/Cas9 technology for potato genome editing. First, we provide strategies for target selection and sgRNA design and describe a Golden Gate-based cloning system to obtain a sgRNA/Cas9-encoding binary vector. We also describe an optimized protocol for ribonucleoprotein (RNP) complex assembly. The binary vector can be used for both Agrobacterium-mediated transformation and transient expression in potato protoplasts, while the RNP complexes are intended to obtain edited potato lines through protoplast transfection and plant regeneration. Finally, we describe procedures to identify the gene-edited potato lines. The methods described here are suitable for potato gene functional analysis and breeding.

Selection of a histidine auxotrophic Marchantia polymorpha strain with an auxotrophic selective marker

Marchantia polymorpha has emerged as a model liverwort species, with molecular tools increasingly available. In the present study, we developed an auxotrophic strain of M. polymorpha and an auxotrophic selective marker gene as new experimental tools for this valuable model system. Using CRISPR (clustered regularly interspaced palindromic repeats)/Cas9-mediated genome editing, we mutated the genomic region for IMIDAZOLEGLYCEROL-PHOSPHATE DEHYDRATASE (IGPD) in M. polymorpha to disrupt the biosynthesis of histidine (igpd). We modified an IGPD gene (IGPDm) with silent mutations, generating a histidine auxotrophic selective marker gene that was not a target of our CRISPR/Cas9-mediated genome editing. The M. polymorpha igpd mutant was a histidine auxotrophic strain, growing only on medium containing histidine. The igpd mutant could be complemented by transformation with the IGPDm gene, indicating that this gene could be used as an auxotrophic selective marker. Using the IGPDm marker in the igpd mutant background, we produced transgenic lines without the need for antibiotic selection. The histidine auxotrophic strain igpd and auxotrophic selective marker IGPDm represent new molecular tools for M. polymorpha research.