Skip to main content

Advertisement

Springer Nature is making SARS-CoV-2 and COVID-19 research free. View research | View latest news | Sign up for updates

De novo sequencing and analysis of the Ulva linza transcriptome to discover putative mechanisms associated with its successful colonization of coastal ecosystems

BMC Genomics volume 13, Article number: 565 (2012) Cite this article

Abstract

Background

The green algal genus Ulva Linnaeus (Ulvaceae, Ulvales, Chlorophyta) is well known for its wide distribution in marine, freshwater, and brackish environments throughout the world. The Ulva species are also highly tolerant of variations in salinity, temperature, and irradiance and are the main cause of green tides, which can have deleterious ecological effects. However, limited genomic information is currently available in this non-model and ecologically important species. Ulva linza is a species that inhabits bedrock in the mid to low intertidal zone, and it is a major contributor to biofouling. Here, we presented the global characterization of the U. linza transcriptome using the Roche GS FLX Titanium platform, with the aim of uncovering the genomic mechanisms underlying rapid and successful colonization of the coastal ecosystems.

Results

De novo assembly of 382,884 reads generated 13,426 contigs with an average length of 1,000 bases. Contiguous sequences were further assembled into 10,784 isotigs with an average length of 1,515 bases. A total of 304,101 reads were nominally identified by BLAST; 4,368 isotigs were functionally annotated with 13,550 GO terms, and 2,404 isotigs having enzyme commission (EC) numbers were assigned to 262 KEGG pathways. When compared with four other full sequenced green algae, 3,457 unique isotigs were found in U. linza and 18 conserved in land plants. In addition, a specific photoprotective mechanism based on both Lhc SR and Psb S proteins and a C4-like carbon-concentrating mechanism were found, which may help U. linza survive stress conditions. At least 19 transporters for essential inorganic nutrients (i.e., nitrogen, phosphorus, and sulphur) were responsible for its ability to take up inorganic nutrients, and at least 25 eukaryotic cytochrome P450s, which is a higher number than that found in other algae, may be related to their strong allelopathy. Multi-origination of the stress related proteins, such as glutamate dehydrogenase, superoxide dismutases, ascorbate peroxidase, catalase and heat-shock proteins, may also contribute to colonization of U. linza under stress conditions.

Conclusions

The transcriptome of U. linza uncovers some potential genomic mechanisms that might explain its ability to rapidly and successfully colonize coastal ecosystems, including the land-specific genes; special photoprotective mechanism based on both Lhc SR and Psb S; development of C4-like carbon-concentrating mechanisms; muti-origin transporters for essential inorganic nutrients; multiple and complex P450s; and glutamate dehydrogenase, superoxide dismutases, ascorbate peroxidase, catalase, and heat-shock proteins that are related to stress resistance.

Background

Green algae originated as much as 1500 million years ago, evolving shortly after the endosymbiotic event that gave rise to early photosynthetic eukaryotes [1]. During the evolutionary history of Earth, green algae have become major players in global energy/biomass production and biogeochemical recycling [2]. The green algal genus Ulva Linnaeus (Ulvaceae, Ulvales, Chlorophyta), which includes the genus formerly known as Enteromorpha[3], is a cosmopolitan intertidal macroalga and includes more than 100 species [4]. It is well known for its distribution in all aquatic habitats including freshwater, brackish environments, marine and fully saline environments throughout the warm temperate and tropical regions of the world [5–7]. This group also possesses features that provide them with a competitive edge for rapid and successful colonization in eutrophic conditions, namely a copious production of reproductive spore bodies; the ability to rapidly take up and store high amounts of inorganic nitrogen; and a wide tolerance to adverse environmental conditions such as temperature, light intensity, salinity, and anoxia [8, 9]. In addition, when these algae grow under eutrophic conditions, they have the competitive advantage of rapid and successful colonization, resulting in a “green tide” [8, 10, 11], which has the following deleterious ecological effects: the uncoupling of biogeochemical cycles in sediments from those in the water column [9], a negative effect on seagrass beds due to shading, disruption of feeding by wading birds [12], the development of a lethal environment due to oxygen deficiency [13] and a shift from a high diversity mixture to low-diversity assemblages of fast growing annuals [14]. Because of their ecological effects, green tides have drawn considerable attention from scientists and governments [8, 9, 12].

Whole genome or transcriptome sequencing could elucidate the genetic basis for this physiological characteristic or exploit this group’s unique genetic resources based on global gene expression [15–17]. Over the last few years, the application of newly developed high-throughput sequencing technologies has allowed hundreds of thousands of high-quality sequence reads to be produced de novo from whole genome or transcriptome templates, thus enabling immediate inroads to genetic studies of organisms for which little or no sequence data exist [18]. For non-model organisms, sequencing and assembling of genomes remains challenging and costly, even when considering high throughput sequencing technology [19]. As an alternative, sequencing transcriptomes is less complex and provides fast and cost-effective access to the gene expression profile of an organism [16, 20]. Transcriptomes present a valuable resource for accelerating gene discovery by expanding gene families [21, 22], improving genome annotation [23], elucidating phylogenetic relationships [24], facilitating breeding programs for both plants and animals by providing SSR and SNP markers [25, 26], and allowing large-scale expression analysis [27, 28] and rapid identification of transcripts involved in specific biological processes [29]. Recent studies have demonstrated the success of 454 de novo sequencing, assembly and analysis of transcriptomes in non-model organisms with no prior genomic resources [30–33]. The generation of such large-scale sequence data will enable functional analyses that were previously limited to model organisms and their rapid application in ecologically important taxa [34].

To date, genomic information about algae is still limited compared to land plants. Moreover, algae are generally far more diverse and evolutionarily divergent than land plants [35]. Currently, nearly complete genomic information is available publicly for only five green microalgae species: Chlamydomonas reinhardtti (Chlorophyceae) [36], Chlorella variabilis (Trebouxiophyceae) [37], Ostreococcus tauri (Prasinophyceae) [38], Micromonas (Prasinophyceae) [39] and Volvox carteri (Chlorophyceae) [40]. To date, there is little genomic data available for the Ulvophyceae, a group that is diverse and evolutionarily divergent from all of the currently sequenced genomes of green microalgae. For the Ulvophyceae group, EST resources have been developed using Sanger sequencing in three species, Acetabularia acetabulum[35], Ulva linza[41] and Ulva prolifera[42]. A comprehensive description of the full complement of genes as a transcriptome, however, is still unavailable. Ulva linza is one of the species found in the mid to low intertidal on bedrock, and it is also a major contributor to biofouling [41]. It must cope with more dramatic changes than land plants because of tidal changes in light intensity, temperature, salinity, and wave action and with the biotic stresses characteristic of dense coastal ecosystems. In this study, we sequenced the transcriptome of U. linza under a wide variety of environmental stress conditions with the aim of uncovering some of potential genomic mechanisms that might explain its ability to rapidly and successfully colonize the coastal ecosystems and form green tides.

Results and Discussion

454 sequencing and de novo assembly of the transcriptome

The normalized cDNA library of U. linza cells grown under the normal and six stress conditions was constructed using SMART technology. The library was subjected to a one-half plate run with the 454 GS FLX Titanium platform, and the transcriptome was assembled from the resulting sequencing reads. Sequencing of the cDNA library generated a total 503,789 raw reads (SRA053388), with an average sequence length of 396 bases (Figure 1, Table 1). After trimming the adapter sequences and removing the short sequences with low complexity and low quality scores, 445,787 high quality reads were obtained, corresponding to 88.5% of the original raw sequences. A total of 382,884 reads were assembled into 13,426 contiguous sequences (contigs), accounting for 85.9% of the assembled reads, while 38,961 reads remained as singletons (reads not assembled into contigs). The size of contigs ranged from 50 to 9,881 bases, with an average length of 1,000 bases. As expected for a randomly fragmented transcriptome, there was a positive relationship between the length of a given contig and the number of reads assembled into that contig. The length distribution and read depth of the contigs are shown in Figure 2, revealing that more than 9,578 contigs are greater than 400 bp, with an average depth of coverage of 14.48 sequences per nucleotide position. Contiguous sequences were further assembled into 10,784 isotigs, which are the putative transcripts constructed using the overlapping contig reads provided as input to the Newbler cDNA assembler. The average length of the isotigs was 1,515 bases, with an N50 of 1,856 bases (50% of the assembled bases were incorporated into isotigs greater than 1,856 bases). The coverage depth for isotigs ranged from 1 to 14, with an average of 1.24 contigs assembled into each isotig. When compared to the transcriptomes of green algae and land plants including C. reinhardtii, V. carteri, P. patens and A. thaliana, we found that a total of 6,519 (82.1%) common genes could be identified in the transcriptome of U. linza.

Figure 1
figure1

Read length distribution of the Ulva linza transcriptome data.

Full size image
Table 1 Summary of the sequencing and assembly
Full size table
Figure 2
figure2

Contig length as a function of the contig read depth.

Full size image

Comparison to other proteomes of green algae

To identify novel features of the U. linza transcriptome, a BLASTX search was used to compare the isotigs to sequences from other organisms belonging to the green algae lineage. In this analysis we included four fully sequenced organisms that diverged at different times during evolution: C. reinhardtii, a model organism for photosynthetic unicellular eukaryotes [36]; O. tauri, which is the smallest free-living eukaryote identified to date and one of the most ancient clade among Chlorophyta with the features of genome compaction and gene family downsizing [38]; C. variabilis, which is used as a model system for studying algal symbioses and algal–viral interactions of freshwater organisms [37]; and Micromonas which offers valuable insights into ecological differentiation and the dynamic nature of early plant evolution [39]. Sequence comparisons revealed that genes shared between U. linza and the four unicellular green microalgae have a consistent frequency and similarity (Figure 3). In contrast, the frequency of unique sequences when compared with different green algal species was clearly different, ranging from about 756 in C. reinhardtii to 2,485 in O. tauri. The 1,729 reads, which have no counterpart in O. tauri but do in C. reinhardtii, were assembled into 1,379 isotigs. Totally 762 isotigs could be annotated, among which at least 441 sequences were also conserved in Volvox, Micromonas and Chlorella (Additional file 1: Table S1). This indicates that O. tauri lost genes compared to other green algae.

Figure 3
figure3

Comparison between the Ulva linza transcriptome and four closely related chlorophyta genomes.

Full size image

In total, 4,081 isotigs in U. linza were found to be conserved in all four of the green algae, thus they might have existed before the species diverged and may therefore play basic roles for their survival (Figure 4). The lineage-specific genes have evolved after the species diverged by gene transfer or gene loss, and may contribute to their specific characteristics and adaptation to their living environments. The analysis also provides us with 3,457 isotigs that are exclusive to U. linza, not found in the other four green algal genomes queried. This corresponds to 32.1% of the total isotig repertoire in U. linza, which is similar to the fraction of unique genes found in Nannochloropis gaditana (30.2%) [43] and C. reinhardtii (35.3%), but is much higher than that of C. variabilis NC64A (18.6%), Micromonas sp. RCC299 (21.0%) and O. tauri (17.5%). Among the 3,457 specific isotigs in U. linza, only 195 could be annotated and 18 of which were conserved in land plants including disease resistance proteins, an ABC transporter, etc. (Additional file 2: Table S2). U. linza lives in intertidal zone of coastal ecosystems, which is a transitional environment from sea water to land, and is therefore exposed to direct sunshine several hours per day. Moreover, the sessile form of life of U. linza cannot escape from stress conditions by swimming deeper, a behavior that is typical of planktonic algae [44]. We predicted that these land-specific genes in U. linza may be related to the colonization of the special coastal ecosystems.

Figure 4
figure4

Venn diagram representation of shared/unique genes of U. linza in comparison with four other sequenced green algae of Chlamydomonas reinhardtii , Chlorella variabilis NC64A, Micromonas sp. RCC299 and Ostreococcus tauri.

Full size image

Functional annotation

All sequences were aligned against the local non-redundant (nr) protein database downloaded from the National Center for Biotechnology Information (NCBI) using the BLASTx algorithm. When the E-value cutoff was set at 10-2, a total of 304,101 reads (72.2% of total) had significant BLAST matches. According to the scores obtained by the BLASTX, the high quality reads were assigned to the taxonomic tree that was representative of Bacteria (10,685; 3.5%) and Eukaryota (218,725; 71.9%). As expected, the taxonomic assignment was well represented by Chlorophyta (90,074; 29.6%) and Embryophyta (11,229; 3.7%) within the Viridiplantae group. Within Chlorophyta, the top four assignments were assigned to Chlorella, Ulva, Volvox and Chlamydomonas.

The isotigs were further annotated with GO terms based on their sequence similarities to known proteins in the UniProt database annotated with GO terms as well as the InterPro and Pfam domains they contain. A total of 4,368 isotigs (40.5%) were functionally annotated with 13,550 GO terms. Among these, 7,192 were assigned at least one GO term in the biological process category, 6,108 in the molecular function category and 3,233 in the cellular component category (Figure 5). In the biological process class (second level GO terms), the majority of the GO terms were grouped into cellular (42.7%) and metabolic (45.2%) processes. However, in the category of molecular function, the vast majority of GO terms were classified into binding (46.7%) and catalytic activities (47.4%). Under the category of cellular components, 32% of all GO terms corresponded to cell parts and organelles.

Figure 5
figure5

Functional categories based on GO annotations.

Full size image

The Kyoto Encyclopedia of Genes and Genomes (KEGG) orthology (KO) is a classification system that provides an alternative functional annotation of genes based on their associated biological pathways (Figure 6). To reconstruct the metabolic pathways involved in U. linza, both reads and isotigs were assigned to the KEGG. In total 56,553 reads and 2,404 isotigs having enzyme commission (EC) numbers were assigned to 262 KEGG pathways. The highest numbers of reads were assigned to pathways on the second level of energy metabolism (25,355), carbohydrate metabolism (17,490), and immune systems (12,391) and on the third level of carbon fixation (9,356), ABC transporters (9,051) and photosynthesis (5,281). In contrast, the highest numbers of isotigs were assigned to pathways on the second level of carbohydrate metabolism (323), translation (311), and amino acid metabolism (254) and on the third level of ribosome (108), spliceosome (91), and purine metabolism (63) (Additional file 3: Table S3). The reason for the difference between reads and isotigs was that only a fraction of reads could be assembled into isotigs, and those unassembled reads may also have functional implications. Read-based KEGG annotation also can be used to estimate the expression abundance of certain KEGG functions. Pathways that possessed the highest number of isotigs were similar to those in C. reinhardtii, V. carteri and P. patens, but different from Arabidopsis thaliana on the third level.

Figure 6
figure6

Metabolic pathway analysis using KEGG.

Full size image

Photoprotection mechanisms in U. linza

In taxa of photosynthetic eukaryotes, the thylakoid membrane-integral light-harvesting complex (LHC) antenna systems, especially LHCII in Photosystem II, play important roles in regulating energy flow to reaction centers [45]. The LHC systems are generally used to harvest and transfer excitation energy into the reaction centers to drive photosynthesis. However, under excess light conditions, they undergo a conformational change and activate a quenching state for energy dissipation, termed non-photochemical quenching (NPQ) [46, 47].

The light-harvesting antenna complex in PSII consists of two types of proteins: major trimer-forming LHCII antenna proteins and minor LHCII proteins. In land plants, major LHCII proteins are encoded by multiple Lhcb1, Lhcb2 and Lhcb3 genes [48], but in green algae, the major LHCII proteins encoded by Lhcbm genes showed great diversity [48–51]. Nineteen Lhcbm-encoding genes have been found in the transcripome of U. linza, and of these, the five genes that had full-length sequences were selected for phylogenetic tree construction. The phylogenetic tree shows that the major LHCII genes have undergone significant divergence in different organisms; this may be an adaptation mechanism to optimize light-harvesting capacity and acclimation to environmental conditions. Three minor PSII antenna have been found in land plants, namely CP26 (Lhcb5), CP29 (Lhcb4), and CP24 (Lhcb6); all of these are located at the interface between the major LHCII proteins and PSII [52, 53]. The minor proteins CP26 and CP29 have been shown to play direct roles in NPQ in land plants [54]. The CP24 protein evolved differently and more recently in land plants as it is absent in all green algae [49, 50]. In U. linza, only the CP26 and CP29 proteins were found as expected. From the phylogenetic tree, it can be seen that CP26 and CP29 are found in all green plant groups examined and evolved distinctly to the major LHCII proteins, which were present prior to diversification of the major lineages (Additional file 4: Figure S1). The conserved nature and distribution of the minor LHCII proteins suggests that they have a significant function.

We also found the co-expressed Lhc SR and Psb S genes, which are photoprotection-related LHC genes. The PsbS protein was recently implicated in NPQ in land plants [55–57]. PsbS was not the site of quenching itself; instead it acts as a pH-dependent trigger to activate the quenching sites by regulating the organization of LHCII antenna system with the PSII core complex [58–61]. However, in green algae, such as Chlamydomonas, no accumulation of the PsbS protein could be detected, although the genes are structurally similar to PsbS [50, 62]. The lack of the PsbS and CP24 proteins in green algae implied that the mechanistic basis of NPQ in these groups is likely different. In fact, the stress-related LHC member LhcSR proteins in green algae have been shown to participate in NPQ, likely indicating an early photoprotective mechanism [63]. Furthermore, the LhcSR3 protein in Chlamydomonas appears to merge both pH-sensing and energy-quenching functions in plants, which are accomplished by PsbS and monomeric Lhcb proteins, respectively [64]. Among the land plants, the mechanism of NPQ was unique in the moss, P. patens, where both the LhcSR and PsbS proteins are active [65]. The mechanism for NPQ was hypothesized to have evolved from being LhcSR-based to being PSBS-based during the transition to land colonization [63, 65]. However, the transcription profiles of both lhc SR and psb S genes in U. linza indicate that the PsbS-based NPQ mechanism emerged before the land colonization. In addition, we found that both the psb S and lhc SR genes were induced by high light (>600 μmol photons m−2 s−1) and the expression patterns of these two genes were consistent (Figure 7), which indicates that both PsbS and LhcSR are active in photoprotection of U. linza. This result also suggests that the photoprotective mechanism in Ulva may differ from that of unicellular green algae. A potential explanation for this difference is that Ulva must cope with dramatic changes in the intensity and spectral quality of light in rocky coastal ecosystems. U. linza is a photophilous alga commonly distributed in shallow and horizontal habitats, like sun-exposed rocky shores. The light irradiance varies more dramatically for the algae than for land plants because of the morning and evening tides.

Figure 7
figure7

Relative mRNA expression of lhc SR and psb S in U. linza exposed to 40-2000 μmol photons m −2 s −1 for 1 h. Data were means ± SD (n=3).

Full size image

Carbon fixation pathway

Besides photoprotection mechanisms, efficient carbon fixation is essential for photosynthesis, and as such likely to facilitate the colonization of coastal ecosystems and the accumulation of high U. linza biomass. Aquatic photosynthetic organisms are exposed to dramatic changes in the supply of dissolved inorganic carbon (Ci; CO2 and/or HCO3-). To acclimate to the Ci-limiting environmental conditions, they have developed a carbon concentrating mechanisms (CCM) to optimize the photosynthesis rate under the CO2-limiting condition [66, 67]. At least one α-carbonic anhydrase (CA) with targeting signals for the chloroplast lumen was found in U. linza. Lumenal-targeted αCAs are well known from C. reinhardtii and higher plants and are essential for growth under ambient CO2 concentrations [68]. Thus, we predicted that the α-CA-based CCM type found in U. linza would show a strong similarity to the CCM in C. reinhardtii[69]. Few genes common to organisms that actively or passively enhance Ci influx in CCMs were found in U. linza. Only a putative soluble protein, LciB, was identified as a candidate for a chloroplast Ci transporter, and no Ci membrane transporter candidate was found. Five low/limited CO2-induced proteins were identified, including two regulatory factors (CCM1 and LCR1) that control the expression of the CO2-related genes.

C4-photosynthesis is an important feature in some higher plants in tropical environments, as it helps them cope well with high light intensities, high temperatures, and dryness [70, 71]. Some intermediate products of the C4 carbon-fixation pathway have been detected in diatoms, brown algae, green algae and red algae [72–77]. The C4-like CCM mechanism in green algae of O. tauri[38] and Micromonas[39] may provide them an ecological advantage in the CO2-limiting conditions of phytoplankton blooms. In Ulva prolifera, coexistence of C3 and C4-like pathway may contribute to its wide distribution and forming blooms [78]. In U. linza, genes encoding all of the enzymes required for C4 photosynthesis were identified implying a putative C4-like CCM mechanism (Additional file 5: Table S4).

Transporters

Transporters are mainly responsible for the acquisition, redistribution and compartmentalization of organic nutrients and inorganic ions, as well as for the efflux of toxic compounds and metabolic end products. A total of 114 transport members were found in the U. linza transcriptome, and the complement resembles that of both green algae and land plants. The ATP-binding cassette (ABC) family (36 members) and the P-type adenosine triphosphatase (ATPase) family (12 members) are large compared to other transporters. One homologue of a plant-specific peroxisomal ABC transporter was found in Ulva, but is absent in the genomes of other green algae. Ulva is known for its ability to take up inorganic nutrients and use them to grow rapidly, which allows it to be successful in nitrogen-rich areas and form the green tides [8, 11]. The transcriptome encodes at least 19 transporters for essential inorganic nutrients (i.e., nitrogen, phosphorus, and sulphur). U. linza obtains nitrogen in at least three different ways, and ammonium is the best nitrogen source for Ulva growth, followed by urea and nitrate [79]. The complete set of genes allowing transport and assimilation of these substrates include six ammonium transporters, three high-affinity nitrate transporters, and two urea transporters. The high number of ammonium transporters indicates that Ulva is a strong competitor for this resource, which is very efficiently pumped from surrounding seawater and is probably the privileged nitrogen source for Ulva. The inorganic nutrient transporters found in Ulva showed multiple origins. The unicellular green algae O. tauri has two prokaryote-like ammonium transporters, and one was found in U. linza (Additional file 6: Figure S2). All of the three phosphate transporters found in U. linza are absent in higher land plants, two occur only in green algae, and the third one is present in moss (Additional file 7: Figure S3). These results suggest that the phosphate transporting mechanism in Ulva may be old and associated with aquatic environments and that it was lost during adaptation to land. Of the three sulfate transporters, two are in the Na+/SO42--family; neither of these is present in higher land plants but one is present in moss. The third one is conserved in both green algae and high plants (Additional file 8: Figure S4).

Cytochrome P450 oxidoreductases

Eukaryotic cytochrome P450s are oxidoreductases, most of which catalyze NADPH- and O2-dependent hydroxylation reactions [80]. Plant P450s participate in a myriad of biochemical pathways, including those devoted to the synthesis of plant products and plant growth regulators [81]. The complexity of the P450 family is reflected in the evolution of the metabolic complexity of organisms. U. linza has at least 25 cytochrome P450s, which is a higher number than that found in other algae (eg. the unicellular green alga Chlamydomonas (13) and brown alga Ectocarpus (12)), but a lower number than that found in land plants (e.g., Physcomitrella (71) and Arabidopsis (286)). Nearly half of the P450s in U. linza are conserved in algae (even cyanobacteria) (Additional file 9: Figure S5). These P450s likely originated in algae and may be related to acclimatization to the aquatic environment. The role of these conserved P450s is poorly understood. A few of members of the CYP97 family are known to be involved in accumulating secondary ketocarotenoids in response to stress in terrestrial plants [82, 83]. Some P450s genes in U. linza are homologous to those in fungus, land plants, and even animals but are absent in unicellular green algae. We predicted that these genes may function in allelopathy in U. linza, thus allowing it to produce toxic substances for defense against herbivory or other phototrophic organisms that compete for light and nutrients [84–86]. U. linza is known to exhibit strong allelopathy: as one of the most common green tide macroalgal species, it has to outcompete other coastal primary producers and establish dominant populations in a short time [87].

Stress-related proteins

Glutamate dehydrogenase (GDH) catalyzes the reversible oxidative deamination of glutamate to α-ketoglutarate and ammonia. In plants, the enzyme can work in either direction depending on environment and stress [88, 89]. It links the two fundamental metabolic processes involving carbon and nitrogen within a cell, and transgenic plants expressing microbial GDHs exhibit improved tolerance to herbicides, water deficits, and pathogen infections [89, 90]. Two GDH genes were found in Arabidopsis, and only one was found in Chlamydomonas. In contrast, eight putative GDH genes were found in U. linza. Seven of these GDH genes have similarities to those found in bacteria, fungus, and even animals (Additional file 10: Figure S6). Lateral gene transfer has been shown to play a significant role in the evolution of the GDH genes [91]. Thus, we predicted that the diverse GDH genes in Ulva may have been obtained recently by gene transfer and function in maintaining a balance between carbon and nitrogen metabolites and enhancing the stress tolerance.

Reactive oxygen species generated during photosynthesis are scavenged by superoxide dismutases (SODs), ascorbate peroxidase (APX), and catalase (CAT) in plants [92–94]. Four putative APX genes were found in U. linza, which is more than has been found in C. reinhardtii (2), V. carteri (3), and P. patens (3) but less than in A. thaliana (6). All four APX genes are conserved in green algae and land plants (Additional file 11: Figure S7). In contrast, the two CAT genes in U. linza show higher similarities to those in fungi and bacteria than to those in green algae, and they have no homologous genes in land plants (Additional file 12: Figure S8). The APX and CAT genes, which are involved in scavenging reactive oxygen, may have evolved in different ways. Two Fe-type and one Mn-type SODs were found in U. linza, but no Cu/Zn-type or Ni-containing SOD was found. However, a Cu/Zn-type SOD was found to be expressed in Ulva fasciata in response to copper stress [95]. We predicted that Fe-type and Mn-type SODs may play dominant roles in response to broad stress conditions, such as high salt, high light, and high temperature conditions. In contrast, the Cu/Zn-type SOD may be related only to heavy metal stresses.

Heat-shock proteins (Hsps) constitute a major class of stress-response proteins that mainly assist in protein refolding under stress conditions. Thirty-two Hsp genes were found in U. linza from four conserved families: Hsp70, Hsp90, Hsp100, and low-molecular-mass Hsps (sHsps). In U. linza, 8 sHsps belong to the Hsp20 family, and these Hsp20 proteins show high similarities to prokaryotes as to green algae and plants (Additional file 13: Table S5; Additional file 14: Figure S9). Similar observations were also made in terrestrial plants, where the sHsps are highly diverse, such as in Arabidopsis, 13 different sHsps are grouped into six classes [96]. In contrast to Hsp20 proteins, most members of the Hsp70, Hsp90 and Hsp100 families found in U. linza were conserved well in eukaryotic algae and plants, with the exception of only three (Additional file 13: Table S5). Compared to other unicellular green algae, number of Hsp90 and Hsp100 proteins was increased in U. linza, especially members of the Hsp90 family (Additional file 15: Table S6). The Hsp90 family is distinct in that most of its known substrates are signal transduction proteins, such as steroid hormone receptors and signaling kinases [97], which may play a key role in signal transduction networks, protein degradation and protein trafficking. Thus, the increasing number of Hsp90 members indicates that U. linza may possess a Hsp/chaperone network that is more complicated than that found in unicellular algae, and it may be related to the acclimation of U. linza to the intricate environment of the intertidal zone.

Conclusions

The green algal genus Ulva Linnaeus (Ulvaceae, Ulvales, Chlorophyta) is well known for its wide distribution in marine, freshwater, and brackish environments throughout the world. The Ulva species are also highly tolerant of variations in salinity, temperature, and irradiance and are the main cause of green tides, which can have deleterious ecological effects. However, only limited genomic information is currently available in this non-model and ecologically important species. The transcriptome of U. linza uncovers some of potential genomic mechanisms that might explain its ability to rapidly and successfully colonize the coastal ecosystems, including the specific photoprotective mechanism based on both Lhc SR and Psb S; development of C4-like carbon-concentrating mechanisms; muti-origin transporters for essential inorganic nutrients; multiple and complex P450s, glutamate dehydrogenase, superoxide dismutases, ascorbate peroxidase, catalase and heat-shock proteins that are related to stress resistance.

Methods

Sampling and culture conditions

The Ulva linza samples were collected in May of 2010 from intertidal zone 20 (35°35´N, 119°30´E) of Zhanqiao Wharf, Qingdao, China. In the laboratory, the intact 21 samples were washed several times with sterile seawater, sterilized with 1% sodium hypochlorite for 2 min, and then rinsed with autoclaved seawater. The sterilized material was then placed into an aquarium containing enriched seawater and which was aerated and maintained at 15°C under a 12:12 h LD photoperiod with 120 μmol photons m−2 s−1 provided by cool white fluorescent tubes.

Library preparation and 454 sequencing

The samples of Ulva linza were subjected to different stress conditions: low temperature (6°C, 2 h), high temperature (42°C, 1 h), high light (400 μmol photons m−2 s−1, 1 h), high salt (93‰, 3 h) and UV-B stress (60 μw cm-2, 3 h). All the treated samples were frozen in liquid nitrogen, and total RNA was extracted using the TRIzol reagent (Invitrogen) according to the manufacturer’s recommendations. The mRNA was purified from the total RNA using the Oligotex mRNA Midi Kit (QIAGEN, Germany). Double-stranded cDNA was then synthesized using the SMART cDNA Library Construction kit (Clontech, USA) by following the manufacturer’s protocol. The cDNA was purified using Qiagen QIAquick PCR purification spin columns. Normalization was performed using the TRIMMER cDNA normalization kit (Evrogen) to decrease the prevalence of abundant transcripts before sequencing, and the quality of cDNA was checked using the Agilent 2100 Bioanalyzer. A half-plate sequencing run was performed at the Tianjin Biochip Corporation following the manufacturer’s protocols. Raw data has been submitted to the Status of the NCBI Sequence Read Archive (http://www.ncbi.nlm.nih.gov/Traces/sra_sub/sub.cgi?).

De novo assembly

The “-cdna” option of the latest release of Roche 454’s Newbler (version 2.5p1) was used to assemble transcriptome. Newbler creates isotigs out of contigs that are consistently connected by a subset of reads. Each isotig corresponds to an alternative transcript, and any contigs or isotigs that share any read overlaps are put into the same isogroup.

Functional annotation

We used local BLASTX to align 454 reads to the GenBank nonredundant protein (nr) database using an E value threshold of 10-2. The unigene sequences were also translated into proteins using ESTScan and the translated protein sequences were then compared to InterPro and pfam domain databases. The gene ontology (GO) terms were assigned to each unigene based on the GO terms annotated to its corresponding homologues in the UniProt database, and based on the InterPro and Pfam domains using interpro2go and pfam2go mapping files provided by the GO website, respectively [98]. Then, the MEGAN v4.0 software [99] was used to assign the annotated reads to KEGG pathways [100]. The number of reads assigned to each KEGG functional category can be used to evaluate gene quantities and expression levels.

Real-time RT-PCR analysis

Real-time RT-PCR approach was employed to analyze the expression of lhc SR and psb S genes under different stress conditions. The primers for real-time RT-PCR analysis in this experiment are shown in Additional file 16: Table S7. The different light intensity treatments were provided at 120, 240, 400, 500, and 600 μmol photons m−2 s−1 and lasted 1 h at 15°C. For each condition, three biological replicates were conducted. The expression of the genes under different conditions was measured by real-time PCR. The qPCR was performed with the ABI StepOne plus Real-Time PCR System (Applied Biosystems, USA) and performed using SYBR Green fluorescence (Takara) according to the manufacturer’s instructions. The RNA levels were expressed relative to those of the 18S rDNA sequence. The thermal profile for real-time PCR was 30 s followed by 40 cycles at 95°C for 5 s, 55°C for 10 s and 72°C for 30 s. Dissociation curve analysis of the amplification products was performed at the end of each PCR reaction to confirm that only one specific PCR product was amplified and detected. The qPCRs were performed in triplicate for each sample. The 2-ΔΔCT method [101] was used to analyze the quantitative real-time PCR data.

References

  1. 1.

    Yoon HS, Hackett J, Ciniglia C, Pinto G, Bhattacharya D: A molecular timeline for the origin of photosynthetic eukaryotes. Mol Biol Evol. 2004, 21: 809-818. 10.1093/molbev/msh075.

  2. 2.

    Grossman AR: Paths toward algal genomics. Plant Physiol. 2005, 137: 410-427. 10.1104/pp.104.053447.

  3. 3.

    Hayden HS, Blomster J, Maggs CA, Silva PC, Stanhope MJ, Waaland JR: Linnaeus was right all along: Ulva and Enteromorpha are not distinct genera. Eur J Phycol. 2003, 38: 277-294. 10.1080/1364253031000136321.

  4. 4.

    Loughnane CJ, McIvor LM, Rindi F, Stengel DB, Guiry MD: Morphology, rbcL phylogeny and distribution of distromatic Ulva (Ulvophyceae, Chlorophyta) in Ireland and southern Britain. Phycologia. 2008, 47: 416-429. 10.2216/PH07-61.1.

  5. 5.

    McAvoy KM, Klug JL: Positive and negative effects of riverine input on the estuarine green alga Ulva intestinalis (syn. Enteromorpha intestinalis)(Linnaeus). Hydrobiologia. 2005, 545: 1-9. 10.1007/s10750-005-1923-5.

  6. 6.

    Shimada S, Yokoyama N, Arai S, Hiraoka M: Phylogeography of the genus Ulva (Ulvophyceae, Chlorophyta), with special reference to the Japanese freshwater and brackish taxa. J Appl Phycol. 2009, 20: 979-989.

  7. 7.

    Hayden HS, Waaland JR: Phylogenetic systematics of the Ulvaceae (Ulvales, Ulvophyceae) using chloroplast and nuclear DNA sequences. J Phycol. 2002, 38: 1200-1212. 10.1046/j.1529-8817.2002.01167.x.

  8. 8.

    Fletcher RL: The occurrence of ‘green tides’: a review. Marine benthic vegetation: recent changes and the effects of eutrophication. Edited by: Schramm W, Nienhuis PH. 1996, Berlin: Springerr-Verlag, 7-43.

  9. 9.

    Valiela I, McClelland J, Hauxwell J, Behr PJ, Hersh D, Foreman K: Macroalgal blooms in shallow estuaries: controls and ecophysiological and ecosystem consequences. Limnol Oceanogr. 1997, 42: 1105-1118. 10.4319/lo.1997.42.5_part_2.1105.

  10. 10.

    Ye N, Zhang X, Mao Y, Liang C, Xu D, Zou J, Zhuang Z, Qing W: ‘Green tides’ are overwhelming the coastline of our blue planet: taking the world’s largest example. Ecol Res. 2011, 26: 477-485. 10.1007/s11284-011-0821-8.

  11. 11.

    Zhang X, Xu D, Mao Y, Li Y, Xue S, Zou J, Lian W, Liang C, Zhuang Z, Wang Q, Ye N: Settlement of vegetative fragments of Ulva prolifera confirmed as an important seed source for succession of a large-scale green tide bloom. Limnol Oceanogr. 2011, 56: 233-242. 10.4319/lo.2011.56.1.0233.

  12. 12.

    Raffaelli DG, Raven JA, Poole LJ: Ecological impact of green macroalgal blooms. Oceanogr Mar Biol. 1998, 36: 97-125.

  13. 13.

    Charlier RH, Morand P, Finkl CW, Thys A: Green tides on the Brittany Coasts. Aplinkos Tyrimai Inžinerija Ir Vadyba. 2007, 3: 52-59.

  14. 14.

    Worm B, Heike K, Sommer U: Algal propagules banks modify competition, consumer and resource control on Baltic rocky shores. Oecologia. 2001, 128: 281-293. 10.1007/s004420100648.

  15. 15.

    Rudd S: Expressed sequence tags: alternative or complement to whole genome sequences?. Trends Plant Sci. 2003, 8: 321-329. 10.1016/S1360-1385(03)00131-6.

  16. 16.

    Bouck A, Vision T: The molecular ecologist’s guide to expressed sequence tags. Mol Ecol. 2007, 16: 907-924.

  17. 17.

    Nagaraj SH, Gasser RB, Ranganathan S: A hitchhiker’s guide to expressed sequence tag (EST) analysis. Brief Bioinform. 2007, 8: 6-21.

  18. 18.

    Riggins CW, Peng Y, Stewart CN, Tranela PJ: Characterization of de novo transcriptome for waterhemp (Amaranthus tuberculatus) using GS-FLX 454 pyrosequencing and its application for studies of herbicide target-site genes. Pest Manag Sci. 2010, 66: 1042-1052. 10.1002/ps.2006.

  19. 19.

    Pop M, Salzberg SL: Bioinformatics challenges of new sequencing technology. Trends Genet. 2008, 24: 142-149. 10.1016/j.tig.2007.12.006.

  20. 20.

    Stinchcombe JR, Hoekstra HE: Combining population genomics and quantitative genetics: finding the genes underlying ecologically important traits. Heredity. 2007, 100: 158-170.

  21. 21.

    Bourdon V, Naef F, Rao PH, Reuter V, Mok SC, Bosl GJ, Koul S, Murty VV, Kucherlapati RS, Chaganti RS: Genomic and expression analysis of the 12p11-p12 amplicon using EST arrays identifies two novel amplified and overexpressed genes. Cancer Res. 2002, 62: 6218-6223.

  22. 22.

    Cheung F, Win J, Lang JM, Hamilton J, Vuong H, Leach JE, Kamoun S, Levesque AC, Tisserat N, Buell CR: Analysis of the Pythium ultimum transcriptome using Sanger and pyrosequencing approaches. BMC Genomics. 2008, 9: 542-10.1186/1471-2164-9-542.

  23. 23.

    Seki M, Narusaka M, Kamiya A, et al: Functional annotation of a full-length Arabidopsis cDNA collection. Science. 2002, 296: 141-145. 10.1126/science.1071006.

  24. 24.

    Nishiyama T, Fujita T, Shin-I T, Seki M, Nishide H, Uchiyama I, Kamiya A, Carninci P, Hayashizaki Y, Shinozaki K, Kohara Y, Hasebe M: Comparative genomics of Physcomitrella patens gametophytic transcriptome and Arabidopsis thaliana: Implication for land plant evolution. Proc Natl Acad Sci USA. 2003, 100: 8007-8012. 10.1073/pnas.0932694100.

  25. 25.

    Gonzalo MJ, Oliver M, Garcia-Mas J, Monfort A, Dolcet-Sanjuan R, Katzir N, Arus P, Monforte AJ: Simple-sequence repeat markers used in merging linkage maps of melon (Cucumis melo L.). Theor Appl Genet. 2005, 110: 802-811. 10.1007/s00122-004-1814-6.

  26. 26.

    Parchman TL, Geist KS, Grahnen JA, Benkman CW, Buerkle CA: Transcriptome sequencing in an ecologically important tree species: assembly, annotation, and marker discovery. BMC Genomics. 2010, 11: 180-196. 10.1186/1471-2164-11-180.

  27. 27.

    Fei Z, Tang X, Alba R, White J, Ronning C, Martin G, Tanksley S, Giovannoni J: Comprehensive EST analysis of tomato and comparative genomics of fruit ripening. Plant J. 2004, 40: 47-59. 10.1111/j.1365-313X.2004.02188.x.

  28. 28.

    Eveland AL, McCarty DR, Koch KE: Transcript profiling by 3'-untranslated region sequencing resolves expression of gene families. Plant Physiol. 2008, 146: 32-44.

  29. 29.

    Rismani-Yazdi H, Haznedaroglu BZ, Bibby K, Peccia J: Transcriptome sequencing and annotation of the microalgae Dunaliella tertiolecta: Pathway description and gene discovery for production of next-generation biofuels. BMC Genomics. 2011, 12: 148-10.1186/1471-2164-12-148.

  30. 30.

    Vera JC, Wheat CW, Fescemyer HW, Frilander MJ, Crawford DL, Hanski I, Marden JH: Rapid transcriptome characterization for a nonmodel organism using 454 pyrosequencing. Mol Ecol. 2008, 17: 1636-1647. 10.1111/j.1365-294X.2008.03666.x.

  31. 31.

    Novaes E, Drost DR, Farmerie WG, Pappas GJ, Grattapaglia D, Sederoff RR, Kirst M: High-throughput gene and SNP discovery in Eucalyptus grandis, an uncharacterized genome. BMC Genomics. 2008, 9: 312-10.1186/1471-2164-9-312.

  32. 32.

    Kristiansson E, Asker N, Forlin L, Larsson DGJ: Characterization of the Zoarces viviparus liver transcriptome using massively parallel pyrosequencing. BMC Genomics. 2009, 10: 345-10.1186/1471-2164-10-345.

  33. 33.

    Meyer E, Aglyamova GV, Wang S, Buchanan-Carter J, Abrego D, Colbourne JK, Willis BL, Matz MV: Sequencing and de novo analysis of a coral larval transcriptome using 454 GSFlx. BMC Genomics. 2009, 10: 219-10.1186/1471-2164-10-219.

  34. 34.

    Wheat CW: Rapidly developing functional genomics in ecological model systems via 454 transcriptome sequencing. Genetica. 2008, 138: 433-451.

  35. 35.

    Henry IM, Wilkinson MD, Hernandez JM, Schwarz-Sommer Z, Grotewold E, Mandoli DF: Comparison of ESTs from juvenile and adult phases of the giant unicellular green alga Acetabularia acetabulum. BMC Plant Biol. 2004, 4: 3-10.1186/1471-2229-4-3.

  36. 36.

    Merchant SS, Prochnik SE, Vallon O, et al: The Chlamydomonas genome reveals the evolution of key animal and plant functions. Science. 2007, 318: 245-250. 10.1126/science.1143609.

  37. 37.

    Blanc G, Duncan G, Agarkova I, et al: The Chlorella variabilis NC64A Genome Reveals Adaptation to Photosymbiosis, Coevolution with Viruses, and Cryptic Sex. Plant cell. 2010, 22: 2943-2955. 10.1105/tpc.110.076406.

  38. 38.

    Derelle E, Ferraz C, Rombauts S, et al: Genome analysis of the smallest free-living eukaryote Ostreococcus tauri unveils many unique features. Proc Natl Acad Sci USA. 2006, 103: 11647-11652. 10.1073/pnas.0604795103.

  39. 39.

    Worden AZ, Lee J-H, Mock T, et al: Green Evolution and Dynamic Adaptations Revealed by Genomes of the Marine Picoeukaryotes Micromonas. Science. 2009, 324: 268-272. 10.1126/science.1167222.

  40. 40.

    Prochnik SE, Umen J, Nedelcu AM, et al: Genomic analysis of organismal complexity in the multicellular green alga Volvox carteri. Science. 2010, 329: 223-226. 10.1126/science.1188800.

  41. 41.

    Stanley MS, Perry RM, Callow JA: Analysis of expressed sequence tags from the green alga Ulva linza (Chlorophyta). J Phycol. 2005, 41: 1219-1226. 10.1111/j.1529-8817.2005.00138.x.

  42. 42.

    Niu J, Hu H, Hu S, Wang G, Peng G, Sun S: Analysis of expressed sequence tags from the Ulva prolifera (Chlorophyta). Chin J Oceanol Limnol. 2010, 28: 26-36. 10.1007/s00343-010-9120-4.

  43. 43.

    Radakovits R, Jinkerson RE, Fuerstenberg SI, Tae H, Settlage RE, Boore JL, Posewitz 1MC: Draft genome sequence and genetic transformation of the oleaginous alga Nannochloropis gaditana. Nat Commun. 2012, 10.1038/ncomms1688.

  44. 44.

    Waters ER: Molecular adaptation and the origin of land plants. Mol Phylogenet Evol. 2003, 29: 456-463. 10.1016/j.ympev.2003.07.018.

  45. 45.

    Neilson JA, Durnford DG: Structural and functional diversification of the light-harvesting complexes in photosynthetic eukaryotes. Photosynth Res. 2001, 106: 57-71.

  46. 46.

    Pascal AA, Liu Z, Broess K, van Oort B, van Amerongen H, Wang C, Horton P, Chang W, Ruban AV: Molecular basis of photoprotection and control of photosynthetic light-harvesting. Nature. 2005, 436: 134-137. 10.1038/nature03795.

  47. 47.

    Ruban AV, Berera R, Ilioaia C, van Stokkum IHM, Kennis JT, Pascal AA, van Amerongen H, Robert B, Horton P, van Grondelle R: Identification of a mechanism of photoprotective energy dissipation in higher plants. Nature. 2007, 450: 575-579. 10.1038/nature06262.

  48. 48.

    Jansson S: The light-harvesting chlorophyll a/b binding proteins. Biochim Biophys Acta. 1994, 1184: 1-19. 10.1016/0005-2728(94)90148-1.

  49. 49.

    Elrad D, Grossman AR: A genome’s-eye view of the light harvesting polypeptides of Chlamydomonas reinhardtii. Curr Genet. 2004, 45: 61-75. 10.1007/s00294-003-0460-x.

  50. 50.

    Koziol AG, Borza T, Ishida K, Keeling P, Lee RW, Durnford DG: Tracing the evolution of the light-harvesting antennae in chlorophyll a/b-containing organisms. Plant Physiol. 2007, 143: 1802-1816. 10.1104/pp.106.092536.

  51. 51.

    Sxi C, Worden AZ, Rodríguez F, Moreau H, Partensky F: New insights into the nature and phylogeny of prasinophyte antenna proteins: Ostreococcus tauri, a case study. Mo Biol Evol. 2005, 11: 2217-2230.

  52. 52.

    Boekema EJ, van Roon H, van Breemen JFL, Dekker JP: Supramolecular organization of Photosystem II and its lightharvesting antenna in partially solubilized Photosystem II membranes. Eur J Biochem. 1999, 266: 444-452. 10.1046/j.1432-1327.1999.00876.x.

  53. 53.

    Yakushevska AE, Keegstra W, Boekema EJ, Dekker JP, Andersson J, Jansson S, Ruban AV, Horton P: The structure of Photosystem II in Arabidopsis: localization of the CP26 and CP29 antenna complexes. Biochemistry. 2003, 42: 608-613. 10.1021/bi027109z.

  54. 54.

    Horton P, Johnson MP, Perez-Bueno ML, Kiss AZ, Ruban AV: Photosynthetic acclimation: Does the dynamic structure and macro-organisation of photosystem II in higher plant grana membranes regulate light harvesting states?. FEBS J. 2008, 275: 1069-1079. 10.1111/j.1742-4658.2008.06263.x.

  55. 55.

    Schultes NP, Peterson RB: Phylogeny-directed structural analysis of the Arabidopsis PsbS protein. Biochem Biophy Res Co. 2007, 355: 464-470. 10.1016/j.bbrc.2007.01.173.

  56. 56.

    Li XP, Björkman O, Shih C, Grossman AR, Rosenquist M, Jansson S, Niyogi KK: A pigment-binding protein essential for regulation of photosynthetic light harvesting. Nature. 2000, 403: 391-395. 10.1038/35000131.

  57. 57.

    Li XP, Müller-Moulé P, Gilmore AM, Niyogi KK: PsbS dependent enhancement of feedback de-excitation protects Photosystem II from photoinhibition. Proc Natl Acad Sci USA. 2002, 99: 15222-15227. 10.1073/pnas.232447699.

  58. 58.

    Li XP, Gilmore AM, Caffari S, Bassi R, Golan T, Kramer D, Niyogi KK: Regulation of photosynthetic light harvesting involves intrathylakoid lumen pH sensing by the PsbS protein. J Biol Chem. 2004, 279: 22866-22874. 10.1074/jbc.M402461200.

  59. 59.

    Bonente G, Howes B, Caffarri S, Smulevich G, Bassi R: The occurrence of the psbS gene product in Chlamydomonas reinhardtii and in other photosynthetic organisms and its correlation with Energy Quenching. Photochem Photobiol. 2008, 84: 1359-1370. 10.1111/j.1751-1097.2008.00456.x.

  60. 60.

    Kiss AZ, Ruban AV, Horton P: The PsbS protein controls the organization of the photosystem II antenna in higher plant thylakoid membranes. J Biol Chem. 2008, 283: 3972-3978.

  61. 61.

    Betterle N, Ballottair M, Zorzan S, et al: Light-induced dissociation of an antenna hetero-oligomer is needed for nonphotochemical quenching induction. J Biol Chem. 2009, 22: 15255-15266.

  62. 62.

    Bonente G, Passarini F, Cazzaniga S, Mancone C, Buia MC, Tripodi M, Bassi R, Caffarri S: The occurrence of the psbS gene product in Chlamydomonas reinhardtii and in other photosynthetic organisms and its correlation with energy quenching. Photochem Photobiol. 2008, 84: 1359-1370. 10.1111/j.1751-1097.2008.00456.x.

  63. 63.

    Peers G, Truong TB, Ostendorf E, Busch A, Elrad D, Grossman AR, Hippler M, Niyogi KK: An ancient light-harvesting protein is critical for the regulation of algal photosynthesis. Nature. 2009, 462: 518-522. 10.1038/nature08587.

  64. 64.

    Bonente G, Ballottari M, Truong TB, Morosinotto T, Ahn TK, et al: Analysis of LhcSR3, a Protein Essential for Feedback De-Excitation in the Green Alga Chlamydomonas reinhardtii. PLoS Biol. 2011, 9: e1000577-10.1371/journal.pbio.1000577.

  65. 65.

    Alboresi A, Gerotto C, Giacometti GM, Bassi R, Morosinotto T: Physcomitrella patens mutants affected on heat dissipation clarify the evolution of photoprotection mechanisms upon land colonization. Proc Natl Acad Sci USA. 2010, 107: 11128-11133. 10.1073/pnas.1002873107.

  66. 66.

    Raven JA, Cockell CS, Rocha CLL: The evolution of inorganic carbon concentrating mechanisms in photosynthesis. Phil Trans R Soc B. 2008, 1504: 2641-2650.

  67. 67.

    Giordano M, Beardall J, Raven JA: CO2 Concentrating mechanisms in algae: mechanisms, evironmental mdulation, and Evolution. Annu Rev of Plant Biol. 2005, 56: 99-131. 10.1146/annurev.arplant.56.032604.144052.

  68. 68.

    Palmqvist K, Yu JW, Badger MR: Carbonic anhydrase activity and inorganic carbon fluxes in low- and high-Ci cells of Chlamydomonas reinhardtii and Scenedesmus obliquus. Physiol Plantarum. 1994, 90: 537-547. 10.1111/j.1399-3054.1994.tb08812.x.

  69. 69.

    Yamano T, Fukuzawa H: Carbon-concentrating mechanism in a green alga Chlamydomonas reinhardtii, revealed by transcriptome analyses. J Basic Microb. 2009, 49: 42-51. 10.1002/jobm.200800352.

  70. 70.

    Sage RF: The evolution of C4 photosynthesis. New Phytol. 2004, 161: 341-370. 10.1111/j.1469-8137.2004.00974.x.

  71. 71.

    Giordano M, Beardall J, Raven JA: CO2 concentrating mechanisms in algae: mechanisms, environmental modulation, and evolution. Annu Rev Plant Biol. 2005, 56: 99-131. 10.1146/annurev.arplant.56.032604.144052.

  72. 72.

    Kremer BP: Metabolic implications of non-photosynthetic carbon fixation in brown macroalgae. Phycologia. 1981, 20: 242-250. 10.2216/i0031-8884-20-3-242.1.

  73. 73.

    Cock JM, et al: The Ectocarpus genome and the independent evolution of multicellularity in brown algae. Nature. 2010, 465: 617-621. 10.1038/nature09016.

  74. 74.

    Kerby NW, Raven JA: Transport and fixation of inorganic carbon by marine algae. Adv Bot Res. 1985, 11: 71-123.

  75. 75.

    Lucas WJ, Berry JA: Inorganic carbon transport in aquatic photosynthetic organisms. Physiol Plant. 1985, 65: 539-543. 10.1111/j.1399-3054.1985.tb08687.x.

  76. 76.

    Reiskind JB, Bowes G: The role of phosphoenolpyruvate carboxykinase in a marine macroalga with C4-like photosynthetic characteristics. Proc Natl Acad Sci USA. 1991, 88: 2883-2887. 10.1073/pnas.88.7.2883.

  77. 77.

    Fan X, Fang Y, Hu S, Wang G: Generation and analysis of 5318 expressed sequence tags from the filamentous sporophyte of Porphyra haitanensis (Rhodophyta). J Phycol. 2007, 43: 1287-1294. 10.1111/j.1529-8817.2007.00415.x.

  78. 78.

    Xu J, Fan X, Zhang X, Xu D, Mou S, et al: Evidence of Coexistence of C3 and C4 Photosynthetic Pathways in a Green-Tide-Forming Alga, Ulva prolifera. PLoS One. 2012, 7 (5): e37438-10.1371/journal.pone.0037438.

  79. 79.

    Brault D, Quéguiner B: Aquaculture-A Biotechnology in Progress. Effect of inorganic and organic nitrogen sources on growth of Ulva gigantean (Kutzing) Bliding. Edited by: Pauw ND, Jaspers E, Ackefors H, Wilkins N. 1989, Bredene: European Aquaculture Society, 425-430.

  80. 80.

    Coon M: Cytochrome P450: nature's most versatile biological catalyst. Annu Rev Pharmacol. 2005, 45: 1-25. 10.1146/annurev.pharmtox.45.120403.100030.

  81. 81.

    Initiative TAG: Analysis of the genome sequence of the flowering plant Arabidopsis thaliana. Nature. 2000, 408: 796-815. 10.1038/35048692.

  82. 82.

    Tian L, Musetti V, Kim J, Magallanes-Lundback M, DellaPenna D: The Arabidopsis LUT1 locus encodes a member of the cytochrome p450 family that is required for carotenoid epsilon-ring hydroxylation activity. Proc Natl Acad Sci USA. 2004, 101: 402-407. 10.1073/pnas.2237237100.

  83. 83.

    Kim J, Smith J, Tian L, Dellapenna D: The evolution and function of carotenoid hydroxylases in Arabidopsis. Plant Cell Physiol. 2009, 50: 463-479. 10.1093/pcp/pcp005.

  84. 84.

    Mulderij G, Mooij WM, Smolders AJP, Van Donk E: Allelopathic inhibition of phytoplankton by exudates from Stratiotes aloides. Aquat Bot. 2005, 82: 284-296. 10.1016/j.aquabot.2005.04.001.

  85. 85.

    Wang Y, Yu ZM, Song XX, Tang XX, Zhang SD: Effects of macroalgae Ulva pertusa (Chlorophyta) and Gracilaria lemaneiformis (Rhodophyta) on growth of four species of bloom-forming dinoflagellates. Aquat Bot. 2007, 86: 139-147. 10.1016/j.aquabot.2006.09.013.

  86. 86.

    Huo YZ, Zhang JH, Xu SN, Tian QT, Zhang YJ, He PM: Effects of seaweed Gracilaria verrucosa on the growth of microalgae: A case study in the laboratory and in an enclosed sea of Hangzhou Bay, China. Harmful Algae. 2011, 10: 411-418. 10.1016/j.hal.2011.02.003.

  87. 87.

    Jin Q, Dong S, Wang C: Allelopathic growth inhibition of Prorocentrum micans (Dinophyta) by Ulva pertusa and Ulva linza (Chlorophyta) in laboratory cultures. Eur J Phycol. 2005, 40: 1-37. 10.1080/09670260400009882.

  88. 88.

    Mungur R, Glass AD, Goodenow DB, Lightfoot DA: Metabolite fingerprinting in transgenic Nicotiana tabacum altered by the Escherichia coli glutamate dehydrogenase gene. J Biomed Biotechnol. 2005, 2: 198-214.

  89. 89.

    Grabowska A, Nowicki M, Kwinta J: Glutamate dehydrogenase of the germinating triticale seeds: gene expression, activity distribution and kinetic characteristics. Acta Physiol Plant. 2011, 33: 1981-1990. 10.1007/s11738-011-0801-1.

  90. 90.

    Lightfoot DA, Bernhardt K, Mungur R, Nolte S, Ameziane R, Colter A, Jones K, Iqbal MJ, Varsa E, Young B: Improved drought tolerance of transgenic Zea mays plants that express the glutamate dehydrogenase gene (gdhA) of E. coli. Euphytica. 2007, 156: 103-116. 10.1007/s10681-007-9357-y.

  91. 91.

    Andersson JO, Roger AJ: Evolution of glutamate dehydrogenase genes: evidence for lateral gene transfer within and between prokaryotes and eukaryotes. BMC Evol Biol. 2003, 3: 14-10.1186/1471-2148-3-14.

  92. 92.

    Asada K, Takahashi M: Production and scavenging of active oxygen in photosynthesis. Photoinhibition. Edited by: Kyle DJ, Osmond CB, Arntzen CJ. 1987, Amsterdam: Elsevier, 227-287.

  93. 93.

    Bowler C, Montagu MV, Inzé D: Superoxide dismutases and stress tolerance. Annu Rev Plant Physiol Plant Mol Biol. 1992, 43: 83-116. 10.1146/annurev.pp.43.060192.000503.

  94. 94.

    Willekens H, Chamnongpol S, Davey M, Schraudner M, Langebartels C, Montagu MV, Inzé D, Camp WV: Catalase is a sink for H2O2 and is indispensable for stress defence in C3 plants. EMBO J. 1997, 16: 4806-4816. 10.1093/emboj/16.16.4806.

  95. 95.

    Wu TM, Lee TM: Regulation of activity and gene expression of antioxidant enzymes in Ulva fasciata Delile (Ulvales, Chlorophyta) in response to excess copper. Phycologia. 2008, 47: 346-360. 10.2216/PH07-77.1.

  96. 96.

    Scharf KD, Seddigue M, Vierling E: The expanding family of Arabidopsis thaliana small heat stress proteins and a new family of proteins containing a-crystallin domains (ACD proteins). Cell Stress Chaperon. 2001, 6: 225-237. 10.1379/1466-1268(2001)006<0225:TEFOAT>2.0.CO;2.

  97. 97.

    Young JC, Moarefi I, Hartl FU: Hsp90: a specialized but essential protein folding tool. J Cell Biol. 2001, 154: 267-273. 10.1083/jcb.200104079.

  98. 98.

    Camon E, Magrane M, Barrell D, Lee V, Dimmer E, Maslen J, Binns D, Harte N, Lopez R, Apweiler R: The Gene Ontology Annotation (GOA) Database: sharingknowledge in Uniprot with Gene Ontology. Nucleic Acids Res. 2004, 32: 262-266. 10.1093/nar/gkh021.

  99. 99.

    Huson DH, Auch A, Qi J, Schuster SC: Megan Analysis of Metagenome Data. Genome Res. 2007, 17: 377-386. 10.1101/gr.5969107.

  100. 100.

    Kanehisa M, Goto S, Furumichi M, Tanabe M, Hirakawa M: KEGG for representation and analysis of molecular networks involving diseases and drugs. Nucleic Acids Res. 2010, 38: 355-360. 10.1093/nar/gkp896.

  101. 101.

    Livak KJ, Schmittgen TD: Analysis of relative gene expression data using real-time quantitative PCR and the 2-ΔΔCTmethod. Methods. 2001, 25: 402-408. 10.1006/meth.2001.1262.

Download references

Acknowledgments

We acknowledge the supercomputing support for the genomics analysis by Beijing Institutes of Life Science, Chinese Academy of Sciences (BIOLS). We appreciate the reviewers very much for their positive and constructive comments and suggestions on our manuscript. This work was supported by the National Natural Science Foundation of China (41176153, 31200187), Hi-Tech Research and Development Program (863) of China (2012AA052103), Shandong Science and Technology plan project (2011GHY11528).

Author information

Affiliations

  1. Yellow Sea Fisheries Research Institute, Chinese Academy of Fishery Sciences, Qingdao, 266071, China
    • Xiaowen Zhang
    • , Naihao Ye
    • , Shanli Mou
    • , Xiao Fan
    • , Jianfang Xu
    • , Dong Xu
    •  & Zhimeng Zhuang
  2. Qingdao University of Science > Technology, Qingdao, 266042, China
    • Chengwei Liang
  3. Key Laboratory of Marine Bioactive Substance, The First Institute of Oceanography, State Oceanic administration (SOA), Qingdao, 266061, China
    • Jianfang Xu
Authors
  1. Search for Xiaowen Zhang in:
  2. Search for Naihao Ye in:
  3. Search for Chengwei Liang in:
  4. Search for Shanli Mou in:
  5. Search for Xiao Fan in:
  6. Search for Jianfang Xu in:
  7. Search for Dong Xu in:
  8. Search for Zhimeng Zhuang in:

Corresponding author

Correspondence to Naihao Ye.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors’ contributions

XZ, NY, and ZZ conceived and designed the study; SM, XF and JX conducted culturing and molecular biology laboratory work; XZ, CL and DX undertook data analyses; all authors contributed to the manuscript. All authors read and approved the final manuscript.

Electronic supplementary material

Additional file 1: Table S1. Annotation of Ulva genes occurring in Chlamydomonas reinhardtii, but not in Ostreococcus tauri. (XML 3 MB)

Additional file 2: Table S2. Eighteen Ulva special isotigs conserved in land plants. (DOC 54 KB)

Additional file 3: Table S3. a. Comparison of Ulva KEGG annotation analyzed by Reads and Isotigs respectively. b. Comparison of KEGG annotation among Ulva linza (isotigs), Chlamydomonas reinhardtii, Volvox carteri, Physcomitrella patens and Arabidopsis thaliana on three different levels. (DOC 51 KB)

Additional file 4: Figure S1. Phylogenetic analysis of Lhcb proteins in six organisms. The phylogenetic tree was constructed by the neighbor-joining algorithm of the MEGA 4.0 program. A total of 1,000 bootstrap replicates were performed. Bootstrap values (1000 replicates) >50% are indicated on the branch. (JPEG 415 KB)

Additional file 5: Table S4. Putative genes encoding the enzymes required for a C4-like CCM. (DOC 40 KB)

Additional file 6: Figure S2. Phylogenetic analysis of putative ammonium transporters in U. linza. The phylogenetic tree was constructed by the neighbor-joining algorithm of the MEGA 4.0 program. A total of 1,000 bootstrap replicates were performed. The prokaryote-like ammonium transporters found in U. linza was showed in red color. (JPEG 529 KB)

Additional file 7: Figure S3. Phylogenetic analysis of putative phosphate transporters in U. linza. The phylogenetic tree was constructed by the neighbor-joining algorithm of the MEGA 4.0 program. A total of 1,000 bootstrap replicates were performed. (JPEG 617 KB)

Additional file 8: Figure S4. Phylogenetic analysis of putative sulfate transporters in U. linza. The phylogenetic tree was constructed by the neighbor-joining algorithm of the MEGA 4.0 program. A total of 1,000 bootstrap replicates were performed. (JPEG 467 KB)

Additional file 9: Figure S5. Phylogenetic analysis of putative P450s in U. linza. The phylogenetic tree was constructed by the neighbor-joining algorithm of the MEGA 4.0 program. A total of 1,000 bootstrap replicates were performed. (JPEG 823 KB)

Additional file 10: Figure S6. Phylogenetic analysis of putative Glutamate dehydrogenases (GDH) in U. linza. The phylogenetic tree was constructed by the neighbor-joining algorithm of the MEGA 4.0 program. A total of 1,000 bootstrap replicates were performed. (JPEG 571 KB)

Additional file 11: Figure S7. Phylogenetic analysis of putative ascorbate peroxidase (APX) genes in U. linza. The phylogenetic tree was constructed by the neighbor-joining algorithm of the MEGA 4.0 program. A total of 1,000 bootstrap replicates were performed. (JPEG 604 KB)

Additional file 12: Figure S8. Phylogenetic analysis of putative catalase (CAT) genes in U. linza. The phylogenetic tree was constructed by the neighbor-joining algorithm of the MEGA 4.0 program. A total of 1,000 bootstrap replicates were performed. (JPEG 597 KB)

Additional file 13: Table S5. Putative Heat-shock proteins found in U. linza. Three genes encoding Hsp70, Hsp90 and Hsp100 that not conserved in green algae and plants were shown in yellow color. (DOC 80 KB)

Additional file 14: Figure S9. Phylogenetic analysis of putative Hsp20 proteins in U. linza. The phylogenetic tree was constructed by the neighbor-joining algorithm of the MEGA 4.0 program. A total of 1,000 bootstrap replicates were performed. (JPEG 771 KB)

Additional file 15: Table S6. Comparison of Heat-shock proteins among Ulva linza, Micromonas sp. RCC299, Ostreococcus tauri, Chlorella variabilis NC64A, Chlamydomonas reinhardtii, Volvox carteri and Arabidopsis thaliana. (DOC 38 KB)

Additional file 16: Table S7. Primers used in the real-time RT-PCR analysis. (DOC 29 KB)

Authors’ original submitted files for images

Rights and permissions

Reprints and Permissions

About this article

Cite this article

Zhang, X., Ye, N., Liang, C. et al. De novo sequencing and analysis of the Ulva linza transcriptome to discover putative mechanisms associated with its successful colonization of coastal ecosystems. BMC Genomics 13, 565 (2012). https://doi.org/10.1186/1471-2164-13-565

Download citation

Keywords

  • Gene Ontology
  • Green Alga
  • Land Plant
  • Coastal Ecosystem
  • Unicellular Green Alga

BMC Genomics

ISSN: 1471-2164

Contact us