The Physcomitrella patens chromosome-scale assembly reveals moss genome structure and evolution

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.

A molecular atlas of plastid and mitochondrial proteins reveals organellar remodeling during plant evolutionary transitions from algae to angiosperms

Algae and plants carry 2 organelles of endosymbiotic origin that have been co-evolving in their host cells for more than a billion years. The biology of plastids and mitochondria can differ significantly across major lineages and organelle changes likely accompanied the adaptation to new ecological niches such as the terrestrial habitat. Based on organelle proteome data and the genomes of 168 phototrophic (Archaeplastida) versus a broad range of 518 non-phototrophic eukaryotes, we screened for changes in plastid and mitochondrial biology across 1 billion years of evolution. Taking into account 331,571 protein families (or orthogroups), we identify 31,625 protein families that are unique to primary plastid-bearing eukaryotes. The 1,906 and 825 protein families are predicted to operate in plastids and mitochondria, respectively. Tracing the evolutionary history of these protein families through evolutionary time uncovers the significant remodeling the organelles experienced from algae to land plants. The analyses of gained orthogroups identifies molecular changes of organelle biology that connect to the diversification of major lineages and facilitated major transitions from chlorophytes en route to the global greening and origin of angiosperms.

Why should we study plant sex chromosomes?

Understanding plant sex chromosomes involves studying interactions between developmental and physiological genetics, genome evolution, and evolutionary ecology. We focus on areas of overlap between these. Ideas about how species with separate sexes (dioecious species, in plant terminology) can evolve are even more relevant to plants than to most animal taxa because dioecy has evolved many times from ancestral functionally hermaphroditic populations, often recently. One aim of studying plant sex chromosomes is to discover how separate males and females evolved from ancestors with no such genetic sex-determining polymorphism, and the diversity in the genetic control of maleness vs femaleness. Different systems share some interesting features, and their differences help to understand why completely sex-linked regions may evolve. In some dioecious plants, the sex-determining genome regions are physically small. In others, regions without crossing over have evolved sometimes extensive regions with properties very similar to those of the familiar animal sex chromosomes. The differences also affect the evolutionary changes possible when the environment (or pollination environment, for angiosperms) changes, as dioecy is an ecologically risky strategy for sessile organisms. Dioecious plants have repeatedly reverted to cosexuality, and hermaphroditic strains of fruit crops such as papaya and grapes are desired by plant breeders. Sex-linked regions are predicted to become enriched in genes with sex differences in expression, especially when higher expression benefits one sex function but harms the other. Such trade-offs may be important for understanding other plant developmental and physiological processes and have direct applications in plant breeding.

Diverse INOSITOL PHOSPHORYLCERAMIDE SYNTHASE mutant alleles of Physcomitrium patens offer new insight into complex sphingolipid metabolism

Sphingolipids are widespread, abundant, and essential lipids in plants and in other eukaryotes. Glycosyl inositol phosphorylceramides (GIPCs) are the most abundant class of plant sphingolipids, and are enriched in the plasma membrane of plant cells. They have been difficult to study due to lethal or pleiotropic mutant phenotypes. To overcome this, we developed a CRISPR/Cas9-based method for generating multiple and varied knockdown and knockout populations of mutants in a given gene of interest in the model moss Physcomitrium patens. This system is uniquely convenient due to the predominantly haploid state of the Physcomitrium life cycle, and totipotency of Physcomitrium protoplasts used for transformation. We used this approach to target the INOSITOL PHOSPHORYLCERAMIDE SYNTHASE (IPCS) gene family, which catalyzes the first, committed step in the synthesis of GIPCs. We isolated knockout single mutants and knockdown higher-order mutants showing a spectrum of deficiencies in GIPC content. Remarkably, we also identified two mutant alleles accumulating inositol phosphorylceramides, the direct products of IPCS activity, and provide our best explanation for this unexpected phenotype. Our approach is broadly applicable for studying essential genes and gene families, and for obtaining unusual lesions within a gene of interest.

Genomes of multicellular algal sisters to land plants illuminate signaling network evolution

Zygnematophyceae are the algal sisters of land plants. Here we sequenced four genomes of filamentous Zygnematophyceae, including chromosome-scale assemblies for three strains of Zygnema circumcarinatum. We inferred traits in the ancestor of Zygnematophyceae and land plants that might have ushered in the conquest of land by plants: expanded genes for signaling cascades, environmental response, and multicellular growth. Zygnematophyceae and land plants share all the major enzymes for cell wall synthesis and remodifications, and gene gains shaped this toolkit. Co-expression network analyses uncover gene cohorts that unite environmental signaling with multicellular developmental programs. Our data shed light on a molecular chassis that balances environmental response and growth modulation across more than 600 million years of streptophyte evolution.

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.

Re-examining meristems through the lens of evo-devo

The concept of the meristem was introduced in 1858 to characterize multicellular, formative, and proliferative tissues that give rise to the entire plant body, based on observations of vascular plants. Although its original definition did not encompass bryophytes, this concept has been used and continuously refined over the past 165 years to describe the diverse apices of all land plants. Here, we re-examine this matter in light of recent evo-devo research and show that, despite displaying high anatomical diversity, land plant meristems are unified by shared genetic control. We also propose a modular view of meristem function and highlight multiple evolutionary mechanisms that are likely to have contributed to the assembly and diversification of the varied meristems during the course of plant evolution.

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.

Drying without dying: A genome database for desiccation-tolerant plants and evolution of desiccation tolerance

Desiccation is typically fatal, but a small number of land plants have evolved vegetative desiccation tolerance (VDT), allowing them to dry without dying through a process called anhydrobiosis. Advances in sequencing technologies have enabled the investigation of genomes for desiccation-tolerant plants over the past decade. However, a dedicated and integrated database for these valuable genomic resources has been lacking. Our prolonged interest in VDT plant genomes motivated us to create the "Drying without Dying" database, which contains a total of 16 VDT-related plant genomes (including 10 mosses) and incorporates 10 genomes that are closely related to VDT plants. The database features bioinformatic tools, such as blast and homologous cluster search, sequence retrieval, Gene Ontology term and metabolic pathway enrichment statistics, expression profiling, co-expression network extraction, and JBrowser exploration for each genome. To demonstrate its utility, we conducted tailored PFAM family statistical analyses, and we discovered that the drought-responsive ABA transporter AWPM-19 family is significantly tandemly duplicated in all bryophytes but rarely so in tracheophytes. Transcriptomic investigations also revealed that response patterns following desiccation diverged between bryophytes and angiosperms. Combined, the analyses provided genomic and transcriptomic evidence supporting a possible divergence and lineage-specific evolution of VDT in plants. The database can be accessed at http://desiccation.novogene.com. We expect this initial release of the "Drying without Dying" plant genome database will facilitate future discovery of VDT genetic resources.

Regulation of developmental gatekeeping and cell fate transition by the calpain protease DEK1 in Physcomitrium patens

Calpains are cysteine proteases that control cell fate transitions whose loss of function causes severe, pleiotropic phenotypes in eukaryotes. Although mainly considered as modulatory proteases, human calpain targets are directed to the N-end rule degradation pathway. Several such targets are transcription factors, hinting at a gene-regulatory role. Here, we analyze the gene-regulatory networks of the moss Physcomitrium patens and characterize the regulons that are misregulated in mutants of the calpain DEFECTIVE KERNEL1 (DEK1). Predicted cleavage patterns of the regulatory hierarchies in five DEK1-controlled subnetworks are consistent with a pleiotropic and regulatory role during cell fate transitions targeting multiple functions. Network structure suggests DEK1-gated sequential transitions between cell fates in 2D-to-3D development. Our method combines comprehensive phenotyping, transcriptomics and data science to dissect phenotypic traits, and our model explains the protease function as a switch gatekeeping cell fate transitions potentially also beyond plant development.

Sequence similarity networks bear out hierarchical relationships of green cytochrome P450

Land plants have diversified enzyme families. One of the most prominent is the cytochrome P450 (CYP or CYP450) family. With over 443,000 CYP proteins sequenced across the tree of life, CYPs are ubiquitous in archaea, bacteria, and eukaryotes. Here, we focused on land plants and algae to study the role of CYP diversification. CYPs, acting as monooxygenases, catalyze hydroxylation reactions crucial for specialized plant metabolic pathways, including detoxification and phytohormone production; the CYPome consists of one enormous superfamily that is divided into clans and families. Their evolutionary history speaks of high substrate promiscuity; radiation and functional diversification have yielded numerous CYP families. To understand the evolutionary relationships within the CYPs, we employed sequence similarity network analyses. We recovered distinct clusters representing different CYP families, reflecting their diversified sequences that we link to the prediction of functionalities. Hierarchical clustering and phylogenetic analysis further elucidated relationships between CYP clans, uncovering their shared deep evolutionary history. We explored the distribution and diversification of CYP subfamilies across plant and algal lineages, uncovering novel candidates and providing insights into the evolution of these enzyme families. This identified unexpected relationships between CYP families, such as the link between CYP82 and CYP74, shedding light on their roles in plant defense signaling pathways. Our approach provides a methodology that brings insights into the emergence of new functions within the CYP450 family, contributing to the evolutionary history of plants and algae. These insights can be further validated and implemented via experimental setups under various external conditions.

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.

Diversified molecular adaptations of inorganic nitrogen assimilation and signaling machineries in plants

Plants evolved sophisticated machineries to monitor levels of external nitrogen supply, respond to nitrogen demand from different tissues and integrate this information for coordinating its assimilation. Although roles of inorganic nitrogen in orchestrating developments have been studied in model plants and crops, systematic understanding of the origin and evolution of its assimilation and signaling machineries remains largely unknown. We expanded taxon samplings of algae and early-diverging land plants, covering all main lineages of Archaeplastida, and reconstructed the evolutionary history of core components involved in inorganic nitrogen assimilation and signaling. Most components associated with inorganic nitrogen assimilation were derived from the ancestral Archaeplastida. Improvements of assimilation machineries by gene duplications and horizontal gene transfers were evident during plant terrestrialization. Clusterization of genes encoding nitrate assimilation proteins might be an adaptive strategy for algae to cope with changeable nitrate availability in different habitats. Green plants evolved complex nitrate signaling machinery that was stepwise improved by domains shuffling and regulation co-option. Our study highlights innovations in inorganic nitrogen assimilation and signaling machineries, ranging from molecular modifications of proteins to genomic rearrangements, which shaped developmental and metabolic adaptations of plants to changeable nutrient availability in environments.

Syntrichia ruralis: emerging model moss genome reveals a conserved and previously unknown regulator of desiccation in flowering plants

Water scarcity, resulting from climate change, poses a significant threat to ecosystems. Syntrichia ruralis, a dryland desiccation-tolerant moss, provides valuable insights into survival of water-limited conditions. We sequenced the genome of S. ruralis, conducted transcriptomic analyses, and performed comparative genomic and transcriptomic analyses with existing genomes and transcriptomes, including with the close relative S. caninervis. We took a genetic approach to characterize the role of an S. ruralis transcription factor, identified in transcriptomic analyses, in Arabidopsis thaliana. The genome was assembled into 12 chromosomes encompassing 21 169 protein-coding genes. Comparative analysis revealed copy number and transcript abundance differences in known desiccation-associated gene families, and highlighted genome-level variation among species that may reflect adaptation to different habitats. A significant number of abscisic acid (ABA)-responsive genes were found to be negatively regulated by a MYB transcription factor (MYB55) that was upstream of the S. ruralis ortholog of ABA-insensitive 3 (ABI3). We determined that this conserved MYB transcription factor, uncharacterized in Arabidopsis, acts as a negative regulator of an ABA-dependent stress response in Arabidopsis. The new genomic resources from this emerging model moss offer novel insights into how plants regulate their responses to water deprivation.

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.

Analysis of Marchantia polymorpha-microorganism interactions: basis for understanding plant-microbe and plant-pathogen interactions

Marchantia polymorpha is a bryophyte gaining significance as a model plant in evolutionary studies in recent years. This is attributed to its small-sequenced genome, standardized transformation methodology, global distribution, and easy and rapid in vitro culturing. As an evolutionary model, M. polymorpha contributes to our understanding of the evolution of plant defensive responses and the associated hormonal signaling pathways. Through its interaction with microorganisms, M. polymorpha serves as a valuable source of knowledge, yielding insights into new microbial species and bioactive compounds. Bibliographic analysis involved collecting, reading, and categorizing documents obtained from the Scopus and Web of Science databases using different search terms. The review was based on 30 articles published between 1995 and 2023, with Japanese and Spanish authors emerging as the most prolific contributors in this field. These articles have been grouped into four main themes: antimicrobial metabolites produced by M. polymorpha; identification and characterization of epiphytic, endophytic, and pathogenic microorganisms; molecular studies of the direct interaction between M. polymorpha and microorganisms; and plant transformation using bacterial vectors. This review highlights the key findings from these articles and identifies potential future research directions.

Phylogenomic insights into the first multicellular streptophyte

Streptophytes are best known as the clade containing the teeming diversity of embryophytes (land plants).1,2,3,4 Next to embryophytes are however a range of freshwater and terrestrial algae that bear important information on the emergence of key traits of land plants. Among these, the Klebsormidiophyceae stand out. Thriving in diverse environments-from mundane (ubiquitous occurrence on tree barks and rocks) to extreme (from the Atacama Desert to the Antarctic)-Klebsormidiophyceae can exhibit filamentous body plans and display remarkable resilience as colonizers of terrestrial habitats.5,6 Currently, the lack of a robust phylogenetic framework for the Klebsormidiophyceae hampers our understanding of the evolutionary history of these key traits. Here, we conducted a phylogenomic analysis utilizing advanced models that can counteract systematic biases. We sequenced 24 new transcriptomes of Klebsormidiophyceae and combined them with 14 previously published genomic and transcriptomic datasets. Using an analysis built on 845 loci and sophisticated mixture models, we establish a phylogenomic framework, dividing the six distinct genera of Klebsormidiophyceae in a novel three-order system, with a deep divergence more than 830 million years ago. Our reconstructions of ancestral states suggest (1) an evolutionary history of multiple transitions between terrestrial-aquatic habitats, with stem Klebsormidiales having conquered land earlier than embryophytes, and (2) that the body plan of the last common ancestor of Klebsormidiophyceae was multicellular, with a high probability that it was filamentous whereas the sarcinoids and unicells in Klebsormidiophyceae are likely derived states. We provide evidence that the first multicellular streptophytes likely lived about a billion years ago.

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.

Shared infection strategy of a fungal pathogen across diverse lineages of land plants, the Fusarium example

Plants engage with a wide variety of microorganisms either in parasitic or mutualistic relationships, which have helped them to adapt to terrestrial ecosystems. Microbial interactions have driven plant evolution and led to the emergence of complex interaction outcomes via suppression of host defenses by evolving pathogens. The evolution of plant-microbe interactions is shaped by conserved host and pathogen gene modules and fast-paced lineage-specific adaptability which determines the interaction outcome. Recent findings from different microbes ranging from bacteria, oomycetes, and fungi suggest recurrent concepts in establishing interactions with evolutionarily distant plant hosts, but also clade-specific adaptation that ultimately contributes to pathogenicity. Here, we revisit some of the latest features that illustrate shared colonization strategies of the fungal pathogen Fusarium oxysporum on distant plant lineages and lineage-specific adaptability of mini-chromosomal units encoding effectors, for shaping host-specific pathogenicity in angiosperms.

A designer synthetic chromosome fragment functions in moss

Rapid advances in DNA synthesis techniques have enabled the assembly and engineering of viral and microbial genomes, presenting new opportunities for synthetic genomics in multicellular eukaryotic organisms. These organisms, characterized by larger genomes, abundant transposons and extensive epigenetic regulation, pose unique challenges. Here we report the in vivo assembly of chromosomal fragments in the moss Physcomitrium patens, producing phenotypically virtually wild-type lines in which one-third of the coding region of a chromosomal arm is replaced by redesigned, chemically synthesized fragments. By eliminating 55.8% of a 155 kb endogenous chromosomal region, we substantially simplified the genome without discernible phenotypic effects, implying that many transposable elements may minimally impact growth. We also introduced other sequence modifications, such as PCRTag incorporation, gene locus swapping and stop codon substitution. Despite these substantial changes, the complex epigenetic landscape was normally established, albeit with some three-dimensional conformation alterations. The synthesis of a partial multicellular eukaryotic chromosome arm lays the foundation for the synthetic moss genome project (SynMoss) and paves the way for genome synthesis in multicellular organisms.

Near telomere-to-telomere genome of the model plant Physcomitrium patens

The model plant Physcomitrium patens has played a pivotal role in enhancing our comprehension of plant evolution and development. However, the current genome harbours numerous regions that remain unfinished and erroneous. To address these issues, we generated an assembly using Oxford Nanopore reads and Hi-C mapping. The assembly incorporates telomeric and centromeric regions, thereby establishing it as a near telomere-to-telomere genome except a region in chromosome 1 that is not fully assembled due to its highly repetitive nature. This near telomere-to-telomere genome resolves the chromosome number at 26 and provides a gap-free genome assembly as well as updated gene models to aid future studies using this model organism.

Telomere Length Variation in Model Bryophytes

The ends of linear chromosomes of most eukaryotes consist of protein-bound DNA arrays called telomeres, which play essential roles in protecting genome integrity. Despite general evolutionary conservation in function, telomeric DNA is known to drastically vary in length and sequence between different eukaryotic lineages. Bryophytes are a group of early diverging land plants that include mosses, liverworts, and hornworts. This group of ancient land plants recently emerged as a new model for important discoveries in genomics and evolutionary biology, as well as for understanding plant adaptations to a terrestrial lifestyle. We measured telomere length in different ecotypes of model bryophyte species, including Physcomitrium patens, Marchantia polymorpha, Ceratodon purpureus, and in Sphagnum isolates. Our data indicate that all analyzed moss and liverwort genotypes have relatively short telomeres. Furthermore, all analyzed ecotypes and isolates of model mosses and liverworts display evidence of substantial natural variation in telomere length. Interestingly, telomere length also differs between male and female strains of the dioecious liverwort M. polymorpha and dioecious moss C. purpureus. Given that bryophytes are extraordinarily well adapted to different ecological niches from polar to tropical environments, our data will contribute to understanding the impact of natural telomere length variation on evolutionary adaptations in this ancient land plant lineage.

Extraordinary preservation of gene collinearity over three hundred million years revealed in homosporous lycophytes

Homosporous lycophytes (Lycopodiaceae) are a deeply diverged lineage in the plant tree of life, having split from heterosporous lycophytes (Selaginella and Isoetes) ~400 Mya. Compared to the heterosporous lineage, Lycopodiaceae has markedly larger genome sizes and remains the last major plant clade for which no chromosome-level assembly has been available. Here, we present chromosomal genome assemblies for two homosporous lycophyte species, the allotetraploid Huperzia asiatica and the diploid Diphasiastrum complanatum. Remarkably, despite that the two species diverged ~350 Mya, around 30% of the genes are still in syntenic blocks. Furthermore, both genomes had undergone independent whole genome duplications, and the resulting intragenomic syntenies have likewise been preserved relatively well. Such slow genome evolution over deep time is in stark contrast to heterosporous lycophytes and is correlated with a decelerated rate of nucleotide substitution. Together, the genomes of H. asiatica and D. complanatum not only fill a crucial gap in the plant genomic landscape but also highlight a potentially meaningful genomic contrast between homosporous and heterosporous species.

Multifactorial analysis of terminator performance on heterologous gene expression in Physcomitrella

KEY MESSAGE

Characterization of Physcomitrella 3'UTRs across different promoters yields endogenous single and double terminators for usage in molecular pharming. The production of recombinant proteins for health applications accounts for a large share of the biopharmaceutical market. While many drugs are produced in microbial and mammalian systems, plants gain more attention as expression hosts to produce eukaryotic proteins. In particular, the good manufacturing practice (GMP)-compliant moss Physcomitrella (Physcomitrium patens) has outstanding features, such as excellent genetic amenability, reproducible bioreactor cultivation, and humanized protein glycosylation patterns. In this study, we selected and characterized novel terminators for their effects on heterologous gene expression. The Physcomitrella genome contains 53,346 unique 3'UTRs (untranslated regions) of which 7964 transcripts contain at least one intron. Over 91% of 3'UTRs exhibit more than one polyadenylation site, indicating the prevalence of alternative polyadenylation in Physcomitrella. Out of all 3'UTRs, 14 terminator candidates were selected and characterized via transient Dual-Luciferase assays, yielding a collection of endogenous terminators performing equally high as established heterologous terminators CaMV35S, AtHSP90, and NOS. High performing candidates were selected for testing as double terminators which impact reporter levels, dependent on terminator identity and positioning. Testing of 3'UTRs among the different promoters NOS, CaMV35S, and PpActin5 showed an increase of more than 1000-fold between promoters PpActin5 and NOS, whereas terminators increased reporter levels by less than tenfold, demonstrating the stronger effect promoters play as compared to terminators. Among selected terminator attributes, the number of polyadenylation sites as well as polyadenylation signals were found to influence terminator performance the most. Our results improve the biotechnology platform Physcomitrella and further our understanding of how terminators influence gene expression in plants in general.

Immunobiodiversity: Conserved and specific immunity across land plants and beyond

Angiosperms represent most plants that humans cultivate, grow, and eat. However, angiosperms are only one of five major land plant lineages. As a whole lineage, plants also include algal groups. All these clades represent a tremendous genetic diversity that can be investigated to reveal the evolutionary history of any given mechanism. In this review, we describe the current model of the plant immune system, discuss its evolution based on the recent literature, and propose future directions for the field. In angiosperms, plant-microbe interactions have been intensively studied, revealing essential cell surface and intracellular immune receptors, as well as metabolic and hormonal defense pathways. Exploring diversity at the genomic and functional levels demonstrates the conservation of these pathways across land plants, some of which are beyond plants. On basis of the conserved mechanisms, lineage-specific variations have occurred, leading to diversified reservoirs of immune mechanisms. In rare cases, this diversity has been harnessed and successfully transferred to other species by integration of wild immune receptors or engineering of novel forms of receptors for improved resistance to pathogens. We propose that exploring further the diversity of immune mechanisms in the whole plant lineage will reveal completely novel sources of resistance to be deployed in crops.

The advent of plant cells in bioreactors

Ever since agriculture started, plants have been bred to obtain better yields, better fruits, or sustainable products under uncertain biotic and abiotic conditions. However, a new way to obtain products from plant cells emerged with the development of recombinant DNA technologies. This led to the possibility of producing exogenous molecules in plants. Furthermore, plant chemodiversity has been the main source of pharmacological molecules, opening a field of plant biotechnology directed to produce high quality plant metabolites. The need for different products by the pharma, cosmetics agriculture and food industry has pushed again to develop new procedures. These include cell production in bioreactors. While plant tissue and cell culture are an established technology, beginning over a hundred years ago, plant cell cultures have shown little impact in biotechnology projects, compared to bacterial, yeasts or animal cells. In this review we address the different types of bioreactors that are currently used for plant cell production and their usage for quality biomolecule production. We make an overview of Nicotiana tabacum, Nicotiana benthamiana, Oryza sativa, Daucus carota, Vitis vinifera and Physcomitrium patens as well-established models for plant cell culture, and some species used to obtain important metabolites, with an insight into the type of bioreactor and production protocols.

Control of stem cell behavior by CLE-JINGASA signaling in the shoot apical meristem in Marchantia polymorpha

Land plants undergo indeterminate growth by the activity of meristems in both gametophyte (haploid) and sporophyte (diploid) generations. In the sporophyte of the flowering plant Arabidopsis thaliana, the apical meristems are located at the shoot and root tips in which a number of regulatory gene homologs are shared for their development, implying deep evolutionary origins. However, little is known about their functional conservation with gametophytic meristems in distantly related land plants such as bryophytes, even though genomic studies have revealed that the subfamily-level diversity of regulatory genes is mostly conserved throughout land plants. Here, we show that a NAM/ATAF/CUC (NAC) domain transcription factor, JINGASA (MpJIN), acts downstream of CLAVATA3 (CLV3)/ESR-related (CLE) peptide signaling and controls stem cell behavior in the gametophytic shoot apical meristem of the liverwort Marchantia polymorpha. In the meristem, strong MpJIN expression was associated with the periclinal cell division at the periphery of the stem cell zone (SCZ), whereas faint MpJIN expression was found at the center of the SCZ. Time course observation indicates that the MpJIN-negative cells are lost from the SCZ and respecified de novo at two separate positions during the dichotomous branching event. Consistently, the induction of MpJIN results in ectopic periclinal cell division in the SCZ and meristem termination. Based on the comparative expression data, we speculate that the function of JIN/FEZ subfamily genes was shared among the shoot apical meristems in the gametophyte and sporophyte generations in early land plants but was lost in certain lineages, including the flowering plant A. thaliana.

Divergent outcomes of genetic conflict on the UV sex chromosomes of Marchantia polymorpha and Ceratodon purpureus

In species with separate sexes, the genome must produce two distinct developmental programs. Sexually dimorphic development may be controlled by either sex-limited loci or biased expression of loci transmitted through both sexes. Variation in the gene content of sex-limited chromosomes demonstrates that eukaryotic species differ markedly in the roles of these two mechanisms in governing sexual dimorphism. The bryophyte model systems Marchantia polymorpha and Ceratodon purpureus provide a particularly striking contrast. Although both species possess a haploid UV sex chromosome system, in which females carry a U chromosome and males carry a V, M. polymorpha relies on biased autosomal expression, while in C. purpureus, sex-linked genes drive dimorphism. Framing these genetic architectures as divergent outcomes of genetic conflict highlights comparative genomic analyses to better understand the evolution of sexual dimorphism.

The ancestral chromatin landscape of land plants

Recent studies have shown that correlations between chromatin modifications and transcription vary among eukaryotes. This is the case for marked differences between the chromatin of the moss Physcomitrium patens and the liverwort Marchantia polymorpha. Mosses and liverworts diverged from hornworts, altogether forming the lineage of bryophytes that shared a common ancestor with land plants. We aimed to describe chromatin in hornworts to establish synapomorphies across bryophytes and approach a definition of the ancestral chromatin organization of land plants. We used genomic methods to define the 3D organization of chromatin and map the chromatin landscape of the model hornwort Anthoceros agrestis. We report that nearly half of the hornwort transposons were associated with facultative heterochromatin and euchromatin and formed the center of topologically associated domains delimited by protein coding genes. Transposons were scattered across autosomes, which contrasted with the dense compartments of constitutive heterochromatin surrounding the centromeres in flowering plants. Most of the features observed in hornworts are also present in liverworts or in mosses but are distinct from flowering plants. Hence, the ancestral genome of bryophytes was likely a patchwork of units of euchromatin interspersed within facultative and constitutive heterochromatin. We propose this genome organization was ancestral to land plants.

Mosses

Often overlooked, these small but otherwise brilliant plants began covering Earth's land masses more than 450 million years ago. They saw the dinosaurs come and go, and they saw us humans coming. Mosses, liverworts and hornworts comprise the bryophytes, the second largest monophyletic clade of land plants (embryophytes), after the vascular plants (tracheophytes). Like all embryophytes, mosses exhibit a haplodiplontic life cycle. This alternation of generations (originally termed Generationswechsel in German) between the haploid gametophyte and the diploid sporophyte implies that every plant genome encodes two distinct ontogenies, in contrast to animal genomes. Contrary to tracheophytes, the haploid gametophyte is the dominant generation in mosses. Haploidy of the major tissues facilitates gene-function annotation via reverse genetics. Nevertheless, the diploid sporophyte of mosses, the spore capsule, is a visible structure unlike the gametophyte of flowering plants, which is largely reduced and embedded in the sporophyte. Visibility of both generations on one plant facilitates the analysis of the alternation of generations, and in a broader sense evo-devo studies. Whereas the conservation of moss morphology over hundreds of millions of years suggests stasis, molecular data reveal fast evolving moss genomes, leaving an enigma for evolutionary biologists. Finally, the extraordinary resilience of mosses may provide lessons for current man-made climate change. In this Primer, we will highlight some of the peculiarities of mosses from historical observations to current genomic data, with an emphasis on their development, reproduction, evolution, biotic interactions, and potential for biotechnology. Mosses from three genera - the living fossil Takakia, the ecosystems engineer Sphagnum, and the model moss Physcomitrella - exemplify the scientific insights and the applications mosses have to offer.

Comparative mutant analyses reveal a novel mechanism of ARF regulation in land plants

A major challenge in plant biology is to understand how the plant hormone auxin regulates diverse transcriptional responses throughout development, in different environments, and in different species. The answer may lie in the specific complement of auxin signaling components in each cell. The balance between activators (class-A AUXIN RESPONSE FACTORS) and repressors (class-B ARFs) is particularly important. It is unclear how this balance is achieved. Through comparative analysis of novel, dominant mutants in maize and the moss Physcomitrium patens , we have discovered a ∼500-million-year-old mechanism of class-B ARF protein level regulation, important in determining cell fate decisions across land plants. Thus, our results add a key piece to the puzzle of how auxin regulates plant development.

Expression divergence of expansin genes drive the heteroblasty in Ceratopteris chingii

BACKGROUND

Sterile-fertile heteroblasty is a common phenomenon observed in ferns, where the leaf shape of a fern sporophyll, responsible for sporangium production, differs from that of a regular trophophyll. However, due to the large size and complexity of most fern genomes, the molecular mechanisms that regulate the formation of these functionally different heteroblasty have remained elusive. To shed light on these mechanisms, we generated a full-length transcriptome of Ceratopteris chingii with PacBio Iso-Seq from five tissue samples. By integrating Illumina-based sequencing short reads, we identified the genes exhibiting the most significant differential expression between sporophylls and trophophylls.

RESULTS

The long reads were assembled, resulting in a total of 24,024 gene models. The differential expressed genes between heteroblasty primarily involved reproduction and cell wall composition, with a particular focus on expansin genes. Reconstructing the phylogeny of expansin genes across 19 plant species, ranging from green algae to seed plants, we identified four ortholog groups for expansins. The observed high expression of expansin genes in the young sporophylls of C. chingii emphasizes their role in the development of heteroblastic leaves. Through gene coexpression analysis, we identified highly divergent expressions of expansin genes both within and between species.

CONCLUSIONS

The specific regulatory interactions and accompanying expression patterns of expansin genes are associated with variations in leaf shapes between sporophylls and trophophylls.

Chromosome-level genome assembly of Niphotrichum japonicum provides new insights into heat stress responses in mosses

With a diversity of approximately 22,000 species, bryophytes (hornworts, liverworts, and mosses) represent a major and diverse lineage of land plants. Bryophytes can thrive in many extreme environments as they can endure the stresses of drought, heat, and cold. The moss Niphotrichum japonicum (Grimmiaceae, Grimmiales) can subsist for extended periods under heat and drought conditions, providing a good candidate for studying the genetic basis underlying such high resilience. Here, we de novo assembled the genome of N. japonicum using Nanopore long reads combined with Hi-C scaffolding technology to anchor the 191.61 Mb assembly into 14 pseudochromosomes. The genome structure of N. japonicum's autosomes is mostly conserved and highly syntenic, in contrast to the sparse and disordered genes present in its sex chromosome. Comparative genomic analysis revealed the presence of 10,019 genes exclusively in N. japonicum. These genes may contribute to the species-specific resilience, as demonstrated by the gene ontology (GO) enrichment. Transcriptome analysis showed that 37.44% (including 3,107 unique genes) of the total annotated genes (26,898) exhibited differential expression as a result of heat-induced stress, and the mechanisms that respond to heat stress are generally conserved across plants. These include the upregulation of HSPs, LEAs, and reactive oxygen species (ROS) scavenging genes, and the downregulation of PPR genes. N. japonicum also appears to have distinctive thermal mechanisms, including species-specific expansion and upregulation of the Self-incomp_S1 gene family, functional divergence of duplicated genes, structural clusters of upregulated genes, and expression piggybacking of hub genes. Overall, our study highlights both shared and species-specific heat tolerance strategies in N. japonicum, providing valuable insights into the heat tolerance mechanism and the evolution of resilient plants.

Towards Plant Synthetic Genomics

Rapid advances in DNA synthesis techniques have allowed the assembly and engineering of viral and microbial genomes. Multicellular eukaryotic organisms, with their larger genomes, abundant transposons, and prevalent epigenetic regulation, present a new frontier to synthetic genomics. Plant synthetic genomics have long been proposed, and exciting progress has been made using the top-down approach. In this perspective, we propose applying bottom-up genome synthesis in multicellular plants, starting from the model moss Physcomitrium patens, in which homologous recombination, DNA delivery, and regeneration are possible, although further optimizations are necessary. We then discuss technical barriers, including genome assembly and plant transformation, associated with synthetic genomics in seed plants.

Identification of SRS transcription factor family in Solanum lycopersicum, and functional characterization of their responses to hormones and abiotic stresses

The SHI RELATED SEQUENCE (SRS) family plays a vital role in the development of multiple plant organs such as floral meristem determinacy, organ morphogenesis, and signal transduction. Nevertheless, there is little understanding of the biological significance of tomato SRS family at this point. Our research identified eight SlSRS family members and classified them into three subfamilies based on phylogenetics, conserved motifs, and characteristic domain analysis. The intraspecies and interspecies collinearity analysis revealed clues of SRS family evolution. Many cis-elements related to hormones, stresses, and plant development can be found in the promoter region of SlSRS genes. All of eight SlSRS proteins were located in the nucleus and possessed transcriptional activity, half of which were transcriptional activators, and the other half were transcriptional repressors. Except for SlSRS1, which showed high transcript accumulation in vegetative organs, most SlSRS genes expressed ubiquitously in all flower organs. In addition, all SlSRS genes could significantly respond to at least four different plant hormones. Further, expression of SlSRS genes were regulated by various abiotic stress conditions. In summary, we systematically analyzed and characterized the SlSRS family, reviewed the expression patterns and preliminarily investigated the protein function, and provided essential information for further functional research of the tomato SRS genes in the determination of reproductive floral organs and the development of plants, and possibly other plants.

Genome and transcriptome of Selaginella kraussiana reveal evolution of root apical meristems in vascular plants

The evolution of roots allowed vascular plants to adapt to land environments. Fossil evidence indicates that roots evolved independently in euphyllophytes (ferns and seed plants) and lycophytes, the two lineages of extant vascular plants. Based on a high-quality genome assembly, mRNA sequencing (mRNA-seq) data, and single-cell RNA-seq data for the lycophyte Selaginella kraussiana, we show that the two root origin events in lycophytes and euphyllophytes adopted partially similar molecular modules in the regulation of root apical meristem (RAM) development. In S. kraussiana, the RAM initiates from the rhizophore primordium guided by auxin and duplicates itself by dichotomous branching. The auxin signaling pathway directly upregulates euAINTEGUMENTAb (SkeuANTb), and then SkeuANTb directly promotes the expression of SkeuANTa and the WUSCHEL-RELATED HOMEOBOX13b (SkWOX13b) for RAM maintenance, partially similar to the molecular pathway involving the euANT-branch PLETHORA (AtPLT) genes and AtWOX5 in root initiation in the seed plant Arabidopsis thaliana. Other molecular modules, e.g., SHORT-ROOT and SCARECROW, also have partially similar expression patterns in the RAMs of S. kraussiana and A. thaliana. Overall, our study not only provides genome and transcriptome tools of S. kraussiana but also indicates the employment of some common molecular modules in RAMs during root origins in lycophytes and euphyllophytes.

Transposable element evolution in plant genome ecosystems

The relationship of transposable elements (TEs) with their host genomes has usually been seen as an arms race between TEs and their host genomes. Consequently, TEs are supposed to amplify by bursts of transposition, when the TE escapes host surveillance, followed by long periods of TE quiescence and efficient host control. Recent data obtained from an increasing number of assembled plant genomes and resequencing population datasets show that TE dynamics is more complex and varies among TE families and their host genomes. This variation ranges from large genomes that accommodate large TE populations to genomes that are very active in TE elimination, and from inconspicuous elements with very low activity to elements with high transposition and elimination rates. The dynamics of each TE family results from a long history of interaction with the host in a genome populated by many other TE families, very much like an evolving ecosystem.

Genome-wide analysis of the GDSL esterase/lipase family genes in Physcomitrium patens and the involvement of GELP31 in spore germination

In plants, the ability to produce hydrophobic substances that would provide protection from dehydration was required for the transition to land. This genome-wide investigation outlines the evolution of GDSL-type esterase/lipase (GELP) proteins in the moss Physcomitrium patens and suggests possible functions of some genes. GELP proteins play roles in the formation of hydrophobic polymers such as cutin and suberin that protect against dehydration and pathogen attack. GELP proteins are also implicated in processes such as pollen development and seed metabolism and germination. The P. patens GELP gene family comprises 48 genes and 14 pseudogenes. Phylogenetic analysis of all P. patens GELP sequences along with vascular plant GELP proteins with reported functions revealed that the P. patens genes clustered within previously identified A, B and C clades. A duplication model predicting the expansion of the GELP gene family within the P. patens lineage was constructed. Expression analysis combined with phylogenetic analysis suggested candidate genes for functions such as defence against pathogens, cutin metabolism, spore development and spore germination. The presence of relatively fewer GELP genes in P. patens may reduce the occurrence of functional redundancy that complicates the characterization of vascular plant GELP genes. Knockout lines of GELP31, which is highly expressed in sporophytes, were constructed. Gelp31 spores contained amorphous oil bodies and germinated late, suggesting (a) role(s) of GELP31 in lipid metabolism in spore development or germination. Future knockout studies of other candidate GELP genes will further elucidate the relationship between expansion of the family and the ability to withstand the harsh land environment.

Insight into gene expression associated with DNA methylation and small RNA in the rhizome-root system of Moso bamboo

Moso bamboo (Phyllostachys edulis), typically a monopodial scattering bamboo, is famous for its rapid growth. The rhizome-root system of Moso bamboo plays a crucial role in its clonal growth and spatial distribution. However, few studies have focused on rhizome-root systems. Here we collected LBs, RTs, and RGFNSs, the most important parts of the rhizome-root system, to study the molecular basis of the rapid growth of Moso bamboo due to epigenetic changes, such as DNA modifications and small RNAs. The angle of the shoot apical meristem of LB gradually decreased with increasing distance from the mother plant, and the methylation levels of LB were much higher than those of RT and RGFNS. 24 nt small RNAs and mCHH exhibited similar distribution patterns in transposable elements, suggesting a potential association between these components. The miRNA abundance of LB gradually increased with increasing distance from the mother plant, and a negative correlation was observed between gene expression levels and mCG and mCHG levels in the gene body. This study paves the way for further exploring the effects of epigenetic factors on the physiology of Moso bamboo.

Environmental gradients reveal stress hubs pre-dating plant terrestrialization

Plant terrestrialization brought forth the land plants (embryophytes). Embryophytes account for most of the biomass on land and evolved from streptophyte algae in a singular event. Recent advances have unravelled the first full genomes of the closest algal relatives of land plants; among the first such species was Mesotaenium endlicherianum. Here we used fine-combed RNA sequencing in tandem with a photophysiological assessment on Mesotaenium exposed to a continuous range of temperature and light cues. Our data establish a grid of 42 different conditions, resulting in 128 transcriptomes and ~1.5 Tbp (~9.9 billion reads) of data to study the combinatory effects of stress response using clustering along gradients. Mesotaenium shares with land plants major hubs in genetic networks underpinning stress response and acclimation. Our data suggest that lipid droplet formation and plastid and cell wall-derived signals have denominated molecular programmes since more than 600 million years of streptophyte evolution-before plants made their first steps on land.

DIENELACTONE HYDROLASE LIKE PROTEIN1 negatively regulates the KAI2-ligand pathway in Marchantia polymorpha

Karrikins are smoke-derived butenolides that induce seed germination and photomorphogenesis in a wide range of plants.1,2,3 KARRIKIN INSENSITIVE2 (KAI2), a paralog of a strigolactone receptor, perceives karrikins or their metabolized products in Arabidopsis thaliana.4,5,6,7 Furthermore, KAI2 is thought to perceive an unidentified plant hormone, called KAI2 ligand (KL).8,9 KL signal is transduced via the interaction between KAI2, MORE AXILLARY GROWTH2 (MAX2), and SUPPRESSOR of MORE AXILLARY GROWTH2 1 LIKE family proteins (SMXLs), followed by the degradation of SMXLs.4,7,10,11,12,13,14 This signaling pathway is conserved both in A. thaliana and the bryophyte Marchantia polymorpha.14 Although the KL signaling pathway is well characterized, the KL metabolism pathways remain poorly understood. Here, we show that DIENELACTONE HYDROLASE LIKE PROTEIN1 (DLP1) is a negative regulator of the KL pathway in M. polymorpha. The KL signal induces DLP1 expression. DLP1 overexpression lines phenocopied the Mpkai2a and Mpmax2 mutants, while dlp1 mutants phenocopied the Mpsmxl mutants. Mutations in the KL signaling genes largely suppressed these phenotypes, indicating that DLP1 acts upstream of the KL signaling pathway, although DLP1 also has KL pathway-independent functions. DLP1 exhibited enzymatic activity toward a potential substrate, suggesting the possibility that DLP1 works through KL inactivation. Investigation of DLP1 homologs in A. thaliana revealed that they do not play a major role in the KL pathway, suggesting different mechanisms for the KL signal regulation. Our findings provide new insights into the regulation of the KL signal in M. polymorpha and the evolution of the KL pathway in land plants.

Long-Distance Electrical and Calcium Signals Evoked by Hydrogen Peroxide in Physcomitrella

Electrical and calcium signals in plants are some of the basic carriers of information that are transmitted over a long distance. Together with reactive oxygen species (ROS) waves, electrical and calcium signals can participate in cell-to-cell signaling, conveying information about different stimuli, e.g. abiotic stress, pathogen infection or mechanical injury. There is no information on the ability of ROS to evoke systemic electrical or calcium signals in the model moss Physcomitrella nor on the relationships between these responses. Here, we show that the external application of hydrogen peroxide (H2O2) evokes electrical signals in the form of long-distance changes in the membrane potential, which transmit through the plant instantly after stimulation. The responses were calcium-dependent since their generation was inhibited by lanthanum, a calcium channel inhibitor (2 mM), and EDTA, a calcium chelator (0.5 mM). The electrical signals were partially dependent on glutamate receptor (GLR) ion channels since knocking-out the GLR genes only slightly reduced the amplitude of the responses. The basal part of the gametophyte, which is rich in protonema cells, was the most sensitive to H2O2. The measurements carried out on the protonema expressing fluorescent calcium biosensor GCaMP3 proved that calcium signals propagated slowly (>5 µm/s) and showed a decrement. We also demonstrate upregulation of a stress-related gene that appears in a distant section of the moss 8 min after the H2O2 treatment. The results help understand the importance of both types of signals in the transmission of information about the appearance of ROS in the plant cell apoplast.

Adaptive evolution of the enigmatic Takakia now facing climate change in Tibet

The most extreme environments are the most vulnerable to transformation under a rapidly changing climate. These ecosystems harbor some of the most specialized species, which will likely suffer the highest extinction rates. We document the steepest temperature increase (2010-2021) on record at altitudes of above 4,000 m, triggering a decline of the relictual and highly adapted moss Takakia lepidozioides. Its de-novo-sequenced genome with 27,467 protein-coding genes includes distinct adaptations to abiotic stresses and comprises the largest number of fast-evolving genes under positive selection. The uplift of the study site in the last 65 million years has resulted in life-threatening UV-B radiation and drastically reduced temperatures, and we detected several of the molecular adaptations of Takakia to these environmental changes. Surprisingly, specific morphological features likely occurred earlier than 165 mya in much warmer environments. Following nearly 400 million years of evolution and resilience, this species is now facing extinction.

Heterologous Production of Artemisinin in Physcomitrium patens by Direct in vivo Assembly of Multiple DNA Fragments

The sesquiterpene lactone compound artemisinin is a natural medicinal product of commercial importance. This Artemisia annua-derived secondary metabolite is well known for its antimalarial activity and has been studied in several other biological assays. However, the major shortcoming in its production and commercialization is its low accumulation in the native plant. Moreover, the chemical synthesis of artemisinin is difficult and expensive due to its complex structure. Hence, an alternative and sustainable production system of artemisinin in a heterologous host is required. Previously, heterologous production of artemisinin was achieved by Agrobacterium-mediated transformation. However, this requires extensive bioengineering of modified Nicotiana plants. Recently, a technique involving direct in vivo assembly of multiple DNA fragments in the moss, P. patens, has been successfully established. We utilized this technique to engineer artemisinin biosynthetic pathway genes into the moss, and artemisinin was obtained without further modifications with high initial production. Here, we provide protocols for establishing moss culture accumulating artemisinin, including culture preparation, transformation method, and compound detection via HS-SPME, UPLC-MRM-MS, and LC-QTOF-MS. The bioengineering of moss opens up a more sustainable, cost effective, and scalable platform not only in artemisinin production but also other high-value specialized metabolites in the future.

The first homosporous lycophyte genome revealed the association between the recent dynamic accumulation of LTR-RTs and genome size variation

The contrasting genome size between homosporous and heterosporous plants is fascinating. Different from the heterosporous seed plants and mainly homosporous ferns, the lycophytes are either heterosporous (Isoetales and Selaginellales) or homosporous (Lycopodiales). Many lycophytes are the resource plants of Huperzine A (HupA) which is invaluable for treating Alzheimer's disease. For the seed-free vascular plants, several high-quality genomes of heterosporous Selaginella, homosporous ferns (maidenhair fern, monkey spider tree fern), and heterosporous ferns (Azolla) have been published and provided important insights into the origin and evolution of early land plants. However, the homosporous lycophyte genome has not been decoded. Here, we assembled the first homosporous lycophyte genome and conducted comparative genomic analyses by applying a reformed pipeline for filtering out non-plant sequences. The obtained genome size of Lycopodium clavatum is 2.30 Gb, distinguished in more than 85% repetitive elements of which 62% is long terminal repeat (LTR). This study disclosed a high birth rate and a low death rate of the LTR-RTs in homosporous lycophytes, but the opposite occurs in heterosporous lycophytes. we propose that the recent activity of LTR-RT is responsible for the immense genome size variation between homosporous and heterosporous lycophytes. By combing Ks analysis with a phylogenetic approach, we discovered two whole genome duplications (WGD). Morover, we identified all the five recognized key enzymes for the HupA biosynthetic pathway in the L. clavatum genome, but found this pathway incomplete in other major lineages of land plants. Overall, this study is of great importance for the medicinal utilization of lycophytes and the decoded genome data will be a key cornerstone to elucidate the evolution and biology of early vascular land plants.

Cis-regulatory elements and transcription factors related to auxin signaling in the streptophyte algae Klebsormidium nitens

The phytohormone auxin affects numerous processes in land plants. The central auxin signaling machinery, called the nuclear auxin pathway, is mediated by its pivotal receptor named TRANSPORT INHIBITOR RESPONSE 1/AUXIN SIGNALING F-BOX (TIR1/AFB). The nuclear auxin pathway is widely conserved in land plants, but auxin also accumulates in various algae. Although auxin affects the growth of several algae, the components that mediate auxin signaling have not been identified. We previously reported that exogenous auxin suppresses cell proliferation in the Klebsormidium nitens that is a member of streptophyte algae, a paraphyletic group sharing the common ancestor with land plants. Although K. nitens lacks TIR1/AFB, auxin affects the expression of numerous genes. Thus, elucidation of the mechanism of auxin-inducible gene expression in K. nitens would provide important insights into the evolution of auxin signaling. Here, we show that some motifs are enriched in the promoter sequences of auxin-inducible genes in K. nitens. We also found that the transcription factor KnRAV activates several auxin-inducible genes and directly binds the promoter of KnLBD1, a representative auxin-inducible gene. We propose that KnRAV has the potential to regulate auxin-responsive gene expression in K. nitens.

Cellulose synthase-like D movement in the plasma membrane requires enzymatic activity

Cellulose Synthase-Like D (CSLD) proteins, important for tip growth and cell division, are known to generate β-1,4-glucan. However, whether they are propelled in the membrane as the glucan chains they produce assemble into microfibrils is unknown. To address this, we endogenously tagged all eight CSLDs in Physcomitrium patens and discovered that they all localize to the apex of tip-growing cells and to the cell plate during cytokinesis. Actin is required to target CSLD to cell tips concomitant with cell expansion, but not to cell plates, which depend on actin and CSLD for structural support. Like Cellulose Synthase (CESA), CSLD requires catalytic activity to move in the plasma membrane. We discovered that CSLD moves significantly faster, with shorter duration and less linear trajectories than CESA. In contrast to CESA, CSLD movement was insensitive to the cellulose synthesis inhibitor isoxaben, suggesting that CSLD and CESA function within different complexes possibly producing structurally distinct cellulose microfibrils.

Clade-D auxin response factors regulate auxin signaling and development in the moss Physcomitrium patens

Auxin response factors (ARFs) are a family of transcription factors that are responsible for regulating gene expression in response to changes in auxin level. The analysis of ARF sequence and activity indicates that there are 2 major groups: activators and repressors. One clade of ARFs, clade-D, is sister to clade-A activating ARFs, but are unique in that they lack a DNA-binding domain. Clade-D ARFs are present in lycophytes and bryophytes but absent in other plant lineages. The transcriptional activity of clade-D ARFs, as well as how they regulate gene expression, is not well understood. Here, we report that clade-D ARFs are transcriptional activators in the model bryophyte Physcomitrium patens and have a major role in the development of this species. Δarfddub protonemata exhibit a delay in filament branching, as well as a delay in the chloronema to caulonema transition. Additionally, leafy gametophore development in Δarfddub lines lags behind wild type. We present evidence that ARFd1 interacts with activating ARFs via their PB1 domains, but not with repressing ARFs. Based on these results, we propose a model in which clade-D ARFs enhance gene expression by interacting with DNA bound clade-A ARFs. Further, we show that ARFd1 must form oligomers for full activity.

Completing the TRB family: newly characterized members show ancient evolutionary origins and distinct localization, yet similar interactions

Telomere repeat binding proteins (TRBs) belong to a family of proteins possessing a Myb-like domain which binds to telomeric repeats. Three members of this family (TRB1, TRB2, TRB3) from Arabidopsis thaliana have already been described as associated with terminal telomeric repeats (telomeres) or short interstitial telomeric repeats in gene promoters (telo-boxes). They are also known to interact with several protein complexes: telomerase, Polycomb repressive complex 2 (PRC2) E(z) subunits and the PEAT complex (PWOs-EPCRs-ARIDs-TRBs). Here we characterize two novel members of the TRB family (TRB4 and TRB5). Our wide phylogenetic analyses have shown that TRB proteins evolved in the plant kingdom after the transition to a terrestrial habitat in Streptophyta, and consequently TRBs diversified in seed plants. TRB4-5 share common TRB motifs while differing in several others and seem to have an earlier phylogenetic origin than TRB1-3. Their common Myb-like domains bind long arrays of telomeric repeats in vitro, and we have determined the minimal recognition motif of all TRBs as one telo-box. Our data indicate that despite the distinct localization patterns of TRB1-3 and TRB4-5 in situ, all members of TRB family mutually interact and also bind to telomerase/PRC2/PEAT complexes. Additionally, we have detected novel interactions between TRB4-5 and EMF2 and VRN2, which are Su(z)12 subunits of PRC2.

Co-action of COP1, SPA and cryptochrome in light signal transduction and photomorphogenesis of the moss Physcomitrium patens

The Arabidopsis COP1/SPA ubiquitin ligase suppresses photomorphogenesis in darkness. In the light, photoreceptors inactivate COP1/SPA to allow a light response. While SPA genes are specific to the green lineage, COP1 also exists in humans. This raises the question of when in evolution plant COP1 acquired the need for SPA accessory proteins. We addressed this question by generating Physcomitrium Ppcop1 mutants and comparing their visible and molecular phenotypes with those of Physcomitrium Ppspa mutants. The phenotype of Ppcop1 nonuple mutants resembles that of Ppspa mutants. Most importantly, both mutants produce green chloroplasts in complete darkness. They also exhibit dwarfed gametophores, disturbed branching of protonemata and absent gravitropism. RNA-sequencing analysis indicates that both mutants undergo weak constitutive light signaling in darkness. PpCOP1 and PpSPA proteins form a complex and they interact via their WD repeat domains with the VP motif of the cryptochrome CCE domain in a blue light-dependent manner. This resembles the interaction of Arabidopsis SPA proteins with Arabidopsis CRY1, and is different from that with Arabidopsis CRY2. Taken together, the data indicate that PpCOP1 and PpSPA act together to regulate growth and development of Physcomitrium. However, in contrast to their Arabidopsis orthologs, PpCOP1 and PpSPA proteins execute only partial suppression of light signaling in darkness. Hence, additional repressors may exist that contribute to the repression of a light response in dark-exposed Physcomitrium.