A snapshot of the Physcomitrella N-terminome reveals N-terminal methylation of organellar proteins.

KEY MESSAGE: Analysis of the N-terminome of Physcomitrella reveals N-terminal monomethylation of nuclear-encoded, mitochondria-localized proteins. Post- or co-translational N-terminal modifications of proteins influence their half-life as well as mediating protein sorting to organelles via cleavable N-terminal sequences that are recognized by the respective translocation machinery. Here, we provide an overview on the current modification state of the N-termini of over 4500 proteins from the model moss Physcomitrella (Physcomitrium patens) using a compilation of 24 N-terminomics datasets. Our data reveal distinct proteoforms and modification states and confirm predicted targeting peptide cleavage sites of 1,144 proteins localized to plastids and the thylakoid lumen, to mitochondria, and to the secretory pathway. In addition, we uncover extended N-terminal methylation of mitochondrial proteins. Moreover, we identified PpNTM1 (P. patens alpha N-terminal protein methyltransferase 1) as a candidate for protein methylation in plastids, mitochondria, and the cytosol. These data can now be used to optimize computational targeting predictors, for customized protein fusions and their targeted localization in biotechnology, and offer novel insights into potential dual targeting of proteins.

Nuclear interactions between the Pseudo-Response Regulator clock proteins and the Multi-Step Phosphorelay mediator Histidine-containing phosphotransfers in the moss Physcomitrium patens.

Pseudo-Response Regulator (PRR) proteins constitute a fundamental set of circadian clock components in plants. PRRs have an amino acid sequence stretch with similarity to the receiver (REC) domain of response regulators (RRs) in the Multi-Step Phosphorelay (MSP). However, it has never been elucidated whether PRRs interact with Histidine-containing Phosphotransfer (HPt) proteins, which transfer a phosphate to RRs. Here, we studied whether PRRs interact with HPts in the moss Physcomitrium patens by the Yeast Two-Hybrid system and Bimolecular Fluorescence Complementation. P. patens PRR1/2/3 interacted with HPt1/2 in the nucleus, but not with HPt3, suggesting that P. patens PRRs function as authentic RRs. We discuss these results in relation to the evolution and diversity of the plant circadian clocks.

Insights into a functional synthetic plant genome.

Synthetic genomics involves the design, assembly, and transfer of artificially synthesized DNA fragments into target hosts to replace the native genome and construct viable forms of life. With advances in DNA synthesis and assembly techniques, the application of synthetic genomics in viruses, bacteria, and yeast has improved our knowledge of genome organization and function. Multicellular eukaryotic organisms are characterized by larger genomes, more complex epigenetic regulation, and widespread transposable elements, making genome synthesis challenging. Recently, the first synthetic multicellular eukaryotic organism was generated in the model plant Physcomitrium patens with a partially synthetic chromosome arm. Here, we introduce the design and assembly principles of moss genome synthesis. We also discuss the remaining technical barriers in the application of synthetic genomics in seed plants.

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.

CLE peptides act via the receptor-like kinase CRINKLY 4 in Physcomitrium patens gametophore development.

The CLAVATA pathway plays a key role in the regulation of multicellular shoot and root meristems in flowering plants. In Arabidopsis, CLAVATA 3-like signaling peptides (CLEs) act via receptor-like kinases CLAVATA 1 and CRINKLY 4 (CR4). In the moss Physcomitrium patens, PpCLAVATA and PpCR4 were previously studied independently and shown to play conserved roles in the regulation of cell proliferation and differentiation. The plant calpain DEFECTIVE KERNEL 1 (DEK1) has been identified as another key regulator of cell division and cell fate in vascular plants and bryophytes. The functional interaction between CLAVATA, CR4, and DEK1 remains unknown. Here, we show that P. patens crinkly4 and dek1 mutants respond differently to CLE peptide treatments suggesting their distinct roles in the CLAVATA pathway. Reduced CLAVATA-mediated suppression of leafy shoot growth in Δcr4 mutants indicates that PpCR4 is involved in CLV3p perception, most likely as a receptor. The CLV3p strongly suppressed leaf vein development in Δcr4 mutants, suggesting that other receptors are involved in these processes and indicating a potential role of PpCR4 in organ sensitization to CLEs.

Mitochondrial respiration is essential for photosynthesis-dependent ATP supply of the plant cytosol.

Plants rely on solar energy to synthesize ATP and NADPH for photosynthetic carbon fixation and all cellular need. Mitochondrial respiration is essential in plants, but this may be due to heterotrophic bottlenecks during plant development or because it is also necessary in photosynthetically active cells. In this study, we examined in vivo changes of cytosolic ATP concentration in response to light, employing a biosensing strategy in the moss Physcomitrium patens and revealing increased cytosolic ATP concentration caused by photosynthetic activity. Plants depleted of respiratory Complex I showed decreased cytosolic ATP accumulation, highlighting a critical role of mitochondrial respiration in light-dependent ATP supply of the cytosol. Consistently, targeting mitochondrial ATP production directly, through the construction of mutants deficient in mitochondrial ATPase (complex V), led to drastic growth reduction, despite only minor alterations in photosynthetic electron transport activity. Since P. patens is photoautotrophic throughout its development, we conclude that heterotrophic bottlenecks cannot account for the indispensable role of mitochondrial respiration in plants. Instead, our results support that mitochondrial respiration is essential for ATP provision to the cytosol in photosynthesizing cells. Mitochondrial respiration provides metabolic integration, ensuring supply of cytosolic ATP essential for supporting plant growth and development.

Single-cell RNA sequencing reveals dynamics of gene expression for 2D elongation and 3D growth in Physcomitrium patens.

The transition from two-dimensional (2D) to 3D growth likely facilitated plants to colonize land, but its heterogeneity is not well understood. In this study, we utilized single-cell RNA sequencing to analyze the moss Physcomitrium patens, whose morphogenesis involves a transition from 2D to 3D growth. We profiled over 17,000 single cells covering all major vegetative tissues, including 2D filaments (chloronema and caulonema) and 3D structures (bud and gametophore). Pseudotime analyses revealed larger numbers of candidate genes that determine cell fates for 2D tip elongation or 3D bud differentiation. Using weighted gene co-expression network analysis, we identified a module that connects ß-type carbonic anhydrases (ßCAs) with auxin. We further validated the cellular expression patterns of ßCAs and demonstrated their roles in 3D gametophore development. Overall, our study provides insights into cellular heterogeneity in a moss and identifies molecular signatures that underpin the 2D-to-3D growth transition at single-cell resolution.

Illuminating the role of the calyptra in sporophyte development.

The study of moss calyptra form and function began almost 250 years ago, but calyptra research has remained a niche endeavor focusing on only a small number of species. Recent advances have focused on calyptra cuticular waxes, which function in dehydration protection of the immature sporophyte apex. The physical presence of the calyptra also plays a role in sporophyte development, potentially via its influence on auxin transport. Progress developing genomic resources for mosses beyond the model Physcomitrium patens, specifically for species with larger calyptrae and taller sporophytes, in combination with advances in CRISPR-Cas9 genome editing will enable the influence of the calyptra on gene expression and the production of RNAs and proteins that coordinate sporophyte development to be explored.

CBP60b clade proteins are prototypical transcription factors mediating immunity.

Transcriptional reprogramming is critical for plant immunity. Several calmodulin (CaM)-binding protein 60 (CBP60) family transcription factors (TFs) in Arabidopsis (Arabidopsis thaliana), including CBP60g, systemic acquired resistance deficient 1 (SARD1), CBP60a, and CBP60b, are critical for and show distinct roles in immunity. However, there are additional CBP60 members whose function is unclear. We report here that Arabidopsis CBP60c-f, 4 uncharacterized CBP60 members, play redundant roles with CBP60b in the transcriptional regulation of immunity responses, whose pCBP60b-driven expression compensates the loss of CBP60b. By contrast, neither CBP60g nor SARD1 is interchangeable with CBP60b, suggesting clade-specific functionalization. We further show that the function of CBP60b clade TFs relies on DNA-binding domains (DBDs) and CaM-binding domains, suggesting that they are downstream components of calcium signaling. Importantly, we demonstrate that CBP60s encoded in earliest land plant lineage Physcomitrium patens and Selaginella moellendorffii are functionally homologous to Arabidopsis CBP60b, suggesting that the CBP60b clade contains the prototype TFs of the CBP60 family. Furthermore, tomato and cucumber CBP60b-like genes rescue the defects of Arabidopsis cbp60b and activate the expression of tomato and cucumber SALICYLIC ACID INDUCTION DEFICIIENT2 (SID2) and ENHANCED DISEASE SUSCEPTIBILITY 1 (EDS1) genes, suggesting that immune response pathways centered on CBP60b are also evolutionarily conserved. Together, these findings suggest that CBP60b clade TFs are functionally conserved in evolution and positively mediate immunity.

NSE5 subunit interacts with distant regions of the SMC arms in the Physcomitrium patens SMC5/6 complex.

Structural maintenance of chromosome (SMC) complexes play roles in cohesion, condensation, replication, transcription, and DNA repair. Their cores are composed of SMC proteins with a unique structure consisting of an ATPase head, long arm, and hinge. SMC complexes form long rod-like structures, which can change to ring-like and elbow-bent conformations upon binding ATP, DNA, and other regulatory factors. These SMC dynamic conformational changes are involved in their loading, translocation, and DNA loop extrusion. Here, we examined the binding and role of the PpNSE5 regulatory factor of Physcomitrium patens PpSMC5/6 complex. We found that the PpNSE5 C-terminal half (aa230-505) is required for binding to its PpNSE6 partner, while the N-terminal half (aa1-230) binds PpSMC subunits. Specifically, the first 71 amino acids of PpNSE5 were required for binding to PpSMC6. Interestingly, the PpNSE5 binding required the PpSMC6 head-proximal joint region and PpSMC5 hinge-proximal arm, suggesting a long distance between binding sites on PpSMC5 and PpSMC6 arms. Therefore, we hypothesize that PpNSE5 either links two antiparallel SMC5/6 complexes or binds one SMC5/6 in elbow-bent conformation, the later model being consistent with the role of NSE5/NSE6 dimer as SMC5/6 loading factor to DNA lesions. In addition, we generated the P. patens Ppnse5KO1 mutant line with an N-terminally truncated version of PpNSE5, which exhibited DNA repair defects while keeping a normal number of rDNA repeats. As the first 71 amino acids of PpNSE5 are required for PpSMC6 binding, our results suggest the role of PpNSE5-PpSMC6 interaction in SMC5/6 loading to DNA lesions.

Blooming balloons: Searching for mechanisms of the inflated calyx.

Studying morphological novelties offers special insights into developmental biology and evolution. The inflated calyx syndrome (ICS) is a largely unrecognized but fascinating feature of flower development, where sepals form balloon-like husks that encapsulate fruits. Despite its independent emergence in many lineages of flowering plants, the genetic and molecular mechanisms of ICS remain unknown. Early studies in the Solanaceae genus Physalis put forth key roles of MADS-box genes in ICS. However, recent work suggests these classical floral identity transcription factors were false leads. With newfound capabilities that allow rapid development of genetic systems through genomics and genome editing, Physalis has re-emerged as the most tractable model species for dissecting ICS. This review revisits current understanding of ICS and highlights how recent advancements enable a reset in the search for genetic and molecular mechanisms using unbiased, systematic approaches.

The dehydration-responsive protein PpFAS1.3 in moss Physcomitrium patens plays a regulatory role in lipid metabolism.

Moss plants appear in the early stages of land colonization and possess varying degrees of dehydration tolerance. In this study, a protein called PpFAS1.3 was identified, which contains a fasciclin 1-like domain and is essential for the moss Physcomitrium patens' response to short-term rapid dehydration. When the FAS1.3 protein was knocked out, leafyshoots showed a significant decrease in tolerance to rapid dehydration, resulting in accelerated water loss and increased membrane leakage. Phylogenetic analysis suggests that PpFAS1.3 and its homologous proteins may have originated from bacteria and are specifically found in non-vascular plants like mosses and liverworts. As a dehydration-related protein, FAS1.3 plays a significant role in regulating lipid metabolism, particularly in the synthesis of free fatty acids (FFA) and the metabolism of two phospholipids, PC and PA. This discovery highlights the close connection between PpFAS1.3 and lipid metabolism, providing new insights into the molecular mechanisms underlying plant adaptation to stresses.

Moss-produced human complement factor H with modified glycans has an extended half-life and improved biological activity.

Most drugs that target the complement system are designed to inhibit the complement pathway at either the proximal or terminal levels. The use of a natural complement regulator such as factor H (FH) could provide a superior treatment option by restoring the balance of an overactive complement system while preserving its normal physiological functions. Until now, the systemic treatment of complement-associated disorders with FH has been deemed unfeasible, primarily due to high production costs, risks related to FH purified from donors' blood, and the challenging expression of recombinant FH in different host systems. We recently demonstrated that a moss-based expression system can produce high yields of properly folded, fully functional, recombinant FH. However, the half-life of the initial variant (CPV-101) was relatively short. Here we show that the same polypeptide with modified glycosylation (CPV-104) achieves a pharmacokinetic profile comparable to that of native FH derived from human serum. The treatment of FH-deficient mice with CPV-104 significantly improved important efficacy parameters such as the normalization of serum C3 levels and the rapid degradation of C3 deposits in the kidney compared to treatment with CPV-101. Furthermore, CPV-104 showed comparable functionality to serum-derived FH in vitro, as well as similar performance in ex vivo assays involving samples from patients with atypical hemolytic uremic syndrome, C3 glomerulopathy and paroxysomal nocturnal hematuria. CPV-104 - the human FH analog expressed in moss - will therefore allow the treatment of complement-associated human diseases by rebalancing instead of inhibiting the complement cascade.

Crossroads of assembling a moss genome: navigating contaminants and horizontal gene transfer in the moss Physcomitrellopsis africana.

The first chromosome-scale reference genome of the rare narrow-endemic African moss Physcomitrellopsis africana (P. africana) is presented here. Assembled from 73 × Oxford Nanopore Technologies (ONT) long reads and 163 × Beijing Genomics Institute (BGI)-seq short reads, the 414 Mb reference comprises 26 chromosomes and 22,925 protein-coding genes [Benchmarking Universal Single-Copy Ortholog (BUSCO) scores: C:94.8% (D:13.9%)]. This genome holds 2 genes that withstood rigorous filtration of microbial contaminants, have no homolog in other land plants, and are thus interpreted as resulting from 2 unique horizontal gene transfers (HGTs) from microbes. Further, P. africana shares 176 of the 273 published HGT candidates identified in Physcomitrium patens (P. patens), but lacks 98 of these, highlighting that perhaps as many as 91 genes were acquired in P. patens in the last 40 million years following its divergence from its common ancestor with P. africana. These observations suggest rather continuous gene gains via HGT followed by potential losses during the diversification of the Funariaceae. Our findings showcase both dynamic flux in plant HGTs over evolutionarily "short" timescales, alongside enduring impacts of successful integrations, like those still functionally maintained in extant P. africana. Furthermore, this study describes the informatic processes employed to distinguish contaminants from candidate HGT events.

Conquering new grounds: plant organellar C-to-U RNA editing factors can be functional in the plant cytosol.

Plant mitochondrial and chloroplast transcripts are subject to numerous events of specific cytidine-to-uridine (C-to-U) RNA editing to correct genetic information. Key protein factors for this process are specific RNA-binding pentatricopeptide repeat (PPR) proteins, which are encoded in the nucleus and post-translationally imported into the two endosymbiotic organelles. Despite hundreds of C-to-U editing sites in the plant organelles, no comparable editing has been found for nucleo-cytosolic mRNAs raising the question why plant RNA editing is restricted to chloroplasts and mitochondria. Here, we addressed this issue in the model moss Physcomitrium patens, where all PPR-type RNA editing factors comprise specific RNA-binding and cytidine deamination functionalities in single proteins. To explore whether organelle-type RNA editing can principally also take place in the plant cytosol, we expressed PPR56, PPR65 and PPR78, three editing factors recently shown to also function in a bacterial setup, together with cytosolic co-transcribed native targets in Physcomitrium. While we obtained unsatisfying results upon their constitutive expression, we found strong cytosolic RNA editing under hormone-inducible expression. Moreover, RNA-Seq analyses revealed varying numbers of up to more than 900 off-targets in other cytosolic transcripts. We conclude that PPR-mediated C-to-U RNA editing is not per se incompatible with the plant cytosol but that its limited target specificity has restricted its occurrence to the much less complex transcriptomes of mitochondria and chloroplast in the course of evolution.

Ethylene controls three-dimensional growth involving reduced auxin levels in the moss Physcomitrium patens.

The conquest of land by plants was concomitant with, and possibly enabled by, the evolution of three-dimensional (3D) growth. The moss Physcomitrium patens provides a model system for elucidating molecular mechanisms in the initiation of 3D growth. Here, we investigate whether the phytohormone ethylene, which is believed to have been a signal before land plant emergence, plays a role in 3D growth regulation in P. patens. We report ethylene controls 3D gametophore formation, based on results from exogenously applied ethylene and genetic manipulation of PpEIN2, which is a central component in the ethylene signaling pathway. Overexpression (OE) of PpEIN2 activates ethylene responses and leads to earlier formation of gametophores with fewer gametophores produced thereafter, phenocopying ethylene-treated wild-type. Conversely, Ppein2 knockout mutants, which are ethylene insensitive, show initially delayed gametophore formation with more gametophores produced later. Furthermore, pharmacological and biochemical analyses reveal auxin levels are decreased in the OE lines but increased in the knockout mutants. Our results suggest that evolutionarily, ethylene and auxin molecular networks were recruited to build the plant body plan in ancestral land plants. This might have played a role in enabling ancient plants to acclimate to the continental surfaces of the planet.

The SpRY Cas9 variant release the PAM sequence constraint for genome editing in the model plant Physcomitrium patens.

Genome editing via CRISPR/Cas has enabled targeted genetic modifications in various species, including plants. The requirement for specific protospacer-adjacent motifs (PAMs) near the target gene, as seen with Cas nucleases like SpCas9, limits its application. PAMless SpCas9 variants, designed with a relaxed PAM requirement, have widened targeting options. However, these so-call PAMless SpCas9 still show variation of editing efficiency depending on the PAM and their efficiency lags behind the native SpCas9. Here we assess the potential of a PAMless SpCas9 variant for genome editing in the model plant Physcomitrium patens. For this purpose, we developed a SpRYCas9i variant, where expression was optimized, and tested its editing efficiency using the APT as a reporter gene. We show that the near PAMless SpRYCas9i effectively recognizes specific PAMs in P. patens that are not or poorly recognized by the native SpCas9. Pattern of mutations found using the SpRYCas9i are similar to the ones found with the SpCas9 and we could not detect off-target activity for the sgRNAs tested in this study. These findings contribute to advancing versatile genome editing techniques in plants.

Genome sequence and cell biological toolbox of the highly regenerative, coenocytic green feather alga Bryopsis.

Green feather algae (Bryopsidales) undergo a unique life cycle in which a single cell repeatedly executes nuclear division without cytokinesis, resulting in the development of a thallus (>100 mm) with characteristic morphology called coenocyte. Bryopsis is a representative coenocytic alga that has exceptionally high regeneration ability: extruded cytoplasm aggregates rapidly in seawater, leading to the formation of protoplasts. However, the genetic basis of the unique cell biology of Bryopsis remains poorly understood. Here, we present a high-quality assembly and annotation of the nuclear genome of Bryopsis sp. (90.7 Mbp, 27 contigs, N50 = 6.7 Mbp, 14 034 protein-coding genes). Comparative genomic analyses indicate that the genes encoding BPL-1/Bryohealin, the aggregation-promoting lectin, are heavily duplicated in Bryopsis, whereas homologous genes are absent in other ulvophyceans, suggesting the basis of regeneration capability of Bryopsis. Bryopsis sp. possesses >30 kinesins but only a single myosin, which differs from other green algae that have multiple types of myosin genes. Consistent with this biased motor toolkit, we observed that the bidirectional motility of chloroplasts in the cytoplasm was dependent on microtubules but not actin in Bryopsis sp. Most genes required for cytokinesis in plants are present in Bryopsis, including those in the SNARE or kinesin superfamily. Nevertheless, a kinesin crucial for cytokinesis initiation in plants (NACK/Kinesin-7II) is hardly expressed in the coenocytic part of the thallus, possibly underlying the lack of cytokinesis in this portion. The present genome sequence lays the foundation for experimental biology in coenocytic macroalgae.

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.