Search

Gene functional descriptions, which are typically derived from sequence similarity to experimentally validated genes in a handful of model species, offer a crucial line of evidence when searching for candidate genes that underlie trait variation. Plant responses to environmental cues, including gene expression regulatory variation, represent important resources for understanding gene function and crucial targets for plant improvement through gene editing and other biotechnologies. However, even after years of effort and numerous large-scale functional characterization studies, biological roles of large proportions of protein coding genes across the plant phylogeny are poorly annotated. Here we describe the Joint Genome Institute (JGI) Plant Gene Atlas, a public and updateable data resource consisting of transcript abundance assays from 2,090 samples derived from 604 tissues or conditions across 18 diverse species. We integrated across these diverse conditions and genotypes by analyzing expression profiles, building gene clusters that exhibited tissue/condition specific expression, and testing for transcriptional modulation in response to environmental queues. For example, we discovered extensive phylogenetically constrained and condition-specific expression profiles across many gene families and genes without any functional annotation. Such conserved expression patterns and other tightly co-expressed gene clusters let us assign expression derived functional descriptions to 64,620 genes with otherwise unknown functions. The ever-expanding Gene Atlas resource is available at JGI Plant Gene Atlas (https://plantgeneatlas.jgi.doe.gov) and Phytozome (https://phytozome-next.jgi.doe.gov), providing bulk access to data and user-specified queries of gene sets. Combined, these web interfaces let users access differentially expressed genes, track orthologs across the Gene Atlas plants, graphically represent co-expressed genes, and visualize gene ontology and pathway enrichments.

Physcomitrella patens is a bryophyte model plant that is often used to study plant evolution and development. Its resources are of great importance for comparative genomics and evo-devo approaches. However, expression data from Physcomitrella patens were so far generated using different gene annotation versions and three different platforms: CombiMatrix and NimbleGen expression microarrays and RNA sequencing. The currently available P. patens expression data are distributed across three tools with different visualization methods to access the data. Here, we introduce an interactive expression atlas, Physcomitrella Expression Atlas Tool (PEATmoss), that unifies publicly available expression data for P. patens and provides multiple visualization methods to query the data in a single web-based tool. Moreover, PEATmoss includes 35 expression experiments not previously available in any other expression atlas. To facilitate gene expression queries across different gene annotation versions, and to access P. patens annotations and related resources, a lookup database and web tool linked to PEATmoss was implemented. PEATmoss can be accessed at https://peatmoss.online.uni-marburg.de

Nuclear-encoded pentatricopeptide repeat (PPR) proteins are key factors for site-specific RNA editing, converting cytidines into uridines in plant mitochondria and chloroplasts. All editing factors in the model moss Physcomitrella patens have a C-terminal DYW domain with similarity to cytidine deaminase. However, numerous editing factors in flowering plants lack such a terminal DYW domain, questioning its immediate role in the pyrimidine base conversion process. Here we further investigate the Physcomitrella DYW-type PPR protein PPR_78, responsible for mitochondrial editing sites cox1eU755SL and rps14eU137SL. Complementation assays with truncated proteins demonstrate that the DYW domain is essential for full PPR_78 editing functionality. The DYW domain can be replaced, however, with its counterpart from another editing factor, PPR_79. The PPR_78 ortholog of the related moss Funaria hygrometrica fully complements the Physcomitrella mutant for editing at both sites, although the editing site in rps14 is lacking in Funaria. Editing factor orthologs in different taxa may thus retain editing capacity for multiple sites despite the absence of editing requirement.

Across the plant kingdom, phytochrome (PHY) photoreceptors play an important role during adaptive and developmental responses to light. In Arabidopsis thaliana, light-activated PHYs accumulate in the nucleus, where they regulate downstream signaling components, such as phytochrome interacting factors (PIFs). PIFs are transcription factors that act as repressors of photomorphogenesis; their inhibition by PHYs leads to substantial changes in gene expression. The nuclear function of PHYs, however, has so far been investigated in only a few non-seed plants. Here, we identified putative target genes of PHY signaling in the moss Physcomitrella patens and found light-regulated genes that are putative orthologs of PIF-controlled genes in Arabidopsis. Phylogenetic analyses revealed that an ancestral PIF-like gene was already present in streptophyte algae, i.e., before the water-to-land transition of plants. The PIF homologs in the genome of P. patens resemble Arabidopsis PIFs in their protein domain structure, molecular properties, and physiological effects, albeit with notable differences in the motif-dependent PHY interaction. Our results suggest that P. patens PIFs are involved in PHY signaling. The PHY-PIF signaling node that relays light signals to target genes has been largely conserved during land plant evolution, with evidence of lineage-specific diversification.

The phytohormone abscisic acid (ABA) is a pivotal regulator of gene expression in response to various environmental stresses such as desiccation, salt and cold causing major changes in plant development and physiology. Here we show that in the moss Physcomitrella patens exogenous application of ABA triggers the formation of vegetative diaspores (brachycytes or brood cells) that enable plant survival in unfavorable environmental conditions. Such diaspores are round-shaped cells characterized by the loss of the central vacuole, due to an increased starch and lipid storage preparing these cells for growth upon suitable environmental conditions. To gain insights into the gene regulation underlying these developmental and physiological changes, we analyzed early transcriptome changes after 30, 60, and 180 min of ABA application and identified 1,030 differentially expressed genes. Among these, several groups can be linked to specific morphological and physiological changes during diaspore formation, such as genes involved in cell wall modifications. Furthermore, almost all members of ABA-dependent signaling and regulation were transcriptionally induced. Network analysis of transcription-associated genes revealed a large overlap of our study with ABA-dependent regulation in response to dehydration, cold stress, and UV-B light, indicating a fundamental function of ABA in diverse stress responses in moss. We also studied the evolutionary conservation of ABA-dependent regulation between moss and the seed plant Arabidopsis thaliana pointing to an early evolution of ABA-mediated stress adaptation during the conquest of the terrestrial habitat by plants.

et al., 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., The work conducted by the US Department of Energy Joint Genome Institute is supported by the Office of Science of the US Department of Energy under Contract No. DE-AC02-05CH11231. Support to RR and SAR by the German Research Foundation (DFG RE 837/10-2), the Excellence Initiative of the German Federal and State Governments (EXC 294), and by the German Federal Ministry of Education and Research (BMBF FRISYS), is highly appreciated. CoGe is supported by the US National Science Foundation under Award Numbers IOS-339156 and IOS-1444490, CyVerse is supported by the U.S. National Science Foundation under Award Numbers DBI-0735191 and DBI-1265383. YK and ACC are grateful for support from the UK Biological Sciences and Biotechnology Research Council (Grant BB/F001797/1). KV acknowledges the Multidisciplinary Research Partnership ‘Bioinformatics: from nucleotides to networks’ Project (no 01MR0410W) of Ghent University. JC is grateful for support from the Spanish Ministerio de Economía y Competitividad (Grant AGL2013-43244-R). RSQ is grateful to Monsanto (St. Louis, MO, USA) for sequencing genomic DNA of P. patens accession Kaskaskia. The 1000 Plants (1 KP) initiative, led by GKSW, is funded by the Alberta Ministry of Innovation and Advanced Education, Alberta Innovates Technology Futures (AITF), Innovates Centres of Research Excellence (iCORE), Musea Ventures, BGI-Shenzhen and China National Genebank (CNGB). TW was supported by EMBO Long-Term Fellowships (ALTF 1166-2011) and by Marie Curie Actions (European Commission EMBOCOFUND2010, GA-2010-267146). The work conducted at PGSB was supported by the German Research Foundation (SFB924) and German Ministry of Education and Research (BMBF, 031A536/de.NBI).

High-throughput RNA sequencing (RNA-seq) has recently become the method of choice to define and analyze transcriptomes. For the model moss Physcomitrella patens, although this method has been used to help analyze specific perturbations, no overall reference dataset has yet been established. In the framework of the Gene Atlas project, the Joint Genome Institute selected P. patens as a flagship genome, opening the way to generate the first comprehensive transcriptome dataset for this moss. The first round of sequencing described here is composed of 99 independent libraries spanning 34 different developmental stages and conditions. Upon dataset quality control and processing through read mapping, 28 509 of the 34 361 v3.3 gene models (83) were detected to be expressed across the samples. Differentially expressed genes (DEGs) were calculated across the dataset to permit perturbation comparisons between conditions. The analysis of the three most distinct and abundant P. patens growth stages – protonema, gametophore and sporophyte – allowed us to define both general transcriptional patterns and stage-specific transcripts. As an example of variation of physico-chemical growth conditions, we detail here the impact of ammonium supplementation under standard growth conditions on the protonemal transcriptome. Finally, the cooperative nature of this project allowed us to analyze inter-laboratory variation, as 13 different laboratories around the world provided samples. We compare differences in the replication of experiments in a single laboratory and between different laboratories. © 2018 The Authors The Plant Journal © 2018 John Wiley Sons Ltd

The moss Physcomitrella patens is used both as an evo-devo model and biotechnological production system for metabolites and pharmaceuticals. Strong in vivo expression of genes of interest is important for production of recombinant proteins, e.g., selectable markers, fluorescent proteins, or enzymes. In this regard, the choice of the promoter sequence as well as codon usage optimization are two important inside factors to consider in order to obtain optimum protein accumulation level. To reliably quantify fluorescence, we transfected protoplasts with promoter:GFP fusion constructs and measured fluorescence intensity of living protoplasts in a plate reader system. We used the red fluorescent protein mCherry under 2x 35S promoter control as second reporter to normalize for different transfection efficiencies. We derived a novel endogenous promoter and compared deletion variants with exogenous promoters. We used different codon-adapted green fluorescent protein (GFP) genes to evaluate the influence of promoter choice and codon optimization on protein accumulation in P. patens, and show that the promoter of the gene of P. patens chlorophyll a/b binding protein lhcsr1 drives expression of GFP in protoplasts significantly (more than twofold) better than the commonly used 2x 35S promoter or the rice actin1 promoter. We identified a shortened 677 bp version of the lhcsr1 promoter that retains full activity in protoplasts. The codon optimized GFP yields significantly (more than twofold) stronger fluorescence signals and thus demonstrates that adjusting codon usage in P. patens can increase expression strength. In combination, new promotor and codon optimized GFP conveyed sixfold increased fluorescence signal.

The moss Physcomitrella patens is an important model organism for studying plant evolution, development, physiology and biotechnology. Here we have generated microarray gene expression data covering the principal developmental stages, culture forms and some environmental/stress conditions. Example analyses of developmental stages and growth conditions as well as abiotic stress treatments demonstrate that (i) growth stage is dominant over culture conditions, (ii) liquid culture is not stressful for the plant, (iii) low pH might aid protoplastation by reduced expression of cell wall structure genes, (iv) largely the same gene pool mediates response to dehydration and rehydration, and (v) AP2/EREBP transcription factors play important roles in stress response reactions. With regard to the AP2 gene family, phylogenetic analysis and comparison with Arabidopsis thaliana shows commonalities as well as uniquely expressed family members under drought, light perturbations and protoplastation. Gene expression profiles for P. patens are available for the scientific community via the easy-to-use tool at https://www.genevestigator.com. By providing large-scale expression profiles, the usability of this model organism is further enhanced, for example by enabling selection of control genes for quantitative real-time PCR. Now, gene expression levels across a broad range of conditions can be accessed online for P. patens.

Rich ecotype collections are used for several plant models to unravel the molecular causes of phenotypic differences, and to investigate the effects of environmental adaption and acclimation. For the model moss Physcomitrella patens collections of accessions are available, and have been used for phylogenetic and taxonomic studies, for example, but few have been investigated further for phenotypic differences. Here, we focus on the Reute accession and provide expression profiling and comparative developmental data for several stages of sporophyte development, as well as information on genetic variation via genomic sequencing. We analysed cross-technology and cross-laboratory data to define a confident set of 15 mature sporophyte-specific genes. We find that the standard laboratory strain Gransden produces fewer sporophytes than Reute or Villersexel, although gametangia develop with the same time course and do not show evident morphological differences. Reute exhibits less genetic variation relative to Gransden than Villersexel, yet we found variation between Gransden and Reute in the expression profiles of several genes, as well as variation hot spots and genes that appear to evolve under positive Darwinian selection. We analyzed expression differences between the ecotypes for selected candidate genes in the GRAS transcription factor family, the chalcone synthase family and in genes involved in cell wall modification that are potentially related to phenotypic differences. We confirm that Reute is a P. patens ecotype, and suggest its use for reverse-genetics studies that involve progression through the life cycle and multiple generations.

The biotechnological production of proteins is currently achieved via expression systems derived from different lineages. In the past years transgenic plants have proven to be able to compete with bacteria or mammalian cell systems. Gene engineering approaches exist to raise yields by controlling mandatory processes in the course of biopharmaceutical protein production. Here we review and discuss the current status and recent improvements of parameters influencing recombinant protein production in transgenic plants. In particular, this review focuses on the so-called inside (mRNA sequence and structure) and outside factors (host and production system/conditions), which are adjustable and allow to optimize protein production via gene engineering.

The Plant Ontology (PO; http://www.plantontology.org/) is a publicly available, collaborative effort to develop and maintain a controlled, structured vocabulary (‘ontology’) of terms to describe plant anatomy, morphology and the stages of plant development. The goals of the PO are to link (annotate) gene expression and phenotype data to plant structures and stages of plant development, using the data model adopted by the Gene Ontology. From its original design covering only rice, maize and Arabidopsis, the scope of the PO has been expanded to include all green plants. The PO was the first multispecies anatomy ontology developed for the annotation of genes and phenotypes. Also, to our knowledge, it was one of the first biological ontologies that provides translations (via synonyms) in non-English languages such as Japanese and Spanish. As of Release #18 (July 2012), there are about 2.2 million annotations linking PO terms to >110,000 unique data objects representing genes or gene models, proteins, RNAs, germplasm and quantitative trait loci (QTLs) from 22 plant species. In this paper, we focus on the plant anatomical entity branch of the PO, describing the organizing principles, resources available to users and examples of how the PO is integrated into other plant genomics databases and web portals. We also provide two examples of comparative analyses, demonstrating how the ontology structure and PO-annotated data can be used to discover the patterns of expression of the LEAFY (LFY) and terpene synthase (TPS) gene homologs.