Search

The moss Physcomitrium (Physcomitrella) patens is a bryophyte that provides genetic information about the adaptation to the life on land by early Embryophytes and is a reference organism for comparative evolutionary studies in plants. Copper is an essential micronutrient for every living organism, its transport across the plasma membrane is achieved by the copper transport protein family COPT/CTR. Two genes related to the COPT family were identified in Physcomitrella patens, PpaCOPT1 and PpaCOPT2. Homology modelling of both proteins showed the presence of three putative transmembrane domains (TMD) and the Mx3M motif, constituting a potential Cu + selectivity filter present in other members of this family. Functional characterization of PpaCOPT1 and PpaCOPT2 in the yeast mutant ctr1Δctr3Δ restored its growth on medium with non-fermentable carbon sources at micromolar Cu concentrations, providing support that these two moss proteins function as high affinity Cu + transporters. Localization of PpaCOPT1 and PpaCOPT2 in yeast cells was observed at the tonoplast and plasma membrane, respectively. The heterologous expression of PpaCOPT2 in tobacco epidermal cells co-localized with the plasma membrane marker. Finally, only PpaCOPT1 was expressed in seven-day old protonema and was influenced by extracellular copper levels. This evidence suggests different roles of PpaCOPT1 and PpaCOPT2 in copper homeostasis in Physcomitrella patens.

The homologs of VASCULAR RELATED NAC-DOMAIN in the peat moss Sphagnum palustre were identified and these transcriptional activity as the VNS family was conserved. In angiosperms, xylem vessel element differentiation is governed by the master regulators VASCULAR RELATED NAC-DOMAIN6 (VND6) and VND7, encoding plant-specific NAC transcription factors. Although vessel elements have not been found in bryophytes, differentiation of the water-conducting hydroid cells in the moss Physcomitrella patens is regulated by VND homologs termed VND-NST-SOMBRERO (VNS) genes. VNS genes are conserved in the land plant lineage, but their functions have not been elucidated outside of angiosperms and P. patens. The peat moss Sphagnum palustre, of class Sphagnopsida in the phylum Bryophyta, does not have hydroids and instead uses hyaline cells with thickened, helical-patterned cell walls and pores to store water in the leaves. Here, we performed whole-transcriptome analysis and de novo assembly using next generation sequencing in S. palustre, obtaining sequences for 68,305 genes. Among them, we identified seven VNS-like genes, SpVNS1-A, SpVNS1-B, SpVNS2-A, SpVNS2-B, SpVNS3-A, SpVNS3-B, and SpVNS4-A. Transient expression of these VNS-like genes, with the exception of SpVNS2-A, in Nicotiana benthamiana leaf cells resulted in ectopic thickening of secondary walls. This result suggests that the transcriptional activity observed in other VNS family members is functionally conserved in the VNS homologs of S. palustre.

Plants utilize a plethora of peptide signals to regulate their immune response. Peptide ligands and their cognate receptors involved in immune signaling share common motifs among many species of vascular plants. However, the origin and evolution of immune peptides is still poorly understood. Here, we searched for genes encoding small secreted peptides in the genomes of three bryophyte lineages-mosses, liverworts and hornworts-that occupy a critical position in the study of land plant evolution. We found that bryophytes shared common predicted small secreted peptides (SSPs) with vascular plants. The number of SSPs is higher in the genomes of mosses than in both the liverwort Marchantia polymorpha and the hornwort Anthoceros sp. The synthetic peptide elicitors-AtPEP and StPEP-specific for vascular plants, triggered ROS production in the protonema of the moss Physcomitrella patens, suggesting the possibility of recognizing peptide ligands from angiosperms by moss receptors. Mass spectrometry analysis of the moss Physcomitrella patens, both the wild type and the Δcerk mutant secretomes, revealed peptides that specifically responded to chitosan treatment, suggesting their role in immune signaling.

This work demonstrates that PpCIPK1, a putative protein kinase, participates in regulating plant salt tolerance in moss Physcomitrella patens. Calcineurin B-Like protein (CBL)-interacting protein kinases (CIPKs) have been reported to be involved in multiple signaling networks and function in plant growth and stress responses, however, their biological functions in non-seed plants have not been well characterized. In this study, we report that PpCIPK1, a putative protein kinase, participates in regulating plant salt tolerance in moss Physcomitrella patens (P. patens). Phylogenetic analysis revealed that PpCIPK1 shared high similarity with its homologs in higher plants. PpCIPK1 transcription level was induced upon salt stress in P. patens. Using homologous recombination, we constructed PpCIPK1 knockout mutant lines (PpCIPK1 KO). Salt sensitivity analysis showed that independent PpCIPK1 KO plants exhibited severe growth inhibition and developmental deficiency of gametophytes under salt stress condition compared to that of wild-type P. patens (WT). Consistently, ionic homeostasis was disrupted in plants due to PpCIPK1 deletion, and high level of H2O2 was accumulated in PpCIPK1 KO than that in WT. Furthermore, PpCIPK1 functions in regulating photosynthetic activity in response to salt stress. Interestingly, we observed that PpCIPK1 could completely rescue the salt-sensitive phenotype of sos2-1 to WT level in Arabidopsis, indicating that AtSOS2 and PpCIPK1 are functionally conserved. In conclusion, our work provides evidence that PpCIPK1 participates in salt tolerance regulation in P. patens.

Evolutionary conserved defense mechanisms present in extant bryophytes and angiosperms, as well as moss-specific defenses are part of the immune response of Physcomitrium patens. Bryophytes and tracheophytes are descendants of early land plants that evolved adaptation mechanisms to cope with different kinds of terrestrial stresses, including drought, variations in temperature and UV radiation, as well as defense mechanisms against microorganisms present in the air and soil. Although great advances have been made on pathogen perception and subsequent defense activation in angiosperms, limited information is available in bryophytes. In this study, a transcriptomic approach uncovered the molecular mechanisms underlying the defense response of the bryophyte Physcomitrium patens (previously Physcomitrella patens) against the important plant pathogen Botrytis cinerea. A total of 3.072 differentially expressed genes were significantly affected during B. cinerea infection, including genes encoding proteins with known functions in angiosperm immunity and involved in pathogen perception, signaling, transcription, hormonal signaling, metabolic pathways such as shikimate and phenylpropanoid, and proteins with diverse role in defense against biotic stress. Similarly as in other plants, B. cinerea infection leads to downregulation of genes involved in photosynthesis and cell cycle progression. These results highlight the existence of evolutionary conserved defense responses to pathogens throughout the green plant lineage, suggesting that they were probably present in the common ancestors of land plants. Moreover, several genes acquired by horizontal transfer from prokaryotes and fungi, and a high number of P. patens-specific orphan genes were differentially expressed during B. cinerea infection, suggesting that they are important players in the moss immune response.

Plants harvest light energy utilized for photosynthesis by light-harvesting complex I and II (LHCI and LHCII) surrounding photosystem I and II (PSI and PSII), respectively. During the evolution of green plants, moss is at an evolutionarily intermediate position from aquatic photosynthetic organisms to land plants, being the first photosynthetic organisms that landed. Here, we report the structure of the PSI–LHCI supercomplex from the moss Physcomitrella patens (Pp) at 3.23 Å resolution solved by cryo-electron microscopy. Our structure revealed that four Lhca subunits are associated with the PSI core in an order of Lhca1–Lhca5–Lhca2–Lhca3. This number is much decreased from 8 to 10, the number of subunits in most green algal PSI–LHCI, but the same as those of land plants. Although Pp PSI–LHCI has a similar structure as PSI–LHCI of land plants, it has Lhca5, instead of Lhca4, in the second position of Lhca, and several differences were found in the arrangement of chlorophylls among green algal, moss, and land plant PSI–LHCI. One chlorophyll, PsaF–Chl 305, which is found in the moss PSI–LHCI, is located at the gap region between the two middle Lhca subunits and the PSI core, and therefore may make the excitation energy transfer from LHCI to the core more efficient than that of land plants. On the other hand, energy-transfer paths at the two side Lhca subunits are relatively conserved. These results provide a structural basis for unravelling the mechanisms of light-energy harvesting and transfer in the moss PSI–LHCI, as well as important clues on the changes of PSI–LHCI after landing.

Plants perceive dynamic light conditions and optimize their growth and development accordingly by regulating gene expression at multiple levels. Alternative splicing (AS), a widespread mechanism in eukaryotes that post-transcriptionally generates two or more messenger RNAs (mRNAs) from the same pre-mRNA, is rapidly controlled by light. However, a detailed mechanism of light-regulated AS is still not clear. In this study, we demonstrate that histone 3 lysine 36 trimethylation (H3K36me3) rapidly and differentially responds to light at specific gene loci with light-regulated intron retention (IR) of their transcripts in the moss Physcomitrella patens. However, the level of H3K36me3 following exposure to light is inversely related to that of IR events. Physcomitrella patens MORF-related gene 1 (PpMRG1), a chromatin adaptor, bound with higher affinity to H3K36me3 in light conditions than in darkness and was differentially targeted to gene loci showing light-responsive IR. Transcriptome analysis indicated that PpMRG1 functions in the regulation of light-mediated AS. Furthermore, PpMRG1 was also involved in red light-mediated phototropic responses. Our results suggest that light regulates histone methylation, which leads to alterations of AS patterns. The chromatin adaptor PpMRG1 potentially participates in light-mediated AS, revealing that chromatin-coupled regulation of pre-mRNA splicing is an important aspect of the plant’s response to environmental changes.

Key message 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. Abstract 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.

Cytosine (DNA) methylation in plants regulates the expression of genes and transposons. While methylation in plant genomes occurs at CG, CHG, and CHH sequence contexts, the comparative roles of the individual methylation contexts remain elusive. Here, we present Physcomitrella patens as the second plant system, besides Arabidopsis thaliana, with viable mutants with an essentially complete loss of methylation in the CG and non-CG contexts. In contrast to A. thaliana, P. patens has more robust CHH methylation, similar CG and CHG methylation levels, and minimal cross-talk between CG and non-CG methylation, making it possible to study context-specific effects independently. Our data found CHH methylation to act in redundancy with symmetric methylation in silencing transposons and to regulate the expression of CG/CHG-depleted transposons. Specific elimination of CG methylation did not dysregulate transposons or genes. In contrast, exclusive removal of non-CG methylation massively up-regulated transposons and genes. In addition, comparing two exclusively but equally CG- or CHG-methylated genomes, we show that CHG methylation acts as a greater transcriptional regulator than CG methylation. These results disentangle the transcriptional roles of CG and non-CG, as well as symmetric and asymmetric methylation in a plant genome, and point to the crucial role of non-CG methylation in genome regulation.

With global climate change, water scarcity threatens whole agro/ecosystems. The desert moss Syntrichia caninervis, an extremophile, offers novel insights into surviving desiccation and heat. The sequenced S. caninervis genome consists of 13 chromosomes containing 16 545 protein-coding genes and 2666 unplaced scaffolds. Syntenic relationships within the S. caninervis and Physcomitrella patens genomes indicate the S. caninervis genome has undergone a single whole genome duplication event (compared to two for P. patens) and evidence suggests chromosomal or segmental losses in the evolutionary history of S. caninervis. The genome contains a large sex chromosome composed primarily of repetitive sequences with a large number of Copia and Gypsy elements. Orthogroup analyses revealed an expansion of ELIP genes encoding proteins important in photoprotection. The transcriptomic response to desiccation identified four structural clusters of novel genes. The genomic resources established for this extremophile offer new perspectives for understanding the evolution of desiccation tolerance in plants.

Homologous genes for the peptidoglycan precursor flippase MurJ, and peptidoglycan hydrolases: lytic transglycosylase MltB, and DD-carboxypeptidase VanY are required for chloroplast division in the moss Physcomitrella patens. The moss Physcomitrella patens is used as a model plant to study plastid peptidoglycan biosynthesis. In bacteria, MurJ flippase transports peptidoglycan precursors from the cytoplasm to the periplasm. In this study, we identified a MurJ homolog (PpMurJ) in the P. patens genome. Bacteria employ peptidoglycan degradation and recycling pathways for cell division. We also searched the P. patens genome for genes homologous to bacterial peptidoglycan hydrolases and identified genes homologous for the lytic transglycosylase mltB, N-acetylglucosaminidase nagZ, and LD-carboxypeptidase ldcA in addition to a putative DD-carboxypeptidase vanY reported previously. Moreover, we found a ß-lactamase-like gene (Pplactamase). GFP fusion proteins with either PpMltB or PpVanY were detected in the chloroplasts, whereas fusion proteins with PpNagZ, PpLdcA, or Pplactamase localized in the cytoplasm. Experiments seeking PpMurJ-GFP fusion proteins failed. PpMurJ gene disruption in P. patens resulted in the appearance of macrochloroplasts in protonemal cells. Compared with the numbers of chloroplasts in wild-type plants (38.9 ± 4.9), PpMltB knockout and PpVanY knockout had lower numbers of chloroplasts (14.3 ± 6.7 and 28.1 ± 5.9, respectively). No differences in chloroplast numbers were observed after PpNagZ, PpLdcA, or Pplactamase single-knockout. Chloroplast numbers in PpMltB/PpVanY double-knockout cells were similar to those in PpMltB single-knockout cells. Zymogram analysis of the recombinant PpMltB protein revealed its peptidoglycan hydrolase activity. Our results imply that PpMurJ, PpMltB and PpVanY play a critical role in chloroplast division in the moss P. patens.

Cytidine-to-uridine RNA editing is a posttranscriptional process in plant organelles, mediated by specific pentatricopeptide repeat (PPR) proteins. In angiosperms, hundreds of sites undergo RNA editing. By contrast, only 13 sites are edited in the moss Physcomitrium (Physcomitrella) patens. Some are conserved between the two species, like the mitochondrial editing site nad5eU598RC. The PPR proteins assigned to this editing site are known in both species: the DYW-type PPR protein PPR79 in P. patens and the E+-type PPR protein CWM1 in Arabidopsis (Arabidopsis thaliana). CWM1 also edits sites ccmCeU463RC, ccmBeU428SL, and nad5eU609VV. Here, we reciprocally expressed the P. patens and Arabidopsis editing factors in the respective other genetic environment. Surprisingly, the P. patens editing factor edited all target sites when expressed in the Arabidopsis cwm1 mutant background, even when carboxy-terminally truncated. Conversely, neither Arabidopsis CWM1 nor CWM1-PPR79 chimeras restored editing in P. patens ppr79 knockout plants. A CWM1-like PPR protein from the early diverging angiosperm macadamia (Macadamia integrifolia) features a complete DYW domain and fully rescued editing of nad5eU598RC when expressed in P. patens. We conclude that (1) the independently evolved P. patens editing factor PPR79 faithfully operates in the more complex Arabidopsis editing system, (2) truncated PPR79 recruits catalytic DYW domains in trans when expressed in Arabidopsis, and (3) the macadamia CWM1-like protein retains the capacity to work in the less complex P. patens editing environment.

Phytochromes are red and far-red light receptors in plants that control growth and development in response to changes in the environment. Light-activated phytochromes enter the nucleus and act on a set of downstream signalling components to regulate gene expression. PHYTOCHROME INTERACTING FACTORs (PIFs) belong to the basic helix-loop-helix family of transcription factors and directly bind to light-activated phytochromes. Potential homologues of PIFs have been identified in ferns, bryophytes and streptophyte algae, and it has been shown that the potential PIF homologues from Physcomitrella patens, PIF1 to PIF4, have PIF function when expressed in Arabidopsis. However, their function in Physcomitrella is still unknown. Seed plant PIFs bind to G-box-containing promoters and, therefore, we searched the Physcomitrella genome for genes that contain G-boxes in their promoter. Here, we show that Physcomitrella PIFs activate these promoters in a G-box-dependent manner, suggesting that they could be direct PIF targets. Furthermore, we generated Physcomitrella pif1, pif2, pif3 and pif4 knock out mutant lines and quantified the expression of potential PIF direct target genes. The expression of these genes was generally reduced in pif mutants compared to the wildtype, but for several genes, the relative induction upon a short light treatment was higher in pif mutants than the wildtype. In contrast, expression of these genes was strongly repressed in continuous light, and pif mutants showed partial downregulation of these genes in the dark. Thus, the overall function of PIFs in light-regulated gene expression might be an ancient property of PIFs.

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.

LBDs, as the plant-specific gene family, play essential roles in lateral organ development, plant regeneration, as well as abiotic stress and pathogen response. However, the number and characteristic of LBD genes in Pyscomitrella patens were still obscure. This study was performed to identify the LBD family gene in moss and to determine the expression profiles of LBDs under the abiotic and pathogen stress. Complete genome sequences and transcriptomes of P. patens were downloaded from the Ensembl plant database. The hidden Markov model-based profile of the conserved LOB domain was submitted as a query to identify all potential LOB domain sequences with HMMER software. Expression profiles of PpLBDs were obtained based on the GEO public database and qRT-PCR analysis. In this study, a total of 31 LBDs were identified in the P. patens genome, divided into two classes based on the presence of the leucine zipper-like coiled-coil motif. A phylogenetic relationship was obtained between 31 proteins from P. patens and 43 proteins from the Arabidopsis thaliana genome, providing insights into their conserved and potential functions. Furthermore, the exon–intron organization of each PpLBD were analyzed. All PpLBD contain the conserved DNA binding motif (CX2CX6CX3C zinc finger-like motif), and were predicted to be located in cell nuclear. The 31 PpLBD genes were unevenly assigned to 18 out of 27 chromosomes based on the physical positions. Among these genes, PpLBD27 was not only remarkably highest expressed in desiccation, but also a susceptible gene to pathogens through jasmonic acid-mediated signaling pathway. Most of PpLBDs were up-regulated with the treatment of mannitol. These results showed they were differentially induced and their potential functions in the environmental stimulus of the early terrestrial colonizers. Despite significant differences in the life cycle in P. patens and flowering plants, their functions involved in abiotic and biotic stress-regulated by LBDs have been identified and appear to be conserved in the two lineages. These results provided a comprehensive analysis of PpLBDs and paved insights into studies aimed at a better understanding of PpLBDs.

Plant colonization of land has been intimately associated with mycorrhizae or mycorrhizae-like fungi. Despite the pivotal role of fungi in plant adaptation, it remains unclear whether and how gene acquisition following fungal interaction might have affected the development of land plants. Here we report a macro2 domain gene in bryophytes that is likely derived from Mucoromycota, a group that includes some mycorrhizae-like fungi found in the earliest land plants. Experimental and transcriptomic evidence suggests that this macro2 domain gene in the moss Physcomitrella patens, PpMACRO2, is important in epigenetic modification, stem cell function, cell reprogramming and other processes. Gene knockout and over-expression of PpMACRO2 significantly change the number and size of gametophores. These findings provide insights into the role of fungal association and the ancestral gene repertoire in the early evolution of land plants.

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

SQUAMOSA PROMOTER BINDING PROTEIN-LIKE (SPL) transcription factors play critical roles in regulating diverse aspects of plant growth and development, including vegetative phase change, plant architecture, anthocyanin accumulation, lateral root growth, etc. In the present study, 15 SPL genes were identified based on the genome data of Codonopsis pilosula, a well-known medicinal plant. Phylogenetic analysis clustered CpSPLs into eight groups (G1-G8) along with SPLs from Arabidopsis thaliana, Solanum lycopersicum, Oryza sativa and Physcomitrella patens. CpSPLs in the same group share similar gene structure and conserved motif composition. Cis-acting elements responding to light, stress and phytohormone widely exist in their promoter regions. Our qRT-PCR results indicated that 15 CpSPLs were differentially expressed in different tissues (root, stem, leaf, flower and calyx), different developmental periods (1, 2 and 3 months after germination) and various conditions (NaCl, MeJA and ABA treatment). Compared with the control, overexpression of CpSPL2 or CpSPL10 significantly promoted not only the growth of hairy roots, but also the accumulation of total saponins and lobetyolin. Our results established a foundation for further investigation of CpSPLs and provided novel insights into their biological functions. As far as we know, this is the first experimental research on gene function in C. pilosula.

eIF4A is a RNA-stimulated ATPase and helicase. Besides its key role in regulating cap-dependent translation initiation in eukaryotes, it also performs specific functions in regulating cell cycle progression, plant growth and abiotic stress tolerance. Flowering plants encode three eIF4A paralogues, eIF4A1, eIF4A2 and eIF4A3 that share conserved sequence motifs but differ in functions. To date, however, no information is available on eIF4A in basal land plants. In this study we report that genome of the moss Physcomitrella patens encodes multiple eIF4A genes. The encoded proteins possess the highly conserved motifs characteristic of the DEAD box helicases. Spatial expression analysis shows these genes to be ubiquitously expressed in all tissue types with Pp3c6_1080V3.1 showing high expression in filamentous protonemata. Targeted deletion of conserved core motifs in Pp3c6_1080V3.1 slowed protonemata growth and resulted in dwarfing of leafy gametophores suggesting a role for Pp3c6_1080V3.1 in regulating cell division/elongation. Rapid and strong induction of Pp3c6_1080V3.1 under salt stress and slow recovery of knockout plants upon exposure to high salt further suggest Pp3c6_1080V3.1 to be involved in stress management in P. patens. Protein-protein interaction studies that show Pp3c6_1080V3.1 to interact with the Physcomitrella heterogenous ribonucleoprotein, LIF2L1, a transcriptional regulator of stress-responsive genes in Arabidopsis. The results presented in this study provide insight into evolutionary conserved functions of eIF4A and shed light on the novel link between eIF4A activities and stress mitigation pathways/RNA metabolic processes in P. patens.

Plants perceive environmental light conditions and optimize their growth and development accordingly by regulating gene activity at multiple levels. Photoreceptors are important for light sensing and downstream gene regulation. Phytochromes, red/far-red light receptors, are believed to regulate light-responsive alternative splicing, but little is known about the underlying mechanism. Alternative splicing is primarily regulated by transacting factors, such as splicing regulators, and by cis-acting elements in precursor mRNA. In the moss Physcomitrella patens, we show that phytochrome 4 (PpPHY4) directly interacts with a splicing regulator, heterogeneous nuclear ribonucleoprotein F1 (PphnRNP-F1), in the nucleus to regulate light-responsive alternative splicing. RNA sequencing analysis revealed that PpPHY4 and PphnRNP-F1 coregulate 70% of intron retention (IR) events in response to red light. A repetitive GAA motif was identified to be an exonic splicing silencer that controls red light-responsive IR. Biochemical studies indicated that PphnRNP-F1 is recruited by the GAA motif to form RNA-protein complexes. Finally, red light elevates PphnRNP-F1 protein levels via PpPHY4, increasing levels of IR. We propose that PpPHY4 and PphnRNP-F1 regulate alternative splicing through an exonic splicing silencer to control splicing machinery activity in response to light.

Due to their high efficiency, specificity, and flexibility, programmable nucleases, such as those of the clustered regularly interspaced short palindromic repeats (CRISPR)/Cas12a (Cpf1) system, have greatly expanded the applicability of editing the genomes of various organisms. Genes from different gene families or genes with redundant functions in the same gene family can be examined by assembling multiple CRISPR RNAs (crRNAs) in a single vector. However, the activity and efficiency of CRISPR/Cas12a in the non-vascular plant Physcomitrella patens are largely unknown. Here, we demonstrate that LbCas12a together with its mature crRNA can target multiple loci simultaneously in P. patens with high efficiency via co-delivery of LbCas12a and a crRNA expression cassette in vivo. The mutation frequencies induced by CRISPR/LbCas12a at a single locus ranged from 26.5 to 100%, with diverse deletions being the most common type of mutation. Our method expands the repertoire of genome editing tools available for P. patens and facilitates the creation of loss-of-function mutants of multiple genes from different gene families.

Plant U-box (PUB) E3 ubiquitin ligases play crucial roles in the plant response to abiotic stress and the phytohormone abscisic acid (ABA) signaling, but little is known about them in bryophytes. Here, a representative U-box armadillo repeat (PUB-ARM) ubiquitin E3 ligase from Antarctic moss Pohlia nutans (PnSAG1), was explored for its role in abiotic stress response in Arabidopsis thaliana and Physcomitrella patens. The expression of PnSAG1 was rapidly induced by exogenous abscisic acid (ABA), salt, cold and drought stresses. PnSAG1 was localized to the cytoplasm and showed E3 ubiquitin ligase activity by in vitro ubiquitination assay. The PnSAG1-overexpressing Arabidopsis enhanced the sensitivity with respect to ABA and salt stress during seed germination and early root growth. Similarly, heterogeneous overexpression of PnSAG1 in P. patens was more sensitive to the salinity and ABA in their gametophyte growth. The analysis by RT-qPCR revealed that the expression of salt stress/ABA-related genes were downregulated in PnSAG1-overexpressing plants after salt treatment. Taken together, our results indicated that PnSAG1 plays a negative role in plant response to ABA and salt stress.

Inactive auxin conjugates are accumulated in plants and hydrolyzed to recover phytohormone action. A family of metallopeptidase orthologues has been conserved in Plantae to help regulate auxin homeostatic levels during growth and development. This hydrolase family was recently traced back to liverwort, the most ancient extant land plant lineage. Liverwort’s auxin hydrolase has little activity against auxin conjugate substrates and does not appear to actively regulate auxin. This finding, along with data that shows moss can synthesize auxin conjugates, led to examining another bryophyte lineage, Physcomitrella patens. We have identified and isolated three M20D hydrolase paralogues from moss. The isolated enzymes strongly recognize and cleave a variety of auxin conjugates, including those of indole butyric and indole propionic acids. These P. patens hydrolases not only appear to be “cryptic”, but they are likely to have derived from soil bacteria through Horizontal Gene Transfer. Additionally, support is presented that the plant-type M20D peptidase family may have been universally lost from mosses after divergence from the common ancestor with liverwort.

Chloroplast gene expression is controlled by numerous nuclear-encoded RNA-binding proteins. Among these, pentatricopeptide repeat (PPR) proteins are known to be key players in post-transcriptional regulation in chloroplasts. However, the functions of many PPR proteins remain unknown. In this study, we characterized the function of a chloroplast-localized P-class PPR protein PpPPR_21 in Physcomitrella patens. Knockout (KO) mutants of PpPPR_21 exhibited reduced protonemata growth and lower photosynthetic activity. Immunoblot analysis and blue-native gel analysis showed a remarkable reduction of the photosystem II (PSII) reaction center protein and poor formation of the PSII supercomplexes in the KO mutants. To assess whether PpPPR_21 is involved in chloroplast gene expression, chloroplast genome-wide microarray analysis and Northern blot hybridization were performed. These analyses indicated that the psbI-ycf12 transcript encoding the low molecular weight subunits of PSII did not accumulate in the KO mutants while other psb transcripts accumulated at similar levels in wild-type and KO mutants. A complemented PpPPR_21KO moss transformed with the cognate full-length PpPPR_21cDNA rescued the level of accumulation of psbI-ycf12 transcript. RNA-binding experiments showed that the recombinant PpPPR_21 bound efficiently to the 5' untranslated and translated regions of psbImRNA. The present study suggests that PpPPR_21 may be essential for the accumulation of a stable psbI-ycf12mRNA.

In vascular (Arabidopsis thaliana) and non-vascular (Physcomitrella patens) plants, PHOSPHATE 1 (PHO1) homologs play important roles in the acquisition and transfer of phosphate. The tomato genome contains six genes (SlPHO1;1-SlPHO1;6) homologous to AtPHO1. The six proteins have typical characteristics of the plant PHO1 family, such as the three Syg1/Pho81/XPRI (SPX) subdomains in the N-terminal portion and one ERD1/XPR1/SYG1 (EXS) domain in the C-terminal portion. Phylogenetic analysis revealed that the SlPHO1 family is subdivided into three clusters. A pairwise comparison indicated that SlPHO1;1 showed the highest level of sequence identity/similarity (67.39/76.21%) to AtPHO1. SlPHO1;1 deletion mutants induced by clustered regularly interspaced short palindromic repeats (CRISPR)/Cas9 displayed typical phenotypes of Pi starvation, such as decreased shoot fresh weight and increased root fresh weight, therefore having a greater root-to-shoot ratio. Mutants also accumulated more anthocyanin and had more soluble Pi content in the root and less in the shoot. These results indicate that SlPHO1;1 plays an important role in Pi transport in the tomato at seedling stage.

Background: MicroRNAs (miRNAs) are groups of small non-protein-coding endogenous single stranded RNAs with approximately 18-24 nucleotides in length. High evolutionary sequence conservation of miRNAs among plant species and availability of powerful computational tools allow identification of new orthologs and paralogs. Methods: New conserved miRNAs in P. patens were found by EST-based homology search approaches. All candidates were screened according to a series of miRNA filtering criteria. Unigene, DFCI Gene Index (PpspGI) databases and psRNATarget algorithm were applied to identify target transcripts using P. patens putative conserved miRNA sequences. Results: Nineteen conserved P. patens miRNAs were identified. The sequences were homologous to known reference plant mature miRNA from 10 miRNA families. They could be folded into the typical miRNA secondary structures. RepeatMasker algorithm demonstrated that ppt-miR2919e and pptmiR1533 had simple sequence repeats in their sequences. Target sites (49 genes) were identified for 7 out of 19 miRNAs. GO and KEGG analysis of targets indicated the involvement of some in important multiple biological and metabolic processes. Conclusion: The majority of the registered miRNAs in databases were predicted by computational approaches while many more have remained unknown. Due to the conserved nature of miRNAs in plant species from closely to distantly related, homology search-based approaches between plants species could lead to the identification of novel miRNAs in other plant species providing baseline information for further search about the biological functions and evolution of miRNAs.

Premise of the study We used the model species Physcomitrella patens to examine chirality in moss gametophores. Chirality is manifested in the direction of consecutive apical cell divisions, cell plate configurations, and deviations of leaf connecting lines from the vertical course. However, the frequencies of chiral configurations of all these processes as well as their mutual dependence-especially in the case of gametophore branching-are not known. Other moss species were checked to determine the universality of our findings. Methods The gametophore structure of Physcomitrella patens grown in the laboratory under controlled conditions was investigated using light microscopy and compared with that of other moss species collected from their natural stands. Key results In all investigated moss species, the tetrahedral apical cell exhibits either clockwise or counterclockwise consecutive divisions, and selection of this directionality in the primary axis is random. It is, however, related to cell plate configuration. If the plate is skewed, leaf-producing segments arising from the apical cell cleavage exhibit circumferential rotation. Three parallel lines connecting the leaves deviate from a vertical course, but always in the same direction as that of leaf initiation; thus, the angular distance between consecutive leaves increases to >120°. Lateral branches are exclusively antidromous. Conclusions Gametophore chiral configuration appears to be useful in resolving problems of moss modular growth and branching. Morphological and anatomical evidence strongly suggests that an unknown direction-sensing mechanism controls the development of moss axial organs. We propose that leaves are responsible for a horizontal gradient of sugar signals that develops along the gametophore circumference, thus influencing branching-unit chirality.

Plant cell walls are highly dynamic structures that are composed predominately of polysaccharides. As such, endogenous carbohydrate active enzymes (CAZymes) are central to the synthesis and subsequent modification of plant cells during morphogenesis. The endo-glucanase 16 (EG16) members constitute a distinct group of plant CAZymes, angiosperm orthologs of which were recently shown to have dual β-glucan/xyloglucan hydrolase activity. Molecular phylogeny indicates that EG16 members comprise a sister clade with a deep evolutionary relationship to the widely studied apoplastic xyloglucan endo-transglycosylases/hydrolases (XTH). A cross-genome survey indicated that EG16 members occur as a single ortholog across species and are widespread in early diverging plants, including the non-vascular bryophytes, for which functional data were previously lacking. Remarkably, enzymological characterization of an EG16 ortholog from the model moss Physcomitrella patens (PpEG16) revealed that EG16 activity and sequence/structure are highly conserved across 500 million years of plant evolution, vis-à-vis orthologs from grapevine and poplar. Ex vivo biomechanical assays demonstrated that the application of EG16 gene products caused abrupt breakage of etiolated hypocotyls rather than slow extension, thereby indicating a mode-of-action distinct from endogenous expansins and microbial endo-glucanases. The biochemical data presented here will inform future genomic, genetic, and physiological studies of EG16 enzymes.

Destabilization of organelle genomes causes organelle dysfunction that appears as abnormal growth in plants and diseases in human. In plants, loss of the bacterial-type homologous recombination repair (HRR) factors RECA and RECG induces organelle genome instability. In this study, we show the landscape of organelle genome instability in Physcomitrella patens HRR knockout mutants by deep sequencing in combination with informatics approaches. Genome-wide maps of rearrangement positions in the organelle genomes, which exhibited prominent mutant-specific patterns, were highly biased in terms of direction and location and often associated with dramatic variation in read depth. The rearrangements were location-dependent and mostly derived from the asymmetric products of microhomology-mediated recombination. Our results provide an overall picture of organelle-specific gross genomic rearrangements in the HRR mutants, and suggest that chloroplasts and mitochondria share common mechanisms for replication-related rearrangements., Masaki Odahara and Kensuke Nakamura et al. use deep paired-end sequencing to examine organellar genome recombination when homologous recombination repair genes are individually knocked out in the moss, Physcomitrella patens. Their results suggest that chloroplasts and mitochondria share a common mechanism for replication-related rearrangements.

Genes encoding ribosomal RNA (rDNA) are essential for cell survival and are particularly sensitive to factors leading to genomic instability. Their repetitive character makes them prone to inappropriate recombinational events arising from collision of transcriptional and replication machineries, resulting in unstable rDNA copy numbers. In this review, we summarize current knowledge on the structure and organization of rDNA, its role in sensing changes in the genome, and its linkage to aging. We also review recent findings on the main factors involved in chromatin assembly and DNA repair in the maintenance of rDNA stability in the model plants Arabidopsis thaliana and the moss Physcomitrella patens, providing a view across the plant evolutionary tree.

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., This work was supported by grants from the Spanish Ministerio de Economía, Industria y Competitividad (AGL2016-78992-R/FEDER) and Ministerio de Ciencia e Innovación (PID2019-106374RB-I00/AEI/10.13039/501100011033) to JMC. The IJPB benefits from the support of the LabEx Saclay Plant Sciences (SPS) (ANR-10-LABX-0040-SPS).

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.

ANGUSTIFOLIA (AN) is a plant-specific subfamily of the CtBP/BARS/AN family, characterized by a plant-specific C-terminal domain of approximately 200 amino acids. Previously, we revealed that double knockout (DKO) lines of Physcomitrium (Physcomitrella) patens ANGUSTIFOLIA genes (PpAN1-1 and PpAN1-2) show defects in gametophore height and the lengths of the seta and foot region of sporophytes, by reduced cell elongation. In addition to two canonical ANs, the genome of P. patens has two atypical ANs without a coding region for a plant-specific C-terminus (PpAN2-1 and PpAN2-2); these were investigated in this study. Similar to PpAN1s, both promoters of the PpAN2 genes were highly active in the stems of haploid gametophores and in the middle-to-basal region of young diploid sporophytes that develop into the seta and foot. Analyses of PpAN2-1/2-2 DKO and PpAN quadruple knockout (QKO) lines implied that these four AN genes have partially redundant functions to regulate cell elongation in their expression regions. Transgenic strains harboring P. patens α-tubulin fused to green fluorescent protein, which were generated from a QKO line, showed that the orientation of the microtubules in the gametophore tips in the PpAN QKO lines was unchanged from the wild-type and PpAN1-1/1-2 DKO plants. In addition to both PpAN2-1 and PpAN2-2, short Arabidopsis AN without the C-terminus of 200 amino acids could rescue the Arabidopsis thaliana an-1 phenotypes, implying AN activity is dependent on the N-terminal regions.

β-carbonic anhydrases, which function in regulating plant growth, C/N status, and stomata number, showed functional redundancy and divergence in Physcomitrella patens. Carbonic anhydrases (CAs) catalyze the interconversion of CO2 and HCO3−. Plants have three evolutionarily unrelated CA families: α-, β-, and γ-CAs. βCAs are abundant in plants and are involved in CO2 assimilation, stress responses, and stomata formation. Recent studies of βCAs have mainly examined C3 or C4 plants, whereas their functions in non-vascular plants are mostly unknown. In this study, phylogenetic analysis revealed that the evolution of βCAs were conserved between subaerial green algae and bryophytes after terrestrialization event, and βCAs from some cyanobacteria might begin evolving for the adaptation of terrestrial environment/habitat. In addition, we investigated the physiological roles of βCAs in the basal land plant Physcomitrella patens. High PpβCA expression levels in different tissues suggest that PpβCAs play important roles in development in P. patens. Plants treated with 1–10 mM NaHCO3 had higher fresh and dry weight, PpβCA expression, total CA activity, and photosynthetic yield (Fv/Fm) compared with water-treated plants. However, treatment with 10 mM NaHCO3 influenced the C/N status. Further study of six Ppβca single-gene mutants revealed that PpβCAs have functional redundancy and divergence in regulating the C/N ratio of plants and stomatal formation. This study provides new insight into the physiological roles of βCAs in basal land plants.

RNA interference (RNAi) enables flexible and dynamic interrogation of entire gene families or essential genes without the need for exogenous proteins, unlike CRISPR-Cas technology. Unfortunately, isolation of plants undergoing potent gene silencing requires laborious design, visual screening, and physical separation for downstream characterization. Here, we developed an adenine phosphoribosyltransferase (APT)-based RNAi technology (APTi) in Physcomitrella patens that improves upon the multiple limitations of current RNAi techniques. APTi exploits the prosurvival output of transiently silencing APT in the presence of 2-fluoroadenine, thereby establishing survival itself as a reporter of RNAi. To maximize the silencing efficacy of gene targets, we created vectors that facilitate insertion of any gene target sequence in tandem with the APT silencing motif. We tested the efficacy of APTi with two gene families, the actin-dependent motor, myosin XI (a,b), and the putative chitin receptor Lyk5 (a,b,c). The APTi approach resulted in a homogenous population of transient P. patens mutants specific for our gene targets with zero surviving background plants within 8 d. The observed mutants directly corresponded to a maximal 93% reduction of myosin XI protein and complete loss of chitin-induced calcium spiking in the Lyk5-RNAi background. The positive selection nature of APTi represents a fundamental improvement in RNAi technology and will contribute to the growing demand for technologies amenable to high-throughput phenotyping.

When DNA double-strand breaks occur, four-stranded DNA structures called Holliday junctions (HJs) form during homologous recombination. Because HJs connect homologous DNA by a covalent link, resolution of HJ is crucial to terminate homologous recombination and segregate the pair of DNA molecules faithfully. We recently identified Monokaryotic Chloroplast1 (MOC1) as a plastid DNA HJ resolvase in algae and plants. Although Cruciform cutting endonuclease1 (CCE1) was identified as a mitochondrial DNA HJ resolvase in yeasts, homologs or other mitochondrial HJ resolvases have not been identified in other eukaryotes. Here, we demonstrate that MOC1 depletion in the green alga Chlamydomonas reinhardtii and the moss Physcomitrella patens induced ectopic recombination between short dispersed repeats in ptDNA. In addition, MOC1 depletion disorganized thylakoid membranes in plastids. In some land plant lineages, such as the moss P. patens, a liverwort and a fern, MOC1 dually targeted to plastids and mitochondria. Moreover, mitochondrial targeting of MOC1 was also predicted in charophyte algae and some land plant species. Besides causing instability of plastid DNA, MOC1 depletion in P. patens induced short dispersed repeat-mediated ectopic recombination in mitochondrial DNA and disorganized cristae in mitochondria. Similar phenotypes in plastids and mitochondria were previously observed in mutants of plastid-targeted (RECA2) and mitochondrion-targeted (RECA1) recombinases, respectively. These results suggest that MOC1 functions in the double-strand break repair in which a recombinase generates HJs and MOC1 resolves HJs in mitochondria of some lineages of algae and plants as well as in plastids in algae and plants.

As a broad-spectrum disease resistance factor, MLO is involved in a variety of biotic and abiotic stress responses in plants. To figure out the structural features, phylogenetic relationships, and expression patterns of MLO genes, we investigated the genome and transcriptome sequencing data of 28 plant species using bioinformatics tools. A total of 197 MLO genes were identified. They possessed 5&ndash, 7 transmembrane domains, but only partially contained a calmodulin-binding domain. A total of 359 polymorphic sites and 142 haplotypes were found in 143 sequences, indicating the rich nucleotide diversity of MLO genes. The MLO genes were unevenly distributed on chromosomes or scaffolds and were mainly located at the ends, forming clusters (24.1% genes), tandem duplicates (5.7%), and segment duplicates (36.2%). The MLO genes could be classified into three groups by phylogenetic analysis. The angiosperm genes were mainly in subgroup IA, Selaginella moellendorffii genes were in subgroup IA and IIIB, Physcomitrella patens genes were in subgroup IB and IIIA, and almost all algae genes were in group II. About half of the MLO genes had homologs within and across species. The Ka/Ks values were all less than 1, varying 0.01&ndash, 0.78, suggesting that purifying selection had occurred in MLO gene evolution. In tomato, RNA-seq data indicated that SlMLO genes were highly expressed in roots, followed by flowers, buds, and leaves, and also regulated by different biotic stresses. qRT&ndash, PCR analysis revealed that SlMLO genes could respond to tomato bacterial wilt, with SlMLO1, SlMLO2, SlMLO4, and SlMLO6 probably involved in the susceptibility response, whereas SlMLO14 and SlMLO16 being the opposite. These results lay a foundation for the isolation and application of related genes in plant disease resistance breeding.

DNA methylation in plants occurs in CG, CHG, and CHH sites. While depletion of CG methylation in transposons is associated with ample transcriptional activation, it was mainly studied in species with limited non-CG methylation that is linked to CG methylation. Here we profiled transcription in the moss plant, Physcomitrella patens, that has robust non-CG methylation with similar symmetrical CG and CHG methylation levels. Separated contextual methylation mechanisms in Physcomitrella patens enabled generation of numerous context-specific hypomethylated mutants. Our transcriptome data show that specific elimination of CG methylation is fully complemented by non-CG methylation. Conversely, exclusive removal of non-CG methylation massively dysregulated genes and transposons. Moreover, CHG methylation silenced transposons stronger than CG methylation. Lastly, we found non-CG methylation as crucial for silencing CG-depleted transposons. These results demonstrate the potency of non-CG methylation in genome regulation and suggest that it evolved due to moderate silencing and/or rapid mutability of methylated CGs.