Large-scale analysis of 73 329 physcomitrella plants transformed with different gene disruption libraries: production parameters and mutant phenotypes

Putative NAD(P)-Binding Rossmann Fold Protein Is Involved in Chitosan-Induced Peroxidase Activity and Lipoxygenase Expression in Physcomitrium patens

Oxidative burst, the rapid production of high levels of reactive oxygen species in response to external stimuli, is an early defense reaction against pathogens. The fungal elicitor chitosan causes an oxidative burst in the moss Physcomitrium patens (formerly Physcomitrella patens), mainly due to the peroxidase enzyme Prx34. To better understand the chitosan responses in P. patens, we conducted a screen of part of a P. patens mutant collection to isolate plants with less peroxidase activity than wild-type (WT) plants after chitosan treatment. We isolated a P. patens mutant that affected the gene encoding NAD(P)-binding Rossmann fold protein (hereafter, Rossmann fold protein). Three Rossmann fold protein-knockout (KO) plants (named Rossmann fold KO lines) were generated and used to assess extracellular peroxidase activity and expression of defense-responsive genes, including alternative oxidase, lipoxygenase (LOX), NADPH oxidase, and peroxidase (Prx34) in response to chitosan treatment. Extracellular (apoplastic) peroxidase activity was significantly lower in Rossmann fold KO lines than in WT plants after chitosan treatments. Expression of the LOX gene in Rossmann fold KO plants was significantly lower before and after chitosan treatment when compared with WT. Peroxidase activity assays together with gene expression analyses suggest that the Rossmann fold protein might be an important component of the signaling pathway leading to oxidative burst and basal expression of the LOX gene in P. patens. [Formula: see text] Copyright © 2023 The Author(s). This is an open access article distributed under the CC BY-NC-ND 4.0 International license.

Autopolyploidization affects transcript patterns and gene targeting frequencies in Physcomitrella

In Physcomitrella, whole-genome duplications affected the expression of about 3.7% of the protein-encoding genes, some of them relevant for DNA repair, resulting in a massively reduced gene-targeting frequency. Qualitative changes in gene expression after an autopolyploidization event, a pure duplication of the whole genome (WGD), might be relevant for a different regulation of molecular mechanisms between angiosperms growing in a life cycle with a dominant diploid sporophytic stage and the haploid-dominant mosses. Whereas angiosperms repair DNA double-strand breaks (DSB) preferentially via non-homologous end joining (NHEJ), in the moss Physcomitrella homologous recombination (HR) is the main DNA-DSB repair pathway. HR facilitates the precise integration of foreign DNA into the genome via gene targeting (GT). Here, we studied the influence of ploidy on gene expression patterns and GT efficiency in Physcomitrella using haploid plants and autodiploid plants, generated via an artificial WGD. Single cells (protoplasts) were transfected with a GT construct and material from different time-points after transfection was analysed by microarrays and SuperSAGE sequencing. In the SuperSAGE data, we detected 3.7% of the Physcomitrella genes as differentially expressed in response to the WGD event. Among the differentially expressed genes involved in DNA-DSB repair was an upregulated gene encoding the X-ray repair cross-complementing protein 4 (XRCC4), a key player in NHEJ. Analysing the GT efficiency, we observed that autodiploid plants were significantly GT suppressed (p < 0.001) attaining only one third of the expected GT rates. Hence, an alteration of global transcript patterns, including genes related to DNA repair, in autodiploid Physcomitrella plants correlated with a drastic suppression of HR.

Chloroplasts require glutathione reductase to balance reactive oxygen species and maintain efficient photosynthesis

Thiol-based redox-regulation is vital for coordinating chloroplast functions depending on illumination and has been throroughly investigated for thioredoxin-dependent processes. In parallel, glutathione reductase (GR) maintains a highly reduced glutathione pool, enabling glutathione-mediated redox buffering. Yet, how the redox cascades of the thioredoxin and glutathione redox machineries integrate metabolic regulation and detoxification of reactive oxygen species remains largely unresolved because null mutants of plastid/mitochondrial GR are embryo-lethal in Arabidopsis thaliana. To investigate whether maintaining a highly reducing stromal glutathione redox potential (E_GSH ) via GR is necessary for functional photosynthesis and plant growth, we created knockout lines of the homologous enzyme in the model moss Physcomitrella patens. In these viable mutant lines, we found decreasing photosynthetic performance and plant growth with increasing light intensities, whereas ascorbate and zeaxanthin/antheraxanthin levels were elevated. By in vivo monitoring stromal E_GSH dynamics, we show that stromal E_GSH is highly reducing in wild-type and clearly responsive to light, whereas an absence of GR leads to a partial glutathione oxidation, which is not rescued by light. By metabolic labelling, we reveal changing protein abundances in the GR knockout plants, pinpointing the adjustment of chloroplast proteostasis and the induction of plastid protein repair and degradation machineries. Our results indicate that the plastid thioredoxin system is not a functional backup for the plastid glutathione redox systems, whereas GR plays a critical role in maintaining efficient photosynthesis.

Physcomitrella patens, a versatile synthetic biology chassis

During three decades the moss Physcomitrella patens has been developed to a superb green cell factory with the first commercial products on the market. In the past three decades the moss P. patens has been developed from an obscure bryophyte to a model organism in basic biology, biotechnology, and synthetic biology. Some of the key features of this system include a wide range of Omics technologies, precise genome-engineering via homologous recombination with yeast-like efficiency, a certified good-manufacturing-practice production in bioreactors, successful upscaling to 500 L wave reactors, excellent homogeneity of protein products, superb product stability from batch-to-batch, and a reliable procedure for cryopreservation of cell lines in a master cell bank. About a dozen human proteins are being produced in P. patens as potential biopharmaceuticals, some of them are not only similar to their animal-produced counterparts, but are real biobetters with superior performance. A moss-made pharmaceutical successfully passed phase 1 clinical trials, a fragrant moss, and a cosmetic moss-product is already on the market, highlighting the economic potential of this synthetic biology chassis. Here, we focus on the features of mosses as versatile cell factories for synthetic biology and their impact on metabolic engineering.

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

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

Highly efficient gene tagging in the bryophyte Physcomitrella patens using the tobacco (Nicotiana tabacum) Tnt1 retrotransposon

Because of its highly efficient homologous recombination, the moss Physcomitrella patens is a model organism particularly suited for reverse genetics, but this inherent characteristic limits forward genetic approaches. Here, we show that the tobacco (Nicotiana tabacum) retrotransposon Tnt1 efficiently transposes in P. patens, being the first retrotransposon from a vascular plant reported to transpose in a bryophyte. Tnt1 has a remarkable preference for insertion into genic regions, which makes it particularly suited for gene mutation. In order to stabilize Tnt1 insertions and make it easier to select for insertional mutants, we have developed a two-component system where a mini-Tnt1 with a retrotransposition selectable marker can only transpose when Tnt1 proteins are co-expressed from a separate expression unit. We present a new tool with which to produce insertional mutants in P. patens in a rapid and straightforward manner that complements the existing molecular and genetic toolkit for this model species.

A single homeobox gene triggers phase transition, embryogenesis and asexual reproduction

Plants characteristically alternate between haploid gametophytic and diploid sporophytic stages. Meiosis and fertilization respectively initiate these two different ontogenies(1). Genes triggering ectopic embryo development on vegetative sporophytic tissues are well described(2,3); however, a genetic control of embryo development from gametophytic tissues remains elusive. Here, in the moss Physcomitrella patens we show that ectopic overexpression of the homeobox gene BELL1 induces embryo formation and subsequently reproductive diploid sporophytes from specific gametophytic cells without fertilization. In line with this, BELL1 loss-of-function mutants have a wild-type phenotype, except that their egg cells are bigger and unable to form embryos. Our results identify BELL1 as a master regulator for the gametophyte-to-sporophyte transition in P. patens and provide mechanistic insights into the evolution of embryos that can generate multicellular diploid sporophytes. This developmental innovation facilitated the colonization of land by plants about 500 million years ago(4) and thus shaped our current ecosystems.

Moss-made pharmaceuticals: from bench to bedside

Over the past two decades, the moss Physcomitrella patens has been developed from scratch to a model species in basic research and in biotechnology. A fully sequenced genome, outstanding possibilities for precise genome-engineering via homologous recombination (knockout moss), a certified GMP production in moss bioreactors, successful upscaling to 500 L wave reactors, excellent homogeneity of protein glycosylation, remarkable batch-to-batch stability and a safe cryopreservation for master cell banking are some of the key features of the moss system. Several human proteins are being produced in this system as potential biopharmaceuticals. Among the products are tumour-directed monoclonal antibodies with enhanced antibody-dependent cytotoxicity (ADCC), vascular endothelial growth factor (VEGF), complement factor H (FH), keratinocyte growth factor (FGF7/KGF), epidermal growth factor (EGF), hepatocyte growth factor (HGF), asialo-erythropoietin (asialo-EPO, AEPO), alpha-galactosidase (aGal) and beta-glucocerebrosidase (GBA). Further, an Env-derived multi-epitope HIV protein as a candidate vaccine was produced, and first steps for a metabolic engineering of P. patens have been made. Some of the recombinant biopharmaceuticals from moss bioreactors are not only similar to those produced in mammalian systems such as CHO cells, but are of superior quality (biobetters). The first moss-made pharmaceutical, aGal to treat Morbus Fabry, is in clinical trials.

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.

Recombination products suggest the frequent occurrence of aberrant gene replacement in the moss Physcomitrella patens

In gene replacement, a variant of gene targeting, transformed DNA integrates into the genome by homologous recombination (HR) to replace resident sequences. Gene replacement in the moss Physcomitrella patens is extremely efficient, but often large amounts of additional DNA are integrated at the target locus. A detailed analysis of recombination junctions of PpCOL2 gene knockout mutants shows that the integrated DNA can be highly rearranged. Our data suggest that the replaced sequences were excised by HR and became integrated back into the genome by non-homologous end-joining (NHEJ). RAD51-mediated strand-invasion and subsequent strand-exchange is central to the two-end invasion pathway, the major gene replacement pathway in yeast. In this pathway, integration is initiated by the free ends of a single replacement vector-derived donor molecule which then integrates as an entity. Gene replacement in P. patens is entirely RAD51-dependent suggesting the existence of a pathway mechanistically similar to two-end invasion. However, invasion of the two ends does not seem to be stringently coordinated in P. patens. Actually, often only one fragment end became integrated by HR, or one-sided integration of two independent donor fragments occurred simultaneously leading to a double-strand break that is subsequently sealed by NHEJ and thus causes the observed rearrangements.

Architecture of the PPR gene family in the moss Physcomitrella patens

Pentatricopeptide repeat (PPR) proteins are widespread in eukaryotes and in particular, include several hundred members in land plants. The majority of PPR proteins are localized in mitochondria and plastids, where they play a crucial role in various aspects of RNA metabolism at the post-transcriptional level in gene expression. However, many of their functions remain to be characterized. In contrast to vascular plants, the moss Physcomitrella patens has only 105 PPR genes. This number may represent a minimum set of PPR proteins required for post-transcriptional regulation in plant organelles. Here, we review the overall structure of the P. patens PPR gene family and the current status of the functional characterization of moss PPR proteins.

Usefulness of Physcomitrella patens for studying plant organogenesis

In this chapter, we review the main organogenesis features and associated regulation processes of the moss Physcomitrella patens (P. patens), the model plant for the Bryophytes. We highlight how the study of this descendant of the earliest plant species that colonized earth, brings useful keys to understand the mechanisms that determine and control both vascular and non vascular plants organogenesis. Despite its simple morphogenesis pattern, P. patens still requires the fine tuning of organogenesis regulators, including hormone signalling, common to the whole plant kingdom, and which study is facilitated by a high number of molecular tools, among which the powerful possibility of gene targeting/replacement. The recent discovery of moss cells reprogramming capacity completes the picture of an excellent model for studying plant organogenesis.

Novel antimalarial drug targets: hope for new antimalarial drugs

Malaria is a major global threat, that results in more than 2 million deaths each year. The treatment of malaria is becoming extremely difficult due to the emergence of drug-resistant parasites, the absence of an effective vaccine, and the spread of insecticide-resistant vectors. Thus, malarial therapy needs new chemotherapeutic approaches leading to the search for new drug targets. Here, we discuss different approaches to identifying novel antimalarial drug targets. We have also given due attention to the existing validated targets with a view to develop novel, rationally designed lead molecules. Some of the important parasite proteins are claimed to be the targets; however, further in vitro or in vivo structure-function studies of such proteins are crucial to validate these proteins as suitable targets. The interactome analysis among apicoplast, mitochondrion and genomic DNA will also be useful in identifying vital pathways or proteins regulating critical pathways for parasite growth and survival, and could be attractive targets. Molecules responsible for parasite invasion to host erythrocytes and ion channels of infected erythrocytes, essential for intra-erythrocyte survival and stage progression of parasites are also becoming attractive targets. This review will discuss and highlight the current understanding regarding the potential antimalarial drug targets, which could be utilized to develop novel antimalarials.

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.

Large-scale phenotyping of transgenic tobacco plants (Nicotiana tabacum) to identify essential leaf functions

Two of the major challenges in functional genomics are to identify genes that play a key role in biological processes, and to elucidate the biological role of the large numbers of genes whose function is poorly characterized or still completely unknown. In this study, a combination of large-scale expressed sequence tag sequencing, high-throughput gene silencing and visual phenotyping was used to identify genes in which partial inhibition of expression leads to marked phenotypic changes, mostly on leaves. Three normalized tobacco (Nicotiana tabacum) cDNA libraries were prepared directly in a binary vector using different tissues of tobacco as an RNA source, randomly sequenced and clustered. The Agrobacterium-tobacco leaf disc transformation system was used to generate sets of antisense or co-suppression transgenic tobacco plants for over 20 000 randomly chosen clones, each representing an independent cluster. After transfer to the glasshouse, transgenic plants were scored visually after 10-14 days for changes in growth, leaf form and chlorosis or necrosis. Putative hits were validated by repeating the transformation. This procedure is more stringent than the analysis of knockout mutants, because it requires that even a partial decrease in expression generates a phenotype. This procedure identified 88 validated gene/phenotype relations. These included several previously characterized gene/phenotype relationships, demonstrating the validity of the approach. For about one-third, a function could be inferred, but a loss-of-function phenotype had not been described previously. Strikingly, almost one-half of the validated genes were poorly annotated, or had no known function. For 77 of these tobacco sequences, a single or small number of potential orthologues were identified in Arabidopsis. The genes for which orthologues were identified in Arabidopsis included about one-half of the genes whose function was completely unknown. Comparison with published gene/phenotype relations for Arabidopsis knockout mutants revealed surprisingly little overlap with the present study. Our results indicate that partial gene silencing identifies novel gene/phenotype relationships, which are distinct from those uncovered by knockout screens. They also show that it is possible to perform these analyses in a crop species in which full genome sequence information is lacking, and subsequently to transfer the information to a reference species in which functional studies can be performed more effectively.

Differential requirements for RAD51 in Physcomitrella patens and Arabidopsis thaliana development and DNA damage repair

RAD51, the eukaryotic homolog of the bacterial RecA recombinase, plays a central role in homologous recombination (HR) in yeast and animals. Loss of RAD51 function causes lethality in vertebrates but not in other animals or in the flowering plant Arabidopsis thaliana, suggesting that RAD51 is vital for highly developed organisms but not for others. Here, we found that loss of RAD51 function in the moss Physcomitrella patens, a plant of less complexity, caused a significant vegetative phenotype, indicating an important function for RAD51 in this organism. Moreover, loss of RAD51 caused marked hypersensitivity to the double-strand break-inducing agent bleomycin in P. patens but not in Arabidopsis. Therefore, HR is used for somatic DNA damage repair in P. patens but not in Arabidopsis. These data imply fundamental differences in the use of recombination pathways between plants. Moreover, these data demonstrate that the importance of RAD51 for viability is independent of taxonomic position or complexity of an organism. The involvement of HR in DNA damage repair in the slowly evolving species P. patens but not in fast-evolving Arabidopsis suggests that the choice of the recombination pathway is related to the speed of evolution in plants.

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.

The mechanism of gene targeting in Physcomitrella patens: homologous recombination, concatenation and multiple integration

The model bryophyte Physcomitrella patens exhibits high frequencies of gene targeting when transformed with DNA constructs containing sequences homologous with genomic loci. ‘Targeted gene replacement’ (TGR) resulting from homologous recombination (HR) between each end of a targeting construct and the targeted locus occurs when either single or multiple targeting vectors are delivered. In the latter instance simultaneous, multiple, independent integration of different transgenes occurs at the targeted loci. In both single gene and ‘batch’ transformations, DNA can also be found to undergo ‘targeted insertion’ (TI), integrating at one end of the targeted locus by HR with one flanking sequence of the vector accompanied by an apparent non-homologous end-joining (NHEJ) event at the other. Untargeted integration at nonhomologous sites also occurs, but at a lower frequency. Molecular analysis of TI at a single locus shows that this occurs as a consequence of concatenation of the transforming DNA, in planta, prior to integration, followed by HR between a single site in the genomic target and two of its repeated homologues in the concatenated vector. This reinforces the view that HR is the major pathway by which transforming DNA is integrated in Physcomitrella.

Moss systems biology en route: phytohormones in Physcomitrella development

The moss Physcomitrella patens has become a powerful model system in modern plant biology. Highly standardized cell culture techniques, as well as the necessary tools for computational biology, functional genomics and proteomics have been established. Large EST collections are available and the complete moss genome will be released soon. A simple body plan and the small number of different cell types in Physcomitrella facilitate the study of developmental processes. In the filamentous juvenile moss tissue, developmental decisions rely on the differentiation of single cells. Developmental steps are controlled by distinct phytohormones and integration of environmental signals. Especially the phytohormones auxin, cytokinin, and abscisic acid have distinct effects on early moss development. In this article, we review current knowledge about phytohormone influences on early moss development in an attempt to fully unravel the complex regulatory signal transduction networks underlying the developmental decisions of single plant cells in a holistic systems biology approach.

Mosses as model systems for the study of metabolism and development

The haploid gametophyte stage of the moss life cycle is amenable to genetic and biochemical studies. Many species can be cultured on simple defined media, where growth is rapid, making them ideal material for metabolic studies. Developmental responses to hormones and to environmental inputs can be studied both at the level of individual cells and in multicellular tissues. The protonemal stage of gametophyte development comprises cell filaments that extend by the serial division of their apical cells, allowing the investigation of the generation and modification of cell polarity and the role of the cytoskeleton in these processes. Molecular techniques including gene inactivation by targeted gene replacement or by RNA interference, together with the nearly completed sequencing of the Physcomitrella patens genome, open the way for detailed study of the functions of genes involved in both development and metabolism.

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

To gain insight into the transcriptome of the well-used plant model system Physcomitrella patens, several EST sequencing projects have been undertaken. We have clustered, assembled, and annotated all publicly available EST and CDS sequences in order to represent the transcriptome of this non-seed plant. Here, we present our fully annotated knowledge resource for the Physcomitrella patens transcriptome, integrating annotation from the production process of the clustered sequences and from a high-quality annotation pipeline developed during this study. Each transcript is represented as an entity containing full annotations and GO term associations. The whole production, filtering, clustering, and annotation process is being modelled and results in seven datasets, representing the annotated Physcomitrella transcriptome from different perspectives. We were able to annotate 63.4 % of the 26 123 virtual transcripts. The transcript archetype, as covered by our clustered data, is compared to a compilation based on all available Physcomitrella full length CDS. The distribution of the gene ontology annotations (GOA) for the virtual transcriptome of Physcomitrella patens demonstrates consistency in the ratios of the core molecular functions among the plant GOA. However, the metabolism subcategory is over-represented in bryophytes as compared to seed plants. This observation can be taken as an indicator for the wealth of alternative metabolic pathways in moss in comparison to spermatophytes. All resources presented in this study have been made available to the scientific community through a suite of user-friendly web interfaces via www.cosmoss.org and form the basis for assembly and annotation of the moss genome, which will be sequenced in 2005.