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

  • Lothar Wissler1,

    Affiliated with

    • Francisco M Codoñer2,

      Affiliated with

      • Jenny Gu1,

        Affiliated with

        • Thorsten BH Reusch3,

          Affiliated with

          • Jeanine L Olsen4,

            Affiliated with

            • Gabriele Procaccini5Email author and

              Affiliated with

              • Erich Bornberg-Bauer1Email author

                Affiliated with

                BMC Evolutionary Biology201111:8

                DOI: 10.1186/1471-2148-11-8

                Received: 2 June 2010

                Accepted: 12 January 2011

                Published: 12 January 2011

                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 [1315]. 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 [3335] 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).
                http://static-content.springer.com/image/art%3A10.1186%2F1471-2148-11-8/MediaObjects/12862_2010_1617_Fig1_HTML.jpg
                Figure 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, 7680] 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

                Branch

                Cluster ID

                Arabidopsis gene description

                p -value

                LCA

                orthomcl1184

                60 S ribosomal protein L14 (RPL14A)

                <0.001

                LCA

                orthomcl3768

                proteasome maturation factor UMP1 family protein

                0.001

                LCA

                orthomcl1461

                annexin 7, calcium-dependent phospholipid binding (ANNAT7)

                0.001

                LCA

                orthomcl1674

                cytochrome c oxidase 6B (COX6B)

                0.002

                LCA

                orthomcl538

                chaperonin 20, calmodulin binding (CPN20)

                0.003

                LCA

                orthomcl4618

                PSAE-1

                0.004

                LCA

                orthomcl5414

                light harvesting complex PSII 5 (LHCB5)

                0.004

                LCA

                orthomcl1707

                ubiquinol-cytochrome C reductase subunit, mitochondrial

                0.005

                LCA

                orthomcl3901

                nuclear encoded CLP protease 5 (CLPP5)

                0.006

                LCA

                orthomcl1048

                ferredoxin 3 (ATFD3)

                0.007

                LCA

                orthomcl4111

                calmodulin binding (PSAN)

                0.007

                LCA

                orthomcl1171

                protease inhibitor/seed storage/lipid transfer protein (LTP)

                0.007

                LCA

                orthomcl1074

                60 S ribosomal protein L37 (RPL37A)

                0.007

                LCA

                orthomcl433

                lipid transfer protein 3, lipid binding (LTP3)

                0.008

                LCA

                orthomcl1625

                fructose-bisphosphate aldolase, putative

                0.008

                LCA

                orthomcl3789

                PHD finger protein-related

                0.009

                LCA

                orthomcl801

                60 S ribosomal protein L9 (RPL90B)

                0.013

                LCA

                orthomcl1693

                UDP-D-apiose/UDP-D-xylose synthase 1 (AXS1)

                0.016

                LCA

                orthomcl1038

                60 S ribosomal protein L18 (RPL18C)

                0.020

                LCA

                orthomcl922

                glutaredoxin 4, metal ion binding (GRX4)

                0.021

                LCA

                orthomcl1070

                60 S ribosomal protein L6 (RPL6B)

                0.022

                LCA

                orthomcl2948

                scorbate peroxidase 4 (APX4)

                0.025

                LCA

                orthomcl1822

                FK506 binding/peptidyl-prolyl cis-trans isomerase (FKBP15-2)

                0.025

                LCA

                orthomcl626

                40 S ribosomal protein S3A (RPS3aB)

                0.025

                LCA

                orthomcl3845

                copper ion bindng/electron carrier (DRT112)

                0.026

                LCA

                orthomcl1565

                40 S ribosomal protein S15 (RPS15C)

                0.028

                LCA

                orthomcl4326

                PQ-loop repeat family protein/transmembrane family protein

                0.041

                Po

                orthomcl469

                copper ion binding/glutamate-ammninoa ligase (ATGSR1)

                <0.001

                Po

                orthomcl1126

                cytidylate kinase/uridylate kinase (PYR6)

                0.001

                Po

                orthomcl4752

                glycine dehydrogenase, decarboxylating (GDCH)

                0.002

                Po

                orthomcl1625

                fructose-bisphosphate aldolase, putative

                0.003

                Po

                orthomcl1125

                40 S ribosomal protein S24 (RPS24B)

                0.005

                Po

                orthomcl1070

                60 S ribosomal protein L6 (RPL6B)

                0.007

                Po

                orthomcl1673

                cytochrome c-2 (CYTC-2)

                0.008

                Po

                orthomcl1473

                C2 domain-containing protein

                0.011

                Po

                orthomcl1450

                fructose-bisphosphate aldolase, putative

                0.014

                Po

                orthomcl824

                mitochondrial ATP synthase g subunit family protein

                0.014

                Po

                orthomcl4197

                enhancer of sos3-1, metal ion binding/protein binding (ENH1)

                0.018

                Po

                orthomcl4326

                PQ-loop repeat family protein/transmembrane family protein

                0.022

                Po

                orthomcl1896

                microsomal glutathione s-transferase, putative

                0.028

                Po

                orthomcl5121

                frostbite 1, NADH dehydrogenase, ubiquinone (FRO1)

                0.028

                Po

                orthomcl2960

                unknown protein

                0.029

                Po

                orthomcl1930

                cornichon family protein

                0.041

                Po

                orthomcl433

                lipid transfer protein 3, lipid binding (LTP3)

                0.043

                Po

                orthomcl1635

                histone H1-3 (HIS1-3)

                0.046

                Zm

                orthomcl2446

                photosystem I subunit L (PSAL)

                0.001

                Zm

                orthomcl1812

                PS II subunit O-2, oxygen-evolving/poly(U) binding (PSBO2)

                0.003

                Zm

                orthomcl538

                chaperonin 20, calmodulin binding (CPN20)

                0.003

                Zm

                orthomcl3901

                nuclear encoded CLP protease 5 (CLPP5)

                0.004

                Zm

                orthomcl433

                lipid transfer protein 3, lipid binding (LTP3)

                0.005

                Zm

                orthomcl414

                structural constituent of ribosome

                0.005

                Zm

                orthomcl4111

                calmodulin binding (PSAN)

                0.007

                Zm

                orthomcl591

                RuBisCO activator (RCA)

                0.008

                Zm

                orthomcl1057

                photosystem II subunit R (PSBR)

                0.009

                Zm

                orthomcl953

                dormancy-associated protein-like 1 (DYL1)

                0.010

                Zm

                orthomcl3260

                malate dehydrogenase, cytosolic, putative

                0.012

                Zm

                orthomcl1450

                fructose-bisphosphate aldolase, putative

                0.014

                Zm

                orthomcl5948

                prefoldin 6, unfolded protein binding (PDF6)

                0.015

                Zm

                orthomcl1565

                40 S ribosomal protein S15 (RPS15C)

                0.017

                Zm

                orthomcl1808

                universal stress protein (USP) family protein

                0.031

                Zm

                orthomcl824

                mitochondrial ATP synthase g subunit family protein

                0.032

                Zm

                orthomcl4705

                chlorophyll binding (LHCA3)

                0.043

                Zm

                orthomcl3845

                copper ion bindng/electron carrier (DRT112)

                0.044

                Zm

                orthomcl4618

                PSAE-1

                0.045

                Zm

                orthomcl3789

                PHD finger protein-related

                0.049

                Orthologous gene clusters with evidence for positive selection in at least one of the tested branches leading to Zostera marina (Zm), Posidonia oceanica (Po), and their last common ancestor (LCA; see Figure 1A). Each cluster was annotated using the TAIR9 functional description of the representative A. thaliana ortholog. P -values represent the significance of positive selection inferred by the branch-site test for positive selection.

                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

                Map.ID

                Map.Title

                total

                Po

                Zm

                LCA

                01100

                Metabolic pathways

                16

                6

                8

                9

                03010

                Ribosome

                8

                2

                1

                7

                00195

                Photosynthesis

                7

                0

                6

                4

                00190

                Oxidative phosphorylation

                4

                2

                1

                2

                00710

                Carbon fixation in photosynthetic organisms

                3

                2

                2

                1

                01061

                Biosynthesis of phenylpropanoids

                3

                2

                2

                1

                01062

                Biosynthesis of terpenoids and steroids

                3

                2

                2

                1

                01063

                Biosynthesis of alkaloids derived from shikimate pathway

                3

                2

                2

                1

                01064

                Biosynthesis of alkaloids derived from ornithine, lysine and nicotinic acid

                3

                2

                2

                1

                01065

                Biosynthesis of alkaloids derived from histidine and purine

                3

                2

                2

                1

                01066

                Biosynthesis of alkaloids derived from terpenoid and polyketide

                3

                2

                2

                1

                01070

                Biosynthesis of plant hormones

                3

                2

                2

                1

                00010

                Glycolysis/Gluconeogenesis

                2

                2

                1

                1

                00030

                Pentose phosphate pathway

                2

                2

                1

                1

                00051

                Fructose and mannose metabolism

                2

                2

                1

                1

                00196

                Photosynthesis - antenna proteins

                2

                0

                1

                1

                00480

                Glutathione metabolism

                2

                1

                0

                1

                00020

                Citrate cycle (TCA cycle)

                1

                0

                1

                0

                00053

                Ascorbate and aldarate metabolism

                1

                0

                0

                1

                00240

                Pyrimidine metabolism

                1

                1

                0

                0

                00250

                Alanine, aspartate and glutamate metabolism

                1

                1

                0

                0

                00330

                Arginine and proline metabolism

                1

                1

                0

                0

                00520

                Amino sugar and nucleotide sugar metabolism

                1

                0

                0

                1

                00620

                Pyruvate metabolism

                1

                0

                1

                0

                00630

                Glyoxylate and dicarboxylate metabolism

                1

                0

                1

                0

                00910

                Nitrogen metabolism

                1

                1

                0

                0

                00980

                Metabolism of xenobiotics by cytochrome P450

                1

                1

                0

                0

                03050

                Proteasome

                1

                0

                0

                1

                For each pathway, described by the map ID and the title, the total number of associated PS genes are shown as well as the number of PSGs in each of the three branches Zm, Po, and LCA (see Figure 1A). Note that a gene can be associated to more than one pathway.

                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 ( http://static-content.springer.com/image/art%3A10.1186%2F1471-2148-11-8/MediaObjects/12862_2010_1617_IEq1_HTML.gif ) as a major source of inorganic carbon for photosynthesis [11, 12]. The ability to utilize http://static-content.springer.com/image/art%3A10.1186%2F1471-2148-11-8/MediaObjects/12862_2010_1617_IEq1_HTML.gif 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 [6063]. 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/​), 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 < 1 e -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).

                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/​) 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.

                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/​. 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/​, and mapped to the ortholog clusters via the Arabidopsis gene id.

                Declarations

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

                Authors’ Affiliations

                (1)
                Evolutionary Bioinformatics, Institute for Evolution and Biodiversity, University of Muenster
                (2)
                IRSI-Caixa Fundation, Hospital Universitari Germans Trias I Pujol
                (3)
                Leibniz-Institut fuer Meereswissenschaften IFM-Geomar
                (4)
                Dept. Marine Benthic Ecology and Evolution, Centre for Ecological and Evolutionary Studies, University of Groningen, Centre for Life Sciences
                (5)
                Stazione Zoologica A Dohrn, Villa Comunale

                References

                1. Lambers H, Chapin FS, Pons TL: Plant Physiological Ecology. Berlin: Springer; 1998.
                2. den Hartog C: The seagrasses of the world. Amsterdam: North-Holland Publishing Company; 1970.
                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.
                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.View Article
                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.View Article
                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.
                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.View Article
                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.View Article
                9. Touchette, Burkholder : Overview of the physiological ecology of carbon metabolism in seagrasses. J Exp Mar Bio Ecol 2000, 250 (1–2) : 169–205.PubMedView Article
                10. Cox PA: Water-pollinated plants. Sci Am 1993, 269 (4) : 50–56.View Article
                11. Beer S, Rehnberg J: The acquisition of inorganic carbon by the seagrass Zostera marina. Aquatic Botany 1997, 56: 277–283.View Article
                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.View Article
                13. Barbour MG: Is any angiosperm an obligate halophyte? Am Midl Nat 1970, 84: 105–120.View Article
                14. McConnaughey BH: Introduction to Marine Biology. St. Louis, Missouri: Mosby Company; 1978.
                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.View Article
                16. Marschner H: Mineral Nutrition of Higher Plants. London: Academic Press; 1995.
                17. Sudhir P, Murthy SDS: Effects of salt stress on basic processes of photosynthesis. Photosynthetica 2004, 42: 481–486.View Article
                18. Philbrick CT, Les DH: Evolution of aquatic angiosperm reproductive systems. BioScience 1996, 46: 813–826.View Article
                19. Kuo J, den Hartog C: Seagrasses: Biology, Ecology and Conservation Seagrass Morphology, Anatomy, and Ultrastructure, Dordrecht. Netherlands: Springer; 2006.
                20. Doolittle RF: Convergent evolution: the need to be explicit. Trends Biochem Sci 1994, 19: 15–18.PubMedView Article
                21. Wood TE, Burke JM, Rieseberg LH: Parallel genotypic adaptation: when evolution repeats itself. Genetica 2005, 123 (1–2) : 157–170.PubMedView Article
                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.PubMedView Article
                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.PubMedView Article
                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.PubMedView Article
                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.PubMedView Article
                26. Zhang J: Parallel adaptive origins of digestive RNases in Asian and African leaf monkeys. Nat Genet 2006, 38 (7) : 819–823.PubMedView Article
                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.PubMedView Article
                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.PubMedView Article
                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.View Article
                30. Hurst LD: The Ka/Ks ratio: diagnosing the form of sequence evolution. Trends Genet 2002, 18 (9) : 486.PubMedView Article
                31. Yang, Bielawski : Statistical methods for detecting molecular adaptation. Trends Ecol Evol 2000, 15 (12) : 496–503.PubMedView Article
                32. Hedges SB, Dudley J, Kumar S: TimeTree: a public knowledge-base of divergence times among organisms. Bioinformatics 2006, 22 (23) : 2971–2972.PubMedView Article
                33. Bremer K: Early Cretaceous lineages of monocot flowering plants. Proc Natl Acad Sci USA 2000, 97 (9) : 4707–4711.PubMedView Article
                34. Janssen T, Bremer K: The age of major monocot groups inferred from 800+ rbcL sequences. Bot J Linn Soc 2004, 146: 385–398.View Article
                35. Anderson C, Janssen T: The Timetree of Life. Monocots, Oxford University Press; 2009.
                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.PubMedView Article
                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.PubMedView Article
                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.PubMedView Article
                39. International Brachypodium Initiative: Genome sequencing and analysis of the model grass Brachypodium distachyon. Nature 2010, 463 (7282) : 763–768.View Article
                40. Arabidopsis Genome Initiative: Analysis of the genome sequence of the flowering plant Arabidopsis thaliana. Nature 2000, 408 (6814) : 796–815.View Article
                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.PubMedView Article
                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.PubMedView Article
                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.PubMedView Article
                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.PubMedView Article
                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.PubMedView Article
                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.PubMedView Article
                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.PubMedView Article
                48. Larkum AWD, Orth RJ, Duarte CM: Seagrasses: Biology, Ecology and Conservation. Springer, Dordrecht, The Netherlands; 2006.
                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.View Article
                50. Shabala S, Cuin TA: Potassium transport and plant salt tolerance. Physiol Plant 2007, 133: 651–669.View Article
                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
                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.PubMedView Article
                53. Sanders D: The salty tale of Arabidopsis. Curr Biol 2000, 10 (13) : R486-R488.PubMedView Article
                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.View Article
                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.PubMedView Article
                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.PubMedView Article
                57. Dabbs ER: Mutational alterations in 50 proteins of the Escherichia coli ribosome. Mol Gen Genet 1978, 165: 73–78.PubMedView Article
                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.PubMedView Article
                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.View Article
                60. Beer S, Shomer-Ilan A, Waisel Y: Carbon metabolism in seagrasses: Patterns of photosynthetic CO2 incorporation. J Exp Bot 1980, 31: 1019–1026.View Article
                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.
                62. Frost-Christensen H, Sand-Jensen K: The quantum efficiency of photosynthesis in macroalgae and submerged angiosperms. Oecologia 1992, 91: 377–384.View Article
                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.View Article
                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.View Article
                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.PubMedView Article
                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.PubMedView Article
                67. Li L, Stoeckert CJ, Roos DS: OrthoMCL: identification of ortholog groups for eukaryotic genomes. Genome Res 2003, 13 (9) : 2178–2189.PubMedView Article
                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.PubMedView Article
                69. Edgar RC: MUSCLE: multiple sequence alignment with high accuracy and high throughput. Nucleic Acids Res 2004, 32 (5) : 1792–1797.PubMedView Article
                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.
                71. Yang Z: PAML 4: phylogenetic analysis by maximum likelihood. Mol Biol Evol 2007, 24 (8) : 1586–1591.PubMedView Article
                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
                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.PubMedView Article
                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.PubMedView Article
                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.
                76. Bremer K: Early Cretaceous lineages of monocot flowering plants. Proc Natl Acad Sci USA 2000, 97 (9) : 4707–4711.PubMedView Article
                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.PubMedView Article
                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.PubMedView Article
                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.PubMedView Article
                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.PubMedView Article
                81. Schneider TD, Stephens RM: Sequence logos: a new way to display consensus sequences. Nucleic Acids Res 1990, 18 (20) : 6097–6100.PubMedView Article

                Copyright

                © Wissler et al; licensee BioMed Central Ltd. 2011

                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.

                Advertisement