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.

A DELAY OF GERMINATION 1 (DOG1)-like protein regulates spore germination in the moss Physcomitrium patens.

DELAY OF GERMINATION 1 is a key regulator of dormancy in flowering plants before seed germination. Bryophytes develop haploid spores with an analogous function to seeds. Here, we investigate whether DOG1 function during germination is conserved between bryophytes and flowering plants and analyse the underlying mechanism of DOG1 action in the moss Physcomitrium patens. Phylogenetic and in silico expression analyses were performed to identify and characterise DOG1 domain-containing genes in P. patens. Germination assays were performed to characterise a Ppdog1-like1 mutant, and replacement with AtDOG1 was carried out. Yeast two-hybrid assays were used to test the interaction of the PpDOG1-like protein with DELLA proteins from P. patens and A. thaliana. P. patens possesses nine DOG1 domain-containing genes. The DOG1-like protein PpDOG1-L1 (Pp3c3_9650) interacts with PpDELLAa and PpDELLAb and the A. thaliana DELLA protein AtRGA in yeast. Protein truncations revealed the DOG1 domain as necessary and sufficient for interaction with PpDELLA proteins. Spores of Ppdog1-l1 mutant germinate faster than wild type, but replacement with AtDOG1 reverses this effect. Our data demonstrate a role for the PpDOG1-LIKE1 protein in moss spore germination, possibly alongside PpDELLAs. This suggests a conserved DOG1 domain function in germination, albeit with differential adaptation of regulatory networks in seed and spore germination.

JGI Plant Gene Atlas: an updateable transcriptome resource to improve functional gene descriptions across the plant kingdom.

Gene functional descriptions offer a crucial line of evidence for candidate genes underlying trait variation. Conversely, plant responses to environmental cues represent important resources to decipher gene function and subsequently provide molecular targets for plant improvement through gene editing. However, biological roles of large proportions of genes across the plant phylogeny are poorly annotated. Here we describe the Joint Genome Institute (JGI) Plant Gene Atlas, an updateable data resource consisting of transcript abundance assays spanning 18 diverse species. To integrate across these diverse genotypes, we analyzed expression profiles, built gene clusters that exhibited tissue/condition specific expression, and tested for transcriptional response to environmental queues. We discovered extensive phylogenetically constrained and condition-specific expression profiles for genes without any previously documented functional annotation. Such conserved expression patterns and tightly co-expressed gene clusters let us assign expression derived additional biological information to 64 495 genes with otherwise unknown functions. The ever-expanding Gene Atlas resource is available at JGI Plant Gene Atlas (https://plantgeneatlas.jgi.doe.gov) and Phytozome (https://phytozome.jgi.doe.gov/), providing bulk access to data and user-specified queries of gene sets. Combined, these web interfaces let users access differentially expressed genes, track orthologs across the Gene Atlas plants, graphically represent co-expressed genes, and visualize gene ontology and pathway enrichments.

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.

Improved prime editing allows for routine predictable gene editing in Physcomitrium patens.

Efficient and precise gene editing is the gold standard of any reverse genetic study. The recently developed prime editing approach, a modified CRISPR/Cas9 [clustered regularly interspaced palindromic repeats (CRISPR)/CRISPR-associated protein] editing method, has reached the precision goal but its editing rate can be improved. We present an improved methodology that allows for routine prime editing in the model plant Physcomitrium patens, whilst exploring potential new prime editing improvements. Using a standardized protoplast transfection procedure, multiple prime editing guide RNA (pegRNA) structural and prime editor variants were evaluated targeting the APT reporter gene through direct plant selection. Together, enhancements of expression of the prime editor, modifications of the 3' extension of the pegRNA, and the addition of synonymous mutation in the reverse transcriptase template sequence of the pegRNA dramatically improve the editing rate without affecting the quality of the edits. Furthermore, we show that prime editing is amenable to edit a gene of interest through indirect selection, as demonstrated by the generation of a Ppdek10 mutant. Additionally, we determine that a plant retrotransposon reverse transcriptase enables prime editing. Finally, we show for the first time the possibility of performing prime editing with two independently coded peptides.

A vertically transmitted amalgavirus is present in certain accessions of the bryophyte Physcomitrium patens.

In the last few years, next-generation sequencing techniques have started to be used to identify new viruses infecting plants. This has allowed to rapidly increase our knowledge on viruses other than those causing symptoms in economically important crops. Here we used this approach to identify a virus infecting Physcomitrium patens that has the typical structure of the double-stranded RNA endogenous viruses of the Amalgaviridae family, which we named Physcomitrium patens amalgavirus 1, or PHPAV1. PHPAV1 is present only in certain accessions of P. patens, where its RNA can be detected throughout the cell cycle of the plant. Our analysis demonstrates that PHPAV1 can be vertically transmitted through both paternal and maternal germlines, in crosses between accessions that contain the virus with accessions that do not contain it. This work suggests that PHPAV1 can replicate in genomic backgrounds different from those that actually contain the virus and opens the door for future studies on virus-host coevolution.

PEATmoss (Physcomitrella Expression Atlas Tool): a unified gene expression atlas for the model plant Physcomitrella patens.

Physcomitrella patens is a bryophyte model plant that is often used to study plant evolution and development. Its resources are of great importance for comparative genomics and evo-devo approaches. However, expression data from Physcomitrella patens were so far generated using different gene annotation versions and three different platforms: CombiMatrix and NimbleGen expression microarrays and RNA sequencing. The currently available P. patens expression data are distributed across three tools with different visualization methods to access the data. Here, we introduce an interactive expression atlas, Physcomitrella Expression Atlas Tool (PEATmoss), that unifies publicly available expression data for P. patens and provides multiple visualization methods to query the data in a single web-based tool. Moreover, PEATmoss includes 35 expression experiments not previously available in any other expression atlas. To facilitate gene expression queries across different gene annotation versions, and to access P. patens annotations and related resources, a lookup database and web tool linked to PEATmoss was implemented. PEATmoss can be accessed at https://peatmoss.online.uni-marburg.de.

DEK1 displays a strong subcellular polarity during Physcomitrella patens 3D growth.

Defective Kernel 1 (DEK1) is genetically at the nexus of the 3D morphogenesis of land plants. We aimed to localize DEK1 in the moss Physcomitrella patens to decipher its function during this process. To detect DEK1 in vivo, we inserted the tdTomato fluorophore into PpDEK1 gene locus. Confocal microscopy coupled with the use of time-gating allowed the precise DEK1 subcellular localization during 3D morphogenesis. DEK1 localization displays a strong polarized signal, as it is restricted to the plasma membrane domain between recently divided cells during the early steps of 3D growth development as well as during the subsequent vegetative growth. The signal furthermore displays a clear developmental pattern because it is only detectable in recently divided and elongating cells. Additionally, DEK1 localization appears to be independent of its calpain domain proteolytic activity. The DEK1 polar subcellular distribution in 3D tissue developing cells defines a functional cellular framework to explain its role in this developmental phase. Also, the observation of DEK1 during spermatogenesis suggests another biological function for this protein in plants. Finally the DEK1-tagged strain generated here provides a biological platform upon which further investigations into 3D developmental processes can be performed.

Characterisation of evolutionarily conserved key players affecting eukaryotic flagellar motility and fertility using a moss model.

Defects in flagella/cilia are often associated with infertility and disease. Motile male gametes (sperm cells) are an ancestral eukaryotic trait that has been lost in several lineages like flowering plants. Here, we made use of a phenotypic male fertility difference between two moss (Physcomitrella patens) ecotypes to explore spermatozoid function. We compare genetic and epigenetic variation as well as expression profiles between the Gransden and Reute ecotype to identify a set of candidate genes associated with moss male infertility. We generated a loss-of-function mutant of a coiled-coil domain containing 39 (ccdc39) gene that is part of the flagellar hydin network. Defects in mammal and algal homologues of this gene coincide with a loss of fertility, demonstrating the evolutionary conservation of flagellar function related to male fertility across kingdoms. The Ppccdc39 mutant resembles the Gransden phenotype in terms of male fertility. Potentially, several somatic (epi-)mutations occurred during prolonged vegetative propagation of Gransden, causing regulatory differences of for example the homeodomain transcription factor BELL1. Probably these somatic changes are causative for the observed male fertility defect. We propose that moss spermatozoids might be employed as an easily accessible system to study male infertility of humans and animals in terms of flagellar structure and movement.

The Physcomitrella patens gene atlas project: large-scale RNA-seq based expression data.

High-throughput RNA sequencing (RNA-seq) has recently become the method of choice to define and analyze transcriptomes. For the model moss Physcomitrella patens, although this method has been used to help analyze specific perturbations, no overall reference dataset has yet been established. In the framework of the Gene Atlas project, the Joint Genome Institute selected P. patens as a flagship genome, opening the way to generate the first comprehensive transcriptome dataset for this moss. The first round of sequencing described here is composed of 99 independent libraries spanning 34 different developmental stages and conditions. Upon dataset quality control and processing through read mapping, 28 509 of the 34 361 v3.3 gene models (83%) were detected to be expressed across the samples. Differentially expressed genes (DEGs) were calculated across the dataset to permit perturbation comparisons between conditions. The analysis of the three most distinct and abundant P. patens growth stages - protonema, gametophore and sporophyte - allowed us to define both general transcriptional patterns and stage-specific transcripts. As an example of variation of physico-chemical growth conditions, we detail here the impact of ammonium supplementation under standard growth conditions on the protonemal transcriptome. Finally, the cooperative nature of this project allowed us to analyze inter-laboratory variation, as 13 different laboratories around the world provided samples. We compare differences in the replication of experiments in a single laboratory and between different laboratories.

The Physcomitrella patens chromosome-scale assembly reveals moss genome structure and evolution.

The draft genome of the moss model, Physcomitrella patens, comprised approximately 2000 unordered scaffolds. In order to enable analyses of genome structure and evolution we generated a chromosome-scale genome assembly using genetic linkage as well as (end) sequencing of long DNA fragments. We find that 57% of the genome comprises transposable elements (TEs), some of which may be actively transposing during the life cycle. Unlike in flowering plant genomes, gene- and TE-rich regions show an overall even distribution along the chromosomes. However, the chromosomes are mono-centric with peaks of a class of Copia elements potentially coinciding with centromeres. Gene body methylation is evident in 5.7% of the protein-coding genes, typically coinciding with low GC and low expression. Some giant virus insertions are transcriptionally active and might protect gametes from viral infection via siRNA mediated silencing. Structure-based detection methods show that the genome evolved via two rounds of whole genome duplications (WGDs), apparently common in mosses but not in liverworts and hornworts. Several hundred genes are present in colinear regions conserved since the last common ancestor of plants. These syntenic regions are enriched for functions related to plant-specific cell growth and tissue organization. The P. patens genome lacks the TE-rich pericentromeric and gene-rich distal regions typical for most flowering plant genomes. More non-seed plant genomes are needed to unravel how plant genomes evolve, and to understand whether the P. patens genome structure is typical for mosses or bryophytes.

Sporophyte Formation and Life Cycle Completion in Moss Requires Heterotrimeric G-Proteins.

In this study, we report the functional characterization of heterotrimeric G-proteins from a nonvascular plant, the moss Physcomitrella patens. In plants, G-proteins have been characterized from only a few angiosperms to date, where their involvement has been shown during regulation of multiple signaling and developmental pathways affecting overall plant fitness. In addition to its unparalleled evolutionary position in the plant lineages, the P. patens genome also codes for a unique assortment of G-protein components, which includes two copies of Gß and Gγ genes, but no canonical Gα Instead, a single gene encoding an extra-large Gα (XLG) protein exists in the P. patens genome. Here, we demonstrate that in P. patens the canonical Gα is biochemically and functionally replaced by an XLG protein, which works in the same genetic pathway as one of the Gß proteins to control its development. Furthermore, the specific G-protein subunits in P. patens are essential for its life cycle completion. Deletion of the genomic locus of PpXLG or PpGß2 results in smaller, slower growing gametophores. Normal reproductive structures develop on these gametophores, but they are unable to form any sporophyte, the only diploid stage in the moss life cycle. Finally, the mutant phenotypes of ΔPpXLG and ΔPpGß2 can be complemented by the homologous genes from Arabidopsis, AtXLG2 and AtAGB1, respectively, suggesting an overall conservation of their function throughout the plant evolution.

The DEK1 Calpain Linker Functions in Three-Dimensional Body Patterning in Physcomitrella patens.

The DEFECTIVE KERNEL1 (DEK1) calpain is a conserved 240-kD key regulator of three-dimensional body patterning in land plants acting via mitotic cell plane positioning. The activity of the cytosolic C-terminal calpain protease is regulated by the membrane-anchored DEK1 MEM, which is connected to the calpain via the 600-amino acid residue Linker. Similar to the calpain and MEM domains, the Linker is highly conserved in the land plant lineage, the similarity dropping sharply compared with orthologous charophyte sequences. Using site-directed mutagenesis, we studied the effect on Physcomitrella patens development by deleting the Linker and two conserved Linker motifs. The results show that removal of the Linker has nearly the same effect as removal of the entire DEK1 gene. In contrast, deletion of the conserved Laminin_G3 (LG3) domain had a milder effect, perturbing leafy gametophore patterning and archegonia development. The LG3 domain from Marchantia polymorpha is fully functional in P. patens, whereas angiosperm sequences are not functional. Deletion of a C-terminal Linker subsegment containing a potential calpain autolytic site severely disturbs gametophore development. Finally, changing one of the three calpain active-site amino acid residues results in the same phenotype as deleting the entire DEK1 gene. Based on the conserved nature of animal and DEK1 calpains, we propose that the DEK1 MEM-Linker complex inactivates the calpain by forcing apart the two calpain subunits carrying the three amino acids of the active site.

The phenotype of the CRINKLY4 deletion mutant of Physcomitrella patens suggests a broad role in developmental regulation in early land plants.

MAIN CONCLUSIONS: Deletion of the ancestral gene of the land plant multigene family of receptor like kinase CR4 in Physcomitrella patens demonstrates involvement in developmental control of gametophytic and sporophytic organs. The CRINKLY4 (CR4) family of receptor kinases in angiosperms consists of three clades, one including CR4, the CR4-related CCR1 and CCR2, a second including CCR3 and CCR4 family members, and a third and more distant clade. In addition to crinkly leaves in maize, which gave rise to the mutant gene name, CR4 is implicated in ovule, embryo, flower and root development in Arabidopsis thaliana. In root tips of the same species the module including a CLAVATA3/ESR-related protein, an Arabidopsis CR4, a CLAVATA1 and a WUSCHEL-related homeobox 5 (CLE40-ACR4-CLV1-WOX5) is implicated in meristem cell regulation. In embryos and shoots, CR4 acts together with A. thaliana MERISTEM LAYER 1 and PROTODERMAL FACTOR 2 to promote A. thaliana epidermis differentiation. Phylogenetic analysis has demonstrated that early land plants, e.g. mosses carry a single ancestral CR4 gene, together with genes encoding the other members of the CLE40-ACR4-CLV1-WOX5 signaling module. Here we show that CR4 serves as a broad regulator of morphogenesis both in gametophyte phyllids, archegonia and in sporophyte epidermis of the moss Physcomitrella patens. The phenotype of the CR4 deletion mutant in moss provides insight into the role of the ancestral CR4 gene as a regulator of development in early land plants.

Prenylation is required for polar cell elongation, cell adhesion, and differentiation in Physcomitrella patens.

Protein prenylation is required for a variety of growth and developmental processes in flowering plants. Here we report the consequences of loss of function of all known prenylation subunits in the moss Physcomitrella patens. As in Arabidopsis, protein farnesyltransferase and protein geranylgeranyltransferase type I are not required for viability. However, protein geranylgeranyltransferase type I activity is required for cell adhesion, polar cell elongation, and cell differentiation. Loss of protein geranylgeranyltransferase activity results in colonies of round, single-celled organisms that resemble unicellular algae. The loss of protein farnesylation is not as severe but also results in polar cell elongation and differentiation defects. The complete loss of Rab geranylgeranyltransferase activity appears to be lethal in P. patens. Labeling with antibodies to cell wall components support the lack of polarity establishment and the undifferentiated state of geranylgeranyltransferase type I mutant plants. Our results show that prenylated proteins play key roles in P. patens development and differentiation processes.

Genetic analysis of DEFECTIVE KERNEL1 loop function in three-dimensional body patterning in Physcomitrella patens.

DEFECTIVE KERNEL1 (DEK1) of higher plants plays an essential role in position-dependent signaling and consists of a large transmembrane domain (MEM) linked to a protease catalytic domain and a regulatory domain. Here, we show that the postulated sensory Loop of the MEM domain plays an important role in the developmental regulation of DEK1 activity in the moss Physcomitrella patens. Compared with P. patens lacking DEK1 (∆dek1), the dek1∆loop mutant correctly positions the division plane in the bud apical cell. In contrast with an early developmental arrest of ∆dek1 buds, dek1∆loop develops aberrant gametophores lacking expanded phyllids resulting from misregulation of mitotic activity. In contrast with the highly conserved sequence of the protease catalytic domain, the Loop is highly variable in land plants. Functionally, the sequence from Marchantia polymorpha fully complements the dek1∆loop phenotype, whereas sequences from maize (Zea mays) and Arabidopsis (Arabidopsis thaliana) give phenotypes with retarded growth and affected phyllid development. Bioinformatic analysis identifies MEM as a member of the Major Facilitator Superfamily, membrane transporters reacting to stimuli from the external environment. Transcriptome analysis comparing wild-type and ∆dek1 tissues identifies an effect on two groups of transcripts connected to dek1 mutant phenotypes: transcripts related to cell wall remodeling and regulation of the AINTEGUMENTA, PLETHORA, and BABY BOOM2 (APB2) and APB3 transcription factors known to regulate bud initiation. Finally, sequence data support the hypothesis that the advanced charophyte algae that evolved into ancestral land plants lost cytosolic calpains, retaining DEK1 as the sole calpain in the evolving land plant lineage.

The catalytic domain CysPc of the DEK1 calpain is functionally conserved in land plants.

DEK1, the single calpain of land plants, is a member of the ancient membrane bound TML-CysPc-C2L calpain family that dates back 1.5 billion years. Here we show that the CysPc-C2L domains of land plant calpains form a separate sub-clade in the DEK1 clade of the phylogenetic tree of plants. The charophycean alga Mesostigma viride DEK1-like gene is clearly divergent from those in land plants, suggesting that a major evolutionary shift in DEK1 occurred during the transition to land plants. Based on genetic complementation of the Arabidopsis thaliana dek1-3 mutant using CysPc-C2L domains of various origins, we show that these two domains have been functionally conserved within land plants for at least 450 million years. This conclusion is based on the observation that the CysPc-C2L domains of DEK1 from the moss Physcomitrella patens complements the A. thaliana dek1-3 mutant phenotype. In contrast, neither the CysPc-C2L domains from M. viride nor chimeric animal-plant calpains complement this mutant. Co-evolution analysis identified differences in the interactions between the CysPc-C2L residues of DEK1 and classical calpains, supporting the view that the two enzymes are regulated by fundamentally different mechanisms. Using the A. thaliana dek1-3 complementation assay, we show that four conserved amino acid residues of two Ca²âº-binding sites in the CysPc domain of classical calpains are conserved in land plants and functionally essential in A. thaliana DEK1.

Defective Kernel 1 (DEK1) is required for three-dimensional growth in Physcomitrella patens.

Orientation of cell division is critical for plant morphogenesis. This is evident in the formation and function of meristems and for morphogenetic transitions. Mosses undergo such transitions: from two-dimensional tip-growing filaments (protonema) to the generation of three-dimensional leaf-like structures (gametophores). The Defective Kernel 1 (DEK1) protein plays a key role in the perception of and/or response to positional cues that specify the formation and function of the epidermal layer in developing seeds of flowering plants. The moss Physcomitrella patens contains the highly conserved DEK1 gene. Using efficient gene targeting, we generated a precise PpDEK1 deletion (∆dek1), which resulted in normal filamentous growth of protonema. Two distinct mutant phenotypes were observed: an excess of buds on the protonema, and abnormal cell divisions in the emerging buds resulting in developmental arrest and the absence of three-dimensional growth. Overexpression of a complete PpDEK1 cDNA, or the calpain domain of PpDEK1 alone, successfully complements both phenotypes. These results in P. patens demonstrate the morphogenetic importance of the DEK1 protein in the control of oriented cell divisions. As it is not for protonema, it will allow dissection of the structure/function relationships of the different domains of DEK1 using gene targeting in null mutant background.

Regulation of developmental gatekeeping and cell fate transition by the calpain protease DEK1 in Physcomitrium patens.

Calpains are cysteine proteases that control cell fate transitions whose loss of function causes severe, pleiotropic phenotypes in eukaryotes. Although mainly considered as modulatory proteases, human calpain targets are directed to the N-end rule degradation pathway. Several such targets are transcription factors, hinting at a gene-regulatory role. Here, we analyze the gene-regulatory networks of the moss Physcomitrium patens and characterize the regulons that are misregulated in mutants of the calpain DEFECTIVE KERNEL1 (DEK1). Predicted cleavage patterns of the regulatory hierarchies in five DEK1-controlled subnetworks are consistent with a pleiotropic and regulatory role during cell fate transitions targeting multiple functions. Network structure suggests DEK1-gated sequential transitions between cell fates in 2D-to-3D development. Our method combines comprehensive phenotyping, transcriptomics and data science to dissect phenotypic traits, and our model explains the protease function as a switch gatekeeping cell fate transitions potentially also beyond plant development.