Email updates

Keep up to date with the latest news and content from BMC Genomics and BioMed Central.

Open Access Highly Accessed Research article

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

Andreas D Zimmer1, Daniel Lang1, Karol Buchta1, Stephane Rombauts89, Tomoaki Nishiyama5, Mitsuyasu Hasebe106, Yves Van de Peer89, Stefan A Rensing237 and Ralf Reski134*

Author Affiliations

1 Plant Biotechnology, Faculty of Biology, University of Freiburg, Schaenzlestrasse 1, 79104, Freiburg, Germany

2 Institute of Biology II, Faculty of Biology, University of Freiburg, Schaenzlestrasse 1, 79104, Freiburg, Germany

3 BIOSS - Centre for Biological Signalling Studies, Freiburg, Germany

4 FRIAS - Freiburg Institute for Advanced Studies, Freiburg, Germany

5 Advanced Science Research Center, Kanazawa University, Kanazawa, Japan

6 National Institute for Basic Biology, Okazaki, Japan

7 Faculty of Biology, University of Marburg, Karl-von-Frisch-Str. 8, 35043, Marburg, Germany

8 Department of Plant Systems Biology, VIB, Technologiepark 927, B-9052, Ghent, Belgium

9 Department of Plant Biotechnology and Bioinformatics, Ghent University, Technologiepark 927, B-9052, Ghent, Belgium

10 School of life Science, The Graduate University for Advanced Studies, Okazaki 444-8585, Japan

For all author emails, please log on.

BMC Genomics 2013, 14:498  doi:10.1186/1471-2164-14-498


The electronic version of this article is the complete one and can be found online at: http://www.biomedcentral.com/1471-2164/14/498


Received:6 March 2013
Accepted:19 July 2013
Published:23 July 2013

© 2013 Zimmer et al.; licensee BioMed Central Ltd.

This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract

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 webcite 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.

Keywords:
Bryophyte; Physcomitrella patens; Genome annotation; Gene structure; Reference genome; Model organism; UTR; Plant evolution; Non-flowering plant; Orphan genes

Background

Given its phylogenetic key position as an early diverging land plant that bridges the gap of about one billion years between the unicellular green algae and flowering plants, the moss Physcomitrella patens (Physcomitrella) unites most of the attributes desirable for a model organism, including a short generation time, small stature, comparatively low morphological complexity, a haplo-dominant life cycle, traceable cell lineage, high growth rate and simplicity of genetic transformation. Combined with the potential for evolutionary-developmental (evo-devo) studies, these traits have become increasingly attractive to a wide range of plant scientists. Over the last two decades, a growing community has established P. patens as a model organism with a well-developed molecular toolbox including the uniquely efficient gene targeting via homologous recombination and comprehensive genomics resources which have been made available early on using the central web service cosmoss.org [1-4]. The moss is also a promising model for green biotechnology [5-9], which allows the production of safe recombinant proteins with eukaryotic post-translational modifications in competitive quantities.

The draft genome sequence of the moss Physcomitrella patens was published in 2008 [10]. The availability of the genomic sequence and the established molecular toolbox provide the ideal foundation for extensive comparative and evo-devo analyses studies. This is reflected in the publication record - a growing body of researchers from all fields has begun to apply Physcomitrella as an additional model organism for comparative studies [11-20]. Various evo-devo studies have demonstrated the ability of Physcomitrella transgenes to act as functional orthologs in cross-species complementation assays using A. thaliana mutant lines (e.g. [21-24]). Additionally, large-scale analyses and cross-kingdom comparisons increasingly utilize the moss as a representative organism for the plant kingdom [25-29]. This ongoing interest, the available resources, the active community, and the moss’ attractive phylogenetic position recently led the U.S. Department of Energy’s Joint Genome Institute (JGI) to select P. patens as a “plant flagship genome” [30].

The quality of a genome annotation is the bottleneck for any form of downstream and comparative analyses. Particularly affected by flaws are the large-scale, high-throughput approaches employed in systems biology [31]. Following its initial V1.0 annotation, the P. patens genome annotation has been iteratively improved. The draft V1 assembly was based on whole-genome shotgun Sanger sequencing at 8.6x clone depth and comprises 2536 V1 scaffolds. This number was reduced to 2106 in the released V1.1 after removal of bacterial contaminations [10]. After an additional round of scaffold filtering, released as V1.2, the genome sequence of the 27 chromosomes is still scattered over 1995 genomic scaffolds [3]. The filtering of the gene catalogue, in particular by removal of transposable elements and other non-protein coding regions, led to the prediction of 27,966 protein-coding genes [3]. The results of a survey conducted on behalf of the JGI at the annual moss meeting 2011 clearly show that for most groups (79%) the initial sequencing and release of the genome was already “essential”. However, many research topics were listed that “would be enabled if a highly complete and accurate reference genome for Physcomitrella was available”, revealing the need for advanced genome annotation. For example, there were cases of well-characterized moss genes that were present in the genomic sequence but were missing from the gene catalogue [32]. In V1.2 only 4515 (~16%) gene models had both 5’ and 3’-UTRs (untranslated regions). Over 23,000 genes missed either 5’-UTR or 3’-UTR annotation and thus were incomplete. Furthermore, functional annotation was only available for 41% of the genes and hardly any of these annotations were backed by traceable experimental evidence. A further shortcoming was that no established and universal means for scientists and curators existed to link (published) knowledge on moss genes to the digital representations scattered across several databases.

Model organism databases (MODs, e.g. Gramene [33], TAIR [34], FlyBase [35]) are integrated, web-accessible resources, which are prerequisite for the success and the quality of a reference genome [36]. They act as general repositories for all kinds of research data and scientific knowledge that is generated by the scientific community. Thus, they provide the necessary infrastructure for researchers working with model species, are the focal point for scientists new or outside the field and provide conceptual interfaces for data exchange with more general data repositories (e.g. NCBI, UniProt [37], Ensembl [38], Phytozome [39], and PLAZA [40]) to enable comparative analyses and ensure overall data quality [41,42].

Experience from various MODs and whole genome sequencing projects shows that automatic annotation without substantial manual curation is not sufficient to ensure data quality and knowledge discovery [43]. An active community is necessary to transfer all available data, especially the biological knowledge covered by the scientific literature, to the genome annotation [44].

Initially, the cosmoss.org resource was set up to provide access to the P. patens virtual transcriptome assemblies and annotation [1] using BLAST services, keyword search and sequence retrieval. Subsequently, it was extended to provide services for splice site prediction [2], for mining gene families of transcription associated proteins (PlanTAPDB; [45]) and to predict dual protein targeting (ATP; [46]). Since the release of the initial P. patens genome assembly, the resource cosmoss.org additionally provides access to the draft genome sequence [10], the genetic map [47], and the filtered genome annotation V1.2 [3]. Moreover, cosmoss.org serves as a platform to coordinate the analysis and annotation of the P. patens genome sequence. As part of this, a wiki and several mailing lists have been set up to report and discuss the results within the community. Additionally, an integrative genome browser serves as a main entrance point for the exploration of the moss genome and annotation. The integrative cosmoss.org browser is based on the Gbrowse software [48] and provides base pair level resolution for large-scale annotation data covering predictions for all different kinds of genomic regions ranging from protein-coding genes, transposable elements and repeats to tRNA, rRNA, miRNAs, and other non-protein coding RNAs. Furthermore, the annotations are linked to the cosmoss.org internal annotation resources as well as to GenBank [49], Pfam [50], miRBase [51] and comparative genomics resources like Phytozome [39] and PLAZA [40].

While a significant part of cosmoss.org is based on our analyses, we are continuously integrating external published data, e.g. sRNAs [52], miRNAs [53] and EST (expressed sequence tag) or short read data from the sequence read archive (SRA [54]) and from collaborators around the world. In addition, the cosmoss.org gene annotation releases are shared and hosted at the NCBI and the comparative plant resources, Phytozome, PLAZA, and PlantGDB [55]. In July 2009, the Physcomitrella community annotation services were transferred from the JGI website to cosmoss.org and the resource now functions as the central annotation repository for the moss P. patens.

Here, we report the complete re-annotation of the P. patens genome assembly V1 (Table 1), demonstrate its utility for comparative analyses and introduce the extensions to the cosmoss.org resource to act as a permanent central genome annotation repository and model organism database involving:

• The improvement of gene structures with specific focus on the incorporation of transcript evidence to cover alternative splice variants and to derive UTRs.

• Complete renewal of the functional annotation.

• Prediction and annotation of non-protein-coding genes.

• Integration of user annotations: Structural and functional annotations.

• Integration of manual annotations: Development of community annotation services.

Table 1. P. patensgenome annotation releases

Results and discussion

Improved structural annotation of the P. patens genome

The cosmoss.org Physcomitrella patens V1.6 genome annotation reported here is the result of iterative rounds of evidence mapping, repeat masking, gene structure prediction, filtering and model selection and harbors annotation of protein-coding genes, transposable and repetitive elements and, for the first time, definition of non-protein-coding loci. The release comprises 32,275 protein-coding genes, 432 tRNA loci, 798 rDNA regions, 229 miRNA precursors (108 families) [51], 213 snRNA genes, and 6 SRP (signal recognition particle) loci. Considering the number of miRNA families, P. patens with 108 families has an intermediate position between the green alga C. reinhardtii (47 families) and the flowering plant A. thaliana (187 families) [51]. Consistent with previous findings [3,10] about half of the genome consists of full length LTR retrotransposons and related fragments including chromodomain-containing gypsy LTR retrotransposons (Tcn1) shared by fungi and non-flowering plants [56].

V1.6 protein-coding gene predictions are based on multiple sources of evidence. Of prime importance are ESTs (Additional file 1: Table A1) from 19 different experimental conditions, tissue types and developmental stages providing a reliable basis for gene structure prediction. The combined transcript evidence was used to train species-specific prediction models using SpliceMachine [57] and EuGène [58], followed by the generation of weighted consensus gene structures using EVidenceModeler [59] and PASA [60]. The prediction procedure was repeated iteratively. Each round involved several filtering steps for separation of non-protein coding, repeat and transposable element-associated genes. The interim versions V1.3-1.5 were not published but are the basis for V1.6 (see Additional file 2: Table A2). The whole protein-coding gene prediction and annotation process is summarized in Additional file 3: Figure A1. The annotation release V1.6 comprises 32,275 loci coding for 38,357 protein-coding transcripts (see annotation releases overview Table 1). 26,722 (~83%) loci are supported by transcript evidence (i.e. EST or full-length cDNA). The average V1.6 gene has a mean length of 2369 bp and a transcript length of 1389 bp. The number of only 1582 unchanged gene models from V1.1 to V1.6 is an excellent indicator of the extent of changes and improvements that led to the current release. While the changes in release V1.2 (Table 1; [3]) were restricted to the removal/filtering of non-protein-coding genes (7972 filtered models), the complete annotation process leading to V1.6 (Additional file 3: Figure A1) resulted in 22,307 (~80% of V1.2) updated models, 8387 (25% of V1.6) new loci and 4077 (15% of V1.2) models which have been removed due to non-protein coding/transposable element origin. Of the new loci, 1535 transcripts (1338 genes) are part of a gene family with at least one additional plant species, and 2196 transcripts (1456 genes) encode at least one Pfam domain. Also, published genes and gene models released by the scientific community were mapped, manually curated and integrated into V1.6.

Additional file 1: Table A1. P. patens transcript evidence.

Format: XLS Size: 32KB Download file

This file can be viewed with: Microsoft Excel ViewerOpen Data

Additional file 2: Table A2. Gene prediction performance on P. patens reference genes.

Format: XLS Size: 33KB Download file

This file can be viewed with: Microsoft Excel ViewerOpen Data

Additional file 3: Figure A1. Annotation process overview chart.

Format: PDF Size: 614KB Download file

This file can be viewed with: Adobe Acrobat ReaderOpen Data

Inspection of UTR annotations indicates that gene structure completeness is much improved in V1.6. The number of protein-coding loci with both 5’ and 3’ UTRs were increased from 4515 to 15,757 transcripts. The median transcript length increased from 987 (V1.2) to 1248 bases (V1.6) which can be explained by a higher percentage of annotated UTR regions in this release. V1.6 gene models only contain UTR annotations if they are fully supported by transcript evidence.

To further assess model quality and completeness we compared predicted protein sequences from V1.2 and V1.6 based on the coverage of their respective closest homolog in A. thaliana (Figure 1). The direct comparison demonstrates V1.6 models generally as more complete than those in V1.2: 34.5% of the V1.6 proteins with a corresponding V1.2 model cover a higher proportion of the length of their best A. thaliana BLASTP hit.

thumbnailFigure 1. A. thaliana best hit (BLASTP) coverage changes from P. patens V1.2 to V1.6 34.5 % of the protein-coding gene models (V1.6) covers better their closest A. thaliana homolog.

Annotation and characterization of alternative splicing

We used the PASA pipeline [60] to study the extent and characteristics of alternative splicing (AS) in the moss. We used GenomeThreader [61] transcript alignments in exchange of PASA’s standard GMAP alignments [62] as GenomeThreader supports alignment of transcript evidence to multiple (nearly identical) loci. This is important for covering nearly identical tandemly-arrayed genes observed in P. patens[3,10] and of other segmental duplications.

In total, the PASA pipeline reported 39,106 genes (subclusters) out of which 6556 (16.8%) show AS (Table 2). Consistent with previous analyses [10], we find alternative acceptors as the most prevalent event of AS (4443 events) in PASA assemblies. Among gene loci with alternate transcripts, intron retention is the most frequent form (~40% of all genes). Analysis based on large-scale EST libraries largely derived from total RNA preparations, undoubtedly include partially spliced intermediates which might lead to an overestimation of the fraction of intron retention variants, however our data reflect the findings observed for other plants where intron retention is the preferred form of AS, followed by an alternative acceptor [63-65]. Based on the Sanger EST data the percentage of genes with evidence for AS in P. patens (21%) is in the range of flowering plants like A. thaliana (20%) and Oryza sativa (30%). Given the fact that recently the use of next generation deep sequencing data resulted in a substantial upward revision of the number of alternatively spliced genes in A. thaliana to at least 42% [66] we expect the true extent of AS to be similar in the moss. This is also supported by the extent of NAGNAG alternative splicing based on NGS data which is in the same range as in flowering plants [65].

Table 2. Summary statistics of alternative splicing in P. patens genes

In V1.6 we have incorporated information on AS for the first time into the released genome annotation. The integration of AS into the annotation process leads to 3500 (10.8%) annotated loci with an average of 2.52 transcripts per locus and a maximum of 11 transcripts in V1.6. 1775 loci (51% of AS transcripts; 5.5% of all genes) have an altered coding sequence (CDS) due to alternative splicing resulting in 2380 distinct proteins. In contrast to A. thaliana, where the analysis of large-scale full-length cDNAs suggested most splicing events to occur outside of coding regions in the 5’-UTR [67], alternative splicing in moss seems to affect UTRs and CDS regions to a similar degree. 2948 alternative transcripts are due to alternative splicing in the UTRs of 1991 loci (56% of AS transcripts; 6.2% of all genes). Extensive alternative splicing of moss 5’-UTRs was observed previously in an individual study of the MDHAR genes [68].

Insights into the evolution of gene structures along the green lineage

The advanced annotation enabled us to address comparative questions on genome and gene structure evolution in plants. A comprehensive overview of genome and protein-coding gene statistics of several Viridiplantae species is compiled in Table 3 and Additional file 4: Table A3. Our complete re-evaluation of the respective genome annotations supports previous observations, but also provides novel insights. The average volvocine green algal gene is about twice the length of the average land plant gene [69]. This trend is not as pronounced in the transcript and CDS lengths, but still the average algal transcript is longer than those of the land plants (Table 3). While the median exon length is similar in all seven Viridiplantae species, Chlamydomonas introns are longer and more numerous than those of land plants ([69]; Table 3).

Table 3. Protein-coding gene statistics of selected Viridiplantae

Additional file 4: Table A3. Protein-coding gene statistics of selected Viridiplantae.

Format: XLS Size: 55KB Download file

This file can be viewed with: Microsoft Excel ViewerOpen Data

Corresponding with its intermediary phylogenetic position, the median intron length of the average P. patens gene is intermediate (Table 3) between that of algae and flowering plants, implicating a trend towards shorter introns during the evolution of land plants. Whereas the median algal 5’-UTR length are in the range of Tracheophyta, the median moss 5’-UTR length is about 100 bp longer. Moss 5’-UTR regions more frequently comprise multiple exons than other plant genomes, with 40% of genes so annotated. Additionally, although there are substantially more predicted transcripts in O. sativa and more transcripts with 5’-UTR introns (8227), the total number of transcripts with multiple 5’-UTR exons in P. patens is nearly as high (7120) (p-value < 2E-238); Figure 2). This will likely increase as further high-resolution transcriptome data are obtained, since only 18,180 moss transcripts currently have annotated 5’UTRs compared with 31,793 transcripts in the rice annotation.

thumbnailFigure 2. Comparison of 5’-UTR intron numbers in Viridiplantae5’-UTR intron number frequencies of selected Viridiplantae genomes. The y-axis labels give the number of transcripts w/o 5’UTR introns in percentage of all transcripts with 5’UTR.

In A. thaliana introns in the 5’-UTRs tend to be longer than introns in the coding sequence and 3’-UTRs [70]. A plot of P. patens intron lengths (Figure 3) reveals a similar distribution. P. patens 5’-UTR introns are significantly longer than introns in the CDS (p-value = 0) and 3’-UTR (p-value <1.03E-193). In particular, the percentage of introns longer than 500 bp is much higher in the 5’-UTRs than in the CDS or in the 3’-UTRs. The mean 5’-UTR intron length for P. patens is 520 bp, whereas it is 264 bp for the CDS and 268 bp for 3’-UTR introns.

thumbnailFigure 3. P. patens V1.6 intron lengths distribution 5’-UTR, CDS, and 3’UTR intron lengths in comparison. The percentage of introns longer than 500 bp is much higher in 5’-UTRs than in CDS and 3’-UTR introns.

Furthermore, in A. thaliana the 5’-UTR introns are preferentially located close to the initiating ATG codon [70]. A comparative distribution plot of the proximate 5’-UTR intron distance to the start of transcription and translation reveals similar distributions in P. patens and A. thaliana (Figure 4). Certainly the closeness of 5’-UTR introns to the initiating ATG is more pronounced in A. thaliana (A. thaliana ~75% < 65 bp; P. patens ~50% < 65 bp; p-value < 9.7E-165). Chung et al. (2006) provided initial evidence for a positive effect of the presence of 5’-UTR introns on gene expression. These findings and previous reports suggest that properties of 5’-UTR introns influence gene expression level, regulation, translation and nonsense-mediated mRNA decay [71,72]. As a consequence, the high percentage of genes with 5’-UTR introns, their unique length, and the fact that ~50% of P. patens loci with annotated alternative transcripts undergo alternative splicing in UTRs, suggests that the moss makes frequent use of this type of gene regulation.

thumbnailFigure 4. Distance to translation and transcription start sites of 5-UTR intron positions Distribution of 5’-UTR positions for P. patens and A. thaliana transcripts in comparison. The closeness of 5’-UTR to the initiating ATG is more pronounced in A. thaliana. While ~75% of introns are closer than 65 bp in A. thaliana only ~50% are in P. patens.

Improved functional annotation of moss proteins

The predicted V1.6 protein sequences were functionally annotated by homology transfer using BLAST2GO [73] and InterproScan [74]. Existing Gene Ontology Annotations (GOA) from previous releases were mapped, curated and integrated. GOA were extended using a novel in-house subcellular localization annotation pipeline called pred2GOA which allows the integration of GOAs from multiple sources by weighted combination of results from experimental evidence, subcellular target predictions and homology-based methods (see Methods for details) to improve GOA accuracy. The resulting non-redundant, transcript-wise P. patens V1.6 GOA contains 66,234 terms (Table 4). While we could increase the number of terms in the GO namespaces “biological process” (+4900 terms) and “cellular component” (+19,878 terms), the number is reduced for “molecular function” (-15,479 terms) in comparison to V1.2. However, comparison of both GOAs reveals a prominent set of protein kinase family proteins (248 sequences) with a very high number of assigned terms (in total 9176) in V1.2. Manual inspection identified these terms as false assignments. Excluding these false-positive terms leads to a general increase of terms and annotated gene products in V1.6 (Table 4). In total, 18,786 loci, i.e. 58% of all predicted protein-coding genes, have been assigned at least one GO term. This is a substantial improvement over V1.2 (41%). While we could increase the number of genes with a term for “biological process” (BP) by 2285 and “molecular function” (MF) by 2702, the number of genes with a “cellular component” (CC) term has increased by 10,065. A significant part of the assigned “cellular component” terms are based on supervised subcellular target predictions (pred2GOA, see Methods). To help scientists in evaluating the quality of annotations, we rely on the assignment of evidence codes [75], which offer a direct interpretation of the support for each GOA. After the automatic processing of terms, most of the annotations (~77%) is still supported only by IEA (Inferred by Electronic Annotation) followed by ISS (Inferred from Sequence or Structural Similarity) with ~23%. As noted earlier, only a few genes have annotation supported by experimental evidence. This now becomes obvious in the GOA - only 18 genes have evidence codes indicative of experimental support (0.04% of all assigned GO terms). With the development of our cosmoss.org community annotation service discussed below we hope to improve this situation and invite scientists to annotate their genes of interest. Building on and extending the GAF format, our databases and web service also offers additional information for every GOA describing the algorithmic or experimental source, references and crosslinks which can be used to further assess the quality of annotations.

Table 4. Comparison of the Gene Ontology (GO) annotation of P. patens V1.2 and V1.6

The most common and useful application of GOA is the comprehensive functional analysis of experimentally defined gene sets by ontology term enrichment analysis [76]. The significance of results obtained by enrichment analysis strongly depends on the quality and depth of the underlying GOA. In order to assess the utility of the V1.6 GOA for exploratory enrichment analysis, we compared the Arabidopsis and Physcomitrella GOAs to test whether we could reproduce lineage-specific expansions previously discovered using phylogenetic and phylogenomics approaches ([10]; Table 5). In comparison to A. thaliana, gene families like flagellum associated dyneins [10,69], the light harvesting (LHC) superfamily including chlorophyll a/b binding proteins and early light-inducible proteins [10,13], histidine kinases and response regulators [10], phenylalanine ammonia lyases (PALs; [77]) and aldehyde dehydrogenases (ALDHs; [78]), are expanded in P. patens. Comparison of the underlying gene lists of the terms identified by GO enrichment analysis with those of the phylogenetic studies reveals that these specific expansions are also reflected by comparing functional annotations using statistical methods.

Table 5. Selected GO categories: P. patens in comparison to A. thaliana

In summary, we improved the GO annotation of P. patens gene products in general and demonstrate that the annotation can serve as a solid basis for comparative and exploratory analyses. Compared to A. thaliana (TAIR10), where only 7% of all genes have no GO annotation, there are still more than 40% of all loci without any GO annotation for P. patens, and those with annotation are mainly inferred by electronic annotation (IEA). While a significant fraction of this 40% probably comprises orphan genes, special focus should be placed on the improvement of the functional annotation of the P. patens protein-coding genes to unravel its full potential as a reference model organism. To facilitate such efforts we developed the cosmoss.org community annotation interface to browse and alter the P. patens annotation (see section genonaut) discussed below.

Gene families

Gene families (clusters) were retrieved by clustering a multitude of sequences of Archaeplastida proteins (Additional file 5: Table A4) using OrthoMCL [79]. Parameters were optimized as described in Methods to target Archaeplastida gene families sensu strictu, i.e. families of genes that evolved by speciation and duplications after the divergence of the red/green lineage from a single gene in the last common ancestor. Multi-gene or superfamilies thus are split across multiple clusters. The number of P. patens loci in clusters is 32,733, while after the subtraction of the annotated number of 32,275 protein-coding loci, 458 fall into multiple clusters. Some of these cases are due to fragmentary or false structure predictions caused by fusion of two or more distinct gene loci into a single locus and will be resolved in future.

Additional file 5: Table A4. Archeaplastida genomes sources.

Format: XLS Size: 38KB Download file

This file can be viewed with: Microsoft Excel ViewerOpen Data

Strikingly, 48% of all P. patens loci are in Physcomitrella-only clusters (Figure 5), which independently supports an analysis based on protein domains where 52% of all P. patens genes have no Pfam domain [80]. Further filtering of these loci including BLASTP against GenPept and EST support indicates that ~22% (7169) of all loci in Physcomitrella-only clusters have no detectable homolog, while at least ~13% (4157) have no homolog but transcript evidence. These genes might represent true orphan genes, representing species- or lineage-specific adaptive innovations or non-protein-coding genes. The remainder falls into 8208 clusters with at least one other Viridiplantae protein (Figure 6). Thus, a prominent part of the P. patens gene complement could not be clustered with other plant proteins. In part, this could be explained by remaining fragmentary gene predictions and non-protein-coding genes. This is supported by the observation that the average transcript length of these P. patens-only singleton genes is less than 500 bp and the average number of exons is only two (total average is four exons, Table 3). However, at least 40% (4503) of these P. patens-only singleton genes have transcript evidence. Approximately two-thirds of all genes are in small clusters with less than three P. patens members. As already observed in the draft genome [10], several gene families are expanded by comparison to flowering plants. In total 7013 OrthoMCL clusters are shared by P. patens and A. thaliana containing 14,830 and 11,926 genes, respectively. Reported expansions in gene families [10], like LHCs (esp. chlorophyll a/b binding proteins), dyneins, histidine kinases and response regulators are also reflected by OrthoMCL cluster sizes (Figure 6). The majority of expanded families are in OrthoMCL clusters with more than five members in moss.

thumbnailFigure 5. How many moss-specific genes are there? BLAST hits of P. patens-only clusters based on our OrthoMCL clustering with selected Viridiplantae genomes against GenPept (rel. 190). P. patens proteins were excluded from GenPept for this analysis.

thumbnailFigure 6. Gene family sizes in P. patens proteins from several Viridiplantae were clustered using OrthoMCL. Depicted are all protein clusters with regard to P. patens and sorted by cluster size. The clusters were subdivided into P. patens only clusters and clusters with at least one other member. Protein families found to be expanded in P. patens in comparison to A. thaliana are listed.

Analysis of intron-loss and gain in the green lineage

Similar to 5’-UTRs, where we observe a striking number of multi-exon regions in moss, the overall number of single-exon transcripts is remarkable. While Chlamydomonas and Volvox contain fewer than 10% single-exon transcripts (Table 3 and Additional file 4: Table A3), this fraction is on average ~19% in other land plants, while moss possesses 23.4%. This may be due either to fragmentary gene predictions and residual non-protein coding genes, or may reflect secondaryintron losses.It has been observed that introns and their positions are highly conserved during land plant evolution [81]. In A. thaliana and O. sativa intron losses outnumber intron gains [82]. Furthermore, there is evidence for secondary intron loss; e.g. a moss sedoheptulose-1,7-bisphosphatase (SBP) gene which lost six out of seven introns [81]. Such intron losses might account for the relatively high percentage of single-exon genes in P. patens. One suggested mechanism for intron loss involves the reverse transcription of an mRNA followed by the (partial) replacement of the genomic DNA copy by an intron-less cDNA via homologous recombination, called retrocopying [83]. An increased rate of intron loss might thus be facilitated by the extraordinarily high rate of DNA repair by homologous recombination in P. patens. This hypothesis is based on models proposing a prominent role of gene conversion and DNA repair in intron loss [84,85].

Based on shared gene families the extent of secondary intron loss can estimated. There are 805 clusters with single-exon moss genes that have multi-exon homologs in other Viridiplantae. About half of them (404) contain also moss paralogs with exons in the same cluster. As some of these might be due to fragmentary gene models, we excluded models from the analysis that are significantly shorter than the median transcript lengths of non-Physcomitrella cluster members. A stringent way to define intron loss is to consider the median number of exons found in genes of other species in the same cluster. Based on these criteria, we find evidence for secondary loss of introns in about 14% (4405) of all P. patens genes (Table 6).

Table 6. Number of introns in Viridiplantae

About 3% (941) of the P. patens genes seem to have lost their introns entirely, i.e. are of putative retrocopy origin, which is in the range of the other land plants (Table 6). The extension of the analysis to Arabidopsis, Rice and Chlamydomonas supports the findings from the comparison of absolute numbers of single-exon genes between algae and land plants described in the previous sections. The alga has significantly less (0.7%) single-exon genes than the three land plants (3-4%) under study. One likely scenario is an increased activity of transposons resulting in a secondary, maybe more recent, intron gain in algae [69,86]. This view is supported by the observation that intron positions are often not conserved between the two algae and the land plants (data not shown). Considering the comparable rate of intron loss in plant gene families, the fact that the total number of multi-exon transcripts in P. patens is similar to vascular plants, that only 60% of the 8979 single exon models are supported by expression evidence, and that more than half of these transcripts are shorter than 500 bp (which is less than half of the mean of P. patens transcripts, 1389 nt), leads us to conclude that a significant number of the predicted single exon genes represent fragmentary predictions, non-protein coding genes or pseudogenes.

Gene family size evolution in Viridiplantae

We compared gene family size distributions along the green lineage by looking at all clusters containing P. patens and at least one other Viridiplantae species (8208 clusters/families; Figure 6). This analysis revealed that protein families with two to six members are more frequent in P. patens than in the other Viridiplantae (Figure 7), supporting the previously hypothesized balanced retention of paralogs that act as pseudoalleles in a haploid-dominant land plant [3]. This is most pronounced in clusters with two members, which amount to 27% in P. patens while being substantially less frequent in the other green organisms. In total, 832 clusters point to an expansion in P. patens in comparison to the other Viridiplantae. 184 of these P. patens gene families are more than twice the size of their largest green lineage counterpart.

thumbnailFigure 7. Comparison of gene family sizes in conserved clusters distribution of genes per cluster and species common to P. patens and at least one other Viridiplantae (8208 cluster).

Among these expansions house-keeping and metabolic gene functions are most prominent, independently supporting analyses of paralogs retained after the proposed whole genome duplication event ~45 million years ago [87] and of the unique presence of identical tandemly-arrayed genes [10]. Prominent examples of this expanded category of genes are abundantly expressed components of multimeric protein complexes like the ribosome and proteasome. In total, 86 of the clusters represent the different abundantly expressed protein components of the ribosome. On average these harbor ~1.5 times more genes in P. patens than in A. thaliana. For the 36 clusters representing structural components of the proteasome we observed on average 1.6 times the gene complement of Arabidopsis in the moss. These expansions are also detectable in the functional comparison of the two species using GO enrichment analysis (Additional file 6: Table A5).

Additional file 6: Table A5. topGO Gene ontology (GO) enrichment analysis between P. patens and A. thaliana.

Format: XLS Size: 1.6MB Download file

This file can be viewed with: Microsoft Excel ViewerOpen Data

In addition, families encoding for smaller complexes and monomers are expanded, including the Light Harvesting Complex II (LHCII) major antenna, Ribulose-1, 5-bisphosphate carboxylase oxygenase (RuBisCO) small subunit, TOC12 (translocase of the outer chloroplast membrane 12) and components of the splicing and translation machineries.

The expansion of house-keeping and metabolic gene functions is also mirrored by the findings of the GO enrichment analysis comparing P. patens and A. thaliana (Table 5 and Additional file 6: Table A5), which revealed an increased complement of genes involved in translation, oxidation-reduction, electron transport, microtubule based movement, glycolysis, and ATP synthesis. Additionally, specific expansions occurred which might represent lineage-/species-specific adaptations (e.g. expansins, MIKC*-type MADS box transcription factors, late embryogenesis abundant (LEA) proteins, early response to dehydration (ERD) proteins, cationic peroxidases) and in enriched GO “biological process” annotations (e.g. chitin catabolic processes, the phosphoenolpyruvate-dependent sugar phosphotransferase system, cell wall macromolecule catabolic processes, phosphatidylinositol-mediated signaling, chromatin assembly or disassembly, cell redox homeostasis, ciliary or flagellar motility and DNA repair). Some of these enriched categories like “ciliary or flagellar motility”, which can be explained by the absence of flagellated sperm in flowering plants, confirm the findings of previous analyses [10,69,88]. The majority of the above listed categories and families (nine out of 13) represents true novel insights which will help to unravel the much-cited qualities of the mosses in coping with abiotic and biotic stressors like their unique ability to repair DNA damage [89-91].

Extension of the cosmoss.org resource to provide a permanent, central model organism database and annotation repository for P. patens

Protein-coding gene models were clustered to loci and the resulting locus definitions were used to derive information-rich locus identifiers (Cosmoss Gene ID; CGI; see Additional file 7: S2 and Additional file 8: Figure A2, and [92] for details).

Additional file 7. Supplementary_text - Codon usage & Protein-coding locus definition.

Format: DOCX Size: 17KB Download fileOpen Data

Additional file 8: Figure A2. Gene ID naming – Cosmoss.org gene identifiers.

Format: TIFF Size: 85KB Download fileOpen Data

With the development of genonaut [93], we have extended cosmoss.org by the capability to annotate Physcomitrella patens genes with regard to gene name, product name, description, and Gene Ontology (GO) terms. The interface which allows searching, browsing and editing of annotation is modular and can be extended to support additional (ontology) annotations as gene features. Traceability of annotations in terms of author and experimental evidence is crucial for quality assessment of information retrieval. Thus, the genonaut interface accepts the alteration of an existing gene description only if the source is specified. Integration of multiple sources into a unique abstraction layer is achieved by assigning unique and permanent Cosmoss Reference IDs (CRID) to author statements from all sources. The highest quality author statements are experimental evidence provided as references to peer-reviewed publications. The easiest way to achieve this is to provide a valid PubMed ID and the system automatically retrieves all relevant information from NCBI PubMed. If the source is a publication that is not tracked in PubMed, a custom reference can be created. If no publication is available, as the “weakest” possible evidence, a note in form of a text comment or web link describing the evidence is required. Besides the references the genonaut interface allows to link to other resources via database cross-references (Dbxref).

The gene products can be further annotated using GO terms. To allow convenient manual GOA, the genonaut interface assists the annotator by allowing the user to browse the appropriate GO namespace by keywords to assign the correct terms. In addition to the mandatory reference for a genonaut annotation, we have integrated the assignment of the GO evidence codes [75]. In this way the quality of each assigned GO term is directly discernible. As traceability is crucial for the maintenance of annotation quality, the genonaut system traces every annotation change using a history system. Thus it is possible to trace changes and possibly revert to a previous state if needed, but more importantly to comprehend the annotation history of every gene and annotation version.

Whereas the annotation browser capabilities are publicly available, the editor functions are restricted to registered cosmoss.org users (December 2012: 228 annotator accounts). The cosmoss.org curator team acts as a superior authority which supervises and validates the user provided annotations by direct personal communication.

Moreover, the genonaut interface provides a starting point to retrieve detailed annotation about P. patens protein-coding genes. Besides the possibility to search and edit the annotations and annotation history, the genonaut interface is linked to the sequence retrieval, genome browser and sequence viewer providing transcript and protein domain annotations.

To support the manual curation of gene structures we have integrated and adapted the Apollo structural gene annotation editor [94] for the cosmoss.org genome browser. Generated user_models (December 2012: 830 manually curated transcripts) are assigned CGIs (Additional files 5 and 5) from the user model namespace extended using the authors username to allow multiple versions per locus (e.g. Pp1s275_35U2 ➔ Pp1s275_35U2__zimmer.1).

Conclusion

Here we describe the complete re-annotation of the P. patens V1 genome assembly comprising structural and functional annotation of protein-coding genes and, for the first time, description of non-protein loci including tRNA, rRNA, miRNA, snoRNA and snRNA loci. Compared to V1.2 the improved structural annotation V1.6 resulted in 8387 additional protein-coding loci, 11,242 more complete genes and only 1582 unaltered gene structures. 70% of the 32,275 protein-coding genes are supported by EST evidence. Nearly half (~49%) of the protein-coding loci in V1.6 are now be considered complete, containing both UTRs. Furthermore, the information-rich cosmoss.org locus IDs also carry information on the chromosomal/scaffold localization and about alternative splicing of transcripts.

We significantly increased the number of genes with functional annotations (58% as compared to 41% in V1.2) in form of GO term annotations (GOA). Our quality assessment of the V1.6 GOA demonstrates sufficient annotation depth to recover results from previous high-quality phylogeny-based approaches using ontology term enrichment analysis. Nevertheless, there are still 41% of all loci without any GO annotation and only 0.04% of all assigned GO terms are supported by direct experimental evidence. Although this is a common phenomenon for most available plant genomes, special focus needs to be placed on the improvement of the functional annotation of the P. patens protein-coding genes in order for it to serve as a reference model organism and plant “flagship”. With the development of the cosmoss.org community annotation services described here allowing users to browse, view and alter functional and structural annotations of moss genes, transfer and exchange of knowledge is greatly facilitated. Including the described extensions of the resource, cosmoss.org now is well-equipped to serve as a permanent, central model organism database and annotation repository for P. patens.

We demonstrate the utility of the provided annotation and resources for the comparative study of plant evolution including the analysis of codon usage (see Additional file 7: S1 and Additional file 9: Figure A3), alternative splicing, gene structure evolution as well as the detection of lineage- and species-specific expansions of gene families and biological processes.

Additional file 9: Figure A3. Codon usage of P. patens and A. thaliana transcripts.

Format: PDF Size: 51KB Download file

This file can be viewed with: Adobe Acrobat ReaderOpen Data

Results from our comparative analyses were mostly consistent with previous observations, but also provided several novel insights. In particular, we found further evidence for intron loss during land plant evolution and secondary intron gain in the alga Chlamydomonas. Investigation of alternative splicing and gene structures revealed a unique complexity of 5’-UTRs in the moss, pointing to the importance of UTRs for the regulation of gene expression in this early diverging land plant.

Our comparative analysis of functional annotations and protein clusters revealed expansions of moss house-keeping and metabolic gene functions as well as hitherto unknown lineage-specific expansions. In total, 832 gene clusters are expanded in P. patens and at least ~13% of all gene loci are orphan genes as they have no homolog in other as yet published genomes. Subsequent functional analysis of this data set will further extend our understanding of the unique capabilities to cope with abiotic and biotic stressors and to efficiently repair DNA damage.

Methods

P. patens reference genes

In total, 137 manually annotated and validated Physcomitrella patens gene structures are in the cosmoss.org genome browser track “Ppref genes”. Some of the genes are directly derived from published genes (GenBank) or provided by the scientific community, but the majority was extended or corrected manually using ESTs and FLcDNAs. In addition, in-house validated sequenced P. patens gene structures were added.

EuGène P. patens gene prediction process

EuGène [58] allows combining various type of evidence including e.g. ESTs, mate pair information, homologous sequences, existing gene predictions and splice site predictions. Our EuGène predictions for P. patens are based on a splice site prediction using SpliceMachine [57] trained on filtered P. patens EST alignments, coding, intronic and intergenic regions, homology evidence (A. thaliana homologs) and are filtered using transposon related sequences and gaps. The optimal parameters were determined on an independent P. patens reference genes set not used for training of EuGène. Two P. patens whole genome EuGène predictions went into generation of the release V1.6. The first contains 37,872 predicted loci and was restricted to generate UTR regions only if transcript evidence is available and the second contains UTR regions predicted ab initio (46,071 loci).

The first EuGène model predicts ~94% of all CDS exons in the reference genes correctly, whereas the second does so for ~95% (results are summarized in Table 3). Both predict 76% of all reference CDS without any error. That implies that the EuGène predictions perform well in predicting the exons but split several loci into two or more distinct genes. While manually inspecting these problematic loci we have noticed that gene models created by the JGI [10], which were not been selected for release V1.1, were often better than the selected model model, or could be used to overcome, or respectively complement, the EuGène predictions. As the method of choice to combine all available P. patens evidence and to further improve the protein-coding gene structures we have used EVidenceModeler [59].

EvidenceModeler (EVM) - weighted consensus gene model predictions

EVM (Haas et al. 2008) combines evidence from different sources into a consensus gene structure prediction. With the possibility to weight and unite the different evidence and optimize their combination, EVM utilizes the different sources by equating the drawbacks of individual sources but also boosting their strong points. For P. patens we have used EVM to find the optimal combination of PASA [60] transcript assemblies, EST alignments and five different whole genome protein-coding gene predictions. The process is also described in Additional file 3: Figure A1. The resulting models were subsequently subjected to PASA to model the UTR regions. As a consequence, all UTRs in release V1.6 are supported by transcript evidence. The utilization of EVM has enabled us to increase the prediction performance on the reference gene set (86.1% of all CDS and 97.3% of all CDS exons are correct; see Additional file 2: Table A2).

Additional gene structure predictions using EuGène

EuGène [58] version 3.4 was adapted and trained for P. patens on the basis of the Ppref genes set (mentioned above) and including a species-specific splice-site prediction and IMM (Interpolated Markov Model) models trained on intergenic, CDS and UTR regions. Gaps in the genomic sequence and repetitive regions, in particularly LTR retrotransposons [10], were masked for the training and predictions. As additional evidence, homologous protein sequences (Swissprot rel. 13.4 and Arabidopsis thaliana TAIR7 [34] homologs) and in particular EST alignments 360,974 from GenomeThreader [61], 118,243 from sim4 [95] and 97,373 from exonerate [96] were used. If available, we also provided EST mate pair information into EuGène (63,945 EST mate pairs). Two whole genome EuGène predictions were used for the consensus model approach leading to V1.6. The training input of these two models was the same except that one model was additionally trained with 5’- and 3’-UTR regions (EuGène MarkovIMM plugin).

Splice site prediction

P. patens EST alignments were loaded into a Bio::DB::Seqfeature database. The exclusion of alternative splicing and bad quality EST alignments led to a distinct, species-specific splice-site training set for SpliceMachine. Splice sites were only taken into account where GenomeThreader [61], sim4 [95] and exonerate [96] did exactly the same EST alignment. Two models, one donor (GT) and one acceptor (AG) model have been generated and have been used for EuGène training.

PASA – transcript assemblies

PASA [61] assemblies were performed as directed on the software homepage [97], with the following modifications: Per default PASA uses GMAP [62] for transcript alignments, however, our evaluation process (data not shown) reveals improved EST alignments using GenomeThreader. PASA offers the possibility to include alignments in GFF3 format via the --IMPORT_CUSTOM_ALIGNMENTS_GFF3 switch. The P. patens GenomeThreader EST alignments were converted into PASA compliant GFF3 format. PASA per default supports only one alignment per EST/transcript. Therefore with regard to duplicated genes (especially tandemly arrayed genes and (near-)identical genes) the corresponding EST alignments were renamed. E.g. the sequence ppsp14d22fl matches to 2 loci in the genome, so the corresponding sequence is duplicated and renamed: ppsp14d22fl and ppsp14d22fl_2. Transcript alignments with less than 90% EST length coverage or less than 95% alignment identity were discarded. The maximum allowed intron length was set to 20,000 nt, based on the longest observed intron supported by Sanger ESTs. Three cycles of PASA annotation loading, annotation comparison, and annotation updates were used to maximize the incorporation of transcript alignments into the transcript assemblies.

Consensus gene predictions – unfiltered V1.6 models

The P. patens JGI AllModels V1.1 (http://ftp.jgi-psf.org/pub/JGI_data/Physcomitrella_patens/v1.1/transcripts.Phypa1_1.AllModels.fasta.gz webcite) and the cosmoss.org EuGène models, described in the previous section, and transcript alignments (PASA, exonerate and sim4) were combined and weighted with EVidenceModeler (EVM). EVM combines all evidence per locus into one consensus gene structure model. EVM was trained on the Ppref genes deducting reference genes used for the training of EuGène. Subsequently the UTR regions were modeled by P. patens PASA transcript alignment assemblies as described on the PASA homepage [97]. In this context models with alternative splicing evidence were generated and incorporated into V1.6.

Functional gene annotation

The unfiltered V1.6 gene models and the corresponding predicted protein sequences, respectively, were subjected to BLAST2GO [73]. The initial BLASTP (e-value cut-off: 1E-4) search was performed against a P. patens V1.1 subtracted GenPept release 172.0 to allow an unbiased functional annotation, independent from the previous entirely automatic V1.1 GOA. For the BLAST2GO annotation step the minimum coverage between a hit and its HSP was set to 40%. The validation step as well as the integration of InterProScan V4.5 (InterPro release v22.0) based GO annotation was used to generate the GOA for the P. patens V1.6 proteins.

GO term annotation was extended by experimental evidence, various subcellular target predictions and homology-based methods using the pred2GOA method described below. Existing functional annotations like gene names, description lines, GO terms and KEGG EC numbers and KO terms were collected from JGI and Kyoto Encyclopedia of Genes and Genomes [98] and combined into a non-redundant database. Existing experimental evidence for subcellular localization was also manually integrated as GOA if available. Additional functional annotations were created using the homology-based methods BLAST2GO, IPRScan [74] and KAAS [99]. These two steps resulted in GO terms with IEA evidence codes. Subcellular localization of the protein sequences was predicted considering the individual gene’s full-length status using several tools: TAPScan [100], MultiLoc [101], WolfPSORT [102], TargetP [103], ChloroP [104], SignalP [105], Prosite KDEL [106], HMMTOP [107] and MEMSAT3 [108]. To combine these predictions with the existing GOA, we developed the algorithm pred2GOA which allows weighted integration of GOA from multiple sources. Resulting predictions were translated into cellular component GO terms and compared to the existing GOA at the GO slim (plants) level. The assignment of GO terms is based on a weighted majority rule consensus. If the underlying gene did not have an annotated UTR and a start codon, at least one of the predictors’ predictions had to be based on more than the N-terminal region of the protein. Resulting GO term assignments were reviewed for consistency and defined as “Inferred from Sequence or Structural Similarity” (ISS). The resulting GO annotation was mapped to GO slim terms using the Blast2GO internal mapping function using the “goslim_plant.obo” ontology subset.

GO enrichment analysis

The enrichment analyses were performed using the Bioconductor package topGO [109]. We used and compared results from both the classical Fisher’s Exact Test and the topGO algorithm “weight01” for test statistics with a p-value cut-off 0.05. The Arabidopsis GOA was downloaded from the TAIR ftp server (TAIR10; http://ftp.arabidopsis.org/Ontologies/Gene_Ontology/ATH_GO_GOSLIM.txt webcite).

Protein clusters and representative gene model selection

Although certainly useful in single gene analysis and the study of AS, the use of all splice-variants in large-scale comparative analysis introduces undesirable complications. Therefore gene catalogues are usually reduced to one representative isoform per locus prior to large-scale analyses. By convention, the representative model is the variant with the lowest splice variant index. A representative gene model per locus was selected by the following procedure: At the first step, all P. patens V1.6 proteins together with the proteins of various sequenced Archeaplastida (Additional file 5: Table A4) were clustered using OrthoMCL V.2 [79]. The inflation parameter was optimized using a set of reference gene families (e.g. [10,15]). The inflation value was set to 1.3. The resultant clusters were aligned with MAFFT ginsi v6.712b [110] and subsequently the uncorrected pairwise distances were calculated with distmat (EMBOSS 6.1.0; [111]). The distmat matrices were used to select the Physcomitrella model as representative for a locus that has the least of all substitutions (minimum distmat distance) to any non-Physcomitrella cluster member. Clusters with only Physcomitrella sequences were subjected to a BLASTP search against GenPept (release 172.0). In this case the representative per locus was set to the Physcomitrella sequence which covers its closest GenPept hit, in terms of alignment length, best. If still no clear representative could be determined, the model with the longest transcript length was chosen as the representative. Due to the fact that the clustering is based on sequence information alone, not all splice variants per locus were grouped into the same cluster in some cases. Thus, the number of loci in clusters slightly exceeds the number of physical loci.

RNA genes

Pre-tRNA genes were predicted by combining results from tRNAscan-SE 1.21 [112] and ARAGORN 1.1 [113]. rRNA loci were predicted using RNAmmer-1.2 [114] and extended by mapping the available SILVA Physcomitrella rRNAs [115] to the genome with BLASTN.

Non-protein coding loci and gene families where determined using Infernal [116] with the RFAM (version 8.1) [117] covariance models and by integration of the miRBase 18 [51]Physcomitrella miRNA classifications and annotations.

P. patens V1.6 protein-coding genes filtering

On the basis of a Bio::DB::SeqFeature database [118] the V1.6 gene models were filtered against the annotated LTR-retrotransposons [10], Repbase (RELEASE 20080801; [119]) using RepeatMasker v3.26 [120] and the non-coding RNA described in the previous section.

Codon usage

The effective number of codons was calculated with CodonW 1.4.4 [121].

Statistics

Fisher’s exact tests and Wilcoxon tests were performed with R 3.0.0. The p-values were corrected for multiple testing using fdr.

Availability of supporting data

The complete annotation and sequence information comprising: Gene structure releases in GFF3 format: gene, transcript, CDS, protein, UTR sequences in FASTA format; Mappings/lookup tables, GO annotations in GAF format are accessible via the cosmoss.org download section [122].

Competing interests

The authors declare that they have no competing interests.

Authors’ contributions

ADZ, SR carried out the protein-coding gene structure predictions and annotation. YP supervised the EuGène training and prediction. DL performed non-protein coding gene annotation. ADZ and DL developed the functional annotation pipeline, performed the functional annotation, and drafted the manuscript. DL and ADZ conceived, supervised and designed the curation process, annotation interfaces and database backends. KB implemented the genonaut web interface. TN and MH generated and contributed FLcDNA sequences, discussed results and helped in writing the manuscript. RR and SAR supervised the study and contributed to writing the manuscript. All authors read and approved the final manuscript.

Acknowledgements

We thank Tomomichi Fujita, Keiko Sakakibara, Atsushi Toyoda, Asao Fujiyama, and Yuji Kohara for the generation of FLcDNA libraries and sequencing. We are indebted to our collaborators at the US DoE JGI and the Moss Genome Consortium. Funding by the German Research Foundation DFG (RE 837/10-2) and by MEXT, JSPS and JST is gratefully acknowledged. The article processing charge was funded by the German Research Foundation (DFG) and the Albert Ludwigs University Freiburg in the funding programme Open Access Publishing.

References

  1. Lang D, Eisinger J, Reski R, Rensing SA: Representation and high-quality annotation of the Physcomitrella patens transcriptome demonstrates a high proportion of proteins involved in metabolism in mosses.

    Plant Biol (Stuttg) 2005, 7:238-250. Publisher Full Text OpenURL

  2. Rensing SA, Fritzowsky D, Lang D, Reski R: Protein encoding genes in an ancient plant: analysis of codon usage, retained genes and splice sites in a moss, Physcomitrella patens.

    BMC Genomics 2005, 6:43. PubMed Abstract | BioMed Central Full Text | PubMed Central Full Text OpenURL

  3. Lang D, Zimmer AD, Rensing SA, Reski R: Exploring plant biodiversity: the Physcomitrella genome and beyond.

    Trends Plant Sci 2008, 13:542-549. PubMed Abstract | Publisher Full Text OpenURL

  4. cosmoss.org http://www.cosmoss.org webcite.

  5. Saidi Y, Finka A, Chakhporanian M, Zryd JP, Schaefer DG, Goloubinoff P: Controlled expression of recombinant proteins in Physcomitrella patens by a conditional heat-shock promoter: a tool for plant research and biotechnology.

    Plant Mol Biol 2005, 59:697-711. PubMed Abstract | Publisher Full Text OpenURL

  6. Lucumi A, Posten C, Pons MN: Image analysis supported moss cell disruption in photo-bioreactors.

    Plant Biol (Stuttg) 2005, 7:276-282. Publisher Full Text OpenURL

  7. Decker EL, Reski R: Moss bioreactors producing improved biopharmaceuticals.

    Curr Opin Biotechnol 2007, 18:393-398. PubMed Abstract | Publisher Full Text OpenURL

  8. Büttner-Mainik A, Parsons J, Jérôme H, Hartmann A, Lamer S, Schaaf A, Schlosser A, Zipfel PF, Reski R, Decker EL: Production of biologically active recombinant human factor H in Physcomitrella.

    Plant Biotechnol J 2011, 9:373-383. PubMed Abstract | Publisher Full Text OpenURL

  9. Parsons J, Altmann F, Arrenberg CK, Koprivova A, Beike AK, Stemmer C, Gorr G, Reski R, Decker EL: Moss-based production of asialo-erythropoietin devoid of Lewis A and other plant-typical carbohydrate determinants.

    Plant Biotechnol J 2012, 10:851-861. PubMed Abstract | Publisher Full Text OpenURL

  10. Rensing SA, Lang D, Zimmer AD, Terry A, Salamov A, Shapiro H, Nishiyama T, Perroud PF, Lindquist EA, Kamisugi Y, Tanahashi T, Sakakibara K, Fujita T, Oishi K, Shin-I T, Kuroki Y, Toyoda A, Suzuki Y, Hashimoto SI, Yamaguchi K, Sugano S, Kohara Y, Fujiyama A, Anterola A, Aoki S, Ashton N, Barbazuk WB, Barker E, Bennetzen JL, Blankenship R, et al.: The Physcomitrella genome reveals evolutionary insights into the conquest of land by plants.

    Science 2008, 319:64-69. PubMed Abstract | Publisher Full Text OpenURL

  11. Hirano K, Nakajima M, Asano K, Nishiyama T, Sakakibara H, Kojima M, Katoh E, Xiang H, Tanahashi T, Hasebe M, Banks JA, Ashikari M, Kitano H, Ueguchi-Tanaka M, Matsuoka M: The GID1-Mediated gibberellin perception mechanism is conserved in the Lycophyte Selaginella moellendorffii but not in the Bryophyte Physcomitrella patens.

    Plant Cell 2007, 26:26. OpenURL

  12. Axtell MJ, Snyder JA, Bartel DP: Common functions for diverse small RNAs of land plants.

    Plant Cell 2007, 19:1750-1769. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  13. Alboresi A, Caffarri S, Nogue F, Bassi R, Morosinotto T: In silico and biochemical analysis of Physcomitrella patens photosynthetic antenna: identification of subunits which evolved upon land adaptation.

    PLoS One 2008, 3:e2033. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  14. Gitzinger M, Parsons J, Reski R, Fussenegger M: Functional cross-kingdom conservation of mammalian and moss (Physcomitrella patens) transcription, translation and secretion machineries.

    Plant Biotechnol J 2009, 7:73-86. PubMed Abstract | Publisher Full Text OpenURL

  15. Paponov IA, Teale W, Lang D, Paponov M, Reski R, Rensing SA, Palme K: The evolution of nuclear auxin signalling.

    BMC Evol Biol 2009, 9:126. PubMed Abstract | BioMed Central Full Text | PubMed Central Full Text OpenURL

  16. Pils B, Heyl A: Unraveling the evolution of cytokinin signaling.

    Plant Physiol 2009, 151:782-791. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  17. Mukherjee K, Brocchieri L, Bürglin TR: A comprehensive classification and evolutionary analysis of plant homeobox genes.

    Mol Biol Evol 2009, 26:2775-2794. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  18. Bartels S, González-Besteiro MA, Lang D, Ulm R: Emerging functions for plant MAP kinase phosphatases.

    Trends Plant Sci 2010, 15:322-329. PubMed Abstract | Publisher Full Text OpenURL

  19. Pitsch NT, Witsch B, Baier M: Comparison of the chloroplast peroxidase system in the chlorophyte Chlamydomonas reinhardtii, the bryophyte Physcomitrella patens, the lycophyte Selaginella moellendorffii and the seed plant Arabidopsis thaliana.

    BMC Plant Biol 2010, 10:133. PubMed Abstract | BioMed Central Full Text | PubMed Central Full Text OpenURL

  20. Aya K, Hiwatashi Y, Kojima M, Sakakibara H, Ueguchi-Tanaka M, Hasebe M, Matsuoka M: The Gibberellin perception system evolved to regulate a pre-existing GAMYB-mediated system during land plant evolution.

    Nat Commun 2011, 2:544. PubMed Abstract | Publisher Full Text OpenURL

  21. Yasumura Y, Moylan EC, Langdale JA: A conserved transcription factor mediates nuclear control of organelle biogenesis in anciently diverged land plants.

    Plant Cell 2005, 17:1894-1907. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  22. Menand B, Yi K, Jouannic S, Hoffmann L, Ryan E, Linstead P, Schaefer DG, Dolan L: An ancient mechanism controls the development of cells with a rooting function in land plants.

    Science 2007, 316:1477-1480. PubMed Abstract | Publisher Full Text OpenURL

  23. Marella HH, Sakata Y, Quatrano RS: Characterization and functional analysis of ABSCISIC ACID INSENSITIVE3-like genes from Physcomitrella patens.

    Plant J 2006, 46:1032-1044. PubMed Abstract | Publisher Full Text OpenURL

  24. Perroud PF, Quatrano RS: BRICK1 is required for apical cell growth in filaments of the moss Physcomitrella patens but not for gametophore morphology.

    Plant Cell 2008, 20:411-422. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  25. Zemach A, McDaniel IE, Silva P, Zilberman D: Genome-wide evolutionary analysis of eukaryotic DNA methylation.

    Science 2010, 328:916-919. PubMed Abstract | Publisher Full Text OpenURL

  26. Heidel AJ, Lawal HM, Felder M, Schilde C, Helps NR, Tunggal B, Rivero F, John U, Schleicher M, Eichinger L, Platzer M, Noegel AA, Schaap P, Glöckner G: Phylogeny-wide analysis of social amoeba genomes highlights ancient origins for complex intercellular communication.

    Genome Res 2011, 21:1882-1891. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  27. Yang X, Tschaplinski TJ, Hurst GB, Jawdy S, Abraham PE, Lankford PK, Adams RM, Shah MB, Hettich RL, Lindquist E, Kalluri UC, Gunter LE, Pennacchio C, Tuskan GA: Discovery and annotation of small proteins using genomics, proteomics, and computational approaches.

    Genome Res 2011, 21:634-641. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  28. Fu Y, Bannach O, Chen H, Teune JH, Schmitz A, Steger G, Xiong L, Barbazuk WB: Alternative splicing of anciently exonized 5S rRNA regulates plant transcription factor TFIIIA.

    Genome Res 2009, 19:913-921. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  29. Tang H, Wang X, Bowers JE, Ming R, Alam M, Paterson AH: Unraveling ancient hexaploidy through multiply-aligned angiosperm gene maps.

    Genome Res 2008, 18:1944-1954. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  30. JGI plant flagship genomes http://hagsc.org/projects/flagship webcite

  31. Van den Berg BHJ, McCarthy FM, Lamont SJ, Burgess SC: Re-annotation is an essential step in systems biology modeling of functional genomics data.

    PLoS One 2010, 5:e10642. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  32. Markmann-Mulisch U, Wendeler E, Zobell O, Schween G, Steinbiss HH, Reiss B: Differential requirements for RAD51 in Physcomitrella patens and Arabidopsis thaliana development and DNA damage repair.

    Plant Cell 2007, 19:3080-3089. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  33. Liang C, Jaiswal P, Hebbard C, Avraham S, Buckler ES, Casstevens T, Hurwitz B, McCouch S, Ni J, Pujar A, Ravenscroft D, Ren L, Spooner W, Tecle I, Thomason J, Tung C, Wei X, Yap I, Youens-Clark K, Ware D, Stein L: Gramene: a growing plant comparative genomics resource.

    Nucleic Acids Res 2008, 36:D947-D953. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  34. Swarbreck D, Wilks C, Lamesch P, Berardini TZ, Garcia-Hernandez M, Foerster H, Li D, Meyer T, Muller R, Ploetz L, Radenbaugh A, Singh S, Swing V, Tissier C, Zhang P, Huala E: The Arabidopsis Information Resource (TAIR): gene structure and function annotation.

    Nucleic Acids Res 2008, 36:D1009-D1014. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  35. Tweedie S, Ashburner M, Falls K, Leyland P, McQuilton P, Marygold S, Millburn G, Osumi-Sutherland D, Schroeder A, Seal R, Zhang H: FlyBase: enhancing Drosophila Gene Ontology annotations.

    Nucleic Acids Res 2009, 37:D555-D559. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  36. Hirschman J, Berardini TZ, Drabkin HJ, Howe D: A MOD(ern) perspective on literature curation.

    Mol Genet Genomics 2010, 283:415-425. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  37. Consortium UP: Reorganizing the protein space at the Universal Protein Resource (UniProt).

    Nucleic Acids Res 2012, 40:D71-D75. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  38. Flicek P, Amode MR, Barrell D, Beal K, Brent S, Carvalho-Silva D, Clapham P, Coates G, Fairley S, Fitzgerald S, Gil L, Gordon L, Hendrix M, Hourlier T, Johnson N, Kähäri AK, Keefe D, Keenan S, Kinsella R, Komorowska M, Koscielny G, Kulesha E, Larsson P, Longden I, McLaren W, Muffato M, Overduin B, Pignatelli M, Pritchard B, Riat HS, et al.: Ensembl 2012.

    Nucleic Acids Res 2012, 40:D84-D90. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  39. Goodstein DM, Shu S, Howson R, Neupane R, Hayes RD, Fazo J, Mitros T, Dirks W, Hellsten U, Putnam N, Rokhsar DS: Phytozome: a comparative platform for green plant genomics.

    Nucleic Acids Res 2012, 40:D1178-D1186. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  40. Van Bel M, Proost S, Wischnitzki E, Movahedi S, Scheerlinck C, Van de Peer Y, Vandepoele K: Dissecting plant genomes with the PLAZA comparative genomics platform.

    Plant Physiol 2012, 158:590-600. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  41. Gaudet P, Lane L, Fey P, Bridge A, Poux S, Auchincloss A, Axelsen K, Braconi Quintaje S, Boutet E, Brown P, Coudert E, Datta RS, de Lima WC, de Oliveira Lima T, Duvaud S, Farriol-Mathis N, Ferro Rojas S, Feuermann M, Gateau A, Hinz U, Hulo C, James J, Jimenez S, Jungo F, Keller G, Lemercier P, Lieberherr D, Moinat M, Nikolskaya A, Pedruzzi I, Rivoire C, Roechert B, Schneider M, Stanley E, Tognolli M, Sjölander K, Bougueleret L, Chisholm RL, Bairoch A: Collaborative annotation of genes and proteins between UniProtKB/Swiss-Prot and dictyBase.

    Database (Oxford) 2009, 2009:bap016. Publisher Full Text OpenURL

  42. Kersey PJ, Lawson D, Birney E, Derwent PS, Haimel M, Herrero J, Keenan S, Kerhornou A, Koscielny G, Kähäri A, Kinsella RJ, Kulesha E, Maheswari U, Megy K, Nuhn M, Proctor G, Staines D, Valentin F, Vilella AJ, Yates A: Ensembl genomes: extending Ensembl across the taxonomic space.

    Nucleic Acids Res 2010, 38:D563-D569. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  43. Howe D, Costanzo M, Fey P, Gojobori T, Hannick L, Hide W, Hill DP, Kania R, Schaeffer M, St Pierre S, Twigger S, White O, Rhee SY: Big data: the future of biocuration.

    Nature 2008, 455:47-50. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  44. Mazumder R, Natale DA, Julio JAE, Yeh LS, Wu CH: Community annotation in biology.

    Biol Direct 2010, 5:12. PubMed Abstract | BioMed Central Full Text | PubMed Central Full Text OpenURL

  45. Richardt S, Lang D, Reski R, Frank W, Rensing SA: PlanTAPDB, a phylogeny-based resource of plant transcription-associated proteins.

    Plant Physiol 2007, 143:1452-1466. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  46. Mitschke J, Fuss J, Blum T, Höglund A, Reski R, Kohlbacher O, Rensing SA: Prediction of dual protein targeting to plant organelles.

    New Phytol 2009, 183:224-235. PubMed Abstract | Publisher Full Text OpenURL

  47. Kamisugi Y, Von Stackelberg M, Lang D, Care M, Reski R, Rensing SA, Cuming AC: A sequence-anchored genetic linkage map for the moss, Physcomitrella patens.

    Plant J 2008, 56:855-866. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  48. Stein LD, Mungall C, Shu S, Caudy M, Mangone M, Day A, Nickerson E, Stajich JE, Harris TW, Arva A, Lewis S: The generic genome browser: a building block for a model organism system database.

    Genome Res 2002, 12:1599-1610. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  49. Benson DA, Karsch-Mizrachi I, Lipman DJ, Ostell J, Sayers EW: GenBank.

    Nucleic Acids Res 2011, 39:D32-D37. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  50. Punta M, Coggill PC, Eberhardt RY, Mistry J, Tate J, Boursnell C, Pang N, Forslund K, Ceric G, Clements J, Heger A, Holm L, Sonnhammer ELL, Eddy SR, Bateman A, Finn RD: The Pfam protein families database.

    Nucleic Acids Res 2012, 40:D290-D301. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  51. Kozomara A, Griffiths-Jones S: miRBase: integrating microRNA annotation and deep-sequencing data.

    Nucleic Acids Res 2011, 39:D152-D157. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  52. Axtell MJ, Jan C, Rajagopalan R, Bartel DP: A two-hit trigger for siRNA biogenesis in plants.

    Cell 2006, 127:565-577. PubMed Abstract | Publisher Full Text OpenURL

  53. Griffiths-Jones S: The microRNA Registry.

    Nucleic Acids Res 2004, 32:D109-D111. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  54. Kodama Y, Shumway M, Leinonen R: The Sequence Read Archive: explosive growth of sequencing data.

    Nucleic Acids Res 2012, 40:D54-D56. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  55. Duvick J, Fu A, Muppirala U, Sabharwal M, Wilkerson MD, Lawrence CJ, Lushbough C, Brendel V: PlantGDB: a resource for comparative plant genomics.

    Nucleic Acids Res 2008, 36:D959-D965. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  56. Novikova O, Smyshlyaev G, Blinov A: Evolutionary genomics revealed interkingdom distribution of Tcn1-like chromodomain-containing Gypsy LTR retrotransposons among fungi and plants.

    BMC Genomics 2010, 11:231. PubMed Abstract | BioMed Central Full Text | PubMed Central Full Text OpenURL

  57. Degroeve S, Saeys Y, De Baets B, Rouze P, Van de Peer Y: SpliceMachine: predicting splice sites from high-dimensional local context representations.

    Bioinformatics 2005, 21:1332-1338. PubMed Abstract | Publisher Full Text OpenURL

  58. Schiex T, Moisan A, Rouzé P: Computational Biology. In EuGène: an Eukaryotic Gene Finder that combines several sources of evidence. 2066 edition. Edited by Gascuel O, Sagot MF. Berlin / Heidelberg: Springer; 2001:111-125. PubMed Abstract | Publisher Full Text OpenURL

  59. Haas BJ, Salzberg SL, Zhu W, Pertea M, Allen JE, Orvis J, White O, Buell CR, Wortman JR: Automated eukaryotic gene structure annotation using EVidenceModeler and the Program to Assemble Spliced Alignments.

    Genome Biol 2008, 9:R7. PubMed Abstract | BioMed Central Full Text | PubMed Central Full Text OpenURL

  60. Haas BJ, Delcher AL, Mount SM, Wortman JR, Smith RK, Hannick LI, Maiti R, Ronning CM, Rusch DB, Town CD, Salzberg SL, White O: Improving the Arabidopsis genome annotation using maximal transcript alignment assemblies.

    Nucleic Acids Res 2003, 31:5654-5666. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  61. Gremme G, Brendel V, Sparks ME, Kurtz S: Engineering a Software Tool for Gene Structure Prediction in Higher Organisms.

    Inf Softw Technol 2005, 47:965-978. Publisher Full Text OpenURL

  62. Wu TD, Watanabe CK: GMAP: a genomic mapping and alignment program for mRNA and EST sequences.

    Bioinformatics 2005, 21:1859-1875. PubMed Abstract | Publisher Full Text OpenURL

  63. Wang BB, Brendel V: Genomewide comparative analysis of alternative splicing in plants.

    Proc Natl Acad Sci USA 2006, 103:7175-7180. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  64. Barbazuk WB, Fu Y, McGinnis KM: Genome-wide analyses of alternative splicing in plants: opportunities and challenges.

    Genome Res 2008, 18:1381-1392. PubMed Abstract | Publisher Full Text OpenURL

  65. Sinha R, Zimmer AD, Bolte K, Lang D, Reski R, Platzer M, Rensing SA, Backofen R: Identification and characterization of NAGNAG alternative splicing in the moss Physcomitrella patens.

    BMC Plant Biol 2010, 10:76. PubMed Abstract | BioMed Central Full Text | PubMed Central Full Text OpenURL

  66. Filichkin SA, Priest HD, Givan SA, Shen R, Bryant DW, Fox SE, Wong WK, Mockler TC: Genome-wide mapping of alternative splicing in Arabidopsis thaliana.

    Genome Res 2010, 20:45-58. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  67. Alexandrov NN, Troukhan ME, Brover VV, Tatarinova T, Flavell RB, Feldmann KA: Features of Arabidopsis genes and genome discovered using full-length cDNAs.

    Plant Mol Biol 2006, 60:69-85. PubMed Abstract | Publisher Full Text OpenURL

  68. Lunde C, Baumann U, Shirley NJ, Drew DP, Fincher GB: Gene structure and expression pattern analysis of three monodehydroascorbate reductase (Mdhar) genes in Physcomitrella patens: implications for the evolution of the MDHAR family in plants.

    Plant Mol Biol 2006, 60:259-275. PubMed Abstract | Publisher Full Text OpenURL

  69. Merchant SS, Prochnik SE, Vallon O, et al.: The Chlamydomonas genome reveals the evolution of key animal and plant functions.

    Science 2007, 318:245-250. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  70. Chung BYW, Simons C, Firth AE, Brown CM, Hellens RP: Effect of 5’UTR introns on gene expression in Arabidopsis thaliana.

    BMC Genomics 2006, 7:120. PubMed Abstract | BioMed Central Full Text | PubMed Central Full Text OpenURL

  71. Weise A, Rodriguez-Franco M, Timm B, Hermann M, Link S, Jost W, Gorr G: Use of Physcomitrella patens actin 5’ regions for high transgene expression: importance of 5’ introns.

    Appl Microbiol Biotechnol 2006, 70:337-345. PubMed Abstract | Publisher Full Text OpenURL

  72. Mignone F, Gissi C, Liuni S, Pesole G: Untranslated regions of mRNAs.

    Genome Biol 2002., 3

    REVIEWS0004

    OpenURL

  73. Conesa A, Götz S, García-Gómez JM, Terol J, Talón M, Robles M: Blast2GO: a universal tool for annotation, visualization and analysis in functional genomics research.

    Bioinformatics 2005, 21:3674-3676. PubMed Abstract | Publisher Full Text OpenURL

  74. Quevillon E, Silventoinen V, Pillai S, Harte N, Mulder N, Apweiler R, Lopez R: InterProScan: protein domains identifier.

    Nucleic Acids Res 2005, 33:W116-W120. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  75. Gene Ontology (GO) evidence codeshttp://www.geneontology.org/GO.evidence.shtml webcite

  76. Huang DW, Sherman BT, Lempicki RA: Bioinformatics enrichment tools: paths toward the comprehensive functional analysis of large gene lists.

    Nucleic Acids Res 2009, 37:1-13. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  77. Wolf L, Rizzini L, Stracke R, Ulm R, Rensing SA: The molecular and physiological responses of Physcomitrella patens to ultraviolet-B radiation.

    Plant Physiol 2010, 153:1123-1134. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  78. Wood AJ, Duff RJ: The aldehyde dehydrogenase (ALDH) gene superfamily of the moss Physcomitrella patens and the algae Chlamydomonas reinhardtii and Ostreococcus tauri.

    Bryologist 2009, 112:1-11. Publisher Full Text OpenURL

  79. Li L, Stoeckert CJ, Roos DS: OrthoMCL: identification of ortholog groups for eukaryotic genomes.

    Genome Res 2003, 13:2178-2189. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  80. Kersting AR, Bornberg-Bauer E, Moore AD, Grath S: Dynamics and adaptive benefits of protein domain emergence and arrangements during plant genome evolution.

    Genome Biol Evol 2012, 4:316-329. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  81. Teich R, Grauvogel C, Petersen J: Intron distribution in Plantae: 500 million years of stasis during land plant evolution.

    Gene 2007, 394:96-104. PubMed Abstract | Publisher Full Text OpenURL

  82. Roy SW, Penny D: Patterns of intron loss and gain in plants: intron loss-dominated evolution and genome-wide comparison of O. sativa and A. thaliana.

    Mol Biol Evol 2007, 24:171-181. PubMed Abstract | Publisher Full Text OpenURL

  83. Frugoli JA, McPeek MA, Thomas TL, McClung CR: Intron loss and gain during evolution of the catalase gene family in angiosperms.

    Genetics 1998, 149:355-365. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  84. Hu K: Intron exclusion and the mystery of intron loss.

    FEBS Lett 2006, 580:6361-6365. PubMed Abstract | Publisher Full Text OpenURL

  85. Roy SW, Gilbert W: The evolution of spliceosomal introns: patterns, puzzles and progress.

    Nat Rev Genet 2006, 7:211-221. PubMed Abstract | Publisher Full Text OpenURL

  86. Funke RP, Kovar JL, Logsdon JM Jr, Corrette-Bennett JC, Straus DR, Weeks DP: Nucleus-encoded, plastid-targeted acetolactate synthase genes in two closely related chlorophytes, Chlamydomonas reihardtii and Volvox carteri: phylogenetic origins and recent insertion of introns.

    Mol Gen Genet 1999, 262:12-21. PubMed Abstract | Publisher Full Text OpenURL

  87. Rensing SA, Ick J, Fawcett JA, Lang D, Zimmer A, De Peer YV, Reski R: An ancient genome duplication contributed to the abundance of metabolic genes in the moss Physcomitrella patens.

    BMC Evol Biol 2007, 7:130. PubMed Abstract | BioMed Central Full Text | PubMed Central Full Text OpenURL

  88. Banks JA, Nishiyama T, Hasebe M, et al.: The Selaginella genome identifies genetic changes associated with the evolution of vascular plants.

    Science 2011, 332:960-963. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  89. Schaefer D, Zryd JP, Knight CD, Cove DJ: Stable transformation of the moss Physcomitrella patens.

    Mol Gen Genet 1991, 226:418-424. PubMed Abstract OpenURL

  90. Hohe A, Reski R: A tool for understanding homologous recombination in plants.

    Plant Cell Rep 2003, 21:1135-1142. PubMed Abstract | Publisher Full Text OpenURL

  91. Kamisugi Y, Schlink K, Rensing SA, Schween G, Von Stackelberg M, Cuming AC, Reski R, Cove DJ: The mechanism of gene targeting in Physcomitrella patens: homologous recombination, concatenation and multiple integration.

    Nucleic Acids Res 2006, 34:6206-6214. OpenURL

  92. P. patens annotation guidelines https://www.cosmoss.org/physcome_project/wiki/Annotation_guidelines webcite

  93. Cosmoss.org gene annotation interface – Genonaut http://www.cosmoss.org/annotation/Genonaut webcite

  94. Ed L, Nomi H, Mark G, Raymond C, Suzanna L: Apollo: a community resource for genome annotation editing.

    Bioinformatics 2009, 25:1836-1837. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  95. Florea L, Hartzell G, Zhang Z, Rubin GM, Miller W: A computer program for aligning a cDNA sequence with a genomic DNA sequence.

    Genome Res 1998, 8:967-974. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  96. Slater GSC, Birney E: Automated generation of heuristics for biological sequence comparison.

    BMC Bioinforma 2005, 6:31. BioMed Central Full Text OpenURL

  97. PASA software homepage http://pasa.sourceforge.net webcite

  98. Ogata H, Goto S, Sato K, Fujibuchi W, Bono H, Kanehisa M: KEGG: Kyoto Encyclopedia of genes and genomes.

    Nucleic Acids Res 1999, 27:29-34. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  99. Moriya Y, Itoh M, Okuda S, Yoshizawa AC, Kanehisa M: KAAS: an automatic genome annotation and pathway reconstruction server.

    Nucleic Acids Res 2007, 35:W182-W185. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  100. Lang D, Weiche B, Timmerhaus G, Richardt S, Riaño-Pachón DM, Corrêa LGG, Reski R, Mueller-Roeber B, Rensing SA: Genome-wide phylogenetic comparative analysis of plant transcriptional regulation: a timeline of loss, gain, expansion and correlation with complexity.

    Genome Biol Evol 2010, 2:488-503. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  101. Höglund A, Donnes P, Blum T, Adolph HW, Kohlbacher O: MultiLoc: prediction of protein subcellular localization using N-terminal targeting sequences, sequence motifs and amino acid composition.

    Bioinformatics 2006, 22:1158-1165. PubMed Abstract | Publisher Full Text OpenURL

  102. Horton P, Park KJ, Obayashi T, Fujita N, Harada H, Adams-Collier CJ, Nakai K: WoLF PSORT: protein localization predictor.

    Nucleic Acids Res 2007, 35:W585-W587. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  103. Emanuelsson O, Nielsen H, Brunak S, Von Heijne G: Predicting subcellular localization of proteins based on their N-terminal amino acid sequence.

    J Mol Biol 2000, 300:1005-1016. PubMed Abstract | Publisher Full Text OpenURL

  104. Emanuelsson O, Nielsen H, Von Heijne G: ChloroP, a neural network-based method for predicting chloroplast transit peptides and their cleavage sites.

    Protein Sci 1999, 8:978-984. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  105. Boden M, Hawkins J: Prediction of subcellular localization using sequence-biased recurrent networks.

    Bioinformatics 2005, 21:2279-2286. PubMed Abstract | Publisher Full Text OpenURL

  106. Hulo N, Bairoch A, Bulliard V, Cerutti L, De Castro E, Langendijk-Genevaux PS, Pagni M, Sigrist CJ: The PROSITE database.

    Nucleic Acids Res 2006, 34:D227-D230. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  107. Tusnády GE, Simon I: The HMMTOP transmembrane topology prediction server.

    Bioinformatics 2001, 17:849-850. PubMed Abstract | Publisher Full Text OpenURL

  108. Jones DT: Improving the accuracy of transmembrane protein topology prediction using evolutionary information.

    Bioinformatics 2007, 23:538-544. PubMed Abstract | Publisher Full Text OpenURL

  109. Alexa A, Rahnenführer J, Lengauer T: Improved scoring of functional groups from gene expression data by decorrelating GO graph structure.

    Bioinformatics 2006, 22:1600-1607. PubMed Abstract | Publisher Full Text OpenURL

  110. Katoh K, Kuma K, Toh H, Miyata T: MAFFT version 5: improvement in accuracy of multiple sequence alignment.

    Nucleic Acids Res 2005, 33:511-518. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  111. Rice P, Longden I, Bleasby A: EMBOSS: the European Molecular Biology Open Software Suite.

    Trends Genet 2000, 16:276-277. PubMed Abstract | Publisher Full Text OpenURL

  112. Schattner P, Brooks AN, Lowe TM: The tRNAscan-SE, snoscan and snoGPS web servers for the detection of tRNAs and snoRNAs.

    Nucleic Acids Res 2005, 33:W686-W689. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  113. Laslett D, Canback B: ARAGORN, a program to detect tRNA genes and tmRNA genes in nucleotide sequences.

    Nucleic Acids Res 2004, 32:11-16. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  114. Lagesen K, Hallin P, Rødland EA, Staerfeldt HH, Rognes T, Ussery DW: RNAmmer: consistent and rapid annotation of ribosomal RNA genes.

    Nucleic Acids Res 2007, 35:3100-3108. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  115. Pruesse E, Quast C, Knittel K, Fuchs BM, Ludwig W, Peplies J, Glöckner FO: SILVA: a comprehensive online resource for quality checked and aligned ribosomal RNA sequence data compatible with ARB.

    Nucleic Acids Res 2007, 35:7188-7196. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  116. Nawrocki EP, Kolbe DL, Eddy SR: Infernal 1.0: inference of RNA alignments.

    Bioinformatics 2009, 25:1335-1337. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  117. Gardner PP, Daub J, Tate J, Moore BL, Osuch IH, Griffiths-Jones S, Finn RD, Nawrocki EP, Kolbe DL, Eddy SR, Bateman A: Rfam: wikipedia, clans and the “decimal” release.

    Nucleic Acids Res 2011, 39:D141-D145. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  118. Stajich JE, Block D, Boulez K, Brenner SE, Chervitz SA, Dagdigian C, Fuellen G, Gilbert JG, Korf I, Lapp H, Lehvaslaiho H, Matsalla C, Mungall CJ, Osborne BI, Pocock MR, Schattner P, Senger M, Stein LD, Stupka E, Wilkinson MD, Birney E: The Bioperl toolkit: Perl modules for the life sciences.

    Genome Res 2002, 12:1611-1618. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  119. Jurka J: Repbase update: a database and an electronic journal of repetitive elements.

    Trends Genet 2000, 16:418-420. PubMed Abstract | Publisher Full Text OpenURL

  120. RepeatMasker Open-3.0 http://www.repeatmasker.org webcite

  121. codonW: http://codonw.sourceforge.net

  122. P. patens V1.6 annotation download section https://www.cosmoss.org/physcome_project/wiki/Genome_Annotation/V1.6#Downloads_genome_annotation_version_1.6 webcite