Tagged mutagenesis and gene-trap in the moss, Physcomitrella patens by shuttle mutagenesis

Δ4-dn-iso-OPDA, a bioactive plant hormone of Marchantia polymorpha

Significant progress has been recently made in our understanding of the evolution of jasmonates biosynthesis and signaling. The bioactive jasmonate activating COI1-JAZ co-receptor differs in bryophytes and vascular plants. Dinor-iso-12-oxo-phytodienoic acid (dn-iso-OPDA) is the bioactive hormone in bryophytes and lycophytes. However, further studies showed that the full activation of hormone signaling in Marchantia polymorpha requires additional unidentified hormones. Δ4-dn-OPDAs were previously identified as novel bioactive jasmonates in M. polymorpha. In this paper, we describe the major bioactive isomer of Δ4-dn-OPDAs as Δ4-dn-iso-OPDA through chemical synthesis, receptor binding assay, and biological activity in M. polymorpha. In addition, we disclosed that Δ4-dn-cis-OPDA is a biosynthetic precursor of Δ4-dn-iso-OPDA. We demonstrated that in planta cis-to-iso conversion of Δ4-dn-cis-OPDA occurs in the biosynthesis of Δ4-dn-iso-OPDA, defining a key biosynthetic step in the chemical evolution of hormone structure. We predict that these findings will facilitate further understanding of the molecular evolution of plant hormone signaling.

Near-infrared imaging of phytochrome-derived autofluorescence in plant nuclei

Capturing images of the nuclear dynamics within live cells is an essential technique for comprehending the intricate biological processes inherent to plant cell nuclei. While various methods exist for imaging nuclei, including combining fluorescent proteins and dyes with microscopy, there is a dearth of commercially available dyes for live-cell imaging. In Arabidopsis thaliana, we discovered that nuclei emit autofluorescence in the near-infrared (NIR) range of the spectrum and devised a non-invasive technique for the visualization of live cell nuclei using this inherent NIR autofluorescence. Our studies demonstrated the capability of the NIR imaging technique to visualize the dynamic behavior of nuclei within primary roots, root hairs, and pollen tubes, which are tissues that harbor a limited number of other organelles displaying autofluorescence. We further demonstrated the applicability of NIR autofluorescence imaging in various other tissues by incorporating fluorescence lifetime imaging techniques. Nuclear autofluorescence was also detected across a wide range of plant species, enabling analyses without the need for transformation. The nuclear autofluorescence in the NIR wavelength range was not observed in animal or yeast cells. Genetic analysis revealed that this autofluorescence was caused by the phytochrome protein. Our studies demonstrated that nuclear autofluorescence imaging can be effectively employed not only in model plants but also for studying nuclei in non-model plant species.

Ethylene controls three-dimensional growth involving reduced auxin levels in the moss Physcomitrium patens

The conquest of land by plants was concomitant with, and possibly enabled by, the evolution of three-dimensional (3D) growth. The moss Physcomitrium patens provides a model system for elucidating molecular mechanisms in the initiation of 3D growth. Here, we investigate whether the phytohormone ethylene, which is believed to have been a signal before land plant emergence, plays a role in 3D growth regulation in P. patens. We report ethylene controls 3D gametophore formation, based on results from exogenously applied ethylene and genetic manipulation of PpEIN2, which is a central component in the ethylene signaling pathway. Overexpression (OE) of PpEIN2 activates ethylene responses and leads to earlier formation of gametophores with fewer gametophores produced thereafter, phenocopying ethylene-treated wild-type. Conversely, Ppein2 knockout mutants, which are ethylene insensitive, show initially delayed gametophore formation with more gametophores produced later. Furthermore, pharmacological and biochemical analyses reveal auxin levels are decreased in the OE lines but increased in the knockout mutants. Our results suggest that evolutionarily, ethylene and auxin molecular networks were recruited to build the plant body plan in ancestral land plants. This might have played a role in enabling ancient plants to acclimate to the continental surfaces of the planet.

Autophagy Promotes Cell Death Induced by Hydrogen Peroxide in Physcomitrium patens

The autophagy-defective mutants (atg5 and atg7) of Physcomitrium patens exhibit strong desiccation tolerance. Here, we examined the effects of H2O2 on wild-type (WT) and autophagy-defective mutants of P. patens, considering that desiccation induces reactive oxygen species (ROS). We found that atg mutants can survive a 30-min treatment with 100 mM H2O2, whereas WT cannot, implying that autophagy promotes cell death induced by H2O2. Concomitant with cell death, vacuole collapse occurred. Intracellular H2O2 levels in both WT and atg5 increased immediately after H2O2 treatment and subsequently reached plateaus, which were higher in WT than in atg5. The ROS scavenger N-acetylcysteine lowered the plateau levels in WT and blocked cell death, suggesting that higher H2O2 plateau caused cell death. The uncoupler of electron transport chain (ETC) carbonyl cyanide m-chlorophenylhydrazone also lowered the H2O2 plateaus, showing that ROS produced in the ETC in mitochondria and/or chloroplasts elevated the H2O2 plateau. The autophagy inhibitor 3-methyladenine lowered the H2O2 plateau and the cell death rate in WT, suggesting that autophagy occurring after H2O2 treatment is involved in the production of ROS. Conversely, the addition of bovine serum albumin, which is endocytosed and supplies amino acids instead of autophagy, elevated the H2O2 plateau in atg5 cells, suggesting that amino acids produced through autophagy promote H2O2 generation. These results clearly show that autophagy causes cell death under certain stress conditions. We propose that autophagy-derived amino acids are catabolized using ETCs in mitochondria and/or chloroplasts and produce H2O2, which in turn promotes the cell death accompanying vacuole collapse.

PpSCARECROW1 (PpSCR1) regulates leaf blade and mid-vein development in Physcomitrium patens

In plants, asymmetric cell divisions result in distinct cell fates forming large and small daughter cells, adding to the cellular diversity in an organ. SCARECROW (SCR), a GRAS domain-containing transcription factor controls asymmetric periclinal cell divisions in flowering plants by governing radial patterning of ground tissue in roots and cell proliferation in leaves. Though SCR homologs are present across land plant lineages, the current understanding of their role in cellular patterning and leaf development is mostly limited to flowering plants. Our phylogenetic analysis identified three SCR homologs in moss Physcomitrium patens, amongst which PpSCR1 showed highest expression in gametophores and its promoter activity was prominent at the mid-vein and the flanking leaf blade cells pointing towards its role in leaf development. Notably, out of the three SCR homologs, only the ppscr1 knock-out lines developed slender leaves with four times narrower leaf blade and three times thicker mid-vein. Detailed histology studies revealed that slender leaf phenotype is either due to the loss of anticlinal cell divisions or failure of periclinal division suppression in the leaf blade. RNA-Seq analyses revealed that genes responsible for cell division and differentiation are expressed differentially in the mutant. PpSCR1 overexpression lines exhibited significantly wider leaf lamina, further reconfirming the role in leaf development. Together, our data suggests that PpSCR1 is involved in the leaf blade and mid-vein development of moss and that its role in the regulation of cell division and proliferation is ancient and conserved among flowering plants and mosses.

Divergence of trafficking and polarization mechanisms for PIN auxin transporters during land plant evolution

The phytohormone auxin, and its directional transport through tissues, plays a fundamental role in the development of higher plants. This polar auxin transport predominantly relies on PIN-FORMED (PIN) auxin exporters. Hence, PIN polarization is crucial for development, but its evolution during the rise of morphological complexity in land plants remains unclear. Here, we performed a cross-species investigation by observing the trafficking and localization of endogenous and exogenous PINs in two bryophytes, Physcomitrium patens and Marchantia polymorpha, and in the flowering plant Arabidopsis thaliana. We confirmed that the GFP fusion did not compromise the auxin export function of all examined PINs by using a radioactive auxin export assay and by observing the phenotypic changes in transgenic bryophytes. Endogenous PINs polarize to filamentous apices, while exogenous Arabidopsis PINs distribute symmetrically on the membrane in both bryophytes. In the Arabidopsis root epidermis, bryophytic PINs have no defined polarity. Pharmacological interference revealed a strong cytoskeletal dependence of bryophytic but not Arabidopsis PIN polarization. The divergence of PIN polarization and trafficking is also observed within the bryophyte clade and between tissues of individual species. These results collectively reveal the divergence of PIN trafficking and polarity mechanisms throughout land plant evolution and the co-evolution of PIN sequence-based and cell-based polarity mechanisms.

RAF-like protein kinases mediate a deeply conserved, rapid auxin response

The plant-signaling molecule auxin triggers fast and slow cellular responses across land plants and algae. The nuclear auxin pathway mediates gene expression and controls growth and development in land plants, but this pathway is absent from algal sister groups. Several components of rapid responses have been identified in Arabidopsis, but it is unknown if these are part of a conserved mechanism. We recently identified a fast, proteome-wide phosphorylation response to auxin. Here, we show that this response occurs across 5 land plant and algal species and converges on a core group of shared targets. We found conserved rapid physiological responses to auxin in the same species and identified rapidly accelerated fibrosarcoma (RAF)-like protein kinases as central mediators of auxin-triggered phosphorylation across species. Genetic analysis connects this kinase to both auxin-triggered protein phosphorylation and rapid cellular response, thus identifying an ancient mechanism for fast auxin responses in the green lineage.

An ABCB transporter regulates anisotropic cell expansion via cuticle deposition in the moss Physcomitrium patens

Anisotropic cell expansion is crucial for the morphogenesis of land plants, as cell migration is restricted by the rigid cell wall. The anisotropy of cell expansion is regulated by mechanisms acting on the deposition or modification of cell wall polysaccharides. Besides the polysaccharide components in the cell wall, a layer of hydrophobic cuticle covers the outer cell wall and is subjected to tensile stress that mechanically restricts cell expansion. However, the molecular machinery that deposits cuticle materials in the appropriate spatiotemporal manner to accommodate cell and tissue expansion remains elusive. Here, we report that PpABCB14, an ATP-binding cassette transporter in the moss Physcomitrium patens, regulates the anisotropy of cell expansion. PpABCB14 localized to expanding regions of leaf cells. Deletion of PpABCB14 resulted in impaired anisotropic cell expansion. Unexpectedly, the cuticle proper was reduced in the mutants, and the cuticular lipid components decreased. Moreover, induced PpABCB14 expression resulted in deformed leaf cells with increased cuticle lipid accumulation on the cell surface. Taken together, PpABCB14 regulates the anisotropy of cell expansion via cuticle deposition, revealing a regulatory mechanism for cell expansion in addition to the mechanisms acting on cell wall polysaccharides.

Peptides from conserved tandem direct repeats of SHORT-LEAF regulate gametophore development in moss P. patens

Tandem direct repeat (TDR)-containing proteins, present across all domains of life, play crucial roles in plant development and defense mechanisms. Previously, we identified that disruption of a bryophyte-specific protein family, SHORT-LEAF (SHLF), possessing the longest reported TDRs, is the cause of the shlf mutant phenotype in Physcomitrium patens. shlf exhibits reduced apical dominance, altered auxin distribution, and 2-fold shorter leaves. However, the molecular role of SHLF was unclear due to the absence of known conserved domains. Through a series of protein domain deletion analyses, here, we demonstrate the importance of the signal peptide and the conserved TDRs and report a minimal functional protein (miniSHLF) containing the N-terminal signal peptide and first two TDRs (N-TDR1-2). We also demonstrate that SHLF behaves as a secretory protein and that the TDRs contribute to a pool of secreted peptides essential for SHLF function. Further, we identified that the mutant secretome lacks SHLF peptides, which are abundant in WT and miniSHLF secretomes. Interestingly, shlf mutants supplemented with the secretome or peptidome from WT or miniSHLF showed complete or partial phenotypic recovery. Transcriptomic and metabolomic analyses revealed that shlf displays an elevated stress response, including high ROS activity and differential accumulation of genes and metabolites involved in the phenylpropanoid pathway, which may affect auxin distribution. The TDR-specific synthetic peptide SHLFpep3 (INIINAPLQGFKIA) also rescued the mutant phenotypes, including the altered auxin distribution, in a dosage-dependent manner and restored the mutant's stress levels. Our study shows that secretory SHLF peptides derived from conserved TDRs regulate moss gametophore development.

Adaptive Optics Microscopy with Wavefront Sensing Based on Neighbor Correlation

Complex structures in living cells and tissues induce wavefront errors when light waves pass through them, and images observed with optical microscopes are undesirably blurred. This problem is especially serious for living plant cells because images are strikingly degraded even within a single cell. Adaptive optics (AO) is expected to be a solution to this problem by correcting such wavefront errors, thus enabling high-resolution imaging. In particular, scene-based AO involves wavefront sensing based on the image correlation between subapertures in a Shack-Hartmann wavefront sensor and thus does not require an intense point light source. However, the complex 3D structures of living cells often cause low correlation between subimages, leading to loss of accuracy in wavefront sensing. This paper proposes a novel method for scene-based sensing using only image correlations between adjacent subapertures. The method can minimize changes between subimages to be correlated and thus prevent inaccuracy in phase estimation. Using an artificial test target mimicking the optical properties of a layer of living plant cells, an imaging performance with a Strehl ratio of approximately 0.5 was confirmed. Upon observation of chloroplast autofluorescence inside living leaf cells of the moss Physcomitrium patens, recovered resolution images were successfully obtained even with complex biological structures. Under bright-field illumination, the proposed method outperformed the conventional method, demonstrating the future potential of this method for label- and damage-free AO microscopy. Several points for improvement in terms of the effect of AO correction are discussed.

Variations in the enzymatic activity of S1-type nucleases results from differences in their active site structures

BACKGROUND

S1-like nucleases are widespread enzymes commonly used in biotechnology and molecular biology. Although it is commonly believed that they are mainly Zn2+-dependent acidic enzymes, we have found that numerous members of this family deviate from this rule. Therefore, in this work, we decided to check how broad is the range of non‑zinc-dependent S1-like nucleases and what is the molecular basis of their activities.

METHODS

S1-like nucleases chosen for analysis were achieved through heterologous expression in appropriate eukaryotic hosts. To characterize nucleases' active-site properties, point mutations were introduced in selected positions. The enzymatic activities of wild-type and mutant nucleases were tested by in-gel nuclease activity assay.

RESULTS

We discovered that S1-like nucleases encoded by non-vascular plants and single-celled protozoa, like their higher plant homologues, exhibit a large variety of catalytic properties. We have shown that these individual properties are determined by specific non-conserved active site residues.

CONCLUSIONS

Our findings demonstrate that mutations that occur during evolution can significantly alter the catalytic properties of S1-like nucleases. As a result, different ions can compete for particular S1-type nucleases' active sites. This phenomenon undermines the existing classification of S1-like nucleases.

GENERAL SIGNIFICANCE

Our findings have numerous implications for applications and understanding the S1-like nucleases' biological functions. For example, new biotechnological applications should take into account their unexpected catalytic properties. Moreover, these results demonstrate that the trinuclear zinc-based model commonly used to characterize the catalytic activities of S1-like nucleases is insufficient to explain the actions of non‑zinc-dependent members of this family.

Protocol: an improved method for inducing sporophyte generation in the model moss Physcomitrium patens under nitrogen starvation

BACKGROUND

Land plants exhibit a haplodiplontic life cycle, whereby multicellular bodies develop in both the haploid and diploid generations. The early-diverging land plants, known as bryophytes, have a haploid-dominant life cycle, in which a short-lived multicellular body in the diploid generation, known as the sporophyte, develops on the maternal haploid gametophyte tissues. The moss Physcomitrium (Physcomitrella) patens has become one of the most powerful model systems in evolutionary plant developmental studies. To induce diploid sporophytes of P. patens, several protocols are implemented. One of the conventional approaches is to grow approximately one-month-old gametophores for another month on Jiffy-7 pellets made from the peat moss that is difficult to fully sterilize. A more efficient method to obtain all tissues throughout the life cycle should accelerate studies of P. patens.

RESULTS

Here, we investigated the effect of nitrogen conditions on the growth and development of P. patens. We provide an improved protocol for the sporophyte induction of P. patens using a BCD-based solid culture medium without Jiffy-7 pellets, based on the finding that the formation of gametangia and subsequent sporophytes is promoted by nitrogen-free growth conditions. The protocol consists of two steps; first, culture the protonemata and gametophores on nitrogen-rich medium under continuous light at 25 °C, and then transfer the gametophores onto nitrogen-free medium under short-day and at 15 °C for sporophyte induction. The protocol enables to shorten the induction period and reduce the culture space.

CONCLUSIONS

Our more efficient and shortened protocol for inducing the formation of sporophytes will contribute to future studies into the fertilization or the diploid sporophyte generation of P. patens.

Kinesin-4 optimizes microtubule orientations for responsive tip growth guidance in moss

Tip-growing cells of, amongst others, plants and fungi secrete wall materials in a highly polarized fashion for fast and efficient colonization of the environment. A polarized microtubule cytoskeleton, in which most microtubule ends are directed toward the growing apex, has been implicated in directing growth. Its organizing principles, in particular regarding maintenance of network unipolarity, have remained elusive. We show that a kinesin-4 protein, hitherto best known for a role in cytokinesis, suppresses encounters between antiparallel microtubules. Without this activity, microtubules hyper-aligned along the growth axis and increasingly grew away from the apex. Cells themselves displayed an overly straight growth path and a delayed gravitropic response. This result revealed conflicting systemic needs for a stable growth direction and an ability to change course in response to extracellular cues. Thus, the use of selective inhibition of microtubule growth at antiparallel overlaps constitutes a new organizing principle within a unipolar microtubule array.

The TOPLESS corepressor regulates developmental switches in the bryophyte Physcomitrium patens that were critical for plant terrestrialisation

The plant-specific TOPLESS (TPL) family of transcriptional corepressors is integral to multiple angiosperm developmental processes. Despite this, we know little about TPL function in other plants. To address this gap, we investigated the roles TPL plays in the bryophyte Physcomitrium patens, which diverged from angiosperms approximately 0.5 billion years ago. Although complete loss of PpTPL function is lethal, transgenic lines with reduced PpTPL activity revealed that PpTPLs are essential for two fundamental developmental switches in this plant: the transitions from basal photosynthetic filaments (chloronemata) to specialised foraging filaments (caulonemata) and from two-dimensional (2D) to three-dimensional (3D) growth. Using a transcriptomics approach, we integrated PpTPL into the regulatory network governing 3D growth and we propose that PpTPLs represent another important class of regulators that are essential for the 2D-to-3D developmental switch. Transcriptomics also revealed a previously unknown role for PpTPL in the regulation of flavonoids. Intriguingly, 3D growth and the formation of caulonemata were crucial innovations that facilitated the colonisation of land by plants, a major transformative event in the history of life on Earth. We conclude that TPL, which existed before the land plants, was co-opted into new developmental pathways, enabling phytoterrestrialisation and the evolution of land plants.

Properties of INDETERMINATE DOMAIN Proteins from Physcomitrium patens: DNA-Binding, Interaction with GRAS Proteins, and Transcriptional Activity

INDETERMINATE DOMAIN (IDD) proteins are plant-specific transcription factors that interact with GRAS proteins, such as DELLA and SHORT ROOT (SHR), to regulate target genes. The combination of IDD and DELLA proteins regulates genes involved in gibberellic acid (GA) synthesis and GA signaling, whereas the combination of IDD with the complex of SHR and SCARECROW, another GRAS protein, regulates genes involved in root tissue formation. Previous bioinformatic research identified seven IDDs, two DELLA, and two SHR genes in Physcomitrium patens, a model organism for non-vascular plants (bryophytes), which lack a GA signaling pathway and roots. In this study, DNA-binding properties and protein-protein interaction of IDDs from P. patens (PpIDD) were analyzed. Our results showed that the DNA-binding properties of PpIDDs were largely conserved between moss and seed plants. Four PpIDDs showed interaction with Arabidopsis DELLA (AtDELLA) proteins but not with PpDELLAs, and one PpIDD showed interaction with PpSHR but not with AtSHR. Moreover, AtIDD10 (JACKDAW) interacted with PpSHR but not with PpDELLAs. Our results indicate that DELLA proteins have modified their structure to interact with IDD proteins during evolution from moss lineage to seed plants, whereas the interaction of IDD and SHR was already present in moss lineage.

Parallel evolution of two AIM24 protein subfamilies and their conserved functions in ER stress tolerance in land plants

Despite decades of efforts in genome sequencing and functional characterization, some important protein families remain poorly understood. In this study, we report the classification, evolution, and functions of the largely uncharacterized AIM24 protein family in plants, including the identification of a novel subfamily. We show that two AIM24 subfamilies (AIM24-A and AIM24-B) are commonly distributed in major plant groups. These two subfamilies not only have modest sequence similarities and different gene structures but also are of independent bacterial ancestry. We performed comparative functional investigations on the two AIM24 subfamilies using three model plants: the moss Physcomitrium patens, the liverwort Marchantia polymorpha, and the flowering plant Arabidopsis thaliana. Intriguingly, despite their significant differences in sequence and gene structure, both AIM24 subfamilies are involved in ER stress tolerance and the unfolded protein response (UPR). In addition, transformation of the AIM24-A gene from P. patens into the AIM24-B null mutant of A. thaliana could at least partially rescue ER stress tolerance and the UPR. We also discuss the role of AIM24 genes in plant development and other cellular activities. This study provides a unique example of parallel evolution in molecular functions and can serve as a foundation for further investigation of the AIM24 family in plants.

Overexpression of Physcomitrium patens cell cycle regulators leads to larger gametophytes

Regulation of cell division is crucial for the development of multicellular organisms, and in plants, this is in part regulated by the D-type cyclins (CYCD) and cyclin-dependent kinase A (CDKA) complex. Cell division regulation in Physcomitrium differs from other plants, by having cell division checks at both the G1 to S and G2 to M transition, controlled by the CYCD1/CDKA2 and CYCD2/CDKA1 complexes, respectively. This led us to hypothesize that upregulation of cell division could be archived in Bryophytes, without the devastating phenotypes observed in Arabidopsis. Overexpressing lines of PpCYCD1, PpCYCD2, PpCDKA1, or PpCDKA2 under Ubiquitin promotor control provided transcriptomic and phenotypical data that confirmed their involvement in the G1 to S or G2 to M transition control. Interestingly, combinatorial overexpression of all four genes produced plants with dominant PpCDKA2 and PpCYCD1 phenotypes and led to plants with twice as large gametophores. No detrimental phenotypes were observed in this line and two of the major carbon sinks in plants, the cell wall and starch, were unaffected by the increased growth rate. These results show that the cell cycle characteristics of P. patens can be manipulated by the ectopic expression of cell cycle regulators.

Using Spinning Disk Microscopy to Observe the Mitotic and Cytokinetic Apparatus in Physcomitrium patens

Protonemata of the moss Physcomitrium patens are ideal structures in which to observe cytoskeletal organization and dynamics. Special care is needed to prepare P. patens cultures for high-resolution microscopy. Here, we describe methods for spinning disk microscopy of dividing P. patens cells expressing sGFP-tubulin and H2B-mCherry, including detailed methods for culturing P. patens.

Horizontally acquired fungal killer protein genes affect cell development in mosses

The moss Physcomitrium patens is crucial for studying plant development and evolution. Although the P. patens genome includes genes acquired from bacteria, fungi and viruses, the functions and evolutionary significance of these acquired genes remain largely unclear. Killer protein 4 (KP4) is a toxin secreted by the phytopathogenic fungus Ustilago maydis that inhibits the growth of sensitive target strains by blocking their calcium uptake. Here, we show that KP4 genes in mosses were acquired from fungi through at least three independent events of horizontal gene transfer. Two paralogous copies of KP4 (PpKP4-1 and PpKP4-2) exist in P. patens. Knockout mutants ppkp4-1 and ppkp4-2 showed cell death at the protonemal stage, and ppkp4-2 also exhibited defects in tip growth. We provide experimental evidence indicating that PpKP4-1/2 affects P. patens protonemal cell development by mediating cytoplasmic calcium and that KP4 genes are functionally conserved between P. patens and fungi. The present study provides additional insights into the role of horizontal gene transfer in land plant development and evolution.

Knockout of DDM1 in Physcomitrium patens disrupts DNA methylation with a minute effect on transposon regulation and development

The Snf2 chromatin remodeler, DECREASE IN DNA METHYLATION 1 (DDM1) facilitates DNA methylation. In flowering plants, DDM1 mediates methylation in heterochromatin, which is targeted primarily by MET1 and CMT methylases and is necessary for silencing transposons and for proper development. DNA methylation mechanisms evolved throughout plant evolution, whereas the role of DDM1 in early terrestrial plants remains elusive. Here, we studied the function of DDM1 in the moss, Physcomitrium (Physcomitrella) patens, which has robust DNA methylation that suppresses transposons and is mediated by a MET1, a CMT, and a DNMT3 methylases. To elucidate the role of DDM1 in P. patens, we have generated a knockout mutant and found DNA methylation to be strongly disrupted at any of its sequence contexts. Symmetric CG and CHG sequences were affected stronger than asymmetric CHH sites. Furthermore, despite their separate targeting mechanisms, CG (MET) and CHG (CMT) methylation were similarly depleted by about 75%. CHH (DNMT3) methylation was overall reduced by about 25%, with an evident hyper-methylation activity within lowly-methylated euchromatic transposon sequences. Despite the strong hypomethylation effect, only a minute number of transposons were transcriptionally activated in Ppddm1. Finally, Ppddm1 was found to develop normally throughout the plant life cycle. These results demonstrate that DNA methylation is strongly dependent on DDM1 in a non-flowering plant; that DDM1 is required for plant-DNMT3 (CHH) methylases, though to a lower extent than for MET1 and CMT enzymes; and that distinct and separate methylation pathways (e.g. MET1-CG and CMT-CHG), can be equally regulated by the chromatin and that DDM1 plays a role in it. Finally, our data suggest that the biological significance of DDM1 in terms of transposon regulation and plant development, is species dependent.

Surface-localized glycoproteins act through class C ARFs to fine-tune gametophore initiation in Physcomitrium patens

Arabinogalactan proteins are functionally diverse cell wall structural glycoproteins that have been implicated in cell wall remodeling, although the mechanistic actions remain elusive. Here, we identify and characterize two AGP glycoproteins, SLEEPING BEAUTY (SB) and SB-like (SBL), that negatively regulate the gametophore bud initiation in Physcomitrium patens by dampening cell wall loosening/softening. Disruption of SB and SBL led to accelerated gametophore formation and altered cell wall compositions. The function of SB is glycosylation dependent and genetically connected with the class C auxin response factor (ARF) transcription factors PpARFC1B and PpARFC2. Transcriptomics profiling showed that SB upregulates PpARFC2, which in turn suppresses a range of cell wall-modifying genes that are required for cell wall loosening/softening. We further show that PpARFC2 binds directly to multiple AuxRE motifs on the cis-regulatory sequences of PECTIN METHYLESTERASE to suppress its expression. Hence, our results demonstrate a mechanism by which the SB modulates the strength of intracellular auxin signaling output, which is necessary to fine-tune the timing of gametophore initials formation.

Functional conservation of an AP2/ERF transcription factor in cuticle formation suggests an important role in the terrestrialization of early land plants

The formation of a hydrophobic cuticle layer on aerial plant parts was a critical innovation for protection from the terrestrial environment during the evolution of land plants. However, little is known about the molecular mechanisms underlying cuticle biogenesis in early terrestrial plants. Here, we report an APETALA2/Ethylene Response Factor (AP2/ERF) transcriptional activator, PpWIN1, involved in cutin and cuticular wax biosynthesis in Physcomitrium patens and Arabidopsis. The transcript levels of PpWIN1 were 2.5-fold higher in gametophores than in the protonema, and increased by approximately 3- to 4.7-fold in the protonema and gametophores under salt and osmotic stresses. PpWIN1 harbouring transcriptional activation activity is localized in the nucleus of tobacco leaf epidermal cells. Δppwin1 knockout mutants displayed a permeable cuticle, increased water loss, and cutin- and wax-deficient phenotypes. In contrast, increased total cutin and wax loads, and decreased water loss rates were observed in PpWIN1-overexpressing Arabidopsis plants. The transcript levels of genes involved in cutin or wax biosynthesis were significantly up-regulated in PpWIN1-overexpressing Arabidopsis lines, indicating that PpWIN1 acts as a transcriptional activator in cuticle biosynthesis. This study suggests that Arabidopsis WIN1/SHN1 orthologs may be functionally conserved from early to vascular land plants.

Three-dimensionally visualized rhizoid system of moss, Physcomitrium patens, by refraction-contrast X-ray micro-computed tomography

Land plants have two types of shoot-supporting systems, root system and rhizoid system, in vascular plants and bryophytes. However, since the evolutionary origin of the systems is different, how much they exploit common systems or distinct systems to architect their structures is largely unknown. To understand the regulatory mechanism of how bryophytes architect the rhizoid system responding to environmental factors, we have developed the methodology to visualize and quantitatively analyze the rhizoid system of the moss, Physcomitrium patens, in 3D. The rhizoids having a diameter of 21.3 µm on the average were visualized by refraction-contrast X-ray micro-computed tomography using coherent X-ray optics available at synchrotron radiation facility SPring-8. Three types of shape (ring-shape, line and black circle) observed in tomographic slices of specimens embedded in paraffin were confirmed to be the rhizoids by optical and electron microscopy. Comprehensive automatic segmentation of the rhizoids, which appeared in three different form types in tomograms, was tested by a method using a Canny edge detector or machine learning. The accuracy of output images was evaluated by comparing with the manually segmented ground truth images using measures such as F1 score and Intersection over Union, revealing that the automatic segmentation using machine learning was more effective than that using the Canny edge detector. Thus, machine learning-based skeletonized 3D model revealed quite dense distribution of rhizoids. We successfully visualized the moss rhizoid system in 3D for the first time.

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

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

Phylogenetic distribution and expression pattern analyses identified a divergent basal body assembly protein involved in land plant spermatogenesis

Land plant spermatozoids commonly possess characteristic structures such as the spline, which consists of a microtubule array, the multilayered structure (MLS) in which the uppermost layer is a continuum of the spline, and multiple flagella. However, the molecular mechanisms underpinning spermatogenesis remain to be elucidated. We successfully identified candidate genes involved in spermatogenesis, deeply divergent BLD10s, by computational analyses combining multiple methods and omics data. We then examined the functions of BLD10s in the liverwort Marchantia polymorpha and the moss Physcomitrium patens. MpBLD10 and PpBLD10 are required for normal basal body (BB) and flagella formation. Mpbld10 mutants exhibited defects in remodeling of the cytoplasm and nucleus during spermatozoid formation, and thus MpBLD10 should be involved in chromatin reorganization and elimination of the cytoplasm during spermiogenesis. We identified orthologs of MpBLD10 and PpBLD10 in diverse Streptophyta and found that MpBLD10 and PpBLD10 are orthologous to BLD10/CEP135 family proteins, which function in BB assembly. However, BLD10s evolved especially quickly in land plants and MpBLD10 might have acquired additional functions in spermatozoid formation through rapid molecular evolution.

An Overview of Pentatricopeptide Repeat (PPR) Proteins in the Moss Physcomitrium patens and Their Role in Organellar Gene Expression

Pentatricopeptide repeat (PPR) proteins are one type of helical repeat protein that are widespread in eukaryotes. In particular, there are several hundred PPR members in flowering plants. The majority of PPR proteins are localized in the plastids and mitochondria, where they play a crucial role in various aspects of RNA metabolism at the post-transcriptional and translational steps during gene expression. Among the early land plants, the moss Physcomitrium (formerly Physcomitrella) patens has at least 107 PPR protein-encoding genes, but most of their functions remain unclear. To elucidate the functions of PPR proteins, a reverse-genetics approach has been applied to P. patens. To date, the molecular functions of 22 PPR proteins were identified as essential factors required for either mRNA processing and stabilization, RNA splicing, or RNA editing. This review examines the P. patens PPR gene family and their current functional characterization. Similarities and a diversity of functions of PPR proteins between P. patens and flowering plants and their roles in the post-transcriptional regulation of organellar gene expression are discussed.

Hydrogen Peroxide Mediates Premature Senescence Caused by Darkness and Inorganic Nitrogen Starvation in Physcomitrium patens

Leaf senescence accompanied by yellowing and Rubisco degradation occurs prematurely in response to various stresses. However, signaling pathways between stress perception and senescence responses are not understood fully, although previous studies suggest the involvement of reactive oxygen species (ROS). While investigating the physiological functions of autophagy in Physcomitrium patens using wild-type (WT) and autophagy-deficient atg5 strains, we found that Physcomitrium colonies senesce prematurely under dark or nitrogen-deficient conditions, with atg5 senescing earlier than WT. In the present study, we measured cellular H2O2, and examined whether H2O2 mediates premature senescence in Physcomitrium colonies. Methyl viologen, an ROS generator, increased cellular H2O2 levels and caused senescence-like symptoms. H2O2 levels were also elevated to the same plateau levels in WT and atg5 under dark or nitrogen-deficient conditions. The ROS scavenger N-acetylcysteine and the ROS source inhibitor carbonyl cyanide m-chlorophenylhydrazone inhibited the increase in H2O2 levels as well as senescence. Upon transfer to a nitrogen-deficient medium, H2O2 levels increased earlier in atg5 than in WT by ~18 h, whereas atg5 yellowed earlier by >2 days. We conclude that the increased H2O2 levels under dark or nitrogen-deficient conditions mediate premature senescence in Physcomitrium but do not explain the different senescence responses of WT and atg5 cells.

Amino Acids Supplied through the Autophagy/Endocytosis Pathway Promote Starch Synthesis in Physcomitrella Protonemal Cells

The physiological implications of autophagy in plant cells have not been fully elucidated. Therefore, we investigated the consequences of autophagy in the moss Physcomitrella by measuring biochemical parameters (fresh and dry weights; starch, amino acid, carbohydrate, and NH3 content) in wild-type (WT) and autophagy-deficient atg5 Physcomitrella cells. We found higher starch levels and a higher net starch synthesis rate in WT cells than in atg5 cells cultured in a glucose-containing culture medium, whereas net starch degradation was similar in the two strains cultured in a glucose-deficient culture medium. Additionally, the treatment of cells with the autophagy inhibitor 3-methyladenine suppressed starch synthesis. Loading bovine serum albumin into atg5 cells through endocytosis, i.e., supplying proteins to vacuoles in the same way as through autophagy, accelerated starch synthesis, whereas loading glutamine through the plasma membrane had no such effect, suggesting that Physcomitrella cells distinguish between different amino acid supply pathways. After net starch synthesis, NH3 levels increased in WT cells, although the change in total amino acid content did not differ between WT and atg5 cells, indicating that autophagy-produced amino acids are oxidized rapidly. We conclude that autophagy promotes starch synthesis in Physcomitrella by supplying the energy obtained by oxidizing autophagy-produced amino acids.

Molecular and physiological responses to desiccation indicate the abscisic acid pathway is conserved in the peat moss, Sphagnum

Mosses of the genus Sphagnum are the main components of peatlands, a major carbon-storing ecosystem. Changes in precipitation patterns are predicted to affect water relations in this ecosystem, but the effect of desiccation on the physiological and molecular processes in Sphagnum is still largely unexplored. Here we show that different Sphagnum species have differential physiological and molecular responses to desiccation but, surprisingly, this is not directly correlated with their position in relation to the water table. In addition, the expression of drought responsive genes is increased upon water withdrawal in all species. This increase in gene expression is accompanied by an increase in abscisic acid (ABA), supporting a role for ABA during desiccation responses in Sphagnum. Not only do ABA levels increase upon desiccation, but Sphagnum plants pre-treated with ABA display increased tolerance to desiccation, suggesting that ABA levels play a functional role in the response. In addition, many of the ABA signalling components are present in Sphagnum and we demonstrate, by complementation in Physcomitrium patens, that Sphagnum ABI3 is functionally conserved. The data presented here, therefore, support a conserved role for ABA in desiccation responses in Sphagnum.

An ancestral function of strigolactones as symbiotic rhizosphere signals

In flowering plants, strigolactones (SLs) have dual functions as hormones that regulate growth and development, and as rhizosphere signaling molecules that induce symbiosis with arbuscular mycorrhizal (AM) fungi. Here, we report the identification of bryosymbiol (BSB), an SL from the bryophyte Marchantia paleacea. BSB is also found in vascular plants, indicating its origin in the common ancestor of land plants. BSB synthesis is enhanced at AM symbiosis permissive conditions and BSB deficient mutants are impaired in AM symbiosis. In contrast, the absence of BSB synthesis has little effect on the growth and gene expression. We show that the introduction of the SL receptor of Arabidopsis renders M. paleacea cells BSB-responsive. These results suggest that BSB is not perceived by M. paleacea cells due to the lack of cognate SL receptors. We propose that SLs originated as AM symbiosis-inducing rhizosphere signaling molecules and were later recruited as plant hormone.

An ancestral function of strigolactones as symbiotic rhizosphere signals

In flowering plants, strigolactones (SLs) have dual functions as hormones that regulate growth and development, and as rhizosphere signaling molecules that induce symbiosis with arbuscular mycorrhizal (AM) fungi. Here, we report the identification of bryosymbiol (BSB), an SL from the bryophyte Marchantia paleacea. BSB is also found in vascular plants, indicating its origin in the common ancestor of land plants. BSB synthesis is enhanced at AM symbiosis permissive conditions and BSB deficient mutants are impaired in AM symbiosis. In contrast, the absence of BSB synthesis has little effect on the growth and gene expression. We show that the introduction of the SL receptor of Arabidopsis renders M. paleacea cells BSB-responsive. These results suggest that BSB is not perceived by M. paleacea cells due to the lack of cognate SL receptors. We propose that SLs originated as AM symbiosis-inducing rhizosphere signaling molecules and were later recruited as plant hormone.

Strigolactones (SLs) regulate angiosperm development and promote symbiosis with arbuscular mycorrhizae. Here the authors show that bryosymbiol, an SL present in bryophytes and angiosperms, promotes AM symbiosis in Marchantia paleacea suggesting an ancestral function of SLs as rhizosphere signals.

The cellular function of ROP GTPase prenylation is important for multicellularity in the moss Physcomitrium patens

A complete picture of how signaling pathways lead to multicellularity is largely unknown. Previously, we generated mutations in a protein prenylation enzyme, GGB, and showed that it is essential for maintaining multicellularity in the moss Physcomitrium patens. Here, we show that ROP GTPases act as downstream factors that are prenylated by GGB and themselves play an important role in the multicellularity of P. patens. We also show that the loss of multicellularity caused by the suppression of GGB or ROP GTPases is due to uncoordinated cell expansion, defects in cell wall integrity and the disturbance of the directional control of cell plate orientation. Expressing prenylatable ROP in the ggb mutant not only rescues multicellularity in protonemata but also results in development of gametophores. Although the prenylation of ROP is important for multicellularity, a higher threshold of active ROP is required for gametophore development. Thus, our results suggest that ROP activation via prenylation by GGB is a key process at both cell and tissue levels, facilitating the developmental transition from one dimension to two dimensions and to three dimensions in P. patens.

Topoisomerase 1α is required for synchronous spermatogenesis in Physcomitrium patens

DNA topoisomerase 1 (TOP1) plays general roles in DNA replication and transcription by regulating DNA topology in land plants and metazoans. TOP1 is also involved in specific developmental events; however, whether TOP1 plays a conserved developmental role among multicellular organisms is unknown. Here, we investigated the developmental roles of TOP1 in the moss Physcomitrium (Physcomitrella) patens with gene targeting, microscopy, 3D image segmentation and crossing experiments. We discovered that the disruption of TOP1α, but not its paralogue TOP1β, leads to a defect in fertilisation and subsequent sporophyte formation in P. patens. In the top1α mutant, the egg cell was functional for fertilisation, while sperm cells were fewer and infertile with disordered structures. We observed that the nuclei volume of wild-type sperm cells synchronously decreases during antheridium development, indicating chromatin condensation towards the compact sperm head. By contrast, the top1α mutant exhibited attenuated cell divisions and asynchronous and defective contraction of the nuclei of sperm cells throughout spermatogenesis. These results indicate that TOP1α is involved in cell division and chromatin condensation during spermatogenesis in P. patens. Our results suggest that the regulation of DNA topology by TOP1 plays a key role in spermatogenesis in both land plants and metazoans.

Abscisic acid switches cell division modes of asymmetric cell division and symmetric cell division in stem cells of protonemal filaments in the moss Physcomitrium patens

Multicellular organisms regulate cell numbers and cell fate by using asymmetric cell division (ACD) and symmetric cell division (SCD) during their development and to adapt to unfavorable environmental conditions. A stem cell self-renews and generates differentiated cells. In plants, various types of cells are produced by ACD or SCD; however, the molecular mechanisms of ACD or SCD and the cell division mode switch are largely unknown. The moss Physcomitrium (Physcomitrella) patens is a suitable model to study plant stem cells due to its simple anatomy. Here, we report the cell division mode switch induced by abscisic acid (ABA) in P. patens. ABA is synthesized in response to abiotic stresses and induces round-shape cells, called brood cells, from cylindrical protonemal cells. Although two daughter cells with distinct sizes were produced by ACD in a protonemal stem cell on ABA-free media, the sizes of two daughter cells became similar with ABA treatment. Actin microfilaments were spatially localized on the apices of apical stem cells in protonemata on ABA-free media, but the polar accumulation was lost under the condition of ABA treatment. Moreover, ABA treatment conferred an identical cell fate to the daughter cells in terms of cell division activity. Collectively, the results indicate ABA may suppress the ACD characteristics but evoke SCD in cells. We also noticed that ABA-induced brood cells not only self-renewed but regenerated protonemal cells when ABA was removed from the media, suggesting that brood cells are novel stem cells that are induced by environmental signals in P. patens.

Characterisation of rapid alkalinisation factors in Physcomitrium patens reveals functional conservation in tip growth

Small signalling peptides are key molecules for cell-to-cell communications in plants. The cysteine-rich signalling peptide, rapid alkalinisation factors (RALFs) family are involved in diverse developmental and stress responses and have expanded considerably during land plant evolution, implying neofunctionalisations in the RALF family. However, the ancestral roles of RALFs when land plant first acquired them remain unknown. Here, we functionally characterised two of the three RALFs in bryophyte Physcomitrium patens using loss-of-function mutants, overexpressors, as well as fluorescent proteins tagged reporter lines. We showed that PpRALF1 and PpRALF2 have overlapping functions in promoting protonema tip growth and elongation, showing a homologous function as the Arabidopsis RALF1 in promoting root hair tip growth. Although both PpRALFs are secreted to the plasma membrane on which PpRALF1 symmetrically localised, PpRALF2 showed a polarised localisation at the growing tip. Notably, proteolytic cleavage of PpRALF1 is necessary for its function. Our data reveal a possible evolutionary origin of the RALF functions and suggest that functional divergence of RALFs is essential to drive complex morphogenesis and to facilitate other novel processes in land plants.

Physcomitrium patens pentatricopeptide repeat protein PpPPR_32 is involved in the accumulation of psaC mRNA encoding the iron sulfur protein of photosystem I

Pentatricopeptide repeat (PPR) proteins are involved in RNA metabolism and also play a role in posttranscriptional regulation during plant organellar gene expression. Although a hundred of PPR proteins exist in the moss Physcomitrium patens, their functions are not fully understood. Here, we report the function of P-class PPR protein PpPPR_32 in P. patens. A transient expression assay using green fluorescent protein demonstrated that the N-terminal region of PpPPR_32 functions as a chloroplast-targeting transit peptide, indicating that PpPPR_32 is localized in chloroplasts. PpPPR_32 knockout (KO) mutants grew autotrophically but with reduced protonema growth and the poor formation of photosystem I (PSI) complexes. Quantitative real-time reverse transcription-polymerase chain reaction and RNA gel blot hybridization analyses revealed a significant reduction in the transcript level of the psaC gene encoding the iron sulfur protein of PSI but no alteration to the transcript levels of other PSI genes. This suggests that PpPPR_32 is specifically involved in the expression level of the psaC gene. Our results indicate that PpPPR_32 is essential for the accumulation of psaC transcript and PSI complexes.

A PSTAIRE-type cyclin-dependent kinase controls light responses in land plants

Light is a critical signal perceived by plants to adapt their growth rate and direction. Although many signaling components have been studied, how plants respond to constantly fluctuating light remains underexplored. Here, we showed that in the moss Physcomitrium (Physcomitrella) patens, the PSTAIRE-type cyclin-dependent kinase PpCDKA is dispensable for growth. Instead, PpCDKA and its homolog in Arabidopsis thaliana control light-induced tropisms and chloroplast movements by probably influencing the cytoskeleton organization independently of the cell cycle. In addition, lower PpCDKA kinase activity was required to elicit light responses relative to cell cycle regulation. Thus, our study suggests that plant CDKAs may have been co-opted to control multiple light responses, and owing to the bistable switch properties of PSTAIRE-type CDKs, the noncanonical functions are widely conserved for eukaryotic environmental adaptation.

Tracking Intercellular Movement of Fluorescent Proteins in Bryophytes

An important approach to investigate intercellular connectivity via plasmodesmata is to visualize and track the movement of fluorescent proteins between cells. The intercellular connectivity is largely controlled by the size exclusion limit of the pores. Over the past few decades, the technique to observe and analyze intercellular movement of a fluorescent protein has been developed mainly in angiosperms such as Arabidopsis thaliana. We recently applied the corresponding system to track the intercellular movement of the fluorescent protein Dendra2 in the moss Physcomitrium (Physcomitrella) patens. The protonemal tissues are particularly suited for observation of the intercellular movement due to the simple organization. Here, we describe a protocol suitable for the analysis of Dendra2 movement between cells in P. patens.

Sensor histidine kinases mediate ABA and osmostress signaling in the moss Physcomitrium patens

To survive fluctuating water availability on land, terrestrial plants must be able to sense water stresses, such as drought and flooding. The plant hormone abscisic acid (ABA) and plant-specific SNF1-related protein kinase 2 (SnRK2) play key roles in plant osmostress responses. We recently reported that, in the moss Physcomitrium patens, ABA and osmostress-dependent SnRK2 activation requires phosphorylation by an upstream RAF-like kinase (ARK). This RAF/SnRK2 module is an evolutionarily conserved mechanism of osmostress signaling in land plants. Surprisingly, ARK is also an ortholog of Arabidopsis CONSTITUTIVE RESPONSE 1 (CTR1), which negatively regulates the ethylene-mediated submergence response of P. patens, indicating a nexus for cross-talk between the two signaling pathways that regulate responses to water availability. However, the mechanism through which the ARK/SnRK2 module is activated in response to water stress remains to be elucidated. Here, we show that a group of ethylene-receptor-related sensor histidine kinases (ETR-HKs) is essential for ABA and osmostress responses in P. patens. The intracellular kinase domain of an ETR-HK from P. patens physically interacts with ARK at the endoplasmic reticulum in planta. Moreover, HK disruptants lack ABA-dependent autophosphorylation of the critical serine residue in the activation loop of ARK, leading to loss of SnRK2 activation in response to ABA and osmostress. Collectively with the notion that ETR-HKs participate in submergence responses, our present data suggest that the HK/ARK module functions as an integration unit for environmental water availability to elicit optimized water stress responses in the moss P. patens.

Genes encoding lipid II flippase MurJ and peptidoglycan hydrolases are required for chloroplast division in the moss Physcomitrella patens

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.

Moss PPR-SMR protein PpPPR_64 influences the expression of a psaA-psaB-rps14 gene cluster and processing of the 23S-4.5S rRNA precursor in chloroplasts

Moss PPR-SMR protein PpPPR_64 is a pTAC2 homolog but is functionally distinct from pTAC2. PpPPR_64 is required for psaA gene expression and its function may have evolved in mosses. The pentatricopeptide repeat (PPR) proteins are key regulatory factors responsible for the control of plant organellar gene expression. A small subset of PPR proteins possess a C-terminal small MutS-related (SMR) domain and have diverse roles in plant organellar biogenesis. However, the function of PPR-SMR proteins is not fully understood. Here, we report the function of PPR-SMR protein PpPPR_64 in the moss Physcomitrium patens. Phylogenetic analysis indicated that PpPPR_64 belongs to the same clade as the Arabidopsis PPR-SMR protein pTAC2. PpPPR_64 knockout (KO) mutants grew autotrophically but with reduced protonemata growth and the poor formation of photosystems' antenna complexes. Quantitative reverse transcription-polymerase chain reaction and RNA gel blot hybridization analyses revealed a significant reduction in transcript levels of the psaA-psaB-rps14 gene cluster but no alteration to transcript levels of most photosynthesis- and non-photosynthesis-related genes. In addition, RNA processing of 23S-4.5S rRNA precursor was impaired in the PpPPR_64 KO mutants. This suggests that PpPPR_64 is specifically involved in the expression level of the psaA-psaB-rps14 gene and in processing of the 23S-4.5S rRNA precursor. Our results indicate that PpPPR_64 is functionally distinct from pTAC2 and is a novel PPR-SMR protein required for proper chloroplast biogenesis in P. patens.

The unique bryophyte-specific repeat-containing protein SHORT-LEAF regulates gametophore development in moss

Convergent evolution of shoot development across plant lineages has prompted numerous comparative genetic studies. Though functional conservation of gene networks governing flowering plant shoot development has been explored in bryophyte gametophore development, the role of bryophyte-specific genes remains unknown. Previously, we have reported Tnt1 insertional mutants of moss defective in gametophore development. Here, we report a mutant (short-leaf; shlf) having two-fold shorter leaves, reduced apical dominance, and low plasmodesmata frequency. UHPLC-MS/MS-based auxin quantification and analysis of soybean (Glycine max) auxin-responsive promoter (GH3:GUS) lines exhibited a striking differential auxin distribution pattern in the mutant gametophore. Whole-genome sequencing and functional characterization of candidate genes revealed that a novel bryophyte-specific gene (SHORT-LEAF; SHLF) is responsible for the shlf phenotype. SHLF represents a unique family of near-perfect tandem direct repeat (TDR)-containing proteins conserved only among mosses and liverworts, as evident from our phylogenetic analysis. Cross-complementation with a Marchantia homolog partially recovered the shlf phenotype, indicating possible functional specialization. The distinctive structure (longest known TDRs), absence of any known conserved domain, localization in the endoplasmic reticulum, and proteolytic cleavage pattern of SHLF imply its function in bryophyte-specific cellular mechanisms. This makes SHLF a potential candidate to study gametophore development and evolutionary adaptations of early land plants.

The d-mannose/l-galactose pathway is the dominant ascorbate biosynthetic route in the moss Physcomitrium patens

Ascorbate is an abundant and indispensable redox compound in plants. Genetic and biochemical studies have established the d-mannose/l-galactose (d-Man/l-Gal) pathway as the predominant ascorbate biosynthetic pathway in streptophytes, while the d-galacturonate (d-GalUA) pathway is found in prasinophytes and euglenoids. Based on the presence of the complete set of genes encoding enzymes involved in the d-Man/l-Gal pathway and an orthologous gene encoding aldonolactonase (ALase) - a key enzyme for the d-GalUA pathway - Physcomitrium patens may possess both pathways. Here, we have characterized the moss ALase as a functional lactonase and evaluated the ascorbate biosynthesis capability of the two pathways using knockout mutants. Physcomitrium patens expresses two ALase paralogs, namely PpALase1 and PpALase2. Kinetic analyses with recombinant enzymes indicated that PpALase1 is a functional enzyme catalyzing the conversion of l-galactonic acid to the final precursor l-galactono-1,4-lactone and that it also reacts with dehydroascorbate as a substrate. Interestingly, mutants lacking PpALase1 (Δal1) showed 1.2-fold higher total ascorbate content than the wild type, and their dehydroascorbate content was increased by 50% compared with that of the wild type. In contrast, the total ascorbate content of mutants lacking PpVTC2-1 (Δvtc2-1) or PpVTC2-2 (Δvtc2-2), which encode the rate-limiting enzyme GDP-l-Gal phosphorylase in the d-Man/l-Gal pathway, was markedly decreased to 46 and 17%, respectively, compared with that of the wild type. Taken together, the dominant ascorbate biosynthetic pathway in P. patens is the d-Man/l-Gal pathway, not the d-GalUA pathway, and PpALase1 may play a significant role in ascorbate metabolism by facilitating dehydroascorbate degradation rather than ascorbate biosynthesis.

Evolution of an assembly factor-based subunit contributed to a novel NDH-PSI supercomplex formation in chloroplasts

Chloroplast NADH dehydrogenase-like (NDH) complex is structurally related to mitochondrial Complex I and forms a supercomplex with two copies of Photosystem I (the NDH-PSI supercomplex) via linker proteins Lhca5 and Lhca6. The latter was acquired relatively recently in a common ancestor of angiosperms. Here we show that NDH-dependent Cyclic Electron Flow 5 (NDF5) is an NDH assembly factor in Arabidopsis. NDF5 initiates the assembly of NDH subunits (PnsB2 and PnsB3) and Lhca6, suggesting that they form a contact site with Lhca6. Our analysis of the NDF5 ortholog in Physcomitrella and angiosperm genomes reveals the subunit PnsB2 to be newly acquired via tandem gene duplication of NDF5 at some point in the evolution of angiosperms. Another Lhca6 contact subunit, PnsB3, has evolved from a protein unrelated to NDH. The structure of the largest photosynthetic electron transport chain complex has become more complicated by acquiring novel subunits and supercomplex formation with PSI.

The chloroplast NDH complex interacts with Photosystem I to form the NDH-PSI supercomplex. Here the authors show that Arabidopsis NDF5 shares a common ancestor with the NDH subunit PnsB2 and acts as an NDH assembly factor initiating the assembly of PnsB2 and the evolutionarily distinct PnsB3.

Red light-regulated interaction of Per-Arnt-Sim histidine kinases with partner histidine-containing phosphotransfer proteins in Physcomitrium patens

Multi-step phosphorelay (MSP) is a broadly distributed signaling system in organisms. In MSP, histidine kinases (HKs) receive various environmental signals and transmit them by autophosphorylation followed by phosphotransfer to partner histidine-containing phosphotransfer proteins (HPts). Previously, we reported that Per-Arnt-Sim (PAS) domain-containing HK1 (PHK1) and PHK2 of the moss Physcomitrium (Physcomitrella) patens repressed red light-induced protonema branching, a critical step in the moss life cycle. In plants, PHK homolog-encoding genes are conserved only in early-diverging lineages such as bryophytes and lycophytes. PHKs-mediated signaling machineries attract attention especially from an evolutionary viewpoint, but they remain uninvestigated. Here, we studied the P. patens PHKs focusing on their subcellular patterns of localization and interaction with HPts. Yeast two-hybrid analysis, a localization assay with a green fluorescent protein, and a bimolecular fluorescence complementation analysis together showed that PHKs are localized and interact with partner HPts mostly in the nucleus, as unprecedented features for plant HKs. Additionally, red light triggered the interactions between PHKs and HPts in the cytoplasm, and light co-repressed the expression of PHK1 and PHK2 as well as genes encoding their partner HPts. Our results emphasize the uniqueness of PHKs-mediated signaling machineries, and functional implications of this uniqueness are discussed.

Metabolic Control of Gametophore Shoot Formation through Arginine in the Moss Physcomitrium patens

Shoot formation is accompanied by active cell proliferation and expansion, requiring that metabolic state adapts to developmental control. Despite the importance of such metabolic reprogramming, it remains unclear how development and metabolism are integrated. Here, we show that disruption of ANGUSTIFOLIA3 orthologs (PpAN3s) compromises gametophore shoot formation in the moss Physcomitrium patens due to defective cell proliferation and expansion. Trans-omics analysis reveals that the downstream activity of PpAN3 is linked to arginine metabolism. Elevating arginine level by chemical treatment leads to stunted gametophores and causes Ppan3 mutant-like transcriptional changes in the wild-type plant. Furthermore, ectopic expression of AtAN3 from Arabidopsis thaliana ameliorates the defective arginine metabolism and promotes gametophore formation in Ppan3 mutants. Together, these findings indicate that arginine metabolism is a key pathway associated with gametophore formation and provide evolutionary insights into the establishment of the shoot system in land plants through the integration of developmental and metabolic processes.

Ultra-deep sequencing reveals dramatic alteration of organellar genomes in Physcomitrella patens due to biased asymmetric recombination

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.

Phytophthora infestans RXLR effector AVR1 disturbs the growth of Physcomitrium patens without affecting Sec5 localization

Plant pathogens often exploit a whole range of effectors to facilitate infection. The RXLR effector AVR1 produced by the oomycete plant pathogen Phytophthora infestans suppresses host defense by targeting Sec5. Sec5 is a subunit of the exocyst, a protein complex that is important for mediating polarized exocytosis during plant development and defense against pathogens. The mechanism by which AVR1 manipulates Sec5 functioning is unknown. In this study, we analyzed the effect of AVR1 on Sec5 localization and functioning in the moss Physcomitrium patens. P. patens has four Sec5 homologs. Two (PpSec5b and PpSec5d) were found to interact with AVR1 in yeast-two-hybrid assays while none of the four showed a positive interaction with AVR1ΔT, a truncated version of AVR1. In P. patens lines carrying β-estradiol inducible AVR1 or AVR1ΔT transgenes, expression of AVR1 or AVR1ΔT caused defects in the development of caulonemal protonema cells and abnormal morphology of chloronema cells. Similar phenotypes were observed in Sec5- or Sec6-silenced P. patens lines, suggesting that both AVR1 and AVR1ΔT affect exocyst functioning in P. patens. With respect to Sec5 localization we found no differences between β-estradiol-treated and untreated transgenic AVR1 lines. Sec5 localizes at the plasma membrane in growing caulonema cells, also during pathogen attack, and its subcellular localization is the same, with or without AVR1 in the vicinity.

Activation of SnRK2 by Raf-like kinase ARK represents a primary mechanism of ABA and abiotic stress responses

The Raf-like protein kinase abscisic acid (ABA) and abiotic stress-responsive Raf-like kinase (ARK) previously identified in the moss Physcomitrium (Physcomitrella) patens acts as an upstream regulator of subgroup III SNF1-related protein kinase2 (SnRK2), the key regulator of ABA and abiotic stress responses. However, the mechanisms underlying activation of ARK by ABA and abiotic stress for the regulation of SnRK2, including the role of ABA receptor-associated group A PP2C (PP2C-A), are not understood. We identified Ser1029 as the phosphorylation site in the activation loop of ARK, which provided a possible mechanism for regulation of its activity. Analysis of transgenic P. patens ark lines expressing ARK-GFP with Ser1029-to-Ala mutation indicated that this replacement causes reductions in ABA-induced gene expression, stress tolerance, and SnRK2 activity. Immunoblot analysis using an anti-phosphopeptide antibody indicated that ABA treatments rapidly stimulate Ser1029 phosphorylation in the wild type (WT). The phosphorylation profile of Ser1029 in ABA-hypersensitive ppabi1 lacking protein phosphatase 2C-A (PP2C-A) was similar to that in the WT, whereas little Ser1029 phosphorylation was observed in ABA-insensitive ark missense mutant lines. Furthermore, newly isolated ppabi1 ark lines showed ABA-insensitive phenotypes similar to those of ark lines. Therefore, ARK is a primary activator of SnRK2, preceding negative regulation by PP2C-A in bryophytes, which provides a prototype mechanism for ABA and abiotic stress responses in plants.