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

Recombinant production of spider silk protein in Physcomitrella photobioreactors

KEY MESSAGE

We report the successful moss-produced recombinant spider silk key protein component containing both the N- and the C-terminal domain. Spider dragline silk stands out as a remarkable biomaterial, representing one of nature's toughest fibres. Its strength rivals that of many synthetic fibres used commercially, rendering it applicable across various industrial and medical domains. However, its widespread utilisation requires cost-effective mass production. Biotechnology presents a promising avenue for achieving this goal, particularly through the production of recombinant dragline silk proteins in transgenic plant systems. This study aimed to assess the feasibility of producing one key protein component of dragline silk, MaSp1, from the western black widow spider, Latrodectus hesperus, the protein LhMaSp1, in the moss Physcomitrella (Physcomitrium patens). Here, we present the successful recombinant production of spider silk protein containing both the N- and C-terminal domains of LhMaSp1 in moss cells. The production of recombinant LhMaSp1 protein in Physcomitrella was performed in shake flasks and in five-litre photobioreactors and the correct synthesis of LhMaSp1 was proven via mass spectrometry. We estimate that the yield of recombinant spider silk protein in Physcomitrella bioreactors is above 0.82 mg/g fresh weight.

Accurate cross-species 5mC detection for Oxford Nanopore sequencing in plants with DeepPlant

Nanopore sequencing enables comprehensive detection of 5-methylcytosine (5mC), particularly in repeat regions. However, CHH methylation detection in plants is limited by the scarcity of high-methylation positive samples, reducing generalization across species. Dorado, the only tool for plant 5mC detection on the R10.4 platform, lacks extensive species testing. Here, we develop DeepPlant, a deep learning model incorporating both Bi-LSTM and Transformer architectures, which significantly improves CHH detection accuracy and performs well for CpG and CHG motifs. We address the scarcity of methylation-positive CHH training samples through screening species with abundant high-methylation CHH sites using bisulfite-sequencing and generate datasets that cover diverse 9-mer motifs for training and testing DeepPlant. Evaluated across nine species, DeepPlant achieves high whole-genome methylation frequency correlations (0.705-0.838) with BS-seq data on CHH, improved by 23.4- 117.6% compared to Dorado. DeepPlant also demonstrates superior single-molecule accuracy and F1 score, offering strong generalization for plant epigenetics research.

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

The plant hormone auxin regulates a wide variety of transcriptional responses depending on the cell type, environment and species. How this diversity is achieved may be related to the specific complement of auxin-signalling components in each cell. The levels of activators (class-A AUXIN RESPONSE FACTORS) and repressors (class-B ARFs) are particularly important. Tight regulation of ARF protein levels is probably key in determining this balance. 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 mediated by proteasome degradation, 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.

Functional genomic perspectives on plant terrestrialization

Plant evolutionary research has made leaps in exploring the deep evolutionary roots of embryophytes. A solid phylogenomic framework was established, allowing evolutionary inferences. Comparative genomic approaches revealed that many genes coding for transcription factors, morphogenetic regulators, specialized metabolic enzymes, phytohormone signaling, and more are not innovations of land plants but have a deep streptophyte algal ancestry. Are these just spurious homologs, or do they actualize traits we deem important in embryophytes? Building on streptophyte algae genome data, current endeavors delve into the functional significance of whole cohorts of homologs by leveraging the power of comparative high-throughput approaches. This ushered in the identification of recurrent themes in function, ultimately providing a functional genomic definition for the toolkit of plant terrestrialization.

Tandem LTR-retrotransposon structures are common and highly polymorphic in plant genomes

BACKGROUND

LTR-retrotransposons (LTR-RT) are a major component of plant genomes and important drivers of genome evolution. Most LTR-RT copies in plant genomes are defective elements found as truncated copies, nested insertions or as part of more complex structures. The recent availability of highly contiguous plant genome assemblies based on long-read sequences now allows to perform detailed characterization of these complex structures and to evaluate their importance for plant genome evolution.

RESULTS

The detailed analysis of two rice loci containing complex LTR-RT structures showed that they consist of tandem arrays of LTR copies sharing internal LTRs. Our analyses suggests that these LTR-RT tandems are the result of a single insertion and not of the recombination of two independent LTR-RT elements. Our results also suggest that gypsy elements may be more prone to form these structures. We show that these structures are highly polymorphic in rice and therefore have the potential to generate genetic variability. We have developed a computational pipeline (IDENTAM) that scans genome sequences and identifies tandem LTR-RT candidates. Using this tool, we have detected 266 tandems in a pangenome built from the genomes of 76 accessions of cultivated and wild rice, showing that tandem LTR-RT structures are frequent and highly polymorphic in rice. Running IDENTAM in the Arabidopsis, almond and cotton genomes showed that LTR-RT tandems are frequent in plant genomes of different size, complexity and ploidy level. The complexity of differentiating intra-element variations at the nucleotide level among haplotypes is very high, and we found that graph-based pangenomic methodologies are appropriate to resolve these structures.

CONCLUSIONS

Our results show that LTR-RT elements can form tandem arrays. These structures are relatively abundant and highly polymorphic in rice and are widespread in the plant kingdom. Future studies will contribute to understanding how these structures originate and whether the variability that they generate has a functional impact.

The Microtubule Cytoskeleton in Bryophytes

Microtubules (MTs) are essential cytoskeletal elements in all eukaryotes, playing critical roles in cell shape, intercellular organization, cell division, and cell motility. The organization of the MT network has undergone significant changes throughout plant evolution. Some MT structures, such as the preprophase band and phragmoplast, are innovations in plant lineages, while others, including the centriole and flagellum, have been lost over time. Bryophytes, consisting of mosses, liverworts, and hornworts, are the earliest land plants and occupy a key phylogenetic position in the evolution of MT organization. In the past two decades, advances in genomics, genetics, and cell imaging technologies have significantly enhanced our understanding of MT organization and function. Two representative species, Physcomitrium patens (moss) and Marchantia polymorph (liverwort), have become established model organisms, and new models for hornworts are emerging. In this review, we summarize the current knowledge of the MT cytoskeleton, drawing from early electron microscopy studies and recent advances in these emerging models. Our aim is to provide a comprehensive overview of the major MT array types and key factors involved in MT organization in bryophytes, offering insights into MT adaptation during plant evolution.

Ancient large-scale gene duplications and diversification in bryophytes illuminate the plant terrestrialization

Large-scale gene duplications (LSGDs) are crucial for evolutionary adaptation and recurrent in vascular plants. However, the role of ancient LSGDs in the terrestrialization and diversification of bryophytes, the second most species-rich group of land plants, remains largely elusive due to limited sampling in bryophytes. Employing the most extensive nuclear gene dataset in bryophytes to date, we reconstructed a time-calibrated phylogenetic tree from 209 species, covering virtually all key bryophyte lineages, for phylogenomic analyses of LSGDs and diversification. We newly identified two ancient LSGDs: one in the most recent common ancestor (MRCA) of extant bryophytes and another in the MRCA of the majority of Jungermanniales s. lato. Duplicated genes from these two LSGDs show significant enrichment in photosynthesis-related processes and structures. Rhizoid-responsive ROOTHAIR DEFECTIVE SIX-LIKE (RSL) genes from ancient LSGDs are present in rhizoidless bryophytes, challenging assumptions about rhizoid absence mechanisms. We highlighted four major diversification rate upshifts, two of which slightly postdated LSGDs, potentially linked to the flourishing of gymnosperms and angiosperms and explaining over 80% of bryophyte diversity. Our findings, supported by extensive bryophyte sampling, highlight the significance of LSGDs in the early terrestrialization and diversification of bryophytes, offering new insights into land plant evolution.

The Marchantia polymorpha pangenome reveals ancient mechanisms of plant adaptation to the environment

Plant adaptation to terrestrial life started 450 million years ago and has played a major role in the evolution of life on Earth. The genetic mechanisms allowing this adaptation to a diversity of terrestrial constraints have been mostly studied by focusing on flowering plants. Here, we gathered a collection of 133 accessions of the model bryophyte Marchantia polymorpha and studied its intraspecific diversity using selection signature analyses, a genome-environment association study and a pangenome. We identified adaptive features, such as peroxidases or nucleotide-binding and leucine-rich repeats (NLRs), also observed in flowering plants, likely inherited from the first land plants. The M. polymorpha pangenome also harbors lineage-specific accessory genes absent from seed plants. We conclude that different land plant lineages still share many elements from the genetic toolkit evolved by their most recent common ancestor to adapt to the terrestrial habitat, refined by lineage-specific polymorphisms and gene family evolution.

Comparative analysis using a chromosome-scale genome assembly for Funaria hygrometrica suggests greater collinearity in mosses than in seed plants

Mosses, the largest lineage of seed-free plants, have smaller and less variable genome sizes than flowering plants. Nevertheless, whether this difference results from divergent genome dynamics is poorly known. Here, we use newly generated chromosome-scale genome assemblies for Funaria hygrometrica and comparative analysis with other moss and seed plant genomes to investigate moss genome dynamics. Although some aspects of moss genome dynamics are seed plant-like, such as the mechanism of genome size change and de novo gain/loss of genes, moss genomes retain higher synteny, and collinearity over evolutionary time than seed plant genomes. Furthermore, transposable elements and genes are more evenly distributed along chromosomes in mosses than in seed plants, a feature shared with other sequenced seed-free plant genomes. Overall, our findings support the hypothesis that large-scale genome structure and dynamics of mosses and seed plants differ. In particular, our data suggest a lower rate of gene order reshuffling along chromosomes in mosses compared to seed plants. We speculate that such lower rate of structural genomic variation and unique chromosome structure in mosses may contribute to their relatively smaller and less variable genome sizes.

Insight into the Characterization of Two Female Suppressor Gene Families: SOFF and SyGI in Plants

BACKGROUND/OBJECTIVES

The Suppressor of Female Function (SOFF) and Shy Girl (SyGI) gene families play vital roles in sex determination in dioecious plants. However, their evolutionary dynamics and functional characteristics remain largely unexplored.

METHODS

Through this study, a systematic bioinformatics analysis of SOFF and SyGI families was performed in plants to explore their evolutionary relationships, gene structures, motif synteny and functional predictions.

RESULTS

Phylogenetic analysis showed that the SOFF family expanded over time and was divided into two subfamilies and seven groups, while SyGI was a smaller family made of compact molecules with three groups. Synteny analysis revealed that 125 duplicated gene pairs were identified in Kiwifruit where WGD/segmental duplication played a major role in duplicating these events. Structural analysis predicted that SOFF genes have a DUF 247 domain with a transmembrane region, while SyGI sequences have an REC-like conserved domain, with a "barrel-shaped" structure consisting of five α-helices and five β-strands. Promoter region analysis highlighted their probable regulatory roles in plant development, hormone signaling and stress responses. Protein interaction analysis exhibited only four SOFF genes with a close interaction with other genes, while SyGI genes had extensive interactions, particularly with cytokinin signal transduction pathways.

CONCLUSIONS

The current study offers a crucial understanding of the molecular evolution and functional characteristics of SOFF and SyGI gene families, providing a foundation for future functional validation and genetic studies on developmental regulation and sex determination in dioecious plants. Also, this research enhances our insight into plant reproductive biology and offers possible targets for breeding and genetic engineering approaches.

Evidence for evolution of a new sex chromosome within the haploid-dominant Marchantiales plant lineage

Sex chromosomes have evolved independently in numerous lineages across the Tree of Life, in both diploid-dominant species, including many animals and plants, and the less studied haploid-dominant plants and algae. Strict genetic sex determination ensures that individuals reproduce by outcrossing. However, species with separate sexes (termed dioecy in diploid plants, and dioicy in haploid plants) may sometimes evolve different sex systems, and become monoicous, with the ability to self-fertilize. Here, we studied dioicy-monoicy transitions in the ancient liverwort haploid-dominant plant lineage, using three telomere-to-telomere gapless chromosome-scale reference genome assemblies from the Ricciaceae group of Marchantiales. Ancestral liverworts are believed to have been dioicous, with U and V chromosomes (chromosome 9) determining femaleness and maleness, respectively. We confirm the finding that monoicy in Ricciocarpos natans evolved from a dioicous ancestor, and most ancestrally U chromosomal genes have been retained on autosomes in this species. We also describe evidence suggesting the possible re-evolution of dioicy in the genus Riccia, with probable de novo establishment of a sex chromosome from an autosome (chromosome 5), and further translocations of genes from the new sex chromosome to autosomes. Our results also indicated that micro-chromosomes are consistent genomic features, and may have evolved independently from sex chromosomes in Ricciocarpos and Riccia lineages.

Time-resolved oxidative signal convergence across the algae-embryophyte divide

The earliest land plants faced a significant challenge in adapting to environmental stressors. Stress on land is unique in its dynamics, entailing swift and drastic changes in light and temperature. While we know that land plants share with their closest streptophyte algal relatives key components of the genetic makeup for dynamic stress responses, their concerted action is little understood. Here, we combine time-course stress profiling using photophysiology, transcriptomics on 2.7 Tbp of data, and metabolite profiling analyses on 270 distinct samples, to study stress kinetics across three 600-million-year-divergent streptophytes. Through co-expression analysis and Granger causal inference we predict a gene regulatory network that retraces a web of ancient signal convergences at ethylene signaling components, osmosensors, and chains of major kinases. These kinase hubs already integrated diverse environmental inputs since before the dawn of plants on land.

Differential GTP-dependent in-vitro polymerization of recombinant Physcomitrella FtsZ proteins

Bacterial cell division and plant chloroplast division require selfassembling Filamentous temperature-sensitive Z (FtsZ) proteins. FtsZ proteins are GTPases sharing structural and biochemical similarities with eukaryotic tubulin. In the moss Physcomitrella, the morphology of the FtsZ polymer networks varies between the different FtsZ isoforms. The underlying mechanism and foundation of the distinct networks is unknown. Here, we investigated the interaction of Physcomitrella FtsZ2-1 with FtsZ1 isoforms via co-immunoprecipitation and mass spectrometry, and found protein-protein interaction in vivo. We tagged FtsZ1-2 and FtsZ2-1 with different fluorophores and expressed both in E. coli, which led to the formation of defined structures within the cells and to an influence on bacterial cell division and morphology. Furthermore, we have optimized the purification protocols for FtsZ1-2 and FtsZ2-1 expressed in E. coli and characterized their GTPase activity and polymerization in vitro. Both FtsZ isoforms showed GTPase activity. Stoichiometric mixing of both proteins led to a significantly increased GTPase activity, indicating a synergistic interaction between them. In light scattering assays, we observed GTP-dependent assembly of FtsZ1-2 and of FtsZ2-1 in a protein concentration dependent manner. Stoichiometric mixing of both proteins resulted in significantly faster polymerization, again indicating a synergistic interaction between them. Under the same conditions used for GTPase and light scattering assays both FtsZ isoforms formed filaments in a GTP-dependent manner as visualized by transmission electron microscopy (TEM). Taken together, our results reveal that Physcomitrella FtsZ1-2 and FtsZ2-1 are functionally different, can synergistically interact in vivo and in vitro, and differ in their properties from FtsZ proteins from bacteria, archaea and vascular plants.

Pan-phylum genomes of hornworts reveal conserved autosomes but dynamic accessory and sex chromosomes

Hornworts, one of the three bryophyte phyla, show some of the deepest divergences in extant land plants, with some families separated by more than 300 million years. Previous hornwort genomes represented only one genus, limiting the ability to infer evolution within hornworts and their early land plant ancestors. Here we report ten new chromosome-scale genomes representing all hornwort families and most of the genera. We found that, despite the deep divergence, synteny was surprisingly conserved across all hornwort genomes, a pattern that might be related to the absence of whole-genome duplication. We further uncovered multiple accessory and putative sex chromosomes that are highly repetitive and CpG methylated. In contrast to autosomes, these chromosomes mostly lack syntenic relationships with one another and are evolutionarily labile. Notable gene retention and losses were identified, including those responsible for flavonoid biosynthesis, stomata patterning and phytohormone reception, which have implications in reconstructing the evolution of early land plants. Together, our pan-phylum genomes revealed an array of conserved and divergent genomic features in hornworts, highlighting the uniqueness of this deeply diverged lineage.

Enhanced sensitivity of TAPscan v4 enables comprehensive analysis of streptophyte transcription factor evolution

Transcription-associated proteins (TAPs) fulfill multiple functions in regulatory and developmental processes and display lineage-specific evolution. TAPscan is a comprehensive and highly reliable tool for genome-wide TAP annotation via domain profiles. Here, we present TAPscan v4, including an updated web interface (https://tapscan.plantcode.cup.uni-freiburg.de/), which enables an in-depth representation of the distribution of 138 TAP families across 678 species from diverse groups of organisms, with a focus on Archaeplastida (plants in the wide sense). With this release, we also make the underlying "Genome Zoo" available, a curated protein data set with scripts and metadata. Eighteen new TAP (sub)families were added as part of the update. Nine of those were gained in the most recent common ancestor of the Streptophyta (comprising streptophyte algae and land plants), or within the streptophyte algae. More than one-third of all detected TAP family gains were identified during the evolution of streptophyte algae, before the emergence of land plants, and are thus likely to have been significant for plant terrestrialization. The TAP complement of the Zygnematophyceae was identified to be the most similar to that of land plants, consistent with the finding that this lineage is sister to land plants. Overall, our data retrace the evolution of streptophyte TAPs, allowing us to pinpoint the regulatory repertoire of the earliest land plants.

The evolutionarily diverged single-stranded DNA-binding proteins SSB1/SSB2 differentially affect the replication, recombination and mutation of organellar genomes in Arabidopsis thaliana

Single-stranded DNA-binding proteins (SSBs) play essential roles in the replication, recombination and repair processes of organellar DNA molecules. In Arabidopsis thaliana, SSBs are encoded by a small family of two genes (SSB1 and SSB2). However, the functional divergence of these two SSB copies in plants remains largely unknown, and detailed studies regarding their roles in the replication and recombination of organellar genomes are still incomplete. In this study, phylogenetic, gene structure and protein motif analyses all suggested that SSB1 and SSB2 probably diverged during the early evolution of seed plants. Based on accurate long-read sequencing results, ssb1 and ssb2 mutants had decreased copy numbers for both mitochondrial DNA (mtDNA) and plastid DNA (ptDNA), accompanied by a slight increase in structural rearrangements mediated by intermediate-sized repeats in mt genome and small-scale variants in both genomes. Our findings provide an important foundation for further investigating the effects of DNA dosage in the regulation of mutation frequencies in plant organellar genomes.

An alternate route for cellulose microfibril biosynthesis in plants

Similar to cellulose synthases (CESAs), cellulose synthase-like D (CSLD) proteins synthesize β-1,4-glucan in plants. CSLDs are important for tip growth and cytokinesis, but it was unknown whether they form membrane complexes in vivo or produce microfibrillar cellulose. We produced viable CESA-deficient mutants of the moss Physcomitrium patens to investigate CSLD function without interfering CESA activity. Microscopy and spectroscopy showed that CESA-deficient mutants synthesize cellulose microfibrils that are indistinguishable from those in vascular plants. Correspondingly, freeze-fracture electron microscopy revealed rosette-shaped particle assemblies in the plasma membrane that are indistinguishable from CESA-containing rosette cellulose synthesis complexes (CSCs). Our data show that proteins other than CESAs, most likely CSLDs, produce cellulose microfibrils in P. patens protonemal filaments. The data suggest that the specialized roles of CSLDs in cytokinesis and tip growth are based on differential expression and different interactions with microtubules and possibly Ca2+, rather than structural differences in the microfibrils they produce.

Sex determination in bryophytes: current state of the art

With the advent of genomic and other omics technologies, the last decades have witnessed a series of steady and important breakthroughs in the understanding of genetic determinants of different reproductive systems in vascular plants and especially on how sexual reproduction shaped their evolution. In contrast, the molecular mechanisms of these fundamental aspects of the biology of bryophytes, a group of non-vascular embryophyte plants sister to all tracheophytes, are still largely obscure. The recent characterization of the sex chromosomes and genetic switches determining sex in bryophytes and emerging approaches for molecular sexing of gametophytes hold great promise for elucidation of the evolutionary history as well as the conservation of this species-rich but understudied group of land plants.

The plant early recombinosome: a high security complex to break DNA during meiosis

KEY MESSAGE

The formacion of numerous unpredictable DNA Double Strand Breaks (DSBs) on chromosomes iniciates meiotic recombination. In this perspective, we propose a 'multi-key lock' model to secure the risky but necesary breaks as well as a 'one per pair of cromatids' model for the topoisomerase-like early recombinosome. During meiosis, homologous chromosomes recombine at few sites of crossing-overs (COs) to ensure correct segregation. The initiation of meiotic recombination involves the formation of DNA double strand breaks (DSBs) during prophase I. Too many DSBs are dangerous for genome integrity: if these DSBs are not properly repaired, it could potentially lead to chromosomal fragmentation. Too few DSBs are also problematic: if the obligate CO cannot form between bivalents, catastrophic unequal segregation of univalents lead to the formation of sterile aneuploid spores. Research on the regulation of the formation of these necessary but risky DSBs has recently advanced in yeast, mammals and plants. DNA DSBs are created by the enzymatic activity of the early recombinosome, a topoisomerase-like complex containing SPO11. This opinion paper reviews recent insights on the regulation of the SPO11 cofactors necessary for the introduction of temporally and spatially controlled DSBs. We propose that a 'multi-key-lock' model for each subunit of the early recombinosome complex is required to secure the formation of DSBs. We also discuss the hypothetical implications that the established topoisomerase-like nature of the SPO11 core-complex can have in creating DSB in only one of the two replicated chromatids of early prophase I meiotic chromosomes. This hypothetical 'one per pair of chromatids' DSB formation model could optimize the faithful repair of the self-inflicted DSBs. Each DSB could use three potential intact homologous DNA sequences as repair template: one from the sister chromatid and the two others from the homologous chromosomes.

MpPUB9, a U-box E3 ubiquitin ligase, acts as a positive regulator by promoting the turnover of MpEXO70.1 under high salinity in Marchantia polymorpha

Marchantia polymorpha, occupying a basal position in the monophyletic assemblage of land plants, displays a notable expansion of plant U-box (PUB) proteins compared with those in animals. We elucidated the roles of MpPUB9 in regulating salt stress tolerance in M. polymorpha. MpPUB9 expression was rapidly induced by high salinity and dehydration. MpPUB9 possessed an intact U-box domain in the N-terminus. MpPUB9-Citrine localized to punctate structures and was peripherally associated with microsomal membranes. Phenotypic analyses demonstrate that the hyponastic and epinastic thallus growth phenotypes, which were induced by the overexpression and suppression of MpPUB9, may provoke salt stress-resistant and -susceptible phenotypes, respectively. MpPUB9 was also found to directly interact with the exocyst protein MpEXO70.1, leading to its ubiquitination. Under high-salinity conditions, though the stability of MpPUB9 was dramatically increased, MpEXO70.1 showed slightly faster turnover rates. Transcriptome analyses showed that salt treatment and the overexpression of MpPUB9 co-upregulated the genes related to the modulation of H2O2 and cell wall organization. Overall, our results suggest that MpPUB9 plays a crucial role in the positive regulation of salt stress tolerance, resulting from its interaction with MpEXO70.1 and modulating turnover of the protein under high-salt conditions via the coordination of UPS with autophagy.

Differential prolyl hydroxylation by six Physcomitrella prolyl-4 hydroxylases

Hydroxylation of prolines to 4-trans-hydroxyproline (Hyp) is mediated by prolyl-4 hydroxylases (P4Hs). In plants, Hyps occur in Hydroxyproline-rich glycoproteins (HRGPs), and are frequently O-glycosylated. While both modifications are important, e.g. for cell wall stability, they are undesired in plant-made pharmaceuticals. Sequence motifs for prolyl-hydroxylation were proposed but did not include data from mosses, such as Physcomitrella. We identified six moss P4Hs by phylogenetic reconstruction. Our analysis of 73 Hyps in 24 secretory proteins from multiple mass spectrometry datasets revealed that prolines near other prolines, alanine, serine, threonine and valine were preferentially hydroxylated. About 95 % of Hyps were predictable with combined established methods. In our data, AOV was the most frequent pattern. A combination of 443 AlphaFold models and MS data with 3000 prolines found Hyps mainly on protein surfaces in disordered regions. Moss-produced human erythropoietin (EPO) exhibited O-glycosylation with arabinose chains on two Hyps. This modification was significantly reduced in a p4h1 knock-out (KO) Physcomitrella mutant. Quantitative proteomics with different p4h mutants revealed specific changes in protein amounts, and a modified prolyl-hydroxylation pattern, suggesting a differential function of the Physcomitrella P4Hs. Quantitative RT-PCR revealed a differential effect of single p4h KOs on the expression of the other five p4h genes, suggesting a partial compensation of the mutation. AlphaFold-Multimer models for Physcomitrella P4H1 and its target EPO peptide superposed with the crystal structure of Chlamydomonas P4H1 suggested significant amino acids in the active centre of the enzyme and revealed differences between P4H1 and the other Physcomitrella P4Hs.

Current and future perspectives for enhancing our understanding of the evolution of plant metabolism

The special issue 'The evolution of plant metabolism' has brought together original research, reviews and opinions that cover various aspects from the full breath of plant metabolism including its interaction with the environment including other species. Here, we briefly summarize these efforts and attempts to extract a consensus opinion of the best manner in which to tackle this subject both now and in the future. This article is part of the theme issue 'The evolution of plant metabolism'.

Conserved carotenoid pigmentation in reproductive organs of Charophyceae

Sexual reproduction in Charophyceae abounds in complex traits. Their gametangia develop as intricate structures, with oogonia spirally surrounded by envelope cells and richly pigmented antheridia. The red-probably protectant-pigmentation of antheridia is conserved across Charophyceae. Chara tomentosa is, however, unique in exhibiting this pigmentation and also in vegetative tissue. Here, we investigated the two sympatric species, C. tomentosa and Chara baltica, and compared their molecular chassis for pigmentation. Using reversed phase C30 high performance liquid chromatography (RP-C30-HPLC), we uncover that the major pigments are β-carotene, δ-carotene and γ-carotene; using headspace solid-phase microextraction coupled to gas chromatography equipped with a mass spectrometer (HS-SPME-GC-MS), we pinpoint that the unusually large carotenoid pool in C. tomentosa gives rise to diverse volatile apocarotenoids, including abundant 6-methyl-5-hepten-2-one. Based on transcriptome analyses, we uncover signatures of the unique biology of Charophycaee and genes for pigment production, including monocyclized carotenoids. The rich carotenoid pool probably serves as a substrate for diverse carotenoid-derived metabolites, signified not only by (i) the volatile apocarotenoids we detected but (ii) the high expression of a gene coding for a cytochrome P450 enzyme related to land plant proteins involved in the biosynthesis of carotenoid-derived hormones. Overall, our data shed light on a key protection strategy of sexual reproduction in the widespread group of macroalgae. The genetic underpinnings of this are shared across hundreds of millions of years of plant and algal evolution. This article is part of the theme issue 'The evolution of plant metabolism'.

Evolution of the regulatory subunits for the heteromeric acetyl-CoA carboxylase

The committed step for de novo fatty acid (FA) synthesis is the ATP-dependent carboxylation of acetyl-coenzyme A catalysed by acetyl-CoA carboxylase (ACCase). In most plants, ACCase is a multi-subunit complex orthologous to prokaryotes. However, unlike prokaryotes, the plant and algal orthologues are comprised both catalytic and additional dedicated regulatory subunits. Novel regulatory subunits, biotin lipoyl attachment domain-containing proteins (BADC) and carboxyltransferase interactors (CTI) (both three-gene families in Arabidopsis) represent new effectors specific to plants and certain algal species. The evolutionary history of these genes in autotrophic eukaryotes remains elusive, making it an ongoing area of research. Analyses of potential protein-protein and co-occurrence interactions, informed by gene network patterns using the STRING database, in Arabidopsis thaliana and Chlamydomonas reinhardtii unveil intricate gene associations with ACCase, suggesting a complex interplay between FA synthesis and other cellular processes. Among both species, a higher number of co-expressed genes was identified in Arabidopsis, indicating a wider potential regulatory network of ACCase in plants. This review investigates the extent to which these genes arose in autotrophic eukaryotes and provides insights into their evolutionary trajectory. This article is part of the theme issue 'The evolution of plant metabolism'.

Evolution of NAC transcription factors from early land plants to domesticated crops

NAC [NO APICAL MERISTEM (NAM), ARABIDOPSIS TRANSCRIPTION ACTIVATOR FACTOR 1/2 (ATAF1/2), and CUP-SHAPED COTYLEDON (CUC2)] transcription factors are key regulators of plant growth, development, and stress responses but were also crucial players during land plant adaptation and crop domestication. Using representative members of green algae, bryophytes, lycophytes, gymnosperms, and angiosperms, we expanded the evolutionary history of NAC transcription factors to unveil the relationships among members of this gene family. We found a massive increase in the number of NAC transcription factors from green algae to lycophytes and an even larger increase in flowering plants. Many of the NAC clades arose later during evolution since we found eudicot- and monocot-specific clades. Cis-elements analysis in NAC promoters showed the presence of abiotic and biotic stress as well as hormonal response elements, which indicate the ancestral function of NAC transcription factor genes in response to environmental stimuli and in plant development. At the transcriptional level, the expression of NAC transcription factors was low or absent in male reproduction, particularly mature pollen, across the plant kingdom. We also identified NAC genes with conserved expression patterns in response to heat stress in Marchantia polymorpha and Oryza sativa. Our study provides further evidence that transcriptional mechanisms associated with stress responses and development emerged early during plant land adaptation and are still conserved in flowering plants and domesticated crops.

Origin, Evolution, and Diversification of the Expansin Family in Plants

The cell wall is a crucial feature that allows ancestral streptophyte green algae to colonize land. Expansin, an extracellular protein that mediates cell wall loosening in a pH-dependent manner, could be a powerful tool for studying cell wall evolution. However, the evolutionary trajectory of the expansin family remains largely unknown. Here, we conducted a comprehensive identification of 2461 expansins across 64 sequenced species, ranging from aquatic algae to terrestrial plants. Expansins originated in chlorophyte algae and may have conferred the ability to loosen cell walls. The four expansin subfamilies originated independently: α-expansin appeared first, followed by β-expansin, and then expansin-like A and expansin-like B, reflecting the evolutionary complexity of plant expansins. Whole genome duplication/segmental duplication and tandem duplication events greatly contributed to expanding the expansin family. Despite notable changes in sequence characteristics, the intron distribution pattern remained relatively conserved among different subfamilies. Phylogenetic analysis divided all the expansins into five clades, with genes from the same subfamily tending to cluster together. Transcriptome data from 16 species across ten lineages and qRT-PCR analysis revealed varying expression patterns of expansin genes, suggesting functional conservation and diversification during evolution. This study enhances our understanding of the evolutionary conservation and dynamics of the expansin family in plants, providing insight into their roles as cell wall-loosening factors.

Reliability of plastid and mitochondrial localisation prediction declines rapidly with the evolutionary distance to the training set increasing

Mitochondria and plastids import thousands of proteins. Their experimental localisation remains a frequent task, but can be resource-intensive and sometimes impossible. Hence, hundreds of studies make use of algorithms that predict a localisation based on a protein's sequence. Their reliability across evolutionary diverse species is unknown. Here, we evaluate the performance of common algorithms (TargetP, Localizer and WoLFPSORT) for four photosynthetic eukaryotes (Arabidopsis thaliana, Zea mays, Physcomitrium patens, and Chlamydomonas reinhardtii) for which experimental plastid and mitochondrial proteome data is available, and 171 eukaryotes using orthology inferences. The match between predictions and experimental data ranges from 75% to as low as 2%. Results worsen as the evolutionary distance between training and query species increases, especially for plant mitochondria for which performance borders on random sampling. Specificity, sensitivity and precision analyses highlight cross-organelle errors and uncover the evolutionary divergence of organelles as the main driver of current performance issues. The results encourage to train the next generation of neural networks on an evolutionary more diverse set of organelle proteins for optimizing performance and reliability.

Plant biology: Mapping meristem morphogenesis in Marchantia

The use of state-of-the-art imaging, underpinned by molecular data, for the first time provides a clear understanding of two fundamental processes in liverworts - the establishment of dorsoventrality and origin of apical meristems. This work opens the door to exploring many new facets of plant morphogenesis.

A snapshot of the Physcomitrella N-terminome reveals N-terminal methylation of organellar proteins

KEY MESSAGE

Analysis of the N-terminome of Physcomitrella reveals N-terminal monomethylation of nuclear-encoded, mitochondria-localized proteins. Post- or co-translational N-terminal modifications of proteins influence their half-life as well as mediating protein sorting to organelles via cleavable N-terminal sequences that are recognized by the respective translocation machinery. Here, we provide an overview on the current modification state of the N-termini of over 4500 proteins from the model moss Physcomitrella (Physcomitrium patens) using a compilation of 24 N-terminomics datasets. Our data reveal distinct proteoforms and modification states and confirm predicted targeting peptide cleavage sites of 1,144 proteins localized to plastids and the thylakoid lumen, to mitochondria, and to the secretory pathway. In addition, we uncover extended N-terminal methylation of mitochondrial proteins. Moreover, we identified PpNTM1 (P. patens alpha N-terminal protein methyltransferase 1) as a candidate for protein methylation in plastids, mitochondria, and the cytosol. These data can now be used to optimize computational targeting predictors, for customized protein fusions and their targeted localization in biotechnology, and offer novel insights into potential dual targeting of proteins.

Insights into a functional synthetic plant genome

Synthetic genomics involves the design, assembly, and transfer of artificially synthesized DNA fragments into target hosts to replace the native genome and construct viable forms of life. With advances in DNA synthesis and assembly techniques, the application of synthetic genomics in viruses, bacteria, and yeast has improved our knowledge of genome organization and function. Multicellular eukaryotic organisms are characterized by larger genomes, more complex epigenetic regulation, and widespread transposable elements, making genome synthesis challenging. Recently, the first synthetic multicellular eukaryotic organism was generated in the model plant Physcomitrium patens with a partially synthetic chromosome arm. Here, we introduce the design and assembly principles of moss genome synthesis. We also discuss the remaining technical barriers in the application of synthetic genomics in seed plants.

Regulation of heterochromatin organization in plants

Heterochromatin is a nuclear area that contains highly condensed and transcriptionally inactive chromatin. Alterations in the organization of heterochromatin are correlated with changes in gene expression and genome stability, which affect various aspects of plant life. Thus, studies of the molecular mechanisms that regulate heterochromatin organization are important for understanding the regulation of plant physiology. Microscopically, heterochromatin can be characterized as chromocenters that are intensely stained with DNA-binding fluorescent dyes. Arabidopsis thaliana exhibits distinctive chromocenters in interphase nuclei, and genetic studies combined with cytological analyses have identified a number of factors that are involved in heterochromatin assembly and organization. In this review, I will summarize the factors involved in the regulation of heterochromatin organization in plants.

Contrasting and conserved roles of NPR pathways in diverged land plant lineages

The NPR proteins function as salicylic acid (SA) receptors in Arabidopsis thaliana. AtNPR1 plays a central role in SA-induced transcriptional reprogramming whereby positively regulates SA-mediated defense. NPRs are found in the genomes of nearly all land plants. However, we know little about the molecular functions and physiological roles of NPRs in most plant species. We conducted phylogenetic and alignment analyses of NPRs from 68 species covering the significant lineages of land plants. To investigate NPR functions in bryophyte lineages, we generated and characterized NPR loss-of-function mutants in the liverwort Marchantia polymorpha. Brassicaceae NPR1-like proteins have characteristically gained or lost functional residues identified in AtNPRs, pointing to the possibility of a unique evolutionary trajectory for the Brassicaceae NPR1-like proteins. We find that the only NPR in M. polymorpha, MpNPR, is not the master regulator of SA-induced transcriptional reprogramming and negatively regulates bacterial resistance in this species. The Mpnpr transcriptome suggested roles of MpNPR in heat and far-red light responses. We identify both Mpnpr and Atnpr1-1 display enhanced thermomorphogenesis. Interspecies complementation analysis indicated that the molecular properties of AtNPR1 and MpNPR are partially conserved. We further show that MpNPR has SA-binding activity. NPRs and NPR-associated pathways have evolved distinctively in diverged land plant lineages to cope with different terrestrial environments.

Genome-wide analysis and identification of the TBL gene family in Eucalyptus grandis

The TRICHOME BIREFRINGENCE-LIKE (TBL) gene encodes a class of proteins related to xylan acetylation, which has been shown to play an important role in plant response to environmental stresses. This gene family has been meticulously investigated in Arabidopsis thaliana, whereas there have been no related reports in Eucalyptus grandis. In this study, we identified 49 TBL genes in E. grandis. A conserved amino acid motif was identified, which plays an important role in the execution of the function of TBL gene family members. The expression of TBL genes was generally upregulated in jasmonic acid-treated experiments, whereas it has been found that jasmonic acid activates the expression of genes involved in the defense functions of the plant body, suggesting that TBL genes play an important function in the response of the plant to stress. The principle of the action of TBL genes is supported by the finding that the xylan acetylation process increases the rigidity of the cell wall of the plant body and thus improves the plant's resistance to stress. The results of this study provide new information about the TBL gene family in E. grandis and will help in the study of the evolution, inheritance, and function of TBL genes in E. grandis, while confirming their functions.

Voice from both sides: a molecular dialogue between transcriptional activators and repressors in seed-to-seedling transition and crop adaptation

High-quality seeds provide valuable nutrients to human society and ensure successful seedling establishment. During maturation, seeds accumulate storage compounds that are required to sustain seedling growth during germination. This review focuses on the epigenetic repression of the embryonic and seed maturation programs in seedlings. We begin with an extensive overview of mutants affecting these processes, illustrating the roles of core proteins and accessory components in the epigenetic machinery by comparing mutants at both phenotypic and molecular levels. We highlight how omics assays help uncover target-specific functional specialization and coordination among various epigenetic mechanisms. Furthermore, we provide an in-depth discussion on the Seed dormancy 4 (Sdr4) transcriptional corepressor family, comparing and contrasting their regulation of seed germination in the dicotyledonous species Arabidopsis and two monocotyledonous crops, rice and wheat. Finally, we compare the similarities in the activation and repression of the embryonic and seed maturation programs through a shared set of cis-regulatory elements and discuss the challenges in applying knowledge largely gained in model species to crops.

Plant Kinesin Repertoires Expand with New Domain Architecture and Contract with the Loss of Flagella

Kinesins are eukaryotic microtubule motor proteins subdivided into conserved families with distinct functional roles. While many kinesin families are widespread in eukaryotes, each organismal lineage maintains a unique kinesin repertoire composed of many families with distinct numbers of genes. Previous genomic surveys indicated that land plant kinesin repertoires differ markedly from other eukaryotes. To determine when repertoires diverged during plant evolution, we performed robust phylogenomic analyses of kinesins in 24 representative plants, two algae, two animals, and one yeast. These analyses show that kinesin repertoires expand and contract coincident with major shifts in the biology of algae and land plants. One kinesin family and five subfamilies, each defined by unique domain architectures, emerged in the green algae. Four of those kinesin groups expanded in ancestors of modern land plants, while six other kinesin groups were lost in the ancestors of pollen-bearing plants. Expansions of different kinesin families and subfamilies occurred in moss and angiosperm lineages. Other kinesin families remained stable and did not expand throughout plant evolution. Collectively these data support a radiation of kinesin domain architectures in algae followed by differential positive and negative selection on kinesins families and subfamilies in different lineages of land plants.

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.

DNA methylation enables recurrent endogenization of giant viruses in an animal relative

5-Methylcytosine (5mC) is a widespread silencing mechanism that controls genomic parasites. In eukaryotes, 5mC has gained complex roles in gene regulation beyond parasite control, yet 5mC has also been lost in many lineages. The causes for 5mC retention and its genomic consequences are still poorly understood. Here, we show that the protist closely related to animals Amoebidium appalachense features both transposon and gene body methylation, a pattern reminiscent of invertebrates and plants. Unexpectedly, hypermethylated genomic regions in Amoebidium derive from viral insertions, including hundreds of endogenized giant viruses, contributing 14% of the proteome. Using a combination of inhibitors and genomic assays, we demonstrate that 5mC silences these giant virus insertions. Moreover, alternative Amoebidium isolates show polymorphic giant virus insertions, highlighting a dynamic process of infection, endogenization, and purging. Our results indicate that 5mC is critical for the controlled coexistence of newly acquired viral DNA into eukaryotic genomes, making Amoebidium a unique model to understand the hybrid origins of eukaryotic DNA.

Crossroads of assembling a moss genome: navigating contaminants and horizontal gene transfer in the moss Physcomitrellopsis africana

The first chromosome-scale reference genome of the rare narrow-endemic African moss Physcomitrellopsis africana (P. africana) is presented here. Assembled from 73 × Oxford Nanopore Technologies (ONT) long reads and 163 × Beijing Genomics Institute (BGI)-seq short reads, the 414 Mb reference comprises 26 chromosomes and 22,925 protein-coding genes [Benchmarking Universal Single-Copy Ortholog (BUSCO) scores: C:94.8% (D:13.9%)]. This genome holds 2 genes that withstood rigorous filtration of microbial contaminants, have no homolog in other land plants, and are thus interpreted as resulting from 2 unique horizontal gene transfers (HGTs) from microbes. Further, P. africana shares 176 of the 273 published HGT candidates identified in Physcomitrium patens (P. patens), but lacks 98 of these, highlighting that perhaps as many as 91 genes were acquired in P. patens in the last 40 million years following its divergence from its common ancestor with P. africana. These observations suggest rather continuous gene gains via HGT followed by potential losses during the diversification of the Funariaceae. Our findings showcase both dynamic flux in plant HGTs over evolutionarily "short" timescales, alongside enduring impacts of successful integrations, like those still functionally maintained in extant P. africana. Furthermore, this study describes the informatic processes employed to distinguish contaminants from candidate HGT events.

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.

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.

Evolutionary assembly of the plant terrestrialization toolkit from protein domains

Land plants (embryophytes) came about in a momentous evolutionary singularity: plant terrestrialization. This event marks not only the conquest of land by plants but also the massive radiation of embryophytes into a diverse array of novel forms and functions. The unique suite of traits present in the earliest land plants is thought to have been ushered in by a burst in genomic novelty. Here, we asked the question of how these bursts were possible. For this, we explored: (i) the initial emergence and (ii) the reshuffling of domains to give rise to hallmark environmental response genes of land plants. We pinpoint that a quarter of the embryophytic genes for stress physiology are specific to the lineage, yet a significant portion of this novelty arises not de novo but from reshuffling and recombining of pre-existing domains. Our data suggest that novel combinations of old genomic substrate shaped the plant terrestrialization toolkit, including hallmark processes in signalling, biotic interactions and specialized metabolism.

Conserved functions of chromatin regulators in basal Archaeplastida

Chromatin is a dynamic network that regulates genome organization and gene expression. Different types of chromatin regulators are highly conserved among Archaeplastida, including unicellular algae, while some chromatin genes are only present in land plant genomes. Here, we review recent advances in understanding the function of conserved chromatin factors in basal land plants and algae. We focus on the role of Polycomb-group genes which mediate H3K27me3-based silencing and play a role in balancing gene dosage and regulating haploid-to-diploid transitions by tissue-specific repression of the transcription factors KNOX and BELL in many representatives of the green lineage. Moreover, H3K27me3 predominantly occupies repetitive elements which can lead to their silencing in a unicellular alga and basal land plants, while it covers mostly protein-coding genes in higher land plants. In addition, we discuss the role of nuclear matrix constituent proteins as putative functional lamin analogs that are highly conserved among land plants and might have an ancestral function in stress response regulation. In summary, our review highlights the importance of studying chromatin regulation in a wide range of organisms in the Archaeplastida.

Designing a synthetic moss genome using GenoDesigner

The de novo synthesis of genomes has made unprecedented progress and achieved milestones, particularly in bacteria and yeast. However, the process of synthesizing a multicellular plant genome has not progressed at the same pace, due to the complexity of multicellular plant genomes, technical difficulties associated with large genome size and structure, and the intricacies of gene regulation and expression in plants. Here we outline the bottom-up design principles for the de novo synthesis of the Physcomitrium patens (that is, earthmoss) genome. To facilitate international collaboration and accessibility, we have developed and launched a public online design platform called GenoDesigner. This platform offers an intuitive graphical interface enabling users to efficiently manipulate extensive genome sequences, even up to the gigabase level. This tool is poised to greatly expedite the synthesis of the P. patens genome, offering an essential reference and roadmap for the synthesis of plant genomes.

Brassinosteroid-dependent phosphorylation of PHOSPHATE STARVATION RESPONSE2 reduces its DNA-binding ability in rice

Brassinosteroids (BRs) are widely used as plant growth regulators in modern agriculture. Understanding how BRs regulate nutrient signaling is crucial for reducing fertilizer usage. Here we elucidate that the central BR signaling inhibitor GSK3/SHAGGY-LIKE KINASE2 (GSK2) interacts directly with and phosphorylates PHOSPHATE STARVATION RESPONSE2 (OsPHR2), the key regulator of phosphate (Pi) signaling, to suppress its transcription factor activity in rice (Oryza sativa). We identify a critical phosphorylation site at serine residue S269 of OsPHR2 and demonstrate that phosphorylation by GSK2 or phosphor-mimic mutation of S269 substantially impairs the DNA-binding activity of OsPHR2, and thus diminishes expression of OsPHR2-induced genes and reduces Pi levels. Like BRs, Pi starvation noticeably induces GSK2 instability. We further show that this site-specific phosphorylation event is conserved in Arabidopsis (Arabidopsis thaliana), but varies among the PHR-family members, being present only in most land plants. These results unveil a distinctive post-transcriptional regulatory mechanism in Pi signaling by which BRs promote Pi acquisition, with a potential contribution to the environmental adaptability of plants during their evolution.

Identification, molecular evolution, codon bias, and expansion analysis of NLP transcription factor family in foxtail millet (Setaria italica L.) and closely related crops

The NODULE-INCEPTION-like protein (NLP) family is a plant-specific transcription factor (TF) family involved in nitrate transport and assimilation in plants, which are essential for improving plant nitrogen use efficiency. Currently, the molecular nature and evolutionary trajectory of NLP genes in the C4 model crop foxtail millet are unknown. Therefore, we performed a comprehensive analysis of NLP and molecular evolution in foxtail millet by scanning the genomes of foxtail millet and representative species of the plant kingdom. We identified seven NLP genes in the foxtail millet genome, all of which are individually and separately distributed on different chromosomes. They were not structurally identical to each other and were mainly expressed on root tissues. We unearthed two key genes (Si5G004100.1 and Si6G248300.1) with a variety of excellent characteristics. Regarding its molecular evolution, we found that NLP genes in Gramineae mainly underwent dispersed duplication, but maize NLP genes were mainly generated via WGD events. Other factors such as base mutations and natural selection have combined to promote the evolution of NLP genes. Intriguingly, the family in plants showed a gradual expansion during evolution with more duplications than losses, contrary to most gene families. In conclusion, this study advances the use of NLP genetic resources and the understanding of molecular evolution in cereals.

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.

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.

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.

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.