Efficient and modular CRISPR-Cas9 vector system for Physcomitrella patens

Comparative mutant analyses reveal a novel mechanism of ARF regulation in land plants

The plant hormone auxin regulates a wide variety of transcriptional responses depending on the cell type, environment and species. How this diversity is achieved may be related to the specific complement of auxin-signalling components in each cell. The levels of activators (class-A AUXIN RESPONSE FACTORS) and repressors (class-B ARFs) are particularly important. Tight regulation of ARF protein levels is probably key in determining this balance. Through comparative analysis of novel, dominant mutants in maize and the moss Physcomitrium patens, we have discovered a ~500-million-year-old mechanism of class-B ARF protein-level regulation mediated by proteasome degradation, important in determining cell fate decisions across land plants. Thus, our results add a key piece to the puzzle of how auxin regulates plant development.

The chloroplast-located HKT transporter plays an important role in fertilization and development in Physcomitrium patens

Cell survival depends on the maintenance of cell homeostasis that involves all the biochemical, genomic and transport processes that take place in all the organelles within a eukaryote cell. In particular, ion homeostasis is required to regulate the membrane potential and solute transport across all membranes, any alteration in these parameters will reflect in the malfunctioning of any organelle, and consequently, in the development of the organism. In plant cells, sodium transporters play a central role in keeping the concentrations of this cation across all membranes under physiological conditions to prevent its toxic effects. HKT transporters are a family of membrane proteins exclusively present in plants, with some homologs present in prokaryotes. HKT transporters have been associated to salt tolerance in plants, retrieving any leak of the cation into the xylem, or removing it from aerial parts, including the flowers, to be transported to the roots along the phloem. This function has been assigned as most of the HKT transporters are located at the plasma membrane. Here, we report the localization of the HKT from Physcomitrium patens to the thylakoid membrane, reminiscent of the prokaryote origin of these family of transporters. Mutation of PpHKT leads to several alterations in the phenotype of the organism, including the lack of sporophyte formation, and changes in expression of many genes. These alterations suggest that the breakdown in chloroplast ion homeostasis triggers a signalling cascade to the nucleus to communicate its status, being important for the moss to complete its life cycle.

Cleavable donor-assisted CRISPR/Cas9 system significantly improves the efficiency of large DNA insertion in Physcomitrium patens

Precise insertion of desired fragments can be achieved by CRISPR/Cas9-based genome editing. However, a decrease in knock-in efficiency has been observed with increasing length of exogenous inserts. In this study, we developed an in vivo cleavable (IVC) donor-assisted CRISPR/Cas9 system to improve efficiency, particularly for larger inserts, in the moss Physcomitrium patens (P. patens). The IVC donor, which contains two Cas9 nuclease recognition sites flanking the homology template, enables the in vivo release of the linear template for homology-directed repair (HDR) when co-delivered with the corresponding CRISPR/Cas9 plasmid into protoplasts. In our experimental framework, two distinct sgRNAs and four different DNA inserts were evaluated. Compared with standard circular donors, IVC donors significantly enhanced the efficiency of CRISPR/Cas9-mediated precise insertion of 5.8, 7.5, and 11.1 kb DNA fragments at the PpPDV2-4 sgRNA target site, improving integration rates from 29.6 to 67.8%, from 15.0 to 72.0%, and from 12.1 to 65.6%, respectively. At an alternative sgRNA2 target site within the Pp6c18_3160 locus, the IVC donor also demonstrated a higher knock-in efficiency for a 7.4 kb fragment compared with the standard circular donor. This IVC donor-assisted CRISPR/Cas9 approach for large fragment knock-in represents a powerful tool for basic research and synthetic biology efforts in moss species. Moreover, this strategy may be potentially applicable to crops that are amenable to protoplast transformation and regeneration, facilitating the improvement of key traits.

An alternate route for cellulose microfibril biosynthesis in plants

Similar to cellulose synthases (CESAs), cellulose synthase-like D (CSLD) proteins synthesize β-1,4-glucan in plants. CSLDs are important for tip growth and cytokinesis, but it was unknown whether they form membrane complexes in vivo or produce microfibrillar cellulose. We produced viable CESA-deficient mutants of the moss Physcomitrium patens to investigate CSLD function without interfering CESA activity. Microscopy and spectroscopy showed that CESA-deficient mutants synthesize cellulose microfibrils that are indistinguishable from those in vascular plants. Correspondingly, freeze-fracture electron microscopy revealed rosette-shaped particle assemblies in the plasma membrane that are indistinguishable from CESA-containing rosette cellulose synthesis complexes (CSCs). Our data show that proteins other than CESAs, most likely CSLDs, produce cellulose microfibrils in P. patens protonemal filaments. The data suggest that the specialized roles of CSLDs in cytokinesis and tip growth are based on differential expression and different interactions with microtubules and possibly Ca2+, rather than structural differences in the microfibrils they produce.

Establishing CRISPR-Cas9 in the sexually dimorphic moss, Ceratodon purpureus

The development of CRISPR technologies provides a powerful tool for understanding the evolution and functionality of essential biological processes. Here we demonstrate successful CRISPR-Cas9 genome editing in the dioecious moss species, Ceratodon purpureus. Using an existing selection system from the distantly related hermaphroditic moss, Physcomitrium patens, we generated knock-outs of the APT reporter gene by employing CRISPR-targeted mutagenesis under expression of native U6 snRNA promoters. Next, we used the native homology-directed repair (HDR) pathway, combined with CRISPR-Cas9, to knock in two reporter genes under expression of an endogenous RPS5A promoter in a newly developed landing site in C. purpureus. Our results show that the molecular tools developed in P. patens can be extended to other mosses across this ecologically important and developmentally variable group. These findings pave the way for precise and powerful experiments aimed at identifying the genetic basis of key functional variation within the bryophytes and between the bryophytes and other land plants.

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.

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.

Comparative mutant analyses reveal a novel mechanism of ARF regulation in land plants

A major challenge in plant biology is to understand how the plant hormone auxin regulates diverse transcriptional responses throughout development, in different environments, and in different species. The answer may lie in the specific complement of auxin signaling components in each cell. The balance between activators (class-A AUXIN RESPONSE FACTORS) and repressors (class-B ARFs) is particularly important. It is unclear how this balance is achieved. Through comparative analysis of novel, dominant mutants in maize and the moss Physcomitrium patens , we have discovered a ∼500-million-year-old mechanism of class-B ARF protein level regulation, important in determining cell fate decisions across land plants. Thus, our results add a key piece to the puzzle of how auxin regulates plant development.

Characterization of the conserved features of the NSE6 subunit of the Physcomitrium patens SMC5/6 complex

Structural maintenance of chromosomes (SMC) complexes are molecular machines ensuring chromatin organization at higher levels. They play direct roles in cohesion, condensation, replication, transcription, and DNA repair. Their cores are composed of long-armed SMC, kleisin, and kleisin-associated subunits. Additional factors, like NSE6 within SMC5/6, bind to SMC core complexes and regulate their activities. In the human HsNSE6/SLF2, we recently identified a new CANIN domain. Here we tracked down its sequence homology to lower plants, selected the bryophyte Physcomitrium patens, and analyzed PpNSE6 protein-protein interactions to explore its conservation in detail. We identified a previously unrecognized core sequence motif conserved from yeasts to humans within the NSE6 CANIN domain. This motif mediates the interaction between NSE6 and its NSE5 partner in yeasts and plants. In addition, the CANIN domain and its preceding PpNSE6 sequences bind both PpSMC5 and PpSMC6 arms. Interestingly, we mapped the PpNSE6-binding site at the PpSMC5 arm right next to the PpNSE2-binding surface. The position of NSE6 at SMC arms suggests its role in the regulation of SMC5/6 dynamics. Consistent with the regulatory role of NSE6 subunits, Ppnse6 mutant lines were viable and sensitive to the DNA-damaging drug bleomycin and lost a large portion of rDNA copies. These moss mutants also exhibited reduced growth and developmental aberrations. Altogether, our data showed the conserved function of the NSE6 subunit and architecture of the SMC5/6 complex across species.

Cellulose synthase-like D movement in the plasma membrane requires enzymatic activity

Cellulose Synthase-Like D (CSLD) proteins, important for tip growth and cell division, are known to generate β-1,4-glucan. However, whether they are propelled in the membrane as the glucan chains they produce assemble into microfibrils is unknown. To address this, we endogenously tagged all eight CSLDs in Physcomitrium patens and discovered that they all localize to the apex of tip-growing cells and to the cell plate during cytokinesis. Actin is required to target CSLD to cell tips concomitant with cell expansion, but not to cell plates, which depend on actin and CSLD for structural support. Like Cellulose Synthase (CESA), CSLD requires catalytic activity to move in the plasma membrane. We discovered that CSLD moves significantly faster, with shorter duration and less linear trajectories than CESA. In contrast to CESA, CSLD movement was insensitive to the cellulose synthesis inhibitor isoxaben, suggesting that CSLD and CESA function within different complexes possibly producing structurally distinct cellulose microfibrils.

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.

A conserved module in the formation of moss midribs and seed plant axillary meristems

Different evolutionary lineages have evolved distinct characteristic body plans and anatomical structures, but their origins are largely elusive. For example, seed plants evolve axillary meristems to enable lateral branching. In moss, the phyllid (leaf) midrib containing specialized cells is responsible for water conduction and support. Midribs function like vascular tissues in flowering plants but may have risen from a different evolutionary path. Here, we demonstrate that midrib formation in the model moss Physcomitrium patens is regulated by orthologs of Arabidopsis LATERAL SUPPRESSOR (LAS), a key regulator of axillary meristem initiation. Midribs are missing in loss-of-function mutants, and ectopic formation of midrib-like structures is induced in overexpression lines. Furthermore, the PpLAS/AtLAS genes have conserved functions in the promotion of cell division in both lineages, which alleviates phenotypes in both Physcomitrium and Arabidopsis las mutants. Our results show that a conserved regulatory module is reused in divergent developmental programs, water-conducting and supporting tissues in moss, and axillary meristem initiation in seed plants.

MS1/MMD1 homologues in the moss Physcomitrium patens are required for male and female gametogenesis

The Arabidopsis Plant HomeoDomain (PHD) proteins AtMS1 and AtMMD1 provide chromatin-mediated transcriptional regulation essential for tapetum-dependent pollen formation. This pollen-based male gametogenesis is a derived trait of seed plants. Male gametogenesis in the common ancestors of land plants is instead likely to have been reminiscent of that in extant bryophytes where flagellated sperms are produced by an elaborate gametophyte generation. Still, also bryophytes possess MS1/MMD1-related PHD proteins. We addressed the function of two MS1/MMD1-homologues in the bryophyte model moss Physcomitrium patens by the generation and analysis of reporter and loss-of-function lines. The two genes are together essential for both male and female fertility by providing functions in the gamete-producing inner cells of antheridia and archegonia. They are furthermore expressed in the diploid sporophyte generation suggesting a function during sporogenesis, a process proposed related by descent to pollen formation in angiosperms. We propose that the moss MS1/MMD1-related regulatory network required for completion of male and female gametogenesis, and possibly for sporogenesis, represent a heritage from ancestral land plants.

Dependence on clade II bHLH transcription factors for nursing of haploid products by tapetal-like cells is conserved between moss sporangia and angiosperm anthers

Clade II basic helix-loop-helix transcription factors (bHLH TFs) are essential for pollen production and tapetal nursing functions in angiosperm anthers. As pollen has been suggested to be related to bryophyte spores by descent, we characterized two Physcomitrium (Physcomitrella) patens clade II bHLH TFs (PpbHLH092 and PpbHLH098), to test if regulation of sporogenous cells and the nursing cells surrounding them is conserved between angiosperm anthers and bryophyte sporangia. We made CRISPR-Cas9 reporter and loss-of-function lines to address the function of PpbHLH092/098. We sectioned and analyzed WT and mutant sporophytes for a comprehensive stage-by-stage comparison of sporangium development. Spore precursors in the P. patens sporangium are surrounded by nursing cells showing striking similarities to tapetal cells in angiosperms. Moss clade II bHLH TFs are essential for the differentiation of these tapetal-like cells and for the production of functional spores. Clade II bHLH TFs provide a conserved role in controlling the sporophytic somatic cells surrounding and nursing the sporogenous cells in both moss sporangia and angiosperm anthers. This supports the hypothesis that such nursing functions in mosses and angiosperms, lineages separated by c. 450 million years, are related by descent.

Genome-wide specificity of plant genome editing by both CRISPR-Cas9 and TALEN

CRISPR and TALENs are efficient systems for gene editing in many organisms including plants. In many cases the CRISPR–Cas or TALEN modules are expressed in the plant cell only transiently. Theoretically, transient expression of the editing modules should limit unexpected effects compared to stable transformation. However, very few studies have measured the off-target and unpredicted effects of editing strategies on the plant genome, and none of them have compared these two major editing systems. We conducted, in Physcomitrium patens, a comprehensive genome-wide investigation of off-target mutations using either a CRISPR–Cas9 or a TALEN strategy. We observed a similar number of differences for the two editing strategies compared to control non-transfected plants, with an average of 8.25 SNVs and 19.5 InDels for the CRISPR-edited plants, and an average of 17.5 SNVs and 32 InDels for the TALEN-edited plants. Interestingly, a comparable number of SNVs and InDels could be detected in the PEG-treated control plants. This shows that except for the on-target modifications, the gene editing tools used in this study did not show a significant off-target activity nor unpredicted effects on the genome, and did not lead to transgene integration. The PEG treatment, a well-established biotechnological method, in itself, was the main source of mutations found in the edited plants.

Cellulose synthesis complexes are homo-oligomeric and hetero-oligomeric in Physcomitrium patens

The common ancestor of seed plants and mosses contained homo-oligomeric cellulose synthesis complexes (CSCs) composed of identical subunits encoded by a single CELLULOSE SYNTHASE (CESA) gene. Seed plants use different CESA isoforms for primary and secondary cell wall deposition. Both primary and secondary CESAs form hetero-oligomeric CSCs that assemble and function in planta only when all the required isoforms are present. The moss Physcomitrium (Physcomitrella) patens has seven CESA genes that can be grouped into two functionally and phylogenetically distinct classes. Previously, we showed that PpCESA3 and/or PpCESA8 (class A) together with PpCESA6 and/or PpCESA7 (class B) form obligate hetero-oligomeric complexes required for normal secondary cell wall deposition. Here, we show that gametophore morphogenesis requires a member of class A, PpCESA5, and is sustained in the absence of other PpCESA isoforms. PpCESA5 also differs from the other class A PpCESAs as it is able to self-interact and does not co-immunoprecipitate with other PpCESA isoforms. These results are consistent with the hypothesis that homo-oligomeric CSCs containing only PpCESA5 subunits synthesize cellulose required for gametophore morphogenesis. Analysis of mutant phenotypes also revealed that, like secondary cell wall deposition, normal protonemal tip growth requires class B isoforms (PpCESA4 or PpCESA10), along with a class A partner (PpCESA3, PpCESA5, or PpCESA8). Thus, P. patens contains both homo-oligomeric and hetero-oligomeric CSCs.

COPII Sec23 proteins form isoform-specific endoplasmic reticulum exit sites with differential effects on polarized growth

Coat Protein complex II (COPII), a coat protein complex that forms vesicles on the endoplasmic reticulum (ER), mediates trafficking to the Golgi. While metazoans have few genes encoding each COPII component, plants have expanded these gene families, leading to the hypothesis that plant COPII has functionally diversified. In the moss Physcomitrium (Physcomitrella) patens, the Sec23/24 gene families are each composed of seven genes. Silencing Sec23/24 revealed isoform-specific contributions to polarized growth, with the closely related Sec23D/E and Sec24C/D essential for protonemal development. Focusing on Sec23, we discovered that Sec23D/E mediate ER-to Golgi transport and are essential for tip growth, with Sec23D localizing to presumptive ER exit sites. In contrast, Sec23A, B, C, F, and G are dispensable and do not quantitatively affect ER-to-Golgi trafficking. However, Δsec23abcfg plants exhibited reduced secretion of plasma membrane cargo. Of the four highly expressed protonemal Sec23 genes, Sec23F/G are members of a divergent Sec23 clade specifically retained in land plants. Notably, Sec23G accumulates on ER-associated foci that are significantly larger, do not overlap with, and are independent of Sec23D. While Sec23D/E form ER exit sites and function as bona fide COPII components essential for tip-growing protonemata, Sec23G and the closely related Sec23F have likely functionally diversified, forming separate and independent ER exit sites and participating in Golgi-independent trafficking pathways.

Targeted Gene Knockouts by Protoplast Transformation in the Moss Physcomitrella patens

Targeted gene knockout is particularly useful for analyzing gene functions in plant growth, signaling, and development. By transforming knockout cassettes consisting of homologous sequences of the target gene into protoplasts, the classical gene targeting method aims to obtain targeted gene replacement, allowing for the characterization of gene functions in vivo. The moss Physcomitrella patens is a known model organism for a high frequency of homologous recombination and thus harbors a remarkable rate of gene targeting. Other moss features, including easy to culture, dominant haploidy phase, and sequenced genome, make gene targeting prevalent in Physcomitrella patens. However, even gene targeting was powerful to generate knockouts, researchers using this method still experienced technical challenges. For example, obtaining a good number of targeted knockouts after protoplast transformation and regeneration disturbed the users. Off-target mutations such as illegitimate random integration mediated by nonhomologous end joining and targeted insertion wherein one junction on-target but the other end off-target is commonly present in the knockouts. Protoplast fusion during transformation and regeneration was also a problem. This review will discuss the advantages and technical challenges of gene targeting. Recently, CRISPR-Cas9 is a revolutionary technology and becoming a hot topic in plant gene editing. In the second part of this review, CRISPR-Cas9 technology will be focused on and compared to gene targeting regarding the practical use in Physcomitrella patens. This review presents an updated perspective of the gene targeting and CRISPR-Cas9 techniques to plant biologists who may consider studying gene functions in the model organism Physcomitrella patens.

Multiplex CRISPR-Cas9 mutagenesis of the phytochrome gene family in Physcomitrium (Physcomitrella) patens

We mutated all seven Physcomitrium (Physcomitrella) patens phytochrome genes using highly-efficient CRISPR-Cas9 procedures. We thereby identified phy5a as the phytochrome primarily responsible for inhibiting gravitropism, proving the utility of the mutant library. The CRISPR-Cas9 system is a powerful tool for genome editing. Here we report highly-efficient multiplex CRISPR-Cas9 editing of the seven-member phytochrome gene family in the model bryophyte Physcomitrium (Physcomitrella) patens. Based on the co-delivery of an improved Cas9 plasmid with multiple sgRNA plasmids and an efficient screening procedure to identify high-order multiple mutants prior to sequencing, we demonstrate successful targeting of all seven PHY genes in a single transfection. We investigated further aspects of the CRISPR methodology in Physcomitrella, including the significance of spacing between paired sgRNA targets and the efficacy of NHEJ and HDR in repairing the chromosome when excising a complete locus. As proof-of-principle, we show that the septuple phy^- mutant remains gravitropic in light, in line with expectations, and on the basis of data from lower order multiplex knockouts conclude that phy5a is the principal phytochrome responsible for inhibiting gravitropism in light. We expect, therefore, that this mutant collection will be valuable for further studies of phytochrome function and that the methods we describe will allow similar approaches to revealing specific functions in other gene families.

Phytohormone biosynthesis and signaling pathways of mosses

Most known phytohormones regulate moss development. We present a comprehensive view of the synthesis and signaling pathways for the most investigated of these compounds in mosses, focusing on the model Physcomitrium patens. The last 50 years of research have shown that most of the known phytohormones are synthesized by the model moss Physcomitrium patens (formerly Physcomitrella patens) and regulate its development, in interaction with responses to biotic and abiotic stresses. Biosynthesis and signaling pathways are best described in P. patens for the three classical hormones auxins, cytokinins and abscisic acid. Furthermore, their roles in almost all steps of development, from early filament growth to gametophore development and sexual reproduction, have been the focus of much research effort over the years. Evidence of hormonal roles exist for ethylene and for CLE signaling peptides, as well as for salicylic acid, although their possible effects on development remain unclear. Production of brassinosteroids by P. patens is still debated, and modes of action for these compounds are even less known. Gibberellin biosynthesis and signaling may have been lost in P. patens, while gibberellin precursors such as ent-kaurene derivatives could be used as signals in a yet to discover pathway. As for jasmonic acid, it is not used per se as a hormone in P. patens, but its precursor OPDA appears to play a corresponding role in defense against abiotic stress. We have tried to gather a comprehensive view of the biosynthesis and signaling pathways for all these compounds in mosses, without forgetting strigolactones, the last class of plant hormones to be reported. Study of the strigolactone response in P. patens points to a novel signaling compound, the KAI2-ligand, which was likely employed as a hormone prior to land plant emergence.

Evolution of plant telomerase RNAs: farther to the past, deeper to the roots

The enormous sequence heterogeneity of telomerase RNA (TR) subunits has thus far complicated their characterization in a wider phylogenetic range. Our recent finding that land plant TRs are, similarly to known ciliate TRs, transcribed by RNA polymerase III and under the control of the type-3 promoter, allowed us to design a novel strategy to characterize TRs in early diverging Viridiplantae taxa, as well as in ciliates and other Diaphoretickes lineages. Starting with the characterization of the upstream sequence element of the type 3 promoter that is conserved in a number of small nuclear RNAs, and the expected minimum TR template region as search features, we identified candidate TRs in selected Diaphoretickes genomes. Homologous TRs were then used to build covariance models to identify TRs in more distant species. Transcripts of the identified TRs were confirmed by transcriptomic data, RT-PCR and Northern hybridization. A templating role for one of our candidates was validated in Physcomitrium patens. Analysis of secondary structure demonstrated a deep conservation of motifs (pseudoknot and template boundary element) observed in all published TRs. These results elucidate the evolution of the earliest eukaryotic TRs, linking the common origin of TRs across Diaphoretickes, and underlying evolutionary transitions in telomere repeats.

Plant PIEZO homologs modulate vacuole morphology during tip growth

In animals, PIEZOs are plasma membrane-localized cation channels involved in diverse mechanosensory processes. We investigated PIEZO function in tip-growing cells in the moss Physcomitrium patens and the flowering plant Arabidopsis thaliana PpPIEZO1 and PpPIEZO2 redundantly contribute to the normal growth, size, and cytoplasmic calcium oscillations of caulonemal cells. Both PpPIEZO1 and PpPIEZO2 localized to vacuolar membranes. Loss-of-function, gain-of-function, and overexpression mutants revealed that moss PIEZO homologs promote increased complexity of vacuolar membranes through tubulation, internalization, and/or fission. Arabidopsis PIEZO1 also localized to the tonoplast and is required for vacuole tubulation in the tips of pollen tubes. We propose that in plant cells the tonoplast has more freedom of movement than the plasma membrane, making it a more effective location for mechanosensory proteins.

Disruption of long-chain base hydroxylation alters growth and impacts sphingolipid synthesis in Physcomitrella patens

Sphingolipids have roles as membrane structural components and as bioactive molecules in plants. In Physcomitrella patens, 4-hydroxysphinganine (phytosphingosine, t18:0) is the predominant sphingolipid long-chain base (LCB). To assess the functional significance of t18:0, CRISPR-Cas9 mutagenesis was used to generate mutant lines lacking the sole SPHINGOID BASE HYDROXYLASE (SBH) gene encoding the hydroxylase responsible for converting sphinganine (d18:0) to t18:0. Total sphingolipid content in sbh protonemata was 2.4-fold higher than in wild-type. Modest changes in glycosyl inositolphosphorylceramide (GIPC) glycosylation patterns occurred. Sphingolipidomic analyses of mutants lacking t18:0 indicated modest alterations in acyl-chain pairing with d18:0 in GIPCs and ceramides, but dramatic alterations in acyl-chain pairing in glucosylceramides, in which 4,8-sphingadienine (d18:2) was the principal LCB. A striking accumulation of free and phosphorylated LCBs accompanied loss of the hydroxylase. The sbh lines exhibited altered morphology, including smaller chloronemal cell size, irregular cell shape, reduced gametophore size, and increased pigmentation. In the presence of the synthetic trihydroxy LCB t17:0, the endogenous sphingolipid content of sbh lines decreased to wild-type levels, and the mutants exhibited phenotypes more similar to wild-type plants. These results demonstrate the importance of sphingolipid content and composition to Physcomitrella growth. They also illuminate similarities in regulating sphingolipid content but differences in regulating sphingolipid species composition between the bryophyte P. patens and angiosperm A. thaliana.

A blueprint for gene function analysis through Base Editing in the model plant Physcomitrium (Physcomitrella) patens

CRISPR-Cas9 has proven to be highly valuable for genome editing in plants, including the model plant Physcomitrium patens. However, the fact that most of the editing events produced using the native Cas9 nuclease correspond to small insertions and deletions is a limitation. CRISPR-Cas9 base editors enable targeted mutation of single nucleotides in eukaryotic genomes and therefore overcome this limitation. Here, we report two programmable base-editing systems to induce precise cytosine or adenine conversions in P. patens. Using cytosine or adenine base editors, site-specific single-base mutations can be achieved with an efficiency up to 55%, without off-target mutations. Using the APT gene as a reporter of editing, we could show that both base editors can be used in simplex or multiplex, allowing for the production of protein variants with multiple amino-acid changes. Finally, we set up a co-editing selection system, named selecting modification of APRT to report gene targeting (SMART), allowing up to 90% efficiency site-specific base editing in P. patens. These two base editors will facilitate gene functional analysis in P. patens, allowing for site-specific editing of a given base through single sgRNA base editing or for in planta evolution of a given gene through the production of randomly mutagenised variants using multiple sgRNA base editing.

SABRE populates ER domains essential for cell plate maturation and cell expansion influencing cell and tissue patterning

SABRE, which is found throughout eukaryotes and was originally identified in plants, mediates cell expansion, division plane orientation, and planar polarity in plants. How and where SABRE mediates these processes remain open questions. We deleted SABRE in Physcomitrium patens, an excellent model for cell biology. SABRE null mutants were stunted, similar to phenotypes in seed plants. Additionally, polarized growing cells were delayed in cytokinesis, sometimes resulting in catastrophic failures. A functional SABRE fluorescent fusion protein localized to dynamic puncta on regions of the endoplasmic reticulum (ER) during interphase and at the cell plate during cell division. Without SABRE, cells accumulated ER aggregates and the ER abnormally buckled along the developing cell plate. Notably, callose deposition was delayed in ∆sabre, and in cells that failed to divide, abnormal callose accumulations formed at the cell plate. Our findings revealed a surprising and fundamental role for the ER in cell plate maturation.

Mutations in the Physcomitrium patens gene encoding Aminodeoxychorismate Synthase confer auxotrophic phenotypes

To facilitate genetic mapping of developmental mutants of Physcomitrium patens, we produced a genetic marker that combines recessive auxotrophy with dominant positive selection. We first identified the gene affected by the pabB4 auxotrophic mutation and then replaced it with a cassette that confers antibiotic resistance. This strain may be used to produce bi-parental somatic hybrids with nearly any other strain.

A Fully Functional ROP Fluorescent Fusion Protein Reveals Roles for This GTPase in Subcellular and Tissue-Level Patterning

Rho of Plants (ROPs) are GTPases that regulate polarity and patterned wall deposition in plants. As these small, globular proteins have many interactors, it has been difficult to ensure that methods to visualize ROP in live cells do not affect ROP function. Here, motivated by work in fission yeast (Schizosaccharomyces pombe), we generated a fluorescent moss (Physcomitrium [Physcomitrella] patens) ROP4 fusion protein by inserting mNeonGreen after Gly-134. Plants harboring tagged ROP4 and no other ROP genes were phenotypically normal. Plants lacking all four ROP genes comprised an unpatterned clump of spherical cells that were unable to form gametophores, demonstrating that ROP is essentially for spatial patterning at the cellular and tissue levels. The functional ROP fusion protein formed a steep gradient at the apical plasma membranes of growing tip cells. ROP also predicted the site of branch formation in the apical cell at the onset of mitosis, which occurs one to two cell cycles before a branch cell emerges. While fluorescence recovery after photobleaching studies demonstrated that ROP dynamics do not depend on the cytoskeleton, acute depolymerization of the cytoskeleton removed ROP from the membrane only in recently divided cells, pointing to a feedback mechanism between the cell cycle, cytoskeleton, and ROP.

High gene space divergence contrasts with frozen vegetative architecture in the moss family Funariaceae

A new paradigm has slowly emerged regarding the diversification of bryophytes, with inferences from molecular data highlighting a dynamic evolution of their genome. However, comparative studies of expressed genes among closely related taxa is so far missing. Here we contrast the dimensions of the vegetative transcriptome of Funaria hygrometrica and Physcomitrium pyriforme against the genome of their relative, Physcomitrium (Physcomitrella) patens. These three species of Funariaceae share highly conserved vegetative bodies, and are partially sympatric, growing on mineral soil in mostly temperate regions. We analyzed the vegetative gametophytic transcriptome of F. hygrometrica and P. pyriforme and mapped short reads, transcripts, and proteins to the genome and gene space of P. patens. Only about half of the transcripts of F. hygrometrica map to their ortholog in P. patens, whereas at least 90% of those of P. pyriforme align to loci in P. patens. Such divergence is unexpected given the high morphological similarity of the gametophyte but reflects the estimated times of divergence of F. hygrometrica and P. pyriforme from P. patens, namely 55 and 20 mya, respectively. The newly sampled transcriptomes bear signatures of at least one, rather ancient, whole genome duplication (WGD), which may be shared with one reported for P. patens. The transcriptomes of F. hygrometrica and P. pyriforme reveal significant contractions or expansions of different gene families. While transcriptomes offer only an incomplete estimate of the gene space, the high number of transcripts obtained suggest a significant divergence in gene sequences, and gene number among the three species, indicative of a rather strong, dynamic genome evolution, shaped in part by whole, partial or localized genome duplication. The gene ontology of their specific and rapidly-evolving protein families, suggests that the evolution of the Funariaceae may have been driven by the diversification of metabolic genes that may optimize the adaptations to environmental conditions, a hypothesis well in line with ecological patterns in the genetic diversity and structure in seed plants.

What can lycophytes teach us about plant evolution and development? Modern perspectives on an ancient lineage

All Evo-Devo studies rely on representative sampling across the tree of interest to elucidate evolutionary trajectories through time. In land plants, genetic resources are well established in model species representing lineages including bryophytes (mosses, liverworts, and hornworts), monilophytes (ferns and allies), and seed plants (gymnosperms and flowering plants), but few resources are available for lycophytes (club mosses, spike mosses, and quillworts). Living lycophytes are a sister group to the euphyllophytes (the fern and seed plant clade), and have retained several ancestral morphological traits despite divergence from a common ancestor of vascular plants around 420 million years ago. This sister relationship offers a unique opportunity to study the conservation of traits such as sporophyte branching, vasculature, and indeterminacy, as well as the convergent evolution of traits such as leaves and roots which have evolved independently in each vascular plant lineage. To elucidate the evolution of vascular development and leaf formation, molecular studies using RNA Seq, quantitative reverse transcription polymerase chain reaction, in situ hybridisation and phylogenetics have revealed the diversification and expression patterns of KNOX, ARP, HD-ZIP, KANADI, and WOX gene families in lycophytes. However, the molecular basis of further trait evolution is not known. Here we describe morphological traits of living lycophytes and their extinct relatives, consider the molecular underpinnings of trait evolution and discuss future research required in lycophytes to understand the key evolutionary innovations enabling the growth and development of all vascular plants.

The Moss Physcomitrium (Physcomitrella) patens: A Model Organism for Non-Seed Plants

Since the discovery two decades ago that transgenes are efficiently integrated into the genome of Physcomitrella patens by homologous recombination, this moss has been a premier model system to study evolutionary developmental biology questions, stem cell reprogramming, and the biology of nonvascular plants. P patens was the first non-seed plant to have its genome sequenced. With this level of genomic information, together with increasing molecular genetic tools, a large number of reverse genetic studies have propelled the use of this model system. A number of technological advances have recently opened the door to forward genetics as well as extremely efficient and precise genome editing in P patens Additionally, careful phylogenetic studies with increased resolution have suggested that P patens emerged from within Physcomitrium Thus, rather than Physcomitrella patens, the species should be named Physcomitrium patens Here we review these advances and describe the areas where P patens has had the most impact on plant biology.

In vivo analysis of formin dynamics in the moss P. patens reveals functional class diversification

Formins are actin regulators critical for diverse processes across eukaryotes. With many formins in plants and animals, it has been challenging to determine formin function in vivo We found that the phylogenetically distinct class I integral membrane formins (denoted For1) from the moss P.patens enrich at sites of membrane turnover, with For1D more tightly associated with the plasma membrane than For1A. To probe formin function, we generated formin-null lines with greatly reduced formin complexity. We found that For1A and For1D help to anchor actin near the cell apex, with For1A contributing to formation of cytosolic actin, while For1D contributes to plasma membrane-associated actin. At the cortex, For1A and For1D localized to motile puncta and differentially impacted actin dynamics. We found that class I cortical formin mobility depended on microtubules and only moderately on actin, whereas class II formin (denoted For2) mobility solely depended on actin. Moreover, cortical For2A tightly correlated with the puncta labeled by the endocytic membrane dye FM4-64, and null mutants in class I formins did not affect uptake of a similar dye, FM1-43, suggesting that class I and II formins are involved in distinct membrane trafficking pathways.

Expanding the CRISPR Toolbox in P. patens Using SpCas9-NG Variant and Application for Gene and Base Editing in Solanaceae Crops

Genome editing has become a major tool for both functional studies and plant breeding in several species. Besides generating knockouts through the classical CRISPR-Cas9 system, recent development of CRISPR base editing holds great and exciting opportunities for the production of gain-of-function mutants. The PAM requirement is a strong limitation for CRISPR technologies such as base editing, because the base substitution mainly occurs in a small edition window. As precise single amino-acid substitution can be responsible for functions associated to some domains or agronomic traits, development of Cas9 variants with relaxed PAM recognition is of upmost importance for gene function analysis and plant breeding. Recently, the SpCas9-NG variant that recognizes the NGN PAM has been successfully tested in plants, mainly in monocotyledon species. In this work, we studied the efficiency of SpCas9-NG in the model moss Physcomitrellapatens and two Solanaceae crops (Solanum lycopersicum and Solanum tuberosum) for both classical CRISPR-generated gene knock-out and cytosine base editing. We showed that the SpCas9-NG greatly expands the scope of genome editing by allowing the targeting of non-canonical NGT and NGA PAMs. The CRISPR toolbox developed in our study opens up new gene function analysis and plant breeding perspectives for model and crop plants.