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.

MicroRNAs (miRNAs) control gene expression in animals and plants. Like another class of small RNAs, siRNAs, they affect gene expression posttranscriptionally. While siRNAs in addition act in transcriptional gene silencing, a role of miRNAs in transcriptional regulation has been less clear. We show here that in moss Physcomitrella patens mutants without a DICER-LIKE1b gene, maturation of miRNAs is normal but cleavage of target RNAs is abolished and levels of these transcripts are drastically reduced. These mutants accumulate miRNA:target-RNA duplexes and show hypermethylation of the genes encoding target RNAs, leading to gene silencing. This pathway occurs also in the wild-type upon hormone treatment. We propose that initiation of epigenetic silencing by DNA methylation depends on the ratio of the miRNA and its target RNA.

Little is known about the division of eukaryotic cell organelles and up to now neither in animals nor in plants has a gene product been shown to mediate this process. A cDNA encoding a homolog of the bacterial cell division protein FtsZ, an ancestral tubulin, was isolated from the eukaryote Physcomitrella patens and used to disrupt efficiently the genomic locus in this terrestrial seedless plant. Seven out of 51 transgenics obtained were knockout plants generated by homologous recombination; they were specifically impeded in plastid division with no detectable effect on mitochondrial division or plant morphology. Implications on the theory of endosymbiosis and on the use of reverse genetics in plants are discussed.

SummaryThe degradation of aromatic compounds follows different biochemical principles in aerobic and anaerobic microorganisms. While aerobes dearomatize and cleave the aromatic ring by oxygenases, facultative anaerobes utilize an ATP-dependent ring reductase for the dearomatization of the activated key intermediate benzoyl-coenzyme A (CoA). In this work, the aromatic metabolism was studied in the obligately anaerobic model organism Geobacter metallireducens . The gene coding for a putative carboxylic acid-CoA ligase was heterologously overexpressed and the gene product was characterized as a highly specific benzoate-CoA ligase catalysing the initial step of benzoate metabolism. However, no evidence for the presence of an ATP-dependent benzoyl-CoA reductase as observed in facultative anaerobes was obtained. In a proteomic approach benzoate-induced proteins were identified; the corresponding genes are organized in two clusters comprising 44 genes. Induction of representative genes during growth on benzoate was confirmed by reverse transcription polymerase chain reaction. The results obtained suggest that benzoate is activated to benzoyl-CoA, which is then reductively dearomatized to cyclohexa-1,5-diene-1-carbonyl-CoA, followed by β β β β -oxidation reactions to acetyl-CoA units, as in facultatively anaerobic bacteria. However, in G. metallireducens the process of reductive benzene ring dearomatization appears to be catalysed by a set of completely different protein components comprising putative molybdenum and selenocysteine containing enzymes.

SummaryBackgroundPlant body plans arise by the activity of meristematic growing tips during development and radiated independently in the gametophyte (n) and sporophyte (2n) stages of the life cycle during evolution. Although auxin and its intercellular transport by PIN family efflux carriers are primary regulators of sporophytic shoot development in flowering plants, the extent of conservation in PIN function within the land plants and the mechanisms regulating bryophyte gametophytic shoot development are largely unknown.ResultsWe have found that treating gametophytic shoots of the moss Physcomitrella patens with exogenous auxins and auxin transport inhibitors disrupts apical function and leaf development. Two plasma membrane-targeted PIN proteins are expressed in leafy shoots, and pin mutants resemble plants treated with auxins or auxin transport inhibitors. PIN-mediated auxin transport regulates apical cell function, leaf initiation, leaf shape, and shoot tropisms in moss gametophytes. pin mutant sporophytes are sometimes branched, reproducing a phenotype only previously seen in the fossil record and in rare natural moss variants.ConclusionsOur results show that PIN-mediated auxin transport is an ancient, conserved regulator of shoot development.

Evolutionary retention of duplicated genes encoding transcription-associated proteins (TAPs, comprising transcription factors and other transcriptional regulators) has been hypothesized to be positively correlated with increasing morphological complexity and paleopolyploidizations, especially within the plant kingdom. Here, we present the most comprehensive set of classification rules for TAPs and its application for genome-wide analyses of plants and algae. Using a dated species tree and phylogenetic comparative (PC) analyses, we define the timeline of TAP loss, gain, and expansion among Viridiplantae and find that two major bursts of gain/expansion occurred, coinciding with the water-to-land transition and the radiation of flowering plants. For the first time, we provide PC proof for the long-standing hypothesis that TAPs are major driving forces behind the evolution of morphological complexity, the latter in Plantae being shaped significantly by polyploidization and subsequent biased paleolog retention. Principal component analysis incorporating the number of TAPs per genome provides an alternate and significant proxy for complexity, ideally suited for PC genomics. Our work lays the ground for further interrogation of the shaping of gene regulatory networks underlying the evolution of organism complexity.

BackgroundThe 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.ResultsHere, 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.ConclusionsComparative 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.