We report the draft genome sequence of the model moss Physcomitrella patens and compare its features with those of flowering plants, from which it is separated by more than 400 million years, and unicellular aquatic algae. This comparison reveals genomic changes concomitant with the evolutionary movement to land, including a general increase in gene family complexity; loss of genes associated with aquatic environments (e.g., flagellar arms); acquisition of genes for tolerating terrestrial stresses (e.g., variation in temperature and water availability); and the development of the auxin and abscisic acid signaling pathways for coordinating multicellular growth and dehydration response. The Physcomitrella genome provides a resource for phylogenetic inferences about gene function and for experimental analysis of plant processes through this plant's unique facility for reverse genetics.

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.

We show in bryophytes that abscisic acid (ABA) pretreatment of moss (Physcomitrella patens) cells confers desiccation tolerance. In angiosperms, both ABA and the transcriptional regulator ABSCISIC ACID INSENSITIVE 3 (ABI3) are required to protect the seed during desiccation. ABA was not able to protect moss cells in stable deletion lines of ABI3 (DeltaPpabi3). Hence, moss has the same functional link between ABA, ABI3, and the desiccation tolerance phenotype that is found in angiosperms. Furthermore, we identified 22 genes that were induced during ABA pretreatment in wild-type lines. When their expression was compared with that of DeltaPpabi3 during ABA pretreatment and immediately after desiccation, a new target of ABI3 action appears to be in the recovery period.

INTRODUCTIONThe moss Physcomitrella patens has been used as an experimental organism for more than 80 years. Within the last 15 years, its use as a model to explore plant functions has increased enormously. The ability to use gene targeting and RNA interference methods to study gene function, the availability of many tools for comparative and functional genomics (including a sequenced and assembled genome, physical and genetic maps, and more than 250,000 expressed sequence tags [ESTs]), and a dominant haploid phase that allows direct forward genetic analysis have all led to a surge of new activity. P. patens can be easily cultured and spends the majority of its life cycle in the haploid state, allowing the application of experimental techniques similar to those used in microbes and yeast. Its development is relatively simple, and it generates only a few tissues that contain a limited number of cell types. Although mosses lack vascular tissue, true roots/stems/leaves, and flowers and seeds, many signaling pathways found in angiosperms are intact in moss. For example, the phytohormones auxin, cytokinin, and abscisic acid, as well as the photomorphogenic pigments phytochrome and cryptochrome, are all interwoven into distinct but overlapping pathways and linked to clear developmental phenotypes. In addition, about one-quarter of the moss genome contains genes with no known function based on sequence motifs, raising the likelihood of successful discovery efforts to identify new and novel gene functions.

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

Summary The sequencing of the Physcomitrella patens genome, combined with the high frequency of gene targeting in this species, makes it ideal for reverse genetic studies. For forward genetic studies, experimental crosses and genetic analysis of progeny are essential.Since P. patens is monoicous, producing both male and female gametes on the same gametophore, and undergoing self-fertilization at a high frequency, the identification of crossed sporophytes is difficult. Usually spores from many sporophytes from a mixed culture must be testsed for the production of recombinant progeny.Here, we describe the use of transgenic lines that express a fluorescent transgene constitutively, to provide a direct visual screen for hybrid sporophytes.We show that segregations in crosses obtained with this technique are as expected, and demonstrate its utility for the study of the rate of outcrossing between three isolates of P. patens.

Orientation of cell division is critical for plant morphogenesis. This is evident in the formation and function of meristems and for morphogenetic transitions. Mosses undergo such transitions: from two-dimensional tip-growing filaments (protonema) to the generation of three-dimensional leaf-like structures (gametophores).The Defective Kernel 1 (DEK1) protein plays a key role in the perception of and/or response to positional cues that specify the formation and function of the epidermal layer in developing seeds of flowering plants. The moss Physcomitrella patens contains the highly conserved DEK1 gene.Using efficient gene targeting, we generated a precise PpDEK1 deletion (Δdek1), which resulted in normal filamentous growth of protonema. Two distinct mutant phenotypes were observed: an excess of buds on the protonema, and abnormal cell divisions in the emerging buds resulting in developmental arrest and the absence of three-dimensional growth. Overexpression of a complete PpDEK1 cDNA, or the calpain domain of PpDEK1 alone, successfully complements both phenotypes.These results in P. patens demonstrate the morphogenetic importance of the DEK1 protein in the control of oriented cell divisions. As it is not for protonema, it will allow dissection of the structure/function relationships of the different domains of DEK1 using gene targeting in null mutant background.

RNAi is a powerful method for generating loss of function mutants, especially for targeting genes belonging to large gene families. We have recently shown that RNAi functions in the moss Physcomitrella patens. We obtained stable lines that show constitutive silencing of a nuclearly localized GFP:GUS fusion protein (NLS:GFP:GUS). However, lines that display silencing of the protein do not necessarily have reduced transcript levels. Therefore, a system has been developed that silences the NLS:GFP:GUS reporter construct at the same time as it silences a gene of interest. Gateway (Invitrogen) recombination cassettes were incorporated into these vectors to facilitate cloning of many different cDNA sequences. In addition, vectors were generated that contain genomic moss DNA sequence information to increase the production of stable moss lines. Transformation with these constructs results in strong silencing within 24 h and is stable for at least a month after transformation. FtsZ2-1, whose loss of function phenotype is known, was incorporated as a test case for analyzing phenotypes. One hundred per cent of regenerating colonies that have silenced GFP exhibit a loss of function FtsZ2-1 phenotype, validating the use of this system to assay phenotypes for plant genes of unknown function.