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.

The mosaic oat genome gives insights into a uniquely healthy cereal crop.

Cultivated oat (Avena sativa L.) is an allohexaploid (AACCDD, 2n = 6x = 42) thought to have been domesticated more than 3,000 years ago while growing as a weed in wheat, emmer and barley fields in Anatolia1,2. Oat has a low carbon footprint, substantial health benefits and the potential to replace animal-based food products. However, the lack of a fully annotated reference genome has hampered efforts to deconvolute its complex evolutionary history and functional gene dynamics. Here we present a high-quality reference genome of A. sativa and close relatives of its diploid (Avena longiglumis, AA, 2n = 14) and tetraploid (Avena insularis, CCDD, 2n = 4x = 28) progenitors. We reveal the mosaic structure of the oat genome, trace large-scale genomic reorganizations in the polyploidization history of oat and illustrate a breeding barrier associated with the genome architecture of oat. We showcase detailed analyses of gene families implicated in human health and nutrition, which adds to the evidence supporting oat safety in gluten-free diets, and we perform mapping-by-sequencing of an agronomic trait related to water-use efficiency. This resource for the Avena genus will help to leverage knowledge from other cereal genomes, improve understanding of basic oat biology and accelerate genomics-assisted breeding and reanalysis of quantitative trait studies.

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.

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.

Viral suppressor of RNA silencing in vascular plants also interferes with the development of the bryophyte Physcomitrella patens.

Plant viruses are important pathogens able to overcome plant defense mechanisms using their viral suppressors of RNA silencing (VSR). Small RNA pathways of bryophytes and vascular plants have significant similarities, but little is known about how viruses interact with mosses. This study elucidated the responses of Physcomitrella patens to two different VSRs. We transformed P. patens plants to express VSR P19 from tomato bushy stunt virus and VSR 2b from cucumber mosaic virus, respectively. RNA sequencing and quantitative PCR were used to detect the effects of VSRs on gene expression. Small RNA (sRNA) sequencing was used to estimate the influences of VSRs on the sRNA pool of P. patens. Expression of either VSR-encoding gene caused developmental disorders in P. patens. The transcripts of four different transcription factors (AP2/erf, EREB-11 and two MYBs) accumulated in the P19 lines. sRNA sequencing revealed that VSR P19 significantly changed the microRNA pool in P. patens. Our results suggest that VSR P19 is functional in P. patens and affects the abundance of specific microRNAs interfering with gene expression. The results open new opportunities for using Physcomitrella as an alternative system to study plant-virus interactions.

Axenic in vitro cultivation of 19 peat moss (Sphagnum L.) species as a resource for basic biology, biotechnology, and paludiculture.

Sphagnum farming can substitute peat with renewable biomass and thus help mitigate climate change. Large volumes of the required founder material can only be supplied sustainably by axenic cultivation in bioreactors. We established axenic in vitro cultures from sporophytes of 19 Sphagnum species collected in Austria, Germany, Latvia, the Netherlands, Russia, and Sweden: S. angustifolium, S. balticum, S. capillifolium, S. centrale, S. compactum, S. cuspidatum, S. fallax, S. fimbriatum, S. fuscum, S. lindbergii, S. medium/divinum, S. palustre, S. papillosum, S. rubellum, S. russowii, S. squarrosum, S. subnitens, S. subfulvum and S. warnstorfii. These species cover five of the six European Sphagnum subgenera; namely, Acutifolia, Cuspidata, Rigida, Sphagnum and Squarrosa. Their growth was measured in suspension cultures, whereas their ploidy was determined by flow cytometry and compared with the genome size of Physcomitrella patens. We identified haploid and diploid Sphagnum species, found that their cells are predominantly arrested in the G1 phase of the cell cycle, and did not find a correlation between plant productivity and ploidy. DNA barcoding was achieved by sequencing introns of the BRK1 genes. With this collection, high-quality founder material for diverse large-scale applications, but also for basic Sphagnum research, is available from the International Moss Stock Center.

RecQ Helicases Function in Development, DNA Repair, and Gene Targeting in Physcomitrella patens.

RecQ DNA helicases are genome surveillance proteins found in all kingdoms of life. They are characterized best in humans, as mutations in RecQ genes lead to developmental abnormalities and diseases. To better understand RecQ functions in plants we concentrated on Arabidopsis thaliana and Physcomitrella patens, the model species predominantly used for studies on DNA repair and gene targeting. Phylogenetic analysis of the six P. patens RecQ genes revealed their orthologs in humans and plants. Because Arabidopsis and P. patens differ in their RecQ4 and RecQ6 genes, reporter and deletion moss mutants were generated and gene functions studied in reciprocal cross-species and cross-kingdom approaches. Both proteins can be found in meristematic moss tissues, although at low levels and with distinct expression patterns. PpRecQ4 is involved in embryogenesis and in subsequent development as demonstrated by sterility of ΔPpRecQ4 mutants and by morphological aberrations. Additionally, ΔPpRecQ4 displays an increased sensitivity to DNA damages and an increased rate of gene targeting. Therefore, we conclude that PpRecQ4 acts as a repressor of recombination. In contrast, PpRecQ6 is not obviously important for moss development or DNA repair but does function as a potent enhancer of gene targeting.

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

The draft genome of the moss model, Physcomitrella patens, comprised approximately 2000 unordered scaffolds. In order to enable analyses of genome structure and evolution we generated a chromosome-scale genome assembly using genetic linkage as well as (end) sequencing of long DNA fragments. We find that 57% of the genome comprises transposable elements (TEs), some of which may be actively transposing during the life cycle. Unlike in flowering plant genomes, gene- and TE-rich regions show an overall even distribution along the chromosomes. However, the chromosomes are mono-centric with peaks of a class of Copia elements potentially coinciding with centromeres. Gene body methylation is evident in 5.7% of the protein-coding genes, typically coinciding with low GC and low expression. Some giant virus insertions are transcriptionally active and might protect gametes from viral infection via siRNA mediated silencing. Structure-based detection methods show that the genome evolved via two rounds of whole genome duplications (WGDs), apparently common in mosses but not in liverworts and hornworts. Several hundred genes are present in colinear regions conserved since the last common ancestor of plants. These syntenic regions are enriched for functions related to plant-specific cell growth and tissue organization. The P. patens genome lacks the TE-rich pericentromeric and gene-rich distal regions typical for most flowering plant genomes. More non-seed plant genomes are needed to unravel how plant genomes evolve, and to understand whether the P. patens genome structure is typical for mosses or bryophytes.

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.

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.

Expression of a human cDNA in moss results in spliced mRNAs and fragmentary protein isoforms.

Production of biopharmaceuticals relies on the expression of mammalian cDNAs in host organisms. Here we show that the expression of a human cDNA in the moss Physcomitrium patens generates the expected full-length and four additional transcripts due to unexpected splicing. This mRNA splicing results in non-functional protein isoforms, cellular misallocation of the proteins and low product yields. We integrated these results together with the results of our analysis of all 32,926 protein-encoding Physcomitrella genes and their 87,533 annotated transcripts in a web application, physCO, for automatized optimization. A thus optimized cDNA results in about twelve times more protein, which correctly localizes to the ER. An analysis of codon preferences of different production hosts suggests that similar effects occur also in non-plant hosts. We anticipate that the use of our methodology will prevent so far undetected mRNA heterosplicing resulting in maximized functional protein amounts for basic biology and biotechnology.

Convergence of sphingolipid desaturation across over 500 million years of plant evolution.

For plants, acclimation to low temperatures is fundamental to survival. This process involves the modification of lipids to maintain membrane fluidity. We previously identified a new cold-induced putative desaturase in Physcomitrium (Physcomitrella) patens. Lipid profiles of null mutants of this gene lack sphingolipids containing monounsaturated C24 fatty acids, classifying the new protein as sphingolipid fatty acid denaturase (PpSFD). PpSFD mutants showed a cold-sensitive phenotype as well as higher susceptibility to the oomycete Pythium, assigning functions in stress tolerance for PpSFD. Ectopic expression of PpSFD in the Atads2.1 (acyl coenzyme A desaturase-like 2) Arabidopsis thaliana mutant functionally complemented its cold-sensitive phenotype. While these two enzymes catalyse a similar reaction, their evolutionary origin is clearly different since AtADS2 is a methyl-end desaturase whereas PpSFD is a cytochrome b5 fusion desaturase. Altogether, we suggest that adjustment of membrane fluidity evolved independently in mosses and seed plants, which diverged more than 500 million years ago.

Selfing in Haploid Plants and Efficacy of Selection: Codon Usage Bias in the Model Moss Physcomitrella patens.

A long-term reduction in effective population size will lead to major shift in genome evolution. In particular, when effective population size is small, genetic drift becomes dominant over natural selection. The onset of self-fertilization is one evolutionary event considerably reducing effective size of populations. Theory predicts that this reduction should be more dramatic in organisms capable for haploid than for diploid selfing. Although theoretically well-grounded, this assertion received mixed experimental support. Here, we test this hypothesis by analyzing synonymous codon usage bias of genes in the model moss Physcomitrella patens frequently undergoing haploid selfing. In line with population genetic theory, we found that the effect of natural selection on synonymous codon usage bias is very weak. Our conclusion is supported by four independent lines of evidence: 1) Very weak or nonsignificant correlation between gene expression and codon usage bias, 2) no increased codon usage bias in more broadly expressed genes, 3) no evidence that codon usage bias would constrain synonymous and nonsynonymous divergence, and 4) predominant role of genetic drift on synonymous codon usage predicted by a model-based analysis. These findings show striking similarity to those observed in AT-rich genomes with weak selection for optimal codon usage and GC content overall. Our finding is in contrast to a previous study reporting adaptive codon usage bias in the moss P. patens.