CLAVATA modulates auxin homeostasis and transport to regulate stem cell identity and plant shape in a moss.

The CLAVATA pathway is a key regulator of stem cell function in the multicellular shoot tips of Arabidopsis, where it acts via the WUSCHEL transcription factor to modulate hormone homeostasis. Broad-scale evolutionary comparisons have shown that CLAVATA is a conserved regulator of land plant stem cell function, but CLAVATA acts independently of WUSCHEL-like (WOX) proteins in bryophytes. The relationship between CLAVATA, hormone homeostasis and the evolution of land plant stem cell functions is unknown. Here we show that in the moss, Physcomitrella (Physcomitrium patens), CLAVATA affects stem cell activity by modulating hormone homeostasis. CLAVATA pathway genes are expressed in the tip cells of filamentous tissues, regulating cell identity, filament branching, plant spread and auxin synthesis. The receptor-like kinase PpRPK2 plays the major role, and Pprpk2 mutants have abnormal responses to cytokinin, auxin and auxin transport inhibition, and show reduced expression of PIN auxin transporters. We propose a model whereby PpRPK2 modulates auxin gradients in filaments to determine stem cell identity and overall plant form. Our data indicate that CLAVATA-mediated auxin homeostasis is a fundamental property of plant stem cell function, probably exhibited by the last shared common ancestor of land plants.

Roles of RAD51 and RTEL1 in telomere and rDNA stability in Physcomitrella patens.

Telomeres and ribosomal RNA genes (rDNA) are essential for cell survival and particularly sensitive to factors affecting genome stability. Here, we examine the role of RAD51 and its antagonist, RTEL1, in the moss Physcomitrella patens. In corresponding mutants, we analyse their sensitivity to DNA damage, the maintenance of telomeres and rDNA, and repair of double-stranded breaks (DSBs) induced by genotoxins with various modes of action. While the loss of RTEL1 results in rapid telomere shortening, concurrent loss of both RAD51 genes has no effect on telomere lengths. We further demonstrate here the linked arrangement of 5S and 45S rRNA genes in P. patens. The spacer between 5S and 18S rRNA genes, especially the region downstream from the transcription start site, shows conspicuous clustering of sites with a high propensity to form quadruplex (G4) structures. Copy numbers of 5S and 18S rDNA are reduced moderately in the pprtel1 mutant, and significantly in the double pprad51-1-2 mutant, with no progression during subsequent cultivation. While reductions in 45S rDNA copy numbers observed in pprtel1 and pprad51-1-2 plants apply also to 5S rDNA, changes in transcript levels are different for 45S and 5S rRNA, indicating their independent transcription by RNA polymerase I and III, respectively. The loss of SOL (Sog One-Like), a transcription factor regulating numerous genes involved in DSB repair, increases the rate of DSB repair in dividing as well as differentiated tissue, and through deactivation of G2/M cell-cycle checkpoint allows the cell-cycle progression manifested as a phenotype resistant to bleomycin.

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 ancestral stomatal patterning module revealed in the non-vascular land plant Physcomitrella patens.

The patterning of stomata plays a vital role in plant development and has emerged as a paradigm for the role of peptide signals in the spatial control of cellular differentiation. Research in Arabidopsis has identified a series of epidermal patterning factors (EPFs), which interact with an array of membrane-localised receptors and associated proteins (encoded by ERECTA and TMM genes) to control stomatal density and distribution. However, although it is well-established that stomata arose very early in the evolution of land plants, until now it has been unclear whether the established angiosperm stomatal patterning system represented by the EPF/TMM/ERECTA module reflects a conserved, universal mechanism in the plant kingdom. Here, we use molecular genetics to show that the moss Physcomitrella patens has conserved homologues of angiosperm EPF, TMM and at least one ERECTA gene that function together to permit the correct patterning of stomata and that, moreover, elements of the module retain function when transferred to Arabidopsis Our data characterise the stomatal patterning system in an evolutionarily distinct branch of plants and support the hypothesis that the EPF/TMM/ERECTA module represents an ancient patterning system.

Genetic Analysis of Physcomitrella patens Identifies ABSCISIC ACID NON-RESPONSIVE, a Regulator of ABA Responses Unique to Basal Land Plants and Required for Desiccation Tolerance.

The anatomically simple plants that first colonized land must have acquired molecular and biochemical adaptations to drought stress. Abscisic acid (ABA) coordinates responses leading to desiccation tolerance in all land plants. We identified ABA nonresponsive mutants in the model bryophyte Physcomitrella patens and genotyped a segregating population to map and identify the ABA NON-RESPONSIVE (ANR) gene encoding a modular protein kinase comprising an N-terminal PAS domain, a central EDR domain, and a C-terminal MAPKKK-like domain. anr mutants fail to accumulate dehydration tolerance-associated gene products in response to drought, ABA, or osmotic stress and do not acquire ABA-dependent desiccation tolerance. The crystal structure of the PAS domain, determined to 1.7-Å resolution, shows a conserved PAS-fold that dimerizes through a weak dimerization interface. Targeted mutagenesis of a conserved tryptophan residue within the PAS domain generates plants with ABA nonresponsive growth and strongly attenuated ABA-responsive gene expression, whereas deleting this domain retains a fully ABA-responsive phenotype. ANR orthologs are found in early-diverging land plant lineages and aquatic algae but are absent from more recently diverged vascular plants. We propose that ANR genes represent an ancestral adaptation that enabled drought stress survival of the first terrestrial colonizers but were lost during land plant evolution.

Giant peroxisomes in a moss (Physcomitrella patens) peroxisomal biogenesis factor 11 mutant.

Peroxisomal biogenesis factor 11 (PEX11) proteins are found in yeasts, mammals and plants, and play a role in peroxisome morphology and regulation of peroxisome division. The moss Physcomitrella patens has six PEX11 isoforms which fall into two subfamilies, similar to those found in monocots and dicots. We carried out targeted gene disruption of the Phypa_PEX11-1 gene and compared the morphological and cellular phenotypes of the wild-type and mutant strains. The mutant grew more slowly and the development of gametophores was retarded. Mutant chloronemal filaments contained large cellular structures which excluded all other cellular organelles. Expression of fluorescent reporter proteins revealed that the mutant strain had greatly enlarged peroxisomes up to 10 µm in diameter. Expression of a vacuolar membrane marker confirmed that the enlarged structures were not vacuoles, or peroxisomes sequestered within vacuoles as a result of pexophagy. Phypa_PEX11 targeted to peroxisome membranes could rescue the knock out phenotype and interacted with Fission1 on the peroxisome membrane. Moss PEX11 functions in peroxisome division similar to PEX11 in other organisms but the mutant phenotype is more extreme and environmentally determined, making P. patens a powerful system in which to address mechanisms of peroxisome proliferation and division.

RAD51B plays an essential role during somatic and meiotic recombination in Physcomitrella.

The eukaryotic RecA homologue Rad51 is a key factor in homologous recombination and recombinational repair. Rad51-like proteins have been identified in yeast (Rad55, Rad57 and Dmc1), plants and vertebrates (RAD51B, RAD51C, RAD51D, XRCC2, XRCC3 and DMC1). RAD51 and DMC1 are the strand-exchange proteins forming a nucleofilament for strand invasion, however, the function of the paralogues in the process of homologous recombination is less clear. In yeast the two Rad51 paralogues, Rad55 and Rad57, have been shown to be involved in somatic and meiotic HR and they are essential to the formation of the Rad51/DNA nucleofilament counterbalancing the anti-recombinase activity of the SRS2 helicase. Here, we examined the role of RAD51B in the model bryophyte Physcomitrella patens. Mutant analysis shows that RAD51B is essential for the maintenance of genome integrity, for resistance to DNA damaging agents and for gene targeting. Furthermore, we set up methods to investigate meiosis in Physcomitrella and we demonstrate that the RAD51B protein is essential for meiotic homologous recombination. Finally, we show that all these functions are independent of the SRS2 anti-recombinase protein, which is in striking contrast to what is found in budding yeast where the RAD51 paralogues are fully dependent on the SRS2 anti-recombinase function.

Conservation of Male Sterility 2 function during spore and pollen wall development supports an evolutionarily early recruitment of a core component in the sporopollenin biosynthetic pathway.

The early evolution of plants required the acquisition of a number of key adaptations to overcome physiological difficulties associated with survival on land. One of these was a tough sporopollenin wall that enclosed reproductive propagules and provided protection from desiccation and UV-B radiation. All land plants possess such walled spores (or their derived homologue, pollen). We took a reverse genetics approach, consisting of knock-out and complementation experiments to test the functional conservation of the sporopollenin-associated gene MALE STERILTY 2 (which is essential for pollen wall development in Arabidopsis thaliana) in the bryophyte Physcomitrella patens. Knock-outs of a putative moss homologue of the A. thaliana MS2 gene, which is highly expressed in the moss sporophyte, led to spores with highly defective walls comparable to that observed in the A. thaliana ms2 mutant, and extremely compromised germination. Conversely, the moss MS2 gene could not rescue the A. thaliana ms2 phenotype. The results presented here suggest that a core component of the biochemical and developmental pathway required for angiosperm pollen wall development was recruited early in land plant evolution but the continued increase in pollen wall complexity observed in angiosperms has been accompanied by divergence in MS2 gene function.

MRE11 and RAD50, but not NBS1, are essential for gene targeting in the moss Physcomitrella patens.

The moss Physcomitrella patens is unique among plant models for the high frequency with which targeted transgene insertion occurs via homologous recombination. Transgene integration is believed to utilize existing machinery for the detection and repair of DNA double-strand breaks (DSBs). We undertook targeted knockout of the Physcomitrella genes encoding components of the principal sensor of DNA DSBs, the MRN complex. Loss of function of PpMRE11 or PpRAD50 strongly and specifically inhibited gene targeting, whilst rates of untargeted transgene integration were relatively unaffected. In contrast, disruption of the PpNBS1 gene retained the wild-type capacity to integrate transforming DNA efficiently at homologous loci. Analysis of the kinetics of DNA-DSB repair in wild-type and mutant plants by single-nucleus agarose gel electrophoresis revealed that bleomycin-induced fragmentation of genomic DNA was repaired at approximately equal rates in each genotype, although both the Ppmre11 and Pprad50 mutants exhibited severely restricted growth and development and enhanced sensitivity to UV-B and bleomycin-induced DNA damage, compared with wild-type and Ppnbs1 plants. This implies that while extensive DNA repair can occur in the absence of a functional MRN complex; this is unsupervised in nature and results in the accumulation of deleterious mutations incompatible with normal growth and development.

The K/HDEL receptor does not recycle but instead acts as a Golgi-gatekeeper.

Accurately measuring the ability of the K/HDEL receptor (ERD2) to retain the ER cargo Amy-HDEL has questioned earlier results on which the popular receptor recycling model is based upon. Here we demonstrate that ERD2 Golgi-retention, rather than fast ER export supports its function. Ligand-induced ERD2 redistribution is only observed when the C-terminus is masked or mutated, compromising the signal that prevents Golgi-to-ER transport of the receptor. Forcing COPI mediated retrograde transport destroys receptor function, but introducing ER-to-Golgi export or cis-Golgi retention signals re-activate ERD2 when its endogenous Golgi-retention signal is masked or deleted. We propose that ERD2 remains fixed as a Golgi gatekeeper, capturing K/HDEL proteins when they arrive and releasing them again into a subdomain for retrograde transport back to the ER. An in vivo ligand:receptor ratio far greater than 100 to 1 strongly supports this model, and the underlying mechanism appears to be extremely conserved across kingdoms.

A minus-end directed kinesin motor directs gravitropism in Physcomitrella patens.

Gravity is a critical environmental factor regulating directional growth and morphogenesis in plants, and gravitropism is the process by which plants perceive and respond to the gravity vector. The cytoskeleton is proposed to play important roles in gravitropism, but the underlying mechanisms are obscure. Here we use genetic screening in Physcomitrella patens, to identify a locus GTRC, that when mutated, reverses the direction of protonemal gravitropism. GTRC encodes a processive minus-end-directed KCHb kinesin, and its N-terminal, C-terminal and motor domains are all essential for transducing the gravity signal. Chimeric analysis between GTRC/KCHb and KCHa reveal a unique role for the N-terminus of GTRC in gravitropism. Further study shows that gravity-triggered normal asymmetric distribution of actin filaments in the tip of protonema is dependent on GTRC. Thus, our work identifies a microtubule-based cellular motor that determines the direction of plant gravitropism via mediating the asymmetric distribution of actin filaments.

The SOS1 transporter of Physcomitrella patens mediates sodium efflux in planta.

• SOS1 is an Na(+)/H(+) antiporter that plays a central role in Na(+) tolerance in land plants. SOS1 mediation of Na(+) efflux has been studied in plasma-membrane vesicles and deduced from the SOS1 suppression of the Na(+) sensitivity of yeast mutants defective in Na(+) -efflux. However, SOS1-mediated Na(+) efflux has not been characterized in either plant or yeast cells. Here, we use Physcomitrella patens to investigate the function of SOS1 in planta. • In P. patens, a nonvascular plant in which the study of ion cellular fluxes is technically simple, the existence of two SOS1 genes suggests that the Na(+) efflux remaining after the deletion of the ENA1 ATPase is mediated by a SOS1 system. Therefore, we cloned the P. patens SOS1 and SOS1B genes (PpSOS1 and PpSOS1B, respectively) and complementary DNAs, and constructed the PpΔsos1 and PpΔena1/PpΔsos1 deletion lines by gene targeting. • Comparison of wild-type, and PpΔsos1 and PpΔena1/PpΔsos1 mutant lines revealed that PpSOS1 is crucial for Na(+) efflux and that the PpΔsos1 line, and especially the PpΔena1/PpΔsos1 lines, showed excessive Na(+) accumulation and Na(+)-triggered cell death. The PpΔsos1 and PpΔena1/PpΔsos1 lines showed impaired high-affinity K(+) uptake. • Our data support the hypothesis that PpSOS1 mediates cellular Na(+) efflux and that PpSOS1 enhances K(+) uptake by an indirect effect.

A sequence-anchored genetic linkage map for the moss, Physcomitrella patens.

The moss Physcomitrella patens is a model for the study of plant cell biology and, by virtue of its basal position in land plant phylogeny, for comparative analysis of the evolution of plant gene function and development. It is ideally suited for 'reverse genetic' analysis by virtue of its outstanding ability to undertake targeted transgene integration by homologous recombination. However, gene identification through mutagenesis and map-based cloning has hitherto not been possible, due to the lack of a genetic linkage map. Using molecular markers [amplified fragment length polymorphisms (AFLP) and simple sequence repeats (SSR)] we have generated genetic linkage maps for Physcomitrella. One hundred and seventy-nine gene-specific SSR markers were mapped in 46 linkage groups, and 1574 polymorphic AFLP markers were identified. Integrating the SSR- and AFLP-based maps generated 31 linkage groups comprising 1420 markers. Anchorage of the integrated linkage map with gene-specific SSR markers coupled with computational prediction of AFLP loci has enabled its correspondence with the newly sequenced Physcomitrella genome. The generation of a linkage map densely populated with molecular markers and anchored to the genome sequence now provides a resource for forward genetic interrogation of the organism and for the development of a pipeline for the map-based cloning of Physcomitrella genes. This will radically enhance the potential of Physcomitrella for determining how gene function has evolved for the acquisition of complex developmental strategies within the plant kingdom.

Erratum: Publisher Correction: SnRK2 protein kinases represent an ancient system in plants for adaptation to a terrestrial environment.

[This corrects the article DOI: 10.1038/s42003-019-0281-1.].

SnRK2 protein kinases represent an ancient system in plants for adaptation to a terrestrial environment.

The SNF1-related protein kinase 2 (SnRK2) family includes key regulators of osmostress and abscisic acid (ABA) responses in angiosperms and can be classified into three subclasses. Subclass III SnRK2s act in the ABA response while ABA-nonresponsive subclass I SnRK2s are regulated through osmostress. Here we report that an ancient subclass III SnRK2-based signalling module including ABA and an upstream Raf-like kinase (ARK) exclusively protects the moss Physcomitrella patens from drought. Subclass III SnRK2s from both Arabidopsis and from the semiterrestrial alga Klebsormidium nitens, which contains all the components of ABA signalling except ABA receptors, complement Physcomitrella snrk2 - mutants, whereas Arabidopsis subclass I SnRK2 cannot. We propose that the earliest land plants developed the ABA/ARK/subclass III SnRK2 signalling module by recruiting ABA to regulate a pre-existing dehydration response and that subsequently a novel subclass I SnRK2 system evolved in vascular plants conferring osmostress protection independently from the ancient system.

The mechanism of gene targeting in Physcomitrella patens: homologous recombination, concatenation and multiple integration.

The model bryophyte Physcomitrella patens exhibits high frequencies of gene targeting when transformed with DNA constructs containing sequences homologous with genomic loci. 'Targeted gene replacement' (TGR) resulting from homologous recombination (HR) between each end of a targeting construct and the targeted locus occurs when either single or multiple targeting vectors are delivered. In the latter instance simultaneous, multiple, independent integration of different transgenes occurs at the targeted loci. In both single gene and 'batch' transformations, DNA can also be found to undergo 'targeted insertion' (TI), integrating at one end of the targeted locus by HR with one flanking sequence of the vector accompanied by an apparent non-homologous end-joining (NHEJ) event at the other. Untargeted integration at nonhomologous sites also occurs, but at a lower frequency. Molecular analysis of TI at a single locus shows that this occurs as a consequence of concatenation of the transforming DNA, in planta, prior to integration, followed by HR between a single site in the genomic target and two of its repeated homologues in the concatenated vector. This reinforces the view that HR is the major pathway by which transforming DNA is integrated in Physcomitrella.

CLAVATA Was a Genetic Novelty for the Morphological Innovation of 3D Growth in Land Plants.

How genes shape diverse plant and animal body forms is a key question in biology. Unlike animal cells, plant cells are confined by rigid cell walls, and cell division plane orientation and growth rather than cell movement determine overall body form. The emergence of plants on land coincided with a new capacity to rotate stem cell divisions through multiple planes, and this enabled three-dimensional (3D) forms to arise from ancestral forms constrained to 2D growth. The genes involved in this evolutionary innovation are largely unknown. The evolution of 3D growth is recapitulated during the development of modern mosses when leafy shoots arise from a filamentous (2D) precursor tissue. Here, we show that a conserved, CLAVATA peptide and receptor-like kinase pathway originated with land plants and orients stem cell division planes during the transition from 2D to 3D growth in a moss, Physcomitrella. We find that this newly identified role for CLAVATA in regulating cell division plane orientation is shared between Physcomitrella and Arabidopsis. We report that roles for CLAVATA in regulating cell proliferation and cell fate are also shared and that CLAVATA-like peptides act via conserved receptor components in Physcomitrella. Our results suggest that CLAVATA was a genetic novelty enabling the morphological innovation of 3D growth in land plants.

Parameters determining the efficiency of gene targeting in the moss Physcomitrella patens.

In the moss Physcomitrella patens, transforming DNA containing homologous sequences integrates predominantly by homologous recombination with its genomic target. A systematic investigation of the parameters that determine gene targeting efficiency shows a direct relationship between homology length and targeting frequency for replacement vectors (a selectable marker flanked by homologous DNA). Overall homology of only 1 kb is sufficient to achieve a 50% yield of targeted transformants. Targeting may occur through homologous recombination in one arm, accompanied by non-homologous end-joining by the other arm of the vector, or by allele replacement following two homologous recombination events. Allele replacement frequency depends on the symmetry of the targeting vector, being proportional to the length of the shorter arm. Allele replacement may involve insertion of multiple copies of the transforming DNA, accompanied by ectopic insertions at non-homologous sites. Single-copy and single insertions at targeted loci (targeted gene replacements, 'TGR') occur with a frequency of 7-20% of all transformants when the minimum requirements for allele replacement are met. Homologous recombination in Physcomitrella is substantially more efficient than in any multicellular eukaryote, recommending it as the outstanding model for the study of homologous recombination in plants.

Towards designer organelles by subverting the peroxisomal import pathway.

The development of 'designer' organelles could be a key strategy to enable foreign pathways to be efficiently controlled within eukaryotic biotechnology. A fundamental component of any such system will be the implementation of a bespoke protein import pathway that can selectively deliver constituent proteins to the new compartment in the presence of existing endogenous trafficking systems. Here we show that the protein-protein interactions that control the peroxisomal protein import pathway can be manipulated to create a pair of interacting partners that still support protein import in moss cells, but are orthogonal to the naturally occurring pathways. In addition to providing a valuable experimental tool to give new insights into peroxisomal protein import, the variant receptor-signal sequence pair forms the basis of a system in which normal peroxisomal function is downregulated and replaced with an alternative pathway, an essential first step in the creation of a designer organelle.Designer organelles could allow the isolation of synthetic biological pathways from endogenous components of the host cell. Here the authors engineer a peroxisomal protein import pathway orthogonal to the naturally occurring system.

The Transcriptional Response to DNA-Double-Strand Breaks in Physcomitrella patens.

The model bryophyte Physcomitrella patens is unique among plants in supporting the generation of mutant alleles by facile homologous recombination-mediated gene targeting (GT). Reasoning that targeted transgene integration occurs through the capture of transforming DNA by the homology-dependent pathway for DNA double-strand break (DNA-DSB) repair, we analysed the genome-wide transcriptomic response to bleomycin-induced DNA damage and generated mutants in candidate DNA repair genes. Massively parallel (Illumina) cDNA sequencing identified potential participants in gene targeting. Transcripts encoding DNA repair proteins active in multiple repair pathways were significantly up-regulated. These included Rad51, CtIP, DNA ligase 1, Replication protein A and ATR in homology-dependent repair, Xrcc4, DNA ligase 4, Ku70 and Ku80 in non-homologous end-joining and Rad1, Tebichi/polymerase theta, PARP in microhomology-mediated end-joining. Differentially regulated cell-cycle components included up-regulated Rad9 and Hus1 DNA-damage-related checkpoint proteins and down-regulated D-type cyclins and B-type CDKs, commensurate with the imposition of a checkpoint at G2 of the cell cycle characteristic of homology-dependent DNA-DSB repair. Candidate genes, including ATP-dependent chromatin remodelling helicases associated with repair and recombination, were knocked out and analysed for growth defects, hypersensitivity to DNA damage and reduced GT efficiency. Targeted knockout of PpCtIP, a cell-cycle activated mediator of homology-dependent DSB resection, resulted in bleomycin-hypersensitivity and greatly reduced GT efficiency.