Representation and high-quality annotation of the Physcomitrella patens transcriptome demonstrates a high proportion of proteins involved in metabolism in mosses

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.

Agrobacterium-mediated Tnt1 mutagenesis of moss protonemal filaments and generation of stable mutants with impaired gametophyte

The gametophyte of moss exhibits a simple body plan, yet its growth is regulated by complex developmental phenomena similar to angiosperms. Because moss can be easily maintained under laboratory conditions, amenable for gene targeting and the availability of genome sequence, P. patens has become an attractive model system for studying evolutionary traits. Until date, there has been no Agrobacterium-mediated Tnt1 mutagenesis protocol for haploid protonemal filaments of moss. Hence, we attempted to use the intact tobacco Tnt1 retrotransposon as a mutagen for P. patens. Bioinformatic analysis of initiator methionyl-tRNA (Met-tRNAi), a critical host factor for Tnt1 transposition process, suggested that it can be explored as a mutagen for bryophytes. Using protonemal filaments and Agrobacterium-mediated transformation, 75 Tnt1 mutants have been generated and cryopreserved. SSAP analysis and TAIL-PCR revealed that Tnt1 is functional in P. patens and has a high-preference for gene and GC-rich regions. In addition, LTR::GUS lines exhibited a basal but tissue-specific inducible expression pattern. Forward genetic screen resulted in 5 novel phenotypes related to hormonal and gravity response, phyllid, and gamete development. SSAP analysis suggests that the Tnt1 insertion pattern is stable under normal growth conditions and the high-frequency phenotypic deviations are possibly due to the combination of haploid explant (protonema) and the choice of mutagen (Tnt1). We demonstrate that Agrobacterium-mediated Tnt1 insertional mutagenesis could generate stable P. patens mutant populations for future forward genetic studies.

Abscisic acid, cold and salt stimulate conserved metabolic regulation in the moss Physcomitrella patens

Salt and cold are major abiotic stresses that have adverse effects on plant growth and development. To cope with these stresses and their detrimental effects plants have evolved several metabolic, biochemical and physiological processes that are mainly triggered and mediated by the plant hormone abscisic acid (ABA). To elucidate the metabolic responses of the moss Physcomitrella patens, which serves as a model plant for abiotic stress adaptation, we performed GC-MS-based metabolic profiling of plants challenged for 5 and 28 h with either salt, cold or ABA. Our results indicate significant changes in the accumulation of several sugars including maltose, isomaltose and trehalose, amino acids including arginine, histidine, ornithine, tryptophan and tyrosine, and organic acids mainly citric acid and malonic acid. The metabolic responses provoked by ABA, cold and salt show considerable similarities. The accumulation of certain metabolites positively correlates with gene expression data whereas some metabolites do not show correlation with cognate transcript abundance. To place our results into an evolutionary context we compared the ABA- and stress-induced metabolic changes in moss to available metabolic profiles of the seed plant Arabidopsis thaliana. We detected considerable conservation between the species, indicating early evolution of stress-associated metabolic adaptations that probably occurred at the plant water-to-land transition.

Identification of Putative Precursor Genes for the Biosynthesis of Cannabinoid-Like Compound in Radula marginata

The liverwort Radula marginata belongs to the bryophyte division of land plants and is a prospective alternate source of cannabinoid-like compounds. However, mechanistic insights into the molecular pathways directing the synthesis of these cannabinoid-like compounds have been hindered due to the lack of genetic information. This prompted us to do deep sequencing, de novo assembly and annotation of R. marginata transcriptome, which resulted in the identification and validation of the genes for cannabinoid biosynthetic pathway. In total, we have identified 11,421 putative genes encoding 1,554 enzymes from 145 biosynthetic pathways. Interestingly, we have identified all the upstream genes of the central precursor of cannabinoid biosynthesis, cannabigerolic acid (CBGA), including its two first intermediates, stilbene acid (SA) and geranyl diphosphate (GPP). Expression of all these genes was validated using quantitative real-time PCR. We have characterized the protein structure of stilbene synthase (STS), which is considered as a homolog of olivetolic acid in R. marginata. Moreover, the metabolomics approach enabled us to identify CBGA-analogous compounds using electrospray ionization mass spectrometry (ESI-MS/MS) and gas chromatography mass spectrometry (GC-MS). Transcriptomic analysis revealed 1085 transcription factors (TF) from 39 families. Comparative analysis showed that six TF families have been uniquely predicted in R. marginata. In addition, the bioinformatics analysis predicted a large number of simple sequence repeats (SSRs) and non-coding RNAs (ncRNAs). Our results collectively provide mechanistic insights into the putative precursor genes for the biosynthesis of cannabinoid-like compounds and a novel transcriptomic resource for R. marginata. The large-scale transcriptomic resource generated in this study would further serve as a reference transcriptome to explore the Radulaceae family.

Pilot satellitome analysis of the model plant, Physcomitrellapatens, revealed a transcribed and high-copy IGS related tandem repeat

Satellite DNA (satDNA) constitutes a substantial part of eukaryotic genomes. In the last decade, it has been shown that satDNA is not an inert part of the genome and its function extends beyond the nuclear membrane. However, the number of model plant species suitable for studying the novel horizons of satDNA functionality is low. Here, we explored the satellitome of the model "basal" plant, Physcomitrellapatens (Hedwig, 1801) Bruch & Schimper, 1849 (moss), which has a number of advantages for deep functional and evolutionary research. Using a newly developed pyTanFinder pipeline (https://github.com/Kirovez/pyTanFinder) coupled with fluorescence in situ hybridization (FISH), we identified five high copy number tandem repeats (TRs) occupying a long DNA array in the moss genome. The nuclear organization study revealed that two TRs had distinct locations in the moss genome, concentrating in the heterochromatin and knob-rDNA like chromatin bodies. Further genomic, epigenetic and transcriptomic analysis showed that one TR, named PpNATR76, was located in the intergenic spacer (IGS) region and transcribed into long non-coding RNAs (lncRNAs). Several specific features of PpNATR76 lncRNAs make them very similar with the recently discovered human lncRNAs, raising a number of questions for future studies. This work provides new resources for functional studies of satellitome in plants using the model organism P.patens, and describes a list of tandem repeats for further analysis.

Plant Carbonic Anhydrases: Structures, Locations, Evolution, and Physiological Roles

Carbonic anhydrases (CAs) are zinc metalloenzymes that catalyze the interconversion of CO2 and HCO3- and are ubiquitous in nature. Higher plants contain three evolutionarily distinct CA families, αCAs, βCAs, and γCAs, where each family is represented by multiple isoforms in all species. Alternative splicing of CA transcripts appears common; consequently, the number of functional CA isoforms in a species may exceed the number of genes. CAs are expressed in numerous plant tissues and in different cellular locations. The most prevalent CAs are those in the chloroplast, cytosol, and mitochondria. This diversity in location is paralleled in the many physiological and biochemical roles that CAs play in plants. In this review, the number and types of CAs in C3, C4, and crassulacean acid metabolism (CAM) plants are considered, and the roles of the α and γCAs are briefly discussed. The remainder of the review focuses on plant βCAs and includes the identification of homologs between species using phylogenetic approaches, a consideration of the inter- and intracellular localization of the proteins, along with the evidence for alternative splice forms. Current understanding of βCA tissue-specific expression patterns and what controls them are reviewed, and the physiological roles for which βCAs have been implicated are presented.

Adaptation Mechanisms in the Evolution of Moss Defenses to Microbes

Bryophytes, including mosses, liverworts and hornworts are early land plants that have evolved key adaptation mechanisms to cope with abiotic stresses and microorganisms. Microbial symbioses facilitated plant colonization of land by enhancing nutrient uptake leading to improved plant growth and fitness. In addition, early land plants acquired novel defense mechanisms to protect plant tissues from pre-existing microbial pathogens. Due to its evolutionary stage linking unicellular green algae to vascular plants, the non-vascular moss Physcomitrella patens is an interesting organism to explore the adaptation mechanisms developed in the evolution of plant defenses to microbes. Cellular and biochemical approaches, gene expression profiles, and functional analysis of genes by targeted gene disruption have revealed that several defense mechanisms against microbial pathogens are conserved between mosses and flowering plants. P. patens perceives pathogen associated molecular patterns by plasma membrane receptor(s) and transduces the signal through a MAP kinase (MAPK) cascade leading to the activation of cell wall associated defenses and expression of genes that encode proteins with different roles in plant resistance. After pathogen assault, P. patens also activates the production of ROS, induces a HR-like reaction and increases levels of some hormones. Furthermore, alternative metabolic pathways are present in P. patens leading to the production of a distinct metabolic scenario than flowering plants that could contribute to defense. P. patens has acquired genes by horizontal transfer from prokaryotes and fungi, and some of them could represent adaptive benefits for resistance to biotic stress. In this review, the current knowledge related to the evolution of plant defense responses against pathogens will be discussed, focusing on the latest advances made in the model plant P. patens.

The moss Funaria hygrometrica has cuticular wax similar to vascular plants, with distinct composition on leafy gametophyte, calyptra and sporophyte capsule surfaces

BACKGROUND AND AIMS

Aerial surfaces of land plants are covered with a waxy cuticle to protect against water loss. The amount and composition of cuticular waxes on moss surfaces had rarely been investigated. Accordingly, the degree of similarity between moss and vascular plant waxes, and between maternal and offspring moss structure waxes is unknown. To resolve these issues, this study aimed at providing a comprehensive analysis of the waxes on the leafy gametophyte, gametophyte calyptra and sporophyte capsule of the moss Funaria hygrometrica

METHODS

Waxes were extracted from the surfaces of leafy gametophytes, gametophyte calyptrae and sporophyte capsules, separated by gas chromatography, identified qualitatively with mass spectrometry, and quantified with flame ionization detection. Diagnostic mass spectral peaks were used to determine the isomer composition of wax esters.

KEY RESULTS

The surfaces of the leafy gametophyte, calyptra and sporophyte capsule of F. hygrometrica were covered with 0·94, 2·0 and 0·44 μg cm(-2) wax, respectively. While each wax mixture was composed of mainly fatty acid alkyl esters, the waxes from maternal and offspring structures had unique compositional markers. β-Hydroxy fatty acid alkyl esters were limited to the leafy gametophyte and calyptra, while alkanes, aldehydes and diol esters were restricted to the sporophyte capsule. Ubiquitous fatty acids, alcohols, fatty acid alkyl esters, aldehydes and alkanes were all found on at least one surface.

CONCLUSIONS

This is the first study to determine wax coverage (μg cm(-2)) on a moss surface, enabling direct comparisons with vascular plants, which were shown to have an equal amount or more wax than F. hygrometrica Wax ester biosynthesis is of particular importance in this species, and the ester-forming enzyme(s) in different parts of the moss may have different substrate preferences. Furthermore, the alkane-forming wax biosynthesis pathway, found widely in vascular plants, is active in the sporophyte capsule, but not in the leafy gametophyte or calyptra. Overall, wax composition and coverage on F. hygrometrica were similar to those reported for some vascular plant species, suggesting that the underlying biosynthetic processes in plants of both lineages were inherited from a common ancestor.

The mitochondrial proteome of the moss Physcomitrella patens

Extant basal land plants are routinely used to trace plant evolution and to track strategies for high abiotic stress resistance. Whereas the structure of mitochondrial genomes and RNA editing are already well studied, mitochondrial proteome research is restricted to a few data sets. While the mitochondrial proteome of the model moss Physcomitrella patens is covered to an estimated 15-25% by proteomic evidence to date, the available data have already provided insights into the evolution of metabolic compartmentation, dual targeting and mitochondrial heterogeneity. This review summarizes the current knowledge about the mitochondrial proteome of P. patens, and gives a perspective on its use as a mitochondrial model system. Its amenability to gene editing, metabolic labelling as well as fluorescence microscopy provides a unique platform to study open questions in mitochondrial biology, such as regulation of protein stability, responses to stress and connectivity to other organelles. Future challenges will include improving the proteomic resources for P. patens, and to link protein inventories and modifications as well as evolutionary differences to the functional level.

CHASE domain-containing receptors play an essential role in the cytokinin response of the moss Physcomitrella patens

While the molecular basis for cytokinin action is quite well understood in flowering plants, little is known about the cytokinin signal transduction in early diverging land plants. The genome of the bryophyte Physcomitrella patens (Hedw.) B.S. encodes three classical cytokinin receptors, the CHASE domain-containing histidine kinases, CHK1, CHK2, and CHK3. In a complementation assay with protoplasts of receptor-deficient Arabidopsis thaliana as well as in cytokinin binding assays, we found evidence that CHK1 and CHK2 receptors can function in cytokinin perception. Using gene targeting, we generated a collection of CHK knockout mutants comprising single (Δchk1, Δchk2, Δchk3), double (Δchk1,2, Δchk1,3, Δchk2,3), and triple (Δchk1,2,3) mutants. Mutants were characterized for their cytokinin response and differentiation capacities. While the wild type did not grow on high doses of cytokinin (1 µM benzyladenine), the Δchk1,2,3 mutant exhibited normal protonema growth. Bud induction assays showed that all three cytokinin receptors contribute to the triggering of budding, albeit to different extents. Furthermore, while the triple mutant showed no response in this bioassay, the remaining mutants displayed budding responses in a diverse manner to different types and concentrations of cytokinins. Determination of cytokinin levels in mutants showed no drastic changes for any of the cytokinins; thus, in contrast to Arabidopsis, revealing only small impacts of cytokinin signaling on homeostasis. In summary, our study provides a first insight into the molecular action of cytokinin in an early diverging land plant and demonstrates that CHK receptors play an essential role in bud induction and gametophore development.

Predicted protein-protein interactions in the moss Physcomitrella patens: a new bioinformatic resource

Background

Physcomitrella patens, a haploid dominant plant, is fast becoming a useful molecular genetics and bioinformatics tool due to its key phylogenetic position as a bryophyte in the post-genomic era. Genome sequences from select reference species were compared bioinformatically to Physcomitrella patens using reciprocal blasts with the InParanoid software package. A reference protein interaction database assembled using MySQL by compiling BioGrid, BIND, DIP, and Intact databases was queried for moss orthologs existing for both interacting partners. This method has been used to successfully predict interactions for a number of angiosperm plants.

Results

The first predicted protein-protein interactome for a bryophyte based on the interolog method contains 67,740 unique interactions from 5,695 different Physcomitrella patens proteins. Most conserved interactions among proteins were those associated with metabolic processes. Over-represented Gene Ontology categories are reported here.

Conclusion

Addition of moss, a plant representative 200 million years diverged from angiosperms to interactomic research greatly expands the possibility of conducting comparative analyses giving tremendous insight into network evolution of land plants. This work helps demonstrate the utility of “guilt-by-association” models for predicting protein interactions, providing provisional roadmaps that can be explored using experimental approaches. Included with this dataset is a method for characterizing subnetworks and investigating specific processes, such as the Calvin-Benson-Bassham cycle.

Electronic supplementary material

The online version of this article (doi:10.1186/s12859-015-0524-1) contains supplementary material, which is available to authorized users.

Tc1-like transposable elements in plant genomes

BACKGROUND

The Tc1/mariner superfamily of transposable elements (TEs) is widespread in animal genomes. Mariner-like elements, which bear a DDD triad catalytic motif, have been identified in a wide range of flowering plant species. However, as the founding member of the superfamily, Tc1-like elements that bear a DD34E triad catalytic motif are only known to unikonts (animals, fungi, and Entamoeba).

RESULTS

Here we report the identification of Tc1-like elements (TLEs) in plant genomes. These elements bear the four terminal nucleotides and the characteristic DD34E triad motif of Tc1 element. The two TLE families (PpTc1, PpTc2) identified in the moss (Physcomitrella patens) genome contain highly similar copies. Multiple copies of PpTc1 are actively transcribed and the transcripts encode intact full length transposase coding sequences. TLEs are also found in angiosperm genome sequence databases of rice (Oryza sativa), dwarf birch (Betula nana), cabbage (Brassica rapa), hemp (Cannabis sativa), barley (Hordium valgare), lettuce (Lactuta sativa), poplar (Populus trichocarpa), pear (Pyrus x bretschneideri), and wheat (Triticum urartu).

CONCLUSIONS

This study extends the occurrence of TLEs to the plant phylum. The elements in the moss genome have amplified recently and may still be capable of transposition. The TLEs are also present in angiosperm genomes, but apparently much less abundant than in moss.

Isopentenyltransferase-1 (IPT1) knockout in Physcomitrella together with phylogenetic analyses of IPTs provide insights into evolution of plant cytokinin biosynthesis

The moss Physcomitrella patens is part of an early divergent clade of land plants utilizing the plant hormone cytokinin for growth control. The rate-limiting step of cytokinin biosynthesis is mediated by isopentenyltransferases (IPTs), found in land plants either as adenylate-IPTs or as tRNA-IPTs. Although a dominant part of cytokinins in flowering plants are synthesized by adenylate-IPTs, the Physcomitrella genome only encodes homologues of tRNA-IPTs. This study therefore looked into the question of whether cytokinins in moss derive from tRNA exclusively. Targeted gene knockout of ipt1 (d|ipt1) along with localization studies revealed that the chloroplast-bound IPT1 was almost exclusively responsible for the A37 prenylation of tRNA in Physcomitrella. Ultra-performance liquid chromatography-tandem mass spectrometry (UPLC-MS/MS)-based cytokinin profiling demonstrated that the total amount of all free cytokinins in tissue was almost unaffected. However, the knockout plants showed increased levels of the N (6) -isopentenyladenine (iP)- and trans-zeatin (tZ)-type cytokinins, considered to provide active forms, while cis-zeatin (cZ)-type cytokinins were reduced. The data provide evidence for an additional and unexpected tRNA-independent cytokinin biosynthetic pathway in moss. Comprehensive phylogenetic analysis indicates a diversification of tRNA-IPT-like genes in bryophytes probably related to additional functions.

Reannotation and extended community resources for the genome of the non-seed plant Physcomitrella patens provide insights into the evolution of plant gene structures and functions

Background

The moss Physcomitrella patens as a model species provides an important reference for early-diverging lineages of plants and the release of the genome in 2008 opened the doors to genome-wide studies. The usability of a reference genome greatly depends on the quality of the annotation and the availability of centralized community resources. Therefore, in the light of accumulating evidence for missing genes, fragmentary gene structures, false annotations and a low rate of functional annotations on the original release, we decided to improve the moss genome annotation.

Results

Here, we report the complete moss genome re-annotation (designated V1.6) incorporating the increased transcript availability from a multitude of developmental stages and tissue types. We demonstrate the utility of the improved P. patens genome annotation for comparative genomics and new extensions to the cosmoss.org resource as a central repository for this plant “flagship” genome. The structural annotation of 32,275 protein-coding genes results in 8387 additional loci including 1456 loci with known protein domains or homologs in Plantae. This is the first release to include information on transcript isoforms, suggesting alternative splicing events for at least 10.8% of the loci. Furthermore, this release now also provides information on non-protein-coding loci. Functional annotations were improved regarding quality and coverage, resulting in 58% annotated loci (previously: 41%) that comprise also 7200 additional loci with GO annotations. Access and manual curation of the functional and structural genome annotation is provided via the http://www.cosmoss.org model organism database.

Conclusions

Comparative analysis of gene structure evolution along the green plant lineage provides novel insights, such as a comparatively high number of loci with 5’-UTR introns in the moss. Comparative analysis of functional annotations reveals expansions of moss house-keeping and metabolic genes and further possibly adaptive, lineage-specific expansions and gains including at least 13% orphan genes.

The plant ontology as a tool for comparative plant anatomy and genomic analyses

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.

Aldehyde dehydrogenase (ALDH) superfamily in plants: gene nomenclature and comparative genomics

In recent years, there has been a significant increase in the number of completely sequenced plant genomes. The comparison of fully sequenced genomes allows for identification of new gene family members, as well as comprehensive analysis of gene family evolution. The aldehyde dehydrogenase (ALDH) gene superfamily comprises a group of enzymes involved in the NAD(+)- or NADP(+)-dependent conversion of various aldehydes to their corresponding carboxylic acids. ALDH enzymes are involved in processing many aldehydes that serve as biogenic intermediates in a wide range of metabolic pathways. In addition, many of these enzymes function as 'aldehyde scavengers' by removing reactive aldehydes generated during the oxidative degradation of lipid membranes, also known as lipid peroxidation. Plants and animals share many ALDH families, and many genes are highly conserved between these two evolutionarily distinct groups. Conversely, both plants and animals also contain unique ALDH genes and families. Herein we carried out genome-wide identification of ALDH genes in a number of plant species-including Arabidopsis thaliana (thale crest), Chlamydomonas reinhardtii (unicellular algae), Oryza sativa (rice), Physcomitrella patens (moss), Vitis vinifera (grapevine) and Zea mays (maize). These data were then combined with previous analysis of Populus trichocarpa (poplar tree), Selaginella moellindorffii (gemmiferous spikemoss), Sorghum bicolor (sorghum) and Volvox carteri (colonial algae) for a comprehensive evolutionary comparison of the plant ALDH superfamily. As a result, newly identified genes can be more easily analyzed and gene names can be assigned according to current nomenclature guidelines; our goal is to clarify previously confusing and conflicting names and classifications that might confound results and prevent accurate comparisons between studies.

Identification of cyclic nucleotide gated channels using regular expressions

Cyclic nucleotide-gated channels (CNGCs) are nonselective cation channels found in plants, animals, and some bacteria. They have a six-transmembrane/one-pore structure, a cytosolic cyclic nucleotide-binding domain, and a cytosolic calmodulin-binding domain. Despite their functional similarities, the plant CNGC family members appear to have different conserved amino acid motifs within corresponding functional domains than animal and bacterial CNGCs do. Here we describe the development and application of methods employing plant CNGC-specific sequence motifs as diagnostic tools to identify novel candidate channels in different plants. These methods are used to evaluate the validity of annotations of putative orthologs of CNGCs from plant genomes. The methods detail how to employ regular expressions of conserved amino acids in functional domains of annotated CNGCs and together with Web tools such as PHI-BLAST and ScanProsite to identify novel candidate CNGCs in species including Physcomitrella patens.

A CELLULOSE SYNTHASE (CESA) gene essential for gametophore morphogenesis in the moss Physcomitrella patens

In seed plants, different groups of orthologous genes encode the CELLULOSE SYNTHASE (CESA) proteins that are responsible for cellulose biosynthesis in primary and secondary cell walls. The seven CESA sequences of the moss Physcomitrella patens (Hedw.) B. S. G. form a monophyletic sister group to seed plant CESAs, consistent with independent CESA diversification and specialization in moss and seed plant lines. The role of PpCESA5 in the development of P. patens was investigated by targeted mutagenesis. The cesa5 knockout lines were tested for cellulose deficiency using carbohydrate-binding module affinity cytochemistry and the morphology of the leafy gametophores was analyzed by 3D reconstruction of confocal images. No defects were identified in the development of the filamentous protonema or in production of bud initials that normally give rise to the leafy gametophores. However, the gametophore buds were cellulose deficient and defects in subsequent cell expansion, cytokinesis, and leaf initiation resulted in the formation of irregular cell clumps instead of leafy shoots. Analysis of the cesa5 knockout phenotype indicates that a biophysical model of organogenesis can be extended to the moss gametophore shoot apical meristem.

Glycoprotein production in moss bioreactors

Complex multimeric recombinant proteins such as therapeutic antibodies require a eukaryotic expression system. Transgenic plants may serve as promising alternatives to the currently favored mammalian cell lines or hybridomas. In contrast to prokaryotic systems, posttranslational modifications of plant and human proteins resemble each other largely, among those, protein N-glycosylation of the complex type. However, a few plant-specific sugar residues may cause immune reactions in humans, representing an obstacle for the broad use of plant-based systems as biopharmaceutical production hosts. The moss Physcomitrella patens represents a flexible tissue-culture system for the contained production and secretion of recombinant biopharmaceuticals in photobioreactors. The recent synthesis of therapeutic proteins as a scFv antibody fragment or the large and heavily modified complement regulator factor H demonstrate the versatility of this expression system. A uniquely efficient gene targeting mechanism can be employed to precisely engineer the glycosylation machinery for recombinant products. In this way, P. patens lines with non-immunogenic optimized glycan structures were created. Therapeutic antibodies produced in these strains exhibited antibody-dependent cellular cytotoxicity superior to the same molecules synthesized in mammalian cell lines.

Metabolite profiling of the moss Physcomitrella patens reveals evolutionary conservation of osmoprotective substances

The moss Physcomitrella patens is suitable for systems biology studies, as it can be grown axenically under standardised conditions in plain mineral medium and comprises only few cell types. We report on metabolite profiling of two major P. patens tissues, filamentous protonema and leafy gametophores, from different culture conditions. A total of 96 compounds were detected, 21 of them as yet unknown in public databases. Protonema and gametophores had distinct metabolic profiles, especially with regard to saccharides, sugar derivates, amino acids, lignin precursors and nitrogen-rich storage compounds. A hydroponic culture was established for P. patens, and was used to apply drought stress under physiological conditions. This treatment led to accumulation of osmoprotectants, such as altrose, maltitol, ascorbic acid and proline. Thus, these osmoprotectants are not unique to seed plants but have evolved at an early phase of the colonization of land by plants.

Evolution of the land plant exocyst complexes

Exocyst is an evolutionarily conserved vesicle tethering complex functioning especially in the last stage of exocytosis. Homologs of its eight canonical subunits - Sec3, Sec5, Sec6, Sec8, Sec10, Sec15, Exo70, and Exo84 - were found also in higher plants and confirmed to form complexes in vivo, and to participate in cell growth including polarized expansion of pollen tubes and root hairs. Here we present results of a phylogenetic study of land plant exocyst subunits encoded by a selection of completely sequenced genomes representing a variety of plant, mostly angiosperm, lineages. According to their evolution histories, plant exocyst subunits can be divided into several groups. The core subunits Sec6, Sec8, and Sec10, together with Sec3 and Sec5, underwent few, if any fixed duplications in the tracheophytes (though they did amplify in the moss Physcomitrella patens), while others form larger families, with the number of paralogs ranging typically from two to eight per genome (Sec15, Exo84) to several dozens per genome (Exo70). Most of the diversity, which can be in some cases traced down to the origins of land plants, can be attributed to the peripheral subunits Exo84 and, in particular, Exo70. As predicted previously, early land plants (including possibly also the Rhyniophytes) encoded three ancestral Exo70 paralogs which further diversified in the course of land plant evolution. Our results imply that plants do not have a single "Exocyst complex" - instead, they appear to possess a diversity of exocyst variants unparalleled among other organisms studied so far. This feature might perhaps be directly related to the demands of building and maintenance of the complicated and spatially diverse structures of the endomembranes and cell surfaces in multicellular land plants.

Ammonium and urea transporter inventory of the selaginella and physcomitrella genomes

Ammonium and urea are important nitrogen sources for autotrophic organisms. Plant genomes encode several families of specific transporters for these molecules, plus other uptake mechanisms such as aquaporins and ABC transporters. Selaginella and Physcomitrella are representatives of lycophytes and bryophytes, respectively, and the recent completion of their genome sequences provided us with an opportunity for comparative genome studies, with special emphasis on the adaptive processes that accompanied the conquest of dry land and the evolution of a vascular system. Our phylogenetic analysis revealed that the number of genes encoding urea transporters underwent a progressive reduction during evolution, eventually down to a single copy in vascular plants. Conversely, no clear evolutionary pattern was found for ammonium transporters, and their number and distribution in families varies between species. In particular Selaginella, similar to rice, favors the AMT2/MEP family of ammonium transporters over the plant-specific AMT1 type. In comparison, Physcomitrella presents several members belonging to both families.

PpASCL, a moss ortholog of anther-specific chalcone synthase-like enzymes, is a hydroxyalkylpyrone synthase involved in an evolutionarily conserved sporopollenin biosynthesis pathway

Sporopollenin is the main constituent of the exine layer of spore and pollen walls. Recently, several Arabidopsis genes, including polyketide synthase A (PKSA), which encodes an anther-specific chalcone synthase-like enzyme (ASCL), have been shown to be involved in sporopollenin biosynthesis. The genome of the moss Physcomitrella patens contains putative orthologs of the Arabidopsis sporopollenin biosynthesis genes. We analyzed available P.patens expressed sequence tag (EST) data for putative moss orthologs of the Arabidopsis genes of sporopollenin biosynthesis and studied the enzymatic properties and reaction mechanism of recombinant PpASCL, the P.patens ortholog of Arabidopsis PKSA. We also generated structure models of PpASCL and Arabidopsis PKSA to study their substrate specificity. Physcomitrella patens orthologs of Arabidopsis genes for sporopollenin biosynthesis were found to be expressed in the sporophyte generation. Similarly to Arabidopsis PKSA, PpASCL condenses hydroxy fatty acyl-CoA esters with malonyl-CoA and produces hydroxyalkyl α-pyrones that probably serve as building blocks of sporopollenin. The ASCL-specific set of Gly-Gly-Ala residues predicted by the models to be located at the floor of the putative active site is proposed to serve as the opening of an acyl-binding tunnel in ASCL. These results suggest that ASCL functions together with other sporophyte-specific enzymes to provide polyhydroxylated precursors of sporopollenin in a pathway common to land plants.

Functional diversification of thylakoidal processing peptidases in Arabidopsis thaliana

Thylakoidal processing peptidase (TPP) is responsible for removing amino-terminal thylakoid-transfer signals from several proteins in the thylakoid lumen. Three TPP isoforms are encoded by the nuclear genome of Arabidopsis thaliana. Previous studies showed that one of them termed plastidic type I signal peptidase 1 (Plsp1) was necessary for processing three thylakoidal proteins and one protein in the chloroplast envelope in vivo. The lack of Plsp1 resulted in seedling lethality, apparently due to disruption of proper thylakoid development. The physiological roles of the other two TPP homologs remain unknown. Here we show that the three A. thaliana TPP isoforms evolved to acquire diverse functions. Phylogenetic analysis revealed that TPP may have originated before the endosymbiotic event, and that there are two groups of TPP in seed plants: one includes Plsp1 and another comprises the other two A. thaliana TPP homologs, which are named as Plsp2A and Plsp2B in this study. The duplication leading to the two groups predates the gymnosperm-angiosperm divergence, and the separation of Plsp2A and Plsp2B occurred after the Malvaceae-Brassicaceae diversification. Quantitative reverse transcription-PCR assay revealed that the two PLSP2 genes were co-expressed in both photosynthetic tissues and roots, whereas the PLSP1 transcript accumulated predominantly in photosynthetic tissues. Both PLSP2 genes were expressed in the aerial parts of the plsp1-null mutant at levels comparable to those in wild-type plants. The seedling-lethal phenotype of the plsp1-null mutant could be rescued by a constitutive expression of Plsp1 cDNA but not by that of Plsp2A or Plsp2B. These results indicate that Plsp1 and Plsp2 evolved to function differently, and that neither of the Plsp2 isoforms is necessary for proper thylakoid development in photosynthetic tissues.

Two novel types of hexokinases in the moss Physcomitrella patens

BACKGROUND

Hexokinase catalyzes the phosphorylation of glucose and fructose, but it is also involved in sugar sensing in both fungi and plants. We have previously described two types of hexokinases in the moss Physcomitrella. Type A, exemplified by PpHxk1, the major hexokinase in Physcomitrella, is a soluble protein that localizes to the chloroplast stroma. Type B, exemplified by PpHxk2, has an N-terminal membrane anchor. Both types are found also in vascular plants, and localize to the chloroplast stroma and mitochondrial membranes, respectively.

RESULTS

We have now characterized all 11 hexokinase encoding genes in Physcomitrella. Based on their N-terminal sequences and intracellular localizations, three of the encoded proteins are type A hexokinases and four are type B hexokinases. One of the type B hexokinases has a splice variant without a membrane anchor, that localizes to the cytosol and the nucleus. However, we also found two new types of hexokinases with no obvious orthologs in vascular plants. Type C, encoded by a single gene, has neither transit peptide nor membrane anchor, and is found in the cytosol and in the nucleus. Type D hexokinases, encoded by three genes, have membrane anchors and localize to mitochondrial membranes, but their sequences differ from those of the type B hexokinases. Interestingly, all moss hexokinases are more similar to each other in overall sequence than to hexokinases from other plants, even though characteristic sequence motifs such as the membrane anchor of the type B hexokinases are highly conserved between moss and vascular plants, indicating a common origin for hexokinases of the same type.

CONCLUSIONS

We conclude that the hexokinase gene family is more diverse in Physcomitrella, encoding two additional types of hexokinases that are absent in vascular plants. In particular, the presence of a cytosolic and nuclear hexokinase (type C) sets Physcomitrella apart from vascular plants, and instead resembles yeast, where all hexokinases localize to the cytosol. The fact that all moss hexokinases are more similar to each other than to hexokinases from vascular plants, even though both type A and type B hexokinases are present in all plants, further suggests that the hexokinase gene family in Physcomitrella has undergone concerted evolution.

In silico comparative analysis of SSR markers in plants

BACKGROUND

The adverse environmental conditions impose extreme limitation to growth and plant development, restricting the genetic potential and reflecting on plant yield losses. The progress obtained by classic plant breeding methods aiming at increasing abiotic stress tolerances have not been enough to cope with increasing food demands. New target genes need to be identified to reach this goal, which requires extensive studies of the related biological mechanisms. Comparative analyses in ancestral plant groups can help to elucidate yet unclear biological processes.

RESULTS

In this study, we surveyed the occurrence patterns of expressed sequence tag-derived microsatellite markers for model plants. A total of 13,133 SSR markers were discovered using the SSRLocator software in non-redundant EST databases made for all eleven species chosen for this study. The dimer motifs are more frequent in lower plant species, such as green algae and mosses, and the trimer motifs are more frequent for the majority of higher plant groups, such as monocots and dicots. With this in silico study we confirm several microsatellite plant survey results made with available bioinformatics tools.

CONCLUSIONS

The comparative studies of EST-SSR markers among all plant lineages is well suited for plant evolution studies as well as for future studies of transferability of molecular markers.

Functional analysis of the group 4 late embryogenesis abundant proteins reveals their relevance in the adaptive response during water deficit in Arabidopsis

Late-Embryogenesis Abundant (LEA) proteins accumulate to high levels during the last stages of seed development, when desiccation tolerance is acquired, and in vegetative and reproductive tissues under water deficit, leading to the hypothesis that these proteins play a role in the adaptation of plants to this stress condition. In this work, we obtained the accumulation patterns of the Arabidopsis (Arabidopsis thaliana) group 4 LEA proteins during different developmental stages and plant organs in response to water deficit. We demonstrate that overexpression of a representative member of this group of proteins confers tolerance to severe drought in Arabidopsis plants. Moreover, we show that deficiency of LEA proteins in this group leads to susceptible phenotypes upon water limitation, during germination, or in mature plants after recovery from severe dehydration. Upon recovery from this stress condition, mutant plants showed a reduced number of floral and axillary buds when compared with wild-type plants. The lack of these proteins also correlates with a reduced seed production under optimal irrigation, supporting a role in fruit and/or seed development. A bioinformatic analysis of group 4 LEA proteins from many plant genera showed that there are two subgroups, originated through ancient gene duplication and a subsequent functional specialization. This study represents, to our knowledge, the first genetic evidence showing that one of the LEA protein groups is directly involved in the adaptive response of higher plants to water deficit, and it provides data indicating that the function of these proteins is not redundant to that of the other LEA proteins.

Regulation of sulfate assimilation in Physcomitrella patens: mosses are different!

Sulfur is an essential nutrient, taken up as sulfate from soil, reduced and incorporated into bioorganic compounds in plant cells. The pathway of sulfate assimilation is highly regulated in a demand-driven manner in seed plants. To test the evolutionary conservation of the regulatory mechanisms, we analyzed regulation of the pathway in the model for basal plants, the moss Physcomitrella patens. While in Arabidopsis the key enzyme of sulfate assimilation, adenosine 5'-phosphosulfate reductase (APR), is feedback repressed by thiols and induced by reduced levels of glutathione, in P. patens such regulation does not occur. The control of the pathway was not moved to other components as these conditions affected neither mRNA accumulation of other genes of sulfate assimilation nor sulfate uptake. Other treatments known to regulate APR, O-acetylserine, cadmium and sulfur deficiency affected APR transcript levels, but not enzyme activity. It appears that the sulfate assimilation pathway in P. patens is much more robust than in seed plants. Thus, the regulatory networks controlling the pathway have probably evolved only later in the evolution of the seed plants after separation of the bryophytes.

Targeted knock-out of a gene encoding sulfite reductase in the moss Physcomitrella patens affects gametophytic and sporophytic development

A key step in sulfate assimilation into cysteine is the reduction of sulfite to sulfide by sulfite reductase (SiR). This enzyme is encoded by three genes in the moss Physcomitrella patens. To obtain a first insight into the roles of the individual isoforms, we deleted the gene encoding the SiR1 isoform in P. patens by homologous recombination and subsequently analysed the DeltaSiR1 mutants. While DeltaSiR1 mutants showed no obvious alteration in sulfur metabolism, their regeneration from protoplasts and their ability to produce mature spores was significantly affected, highlighting an unexpected link between moss sulfate assimilation and development, that is yet to be characterized.

Genome-wide analysis of the chalcone synthase superfamily genes of Physcomitrella patens

Enzymes of the chalcone synthase (CHS) superfamily catalyze the production of a variety of secondary metabolites in bacteria, fungi and plants. Some of these metabolites have played important roles during the early evolution of land plants by providing protection from various environmental assaults including UV irradiation. The genome of the moss, Physcomitrella patens, contains at least 17 putative CHS superfamily genes. Three of these genes (PpCHS2b, PpCHS3 and PpCHS5) exist in multiple copies and all have corresponding ESTs. PpCHS11 and probably also PpCHS9 encode non-CHS enzymes, while PpCHS10 appears to be an ortholog of plant genes encoding anther-specific CHS-like enzymes. It was inferred from the genomic locations of genes comprising it that the moss CHS superfamily expanded through tandem and segmental duplication events. Inferred exon-intron architectures and results from phylogenetic analysis of representative CHS superfamily genes of P. patens and other plants showed that intron gain and loss occurred several times during evolution of this gene superfamily. A high proportion of P. patens CHS genes (7 of 14 genes for which the full sequence is known and probably 3 additional genes) are intronless, prompting speculation that CHS gene duplication via retrotransposition has occurred at least twice in the moss lineage. Analyses of sequence similarities, catalytic motifs and EST data indicated that a surprisingly large number (as many as 13) of the moss CHS superfamily genes probably encode active CHS. EST distribution data and different light responsiveness observed with selected genes provide evidence for their differential regulation. Observed diversity within the moss CHS superfamily and amenability to gene manipulation make Physcomitrella a highly suitable model system for studying expansion and functional diversification of the plant CHS superfamily of genes.

Microarray analysis of the moss Physcomitrella patens reveals evolutionarily conserved transcriptional regulation of salt stress and abscisic acid signalling

Regulatory networks of salt stress and abscisic acid (ABA) responses have previously been analyzed in seed plants. Here, we report microarray expression profiles of 439 genes encoding transcription-associated proteins (TAPs) in response to salt stress and ABA in the salt-tolerant moss Physcomitrella patens. Fourteen and 56 TAP genes were differentially expressed within 60 min of NaCl and ABA treatment, respectively, indicating that these responses are regulated at the transcriptional level. Overlapping expression profiles, as well as the up-regulation of ABA biosynthesis genes, suggest that ABA mediates the salt stress responses in P. patens. Comparison to public gene expression data of Arabidopsis thaliana and phylogenetic analyses suggest that the role of DREB-like, Dof, and bHLH TAPs in salt stress responses have been conserved during embryophyte evolution, and that the function of ABI3-like, bZIP, HAP3, and CO-like TAPs in seed development and flowering emerged from pre-existing ABA and light signalling pathways.

Proteomic analysis of the cold stress response in the moss, Physcomitrella patens

Cold stress has adverse effects on plant growth and development. Plants respond and acclimate to cold stress through various biochemical and physiological processes, thereby acquiring stress tolerance. To better understand the basis for tolerance, we carried out a proteomic study in the model moss, Physcomitrella patens, characterizing gametophore proteins with 2-DE and mass spectroscopy. Following exposure to 0 degrees C for up to 3 days, out of the more than 1000 protein spots reproducibly resolved, only 45 changed in abundance by at least 1.5-fold. Of these, 35 were identified by tryptic digestion and mass spectroscopy. Photosynthetic proteins decreased, whereas many catabolic proteins increased. In addition, cold stress up-regulated a variety of signaling, cytoskeleton, and defense proteins and few proteins in these classes were down-regulated. Up-regulated proteins include the 14-3-3-like protein, actin, HSP70s, lipoxygenases, and cytochrome P450 proteins. These results point to pathways that are important for the mechanism of cold stress response in P. patens and by extension to the entire plant kingdom.

Pythium infection activates conserved plant defense responses in mosses

The moss Physcomitrella patens (P. patens) is a useful model to study abiotic stress responses since it is highly tolerant to drought, salt and osmotic stress. However, very little is known about the defense mechanisms activated in this moss after pathogen assault. In this study, we show that P. patens activated multiple and similar responses against Pythium irregulare and Pythium debaryanum, including the reinforcement of the cell wall, induction of the defense genes CHS, LOX and PAL, and accumulation of the signaling molecules jasmonic acid (JA) and its precursor 12-oxo-phytodienoic acid (OPDA). However, theses responses were not sufficient and infection could not be prevented leading to hyphae colonization of moss tissues and plant decay. Pythium infection induced reactive oxygen species production and caused cell death of moss tissues. Taken together, these data indicate that Pythium infection activates in P. patens common responses to those previously characterized in flowering plants. Microscopic analysis also revealed intracellular relocation of chloroplasts in Pythium-infected tissues toward the infection site. In addition, OPDA, JA and its methyl ester methyl jasmonate induced the expression of PAL. Our results show for the first time JA and OPDA accumulation in a moss and suggest that this defense pathway is functional and has been maintained during the evolution of plants.

Exploring plant biodiversity: the Physcomitrella genome and beyond

For decades, plant molecular biology has focused on only a few angiosperm species. Recently, the approximately 500mega base pairs (Mb) of the haploid Physcomitrella patens genome were sequenced and annotated. Mosses such as P. patens occupy a key evolutionary position halfway between green algae and flowering plants. This draft genome, in comparison to existing genome data from other plants, allows evolutionary insights into the conquest of land by plants and the molecular biodiversity that land plants exhibit. As a model organism, P. patens provides a well-developed molecular toolbox, including efficient gene targeting in combination with the morphologically simple moss tissues. We describe current as well as future tools for P. patens research and the prospects they offer for plant research in general.

Identification of a 4-coumarate:CoA ligase gene family in the moss, Physcomitrella patens

Since the early evolution of land plants from primitive green algae, phenylpropanoid compounds have played an important role. In the biosynthesis of phenylpropanoids, 4-coumarate:CoA ligase (4CL; EC 6.2.1.12) has a pivotal role at the divergence point from general phenylpropanoid metabolism to several major branch pathways. Although higher plant 4CLs have been extensively studied, little information is available on the enzymes from bryophytes. In Physcomitrella patens, we have identified a 4CL gene family consisting of four members, taking advantage of the available EST sequences and a draft sequence of the P. patens genome. The encoded proteins of three of the genes display similar substrate utilization profiles with highest catalytic efficiency towards 4-coumarate. Interestingly, the efficiency with cinnamate as substrate is in the same range as with caffeate and ferulate. The deduced proteins of the four genes share sequence identities between 78% and 86%. The intron/exon structures are pair wise similar. Pp4CL2 and Pp4CL3 each consists of four exons and three introns, whereas Pp4CL1 and Pp4CL4 are characterized each by five exons and four introns. Pp4CL1, Pp4CL2 and Pp4CL3 are expressed in both gametophore and protonema tissue of P. patens, unlike Pp4CL4 whose expression could not be demonstrated under the conditions employed. Phylogenetic analysis suggests an early evolutionary divergence of Pp4CL gene family members. Using Streptomyces coelicolor cinnamate:CoA ligase (ScCCL) as an outgroup, the P. patens 4CLs are clearly separated from the spermatophyte proteins, but are intercalated between the angiosperm 4CL class I and class II. A comparison of three P. patens subspecies from diverse geographical locations shows high sequence identities for the four 4CL isoforms.

Riboswitch control of gene expression in plants by splicing and alternative 3' end processing of mRNAs

The most widespread riboswitch class, found in organisms from all three domains of life, is responsive to the vitamin B(1) derivative thiamin pyrophosphate (TPP). We have established that a TPP-sensing riboswitch is present in the 3' untranslated region (UTR) of the thiamin biosynthetic gene THIC of all plant species examined. The THIC TPP riboswitch controls the formation of transcripts with alternative 3' UTR lengths, which affect mRNA accumulation and protein production. We demonstrate that riboswitch-mediated regulation of alternative 3' end processing is critical for TPP-dependent feedback control of THIC expression. Our data reveal a mechanism whereby metabolite-dependent alteration of RNA folding controls splicing and alternative 3' end processing of mRNAs. These findings highlight the importance of metabolite sensing by riboswitches in plants and further reveal the significance of alternative 3' end processing as a mechanism of gene control in eukaryotes.

Erwinia carotovora elicitors and Botrytis cinerea activate defense responses in Physcomitrella patens

BACKGROUND

Vascular plants respond to pathogens by activating a diverse array of defense mechanisms. Studies with these plants have provided a wealth of information on pathogen recognition, signal transduction and the activation of defense responses. However, very little is known about the infection and defense responses of the bryophyte, Physcomitrella patens, to well-studied phytopathogens. The purpose of this study was to determine: i) whether two representative broad host range pathogens, Erwinia carotovora ssp. carotovora (E.c. carotovora) and Botrytis cinerea (B. cinerea), could infect Physcomitrella, and ii) whether B. cinerea, elicitors of a harpin (HrpN) producing E.c. carotovora strain (SCC1) or a HrpN-negative strain (SCC3193), could cause disease symptoms and induce defense responses in Physcomitrella.

RESULTS

B. cinerea and E.c. carotovora were found to readily infect Physcomitrella gametophytic tissues and cause disease symptoms. Treatments with B. cinerea spores or cell-free culture filtrates from E.c. carotovoraSCC1 (CF(SCC1)), resulted in disease development with severe maceration of Physcomitrella tissues, while CF(SCC3193) produced only mild maceration. Although increased cell death was observed with either the CFs or B. cinerea, the occurrence of cytoplasmic shrinkage was only visible in Evans blue stained protonemal cells treated with CF(SCC1) or inoculated with B. cinerea. Most cells showing cytoplasmic shrinkage accumulated autofluorescent compounds and brown chloroplasts were evident in a high proportion of these cells. CF treatments and B. cinerea inoculation induced the expression of the defense-related genes: PR-1, PAL, CHS and LOX.

CONCLUSION

B. cinerea and E.c. carotovora elicitors induce a defense response in Physcomitrella, as evidenced by enhanced expression of conserved plant defense-related genes. Since cytoplasmic shrinkage is the most common morphological change observed in plant PCD, and that harpins and B. cinerea induce this type of cell death in vascular plants, our results suggest that E.c. carotovora CFSCC1 containing HrpN and B. cinerea could also induce this type of cell death in Physcomitrella. Our studies thus establish Physcomitrella as an experimental host for investigation of plant-pathogen interactions and B. cinerea and elicitors of E.c. carotovora as promising tools for understanding the mechanisms involved in defense responses and in pathogen-mediated cell death in this simple land plant.

Sulfate assimilation in basal land plants - what does genomic sequencing tell us?

Sulfate assimilation is a pathway providing reduced sulfur for the synthesis of cysteine, methionine, co-enzymes such as iron-sulfur centres, thiamine, lipoic acid, or Coenzyme A, and many secondary metabolites, e.g., glucosinolates or alliins. The pathway is relatively well understood in flowering plants, but very little information exists on sulfate assimilation in basal land plants. Since the finding of a putative 3'-phosphoadenosine 5'-phosphosulfate reductase in PHYSCOMITRELLA PATENS, an enigmatic enzyme thought to exist in fungi and some bacteria only, it has been evident that sulfur metabolism in lower plants may substantially differ from seed plant models. The genomic sequencing of two basal plant species, the Bryophyte PHYSCOMITRELLA PATENS, and the Lycophyte SELAGINELLA MOELLENDORFFII, opens up the possibility to search for differences between lower and higher plants at the genomic level. Here we describe the similarities and differences in the organisation of the sulfate assimilation pathway between basal and advanced land plants derived from genome comparisons of these two species with ARABIDOPSIS THALIANA and ORYZA SATIVA, two seed plants with sequenced genomes. We found differences in the number of genes encoding sulfate transporters, adenosine 5'-phosphosulfate reductase, and sulfite reductase between the lower and higher plants. The consequences for regulation of the pathway and evolution of sulfate assimilation in plants are discussed.

Current achievements in the production of complex biopharmaceuticals with moss bioreactors

Transgenic plants are promising alternatives for the low-cost and safe pathogen-free production of complex recombinant pharmaceutical proteins (molecular farming). Plants as higher eukaryotes perform posttranslational modifications similar to those of mammalian cells. However, plant-specific protein N-glycosylation was shown to be immunogenic, a fact that represents a drawback for many plant systems in biopharmaceutical production. The moss Physcomitrella patens offers unique properties as a contained system for protein production. It is grown in the predominant haploid gametophytic stage as tissue suspension cultures in photobioreactors. Efficient secretory signals and a transient transfection system allow the secretion of freshly synthesized proteins to the surrounding medium. The key advantage of Physcomitrella compared to other plant systems is the feasibility of targeted gene replacements. By this means, moss strains with non-immunogenic humanized glycan patterns were created. Here we present an overview of the relevant aspects for establishing moss as a production system for recombinant biopharmaceuticals.

An ancient genome duplication contributed to the abundance of metabolic genes in the moss Physcomitrella patens

Background:

Analyses of complete genomes and large collections of gene transcripts have shown that most, if not all seed plants have undergone one or more genome duplications in their evolutionary past.

Results:

In this study, based on a large collection of EST sequences, we provide evidence that the haploid moss Physcomitrella patens is a paleopolyploid as well. Based on the construction of linearized phylogenetic trees we infer the genome duplication to have occurred between 30 and 60 million years ago. Gene Ontology and pathway association of the duplicated genes in P. patens reveal different biases of gene retention compared with seed plants.

Conclusion:

Metabolic genes seem to have been retained in excess following the genome duplication in P. patens. This might, at least partly, explain the versatility of metabolism, as described for P. patens and other mosses, in comparison to other land plants.

Comparative analysis of the SBP-box gene families in P. patens and seed plants

To come to a better understanding of the evolution and function of the SBP-box transcription factor family in plants, we identified, isolated and characterized 13 of its members from the moss Physcomitrella patens. For the majority of the moss SBP-box genes, clear orthologous relationships with family members of flowering plants could be established by phylogenetic analysis based on the conserved DNA-binding SBP-domain, as well as additional synapomorphic molecular characters. The P. patens SBP-box genes cluster in four separable groups. One of these consists exclusively of moss genes; the three others are shared with family members of Arabidopsis and rice. Besides the family defining DNA-binding SBP-domain, other features can be found conserved between moss and other plant SBP-domain proteins. An AHA-like motif conserved from the unicellular alga Chlamydomonas reinhardtii to flowering plants, was found able to promote transcription in a heterologous yeast system. The conservation of a functional microRNA response element in the mRNA of three of the moss SBP-box genes supports the idea of an ancient origin of microRNA dependent regulation of SBP-box gene family members. As our current knowledge concerning the roles of SBP-box genes in plant development is scarce and the model system P. patens allows targeted mutation, the material we isolated and characterized will be helpful to generate the mutant phenotypes necessary to further elucidate these roles.

Dating the early evolution of plants: detection and molecular clock analyses of orthologs

Orthologs generally are under selective pressure against loss of function, while paralogs usually accumulate mutations and finally die or deviate in terms of function or regulation. Most ortholog detection methods contaminate the resulting datasets with a substantial amount of paralogs. Therefore we aimed to implement a straightforward method that allows the detection of ortholog clusters with a reduced amount of paralogs from completely sequenced genomes. The described cross-species expansion of the reciprocal best BLAST hit method is a time-effective method for ortholog detection, which results in 68% truly orthologous clusters and the procedure specifically enriches single-copy orthologs. The detection of true orthologs can provide a phylogenetic toolkit to better understand evolutionary processes. In a study across six photosynthetic eukaryotes, nuclear genes of putative mitochondrial origin were shown to be over-represented among single copy orthologs. These orthologs are involved in fundamental biological processes like amino acid metabolism or translation. Molecular clock analyses based on this dataset yielded divergence time estimates for the red/green algae (1,142 MYA), green algae/land plant (725 MYA), mosses/seed plant (496 MYA), gymno-/angiosperm (385 MYA) and monocotyledons/core eudicotyledons (301 MYA) divergence times.

Toward a property/function relationship for metallothioneins: Histidine coordination and unusual cluster composition in a zinc-metallothionein from plants

Early cysteine labeled (E(C)) proteins are plant metallothioneins, which were first identified in wheat embryos and are thought to be seed-specific. An exhaustive analysis of expressed sequence tag (EST) entries reveals that homologs are expressed in embryos of both classes of flowering plants (monocotyledons and dicotyledons), but also occur in conifers (gymnosperms) and seed-free spike moss (lycophyta). Mass spectrometric and elemental analysis results indicate that, contrary to the widely propagated number of five, E(C) binds predominantly six zinc ions in at least two zinc-thiolate clusters. 1H and 111Cd NMR experiments suggest that, in contrast to the majority of previously characterized metallothioneins, two conserved histidine residues participate in metal binding. The collected data is consistent with the presence of clusters unprecedented in metallothioneins so far. This novel cluster composition is accompanied by metal-binding properties that are substantially different from other metallothioneins; thus wheat E(C) binds zinc less strongly than either mammalian or cyanobacterial MTs. The unique biochemical properties of wheat E(C) render it ideally suited for a role in zinc donation to nascent proteins during seed development, a role that has been suggested based on the fact that E(C) is induced by the plant hormone abscisic acid, but not by heavy metals. Our results provide a step further toward developing a property/function relationship for metallothioneins.

The putative moss 3'-phosphoadenosine-5'-phosphosulfate reductase is a novel form of adenosine-5'-phosphosulfate reductase without an iron-sulfur cluster

Sulfate assimilation provides reduced sulfur for synthesis of the amino acids cysteine and methionine and for a range of other metabolites. Sulfate has to be activated prior to reduction by adenylation to adenosine 5'-phosphosulfate (APS). In plants, algae, and many bacteria, this compound is reduced to sulfite by APS reductase (APR); in fungi and some cyanobacteria and gamma-proteobacteria, a second activation step, phosphorylation to 3'-phosphoadenosine 5'-phosphosulfate (PAPS), is necessary before reduction to sulfite by PAPS reductase (PAPR). We found previously that the moss Physcomitrella patens is unique among these organisms in possessing orthologs of both APR and PAPR genes (Koprivova, A., Meyer, A. J., Schween, G., Herschbach, C., Reski, R., and Kopriva, S. (2002) J. Biol. Chem. 277, 32195-32201). To assess the function of the two enzymes, we compared their biochemical properties by analysis of purified recombinant proteins. APR from Physcomitrella is very similar to the well characterized APRs from seed plants. On the other hand, we found that the putative PAPR preferentially reduces APS. Sequence analysis, analysis of UV-visible spectra, and determination of iron revealed that this new APR, named PpAPR-B, does not contain the FeS cluster, which was previously believed to determine the substrate specificity of the otherwise relatively similar enzymes. The lack of the FeS cluster in PpAPR-B catalysis is connected with a lower turnover rate but higher stability of the protein. These findings show that APS reduction without the FeS cluster is possible and that plant sulfate assimilation is predominantly dependent on reduction of APS.

PlanTAPDB, a phylogeny-based resource of plant transcription-associated proteins

Diversification of transcription-associated protein (TAP) families during land plant evolution is a key process yielding increased complexity of plant life. Understanding the evolutionary relationships between these genes is crucial to gain insight into plant evolution. We have determined a substantial set of TAPs that are focused on, but not limited to, land plants using PSI-BLAST searches and subsequent filtering and clustering steps. Phylogenies were created in an automated way using a combination of distance and maximum likelihood methods. Comparison of the data to previously published work confirmed their accuracy and usefulness for the majority of gene families. Evidence is presented that the flowering plant apical stem cell regulator WUSCHEL evolved from an ancestral homeobox gene that was already present after the water-to-land transition. The presence of distinct expanded gene families, such as COP1 and HIT in moss, is discussed within the evolutionary backdrop. Comparative analyses revealed that almost all angiosperm transcription factor families were already present in the earliest land plants, whereas many are missing among unicellular algae. A global analysis not only of transcription factors but also of transcriptional regulators and novel putative families is presented. A wealth of data about plant TAP families and all data accrued throughout their automated detection and analysis are made available via the PlanTAPDB Web interface. Evolutionary relationships of these genes are readily accessible to the nonexpert at a mouse-click. Initial analyses of selected gene families revealed that PlanTAPDB can easily be exerted for knowledge discovery.

Evidence for the rapid expansion of microRNA-mediated regulation in early land plant evolution

BACKGROUND

MicroRNAs (miRNAs) are regulatory RNA molecules that are specified by their mode of action, the structure of primary transcripts, and their typical size of 20-24 nucleotides. Frequently, not only single miRNAs but whole families of closely related miRNAs have been found in animals and plants. Some families are widely conserved among different plant taxa. Hence, it is evident that these conserved miRNAs are of ancient origin and indicate essential functions that have been preserved over long evolutionary time scales. In contrast, other miRNAs seem to be species-specific and consequently must possess very distinct functions. Thus, the analysis of an early-branching species provides a window into the early evolution of fundamental regulatory processes in plants.

RESULTS

Based on a combined experimental-computational approach, we report on the identification of 48 novel miRNAs and their putative targets in the moss Physcomitrella patens. From these, 18 miRNAs and two targets were verified in independent experiments. As a result of our study, the number of known miRNAs in Physcomitrella has been raised to 78. Functional assignments to mRNAs targeted by these miRNAs revealed a bias towards genes that are involved in regulation, cell wall biosynthesis and defense. Eight miRNAs were detected with different expression in protonema and gametophore tissue. The miRNAs 1-50 and 2-51 are located on a shared precursor that are separated by only one nucleotide and become processed in a tissue-specific way.

CONCLUSION

Our data provide evidence for a surprisingly diverse and complex miRNA population in Physcomitrella. Thus, the number and function of miRNAs must have significantly expanded during the evolution of early land plants. As we have described here within, the coupled maturation of two miRNAs from a shared precursor has not been previously identified in plants.

Heterologous expression of cDNAs encoding monodehydroascorbate reductases from the moss, Physcomitrella patens and characterization of the expressed enzymes

Monodehydroascorbate reductase (MDHAR; EC 1.6.5.4) catalyses the reduction of the monodehydroascorbate (MDHA) radical to ascorbate, using NADH or NADPH as an electron donor, and is believed to be involved in maintaining the reactive oxygen scavenging capability of plant cells. This key enzyme in the ascorbate-glutathione cycle has been studied here in the moss Physcomitrella patens, which is tolerant to a range of abiotic stresses and is increasingly used as a model plant. In the present study, three cDNAs encoding different MDHAR isoforms of 47 kDa were identified in P. patens, and found to exhibit enzymic characteristics similar to MDHARs in vascular plants despite low-sequence identity and a distant evolutionary relationship between the species. The three cDNAs for the P. patens MDHAR enzymes were expressed in Escherichia coli and the active enzymes were purified and characterized. Each recombinant protein displayed an absorbance spectrum typical of flavoenzymes and contained a single non-covalently bound FAD coenzyme molecule. The Km and kcat values for the heterologously expressed PpMDHAR enzymes ranged from 8 to 18 microM and 120-130 s(-1), respectively, using NADH as the electron donor. The Km values were at least an order of magnitude higher for NADPH. The Km values for the MDHA radical were approximately 0.5-1.0 microM for each of the purified enzymes, and further kinetic analyses indicated that PpMDHARs follow a 'ping-pong' kinetic mechanism. In contrast to previously published data, site-directed mutagenesis indicated that the conserved cysteine residue is not directly involved in the reduction of MDHA.

XET activity is found near sites of growth and cell elongation in bryophytes and some green algae: new insights into the evolution of primary cell wall elongation

BACKGROUND AND AIMS

In angiosperms xyloglucan endotransglucosylase (XET)/hydrolase (XTH) is involved in reorganization of the cell wall during growth and development. The location of oligo-xyloglucan transglucosylation activity and the presence of XTH expressed sequence tags (ESTs) in the earliest diverging extant plants, i.e. in bryophytes and algae, down to the Phaeophyta was examined. The results provide information on the presence of an XET growth mechanism in bryophytes and algae and contribute to the understanding of the evolution of cell wall elongation in general.

METHODS

Representatives of the different plant lineages were pressed onto an XET test paper and assayed. XET or XET-related activity was visualized as the incorporation of fluorescent signal. The Physcomitrella genome database was screened for the presence of XTHs. In addition, using the 3' RACE technique searches were made for the presence of possible XTH ESTs in the Charophyta.

KEY RESULTS

XET activity was found in the three major divisions of bryophytes at sites corresponding to growing regions. In the Physcomitrella genome two putative XTH-encoding cDNA sequences were identified that contain all domains crucial for XET activity. Furthermore, XET activity was located at the sites of growth in Chara (Charophyta) and Ulva (Chlorophyta) and a putative XTH ancestral enzyme in Chara was identified. No XET activity was identified in the Rhodophyta or Phaeophyta.

CONCLUSIONS

XET activity was shown to be present in all major groups of green plants. These data suggest that an XET-related growth mechanism originated before the evolutionary divergence of the Chlorobionta and open new insights in the evolution of the mechanisms of primary cell wall expansion.