Email updates

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

Open Access Research article

Back to the sea twice: identifying candidate plant genes for molecular evolution to marine life

Lothar Wissler1, Francisco M Codoñer2, Jenny Gu1, Thorsten BH Reusch3, Jeanine L Olsen4, Gabriele Procaccini5* and Erich Bornberg-Bauer1*

Author Affiliations

1 Evolutionary Bioinformatics, Institute for Evolution and Biodiversity, University of Muenster, Huefferstrasse 1, D48149 Muenster, Germany

2 IRSI-Caixa Fundation, Hospital Universitari Germans Trias I Pujol, Crta de Canyet s/n 08916 Badalona, Spain

3 Leibniz-Institut fuer Meereswissenschaften IFM-Geomar, Duesternbrooker Weg 20, D24105 Kiel, Germany

4 Dept. Marine Benthic Ecology and Evolution, Centre for Ecological and Evolutionary Studies, University of Groningen, Centre for Life Sciences, Nijenborgh 7, 9747 AG Groningen, The Netherlands

5 Stazione Zoologica A Dohrn, Villa Comunale, 80121 Naples, Italy

For all author emails, please log on.

BMC Evolutionary Biology 2011, 11:8  doi:10.1186/1471-2148-11-8


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


Received:2 June 2010
Accepted:12 January 2011
Published:12 January 2011

© 2011 Wissler 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

Seagrasses are a polyphyletic group of monocotyledonous angiosperms that have adapted to a completely submerged lifestyle in marine waters. Here, we exploit two collections of expressed sequence tags (ESTs) of two wide-spread and ecologically important seagrass species, the Mediterranean seagrass Posidonia oceanica (L.) Delile and the eelgrass Zostera marina L., which have independently evolved from aquatic ancestors. This replicated, yet independent evolutionary history facilitates the identification of traits that may have evolved in parallel and are possible instrumental candidates for adaptation to a marine habitat.

Results

In our study, we provide the first quantitative perspective on molecular adaptations in two seagrass species. By constructing orthologous gene clusters shared between two seagrasses (Z. marina and P. oceanica) and eight distantly related terrestrial angiosperm species, 51 genes could be identified with detection of positive selection along the seagrass branches of the phylogenetic tree. Characterization of these positively selected genes using KEGG pathways and the Gene Ontology uncovered that these genes are mostly involved in translation, metabolism, and photosynthesis.

Conclusions

These results provide first insights into which seagrass genes have diverged from their terrestrial counterparts via an initial aquatic stage characteristic of the order and to the derived fully-marine stage characteristic of seagrasses. We discuss how adaptive changes in these processes may have contributed to the evolution towards an aquatic and marine existence.

Background

Lambers and co-authors summarized the uniqueness of seagrasses as follows: "Aquatic angiosperms are perhaps comparable to whales: They returned to the water, preserving some features of terrestrial organisms" [1]. The monocotyledonous seagrasses represent, in fact, a polyphyletic group of plants that can live underwater in fully marine environments. At least three independent seagrass lineages, but no other angiosperm species, have evolved to a life in the marine environment [2,3].

Seagrasses consist of about 60 species, most of which superficially resemble terrestrial grasses of the family Poaceae in that they have long, narrow leaves and grow in large meadows. Seagrasses belong to the order of Alismatales which includes 11 families of aquatic-freshwater species and 4 families that are fully marine. The marine families include the Posidoniaceae, Zosteraceae, Hydrocharitaceae, and Cymodoceaceae, and have originated in the Cretaceous period [2]. Phylogenetic analysis of members of the entire order, based on the plastid gene encoding for RuBisCO large subunit [4], indicates that the return into the sea occurred at least three times independently through parallel evolution from a common aquatic-freshwater ancestor of terrestrial origin.

Living submerged in an aqueous environment poses many challenges requiring physiological and morphological adaptations that are distinctive from terrestrial angiosperms. For example, the photosynthetic apparatus needs to be modulated to accommodate the changes in light attenuation through the water depth [5]. Consequently, the overall light intensity is decreased and the wavelength composition of sunlight reaching underwater plants is different. Accordingly, seagrasses have one of the highest light requirements among angiosperms [6,7]. One factor contributing to these high light requirements is the reducing sediments to which seagrasses are rooted. These sediments challenge seagrass root tissues with anaerobic conditions since marine sediments are often oxygen deficient. When the internal transport of oxygen from shoot to root tissues is not sufficient, seagrasses can be forced to resort to fermentative metabolism [8,9]. Submergence also exposes organisms to the forces of wave action and tidal currents that effects reproductive functions and reduces the availability of carbon dioxide (CO2). Consequently, seagrasses have evolved to propagate via hydrophilous pollination [10] and rely on carbonic acid and bicarbonate instead of CO2 [11,12]. Specific to marine environments, seagrasses are often exposed to high salt levels and short-term salinity fluctuations in the coastal and estuarine system [13-15]. Increased levels of sodium (Na+) are known to be toxic, partly due to the fact that both Na+ and potassium (K+) have very similar physicochemical properties. Key metabolic processes in the cytoplasm such as enzymatic reactions, protein synthesis, and ribosome functions rely on K+ as a co-factor [16]. An increased level of Na+ creates a competing environment for K+ binding sites and thus decreases efficiency of these processes. Moreover, detrimental effects can propagate from the cytoplasmic compartment into the chloroplasts, leading to a decreased efficiency of photosynthesis which in turn impairs growth [17].

Strikingly, despite their independent evolutionary routes, seagrasses from the three different lineages have evolved many similar morphologies, life history strategies, and breeding systems [3,18]. This indicates that the aquatic habitat imposes novel selection forces that can lead to parallel evolution. For instance, most seagrass species share a secondarily simplified morphology which includes horizontal rhizomes and strap-like leaves originating from a basal meristem. Additionally, seagrasses have been found to share morphological traits that distinguish them from terrestrial plants such as reduced stamen and corolla, and elongated pollen without exine walls [19]. Except for the genus Enhalus with above-surface pollination, all of the 60 seagrass species exhibit true sub-aqueous pollination by means of filiforme pollen (hydrophily; [10]). This adaptation to a marine habitat is thus an example of morphological parallel evolution [20,21].

Identifying genes and cellular processes that may have adaptive contributions to submerged fully marine habitats is therefore of particular interest. By comparing a group of marine angiosperms to terrestrial angiosperms, consequences of specific selection pressures and molecular adaptations can be uncovered. In general, such phenotypic changes can be caused by both changes in gene expression and the primary sequence of encoded proteins. Protein sequences can be strongly conserved whereas changes in their expression pattern can be adaptive (e.g. [22,23]). Conversely, changes in the coding sequence of genes can modify protein structure, function, and efficiency, and therefore can be used to identify evidence for parallel or convergent evolution as successfully demonstrated in recent studies for sequences in plants [24,25], monkeys [26], and fish [27,28].

In this study, the molecular evolution of an identified set of orthologous genes through changes in the coding sequences is investigated to identify candidate genes that may be involved in morphological and physiological adaptations of seagrasses. Gene expression changes as a second mechanism of phenotypic adaptation will not be addressed due to the limitation of the current dataset, although intra-specific analysis of EST libraries between heat-stressed and control treated Zostera marina has previously been conducted [29]. Comparing orthologous gene sequences of two seagrasses and eight terrestrial angiosperm species allows for the inference of sequence evolution and the statistical assessment of synonymous (dS) and non-synonymous (dN) substitution rates, providing insights into molecular adaptation [30,31] of seagrasses. We use EST libraries which were recently developed for two important seagrass species, the Mediterranean seagrass Posidonia oceanica (L.) Delile and the temperate species Zostera marina L. (eelgrass). These seagrass species are two representatives of three currently recognized independent seagrass lineages (Posidoniaceae and Zosteraceae) [4]. Using a molecular evolution approach, the positive selection (dN/dS > 1) of genes along branches leading to each seagrass species was investigated to identify candidate genes in which adaptations allowed for the transition from a terrestrial to an aquatic - and ultimately marine - lifestyle. Estimates of evolutionary distances can be obtained from the timetree database [32], which lists molecular sequence studies that determined that the two seagrass species Z. marina and P. oceanica split 72.5 to 75 million years ago [33-35] and their evolutionary distance to the terrestrial monocots used in this study is estimated at 131 million years [35].

Results

Construction of the dataset

The aim of this study was to investigate the molecular evolution of genes shared between seagrasses and terrestrial angiosperms following the split from the aquatic, last common ancestor (LCA) of seagrasses from the terrestrial monocots. In order to represent two independent seagrass lineages, expressed sequence tag (EST) data were used for Zostera marina and Posidonia oceanica. Orthologous sequences of the two seagrasses were compared to eight terrestrial angiosperm species with a balanced representation of monocot and dicot clades: four monocots including Zea mays [36], Sorghum bicolor [37], Oryza sativa [38] and Brachypodium distachyon [39]; and four dicots including Arabidopsis thaliana [40], Populus trichocarpa [41], Medicago truncatula [42], and Vitis vinifera [43]. Using sequences from all species, orthologous gene clusters (with one sequence per species) could be constructed for 189 genes. The genomes of the moss Physcomitrella patens [44] and green alga Chlamydomonas reinhardtii [45] were not included in this analysis as these species have split from higher plants roughly 600 and 900 million years ago, respectively. Evolutionary distances of this magnitude would have prevented accurate estimations of mutation rates.

Detection of positive selection after the seagrass splits from terrestrial monocots

Using a maximum likelihood framework, the sequence evolution of each gene was evaluated along the species tree (Figure 1A) by estimating the ratio (ω) between the rates of non-synonymous (dN) and synonymous substitutions (dS) in the coding sequence. The parameters used for the analysis were set such that only the alternative hypothesis allows for positive selection in the foreground branch, and a likelihood ratio test (LRT) can determine whether or not the alternative model is a significantly betterfit to the observed sequence alignment than the null model. To each orthologous gene cluster, the branch site test for positive selection in CODEML (test 2, [46]) was applied, using the evolutionary model that allows for a varying ω within the alignment and thus is sensitive towards positive selection limited to a very small number of sites. Testing for positive selection includes running CODEML twice, both with model A (model = 2; NSsites = 2) but with different constraints for the site classes (see Materials and Methods). Three branches abbreviated Po, Zm, and LCA (see Figure 1A) were used as foreground branches in separate tests to identify positive selection after the split of the two seagrass lineages from the terrestrial monocots. For each model, a likelihood score was obtained and a LRT was performed to test for positive selection with p < 0.05. Separate testing of the three branches allowed for rare cases where a gene was inferred to be positively selected in more than one branch. Accordingly, this approach uncovered 65 cases across 51 genes, where the branch-site test for positive selection was significant at least once for the three tested foreground branches (Table 1, p < 0.05).

thumbnailFigure 1. Identification of positively selected genes in two seagrass species and their last common ancestor. (A) Phylogenetic tree of ten plant species among which the molecular evolution of orthologous gene sequences has been analyzed. Positive selection in seagrass evolution has been tested for each of the three highlighted branches Po, Zm and LCA. Divergence times have been obtained from [34,76-80] and the timetree database [32]. (B) Term cloud of over-represented GeneOntology (GO) terms of positively selected genes compared to all tested genes. For each of the three tested branches, enriched GO terms were determined using all other tested genes as a reference as indicated by the different colors. The size of the GO terms is proportional to the p-value obtained in the enrichment test. This procedure creates a representation similar to sequence logos [81], showing enriched annotation terms instead of sequence conservation patterns. A tabular representation of the enriched GO terms can be found in Additional File 5.

Table 1. Genes with evidence for positive selection in seagrasses

Annotation of positively selected genes (PSGs)

Among the 189 tested genes, 51 genes were identified as positively selected genes (PSGs). Using KEGG pathway information, 30 of the 51 PSGs could be associated to at least one pathway. Metabolic pathways, ribosomes, and photosynthesis showed the highest number of associated genes (Table 2), indicating that several components of these pathways have acquired sequence changes after the split of the common ancestor of seagrasses from the terrestrial monocots 130 MYA. For 27 of the 51 PSGs, positive selection was inferred in the branch leading to the last common ancestor of the two seagrass species (branch LCA, Figure 1A). Signals of positive selection in the LCA branch reflect either adaptation before the split of the two seagrass lineages or parallel evolution after their split. In the LCA branch, positive selection has been inferred mostly for ribosomal and metabolic genes (Table 2). Over-representation analysis of GO terms associated with PSGs in the LCA branch revealed only two functional gene classes significantly enriched, including proteins interacting with calmodulin, and proteins located in the thylakoid lumen (Figure 1B).

Table 2. KEGG pathways that are associated to PSGs

Potential lineage-specific positive selection was also detected for 18 and 20 PSGs in Posidonia and Zostera, respectively (Figure 1A, Table 2). Over-representation analysis using the Gene Ontology (GO) revealed that, within this limited sample, positive selection has acted on different functional classes between the three branches under investigation (Figure 1B). In the Zostera lineage, 6 PSGs were identified to be involved in the photosynthesis pathway (ID: 00195), whereas none of these were observed in the Posidonia lineage, suggesting that parts of the two photosystems and the light reaction have undergone Zostera-specfic adaptation (Table 2). Additionally, GO annotation indicates that in Zostera, positive selection has acted on genes responsive to abiotic stimuli and cold (Figure 1B). PSGs in Posidonia were identified to be mostly involved in metabolic processes and biosynthetic pathways. Together, these findings indicate that the two seagrass lineages have diverged substantially on a molecular level despite a seemingly similar habitat. Nonetheless, many signals of positive selection found in the LCA branch also indicate adaptive traits shared by both lineages. These PSGs may have evolved either in their last common ancestor or in parallel after their split.

Discussion

Positively selected genes associated with central biological pathways

Positive selection for 51 genes was detected after the split from terrestrial monocots based on a maximum likelihood approach. Theoretical models based on confirmed biological data have suggested that molecular adaptation is realized to different extents across the proteome and depends on the functional role of each individual protein [47]. In the present analysis, many of the identified PSGs are involved in the central biological pathways of translation, photosynthesis, and glycolysis. These adaptations are possibly associated to the above mentioned Na+ toxicity which seagrasses have likely experienced during their evolution towards the marine environment. To this respect, molecular adaptation of key cellular processes known to be sensitive towards increased ionic levels such as photosynthesis, translation, and selected metabolic enzymes are expected. Considering the importance of these processes for the survival of an organism over short and evolutionary time scales, it is not surprising to identify strong selection pressure shaping genes which increase salt tolerance.

The available dataset allowed only for the investigation of 189 orthologous clusters, equivalent to ~1% of the A. thaliana genome. Since orthologous clusters include only ESTs from both seagrasses, the presented dataset is not an unbiased sample of the genome and is probably enriched for genes that show significant expression levels in both seagrass species. Nevertheless, the presented analysis was able to provide significant partial insights into the molecular evolution of seagrasses. While the limited size of the current dataset leaves room for further investigations, the well described ecology of seagrasses can be utilized to discuss how these PSGs may have contributed to seagrass adaptation to the marine environment.

Molecular evolution for salt tolerance

A number of terrestrial lineages of plants have evolved into aquatic-freshwater hydrophytes and a number of morphological features are shared by both hydrophytes and seagrasses, e.g., the presence of a diffusive boundary layer around the leaves, a photosynthetic epidermis, loss of stomata and the development of aerenchyma (reviewed in [48]). Physiologically, however, seagrasses must cope with high ion concentrations, inefficient carbon uptake and other physical coping mechanisms that are still poorly understood. One of the questions that has to remain open is how exactly do seagrasses deal with the high salinity of the ocean. Seagrasses have been found to harbor increased intracellular levels of Na+ and K+ as compared to terrestrial angiosperm species [49] as well as to other aquatic angiosperms [48]. In general, salt-tolerant plants compensate osmotic and ionic imbalances with increased K+ import and the accumulation of compatible solutes [50,51]. However, genes that are known to facilitate salt tolerance such as the SOS pathway [52,53] were absent from the orthologous gene clusters and could therefore not be investigated. Thus, the mechanism by which seagrasses achieve either a tolerance of higher salinity levels or employ active mechanisms to decrease intracellular Na+ deserved further investigation with more comprehensive sequence and additional expression datasets.

PSG Group 1: Glycolysis

With two fructose-bisphosphate aldolase enzymes and a malate dehydrogenase, the list of PSGs contains three enzymes of the glycolysis pathway. This observation may be particularly significant due to the challenges imposed by the O2 sink created by the reductive sediment leading to compensation by internal transport of oxygen from shoot to root tissues during the day cycle, as mentioned above. In darkness, seagrasses can even be forced to switch to fermentative metabolism. In P. oceanica, malate has previously been shown to accumulate as a consequence of anoxic conditions [54]. Hence, the positive selection of these three glycolysis genes may be associated with seagrass-specific adaptation to anaerobiosis.

PSG Group 2: Ribosomal Genes

Ten PSGs were found to be ribosomal proteins involved in translation. From an evolutionary point of view, translation is an ancient cellular process, and high selection pressure is expected to act against deleterious mutations, as ribosome functioning affects virtually all cellular processes. In A. thaliana, on average four gene copies encode for any of the approximately 80 different ribosomal proteins [55,56]. This redundancy may reflect the importance of maintaining highly productive translation and protein synthesis. At least three scenarios can explain the seemingly high number of PSGs with ribosomal function: (1) Since translation is salt-sensitive, one can hypothesize that these changes reflect salt tolerance adaptations. The vast majority of signals of positive selection in ribosomes were inferred in the LCA branch so that these changes are shared by both seagrass species. Ultimately, signals of positive selection in ribosomes could be one of the traits that allowed the transition to the marine lifestyle. (2) As ribosomes consist of a multitude of subunits, changes in only a few proteins could cause compensatory mutations in other ribosomal proteins to maintain structure and function of the ribosomal complex. Such compensatory mutations were shown to occur in an E. coli mutant [57], and would increase the number of observed changes and overestimate the number of "adaptive changes". (3) Acquisition of non-ribosomal functions could explain sequence changes in these proteins without them being adaptive in the context of ribosomal functioning. In the primate ribosomal protein S4, positive selection has been shown to occur after gene duplication [58]. Andrés et al. [58] concluded that one gene copy has acquired a non-ribosomal function with 2 to 6 amino acid substitutions identified as positively selected sites. The three presented scenarios are not mutually exclusive and ultimately, more experiments will be required to reveal the nature of the inferred sequence changes.

PSG Group 3: Photosynthesis and carbon fixation

Seven PSGs were related to the photosynthetic pathway and may reflect adaptations to new conditions of carbon fixation and photosynthesis that seagrasses had to face after their split from a terrestrial ancestor. Fixation of CO2 is expected to be more difficult for seagrasses since seawater contains very little dissolved carbon dioxide. While CO2 can readily diffuse from the air through the stomata to the mesophyll cells in terrestrial plants, aquatic plants often have limited CO2 diffusion rates [1]. Factors contributing to slow CO2 diffusion in aquatic plants (and especially in seagrasses) are thick boundary layers around the leaves that are sometimes amplified by the presence of unicellular or multicellular photosynthetic epiphytes that compete for CO2 [59], and the low rate of CO2 transport in water. The two seagrass species under investigation, Z. marina and P. oceanica, are known to utilize bicarbonate (<a onClick="popup('http://www.biomedcentral.com/1471-2148/11/8/mathml/M1','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2148/11/8/mathml/M1">View MathML</a>) as a major source of inorganic carbon for photosynthesis [11,12]. The ability to utilize <a onClick="popup('http://www.biomedcentral.com/1471-2148/11/8/mathml/M1','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2148/11/8/mathml/M1">View MathML</a> could be one of the traits evolved in the LCA branch. In contrast, a set of signals of positive selection specific to the Zostera lineage could relate to the biochemical mechanism used in carbon fixation. Seagrasses have long been regarded as C3 plants, but physiological measurements have gathered indications that several seagrass species, including Z. marina, are C3-C4 intermediates or have various carbon-concentrating mechanisms to aid the RuBisCO enzyme in carbon acquisition [60-63]. Finally, seagrasses are able to activate different mechanisms to cope with conditions of light-limitation and shifted light spectrum [6,7] through long-lasting metabolic adjustments including down-regulation of RuBisCO, enhanced proteolysis [64] and putative changes in the antenna complex. These various unique characteristics of seagrasses are further supported by our results.

Conclusions

We have undertaken the first step in systematically unraveling the molecular basis of seagrass evolution from terrestrial ancestors to a fully marine lifestyle. Only genes that were contained in the available seagrass EST collections could be analyzed in this study. Consequently, the current dataset of orthologous gene clusters for 10 angiosperm species is biased and limited in size. Nevertheless, this study has shed light on the molecular evolution of seagrass genes expressed under native conditions in root and leaf tissues. 51 genes showed evidence for positive selection in seagrass branches indicating that photosynthesis, a few metabolic pathways, and ribosomes have strongly diverged after the split of the common ancestor of seagrasses from terrestrial monocots. Further studies will need to address the following questions: (1) How seagrasses have acquired osmoregulatory capacity to tolerate high salinities, (2) how CO2 is fixated, (3) how their photosynthetic apparatus has evolved for under water light harvesting, and (4) under what conditions anaerobiosis takes place. In this regard, comparisons with the aquatic members of the Alismatales will be necessary to distinguish between more general adaptations to the aquatic environment and those that are marine-specific. Finally, the completion of the Zostera marina genome project, currently under way at the Joint Genome Institute (http://www.jgi.doe.gov/ webcite), will be a milestone in providing more comprehensive datasets in the near future to further our understanding of evolution and adaptation of seagrasses and their aquatic relatives.

Methods

Sequence data

Gene sequences from ten angiosperm species were compared to identify genes with signs of positive selection in seagrasses. The two seagrass species Zostera marina and Posidonia oceanica were represented by expressed sequence tag data. Protein-coding sequences from the genomes of eight terrestrial angiosperm species were used to contrast the seagrass sequences. These species included Zea mays, Sorghum bicolor, Oryza sativa, Brachypodium distachyon, Arabidopsis thaliana, Populus trichocarpa, Medicago truncatula, and Vitis vinifera. In the seagrass ESTs representing putative transcript sequences, open reading frames (ORFs) were predicted based on significant BLASTX matches [65] to protein sequences of the other eight angiosperm species (E < 1e-5). Two sequence datasets were constructed: one containing protein sequences, and another one for the protein coding sequences (CDS). For more information on the EST sequences and how the libraries were built can be found in [66].

Orthologous gene clusters

Using the protein sequences of the ten species, orthologous gene clusters were constructed with OrthoMCL [67] using default settings. Only clusters with at least one sequence per species were used in our analysis. If more than one sequence of any species was contained in an OrthoMCL cluster, all sequences of that species were removed except for the one sequence that showed the highest similarity to all other sequences of the cluster as assessed with T-Coffee [68]. For each 1:1 ortholog cluster (see Additional file 1), multiple sequence alignments (MSAs) of the protein and coding sequences were constructed. First, protein sequences were aligned with MUSCLE [69] (see Additional file 2). Second, PAL2NAL [70] was applied to align the CDS codon-wise, guided by the protein MSA as a reference (see Additional file 3).

Additional file 1. OrthoMCL cluster composition after clustering two seagrass EST datasets and 8 full angiosperm genomes.

Format: CSV Size: 36KB Download fileOpen Data

Additional file 2. Aligned proteins of each OrthoMCL cluster, produced with MUSCLE.

Format: GZ Size: 158KB Download fileOpen Data

Additional file 3. Aligned nucleotide (CDS) sequences of each OrthoMCL cluster, produced with PAL2NAL.

Format: GZ Size: 334KB Download fileOpen Data

Test for positive selection

CODEML from the PAML package [71] (v4.3) was used to identify genes under positive selection using a codon-based maximum likelihood method [72]. The phylogenetic relationships between the 10 tested taxa were obtained from NCBI Taxonomy (http://www.ncbi.nlm.nih.gov/Taxonomy/ webcite) and used as reference tree. To test a foreground branch for positive selection, CODEML was run twice, both with model A (model = 2; NSsites = 2) but with different constraints for the site classes as described for test 2 [46]. Model A requires branches in the tree to belong to either foreground or background branch category, where only foreground lineages are allowed to have experienced positive selection (ω > 0). Four classes of sites are assumed in model A: (1) class 0 codons are conserved through the tree, with 0 ¡ ω0 < 1; (2) class 1 codons evolve neutrally with ω1 = 1; (3) class 2a and (4) class 2b codons differ in their selection mode between foreground and background branches. In background branches, 2a codons are conserved with 0 < ω 0 < 1, and 2b codons are neutral with ω1 = 1. In foreground branches and the null hypothesis run, 2a and 2b codons evolve neutrally with fixed ω2 = 1. In foreground branches under the alternative hypothesis, 2a and 2b codons are positively selected with each ω2 > 1. In separate runs, each of the three branches Zm, Po, and LCA were marked as foreground branches and the branch site test for positive selection was applied (see Additional file 4). Positive selection was inferred if the LRT between the scores of the models corresponding to the null and the alternative hypothesis was < 0.05. The p-values were not adjusted for multiple testing for two reasons. First, the presented dataset is relatively small, and given a 5% error rate, only about 3 false positives are to be expected among the 65 significant cases of positive selection. Second, lowering the p-value cutoff makes the test for positive selection a lot more conservative, dismissing genes where positive selection is limited to a very small number of residues.

Additional file 4. CODEML output of the branch-site test for positive selection for each of the three tested seagrass branches.

Format: GZ Size: 5.9MB Download fileOpen Data

GeneOntology and KEGG pathway annotation

The A. thaliana ortholog of each cluster was used to associate additional annotation to the whole ortholog cluster. First, GeneOntology (GO) annotation [73] was obtained for Arabidopsis thaliana from http://www.geneontology.org/ webcite. Both filtered and unfiltered gene associations of A. thaliana (8 Dec 2009 version) were pooled. From these pooled annotations, only non-redundant mappings to genes from the newest Arabidopsis genome release (TAIR9) were kept. Based on the Arabidopsis ortholog contained in each cluster, GO terms were mapped to PSGs. The R package topGO [74] was used to test enrichment of GO annotation terms in these PSGs, using all tested orthologous clusters as reference (see Additional file 5). Enrichment was assessed by Fisher Exact tests as implemented in topGO's classic algorithm treating each GO term as an independent unit. Second, A. thaliana KEGG pathway annotation [75] was obtained from ftp://ftp.genome.jp/pub/kegg/genes/organisms/ath/ webcite, and mapped to the ortholog clusters via the Arabidopsis gene id.

Additional file 5. Results of the GeneOntology (GO) enrichment analysis, testing the PSGs of each branch against all genes tested for positive selection.

Format: CSV Size: 3KB Download fileOpen Data

Authors' contributions

This study was conceived by TBHR and GP and designed by LW, EBB and FMC. LW performed all computations with help from FMC. Data were interpreted by LW, EBB and FMC. The paper was written by LW, JG, GP, JLO, and EBB, and corrected and approved by all authors.

Acknowledgements

We thank Christian Fufezan for helpful discussion. We thank Andrew Moore for his help constructing the term cloud in Figure 1B. LW and EBB acknowledge support from the Volkswagen Foundation. GP and JLO acknowledge funding from FP6 NoE Marine Genomics Europe (EC contract reference: GOCE-CT-2004505403).

References

  1. Lambers H, Chapin FS, Pons TL: Plant Physiological Ecology. Berlin: Springer; 1998. OpenURL

  2. den Hartog C: The seagrasses of the world. Amsterdam: North-Holland Publishing Company; 1970. OpenURL

  3. Waycott M, Procaccini G, Les DH, Reusch TBH: Seagrasses: Biology, Ecology and Conservation. Seagrass Evolution, Ecology and Conservation: A Genetic Perspective. Berlin: Springer; 2006. OpenURL

  4. Les DH, Cleland MA, Waycott M: Phylogenetic Studies in Alismatidae, II: Evolution of Marine Angiosperms (Seagrasses) and Hydrophily.

    Systematic Botany 1997, 22(3):443-463. Publisher Full Text OpenURL

  5. Dalla Via J, Sturmbauer C, Schönweger G, Sötz E, Mathekowitsch S, Stifter M, Rieger R: Light gradients and meadow structure in Posidonia oceanica: ecomorphological and functional correlates.

    Marine Ecology Progress Series 1998, 163:267-278. Publisher Full Text OpenURL

  6. Orth RJ, Carruthers TJB, Dennison WC, Duarte CM, Fourqurean JW, Jr KLH, Hughes AR, Olyarnik S, Williams SL, Kendrick GA, Kenworthy WJ, Short FT, Waycott M: A Global Crisis for Seagrass Ecosystems.

    BioOne 2006, 56(12):987-996. OpenURL

  7. Dennison WC, Orth RJ, Moore KA, Stevenson JC, Carter V, Kollar S, Bergstrom PW, Batiuk RA: Assessing water quality with submersed aquatic vegetation.

    BioScience 1993, 43(2):86-94. Publisher Full Text OpenURL

  8. Terrados J, Duarte CM, Kamp-Nielsen L, Agawin NSR, Gacia E, Lacap D, Fortes MD, Borum J, Lubanski M, Greve T: Are seagrass growth and survival constrained by the reducing conditions of the sediment?

    Aquatic Botany 1999, 65(1-4):175-197. Publisher Full Text OpenURL

  9. Touchette , Burkholder : Overview of the physiological ecology of carbon metabolism in seagrasses.

    J Exp Mar Bio Ecol 2000, 250(1-2):169-205. PubMed Abstract | Publisher Full Text OpenURL

  10. Cox PA: Water-pollinated plants.

    Sci Am 1993, 269(4):50-56. PubMed Abstract | Publisher Full Text OpenURL

  11. Beer S, Rehnberg J: The acquisition of inorganic carbon by the seagrass Zostera marina.

    Aquatic Botany 1997, 56:277-283. Publisher Full Text OpenURL

  12. Invers O, Pérez M, Romero J: Bicarbonate utilization in seagrass photosynthesis: role of carbonic anhydrase in Posidonia oceanica (L.) Delile and Cymodocea nodosa (Ucria) Ascherson.

    J Exp Mar Bio Ecol 1999, 235:125-133. Publisher Full Text OpenURL

  13. Barbour MG: Is any angiosperm an obligate halophyte?

    Am Midl Nat 1970, 84:105-120. Publisher Full Text OpenURL

  14. McConnaughey BH: Introduction to Marine Biology. St. Louis, Missouri: Mosby Company; 1978. OpenURL

  15. Walker DI, McComb AJ: Salinity response of the seagrass Amphibolis antarctica (Labill.) Sonder et Aschers: an experimental validation of field results.

    Aquatic Botany 1990, 36:359-366. Publisher Full Text OpenURL

  16. Marschner H: Mineral Nutrition of Higher Plants. London: Academic Press; 1995. OpenURL

  17. Sudhir P, Murthy SDS: Effects of salt stress on basic processes of photosynthesis.

    Photosynthetica 2004, 42:481-486. Publisher Full Text OpenURL

  18. Philbrick CT, Les DH: Evolution of aquatic angiosperm reproductive systems.

    BioScience 1996, 46:813-826. Publisher Full Text OpenURL

  19. Kuo J, den Hartog C: Seagrasses: Biology, Ecology and Conservation Seagrass Morphology, Anatomy, and Ultrastructure, Dordrecht. Netherlands: Springer; 2006. OpenURL

  20. Doolittle RF: Convergent evolution: the need to be explicit.

    Trends Biochem Sci 1994, 19:15-18. PubMed Abstract | Publisher Full Text OpenURL

  21. Wood TE, Burke JM, Rieseberg LH: Parallel genotypic adaptation: when evolution repeats itself.

    Genetica 2005, 123(1-2):157-170. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  22. Holloway AK, Lawniczak MKN, Mezey JG, Begun DJ, Jones CD: Adaptive gene expression divergence inferred from population genomics.

    PLoS Genet 2007, 3(10):2007-2013. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  23. Fraser HB, Moses AM, Schadt EE: Evidence for widespread adaptive evolution of gene expression in budding yeast.

    Proc Natl Acad Sci USA 2010, 107(7):2977-2982. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  24. Christin PA, Salamin N, Savolainen V, Duvall MR, Besnard G: C4 Photosynthesis evolved in grasses via parallel adaptive genetic changes.

    Curr Biol 2007, 17(14):1241-1247. PubMed Abstract | Publisher Full Text OpenURL

  25. Christin PA, Salamin N, Muasya AM, Roalson EH, Russier F, Besnard G: Evolutionary switch and genetic convergence on rbcL following the evolution of C4 photosynthesis.

    Mol Biol Evol 2008, 25(11):2361-2368. PubMed Abstract | Publisher Full Text OpenURL

  26. Zhang J: Parallel adaptive origins of digestive RNases in Asian and African leaf monkeys.

    Nat Genet 2006, 38(7):819-823. PubMed Abstract | Publisher Full Text OpenURL

  27. Protas ME, Hersey C, Kochanek D, Zhou Y, Wilkens H, Jeffey WR, Zon LI, Borowsky R, Tabin CJ: Genetic analysis of cavefish reveals molecular convergence in the evolution of albinism.

    Nat Genet 2006, 38:107-111. PubMed Abstract | Publisher Full Text OpenURL

  28. Zakon HH, Lu Y, Zwickl DJ, Hillis DM: Sodium channel genes and the evolution of diversity in communication signals of electric fishes: convergent molecular evolution.

    Proc Natl Acad Sci USA 2006, 103(10):3675-3680. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  29. Reusch TBH, Veron AS, Preuss C, Weiner J, Wissler L, Beck A, Klages S, Kube M, Reinhardt R, Bornberg-Bauer E: Comparative analysis of Expressed Sequence Tag (EST) libraries in the seagrass Zostera marina subjected to temperature stress.

    Mar Biotechnol (NY) 2008, 10(3):297-309. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  30. Hurst LD: The Ka/Ks ratio: diagnosing the form of sequence evolution.

    Trends Genet 2002, 18(9):486. PubMed Abstract | Publisher Full Text OpenURL

  31. Yang , Bielawski : Statistical methods for detecting molecular adaptation.

    Trends Ecol Evol 2000, 15(12):496-503. PubMed Abstract | Publisher Full Text OpenURL

  32. Hedges SB, Dudley J, Kumar S: TimeTree: a public knowledge-base of divergence times among organisms.

    Bioinformatics 2006, 22(23):2971-2972. PubMed Abstract | Publisher Full Text OpenURL

  33. Bremer K: Early Cretaceous lineages of monocot flowering plants.

    Proc Natl Acad Sci USA 2000, 97(9):4707-4711. PubMed Abstract | Publisher Full Text OpenURL

  34. Janssen T, Bremer K: The age of major monocot groups inferred from 800+ rbcL sequences.

    Bot J Linn Soc 2004, 146:385-398. Publisher Full Text OpenURL

  35. Anderson C, Janssen T: The Timetree of Life. Monocots, Oxford University Press; 2009. OpenURL

  36. Schnable PS, Ware D, Fulton RS, Stein JC, Wei F, Pasternak S, Liang C, Zhang J, Fulton L, Graves TA, Minx P, Reily AD, Courtney L, Kruchowski SS, Tomlinson C, Strong C, Delehaunty K, Fronick C, Courtney B, Rock SM, Belter E, Du F, Kim K, Abbott RM, Cotton M, Levy A, Marchetto P, Ochoa K, Jackson SM, Gillam B, Chen W, Yan L, Higginbotham J, Cardenas M, Waligorski J, Applebaum E, Phelps L, Falcone J, Kanchi K, Thane T, Scimone A, Thane N, Henke J, Wang T, Ruppert J, Shah N, Rotter K, Hodges J, Ingenthron E, Cordes M, Kohlberg S, Sgro J, Delgado B, Mead K, Chinwalla A, Leonard S, Crouse K, Collura K, Kudrna D, Currie J, He R, Angelova A, Rajasekar S, Mueller T, Lomeli R, Scara G, Ko A, Delaney K, Wissotski M, Lopez G, Campos D, Braidotti M, Ashley E, Golser W, Kim H, Lee S, Lin J, Dujmic Z, Kim W, Talag J, Zuccolo A, Fan C, Sebastian A, Kramer M, Spiegel L, Nascimento L, Zutavern T, Miller B, Ambroise C, Muller S, Spooner W, Narechania A, Ren L, Wei S, Kumari S, Faga B, Levy MJ, McMahan L, Buren PV, Vaughn MW, Ying K, Yeh CT, Emrich SJ, Jia Y, Kalyanaraman A, Hsia AP, Barbazuk WB, Baucom RS, Brutnell TP, Carpita NC, Chaparro C, Chia JM, Deragon JM, Estill JC, Fu Y, Jeddeloh JA, Han Y, Lee H, Li P, Lisch DR, Liu S, Liu Z, Nagel DH, McCann MC, SanMiguel P, Myers AM, Nettleton D, Nguyen J, Penning BW, Ponnala L, Schneider KL, Schwartz DC, Sharma A, Soderlund C, Springer NM, Sun Q, Wang H, Waterman M, Westerman R, Wolfgruber TK, Yang L, Yu Y, Zhang L, Zhou S, Zhu Q, Bennetzen JL, Dawe RK, Jiang J, Jiang N, Presting GG, Wessler SR, Aluru S, Martienssen RA, Clifton SW, McCombie WR, Wing RA, Wilson RK: The B73 maize genome: complexity, diversity, and dynamics.

    Science 2009, 326(5956):1112-1115. PubMed Abstract | Publisher Full Text OpenURL

  37. Paterson AH, Bowers JE, Bruggmann R, Dubchak I, Grimwood J, Gundlach H, Haberer G, Hellsten U, Mitros T, Poliakov A, Schmutz J, Spannagl M, Tang H, Wang X, Wicker T, Bharti AK, Chapman J, Feltus FA, Gowik U, Grigoriev IV, Lyons E, Maher CA, Martis M, Narechania A, Otillar RP, Penning BW, Salamov AA, Wang Y, Zhang L, Carpita NC, Freeling M, Gingle AR, Hash CT, Keller B, Klein P, Kresovich S, McCann MC, Ming R, Peterson DG, ur Rahman M, Ware D, Westhoff P, Mayer KFX, Messing J, Rokhsar DS: The Sorghum bicolor genome and the diversification of grasses.

    Nature 2009, 457(7229):551-556. PubMed Abstract | Publisher Full Text OpenURL

  38. Yu J, Hu S, Wang J, Wong GKS, Li S, Liu B, Deng Y, Dai L, Zhou Y, Zhang X, Cao M, Liu J, Sun J, Tang J, Chen Y, Huang X, Lin W, Ye C, Tong W, Cong L, Geng J, Han Y, Li L, Li W, Hu G, Huang X, Li W, Li J, Liu Z, Li L, Liu J, Qi Q, Liu J, Li L, Li T, Wang X, Lu H, Wu T, Zhu M, Ni P, Han H, Dong W, Ren X, Feng X, Cui P, Li X, Wang H, Xu X, Zhai W, Xu Z, Zhang J, He S, Zhang J, Xu J, Zhang K, Zheng X, Dong J, Zeng W, Tao L, Ye J, Tan J, Ren X, Chen X, He J, Liu D, Tian W, Tian C, Xia H, Bao Q, Li G, Gao H, Cao T, Wang J, Zhao W, Li P, Chen W, Wang X, Zhang Y, Hu J, Wang J, Liu S, Yang J, Zhang G, Xiong Y, Li Z, Mao L, Zhou C, Zhu Z, Chen R, Hao B, Zheng W, Chen S, Guo W, Li G, Liu S, Tao M, Wang J, Zhu L, Yuan L, Yang H: A draft sequence of the rice genome (Oryza sativa L. ssp. indica).

    Science 2002, 296(5565):79-92. PubMed Abstract | Publisher Full Text OpenURL

  39. International Brachypodium Initiative: Genome sequencing and analysis of the model grass Brachypodium distachyon.

    Nature 2010, 463(7282):763-768. PubMed Abstract | Publisher Full Text OpenURL

  40. Arabidopsis Genome Initiative: Analysis of the genome sequence of the flowering plant Arabidopsis thaliana.

    Nature 2000, 408(6814):796-815. PubMed Abstract | Publisher Full Text OpenURL

  41. Tuskan GA, Difazio S, Jansson S, Bohlmann J, Grigoriev I, Hellsten U, Putnam N, Ralph S, Rombauts S, Salamov A, Schein J, Sterck L, Aerts A, Bhalerao RR, Bhalerao RP, Blaudez D, Boerjan W, Brun A, Brunner A, Busov V, Campbell M, Carlson J, Chalot M, Chapman J, Chen GL, Cooper D, Coutinho PM, Couturier J, Covert S, Cronk Q, Cunningham R, Davis J, Degroeve S, Dejardin A, Depamphilis C, Detter J, Dirks B, Dubchak I, Duplessis S, Ehlting J, Ellis B, Gendler K, Goodstein D, Gribskov M, Grimwood J, Groover A, Gunter L, Hamberger B, Heinze B, Helariutta Y, Henrissat B, Holligan D, Holt R, Huang W, Islam-Faridi N, Jones S, Jones-Rhoades M, Jorgensen R, Joshi C, Kangasjärvi J, Karlsson J, Kelleher C, Kirkpatrick R, Kirst M, Kohler A, Kalluri U, Larimer F, Leebens-Mack J, Leple JC, Locascio P, Lou Y, Lucas S, Martin F, Montanini B, Napoli C, Nelson DR, Nelson C, Nieminen K, Nilsson O, Pereda V, Peter G, Philippe R, Pilate G, Poliakov A, Razumovskaya J, Richardson P, Rinaldi C, Ritland K, Rouzé P, Ryaboy D, Schmutz J, Schrader J, Segerman B, Shin H, Siddiqui A, Sterky F, Terry A, Tsai CJ, Uberbacher E, Unneberg P, Vahala J, Wall K, Wessler S, Yang G, Yin T, Douglas C, Marra M, Sandberg G, de Peer YV, Rokhsar D: The genome of black cottonwood, Populus trichocarpa (Torr. & Gray).

    Science 2006, 313(5793):1596-1604. PubMed Abstract | Publisher Full Text OpenURL

  42. Young ND, Cannon SB, Sato S, Kim D, Cook DR, Town CD, Roe BA, Tabata S: Sequencing the genespaces of Medicago truncatula and Lotus japonicus.

    Plant Physiol 2005, 137(4):1174-1181. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  43. Jaillon O, Aury JM, Noel B, Policriti A, Clepet C, Casagrande A, Choisne N, Aubourg S, Vitulo N, Jubin C, Vezzi A, Legeai F, Hugueney P, Dasilva C, Horner D, Mica E, Jublot D, Poulain J, Bruyére C, Billault A, Segurens B, Gouyvenoux M, Ugarte E, Cattonaro F, Anthouard V, Vico V, Fabbro CD, Alaux M, Gaspero GD, Dumas V, Felice N, Paillard S, Juman I, Moroldo M, Scalabrin S, Canaguier A, Clainche IL, Malacrida G, Durand E, Pesole G, Laucou V, Chatelet P, Merdinoglu D, Delledonne M, Pezzotti M, Lecharny A, Scarpelli C, Artiguenave F, Pé ME, Valle G, Morgante M, Caboche M, Adam-Blondon AF, Weissenbach J, Quétier F, Wincker P, for Grapevine Genome Characterization FIPC: The grapevine genome sequence suggests ancestral hexaploidization in major angiosperm phyla.

    Nature 2007, 449(7161):463-467. PubMed Abstract | Publisher Full Text OpenURL

  44. 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, Cho SH, Dutcher SK, Estelle M, Fawcett JA, Gundlach H, Hanada K, Heyl A, Hicks KA, Hughes J, Lohr M, Mayer K, Melkozernov A, Murata T, Nelson DR, Pils B, Prigge M, Reiss B, Renner T, Rombauts S, Rushton PJ, Sanderfoot A, Schween G, Shiu SH, Stueber K, Theodoulou FL, Tu H, de Peer YV, Verrier PJ, Waters E, Wood A, Yang L, Cove D, Cuming AC, Hasebe M, Lucas S, Mishler BD, Reski R, Grigoriev IV, Quatrano RS, Boore JL: The Physcomitrella genome reveals evolutionary insights into the conquest of land by plants.

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

  45. Merchant SS, Prochnik SE, Vallon O, Harris EH, Karpowicz SJ, Witman GB, Terry A, Salamov A, Fritz-Laylin LK, Maréchal-Drouard L, Marshall WF, Qu LH, Nelson DR, Sanderfoot AA, Spalding MH, Kapitonov VV, Ren Q, Ferris P, Lindquist E, Shapiro H, Lucas SM, Grimwood J, Schmutz J, Cardol P, Cerutti H, Chanfreau G, Chen CL, Cognat V, Croft MT, Dent R, Dutcher S, Fernández E, Fukuzawa H, González-Ballester D, González-Halphen D, Hallmann A, Hanikenne M, Hippler M, Inwood W, Jabbari K, Kalanon M, Kuras R, Lefebvre PA, Lemaire SD, Lobanov AV, Lohr M, Manuell A, Meier I, Mets L, Mittag M, Mittelmeier T, Moroney JV, Moseley J, Napoli C, Nedelcu AM, Niyogi K, Novoselov SV, Paulsen IT, Pazour G, Purton S, Ral JP, no Pachón DMR, Riekhof W, Rymarquis L, Schroda M, Stern D, Umen J, Willows R, Wilson N, Zimmer SL, Allmer J, Balk J, Bisova K, Chen CJ, Elias M, Gendler K, Hauser C, Lamb MR, Ledford H, Long JC, Minagawa J, Page MD, Pan J, Pootakham W, Roje S, Rose A, Stahlberg E, Terauchi AM, Yang P, Ball S, Bowler C, Dieckmann CL, Gladyshev VN, Green P, Jorgensen R, May eld S, Mueller-Roeber B, Rajamani S, Sayre RT, Brokstein P, Dubchak I, Goodstein D, Hornick L, Huang YW, Jhaveri J, Luo Y, Mart nez D, Ngau WCA, Otillar B, Poliakov A, Porter A, Szajkowski L, Werner G, Zhou K, Grigoriev IV, Rokhsar DS, Grossman AR: The Chlamydomonas genome reveals the evolution of key animal and plant functions.

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

  46. Zhang J, Nielsen R, Yang Z: Evaluation of an improved branch-site likelihood method for detecting positive selection at the molecular level.

    Mol Biol Evol 2005, 22(12):2472-2479. PubMed Abstract | Publisher Full Text OpenURL

  47. Gu J, Hilser VJ: Sequence-based analysis of protein energy landscapes reveals nonuniform thermal adaptation within the proteome.

    Mol Biol Evol 2009, 26(10):2217-2227. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  48. Larkum AWD, Orth RJ, Duarte CM: Seagrasses: Biology, Ecology and Conservation. Springer, Dordrecht, The Netherlands; 2006. OpenURL

  49. Touchette BW: Seagrass-salinity interactions: physiological mechanisms used by submersed marine angiosperms for a life at sea.

    J Exp Mar Bio Ecol 2007, 350:194-215. Publisher Full Text OpenURL

  50. Shabala S, Cuin TA: Potassium transport and plant salt tolerance.

    Physiol Plant 2007, 133:651-669. Publisher Full Text OpenURL

  51. Sanchez DH, Siahpoosh MR, Roessner U, Udvardi M, Kopka J: Plant metabolomics reveals conserved and divergent metabolic responses to salinity.

    Physiol Plant 2008, 132(2):209-219. PubMed Abstract | Publisher Full Text OpenURL

  52. Hasegawa PM, Bressan RA, Zhu JK, Bohnert HJ: Plant cellulcar and molecular responses to high salinity.

    Annu Rev Plant Physiol Plant Mol Biol 2000, 51:463-499. PubMed Abstract | Publisher Full Text OpenURL

  53. Sanders D: The salty tale of Arabidopsis.

    Curr Biol 2000, 10(13):R486-R488. PubMed Abstract | Publisher Full Text OpenURL

  54. Pérez M, Invers O, Ruiz JM, Frederiksen MS, Holmer M: Physiological responses of the seagrass Posidonia oceanica to elevated organic matter content in sediments: An experimental assessment.

    J Exp Mar Bio Ecol 2007, 344(2):149-160. OpenURL

  55. Barakat A, Szick-Miranda K, Chang IF, Guyot R, Blanc G, Cooke R, Delseny M, Bailey-Serres J: The organization of cytoplasmic ribosomal protein genes in the Arabidopsis genome.

    Plant Physiol 2001, 127(2):398-415. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  56. Giavalisco P, Wilson D, Kreitler T, Lehrach H, Klose J, Gobom J, Fucini P: High heterogeneity within the ribosomal proteins of the Arabidopsis thaliana 80 S ribosome.

    Plant Mol Biol 2005, 57(4):577-591. PubMed Abstract | Publisher Full Text OpenURL

  57. Dabbs ER: Mutational alterations in 50 proteins of the Escherichia coli ribosome.

    Mol Gen Genet 1978, 165:73-78. PubMed Abstract | Publisher Full Text OpenURL

  58. Andrés O, Kellermann T, López-Giráldez F, Rozas J, Domingo-Roura X, Bosch M: RPS4Y gene family evolution in primates.

    BMC Evol Biol 2008, 8:142. PubMed Abstract | PubMed Central Full Text OpenURL

  59. Hughes AR, Bando KJ, Rodriguez LF, Williams SL: Relative effects of grazers and nutrients on seagrasses: a meta-analysis approach.

    Marine Ecology Progress Series 2004, 282:87-99. Publisher Full Text OpenURL

  60. Beer S, Shomer-Ilan A, Waisel Y: Carbon metabolism in seagrasses: Patterns of photosynthetic CO2 incorporation.

    J Exp Bot 1980, 31:1019-1026. Publisher Full Text OpenURL

  61. Waghmode A, Joshi G: Significance of phosphoglycollate phosphatase and 3-phosphoglycerate phosphatase in photosynthetic carbon assimilation in some marine plants.

    Photosynthetica 1983, 17:193-197. OpenURL

  62. Frost-Christensen H, Sand-Jensen K: The quantum efficiency of photosynthesis in macroalgae and submerged angiosperms.

    Oecologia 1992, 91:377-384. Publisher Full Text OpenURL

  63. Reiskind J, Ginkel TV, Browes G: Evidence that inducible C4-type photosynthesis is a chloroplastic CO2-concentrating mechanism in Hydrilla, a submersed monocot.

    Plant Cell Environ 1997, 20:211-220. Publisher Full Text OpenURL

  64. Mazzuca S, Spadafora A, Filadoro D, Vannini C, Marsoni M, Cozza R, Bracale M, Pangaro T, Innocenti AM: Seagrass light acclimation: 2-DE protein analysis in Posidonia leaves grown in chronic low light conditions.

    J Exp Mar Bio Ecol 2009, 374:13-122. Publisher Full Text OpenURL

  65. Altschul SF, Madden TL, Schäffer AA, Zhang J, Zhang Z, Miller W, Lipman DJ: Gapped BLAST and PSI-BLAST: a new generation of protein database search programs.

    Nucleic Acids Res 1997, 25(17):3389-3402. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  66. Wissler L, Dattolo E, Moore AD, Reusch TBH, Olsen JL, Migliaccio M, Bornberg-Bauer E, Procaccini G: Dr. Zompo: an online data repository for Zostera marina and Posidonia oceanica ESTs.

    Database 2009, 2009:bap009. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

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

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

  68. Notredame C, Higgins DG, Heringa J: T-Coffee: A novel method for fast and accurate multiple sequence alignment.

    J Mol Biol 2000, 302:205-217. PubMed Abstract | Publisher Full Text OpenURL

  69. Edgar RC: MUSCLE: multiple sequence alignment with high accuracy and high throughput.

    Nucleic Acids Res 2004, 32(5):1792-1797. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  70. Suyama M, Torrents D, Bork P: PAL2NAL: robust conversion of protein sequence alignments into the corresponding codon alignments.

    Nucleic Acids Res 2006, (34 Web Server):W609-W612. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  71. Yang Z: PAML 4: phylogenetic analysis by maximum likelihood.

    Mol Biol Evol 2007, 24(8):1586-1591. PubMed Abstract | Publisher Full Text OpenURL

  72. Goldman N, Yang Z: A codon-based model of nucleotide substitution for protein-coding DNA sequences.

    Mol Biol Evol 1994, 11(5):725-736. PubMed Abstract | Publisher Full Text OpenURL

  73. Ashburner M, Ball CA, Blake JA, Botstein D, Butler H, Cherry JM, Davis AP, Dolinski K, Dwight SS, Eppig JT, Harris MA, Hill DP, Issel-Tarver L, Kasarskis A, Lewis S, Matese JC, Richardson JE, Ringwald M, Rubin GM, Sherlock G: Gene ontology: tool for the unification of biology. The Gene Ontology Consortium.

    Nat Genet 2000, 25:25-29. PubMed Abstract | Publisher Full Text OpenURL

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

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

  75. Kanehisa M, Araki M, Goto S, Hattori M, Hirakawa M, Itoh M, Katayama T, Kawashima S, Okuda S, Tokimatsu T, Yamanishi Y: KEGG for linking genomes to life and the environment.

    Nucleic Acids Res 2008, (36 Database):D480-D484. PubMed Abstract | PubMed Central Full Text OpenURL

  76. Bremer K: Early Cretaceous lineages of monocot flowering plants.

    Proc Natl Acad Sci USA 2000, 97(9):4707-4711. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  77. Chalupska D, Lee HY, Faris JD, Evrard A, Chalhoub B, Haselkorn R, Gornicki P: Acc homoeoloci and the evolution of wheat genomes.

    Proc Natl Acad Sci USA 2008, 105(28):9691-9696. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  78. Chaw SM, Chang CC, Chen HL, Li WH: Dating the monocot-dicot divergence and the origin of core eudicots using whole chloroplast genomes.

    J Mol Evol 2004, 58(4):424-441. PubMed Abstract | Publisher Full Text OpenURL

  79. Gaut BS, Doebley JF: DNA sequence evidence for the segmental allotetraploid origin of maize.

    Proc Natl Acad Sci USA 1997, 94(13):6809-6814. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  80. Paterson AH, Bowers JE, Chapman BA: Ancient polyploidization predating divergence of the cereals, and its consequences for comparative genomics.

    Proc Natl Acad Sci USA 2004, 101(26):9903-9908. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  81. Schneider TD, Stephens RM: Sequence logos: a new way to display consensus sequences.

    Nucleic Acids Res 1990, 18(20):6097-6100. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL