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

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.

The Identification of SQS/SQE/OSC Gene Families in Regulating the Biosynthesis of Triterpenes in Potentilla anserina

The tuberous roots of Potentilla anserina (Pan) are an edible and medicinal resource in Qinghai-Tibetan Plateau, China. The triterpenoids from tuberous roots have shown promising anti-cancer, hepatoprotective, and anti-inflammatory properties. In this study, we carried out phylogenetic analysis of squalene synthases (SQSs), squalene epoxidases (SQEs), and oxidosqualene cyclases (OSCs) in the pathway of triterpenes. In total, 6, 26, and 20 genes of SQSs, SQEs, and OSCs were retrieved from the genome of Pan, respectively. Moreover, 6 SQSs and 25 SQEs genes expressed in two sub-genomes (A and B) of Pan. SQSs were not expanded after whole-genome duplication (WGD), and the duplicated genes were detected in SQEs. Twenty OSCs were divided into two clades of cycloartenol synthases (CASs) and β-amyrin synthases (β-ASs) by a phylogenetic tree, characterized with gene duplication and evolutionary divergence. We speculated that β-ASs and CASs may participate in triterpenes synthesis. The data presented act as valuable references for future studies on the triterpene synthetic pathway of Pan.

Genome-wide analysis of bromodomain gene family in Arabidopsis and rice

The bromodomain-containing proteins (BRD-proteins) belongs to family of 'epigenetic mark readers', integral to epigenetic regulation. The BRD-members contain a conserved 'bromodomain' (BRD/BRD-fold: interacts with acetylated-lysine in histones), and several additional domains, making them structurally/functionally diverse. Like animals, plants also contain multiple Brd-homologs, however the extent of their diversity and impact of molecular events (genomic duplications, alternative splicing, AS) therein, is relatively less explored. The present genome-wide analysis of Brd-gene families of Arabidopsis thaliana and Oryza sativa showed extensive diversity in structure of genes/proteins, regulatory elements, expression pattern, domains/motifs, and the bromodomain (w.r.t. length, sequence, location) among the Brd-members. Orthology analysis identified thirteen ortholog groups (OGs), three paralog groups (PGs) and four singleton members (STs). While more than 40% Brd-genes were affected by genomic duplication events in both plants, AS-events affected 60% A. thaliana and 41% O. sativa genes. These molecular events affected various regions (promoters, untranslated regions, exons) of different Brd-members with potential impact on expression and/or structure-function characteristics. RNA-Seq data analysis indicated differences in tissue-specificity and stress response of Brd-members. Analysis by RT-qPCR revealed differential abundance and salt stress response of duplicate A. thaliana and O. sativa Brd-genes. Further analysis of AtBrd gene, AtBrdPG1b showed salinity-induced modulation of splicing pattern. Bromodomain (BRD)-region based phylogenetic analysis placed the A. thaliana and O. sativa homologs into clusters/sub-clusters, mostly consistent with ortholog/paralog groups. The bromodomain-region displayed several conserved signatures in key BRD-fold elements (α-helices, loops), along with variations (1-20 sites) and indels among the BRD-duplicates. Homology modeling and superposition identified structural variations in BRD-folds of divergent and duplicate BRD-members, which might affect their interaction with the chromatin histones, and associated functions. The study also showed contribution of various duplication events in Brd-gene family expansion among diverse plants, including several monocot and dicot plant species.

Genome-Wide Characterization and Expression Profiling of ABA Biosynthesis Genes in a Desert Moss Syntrichia caninervis

Syntrichia caninervis can survive under 80-90% protoplasmic water losses, and it is a model plant in desiccation tolerance research. A previous study has revealed that S. caninervis would accumulate ABA under dehydration stress, while the ABA biosynthesis genes in S. caninervis are still unknown. This study identified one ScABA1, two ScABA4s, five ScNCEDs, twenty-nine ScABA2s, one ScABA3, and four ScAAOs genes, indicating that the ABA biosynthesis genes were complete in S. caninervis. Gene location analysis showed that the ABA biosynthesis genes were evenly distributed in chromosomes but were not allocated to sex chromosomes. Collinear analysis revealed that ScABA1, ScNCED, and ScABA2 had homologous genes in Physcomitrella patens. RT-qPCR detection found that all of the ABA biosynthesis genes responded to abiotic stress; it further indicated that ABA plays an important role in S. caninervis. Moreover, the ABA biosynthesis genes in 19 representative plants were compared to study their phylogenetic and conserved motifs; the results suggested that the ABA biosynthesis genes were closely associated with plant taxa, but these genes had the same conserved domain in each plant. In contrast, there is a huge variation in the exon number between different plant taxa; it revealed that ABA biosynthesis gene structures are closely related to plant taxa. Above all, this study provides strong evidence demonstrating that ABA biosynthesis genes were conserved in the plant kingdom and deepens our understanding of the evolution of the phytohormone ABA.

A Physcomitrella PIN protein acts in spermatogenesis and sporophyte retention

The auxin efflux PIN-FORMED (PIN) proteins are conserved in all land plants and important players in plant development. In the moss Physcomitrella (Physcomitrium patens), three canonical PINs (PpPINA-C) are expressed in the leafy shoot (gametophore). PpPINA and PpPINB show functional activity in vegetative growth and sporophyte development. Here, we examined the role of PpPINC in the life cycle of Physcomitrella. We established reporter and knockout lines for PpPINC and analysed vegetative and reproductive tissues using microscopy and transcriptomic sequencing of moss gametangia. PpPINC is expressed in immature leaves, mature gametangia and during sporophyte development. The sperm cells (spermatozoids) of pinC knockout mutants exhibit increased motility and an altered flagella phenotype. Furthermore, the pinC mutants have a higher portion of differentially expressed genes related to spermatogenesis, increased fertility and an increased abortion rate of premeiotic sporophytes. Here, we show that PpPINC is important for spermatogenesis and sporophyte retention. We propose an evolutionary conserved way of polar growth during early moss embryo development and sporophyte attachment to the gametophore while suggesting the mechanical function in sporophyte retention of a ring structure, the Lorch ring.

An improved assembly of the "Cascade" hop (Humulus lupulus) genome uncovers signatures of molecular evolution and refines time of divergence estimates for the Cannabaceae family

We present a chromosome-level assembly of the Cascade hop (Humulus lupulus L. var. lupulus) genome. The hop genome is large (2.8 Gb) and complex, and early attempts at assembly were fragmented. Recent advances have made assembly of the hop genome more tractable, transforming the extent of investigation that can occur. The chromosome-level assembly of Cascade was developed by scaffolding the previously reported Cascade assembly generated with PacBio long-read sequencing and polishing with Illumina short-read DNA sequencing. We developed gene models and repeat annotations and used a controlled bi-parental mapping population to identify significant sex-associated markers. We assessed molecular evolution in gene sequences, gene family expansion and contraction, and time of divergence from Cannabis sativa and other closely related plant species using Bayesian inference. We identified the putative sex chromosome in the female genome based on significant sex-associated markers from the bi-parental mapping population. While the estimate of repeat content (~64%) is similar to the estimate for the hemp genome, syntenic blocks in hop contain a greater percentage of LTRs. Hop is enriched for disease resistance-associated genes in syntenic gene blocks and expanded gene families. The Cascade chromosome-level assembly will inform cultivation strategies and serve to deepen our understanding of the hop genomic landscape, benefiting hop researchers and the Cannabaceae genomics community.

Chromosome-level genomes of multicellular algal sisters to land plants illuminate signaling network evolution

The filamentous and unicellular algae of the class Zygnematophyceae are the closest algal relatives of land plants. Inferring the properties of the last common ancestor shared by these algae and land plants allows us to identify decisive traits that enabled the conquest of land by plants. We sequenced four genomes of filamentous Zygnematophyceae (three strains of Zygnema circumcarinatum and one strain of Z. cylindricum) and generated chromosome-scale assemblies for all strains of the emerging model system Z. circumcarinatum. Comparative genomic analyses reveal expanded genes for signaling cascades, environmental response, and intracellular trafficking that we associate with multicellularity. Gene family analyses suggest that Zygnematophyceae share all the major enzymes with land plants for cell wall polysaccharide synthesis, degradation, and modifications; most of the enzymes for cell wall innovations, especially for polysaccharide backbone synthesis, were gained more than 700 million years ago. In Zygnematophyceae, these enzyme families expanded, forming co-expressed modules. Transcriptomic profiling of over 19 growth conditions combined with co-expression network analyses uncover cohorts of genes 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.

Newly identified sex chromosomes in the Sphagnum (peat moss) genome alter carbon sequestration and ecosystem dynamics

Peatlands are crucial sinks for atmospheric carbon but are critically threatened due to warming climates. Sphagnum (peat moss) species are keystone members of peatland communities where they actively engineer hyperacidic conditions, which improves their competitive advantage and accelerates ecosystem-level carbon sequestration. To dissect the molecular and physiological sources of this unique biology, we generated chromosome-scale genomes of two Sphagnum species: S. divinum and S. angustifolium. Sphagnum genomes show no gene colinearity with any other reference genome to date, demonstrating that Sphagnum represents an unsampled lineage of land plant evolution. The genomes also revealed an average recombination rate an order of magnitude higher than vascular land plants and short putative U/V sex chromosomes. These newly described sex chromosomes interact with autosomal loci that significantly impact growth across diverse pH conditions. This discovery demonstrates that the ability of Sphagnum to sequester carbon in acidic peat bogs is mediated by interactions between sex, autosomes and environment.

Crossroads in the evolution of plant specialized metabolism

The monophyletic group of embryophytes (land plants) stands out among photosynthetic eukaryotes: they are the sole constituents of the macroscopic flora on land. In their entirety, embryophytes account for the majority of the biomass on land and constitute an astounding biodiversity. What allowed for the massive radiation of this particular lineage? One of the defining features of all land plants is the production of an array of specialized metabolites. The compounds that the specialized metabolic pathways of embryophytes produce have diverse functions, ranging from superabundant structural polymers and compounds that ward off abiotic and biotic challenges, to signaling molecules whose abundance is measured at the nanomolar scale. These specialized metabolites govern the growth, development, and physiology of land plants-including their response to the environment. Hence, specialized metabolites define the biology of land plants as we know it. And they were likely a foundation for their success. It is thus intriguing to find that the closest algal relatives of land plants, freshwater organisms from the grade of streptophyte algae, possess homologs for key enzymes of specialized metabolic pathways known from land plants. Indeed, some studies suggest that signature metabolites emerging from these pathways can be found in streptophyte algae. Here we synthesize the current understanding of which routes of the specialized metabolism of embryophytes can be traced to a time before plants had conquered land.

RALF peptides modulate immune response in the moss Physcomitrium patens

BACKGROUND

RAPID ALKALINIZATION FACTOR (RALFs) are cysteine-rich peptides that regulate multiple physiological processes in plants. This peptide family has considerably expanded during land plant evolution, but the role of ancient RALFs in modulating stress responses is unknown.Results: Here, we used the moss Physcomitrium patens as a model to gain insight into the role of RALF peptides in the coordination of plant growth and stress response in non-vascular plants. The quantitative proteomic analysis revealed concerted downregulation of M6 metalloprotease and some membrane proteins, including those involved in stress response, in PpRALF1, 2 and 3 knockout (KO) lines. The subsequent analysis revealed the role of PpRALF3 in growth regulation under abiotic and biotic stress conditions, implying the importance of RALFs in responding to various adverse conditions in bryophytes. We found that knockout of the PpRALF2 and PpRALF3 genes resulted in increased resistance to bacterial and fungal phytopathogens, Pectobacterium carotovorum and Fusarium solani, suggesting the role of these peptides in negative regulation of the immune response in P. patens. Comparing the transcriptomes of PpRALF3 KO and wild-type plants infected by F. solani showed that the regulation of genes in the phenylpropanoid pathway and those involved in cell wall modification and biogenesis was different in these two genotypes.

CONCLUSION

Thus, our study sheds light on the function of the previously uncharacterized PpRALF3 peptide and gives a clue to the ancestral functions of RALF peptides in plant stress response.

A deeply conserved protease, acylamino acid-releasing enzyme (AARE), acts in ageing in Physcomitrella and Arabidopsis

Reactive oxygen species (ROS) are constant by-products of aerobic life. In excess, ROS lead to cytotoxic protein aggregates, which are a hallmark of ageing in animals and linked to age-related pathologies in humans. Acylamino acid-releasing enzymes (AARE) are bifunctional serine proteases, acting on oxidized proteins. AARE are found in all domains of life, albeit under different names, such as acylpeptide hydrolase (APEH/ACPH), acylaminoacyl peptidase (AAP), or oxidized protein hydrolase (OPH). In humans, AARE malfunction is associated with age-related pathologies, while their function in plants is less clear. Here, we provide a detailed analysis of AARE genes in the plant lineage and an in-depth analysis of AARE localization and function in the moss Physcomitrella and the angiosperm Arabidopsis. AARE loss-of-function mutants have not been described for any organism so far. We generated and analysed such mutants and describe a connection between AARE function, aggregation of oxidized proteins and plant ageing, including accelerated developmental progression and reduced life span. Our findings complement similar findings in animals and humans, and suggest a unified concept of ageing may exist in different life forms.

WOX11: the founder of plant organ regeneration

De novo organ regeneration is the process in which adventitious roots or shoots regenerate from detached or wounded organs. De novo organ regeneration can occur either in natural conditions, e.g. adventitious root regeneration from the wounded sites of detached leaves or stems, or in in-vitro tissue culture, e.g. organ regeneration from callus. In this review, we summarize recent advances in research on the molecular mechanism of de novo organ regeneration, focusing on the role of the WUSCHEL-RELATED HOMEOBOX11 (WOX11) gene in the model plant Arabidopsis thaliana. WOX11 is a direct target of the auxin signaling pathway, and it is expressed in, and regulates the establishment of, the founder cell during de novo root regeneration and callus formation. WOX11 activates the expression of its target genes to initiate root and callus primordia. Therefore, WOX11 links upstream auxin signaling to downstream cell fate transition during regeneration. We also discuss the role of WOX11 in diverse species and its evolution in plants.

Genome-wide signatures of adaptation to extreme environments in red algae

The high temperature, acidity, and heavy metal-rich environments associated with hot springs have a major impact on biological processes in resident cells. One group of photosynthetic eukaryotes, the Cyanidiophyceae (Rhodophyta), has successfully thrived in hot springs and associated sites worldwide for more than 1 billion years. Here, we analyze chromosome-level assemblies from three representative Cyanidiophyceae species to study environmental adaptation at the genomic level. We find that subtelomeric gene duplication of functional genes and loss of canonical eukaryotic traits played a major role in environmental adaptation, in addition to horizontal gene transfer events. Shared responses to environmental stress exist in Cyanidiales and Galdieriales, however, most of the adaptive genes (e.g., for arsenic detoxification) evolved independently in these lineages. Our results underline the power of local selection to shape eukaryotic genomes that may face vastly different stresses in adjacent, extreme microhabitats.

PEATmoss: A Gene Expression Atlas for Bryophytes

PEATmoss is an interactive gene expression atlas for bryophytes, which originally unified Physcomitrium patens RNA-seq and microarray expression data from multiple gene annotation versions. This atlas includes more than 100 experiments of P. patens, is expanding to host Anthoceros agrestis and Marchantia polymorpha, and aims to host data from more species in the future. PEATmoss has multiple visualization methods and tools for data downloading and is connected to the Physcomitrium patens Gene Model Lookup DB (PpGML DB), which links P. patens genes to annotations and resources from several databases and contains tools for gene version lookup and sequence and annotation extraction. Among the new features available in PEATmoss are dataset privacy control, multispecies menu, interactive color scale, co-expression network visualization, and replicate data downloading.

Dirigent protein subfamily function and structure in terrestrial plant phenol metabolism

Aquatic plant transition to land, and subsequent terrestrial plant species diversification, was accompanied by the emergence and massive elaboration of plant phenol chemo-diversity. Concomitantly, dirigent protein (DP) and dirigent-like protein subfamilies, derived from large multigene families, emerged and became extensively diversified. DP biochemical functions as gateway entry points into new and diverse plant phenol skeletal types then markedly expanded. DPs have at least eight non-uniformly distributed subfamilies, with different DP subfamily members of known biochemical/physiological function now implicated as gateway entries to lignan, lignin, aromatic diterpenoid, pterocarpan and isoflavene pathways. While some other DP subfamily members have jacalin domains, both these and indeed the majority of DPs throughout the plant kingdom await discovery of their biochemical roles. Methods and approaches were developed to discover DP biochemical function as gateway entry points to distinct plant phenol skeletal types in land plants. Various DP 3D X-ray structural determinations enabled structure-based comparative sequence analysis and modeling to understand similarities and differences among the different DP subfamilies. We consider that the core DP β-barrel fold and associated characteristics are likely common to all DPs, with several residues conserved and nearly invariant. There is also considerable variation in residue composition and topography of the putative substrate binding pockets, as well as substantial differences in several loops, such as the β1-β2 loop. All DPs likely bind and stabilize quinone methide intermediates, while guiding distinctive regio- and/or stereo-chemical entry into Nature's chemo-diverse land plant phenol metabolic classes.

Phylotranscriptomics of liverworts: revisiting the backbone phylogeny and ancestral gene duplications

BACKGROUND AND AIMS

With some 7300 extant species, liverworts (Marchantiophyta) represent one of the major land plant lineages. The backbone relationships, such as the phylogenetic position of Ptilidiales, and the occurrence and timing of whole-genome duplications, are still contentious.

METHODS

Based on analyses of the newly generated transcriptome data for 38 liverworts and complemented with those publicly available, we reconstructed the evolutionary history of liverworts and inferred gene duplication events along the 55 taxon liverwort species tree.

KEY RESULTS

Our phylogenomic study provided an ordinal-level liverwort nuclear phylogeny and identified extensive gene tree conflicts and cyto-nuclear incongruences. Gene duplication analyses based on integrated phylogenomics and Ks distributions indicated no evidence of whole-genome duplication events along the backbone phylogeny of liverworts.

CONCLUSIONS

With a broadened sampling of liverwort transcriptomes, we re-evaluated the backbone phylogeny of liverworts, and provided evidence for ancient hybridizations followed by incomplete lineage sorting that shaped the deep evolutionary history of liverworts. The lack of whole-genome duplication during the deep evolution of liverworts indicates that liverworts might represent one of the few major embryophyte lineages whose evolution was not driven by whole-genome duplications.

Autophagy modulates growth and development in the moss Physcomitrium patens

Physcomitrium patens apical growing protonemal cells have the singularity that they continue to undergo cell divisions as the plant develops. This feature provides a valuable tool to study autophagy in the context of a multicellular apical growing tissue coupled to development. Herein, we showed that the core autophagy machinery is present in the moss P. patens, and characterized the 2D and 3D growth and development of atg5 and atg7 loss-of-function mutants under optimal and nutrient-deprived conditions. Our results showed that 2D growth of the different morphological and functional protonemata apical growing cells, chloronema and caulonema, is differentially modulated by this process. These differences depend on the protonema cell type and position along the protonemal filament, and growth condition. As a global plant response, the absence of autophagy favors the spread of the colony through protonemata growth at the expense of a reduction of the 3D growth, such as the buds and gametophore development, and thus the adult gametophytic and reproductive phases. Altogether this study provides valuable information suggesting that autophagy has roles during apical growth with differential responses within the cell types of the same tissue and contributes to life cycle progression and thus the growth and development of the 2D and 3D tissues of P. patens.

A phylogenomically informed five-order system for the closest relatives of land plants

The evolution of streptophytes had a profound impact on life on Earth. They brought forth those photosynthetic eukaryotes that today dominate the macroscopic flora: the land plants (Embryophyta).¹ There is convincing evidence that the unicellular/filamentous Zygnematophyceae-and not the morphologically more elaborate Coleochaetophyceae or Charophyceae-are the closest algal relatives of land plants.^2-6 Despite the species richness (>4,000), wide distribution, and key evolutionary position of the zygnematophytes, their internal phylogeny remains largely unresolved.^7,8 There are also putative zygnematophytes with interesting body plan modifications (e.g., filamentous growth) whose phylogenetic affiliations remain unknown. Here, we studied a filamentous green alga (strain MZCH580) from an Austrian peat bog with central or parietal chloroplasts that lack discernible pyrenoids. It represents Mougeotiopsis calospora PALLA, an enigmatic alga that was described more than 120 years ago⁹ but never subjected to molecular analyses. We generated transcriptomic data of M. calospora strain MZCH580 and conducted comprehensive phylogenomic analyses (326 nuclear loci) for 46 taxonomically diverse zygnematophytes. Strain MZCH580 falls in a deep-branching zygnematophycean clade together with some unicellular species and thus represents a formerly unknown zygnematophycean lineage with filamentous growth. Our well-supported phylogenomic tree lets us propose a new five-order system for the Zygnematophyceae and provides evidence for at least five independent origins of true filamentous growth in the closest algal relatives of land plants. This phylogeny provides a robust and comprehensive framework for performing comparative analyses and inferring the evolution of cellular traits and body plans in the closest relatives of land plants.

Endogenous giant viruses contribute to intraspecies genomic variability in the model green alga Chlamydomonas reinhardtii

Chlamydomonas reinhardtii is a unicellular eukaryotic alga that has been studied as a model organism for decades. Despite an extensive history as a model system, phylogenetic and genetic characteristics of viruses infecting this alga have remained elusive. We analyzed high-throughput genome sequence data of C. reinhardtii field isolates, and in six we discovered sequences belonging to endogenous giant viruses that reach up to several 100 kb in length. In addition, we have also discovered the entire genome of a closely related giant virus that is endogenized within the genome of Chlamydomonas incerta, the closest sequenced relative of C. reinhardtii. Endogenous giant viruses add hundreds of new gene families to the host strains, highlighting their contribution to the pangenome dynamics and interstrain genomic variability of C. reinhardtii. Our findings suggest that the endogenization of giant viruses may have important implications for structuring the population dynamics and ecology of protists in the environment.

From Plant to Yeast-Advances in Biosynthesis of Artemisinin

Malaria is a life-threatening disease. Artemisinin-based combination therapy (ACT) is the preferred choice for malaria treatment recommended by the World Health Organization. At present, the main source of artemisinin is extracted from Artemisia annua; however, the artemisinin content in A. annua is only 0.1-1%, which cannot meet global demand. Meanwhile, the chemical synthesis of artemisinin has disadvantages such as complicated steps, high cost and low yield. Therefore, the application of the synthetic biology approach to produce artemisinin in vivo has magnificent prospects. In this review, the biosynthesis pathway of artemisinin was summarized. Then we discussed the advances in the heterologous biosynthesis of artemisinin using microorganisms (Escherichia coli and Saccharomyces cerevisiae) as chassis cells. With yeast as the cell factory, the production of artemisinin was transferred from plant to yeast. Through the optimization of the fermentation process, the yield of artemisinic acid reached 25 g/L, thereby producing the semi-synthesis of artemisinin. Moreover, we reviewed the genetic engineering in A. annua to improve the artemisinin content, which included overexpressing artemisinin biosynthesis pathway genes, blocking key genes in competitive pathways, and regulating the expression of transcription factors related to artemisinin biosynthesis. Finally, the research progress of artemisinin production in other plants (Nicotiana, Physcomitrella, etc.) was discussed. The current advances in artemisinin biosynthesis may help lay the foundation for the remarkable up-regulation of artemisinin production in A. annua through gene editing or molecular design breeding in the future.

Characterization and evolutionary diversification of the phospholipase D gene family in mosses

Plant phospholipase D (PLD) exerts important roles in various biological processes, such as intracellular signaling and morphological development. Our knowledge about early land plant PLDs is still underdeveloped. In this study, we identified 84 PLD genes in six mosses, i.e., Physcomitrella patens, Ceratodon purpureus, Fontinalis antipyretica, Pleurozium schreberi, Sphagnum magellanicum, and Sphagnum fallax. These PLDs were classified into four clades (I–IV). We showed that PLD underwent rapid expansion in mosses. A total of six conserved domains and two core HKD motifs were detected. Structure analysis uncovered that the moss PLDs from within a clade generally exhibited similar exon-intron organization. Cis-elements prediction and expression analyses indicated that P. patens PLDs had key roles in stress responsiveness and plant development. Particularly, about half of the P. patens PLDs (e.g., PpPLD1, PpPLD2, and PpPLD5) were differentially expressed under biotic and abiotic stresses. We also determined the expression pattern of P. patens PLD genes in various tissues and at different stages of development. Although the moss, clubmoss, liverwort, and fern PLDs evolved largely under functional constraints, we found episodic positive selection in the moss PLDs, e.g., C. purpureus PLD2 and P. patens PLD11. We infer that the evolutionary force acting on the PLDs may have facilitated moss colonization of land. Our work provides valuable insights into the diversification of moss PLD genes, and can be used for future studies of their functions.

Genome-wide characterization of C2H2 zinc-finger gene family provides insight into the mechanisms and evolution of the dehydration-rehydration responses in Physcomitrium and Arabidopsis

Dehydration tolerance is a vital factor for land plant evolution and world agricultural production. Numerous studies enlightened that the plant-specific C2H2-type zinc-finger proteins (C2H2-ZFPs) as master regulators played pivotal roles in the abiotic stress responses of plants. However, a comprehensive understanding of the evolution of C2H2-ZFPs in terrestrial plants and its regulatory mechanism in dehydration and rehydration response remains a mystery. In this study, the genome-wide identification of C2H2-ZFP genes revealed 549 homologs in the representatives of terrestrial plant lineages from liverwort to angiosperms. Based on the characteristics of the conserved C2H2-ZF domains, four major C2H2-ZF types (M-, Z-, Q-, and D-type) were identified in the C2H2-ZFPs, with the dominants of M-type in all selected species and followed by Z-type in non-seed plants and Q-type in seed plants, respectively. Phylogenetic analyses of the identified C2H2-ZFPs supported four major groups in the land plant representatives, among which the members from the desiccation-tolerant Physcomitrium patens and the dehydration-sensitive Arabidopsis thaliana displayed different topological relationships in the phylogenies reconstructed for a single species. C2H2-ZFPs clustered in the same subclades shared similar features in their conserved domains and gene structures. Approximately, 81% of the C2H2-ZFP promoters of all 549 identified C2H2-ZFPs harbored the conserved ABA-responsive elements (ABREs) and/or dehydration-responsive elements (DREs). Comparative transcriptomic analyses showed that 50 PpZFPs and 56 AtZFPs significantly changed their transcripts abundance. Interestingly, most of the dehydration- and rehydration-responsive PpZPFs and AtZFPs had been predicted to contain the ABRE and DRE elements in their promoter regions and with over half of which phylogenetically belonging to group III. The differences in the expression patterns of C2H2-ZFPs in responses to dehydration and rehydration between P. patens and A. thaliana reflected their different strategies to adapt to dehydration. The identified candidate PpZFPs were specifically induced by moderate dehydration and reached the peak transcript abundance in severe dehydration. Our study lays the foundations for further functional investigation of C2H2-ZFPs in dehydration responses from an evolutionary perspective in land plants. The findings will provide us with genetic resources and potential targets for drought tolerance breeding in crops and beyond.

A seed-like proteome in oil-rich tubers

There are numerous examples of plant organs or developmental stages that are desiccation-tolerant and can withstand extended periods of severe water loss. One prime example are seeds and pollen of many spermatophytes. However, in some plants, also vegetative organs can be desiccation-tolerant. One example are the tubers of yellow nutsedge (Cyperus esculentus), which also store large amounts of lipids similar to seeds. Interestingly, the closest known relative, purple nutsedge (Cyperus rotundus), generates tubers that do not accumulate oil and are not desiccation-tolerant. We generated nanoLC-MS/MS-based proteomes of yellow nutsedge in five replicates of four stages of tuber development and compared them to the proteomes of roots and leaves, yielding 2257 distinct protein groups. Our data reveal a striking upregulation of hallmark proteins of seeds in the tubers. A deeper comparison to the tuber proteome of the close relative purple nutsedge (C. rotundus) and a previously published proteome of Arabidopsis seeds and seedlings indicates that indeed a seed-like proteome was found in yellow but not purple nutsedge. This was further supported by an analysis of the proteome of a lipid droplet-enriched fraction of yellow nutsedge, which also displayed seed-like characteristics. One reason for the differences between the two nutsedge species might be the expression of certain transcription factors homologous to ABSCISIC ACID INSENSITIVE3, WRINKLED1, and LEAFY COTYLEDON1 that drive gene expression in Arabidopsis seed embryos.

Expansion and contraction of small RNA and methylation machinery throughout plant evolution

The revolution in sequencing has created a wealth of plant genomes that can be mined to understand the evolution of biological complexity. Complexity is often driven by gene duplication, which allows paralogs to specialize in an activity of the ancestral gene or acquire novel functions. Angiosperms encode a variety of gene silencing pathways that share related machinery for small RNA biosynthesis and function. Recent phylogenetic analysis of these gene families plots the expansion, specialization, and occasional contraction of this core machinery. This analysis reveals the ancient origin of RNA-directed DNA Methylation in early land plants, or possibly their algal ancestors, as well as ongoing duplications that evolve novel small RNA pathways.

Core circadian clock and light signaling genes brought into genetic linkage across the green lineage

The circadian clock is conserved at both the level of transcriptional networks as well as core genes in plants, ensuring that biological processes are phased to the correct time of day. In the model plant Arabidopsis (Arabidopsis thaliana), the core circadian SHAQKYF-type-MYB (sMYB) genes CIRCADIAN CLOCK ASSOCIATED 1 (CCA1) and REVEILLE (RVE4) show genetic linkage with PSEUDO-RESPONSE REGULATOR 9 (PRR9) and PRR7, respectively. Leveraging chromosome-resolved plant genomes and syntenic ortholog analysis enabled tracing this genetic linkage back to Amborella trichopoda, a sister lineage to the angiosperm, and identifying an additional evolutionarily conserved genetic linkage in light signaling genes. The LHY/CCA1-PRR5/9, RVE4/8-PRR3/7, and PIF3-PHYA genetic linkages emerged in the bryophyte lineage and progressively moved within several genes of each other across an array of angiosperm families representing distinct whole-genome duplication and fractionation events. Soybean (Glycine max) maintained all but two genetic linkages, and expression analysis revealed the PIF3-PHYA linkage overlapping with the E4 maturity group locus was the only pair to robustly cycle with an evening phase, in contrast to the sMYB-PRR morning and midday phase. While most monocots maintain the genetic linkages, they have been lost in the economically important grasses (Poaceae), such as maize (Zea mays), where the genes have been fractionated to separate chromosomes and presence/absence variation results in the segregation of PRR7 paralogs across heterotic groups. The environmental robustness model is put forward, suggesting that evolutionarily conserved genetic linkages ensure superior microhabitat pollinator synchrony, while wide-hybrids or unlinking the genes, as seen in the grasses, result in heterosis, adaptation, and colonization of new ecological niches.

The renaissance and enlightenment of Marchantia as a model system

The liverwort Marchantia polymorpha has been utilized as a model for biological studies since the 18th century. In the past few decades, there has been a Renaissance in its utilization in genomic and genetic approaches to investigating physiological, developmental, and evolutionary aspects of land plant biology. The reasons for its adoption are similar to those of other genetic models, e.g. simple cultivation, ready access via its worldwide distribution, ease of crossing, facile genetics, and more recently, efficient transformation, genome editing, and genomic resources. The haploid gametophyte dominant life cycle of M. polymorpha is conducive to forward genetic approaches. The lack of ancient whole-genome duplications within liverworts facilitates reverse genetic approaches, and possibly related to this genomic stability, liverworts possess sex chromosomes that evolved in the ancestral liverwort. As a representative of one of the three bryophyte lineages, its phylogenetic position allows comparative approaches to provide insights into ancestral land plants. Given the karyotype and genome stability within liverworts, the resources developed for M. polymorpha have facilitated the development of related species as models for biological processes lacking in M. polymorpha.

The moss-specific transcription factor PpERF24 positively modulates immunity against fungal pathogens in Physcomitrium patens

APETALA2/ethylene response factors (AP2/ERFs) transcription factors (TFs) have greatly expanded in land plants compared to algae. In angiosperms, AP2/ERFs play important regulatory functions in plant defenses against pathogens and abiotic stress by controlling the expression of target genes. In the moss Physcomitrium patens, a high number of members of the ERF family are induced during pathogen infection, suggesting that they are important regulators in bryophyte immunity. In the current study, we analyzed a P. patens pathogen-inducible ERF family member designated as PpERF24. Orthologs of PpERF24 were only found in other mosses, while they were absent in the bryophytes Marchantia polymorpha and Anthoceros agrestis, the vascular plant Selaginella moellendorffii, and angiosperms. We show that PpERF24 belongs to a moss-specific clade with distinctive amino acids features in the AP2 domain that binds to the DNA. Interestingly, all P. patens members of the PpERF24 subclade are induced by fungal pathogens. The function of PpERF24 during plant immunity was assessed by an overexpression approach and transcriptomic analysis. Overexpressing lines showed increased defenses to infection by the fungal pathogens Botrytis cinerea and Colletotrichum gloeosporioides evidenced by reduced cellular damage and fungal biomass compared to wild-type plants. Transcriptomic and RT-qPCR analysis revealed that PpERF24 positively regulates the expression levels of defense genes involved in transcriptional regulation, phenylpropanoid and jasmonate pathways, oxidative burst and pathogenesis-related (PR) genes. These findings give novel insights into potential mechanism by which PpERF24 increases plant defenses against several pathogens by regulating important players in plant immunity.

Submergence of the filamentous Zygnematophyceae Mougeotia induces differential gene expression patterns associated with core metabolism and photosynthesis

The streptophyte algal class Zygnematophyceae is the closest algal sister lineage to land plants. In nature, Zygnematophyceae can grow in both terrestrial and freshwater habitats and how they do this is an important unanswered question. Here, we studied what happens to the zygnematophyceaen alga Mougeotia sp., which usually occurs in permanent and temporary freshwater bodies, when it is shifted to liquid growth conditions after growth on a solid substrate. Using global differential gene expression profiling, we identified changes in the core metabolism of the organism interlinked with photosynthesis; the latter went hand in hand with measurable impact on the photophysiology as assessed via pulse amplitude modulation (PAM) fluorometry. Our data reveal a pronounced change in the overall physiology of the alga after submergence and pinpoint candidate genes that play a role. These results provide insight into the importance of photophysiological readjustment when filamentous Zygnematophyceae transition between terrestrial and aquatic habitats.

Plant peptidoglycan precursor biosynthesis: Conservation between moss chloroplasts and Gram-negative bacteria

Accumulating evidence suggests that peptidoglycan, consistent with a bacterial cell wall, is synthesized around the chloroplasts of many photosynthetic eukaryotes, from glaucophyte algae to early-diverging land plants including pteridophyte ferns, but the biosynthetic pathway has not been demonstrated. Here, we employed mass spectrometry and enzymology in a two-fold approach to characterize the synthesis of peptidoglycan in chloroplasts of the moss Physcomitrium (Physcomitrella) patens. To drive the accumulation of peptidoglycan pathway intermediates, P. patens was cultured with the antibiotics fosfomycin, D-cycloserine, and carbenicillin, which inhibit key peptidoglycan pathway proteins in bacteria. Mass spectrometry of the trichloroacetic acid-extracted moss metabolome revealed elevated levels of five of the predicted intermediates from uridine diphosphate N-acetylglucosamine (UDP-GlcNAc) through the uridine diphosphate N-acetylmuramic acid (UDP-MurNAc)-D,L-diaminopimelate (DAP)-pentapeptide. Most Gram-negative bacteria, including cyanobacteria, incorporate meso-diaminopimelic acid (D,L-DAP) into the third residue of the stem peptide of peptidoglycan, as opposed to L-lysine, typical of most Gram-positive bacteria. To establish the specificity of D,L-DAP incorporation into the P. patens precursors, we analyzed the recombinant protein UDP-N-acetylmuramoyl-L-alanyl-D-glutamate-2,6-diaminopimelate ligase (MurE) from both P. patens and the cyanobacterium Anabaena sp. (Nostoc sp. strain PCC 7120). Both ligases incorporated D,L-DAP in almost complete preference to L-Lys, consistent with the mass spectrophotometric data, with catalytic efficiencies similar to previously documented Gram-negative bacterial MurE ligases. We discuss how these data accord with the conservation of active site residues common to DL-DAP-incorporating bacterial MurE ligases and of the probability of a horizontal gene transfer event within the plant peptidoglycan pathway.

Evolution of the DEHYDRATION-RESPONSIVE ELEMENT-BINDING PROTEIN subfamily in green plants

Adapting to unfavorable environments is a necessary step in plant terrestrialization and radiation. The dehydration-responsive element-binding (DREB) protein subfamily plays a pivotal role in plant abiotic stress regulation. However, relationships between the origin and expansion of the DREB subfamily and adaptive evolution of land plants are still being elucidated. Here, we constructed the evolutionary history of the DREB subfamily by compiling APETALA2/ethylene-responsive element-binding protein superfamily genes from 169 representative species of green plants. Through extensive phylogenetic analyses and comparative genomic analysis, our results revealed that the DREB subfamily diverged from the ethylene-responsive factor (ERF) subfamily in the common ancestor of Zygnemophyceae and Embryophyta during the colonization of land by plants, followed by expansions to form three different ancient archetypal genes in Zygnemophyceae species, designated as groups archetype-I, archetype-II/III, and archetype-IV. Four large-scale expansions paralleling the evolution of land plants led to the nine-subgroup divergence of group archetype-II/III in angiosperms, and five whole-genome duplications during Brassicaceae and Poaceae radiation shaped the diversity of subgroup IIb-1. We identified a Poaceae-specific gene in subgroup IIb-1, ERF014, remaining in a Poaceae-specific microsynteny block and co-evolving with a small heat shock protein cluster. Expression analyses demonstrated that heat acclimation may have driven the neofunctionalization of ERF014s in Pooideae by engaging in the conserved heat-responsive module in Poaceae. This study provides insights into lineage-specific expansion and neofunctionalization in the DREB subfamily, together with evolutionary information valuable for future functional studies of plant stress biology.

Viral Complexity

Although traditionally viewed as streamlined and simple, discoveries over the last century have revealed that viruses can exhibit surprisingly complex physical structures, genomic organization, ecological interactions, and evolutionary histories. Viruses can have physical dimensions and genome lengths that exceed many cellular lineages, and their infection strategies can involve a remarkable level of physiological remodeling of their host cells. Virus–virus communication and widespread forms of hyperparasitism have been shown to be common in the virosphere, demonstrating that dynamic ecological interactions often shape their success. And the evolutionary histories of viruses are often fraught with complexities, with chimeric genomes including genes derived from numerous distinct sources or evolved de novo. Here we will discuss many aspects of this viral complexity, with particular emphasis on large DNA viruses, and provide an outlook for future research.

Measuring the Invisible: The Sequences Causal of Genome Size Differences in Eyebrights (Euphrasia) Revealed by k-mers

Genome size variation within plant taxa is due to presence/absence variation, which may affect low-copy sequences or genomic repeats of various frequency classes. However, identifying the sequences underpinning genome size variation is challenging because genome assemblies commonly contain collapsed representations of repetitive sequences and because genome skimming studies by design miss low-copy number sequences. Here, we take a novel approach based on k-mers, short sub-sequences of equal length k, generated from whole-genome sequencing data of diploid eyebrights (Euphrasia), a group of plants that have considerable genome size variation within a ploidy level. We compare k-mer inventories within and between closely related species, and quantify the contribution of different copy number classes to genome size differences. We further match high-copy number k-mers to specific repeat types as retrieved from the RepeatExplorer2 pipeline. We find genome size differences of up to 230Mbp, equivalent to more than 20% genome size variation. The largest contributions to these differences come from rDNA sequences, a 145-nt genomic satellite and a repeat associated with an Angela transposable element. We also find size differences in the low-copy number class (copy number ≤ 10×) of up to 27 Mbp, possibly indicating differences in gene space between our samples. We demonstrate that it is possible to pinpoint the sequences causing genome size variation within species without the use of a reference genome. Such sequences can serve as targets for future cytogenetic studies. We also show that studies of genome size variation should go beyond repeats if they aim to characterise the full range of genomic variants. To allow future work with other taxonomic groups, we share our k-mer analysis pipeline, which is straightforward to run, relying largely on standard GNU command line tools.

An overview of bioinformatics, genomics, and transcriptomics resources for bryophytes

Bryophytes are useful models for the study of plant evolution, development, plant-fungal symbiosis, stress responses, and gametogenesis. Additionally, their dominant haploid gametophytic phase makes them great models for functional genomics research, allowing straightforward genome editing and gene knockout via CRISPR or homologous recombination. Until 2016, however, the only bryophyte genome sequence published was that of Physcomitrium patens. Throughout recent years, several other bryophyte genomes and transcriptome datasets became available, enabling better comparative genomics in evolutionary studies. The increase in the number of bryophyte genome and transcriptome resources available has yielded a plethora of annotations, databases, and bioinformatics tools to access the new data, which covers the large diversity of this clade and whose biology comprises features such as association with arbuscular mycorrhiza fungi, sex chromosomes, low gene redundancy, or loss of RNA editing genes for organellar transcripts. Here we provide a guide to resources available for bryophytes with regards to genome and transcriptome databases and bioinformatics tools.

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

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

Diversity, phylogeny, and adaptation of bryophytes: insights from genomic and transcriptomic data

Bryophytes including mosses, liverworts, and hornworts are among the earliest land plants, and occupy a crucial phylogenetic position to aid in the understanding of plant terrestrialization. Despite their small size and simple structure, bryophytes are the second largest group of extant land plants. They live ubiquitously in various habitats and are highly diversified, with adaptive strategies to modern ecosystems on Earth. More and more genomes and transcriptomes have been assembled to address fundamental questions in plant biology. Here, we review recent advances in bryophytes associated with diversity, phylogeny, and ecological adaptation. Phylogenomic studies have provided increasing supports for the monophyly of bryophytes, with hornworts sister to the Setaphyta clade including liverworts and mosses. Further comparative genomic analyses revealed that multiple whole-genome duplications might have contributed to the species richness and morphological diversity in mosses. We highlight that the biological changes through gene gain or neofunctionalization that primarily evolved in bryophytes have facilitated the adaptation to early land environments; among the strategies to adapt to modern ecosystems in bryophytes, desiccation tolerance is the most remarkable. More genomic information for bryophytes would shed light on key mechanisms for the ecological success of these 'dwarfs' in the plant kingdom.

Apical dominance control by TAR-YUC-mediated auxin biosynthesis is a deep homology of land plants

A key aim in biology is to identify which genetic changes contributed to the evolution of form through time. Apical dominance, the inhibitory effect exerted by shoot apices on the initiation or outgrowth of distant lateral buds, is a major regulatory mechanism of plant form.¹ Nearly a century of studies in the sporophyte of flowering plants have established the phytohormone auxin as a front-runner in the search for key factors controlling apical dominance,^2,3 identifying critical roles for long-range polar auxin transport and local auxin biosynthesis in modulating shoot branching.^4-10 A capacity for lateral branching evolved by convergence in the gametophytic shoot of mosses and primed its diversification;¹¹ however, polar auxin transport is relatively unimportant in this developmental process,¹² the contribution of auxin biosynthesis genes has not been assessed, and more generally, the extent of conservation in apical dominance regulation within the land plants remains largely unknown. To fill this knowledge gap, we sought to identify genetic determinants of apical dominance in the moss Physcomitrium patens. Here, we show that leafy shoot apex decapitation releases apical dominance through massive and rapid transcriptional reprogramming of auxin-responsive genes and altering auxin biosynthesis gene activity. We pinpoint a subset of P. patens TRYPTOPHAN AMINO-TRANSFERASE (TAR) and YUCCA FLAVIN MONOOXYGENASE-LIKE (YUC) auxin biosynthesis genes expressed in the main and lateral shoot apices and show that they are essential for coordinating branch initiation and outgrowth. Our results demonstrate that local auxin biosynthesis acts as a pivotal regulator of apical dominance in moss and constitutes a shared mechanism underpinning shoot architecture control in land plants.

Evolution of meristem zonation by CLE gene duplication in land plants

In angiosperms, a negative feedback pathway involving CLAVATA3 (CLV3) peptide and WUSCHEL transcription factor maintains the stem-cell population in the shoot apical meristem and is central for continued shoot growth and organogenesis. An intriguing question is how this cell-signalling system was established during the evolution of land plants. On the basis of two recent studies on CLV3/ESR-related (CLE) genes, this paper proposes a model for the evolution of meristem zonation. The model suggests that a stem-cell-limiting CLV3 pathway is derived from stem-cell-promoting CLE pathways conserved in land pants by gene duplication in the angiosperm lineage. The model can be examined in the future by genomic and developmental studies on diverse plant species.

The genomic basis of the plant island syndrome in Darwin's giant daisies

The repeated, rapid and often pronounced patterns of evolutionary divergence observed in insular plants, or the ‘plant island syndrome’, include changes in leaf phenotypes, growth, as well as the acquisition of a perennial lifestyle. Here, we sequence and describe the genome of the critically endangered, Galápagos-endemic species Scalesia atractyloides Arnot., obtaining a chromosome-resolved, 3.2-Gbp assembly containing 43,093 candidate gene models. Using a combination of fossil transposable elements, k-mer spectra analyses and orthologue assignment, we identify the two ancestral genomes, and date their divergence and the polyploidization event, concluding that the ancestor of all extant Scalesia species was an allotetraploid. There are a comparable number of genes and transposable elements across the two subgenomes, and while their synteny has been mostly conserved, we find multiple inversions that may have facilitated adaptation. We identify clear signatures of selection across genes associated with vascular development, growth, adaptation to salinity and flowering time, thus finding compelling evidence for a genomic basis of the island syndrome in one of Darwin’s giant daisies.

Many island plant species share a syndrome of characteristic phenotype and life history. Cerca et al. find the genomic basis of the plant island syndrome in one of Darwin’s giant daisies, while separating ancestral genomes in a chromosome-resolved polyploid assembly.