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.

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.

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.

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.

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.

Prime Editing in the model plant Physcomitrium patens and its potential in the tetraploid potato

Since its discovery and first applications for genome editing in plants, the clustered regularly interspaced short palindromic repeats (CRISPR)-Cas9 technology has revolutionized plant research and precision crop breeding. Although the classical CRISPR-Cas9 system is a highly efficient tool for disruptive targeted mutagenesis, this system is mostly inefficient for the introduction of precise and predictable nucleotide substitutions. Recently, Prime Editing technology has been developed, allowing the simultaneous generation of nucleotide transitions and transversions but also short defined indels. In this study, we report on the successful use of Prime Editing in two plants of interest: the plant model Physcomitrium patens and the tetraploid and highly heterozygous potato (Solanum tuberosum). In both cases editing rates were lower than with other CRISPR-Cas9 based techniques, but we were able to successfully introduce nucleotide transversions into targeted genes, a unique feature of Prime Editing. Additionally, the analysis of potential off-target mutation sites in P. patens suggested very high targeting fidelity in this organism. The present work paves the way for the use Prime Editing in Physcomitrium patens and potato, however highlighting the limitations that need to be overcome for more efficient precision plant breeding.

Physcomitrium patens Protoplasting and Protoplast Transfection

Protoplast production with the moss Physcomitrium (Physcomitrella) patens has a long and successful history. As a tool, it has not only been the base of reverse genetic studies covering research fields as diverse as development, metabolism, or gene network regulation but also allowed its development as a bioengineering platform for protein production. We present here a standardized protocol for protoplast production from Physcomitrium (Physcomitrella) patens protonemata. Additionally, we detail procedures for their transfection, their plating for optimal regeneration, and three alternative selection approaches. To improve the consistency of protoplast regeneration, we describe a new option for protoplast embedding. The use of an alginate matrix to regenerate moss protoplast alleviates the use of warm agarized medium. Thus, it optimizes transformed protoplast survival without any morphological detrimental effect or impact on transfection efficiency.

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.

The nuclear GUCT domain-containing DEAD-box RNA helicases govern gametophytic and sporophytic development in Physcomitrium patens

KEY MESSAGE

In Physcomitrium patens, PpRH1/PpRH2 are GUCT-domain-containing DEAD-BOX RNA helicases localize to the nucleus. They are implicated in cell and tissue development in all stages of the moss life cycle.

ABSTRACT

The DEAD-box-containing RNA helicase family encompasses a large and functionally important group of enzymes involved in cellular processes committed to the metabolism of RNA, including its transcription, processing, transport, translation and decay. Studies indicate this protein family has implied roles in plant vegetative and reproductive developmental processes as well as response to environmental stresses such has cold and high salinity. We focus here on a small conserved sub-group of GUCT domain-containing RNA helicase in the moss Physcomitrium patens. Phylogenetic analysis shows that RNA helicases containing the GUCT domain form a distinct conserved clade across the green lineage. In this clade, the P. patens genome possesses two closely related paralogues RNA helicases predicted to be nuclear, PpRH1 and PpRH2. Using in-locus gene fluorescent tagging we show that PpRH1 is localized to the nucleus in protonema. Analysis of PpRH1 and PpRH2 deletions, individually and together, indicates their potential roles in protonema, gametophore and sporophyte cellular and tissue development in P. patens. Additionally, the ultrastructural analysis of phyllid chloroplasts in Δrh2 and Δrh1/2 shows distinct starch granule accumulation under standard growth conditions associated with changes in photosynthetic activity parameters. We could not detect effects of either temperature or stress on protonema growth or PpRH1 and PpRH2 expression. Together, these results suggest that nuclear GUCT-containing RNA helicases play a role primarily in developmental processes directly or indirectly linked to photosynthesis activity in the moss P. patens.

SUPPLEMENTARY INFORMATION

The online version contains supplementary material available at 10.1007/s11103-021-01152-w.

Challenges of in vivo protein localization in plants seen through the DEK1 protein lens

During the last 25 y, fluorescent protein tagging has become a tool of choice to investigate protein function in a cellular context. The information gathered with this approach is not only providing insights into protein subcellular localization but also allows contextualizing protein function in multicellular settings. Here we illustrate the power of this method by commenting on the recent successful localization of the large membrane DEK1 protein during three-dimensional body formation in the moss Physcomitrella patens. But as many approaches, protein tagging is not exempt of caveats. The multiple infructuous (failed) attempts to detect DEK1 using a fluorescent protein tag present a good overview of such potential problems. Here we discuss the insertion of different fluorescent proteins at different positions in the PpDEK1 protein and the resulting unintended range of mutant phenotypes. Albeit none of these mutants generated a detectable fluorescent signal they can still provide interesting biological information about DEK1 function.

Single Nucleotide Polymorphism Charting of P. patens Reveals Accumulation of Somatic Mutations During in vitro Culture on the Scale of Natural Variation by Selfing

Introduction: Physcomitrium patens (Hedw.) Mitten (previously known as Physcomitrella patens) was collected by H.L.K. Whitehouse in Gransden Wood (Huntingdonshire, United Kingdom) in 1962 and distributed across the globe starting in 1974. Hence, the Gransden accession has been cultured in vitro in laboratories for half a century. Today, there are more than 13 different pedigrees derived from the original accession. Additionally, accessions from other sites worldwide were collected during the last decades. Methods and Results: In this study, 250 high throughput RNA sequencing (RNA-seq) samples and 25 gDNA samples were used to detect single nucleotide polymorphisms (SNPs). Analyses were performed using five different P. patens accessions and 13 different Gransden pedigrees. SNPs were overlaid with metadata and known phenotypic variations. Unique SNPs defining Gransden pedigrees and accessions were identified and experimentally confirmed. They can be successfully employed for PCR-based identification. Conclusion: We show independent mutations in different Gransden laboratory pedigrees, demonstrating that somatic mutations occur and accumulate during in vitro culture. The frequency of such mutations is similar to those observed in naturally occurring populations. We present evidence that vegetative propagation leads to accumulation of deleterious mutations, and that sexual reproduction purges those. Unique SNP sets for five different P. patens accessions were isolated and can be used to determine individual accessions as well as Gransden pedigrees. Based on that, laboratory methods to easily determine P. patens accessions and Gransden pedigrees are presented.

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.

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.

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.

Analysis of the Localization of Fluorescent PpROP1 and PpROP-GEF4 Fusion Proteins in Moss Protonemata Based on Genomic "Knock-In" and Estradiol-Titratable Expression

Tip growth of pollen tubes, root hairs, and apical cells of moss protonemata is controlled by ROP (Rho of plants) GTPases, which were shown to accumulate at the apical plasma membrane of these cells. However, most ROP localization patterns reported in the literature are based on fluorescent protein tagging and need to be interpreted with caution, as ROP fusion proteins were generally overexpressed at undefined levels, in many cases without assessing effects on tip growth. ROP-GEFs, important regulators of ROP activity, were also described to accumulate at the apical plasma membrane during tip growth. However, to date only the localization of fluorescent ROP-GEF fusion proteins strongly overexpressed using highly active promoters have been investigated. Here, the intracellular distributions of fluorescent PpROP1 and PpROP-GEF4 fusion proteins expressed at essentially endogenous levels in apical cells of Physcomitrella patens "knock-in" protonemata were analyzed. Whereas PpROP-GEF4 was found to associate with a small apical plasma membrane domain, PpROP1 expression was below the detection limit. Estradiol-titratable expression of a fluorescent PpROP1 fusion protein at the lowest detectable level, at which plant development was only marginally affected, was therefore employed to show that PpROP1 also accumulates at the apical plasma membrane, although within a substantially larger domain. Interestingly, RNA-Seq data indicated that the majority of all genes active in protonemata are expressed at lower levels than PpROP1, suggesting that estradiol-titratable expression may represent an important alternative to "knock-in" based analysis of the intracellular distribution of fluorescent fusion proteins in protonemal cells.

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.

An intragenic mutagenesis strategy in Physcomitrella patens to preserve intron splicing

Gene targeting is a powerful reverse genetics technique for site-specific genome modification. Intrinsic homologous recombination in the moss Physcomitrella patens permits highly effective gene targeting, a characteristic that makes this organism a valuable model for functional genetics. Functional characterization of domains located within a multi-domain protein depends on the ability to generate mutants harboring genetic modifications at internal gene positions while maintaining the reading-frames of the flanking exons. In this study, we designed and evaluated different gene targeting constructs for targeted gene manipulation of sequences corresponding to internal domains of the DEFECTIVE KERNEL1 protein in Physcomitrella patens. Our results show that gene targeting-associated mutagenesis of introns can have adverse effects on splicing, corrupting the normal reading frame of the transcript. We show that successful genetic modification of internal sequences of multi-exon genes depends on gene-targeting strategies which insert the selection marker cassette into the 5' end of the intron and preserve the nucleotide sequence of the targeted intron.

Combination of the Endogenous lhcsr1 Promoter and Codon Usage Optimization Boosts Protein Expression in the Moss Physcomitrella patens

The moss Physcomitrella patens is used both as an evo-devo model and biotechnological production system for metabolites and pharmaceuticals. Strong in vivo expression of genes of interest is important for production of recombinant proteins, e.g., selectable markers, fluorescent proteins, or enzymes. In this regard, the choice of the promoter sequence as well as codon usage optimization are two important inside factors to consider in order to obtain optimum protein accumulation level. To reliably quantify fluorescence, we transfected protoplasts with promoter:GFP fusion constructs and measured fluorescence intensity of living protoplasts in a plate reader system. We used the red fluorescent protein mCherry under 2x 35S promoter control as second reporter to normalize for different transfection efficiencies. We derived a novel endogenous promoter and compared deletion variants with exogenous promoters. We used different codon-adapted green fluorescent protein (GFP) genes to evaluate the influence of promoter choice and codon optimization on protein accumulation in P. patens, and show that the promoter of the gene of P. patens chlorophyll a/b binding protein lhcsr1 drives expression of GFP in protoplasts significantly (more than twofold) better than the commonly used 2x 35S promoter or the rice actin1 promoter. We identified a shortened 677 bp version of the lhcsr1 promoter that retains full activity in protoplasts. The codon optimized GFP yields significantly (more than twofold) stronger fluorescence signals and thus demonstrates that adjusting codon usage in P. patens can increase expression strength. In combination, new promotor and codon optimized GFP conveyed sixfold increased fluorescence signal.

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.

Ancient trans-Acting siRNAs Confer Robustness and Sensitivity onto the Auxin Response

Novel developmental programs often evolve via cooption of existing genetic networks. To understand this process, we explored cooption of the TAS3 tasiRNA pathway in the moss Physcomitrella patens. We find an ancestral function for this repeatedly redeployed pathway in the spatial regulation of a conserved set of Auxin Response Factors. In moss, this results in stochastic patterning of the filamentous protonemal tissue. Through modeling and experimentation, we demonstrate that tasiRNA regulation confers sensitivity and robustness onto the auxin response. Increased auxin sensitivity parallels increased developmental sensitivity to nitrogen, a key environmental signal. We propose that the properties lent to the auxin response network, along with the ability to stochastically modulate development in response to environmental cues, have contributed to repeated cooption of the tasiRNA-ARF module during evolution. The signaling properties of a genetic network, and not just its developmental output, are thus critical to understanding evolution of multicellular forms.

Endogenous Small-Noncoding RNAs and Potential Functions in Desiccation Tolerance in Physcomitrella Patens

Early land plants like moss Physcomitrella patens have developed remarkable drought tolerance. Phytohormone abscisic acid (ABA) protects seeds during water stress by activating genes through transcription factors such as ABSCISIC ACID INSENSITIVE (ABI3). Small noncoding RNA (sncRNA), including microRNAs (miRNAs) and endogenous small-interfering RNAs (endo-siRNAs), are key gene regulators in eukaryotes, playing critical roles in stress tolerance in plants. Combining next-generation sequencing and computational analysis, we profiled and characterized sncRNA species from two ABI3 deletion mutants and the wild type P. patens that were subject to ABA treatment in dehydration and rehydration stages. Small RNA profiling using deep sequencing helped identify 22 novel miRNAs and 6 genomic loci producing trans-acting siRNAs (ta-siRNAs) including TAS3a to TAS3e and TAS6. Data from degradome profiling showed that ABI3 genes (ABI3a/b/c) are potentially regulated by the plant-specific miR536 and that other ABA-relevant genes are regulated by miRNAs and ta-siRNAs. We also observed broad variations of miRNAs and ta-siRNAs expression across different stages, suggesting that they could potentially influence desiccation tolerance. This study provided evidence on the potential roles of sncRNA in mediating desiccation-responsive pathways in early land plants.

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

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.

DEK1; missing piece in puzzle of plant development

Patterning of land plant bodies is determined by positioning of cell walls. A crucial event in land plant evolution was the ability to utilize spatial information to direct cell wall deposition. Recent studies of DEK1 in Physcomitrella patens support a role for DEK1 in position dependent cell wall orientation.

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.

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.

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.

An experimental method to facilitate the identification of hybrid sporophytes in the moss Physcomitrella patens using fluorescent tagged lines

• The sequencing of the Physcomitrella patens genome, combined with the high frequency of gene targeting in this species, makes it ideal for reverse genetic studies. For forward genetic studies, experimental crosses and genetic analysis of progeny are essential. • Since P. patens is monoicous, producing both male and female gametes on the same gametophore, and undergoing self-fertilization at a high frequency, the identification of crossed sporophytes is difficult. Usually spores from many sporophytes from a mixed culture must be tested for the production of recombinant progeny. • Here, we describe the use of transgenic lines that express a fluorescent transgene constitutively, to provide a direct visual screen for hybrid sporophytes. • We show that segregations in crosses obtained with this technique are as expected, and demonstrate its utility for the study of the rate of outcrossing between three isolates of P. patens.

Protection of Telomeres 1 is required for telomere integrity in the moss Physcomitrella patens

In vertebrates, the single-stranded telomeric DNA binding protein Protection of Telomeres 1 (POT1) shields chromosome ends and prevents them from eliciting a DNA damage response. By contrast, Arabidopsis thaliana encodes two divergent full-length POT1 paralogs that do not exhibit telomeric DNA binding in vitro and have evolved to mediate telomerase regulation instead of chromosome end protection. To further investigate the role of POT1 in plants, we established the moss Physcomitrella patens as a new model for telomere biology and a counterpoint to Arabidopsis. The sequence and architecture of the telomere tract is similar in P. patens and Arabidopsis, but P. patens harbors only a single-copy POT1 gene. Unlike At POT1 proteins, Pp POT1 efficiently bound single-stranded telomeric DNA in vitro. Deletion of the P. patens POT1 gene resulted in the rapid onset of severe developmental defects and sterility. Although telomerase activity levels were unperturbed, telomeres were substantially shortened, harbored extended G-overhangs, and engaged in end-to-end fusions. We conclude that the telomere capping function of POT1 is conserved in early diverging land plants but is subsequently lost in Arabidopsis.

Transformation of the moss Physcomitrella patens using T-DNA mutagenesis

In this protocol, the transformation of moss (Physcomitrella patens) protoplasts is performed via Agrobacterium-mediated transfer of T-DNA. Protoplasts are incubated with Agrobacterium and acetoseringone in regeneration medium. They are then washed and plated on antibiotic-containing medium to select for T-DNA insertion in stable transformants. The transformation rate for this protocol is typically 10−4 (expressed as the frequency of stable transformants among regenerants surviving the transformation procedure).

Transformation of moss Physcomitrella patens gametophytes using a biolistic projectile delivery system

This method makes use of biolistic bombardment to propel gold particles coated with DNA into protonemal tissue of the moss Physcomitrella patens. It is especially suitable for transient gene expression studies, but it can be used to obtain stable transformants.

Chemical and UV mutagenesis of spores and protonemal tissue from the moss Physcomitrella patens

This protocol describes how to mutagenize spores and protonemal tissue from moss (Physcomitrella patens) using chemicals or ultraviolet (UV) light. Spores are mutagenized using the alkylating agents N-methyl-N′-nitro-N-nitrosoguanidine (NTG) and ethyl methanesulfonate (EMS), and protonemal tissue is mutagenized with NTG and UV light. Compared to alkylating agents, UV is less effective as a mutagen, but it may be advantageous because it is less hazardous and may not lead to clustered lesions.

Transformation of the moss Physcomitrella patens using direct DNA uptake by protoplasts

This protocol describes how to transform moss (Physcomitrella patens) protoplasts using polyethylene glycol (PEG)-mediated DNA uptake. The transformation rates for direct uptake by protoplasts of DNA with and without genomic sequence (a targeting construct) are typically 10−5 and 10−3, respectively. (These are the frequencies of stable transformants among regenerants surviving the transformation procedure.)

Somatic hybridization in the moss Physcomitrella patens using PEG-induced protoplast fusion

As an alternative to sexual crossing, protoplasts from two strains of moss (Physcomitrella patens) can be hybridized using polyethylene glycol (PEG). Although the efficiency is low, it requires no sophisticated apparatus. Hybrids are readily obtained using complementary auxotrophic mutants or strains with transgenic antibiotic resistance markers. It is now routine to obtain hybrids using transgenic strains that are hygromycin- or G418-resistant by selecting hybrids on medium containing both antibiotics.

BRICK1 is required for apical cell growth in filaments of the moss Physcomitrella patens but not for gametophore morphology

When BRK1, a member of the Wave/SCAR complex, is deleted in Physcomitrella patens (Deltabrk1), we report a striking reduction of filament growth resulting in smaller and fewer cells with misplaced cross walls compared with the normal protonemal cells. Using an inducible green fluorescent protein-talin to detect actin in living tissue, a characteristic broad accumulation of actin is observed at the tip of wild-type apical cells, whereas in Deltabrk1, smaller, more distinct foci of actin are present. Insertion of brk1-yfp into Deltabrk1 rescues the mutant phenotype and results in BRK1 being localized only in the tip of apical cells, the exclusive site of cell extension and division in the filament. Like BRK1, ARPC4 and At RABA4d are normally localized at the tip of apical cells and their localization is correlated with rapid tip growth in filaments. However, neither marker accumulates in apical cells of Deltabrk1 filaments. Although the Deltabrk1 phenotypes in protonema are severe, the leafy shoots or gametophores are normally shaped but stunted. These and other results suggest that BRK1 functions directly or indirectly in the selective accumulation/stabilization of actin and other proteins required for polar cell growth of filaments but not for the basic structure of the gametophore.

Physcomitrella patens: mosses enter the genomic age

The sequenced genome of the moss Physcomitrella patens provides a powerful tool for comparative analyses of land plant genomes. In parallel, several tools for studying gene function have been developed in P. patens, including RNA interference, inducible promoters and gene targeting, a unique attribute of this plant system. The results of these initiatives are now being realized. For example, transcriptomic analyses illustrate commonalities among plant lineages in gene content, structure, and regulation. Transgenic studies show that the regulatory factors ABSCISIC ACID INSENSITIVE3 (ABI3) and LEAFY (LFY) have molecular functions that are conserved between moss and angiosperms, in spite of the fact that they function in non-homologous tissues. Future work in P. patens will contribute to our understanding of the molecular basis of plant development and evolution.

Identification and functional characterization of the moss Physcomitrella patens delta5-desaturase gene involved in arachidonic and eicosapentaenoic acid biosynthesis

The moss Physcomitrella patens contains high levels of arachidonic acid and lesser amounts of eicosapentaenoic acid. Here we report the identification and characterization of a delta5-desaturase from P. patens that is associated with the synthesis of these fatty acids. A full-length cDNA for this desaturase was identified by data base searches based on homology to sequences of known delta5-desaturase cDNAs from fungal and algal species. The resulting P. patens cDNA encodes a 480-amino acid polypeptide that contains a predicted N-terminal cytochrome b5-like domain as well as three histidine-rich domains. Expression of the enzyme in Saccharomyces cerevisiae resulted in the production of the delta5-containing fatty acid arachidonic acid in cells that were provided di-homo-gamma-linolenic acid. In addition, the expressed enzyme generated delta5-desaturation products with the C20 substrates omega-6 eicosadienoic and omega-3 eicosatrienoic acids, but no products were detected with the C18 fatty acid linoleic and alpha-linolenic acids or with the C22 fatty acid adrenic and docosapentaenoic acids. When the corresponding P. patens genomic sequence was disrupted by replacement through homologous recombination, a dramatic alteration in the fatty acid composition was observed, i.e. an increase in di-homo-gamma-linolenic and eicosatetraenoic acids accompanied by a concomitant disappearance of the delta5-fatty acid arachidonic and eicosapentaenoic acids. In addition, overexpression of the P. patens cDNA in protoplasts isolated from a disrupted line resulted in the restoration of arachidonic acid synthesis.

The role of ARPC4 in tip growth and alignment of the polar axis in filaments of Physcomitrella patens

When the actin related protein 2/3 (Arp2/3) complex member arpc4 was deleted in Physcomitrella patens (moss), the resulting null mutant (Deltaarpc4) was viable and revealed no gross changes during morphogenesis of filaments into gametophores. However, we observed a striking reduction of filamentous tip growth, resulting in smaller, denser colonies. Although polar responses of Deltaarpc4 filaments to unilateral white light were unaffected, these mutant filaments were defective in their response to polarized white light. These observations strongly suggest a specific role of the Arp2/3 complex as a downstream target for signals regulating oriented tip growth. Insertion of YFP-ARPC4 into Deltaarpc4 rescued the mutant phenotypes and localized ARPC4 exclusively to the tip cell of filaments, the site of actin dynamics and polarized growth. The ability of Deltaarpc4 to perform some but not all cellular responses will allow the study of its function in orientation of tip growth in response to directional cues (e.g. light) in a viable but mutated background.

Molecular cloning, spatial and temporal characterization of spinach SOGA1 cDNA, encoding an alpha subunit of G protein

Heterotrimeric G proteins are an important component of signal transduction pathway in animals. Although these proteins have been described in plants, their exact function and action mode are not clearly defined. In order to analyze the relationship between these proteins and the transduction of light signals in spinach, we have isolated by 5' and 3' RACE-PCR a 1660bp cDNA clone called SOGA1. This codes for a 383aa protein, which reveals a very strong homology with other plant Galpha subunit sequences. Genomic analysis suggested that SOGA1 belonged to a small multiple gene family. Northern blots and in-situ hybridization analyses showed that SOGA1 transcripts accumulate in all organs tested with a specific high level associated with the apex, roots and hypocotyls. Finally, a time-course analysis performed on the green tissues showed that accumulation of SOGA1 transcripts follows a circadian rhythm. However, in-situ hybridization analysis of the apex suggested the opposite behavior, while no variation was observed in the hypocotyl.

Guanosine triphosphate-binding proteins on the plasmalemma of spinach leaf cells

The molecular mechanism of light perception through phytochrome is not well understood. This red-light photosensor has been implicated in various physiological processes, including the photoinduction of flowering. A few recent studies have shown that phytochrome initiates signal transduction chains via guanosine triphosphate (GTP)-binding proteins (G-proteins). We show here by different approaches that G-proteins exist in spinach (Spinacia oleracea L. cv. Nobel). Binding of GTP on the plasmalemma has been partially characterized and its possible regulation by red light examined by in-vitro assays. These experiments indicate a clear regulation of GTP binding by red light and also by Mastoparan. At least three G-proteins or protein subunits were found to be associated with the plasmalemma of leaf cells. The use of an antibody raised against an animal Gβ subunit confirmed the presence of heterotrimeric G-proteins. Separation of a crude membrane extract by free-flow electrophoresis also showed that some G-proteins could exist on the tonoplast.