EST assembly supported by a draft genome sequence: an analysis of the Chlamydomonas reinhardtii transcriptome

Omics Application of Bio-Hydrogen Production Through Green Alga Chlamydomonas reinhardtii

This article summarizes the current knowledge regarding omics approaches, which include genomics, transcriptomics, proteomics and metabolomics, in the context of bio-hydrogen production in Chlamydomonas reinhardtii. In this paper, critical genes (HydA1, Hyd A2, Sulp, Tla1, Sta7, PFL1) involved in H₂ metabolism were identified and analyzed for their function in H₂ accumulation. Furthermore, the advantages of gene microarrays and RNA-seq were compared, as well as their applications in transcriptomic analysis of H₂ production. Moreover, as a useful tool, proteomic analysis could identify different proteins that participate in H₂ metabolism. This review provides fundamental theory and an experimental basis for H₂ production, and further research effort is needed in this field.

MicroRNAs modulate adaption to multiple abiotic stresses in Chlamydomonas reinhardtii

MicroRNAs play an important role in abiotic stress responses in higher plants and animals, but their role in stress adaptation in algae remains unknown. In this study, the expression of identified and putative miRNAs in Chlamydomonas reinhardtii was assessed using quantitative polymerase chain reaction; some of the miRNAs (Cre-miR906-3p) were up-regulated, whereas others (Cre-miR910) were down-regulated when the species was subjected to multiple abiotic stresses. With degradome sequencing data, we also identified ATP4 (the d-subunit of ATP synthase) and NCR2 (NADPH: cytochrome P450 reductase) as one of the several targets of Cre-miR906-3p and Cre-miR910, respectively. Q-PCR data indicated that ATP4, which was expressed inversely in relation to Cre-miR906-3p under stress conditions. Overexpressing of Cre-miR906-3p enhanced resistance to multiple stresses; conversely, overexpressing of ATP4 produced the opposite effect. These data of Q-PCR, degradome sequencing and adaptation of overexpressing lines indicated that Cre-miR906-3p and its target ATP4 were a part of the same pathway for stress adaptation. We found that Cre-miR910 and its target NCR2 were also a part of this pathway. Overexpressing of Cre-miR910 decreased, whereas that of NCR2 increased the adaption to multiple stresses. Our findings suggest that the two classes of miRNAs synergistically mediate stress adaptation in algae.

Do housekeeping genes exist?

The searching of human housekeeping (HK) genes has been a long quest since the emergence of transcriptomics, and is instrumental for us to understand the structure of genome and the fundamentals of biological processes. The resolved genes are frequently used in evolution studies and as normalization standards in quantitative gene-expression analysis. Within the past 20 years, more than a dozen HK-gene studies have been conducted, yet none of them sampled human tissues completely. We believe an integration of these results will help remove false positive genes owing to the inadequate sampling. Surprisingly, we only find one common gene across 15 examined HK-gene datasets comprising 187 different tissue and cell types. Our subsequent analyses suggest that it might not be appropriate to rigidly define HK genes as expressed in all tissue types that have diverse developmental, physiological, and pathological states. It might be beneficial to use more robustly identified HK functions for filtering criteria, in which the representing genes can be a subset of genome. These genes are not necessarily the same, and perhaps need not to be the same, everywhere in our body.

Genomics of Volvocine Algae

Volvocine algae are a group of chlorophytes that together comprise a unique model for evolutionary and developmental biology. The species Chlamydomonas reinhardtii and Volvox carteri represent extremes in morphological diversity within the Volvocine clade. Chlamydomonas is unicellular and reflects the ancestral state of the group, while Volvox is multicellular and has evolved numerous innovations including germ-soma differentiation, sexual dimorphism, and complex morphogenetic patterning. The Chlamydomonas genome sequence has shed light on several areas of eukaryotic cell biology, metabolism and evolution, while the Volvox genome sequence has enabled a comparison with Chlamydomonas that reveals some of the underlying changes that enabled its transition to multicellularity, but also underscores the subtlety of this transition. Many of the tools and resources are in place to further develop Volvocine algae as a model for evolutionary genomics.

ANALYSIS OF EXPRESSED SEQUENCE TAGS FROM THE GREEN ALGA DUNALIELLA SALINA (CHLOROPHYTA)(1)

The unicellular green alga Dunaliella salina (Dunal) Teodor. is a novel model photosynthetic eukaryote for studying photosystems, high salinity acclimation, and carotenoid accumulation. In spite of such significance, there have been limited studies on the Dunaliella genome transcriptome and proteome. To further investigate D. salina, a cDNA library was constructed and sequenced. Here, we present the analysis of the 2,282 expressed sequence tags (ESTs) generated together with 3,990 ESTs from dbEST. A total of 4,148 unique sequences (UniSeqs) were identified, of which 56.1% had sequence similarity with Uniprot entries, suggesting that a large number of unique genes may be harbored by Dunaliella. Additionally, protein family domains were identified to further characterize these sequences. Then, we also compared EST sequences with different complete eukaryotic genomes from several animals, plants, and fungi. We observed notable differences between D. salina and other organisms. This EST collection and its annotation provided a significant resource for basic and applied research on D. salina and laid the foundation for a systematic analysis of the transcriptome basis of green algae development and diversification.

Construction and evaluation of a whole genome microarray of Chlamydomonas reinhardtii

Background

Chlamydomonas reinhardtii is widely accepted as a model organism regarding photosynthesis, circadian rhythm, cell mobility, phototaxis, and biotechnology. The complete annotation of the genome allows transcriptomic studies, however a new microarray platform was needed. Based on the completed annotation of Chlamydomonas reinhardtii a new microarray on an Agilent platform was designed using an extended JGI 3.1 genome data set which included 15000 transcript models.

Results

In total 44000 probes were determined (3 independent probes per transcript model) covering 93% of the transcriptome. Alignment studies with the recently published AUGUSTUS 10.2 annotation confirmed 11000 transcript models resulting in a very good coverage of 70% of the transcriptome (17000). Following the estimation of 10000 predicted genes in Chlamydomonas reinhardtii our new microarray, nevertheless, covers the expected genome by 90-95%.

Conclusions

To demonstrate the capabilities of the new microarray, we analyzed transcript levels for cultures grown under nitrogen as well as sulfate limitation, and compared the results with recently published microarray and RNA-seq data. We could thereby confirm previous results derived from data on nutrient-starvation induced gene expression of a group of genes related to protein transport and adaptation of the metabolism as well as genes related to efficient light harvesting, light energy distribution and photosynthetic electron transport.

Phylogenomic analysis of the Chlamydomonas genome unmasks proteins potentially involved in photosynthetic function and regulation

Chlamydomonas reinhardtii, a unicellular green alga, has been exploited as a reference organism for identifying proteins and activities associated with the photosynthetic apparatus and the functioning of chloroplasts. Recently, the full genome sequence of Chlamydomonas was generated and a set of gene models, representing all genes on the genome, was developed. Using these gene models, and gene models developed for the genomes of other organisms, a phylogenomic, comparative analysis was performed to identify proteins encoded on the Chlamydomonas genome which were likely involved in chloroplast functions (or specifically associated with the green algal lineage); this set of proteins has been designated the GreenCut. Further analyses of those GreenCut proteins with uncharacterized functions and the generation of mutant strains aberrant for these proteins are beginning to unmask new layers of functionality/regulation that are integrated into the workings of the photosynthetic apparatus.

Critical role of 7,8-didemethyl-8-hydroxy-5-deazariboflavin for photoreactivation in Chlamydomonas reinhardtii

DNA photolyases use two noncovalently bound chromophores to catalyze photoreactivation, the blue light-dependent repair of DNA that has been damaged by ultraviolet light. FAD is the catalytic chromophore for all photolyases and is essential for photoreactivation. The identity of the second chromophore is often 7,8-didemethyl-8-hydroxy-5-deazariboflavin (FO). Under standard light conditions, the second chromophore is considered nonessential for photoreactivation because DNA photolyase bound to only FAD is sufficient to catalyze the repair of UV-damaged DNA. phr1 is a photoreactivation-deficient strain of Chlamydomonas. In this work, the PHR1 gene of Chlamydomonas was cloned through molecular mapping and shown to encode a protein similar to known FO synthases. Additional results revealed that the phr1 strain was deficient in an FO-like molecule and that this deficiency, as well as the phr1 photoreactivation deficiency, could be rescued by transformation with DNA constructs containing the PHR1 gene. Furthermore, expression of a PHR1 cDNA in Escherichia coli produced a protein that generated a molecule with characteristics similar to FO. Together, these results indicate that the Chlamydomonas PHR1 gene encodes an FO synthase and that optimal photoreactivation in Chlamydomonas requires FO, a molecule known to serve as a second chromophore for DNA photolyases.

Enzymatic properties of the ferredoxin-dependent nitrite reductase from Chlamydomonas reinhardtii. Evidence for hydroxylamine as a late intermediate in ammonia production

The ferredoxin-dependent nitrite reductase from the green alga Chlamydomonas reinhardtii has been cloned, expressed in Escherichia coli as a His-tagged recombinant protein, and purified to homogeneity. The spectra, kinetic properties and substrate-binding parameters of the C. reinhardtii enzyme are quite similar to those of the ferredoxin-dependent spinach chloroplast nitrite reductase. Computer modeling, based on the published structure of spinach nitrite reductase, predicts that the structure of C. reinhardtii nitrite reductase will be similar to that of the spinach enzyme. Chemical modification studies and the ionic-strength dependence of the enzyme's ability to interact with ferredoxin are consistent with the involvement of arginine and lysine residues on C. reinhardtii nitrite reductase in electrostatically-stabilized binding to ferredoxin. The C. reinhardtii enzyme has been used to demonstrate that hydroxylamine can serve as an electron-accepting substrate for the enzyme and that the product of hydroxylamine reduction is ammonia, providing the first experimental evidence for the hypothesis that hydroxylamine, bound to the enzyme, can serve as a late intermediate during the reduction of nitrite to ammonia catalyzed by the enzyme.

The transcriptome of the colonial marine hydroid Hydractinia echinata

An increasing amount of expressed sequence tag (EST) and genomic data, predominantly for the cnidarians Acropora, Hydra and Nematostella, reveals that cnidarians have a high genomic complexity, despite being one of the morphologically simplest multicellular animals. Considering the diversity of cnidarians, we performed an EST project on the hydroid Hydractinia echinata, to contribute towards a broader coverage of this phylum. After random sequencing of almost 9000 clones, EST characterization revealed a broad diversity in gene content. Corroborating observations in other cnidarians, Hydractinia sequences exhibited a higher sequence similarity to vertebrates than to ecdysozoan invertebrates. A significant number of sequences were hitherto undescribed in metazoans, suggesting that these may be either cnidarian innovations or ancient genes lost in the bilaterian genomes analysed so far. However, we cannot rule out some degree of contamination from commensal bacteria. The identification of unique Hydractinia sequences emphasizes that the acquired genomic information generated so far is not large enough to be representative of the highly diverse cnidarian phylum. Finally, a database was created to store all the acquired information (http://www.mchips.org/hydractinia_echinata.html).

ESTPiper--a web-based analysis pipeline for expressed sequence tags

BACKGROUND

EST sequencing projects are increasing in scale and scope as the genome sequencing technologies migrate from core sequencing centers to individual research laboratories. Effectively, generating EST data is no longer a bottleneck for investigators. However, processing large amounts of EST data remains a non-trivial challenge for many. Web-based EST analysis tools are proving to be the most convenient option for biologists when performing their analysis, so these tools must continuously improve on their utility to keep in step with the growing needs of research communities. We have developed a web-based EST analysis pipeline called ESTPiper, which streamlines typical large-scale EST analysis components.

RESULTS

The intuitive web interface guides users through each step of base calling, data cleaning, assembly, genome alignment, annotation, analysis of gene ontology (GO), and microarray oligonucleotide probe design. Each step is modularized. Therefore, a user can execute them separately or together in batch mode. In addition, the user has control over the parameters used by the underlying programs. Extensive documentation of ESTPiper's functionality is embedded throughout the web site to facilitate understanding of the required input and interpretation of the computational results. The user can also download intermediate results and port files to separate programs for further analysis. In addition, our server provides a time-stamped description of the run history for reproducibility. The pipeline can also be installed locally, allowing researchers to modify ESTPiper to suit their own needs.

CONCLUSION

ESTPiper streamlines the typical process of EST analysis. The pipeline was initially designed in part to support the Daphnia pulex cDNA sequencing project. A web server hosting ESTPiper is provided at http://estpiper.cgb.indiana.edu/ to now support projects of all size. The software is also freely available from the authors for local installations.

A proteomic survey of Chlamydomonas reinhardtii mitochondria sheds new light on the metabolic plasticity of the organelle and on the nature of the alpha-proteobacterial mitochondrial ancestor

Mitochondria play a key role in the life and death of eukaryotic cells, yet the full spectrum of mitochondrial functions is far from being fully understood, especially in photosynthetic organisms. To advance our understanding of mitochondrial functions in a photosynthetic cell, an extensive proteomic survey of Percoll-purified mitochondria from the metabolically versatile, hydrogen-producing green alga Chlamydomonas reinhardtii was performed. Different fractions of purified mitochondria from Chlamydomonas cells grown under aerobic conditions were analyzed by nano-liquid chromatography-electrospray ionization-mass spectrometry after protein separation on sodium dodecyl sulfate polyacrylamide gel electrophoresis or on blue-native polyacrylamide gel electrophoresis. Of the 496 nonredundant proteins identified, 149 are known or predicted to reside in other cellular compartments and were thus excluded from the molecular and evolutionary analyses of the Chlamydomonas proteome. The mitochondrial proteome of the photosynthetic alga reveals important lineage-specific differences with other mitochondrial proteomes, reflecting the high metabolic diversity of the organelle. Some mitochondrial metabolic pathways in Chlamydomonas appear to combine typical mitochondrial enzymes and bacterial-type ones, whereas others are unknown among mitochondriate eukaryotes. The comparison of the Chlamydomonas proteins to their identifiable homologs predicted from 354 sequenced genomes indicated that Arabidopsis is the most closely related nonalgal eukaryote. Furthermore, this phylogenomic analysis shows that free-living alpha-proteobacteria from the metabolically versatile orders Rhizobiales and Rhodobacterales better reflect the gene content of the ancestor of the chlorophyte mitochondria than parasitic alpha-proteobacteria with reduced and specialized genomes.

Update of the Diatom EST Database: a new tool for digital transcriptomics

The Diatom Expressed Sequence Tag (EST) Database was constructed to provide integral access to ESTs from these ecologically and evolutionarily interesting microalgae. It has now been updated with 130 000 Phaeodactylum tricornutum ESTs from 16 cDNA libraries and 77 000 Thalassiosira pseudonana ESTs from seven libraries, derived from cells grown in different nutrient and stress regimes. The updated relational database incorporates results from statistical analyses such as log-likelihood ratios and hierarchical clustering, which help to identify differentially expressed genes under different conditions, and allow similarities in gene expression in different libraries to be investigated in a functional context. The database also incorporates links to the recently sequenced genomes of P. tricornutum and T. pseudonana, enabling an easy cross-talk between the expression pattern of diatom orthologs and the genome browsers. These improvements will facilitate exploration of diatom responses to conditions of ecological relevance and will aid gene function identification of diatom-specific genes and in silico gene prediction in this largely unexplored class of eukaryotes. The updated Diatom EST Database is available at http://www.biologie.ens.fr/diatomics/EST3.

Expressed sequence tags with cDNA termini: previously overlooked resources for gene annotation and transcriptome exploration in Chlamydomonas reinhardtii

Many of Chlamydomonas reinhardtii expressed sequence tags (ESTs) in GenBank dbEST and community EST assemblies were either over- or undertrimmed in terms of their cDNA termini, which are defined as the diagnostic sequence elements that delineate 3'/5' ends of mRNA transcripts. Overtrimming represents a loss of directional, positional, and structural information of transcript ends whereas undertrimming causes unclean spurious sequences retained in ESTs that exert deleterious impacts on downstream EST-based applications. We examined 309,278 raw EST sequencing trace files of C. reinhardtii and found that only 57% had cDNA termini that matched the expected structures specified in their cDNA library constructions while satisfying our minimum length requirement for their final clean sequences. Using GMAP, 156,963 individual ESTs were mapped to the genome successfully, with their in silico-verified cDNA termini anchored to the genome. Our data analysis suggested strong macro- and microheterogeneity of 3'/5' end positions of individual transcripts derived from the same genes in C. reinhardtii. This work annotating differential ends of individual transcripts in the draft genome presents the research community with a new stream of data that will facilitate accurate determination of gene structures, genome annotation, and exploration of the transcriptome and mRNA metabolism in C. reinhardtii.

High quality metabolomic data for Chlamydomonas reinhardtii

The green eukaryote alga Chlamydomonas reinhardtii is a unicellular model to study control of metabolism in a photosynthetic organism. We here present method improvements for metabolite profiling based on GC-TOF mass spectrometry focusing on three parameters: quenching and cell disruption, extract solvent composition and metabolite annotation. These improvements facilitate using smaller cell numbers and hence, smaller culture volumes which enable faster and more precise sampling techniques that eventually lead to a higher number of samples that can be processed, e.g. for time course experiments. Quenching of metabolism was achieved by mixing 1 ml of culture to 1 ml of -70 degrees C cold 70% methanol. After centrifugation, cells were lyophilized and disrupted by milling using 2-6E6 lyophilized cells, around 500-fold less than previously reported. Glass beads were compared to metal balls for milling, and five different extraction solvents were tested. Additionally, all peaks were annotated in an automated way using the GC-TOF database BinBase instead of manual investigation of a single reference chromatogram. Median precision of analysis was used to decide for the eventual procedure which was applied to a proof-of-principle study of time dependent changes of metabolism under standard conditions.

The Chlamydomonas genome reveals the evolution of key animal and plant functions

Chlamydomonas reinhardtii is a unicellular green alga whose lineage diverged from land plants over 1 billion years ago. It is a model system for studying chloroplast-based photosynthesis, as well as the structure, assembly, and function of eukaryotic flagella (cilia), which were inherited from the common ancestor of plants and animals, but lost in land plants. We sequenced the approximately 120-megabase nuclear genome of Chlamydomonas and performed comparative phylogenomic analyses, identifying genes encoding uncharacterized proteins that are likely associated with the function and biogenesis of chloroplasts or eukaryotic flagella. Analyses of the Chlamydomonas genome advance our understanding of the ancestral eukaryotic cell, reveal previously unknown genes associated with photosynthetic and flagellar functions, and establish links between ciliopathy and the composition and function of flagella.

The Chlamydomonas genome reveals the evolution of key animal and plant functions

Chlamydomonas reinhardtii is a unicellular green alga whose lineage diverged from land plants over 1 billion years ago. It is a model system for studying chloroplast-based photosynthesis, as well as the structure, assembly, and function of eukaryotic flagella (cilia), which were inherited from the common ancestor of plants and animals, but lost in land plants. We sequenced the approximately 120-megabase nuclear genome of Chlamydomonas and performed comparative phylogenomic analyses, identifying genes encoding uncharacterized proteins that are likely associated with the function and biogenesis of chloroplasts or eukaryotic flagella. Analyses of the Chlamydomonas genome advance our understanding of the ancestral eukaryotic cell, reveal previously unknown genes associated with photosynthetic and flagellar functions, and establish links between ciliopathy and the composition and function of flagella.