DIX Domain Polymerization Drives Assembly of Plant Cell Polarity Complexes.

Cell polarity is fundamental for tissue morphogenesis in multicellular organisms. Plants and animals evolved multicellularity independently, and it is unknown whether their polarity systems are derived from a single-celled ancestor. Planar polarity in animals is conferred by Wnt signaling, an ancient signaling pathway transduced by Dishevelled, which assembles signalosomes by dynamic head-to-tail DIX domain polymerization. In contrast, polarity-determining pathways in plants are elusive. We recently discovered Arabidopsis SOSEKI proteins, which exhibit polar localization throughout development. Here, we identify SOSEKI as ancient polar proteins across land plants. Concentration-dependent polymerization via a bona fide DIX domain allows these to recruit ANGUSTIFOLIA to polar sites, similar to the polymerization-dependent recruitment of signaling effectors by Dishevelled. Cross-kingdom domain swaps reveal functional equivalence of animal and plant DIX domains. We trace DIX domains to unicellular eukaryotes and thus show that DIX-dependent polymerization is an ancient mechanism conserved between kingdoms and central to polarity proteins.

SUPPRESSOR OF MAX2 1-LIKE (SMXL) homologs are MAX2-dependent repressors of Physcomitrium patens growth.

SUPPRESSOR OF MAX2 (SMAX)1-LIKE (SMXL) proteins are a plant-specific clade of type I HSP100/Clp-ATPases. SMXL genes are present in virtually all land plant genomes. However, they have mainly been studied in angiosperms. In Arabidopsis (Arabidopsis thaliana), 3 functional SMXL subclades have been identified: SMAX1/SMXL2, SMXL345, and SMXL678. Of these, 2 subclades ensure endogenous phytohormone signal transduction. SMAX1/SMXL2 proteins are involved in KAI2 ligand (KL) signaling, while SMXL678 proteins are involved in strigolactone (SL) signaling. Many questions remain regarding the mode of action of these proteins, as well as their ancestral roles. We addressed these questions by investigating the functions of the 4 SMXL genes in the moss Physcomitrium patens. We demonstrate that PpSMXL proteins are involved in the conserved ancestral MAX2-dependent KL signaling pathway and negatively regulate growth. However, PpSMXL proteins expressed in Arabidopsis cannot replace SMAX1 or SMXL2 function in KL signaling, whereas they can functionally replace SMXL4 and SMXL5 and restore root growth. Therefore, the molecular functions of SMXL proteins are conserved, but their interaction networks are not. Moreover, the PpSMXLC/D clade positively regulates SL signal transduction in P. patens. Overall, our data reveal that SMXL proteins in moss mediate crosstalk between the SL and KL signaling pathways.

Conservation of the moss RNA editing factor PPR78 despite the loss of its known cytidine-to-uridine editing sites is explained by a hidden extra target.

Cytidine (C)-to-uridine (U) RNA editing in plant organelles relies on specific RNA-binding pentatricopeptide repeat (PPR) proteins. In the moss Physcomitrium patens, all such RNA editing factors feature a C-terminal DYW domain that acts as the cytidine deaminase for C-to-U conversion. PPR78 of Physcomitrium targets 2 mitochondrial editing sites, cox1eU755SL and rps14eU137SL. Remarkably, the latter is edited to highly variable degrees in different mosses. Here, we aimed to unravel the coevolution of PPR78 and its 2 target sites in mosses. Heterologous complementation in a Physcomitrium knockout line revealed that the variable editing of rps14eU137SL depends on the PPR arrays of different PPR78 orthologues but not their C-terminal domains. Intriguingly, PPR78 has remained conserved despite the simultaneous loss of editing at both known targets among Hypnales (feather mosses), suggesting it serves an additional function. Using a recently established RNA editing assay in Escherichia coli, we confirmed site-specific RNA editing by PPR78 in the bacterium and identified 4 additional off-targets in the bacterial transcriptome. Based on conservation profiles, we predicted ccmFNeU1465RC as a candidate editing target of PPR78 in moss mitochondrial transcriptomes. We confirmed editing at this site in several mosses and verified that PPR78 targets ccmFNeU1465RC in the bacterial editing system, explaining the conservation and functional adaptation of PPR78 during moss evolution.

Phylogeny-linked occurrence of ribosome stalling on the mRNAs of Arabidopsis unfolded protein response factor bZIP60 orthologs in divergent plant species.

The bZIP60, XBP1 and HAC1 mRNAs encode transcription factors that mediate the unfolded protein response (UPR) in plants, animals and yeasts, respectively. Upon UPR, these mRNAs undergo unconventional cytoplasmic splicing on the endoplasmic reticulum (ER) to produce active transcription factors. Although cytoplasmic splicing is conserved, the ER targeting mechanism differs between XBP1 and HAC1. The ER targeting of HAC1 mRNA occurs before translation, whereas that of XBP1 mRNA involves a ribosome-nascent chain complex that is stalled when a hydrophobic peptide emerges from the ribosome; the corresponding mechanism is unknown for bZIP60. Here, we analyzed ribosome stalling on bZIP60 orthologs of plants. Using a cell-free translation system, we detected nascent peptide-mediated ribosome stalling during the translation elongation of the mRNAs of Arabidopsis, rice and Physcomitrium (moss) orthologs, and the termination-step stalling in the Selaginella (lycopod) ortholog, all of which occurred ∼50 amino acids downstream of a hydrophobic region. Transfection experiments showed that ribosome stalling contributes to cytoplasmic splicing in bZIP60u orthologs of Arabidopsis and Selaginella. In contrast, ribosome stalling was undetectable for liverwort, Klebsormidium (basal land plant), and green algae orthologs. This study highlights the evolutionary diversity of ribosome stalling and its contribution to ER targeting in plants.

Beyond a PPR-RNA recognition code: Many aspects matter for the multi-targeting properties of RNA editing factor PPR56.

The mitochondrial C-to-U RNA editing factor PPR56 of the moss Physcomitrium patens is an RNA-binding pentatricopeptide repeat protein equipped with a terminal DYW-type cytidine deaminase domain. Transferred into Escherichia coli, PPR56 works faithfully on its two native RNA editing targets, nad3eU230SL and nad4eU272SL, and also converts cytidines into uridines at over 100 off-targets in the bacterial transcriptome. Accordingly, PPR56 is attractive for detailed mechanistic studies in the heterologous bacterial setup, allowing for scoring differential RNA editing activities of many target and protein variants in reasonable time. Here, we report (i) on the effects of numerous individual and combined PPR56 protein and target modifications, (ii) on the spectrum of off-target C-to-U editing in the bacterial background transcriptome for PPR56 and two variants engineered for target re-direction and (iii) on combinations of targets in tandem or separately at the 5'- and 3'-ends of large mRNAs. The latter experimentation finds enhancement of RNA editing at weak targets in many cases, including cox3eU290SF as a new candidate mitogenome target. We conclude that C-to-U RNA editing can be much enhanced by transcript features also outside the region ultimately targeted by PPRs of a plant editing factor, possibly facilitated by its enrichment or scanning along transcripts.

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

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

Multilevel analysis between Physcomitrium patens and Mortierellaceae endophytes explores potential long-standing interaction among land plants and fungi.

The model moss species Physcomitrium patens has long been used for studying divergence of land plants spanning from bryophytes to angiosperms. In addition to its phylogenetic relationships, the limited number of differential tissues, and comparable morphology to the earliest embryophytes provide a system to represent basic plant architecture. Based on plant-fungal interactions today, it is hypothesized these kingdoms have a long-standing relationship, predating plant terrestrialization. Mortierellaceae have origins diverging from other land fungi paralleling bryophyte divergence, are related to arbuscular mycorrhizal fungi but are free-living, observed to interact with plants, and can be found in moss microbiomes globally. Due to their parallel origins, we assess here how two Mortierellaceae species, Linnemannia elongata and Benniella erionia, interact with P. patens in coculture. We also assess how Mollicute-related or Burkholderia-related endobacterial symbionts (MRE or BRE) of these fungi impact plant response. Coculture interactions are investigated through high-throughput phenomics, microscopy, RNA-sequencing, differential expression profiling, gene ontology enrichment, and comparisons among 99 other P. patens transcriptomic studies. Here we present new high-throughput approaches for measuring P. patens growth, identify novel expression of over 800 genes that are not expressed on traditional agar media, identify subtle interactions between P. patens and Mortierellaceae, and observe changes to plant-fungal interactions dependent on whether MRE or BRE are present. Our study provides insights into how plants and fungal partners may have interacted based on their communications observed today as well as identifying L. elongata and B. erionia as modern fungal endophytes with P. patens.

Genome sequence and cell biological toolbox of the highly regenerative, coenocytic green feather alga Bryopsis.

Green feather algae (Bryopsidales) undergo a unique life cycle in which a single cell repeatedly executes nuclear division without cytokinesis, resulting in the development of a thallus (>100 mm) with characteristic morphology called coenocyte. Bryopsis is a representative coenocytic alga that has exceptionally high regeneration ability: extruded cytoplasm aggregates rapidly in seawater, leading to the formation of protoplasts. However, the genetic basis of the unique cell biology of Bryopsis remains poorly understood. Here, we present a high-quality assembly and annotation of the nuclear genome of Bryopsis sp. (90.7 Mbp, 27 contigs, N50 = 6.7 Mbp, 14 034 protein-coding genes). Comparative genomic analyses indicate that the genes encoding BPL-1/Bryohealin, the aggregation-promoting lectin, are heavily duplicated in Bryopsis, whereas homologous genes are absent in other ulvophyceans, suggesting the basis of regeneration capability of Bryopsis. Bryopsis sp. possesses >30 kinesins but only a single myosin, which differs from other green algae that have multiple types of myosin genes. Consistent with this biased motor toolkit, we observed that the bidirectional motility of chloroplasts in the cytoplasm was dependent on microtubules but not actin in Bryopsis sp. Most genes required for cytokinesis in plants are present in Bryopsis, including those in the SNARE or kinesin superfamily. Nevertheless, a kinesin crucial for cytokinesis initiation in plants (NACK/Kinesin-7II) is hardly expressed in the coenocytic part of the thallus, possibly underlying the lack of cytokinesis in this portion. The present genome sequence lays the foundation for experimental biology in coenocytic macroalgae.

Proteome plasticity during Physcomitrium patens spore germination - from the desiccated phase to heterotrophic growth and reconstitution of photoautotrophy.

The establishment of moss spores is considered a milestone in plant evolution. They harbor protein networks underpinning desiccation tolerance and accumulation of storage compounds that can be found already in algae and that are also utilized in seeds and pollen. Furthermore, germinating spores must produce proteins that drive the transition through heterotrophic growth to the autotrophic plant. To get insight into the plasticity of this proteome, we investigated it at five timepoints of moss (Physcomitrium patens) spore germination and in protonemata and gametophores. The comparison to previously published Arabidopsis proteome data of seedling establishment showed that not only the proteomes of spores and seeds are functionally related, but also the proteomes of germinating spores and young seedlings. We observed similarities with regard to desiccation tolerance, lipid droplet proteome composition, control of dormancy, and ß-oxidation and the glyoxylate cycle. However, there were also striking differences. For example, spores lacked any obvious storage proteins. Furthermore, we did not detect homologs to the main triacylglycerol lipase in Arabidopsis seeds, SUGAR DEPENDENT1. Instead, we discovered a triacylglycerol lipase of the oil body lipase family and a lipoxygenase as being the overall most abundant proteins in spores. This finding indicates an alternative pathway for triacylglycerol degradation via oxylipin intermediates in the moss. The comparison of spores to Nicotiana tabacum pollen indicated similarities for example in regards to resistance to desiccation and hypoxia, but the overall developmental pattern did not align as in the case of seedling establishment and spore germination.

Conquering new grounds: plant organellar C-to-U RNA editing factors can be functional in the plant cytosol.

Plant mitochondrial and chloroplast transcripts are subject to numerous events of specific cytidine-to-uridine (C-to-U) RNA editing to correct genetic information. Key protein factors for this process are specific RNA-binding pentatricopeptide repeat (PPR) proteins, which are encoded in the nucleus and post-translationally imported into the two endosymbiotic organelles. Despite hundreds of C-to-U editing sites in the plant organelles, no comparable editing has been found for nucleo-cytosolic mRNAs raising the question why plant RNA editing is restricted to chloroplasts and mitochondria. Here, we addressed this issue in the model moss Physcomitrium patens, where all PPR-type RNA editing factors comprise specific RNA-binding and cytidine deamination functionalities in single proteins. To explore whether organelle-type RNA editing can principally also take place in the plant cytosol, we expressed PPR56, PPR65 and PPR78, three editing factors recently shown to also function in a bacterial setup, together with cytosolic co-transcribed native targets in Physcomitrium. While we obtained unsatisfying results upon their constitutive expression, we found strong cytosolic RNA editing under hormone-inducible expression. Moreover, RNA-Seq analyses revealed varying numbers of up to more than 900 off-targets in other cytosolic transcripts. We conclude that PPR-mediated C-to-U RNA editing is not per se incompatible with the plant cytosol but that its limited target specificity has restricted its occurrence to the much less complex transcriptomes of mitochondria and chloroplast in the course of evolution.

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.

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.

The bryophytes Physcomitrium patens and Marchantia polymorpha as model systems for studying evolutionary cell and developmental biology in plants.

Bryophytes are nonvascular spore-forming plants. Unlike in flowering plants, the gametophyte (haploid) generation of bryophytes dominates the sporophyte (diploid) generation. A comparison of bryophytes with flowering plants allows us to answer some fundamental questions raised in evolutionary cell and developmental biology. The moss Physcomitrium patens was the first bryophyte with a sequenced genome. Many cell and developmental studies have been conducted in this species using gene targeting by homologous recombination. The liverwort Marchantia polymorpha has recently emerged as an excellent model system with low genomic redundancy in most of its regulatory pathways. With the development of molecular genetic tools such as efficient genome editing, both P. patens and M. polymorpha have provided many valuable insights. Here, we review these advances with a special focus on polarity formation at the cell and tissue levels. We examine current knowledge regarding the cellular mechanisms of polarized cell elongation and cell division, including symmetric and asymmetric cell division. We also examine the role of polar auxin transport in mosses and liverworts. Finally, we discuss the future of evolutionary cell and developmental biological studies in plants.

COPII Sec23 proteins form isoform-specific endoplasmic reticulum exit sites with differential effects on polarized growth.

Coat Protein complex II (COPII), a coat protein complex that forms vesicles on the endoplasmic reticulum (ER), mediates trafficking to the Golgi. While metazoans have few genes encoding each COPII component, plants have expanded these gene families, leading to the hypothesis that plant COPII has functionally diversified. In the moss Physcomitrium (Physcomitrella) patens, the Sec23/24 gene families are each composed of seven genes. Silencing Sec23/24 revealed isoform-specific contributions to polarized growth, with the closely related Sec23D/E and Sec24C/D essential for protonemal development. Focusing on Sec23, we discovered that Sec23D/E mediate ER-to Golgi transport and are essential for tip growth, with Sec23D localizing to presumptive ER exit sites. In contrast, Sec23A, B, C, F, and G are dispensable and do not quantitatively affect ER-to-Golgi trafficking. However, Δsec23abcfg plants exhibited reduced secretion of plasma membrane cargo. Of the four highly expressed protonemal Sec23 genes, Sec23F/G are members of a divergent Sec23 clade specifically retained in land plants. Notably, Sec23G accumulates on ER-associated foci that are significantly larger, do not overlap with, and are independent of Sec23D. While Sec23D/E form ER exit sites and function as bona fide COPII components essential for tip-growing protonemata, Sec23G and the closely related Sec23F have likely functionally diversified, forming separate and independent ER exit sites and participating in Golgi-independent trafficking pathways.

Transcriptional control of gene expression by microRNAs.

MicroRNAs (miRNAs) control gene expression in animals and plants. Like another class of small RNAs, siRNAs, they affect gene expression posttranscriptionally. While siRNAs in addition act in transcriptional gene silencing, a role of miRNAs in transcriptional regulation has been less clear. We show here that in moss Physcomitrella patens mutants without a DICER-LIKE1b gene, maturation of miRNAs is normal but cleavage of target RNAs is abolished and levels of these transcripts are drastically reduced. These mutants accumulate miRNA:target-RNA duplexes and show hypermethylation of the genes encoding target RNAs, leading to gene silencing. This pathway occurs also in the wild-type upon hormone treatment. We propose that initiation of epigenetic silencing by DNA methylation depends on the ratio of the miRNA and its target RNA.

Cytokinin-CLAVATA cross-talk is an ancient mechanism regulating shoot meristem homeostasis in land plants.

SignificancePlants grow from their tips. The gametophore (shoot-like organ) tip of the moss Physcomitrium patens is a single cell that performs the same functions as those of multicellular flowering plants, producing the cells that make leaves and regenerating new stem cells to maintain the shoot tip. Several pathways, including CLAVATA and cytokinin hormonal signaling, regulate stem cell abundance in flowering plants and in mosses, although the mechanisms whereby these pathways regulate stem cell abundance and their conservation between these plant lineages is poorly understood. Using moss, we investigated how PpCLAVATA and cytokinin signaling interact. Overall, we found evidence that PpCLAVATA and cytokinin signaling interact similarly in moss and flowering plants, despite their distinct anatomies, life cycles, and evolutionary distance.

Regulatory ligand binding in plant chalcone isomerase-like (CHIL) proteins.

Chalcone isomerase-like (CHIL) protein is a noncatalytic protein that enhances flavonoid content in green plants by serving as a metabolite binder and a rectifier of chalcone synthase (CHS). Rectification of CHS catalysis occurs through direct protein-protein interactions between CHIL and CHS, which alter CHS kinetics and product profiles, favoring naringenin chalcone (NC) production. These discoveries raise questions about how CHIL proteins interact structurally with metabolites and how CHIL-ligand interactions affect interactions with CHS. Using differential scanning fluorimetry on a CHIL protein from Vitis vinifera (VvCHIL), we report that positive thermostability effects are induced by the binding of NC, and negative thermostability effects are induced by the binding of naringenin. NC further causes positive changes to CHIL-CHS binding, whereas naringenin causes negative changes to VvCHIL-CHS binding. These results suggest that CHILs may act as sensors for ligand-mediated pathway feedback by influencing CHS function. The protein X-ray crystal structure of VvCHIL compared with the protein X-ray crystal structure of a CHIL from Physcomitrella patens reveals key amino acid differences at a ligand-binding site of VvCHIL that can be substituted to nullify the destabilizing effect caused by naringenin. Together, these results support a role for CHIL proteins as metabolite sensors that modulate the committed step of the flavonoid pathway.

N6 -Methyladenosine modification of mRNA contributes to the transition from 2D to 3D growth in the moss Physcomitrium patens.

Plants colonized the land approximately 470 million years ago, coinciding with the development of apical cells that divide in three planes. The molecular mechanisms that underly the development of the 3D growth pattern are poorly understood, mainly because 3D growth in seed plants starts during embryo development. In contrast, the transition from 2D to 3D growth in the moss Physcomitrium patens has been widely studied, and it involves a large turnover of the transcriptome to allow the establishment of stage-specific transcripts that facilitate this developmental transition. N6 -Methyladenosine (m6 A) is the most abundant, dynamic and conserved internal nucleotide modification present on eukaryotic mRNA and serves as a layer of post-transcriptional regulation directly affecting several cellular processes and developmental pathways in many organisms. In Arabidopsis, m6 A has been reported to be essential for organ growth and determination, embryo development and responses to environmental signals. In this study, we identified the main genes of the m6 A methyltransferase complex (MTC), MTA, MTB and FIP37, in P. patens and demonstrate that their inactivation leads to the loss of m6 A in mRNA, a delay in the formation of gametophore buds and defects in spore development. Genome-wide analysis revealed several transcripts affected in the Ppmta background. We demonstrate that the PpAPB1-PpAPB4 transcripts, encoding central factors orchestrating the transition from 2D to 3D growth in P. patens, are modified by m6 A, whereas in the Ppmta mutant the lack of the m6 A marker is associated with a corresponding decrease in transcript accumulation. Overall, we suggest that m6 A is essential to enable the proper accumulation of these and other bud-specific transcripts directing the turnover of stage-specific transcriptomes, and thus promoting the transition from protonema to gametophore buds in 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.

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.