Bio-crude transcriptomics: Gene discovery and metabolic network reconstruction for the biosynthesis of the terpenome of the hydrocarbon oil-producing green alga, Botryococcus braunii race B (Showa)*
© Molnár et al.; licensee BioMed Central Ltd. 2012
Received: 13 April 2012
Accepted: 19 October 2012
Published: 30 October 2012
Microalgae hold promise for yielding a biofuel feedstock that is sustainable, carbon-neutral, distributed, and only minimally disruptive for the production of food and feed by traditional agriculture. Amongst oleaginous eukaryotic algae, the B race of Botryococcus braunii is unique in that it produces large amounts of liquid hydrocarbons of terpenoid origin. These are comparable to fossil crude oil, and are sequestered outside the cells in a communal extracellular polymeric matrix material. Biosynthetic engineering of terpenoid bio-crude production requires identification of genes and reconstruction of metabolic pathways responsible for production of both hydrocarbons and other metabolites of the alga that compete for photosynthetic carbon and energy.
A de novo assembly of 1,334,609 next-generation pyrosequencing reads form the Showa strain of the B race of B. braunii yielded a transcriptomic database of 46,422 contigs with an average length of 756 bp. Contigs were annotated with pathway, ontology, and protein domain identifiers. Manual curation allowed the reconstruction of pathways that produce terpenoid liquid hydrocarbons from primary metabolites, and pathways that divert photosynthetic carbon into tetraterpenoid carotenoids, diterpenoids, and the prenyl chains of meroterpenoid quinones and chlorophyll. Inventories of machine-assembled contigs are also presented for reconstructed pathways for the biosynthesis of competing storage compounds including triacylglycerol and starch. Regeneration of S-adenosylmethionine, and the extracellular localization of the hydrocarbon oils by active transport and possibly autophagy are also investigated.
The construction of an annotated transcriptomic database, publicly available in a web-based data depository and annotation tool, provides a foundation for metabolic pathway and network reconstruction, and facilitates further omics studies in the absence of a genome sequence for the Showa strain of B. braunii, race B. Further, the transcriptome database empowers future biosynthetic engineering approaches for strain improvement and the transfer of desirable traits to heterologous hosts.
Excessive reliance on fossil hydrocarbons for world energy and synthetic chemistry needs has led to environmental degradation such as global warming, economic imbalances, and their associated national and geopolitical risks for producing and exporting nations alike. The problem is most acute in the liquid fuel sector (itself responsible for about two thirds of the global energy demand, ), where renewable sources of energy made the least impact up until now. Microbial conversion of feedstocks (simple sugars, starches and other polysaccharides, and total biomass) to biofuels (primarily ethanol as of now, but preferably other alcohols and fatty acid esters in the future) promises a renewable and potentially globally distributed source of transportation fuel. “Photosynthetic biofuels” directly link the biosynthesis of energy-rich, storable, transportation-friendly, and fuel infrastructure-compatible metabolites to photosynthesis, using land plants, eukaryotic microalgae or cyanobacteria as cell factories to yield a sustainable and potentially carbon-neutral source of fuels. Microalgae may be especially advantageous as they can be grown on marginal or even non-arable land, and may use water sources not directly utilizable in traditional agriculture [2–5]. In addition, microalgae grow faster, have higher photosynthetic productivity [1, 5], and accumulate biofuel feedstocks to a much higher percentage of their total biomass than land plants [6, 7]. These advantages translate to 15–60 times higher annual oil productivity as compared to soybean, the main US oil crop, grown on the same area of land [2, 4, 5].
Oleaginous algae accumulate storage carbon and energy in the form of neutral lipids, mainly triacylglycerols (TAG), which need to be extracted from intracellular oil bodies by disrupting the cells. TAG conversion to biofuel crude usually involves transesterification of the constituent fatty acids with alcohols before the resulting bio-crude is refined to transportation fuels . In contrast, the cosmopolitan green colonial microalga Botryococcus braunii (Chlorophyta, Trebouxiophyceae) stores photosynthetic carbon in the form of liquid hydrocarbons which need no chemical conversion to provide biofuel crude. This ready-made bio-crude is compatible with, and can be directly processed by the existing petrochemical refinery and distribution infrastructure to yield jet fuel, gasoline, and diesel with little coke formation . Further, this liquid hydrocarbon bio-crude is stored in an extracellular polymeric matrix material surrounding the individual cells of the colony . This may permit relatively mild extraction procedures that use nondisruptive, nontoxic solvents, thus allowing the recycling of “de-fatted” algal biomass (“milking” ).
B. braunii strains belong to three races defined by the chemical nature of their liquid hydrocarbon products. Race A strains accumulate C23-C33 alkadienes and alkatrienes derived from very long chain fatty acids (VLCFA) [10, 11]. Race B and race L strains produce triterpenoids (C30-C37 botryococcenes and methylated squalenes) or tetraterpenoids (C40 lycopadiene), respectively [10, 11]. B race strains of B. braunii typically accumulate hydrocarbons at 30-40% of their dry cell weight, although their hydrocarbon content can reach as high as 86% . Hydrocarbon accumulation apparently does not require specific metabolic triggers like nitrogen starvation as seen in TAG-accumulating microalgae [9–11]. Hydrocarbons originating from all races of B. braunii have repeatedly been isolated from fossil crude oils and coal deposits in significant amounts, indicating a possible role for this alga in the formation of our current petroleum reserves [9, 13] (and references therein).
The genomes of about a dozen photosynthetic algae, including model organisms such as Chlamydomonas reinhardtii, Chlorella variabilis and Volvox carteri have been sequenced and annotated (http://genome.jgi-psf.org) to gain insight into unique aspects of algal biology [14–16]. These investigations also led to the reconstructions of metabolic pathways and networks for selected biological processes of interest, including those that govern lipid production [17–19]. Multiple homologs of enzymes involved in triacylglycerol biosynthesis, such as diacylglcyerol acyltransferases, phospholipid diacylglycerol acyltransferase and phosphatidate phosphatases, have been identified in algae [19, 20]. Next-generation sequencing of transcripts from algal cultures grown under various conditions was used to discover regulated genes that are critical in key biochemical processes and biological pathways [21, 22]. Next-generation sequencing has also been utilized to assemble the complete transcriptome of the microalgae Dunaliella tertiolecta and Chlorella vulgaris, and to discover pathways involved in TAG biosynthesis [23, 24].
Functional annotation, the process of assigning biological meaning to genomic and transcriptomic data, is an important step in extracting useful information from genomic and transcriptomic projects. The annotation process involves assigning biological pathway, ontology, and protein domain data to genes and transcripts. In the case of de novo assembly of genomes and transcriptomes, sequence similarity is a common basis for assigning annotations . Sequences from primary annotation databases, such as KEGG , Gene Ontology , and Pfam , can be used to assign biological function to the assembled transcripts, and further manual curation can be used to validate and extend the annotations derived from sequence similarity. To provide the most utility from a collection of functional annotations, data mining environments are commonly used to facilitate the analysis of large-scale datasets and to provide a public repository of functional data. Several such tools are available for algae, including the Algal Functional Annotation Tool, which allows visualization of sets of transcripts on pathway maps and provides for statistically rigorous functional analysis of large sets of transcripts .
Results and discussion
Biomass sample collection and RNA isolation from the Showa strain of B. braunii
The advent of massively parallel high throughput cDNA sequencing techniques now allows the cost-effective de novo assembly and analysis of the transcriptomes of organisms with still unsequenced genomes [23, 24]. The genome size of B. braunii Showa has recently been measured at 166 Mbp, significantly larger than that of the largest sequenced Chlorophyta alga, V. carteri at 138 Mbp . This large genome, and the long and frequent repeat regions present in the B. braunii genome (A. Koppisch, personal communication) make the sequencing and analysis of this genome challenging. Although a transcriptome sequence database cannot capture untranscribed genomic regions (for example promoters) and does not reflect post-transcriptional regulation, it can still serve as a useful tool for gene discovery and metabolic pathway and network reconstruction, and will inform further proteomic and genomic analyses.
As a first step to construct such a database, we isolated total RNA from seven time points (days 0, 3, 5, 8, 14, 18, and 22) during the four week culture cycle of the Showa strain of B. braunii, race B. These time points are centered on the early portion of the culture cycle (days 0–8) because previous studies have shown that the biosynthesis and accumulation rate of botryococcenes is maximal in this period [35, 37, 38]. Additionally, both enzyme activity and gene expression associated with botryococcene biosynthesis have been shown to be maximal during these early time points [35, 37, 38]. Because of this, we hoped that the resulting transcriptomic database would be enriched for transcripts related to liquid hydrocarbon biosynthesis.
Total RNA from each time point was purified using TRIzol and LiCl precipitation to eliminate co-purifying polysaccharides. Each sample was analyzed for protein and polysaccharide contaminations ( Additional file 1: Table S1) and 5-μg aliquots from the samples with the most RNA (days 0, 3, and 5) were analyzed on a denaturing agarose gel ( Additional file 1: Figure S1A). To further confirm successful isolation of good quality RNA, RT-PCR was carried out to amplify cDNA fragments for several B. braunii Showa genes (Additional file 1: Figure S2), including squalene synthase (SS)  and squalene synthase-like-1 (SSL-1) . All of the RNA samples were then pooled into a single sample, treated with DNAse (Additional file 1: Figure S1B), and analyzed for any remaining RNase contamination (Additional file 1: Figure S1C). These analyses indicated that the isolated RNA was of high quality and did not contain active nucleases. The pooled RNA sample was submitted to the Department of Energy Joint Genome Institute (JGI) for transcriptome sequencing.
De novo assembly of the B. braunii Showa transcriptome
454 pyrosequencing yielded 1,334,609 reads (620 Mb of data) representing the ESTs from mRNA isolated from whole cells of a near-axenic B. braunii Showa culture  as described above. The reads have been deposited by the JGI into two publicly available Sequence Read Archive accessions, SRX028986 and SRX028987. We assembled these reads into 46,422 contigs with an average length of 756 bp (Additional file 1: Figure S3) using a multi-step, recursive sequence assembly protocol [39, 40] as described in the Methods. This new transcriptome assembly provides significantly improved coverage over that of Watanabe et al. for a different B race strain of B. braunii (27,427 contigs with an average length of 267 bp for strain BOT-22 ). Contig coverage in the Showa transcriptome assembly spans three orders of magnitude, from 1X to 8,231X. The most highly represented transcripts either encode proteins related to photosynthesis, or do not show significant similarity to sequences in the GenBank non-redundant database (E-value < 1e-5, Additional file 1: Table S2). To benchmark the quality of our assembly, we compared the B. braunii Showa transcriptome to the core set of 458 conserved proteins (CEGs) that occur in a wide variety of eukaryotes . Using the Core Eukaryotic Genes Mapping Approach (CEGMA) algorithm [41, 42], we recovered 451 out of the 458 core proteins (98.4%, E value cutoff ≤ 1e-5), with 325 out of the 458 CEGs (71.0%) yielding alignments whose lengths exceed 60% of either the CEG or the contig sequence.
Functional annotation, web-based annotation tool and data depository
Functional annotations were assigned to all unique transcript sequences using a previously described annotation pipeline . The Kyoto Encyclopedia of Genes and Genomes (KEGG) , MetaCyc , Reactome , Panther , and Pfam  annotation databases were chosen to provide biological pathway, ontology, and protein domain annotations. Orthologous proteins from C. reinhardtii (a model green alga) and Arabidopsis thaliana (thale cress) were also used to infer ontology identifiers from their respective Gene Ontology and MapMan Ontology sets  (Additional file 1: Table S3). Functional annotations were assigned to 20,906 transcripts (45%), while 8,575 sequences yielded significant BlastX hits (E value <1e-5) against the non-redundant protein database of GenBank. Species in the green plant lineage were the most frequent sources from which functional annotations were derived. Top-hit analysis of KEGG protein alignments producing pathway annotations confirmed the algal character of the transcripts (Additional file 1: Table S4). Eleven of the top 15 organisms are algae or land plants, with the two top organisms, V. carteri and C. reinhardtii, together contributing 25% of all KEGG pathway assignments. An analysis of the contigs by the Metagenomics RAST server  found top Blast hits to proteins from fungi (Ascomycota and Basidiomycota, 19.8% of the contigs), animals (Chordata, Arthropoda, and Nematoda, 10.7% of the contigs), and bacteria (Firmicutes, Proteobacteria, and Actinobacteria, 10.1% of the contigs). However, this analysis likely provides a highly inflated estimate of non-Botryococcal transcripts in our database. While some of these transcripts may indeed originate from contaminating organisms reflecting the non-axenic nature of the algal culture and/or sample handling mistakes introduced during the sequencing process, others still likely represent genuine B. braunii transcripts. These transcripts may highlight gaps and biases in the GenBank database, reflect localized high similarities to phylogenetically distant homologues, originate from sequencing/assembly mistakes, or represent genes from recent horizontal gene transfer.
In order to facilitate the exploration of the data and to expedite future B. braunii omics and functional genetic studies, the sequences and their annotations have been deposited into the Algal Functional Annotation Tool  and are publicly accessible from a web-based portal at http://pathways-pellegrini.mcdb.ucla.edu/botryo1. A combined view of annotations from all primary databases for any particular transcript can be accessed by the transcript ID, and provides pathway, ontology, and protein domain data alongside primary sequence data. Functional enrichment testing and dynamic pathway visualization may be performed for lists of transcripts from within the annotation tool. Transcripts may also be looked up by biological function using a keyword search tool that returns lists of transcripts with annotations matching a keyword or a phrase. Lastly, a pathway browser tool allows visualization of transcripts for any KEGG pathway of interest. The public repository of the B. braunii Showa transcriptome assembly and the utilities provided to query the data will furnish a platform to update annotations as these are made available by future studies and characterization of additional strains.
Manual curation and pathway reconstruction of the terpenome
Biosynthesis of terpenes yields a large variety of essential primary metabolites and specialized secondary metabolites in plants. Terpene biosynthesis provides membrane sterols, phytohormones, carrier molecules for N-glycan biosynthesis, pigments and antioxidants, volatile oils, aroma compounds and resins, and various toxins. Terpene biosynthesis also provides the prenyl side chains of photosynthetic pigments and meroterpenoid carrier molecules in oxidative phosphorylation, as well as polyprenyl compounds for the prenylation of proteins [49, 50]. Crucially, it also underpins the production of extracellular liquid hydrocarbons and contributes to the polymeric extracellular matrix materials in race B strains of B. braunii[11, 51]. C30-C37 triterpenes, in the form of methylated, oxidized, and cyclized botryococcenes, as well as methylated squalene, are initially synthesized inside the cells and at least botryococcenes can be found in intracellular oil bodies . However, the majority (95%) of the botryococcenes are deposited into the colony extracellular matrix . These oils that may be described as “bio-crude” are of prime interest for the biofuel industry as petroleum replacements.
To derive a comprehensive picture of the biosynthesis of the terpenome of B. braunii Showa, contigs in our databases that encode proteins for metabolic pathways yielding terpenes and their precursors (Figure 1) were collected and supplemented by further contigs identified by targeted Blast searches and Pfam keyword inquiries, as described in the Methods section. Contigs were extended and assembled into isotigs, and sequence discrepancies and potential frameshifts in the derived protein models were resolved using manual alignments and scrutiny of the assemblies. Manually curated transcript sequences (contigs and isotigs) are referred heretofore as “curated contigs”, while transcripts assembled without human supervision are simply called “contigs”. Contigs and curated contigs were considered as potentially originating from a different organism of the non-axenic culture if their deduced protein products show the highest similarities to fungal, animal or bacterial proteins with no or very low similarities to plant proteins, and if at the same time their codon usage significantly differs from that of B. braunii (http://www.kazusa.or.jp/codon/index.html). Subcellular localization of the deduced proteins was also predicted using TargetP (http://www.cbs.dtu.dk/services/TargetP/).
Based on databank similarities and the prediction of facile translational start and stop codons, full-size protein models were derived for approximately 30% of the manually curated transcripts. A distinguishing feature of B. braunii Showa transcripts is their surprisingly long (>1-2 kb) 3’ untranslated regions (UTRs), as already noted by Okada et al.. Machine-assembled contigs with moderate to high sequence coverage, apparently derived from these long 3’ UTRs, are abundant in the transcriptome database. These may be partially responsible for the relatively large number of contigs in the dataset, and are also in apparent congruence with the large genome size of B. braunii[36, 52].
Transcript abundances were approximated from sequence coverage, i.e. the number of primary sequence reads for a particular contig divided by the length of that contig. True estimation of expression levels would require the utilization of other techniques like qRT-PCR, DNA microarray analysis, or proteomics [18, 21, 22].
Biosynthesis of terpene precursors
In contrast, a complete contingent of deduced enzymes for the MEP/DOXP pathway [49, 55, 56] is well represented in the B. braunii Showa transcriptome. The MEP/DOXP pathway uses D-glyceraldehyde 3-phosphate and pyruvate as its metabolic input. Photosynthetic 3-phospho-D-glycerate (1) may be converted to D-glyceraldehyde 3-phosphate (3) via 3-phospho-D-glyceroyl phosphate (2) by the collective action of phosphoglycerate kinase (PGK, E.C. 188.8.131.52) and glyceraldehyde phosphate dehydrogenase (GAPDH, Figure 4, Additional file 1: Table S6). PGK is represented in our transcriptome by three nonredundant curated contigs for three presumed isoenzymes, one of which may have originated from a fungal cohabitant. Curated contigs for one predicted isoform of the NADP+-dependent GAPDH (E.C. 184.108.40.206, Benson-Calvin cycle) and five inferred isoforms (three of them probably of fungal origin) of the NAD+-dependent GAPDH (E.C. 220.127.116.11, glycolysis and gluconeogenesis) were also identified. During glycolysis, D-glyceraldehyde 3-phosphate (3) may also be converted to (1) by the non-phosphorylating NADP+-dependent glyceraldehyde 3-phosphate dehydrogenase (GAPN, E.C. 18.104.22.168), encoded by a curated contig with moderate coverage (25–99 reads/kb) in the dataset (Additional file 1: Table S6). TargetP predictions for the subcellular localization of these enzymes provided support for the mitochondrial targeting of the putative NAD+-dependent GAPDH (E.C. 22.214.171.124) isoforms, and cytosolic localization for the PGK isozymes. Some transcripts for inferred PGK and NAD+-dependent GAPDH (E.C. 126.96.36.199) are present at very high abundance (>250 reads/kb) in the transcriptome, suggesting a high flux involving D-glyceraldehyde 3-phosphate in the glycolysis/gluconeogenesis pathways in B. braunii Showa. Contigs encoding short fragments of PGK and GAPDH orthologs have recently been described from the race B B. braunii strain BOT-22 . These show 89-99% identity at the amino acid level with the corresponding enzymes from B. braunii Showa predicted in this study.
Pyruvate (6) is formed from (1) in three steps by phosphoglycerate mutase (PGAM, E.C. 188.8.131.52), phosphopyruvate hydratase (ENO, E.C. 184.108.40.206) and pyruvate kinase (PK, E.C. 220.127.116.11) during glycolysis/gluconeogenesis (Figure 4). All these enzymes are encoded in the transcriptome as multiple deduced isozymes (Additional file 1: Table S6). Curated contigs 10955, 16949, 23205, and 41366 show the highest similarities to individual domains of PKs with four similar catalytic domains each, thus probably representing a single multifunctional enzyme. All these curated contigs have moderate coverage. Curated contig 43373, encoding a predicted cytoplasmic ENO, is the only exception and is represented by very abundant ESTs. Two contigs from B. braunii BOT-22 (FX085139 and FX085140) encoding short regions of PK  show 94% amino acid identity to two inferred PK isozymes encoded in our dataset (curated contig 10955 and 41736, respectively).
Four curated contigs code for three predicted isozymes of the first enzyme of the MEP/DOXP pathway, the thiamine diphosphate-dependent 1-deoxy-D-xylulose 5-phosphate synthase (DXS, E.C. 18.104.22.168) that produces 1-deoxy-D-xylulose 5-phosphate (7) from (3) and (6) (Figure 4, Additional file 1: Table S6). Multiple isoenzymes of DXS are routinely found in land plants , and clade into three phylogenetically distinct families . Constitutively expressed DXS isozymes of these plants produce precursors for essential terpenoids, while certain inducible DXS isozymes specialize in stress response and ecological interactions with symbionts or pathogens . In contrast, genomic evidence shows that strains of green algae harbour only a single DXS each. These proteins form a sister clade to the three DXS clades of land plants [35, 49]. The Okada group has recently reported the cloning and biochemical characterization of three isozymes of DXS from B. braunii Showa . Interestingly, these three isoenzymes all fall into the basal clade for Chlorophyta DXSs, thus representing paralogous sequences resulting from gene duplications within the green algal lineage . All three isozymes were found to be expressed simultaneously, and show similar kinetic parameters except for a higher temperature tolerance for DXS-III . Our data also indicate similar, moderate transcript abundances for DXS-I (curated contig 10163, 66.6 reads/kb) and DXS-II (curated contig 42027, 93.7 reads/kb), with a somewhat lower abundance for DXS-III-related ESTs (curated contigs 07667 and 11032, 4.0 reads/kb and 29.4 reads/kb, respectively). Short curated contigs representing orthologs of DXS-II (FX085276 and FX085277) and DXS-III (FX085274 and FX085275) were also identified in the preliminary transcriptomic analysis of the B. braunii race B strain BOT-22 , with amino acid identities to the Showa enzymes in the 72-84% range ( Additional file 1: Table S6). All three DXS isozymes of the Showa strain have been found to contain chloroplast targeting sequences at their N-termini , in agreement with our TargetP predictions ( Additional file 1: Table S6). DXS has been described as one of the rate-limiting steps of the MEP/DOXP pathway in plants [49, 57], thus the expression of three isoforms of this enzyme in B. braunii Showa might provide an increased metabolic flux for the production of terpenoid precursors. Alternatively, each DXS isoform might be associated with the production of a specific class of terpenoids. Considering the absence of the cytoplasmic mevalonate pathway as an alternative to furnish isoprene precursors, evolutionary optimization of the DXS step by repeated gene duplications providing parallel capacity for the production of (7) might have been beneficial for B. braunii race B strains.
1-deoxy-D-xylulose 5-phosphate (7) is converted to 2-C-methyl-D-erythritol 4-phosphate (8, Figure 4) by the NADPH-dependent enzyme 1-deoxy-D-xylulose-5-phosphate reductoisomerase (DXR, E.C. 22.214.171.1247). A single deduced isozyme of DXR is encoded in our dataset by a curated contig with moderate coverage (Additional file 1: Table S6). The DXR-catalyzed reaction is the first committed step towards the production of isoprene precursors since (7) is also utilized for the biosynthesis of thiamine diphosphate and pyridoxal phosphate . The rest of the pathway involves the CTP-dependent conversion of (8) to 4-diphosphocytidyl-2-C-methyl-D-erythritol (9) by IspD (2-C-methyl-D-erythritol 4-phosphate cytidylyltransferase, E.C. 126.96.36.199), and the ATP-dependent phosphorylation of (9) to 4-diphosphocytidyl-2-C-methyl-D-erythritol 2-phosphate (10) by IspE (4-(cytidine 5’-diphospho)-2-C-methyl-D-erythritol kinase, E.C. 188.8.131.52). Curated contigs for these two inferred enzymes are present in the Showa transcriptome at moderate coverage. Formation of 2-C-methyl-D-erythritol 2,4-cyclodiphosphate (11) from (10) by IspF (2-C-methyl-D-erythritol 2,4-cyclodiphosphate synthase, E.C. 184.108.40.206) releases CMP, followed by the two-electron reduction of (11) to (E)-4-hydroxy-3-methylbut-2-en-1-yl diphosphate (12) by the [4Fe-4S] enzyme IspG (4-hydroxy-3-methylbut-2-enyl diphosphate synthase, E.C. 220.127.116.11). The last step of the MEP/DOXP pathway is the formation of both isopentenyl diphosphate (IPP, 13) and dimethylallyl diphosphate (DMAPP, 14) from (12) by another [4Fe-4S] enzyme, IspH (also known as LytB, 4-hydroxy-3-methylbut-2-enyl diphosphate reductase, E.C. 18.104.22.168). Both IspG and IspH accept electrons from a ferredoxin. These electrons may originate directly from the photo-oxidation of water during photosynthetic conditions in the chloroplast without the involvement of reducing cofactors, while a ferredoxin reductase is required in the dark to channel electrons from cellular pools of NADPH . The branching reaction catalysed by IspH is in stark contrast to the MVA pathway that yields IPP (13) exclusively , which has to be later isomerised to DMAPP (14) by Idi (isopentenyl-diphosphate delta-isomerase, E.C. 22.214.171.124). Each of the predicted enzymes for the downstream half of the MEP/DOXP pathway (IspF and onwards) are encoded by single nonredundant curated contigs with high to very high sequence coverage in the B. braunii Showa transcriptome (200.7, 242.9, and 426.6 for IspF, IspG and IspH, respectively), indicating vigorous transcription and perhaps robust metabolic flow through these enzymes. This is in contrast to some plant systems where the same enzymes were found to be rate-limiting [49, 57]. Two presumed isozymes of Idi are encoded by two nonredundant curated contigs with only low coverage in the Showa transcriptome. Another curated contig with moderate coverage for a putative fungal IPP isomerase was also identified (34876, Additional file 1: Table S6). Although Idi is generally present in organisms utilizing the MEP/DOXP pathway for terpenoid precursor biosynthesis, it has not been found strictly essential and plays only a supplementary role in optimizing IPP and DMAPP ratios . Type II IPP isomerases, detected in Streptomyces spp. and in Synechocystis spp. [60, 61], were not found in the Showa transcriptome. TargetP predictions for all MEP/DOXP enzymes suggest chloroplast targeting, with the exception of IspH where both chloroplast and mitochondrial targeting seems equally plausible (Additional file 1: Table S6).
The cyanobacterium Synechocystis sp. PCC6803 contains no MVA pathway but features a full complement of the MEP/DOXP pathway. Nevertheless, neither pyruvate (6) nor 1-deoxy-D-xylulose 5-phosphate (7) has been observed to stimulate IPP biosynthesis in cell extracts. Similarly, fosmidomycin has not been seen to inhibit DXR by reducing IPP biosynthesis in vitro or by reducing cellular growth in vivo. On the other hand, metabolites of the reductive pentose phosphate cycle, especially D-xylulose 5-phosphate, increased IPP formation even in the presence of the DXR inhibitor fosmidomycin . While the immediate entry point of the pentose phosphate cycle metabolites into the MEP/DOXP pathway is currently not known, it is assumed to be downstream of DXR (Figure 4). A survey of the Showa transcriptome identified curated contigs encoding putative chloroplast-targeted RPE (ribulose phosphate 3-epimerase, E.C. 126.96.36.199) and thiamine diphosphate-dependent TKTL (transketolase, E.C. 188.8.131.52) enzymes that yield D-xylulose 5-phosphate in the pentose phosphate cycle (Additional file 1: Table S6). Short contigs that encode fragments of the TKTL enzyme have recently been identified from B. braunii BOT-22 (FX085315 and FX085314) , and show 87-90% amino acid identities with the enzyme encoded by curated contig 32329 of the Showa transcriptome. Transcripts of putative fungal origin were also detected for both RPE and TKTL. Thiamine diphosphate-dependent phosphoketolase (XFP, E.C. 184.108.40.206) was represented only by a single curated contig of potential fungal origin. This enzyme generates (3) and acetyl phosphate from D-xylulose 5-phosphate and inorganic phosphate. Curated contig 30447 for RPE has moderate sequence coverage at 55.8 reads/kb, while reads for curated contig 32329 for TKTL are highly abundant at 436.6 reads/kb.
While the importance, or even the dominance, of the pentose phosphate cycle for isoprene precursor biosynthesis has been speculated upon for the BOT-22 and BOT-70 strains of B. braunii race B [31, 32], these recent studies employed limited transcriptomic and EST datasets. Our data support the presence of a fully functional MEP/DOXP pathway in B. braunii Showa, with multiple paralogous DXSs with moderate EST coverage for each, providing a reasonable entry for isoprene precursor biosynthesis. The sequence coverage of the inferred downstream half of the MEP/DOXP pathway enzymes (IspF and onwards) is higher than that of the upstream half of the pathway, which would be consistent with an anaplerotic feed of metabolites from the pentose phosphate cycle. Future studies of the dynamics of transcription of the MEP/DOXP pathway and the pentose phosphate cycle enzymes by e.g. qRT-PCR, in relation to the age of the culture and hydrocarbon production levels would be necessary to shed more light on this issue. Investigating the inhibition by fosmidomycin of B. braunii cultures in vivo and/or the isolated Showa DXR enzyme in vitro may also provide further clues to map the flux of metabolites through these alternative biosynthetic routes.
Terpene backbone biosynthesis
C10 geranyl diphosphate (15, Figure 5) is synthesized by geranyl diphosphate synthase (GDPS, E.C. 220.127.116.11) from (13) and (14). 15 yields monoterpenes in plants, and serves as a precursor for solanesyl diphosphate synthase (SDPS, E.C. 18.104.22.168) that affords all-trans-nonaprenyl diphosphate (18) for the polyprenyl chains of terpenoid quinones in the mitochondria (ubiquinones) and in the chloroplast (plastoquinones). Both GDPS and SDPS are predicted to be present in the Showa transcriptome as a single nonredundant curated contig each with low sequence coverage (Additional file 1: Table S7).
Farnesyl diphosphate synthase (FDPS, E.C. 22.214.171.124) generates C15 farnesyl diphosphate (16) in two reaction steps via 15. 16 is the precursor for the sesquiterpenes and the triterpenoids (including squalene and phytosterols), and provides the farnesyl side chains for post-translational modification of proteins. 16 is also the substrate for decaprenyl diphosphate synthase (PDSS1, E.C. 126.96.36.199). This enzyme produces all-trans-decaprenyl diphosphate (19) for the side chain of the mitochondrial electron carrier ubiquinone-10 (coenzyme Q10). The cis- prenyltransferase dehydrodolichyl diphosphate synthase (DHDDS, E.C. 2.5.1.-) also uses 16 as its substrate to generate dehydrodolichyl diphosphates (di-trans, poly-cis-polyprenyl diphosphate, 20) that serve as precursors to the glycosyl carrier lipid dolichol for N-glycan biosynthesis. Crucially, 16 is also the precursor for the liquid hydrocarbon triterpenoid botryococcenes and methylated squalenes, and the cell wall ether lipids and some of the matrix polymers of B. braunii Showa [9, 11]. Two nonredundant curated contigs for two putative isozymes of FDPS with 72% amino acid identity were identified in the Showa transcriptome, both with moderate sequence coverage (Additional file 1: Table S7). The FDPS isozyme encoded by curated contig 15137 is predicted to be localized outside the chloroplast, the mitochondrion, or the secretory apparatus of the cell. A third curated contig for a presumed FDPS with extremely low sequence coverage may stem from a fungal source (Additional file 1: Table S7). Deduced DHDDS isozymes are encoded by four curated contigs of low sequence coverage, one of which is potentially of fungal origin. A single curated contig with low sequence coverage encodes a putative PDSS1 that may also have been derived from a fungal cohabitant (Additional file 1: Table S7).
Geranylgeranyl diphosphate (17) is produced by geranylgeranyl diphosphate synthase (GGDPS, E.C. 188.8.131.52), encoded in the Showa transcriptome by only a single nonredundant curated contig with moderate sequence coverage (Additional file 1: Table S7). In higher plants, this enzyme may be present as three isoenzymes - one cytosolic, one plastidic and one mitochondrial . TargetP prediction supports mitochondrial localization for the single GGDPS identified in the Showa transcriptome. The expected plastidic GGDPS might be encoded by the chloroplast genome and thus likely missed by our transcriptome database. Alternatively, the other isoforms may arise from multiple targeting, or simply represent gaps in the current transcriptome database due to lower levels of expression or sequencing/assembly artifacts. Plant GGDPS enzymes may initiate synthesis of 17 from 14 (Figure 5), but the most effective substrate of GGDPSs from animals and fungi is 16[49, 65]. 17 serves as the precursor for the C20 diterpenes including gibberellins, the phytyl side chains of chlorophyll, phylloquinone and tocopherol (vitamin E), and the geranylgeranyl chains of prenylated proteins. 17 is also the precursor for tetraterpenoid carotenoids, including lycopene, lutein, canthaxanthin and others isolated from race L and race B B. braunii strains [11, 66].
Liquid triterpenoid hydrocarbon biosynthesis
In addition to SQSs and phytoene synthases (CrtB, see next section), genome sequence surveys of algae often identify a gene encoding a putative, uncharacterized protein with the head-to-head trans-isoprenyl diphosphate synthase fold. These predicted proteins (Class 1 isoprenoid biosynthesis-related proteins [ISR]), form a clade that is distinct from SQSs and CrtB enzymes. ISRs had been hypothesized earlier to constitute the then-unknown botryococcene synthase , but the recent in vitro studies from the Chappell laboratory, as described above, conclusively showed that these proteins are not necessary for botryococcene biosynthesis . A single nonredundant curated contig with low sequence coverage represents this cryptic ISR enzyme in the Showa transcriptome ( Additional file 1: Table S8).
Hydroxylation, epoxidation, and formation of O-containing heterocycles (for example compounds 31, 32, and 33, Figure 7) increase the complexity of race B liquid and polymer hydrocarbons. These oxidized methylsqualenes, together with carotenoids and very long chain fatty acids (and to a smaller degree oxidized botryococcenes) also support the biosynthesis of ether lipids like braunixanthin 1 (34) and matrix polymers (35) as precursors [9–11, 51, 77]. Candidates for the introduction of the oxygen functionality into linear triterpene structures might be enzymes similar to the flavoprotein squalene monooxygenase (SQLE, E.C. 184.108.40.206) that produces (3S)-2,3-epoxy-2,3-dihydrosqualene from 22 by incorporation of molecular oxygen. Enzymes similar to SQLE might also accept tetramethylsqualene 30 as their substrate and catalyse epoxide formation towards the centre of the long terpene chain, yielding compounds like diepoxy-tetramethylsqualene (32). Reducing equivalents are channelled to these enzymes by NADPH-dependent hemoprotein reductases. The Showa transcriptome contains seven nonredundant curated contigs with very low to low sequence coverage (3.4 to 20.2 reads/kb, Additional file 1: Table S8) encoding SQLE-like enzymes: these putative enzymes are candidates to channel squalenes and maybe botryococcenes towards the production of extracellular matrix materials.
Biosynthesis of other terpenoids
The terpenome of B. braunii Showa includes meroterpenoid quinones, the side chain of chlorophyll, diterpenoid gibberellins, triterpenoid phytosterols, tetraterpenoid carotenoids, polyprenyl carrier molecules (dolichol, see above), and the prenyl chains of proteins. All these primary and secondary metabolites draw on the common isoprene (IPP and DMAPP) pool generated by the chloroplast-based MEP/DOXP pathway, and divert these precursors away from the production of liquid and matrix hydrocarbons. We have curated and catalogued contigs representing putative enzymes involved in these competing pathways.
The prominent liquid hydrocarbon trans,trans- lycopadiene (46, Figure 9) of race L strains of B. braunii is probably biosynthesized by reduction of an acyclic carotenoid such as 15-cis- phytoene (25) or lycopene (42) . Lycopadiene (46) is considered a biomarker for race L strains of B. braunii as there are no reports for the occurrence of this compound in race B strains . Indeed, we could not identify plausible tetraterpene reductases in the Showa transcriptome for the biosynthesis of 46.
The diterpenoid growth hormones gibberellic acids  are derived from geranylgeranyl diphosphate (17) by cyclization followed by multiple oxidations in land plants. There are no unequivocal reports on the production of gibberellins in green algae , nor have genes with high similarity to key gibberellin production pathways been located in genome sequences . Fittingly, the first enzyme of the gibberellic acid pathway that generally harbours both ent- copalyl diphosphate synthase (E.C. 220.127.116.11) and ent- kaurene synthase (E.C. 18.104.22.168) activities could not be identified in the B. braunii Showa transcriptome. On the other hand, curated contigs that may encode enzymes for the oxidative processing of the tetracyclic intermediate ent-kaurene to gibberellin A4 are present in the transcriptome at low sequence coverage (Additional file 1: Table S12).
In addition to the biosynthesis of many primary metabolites of terpenoid or other origin, the production of large amounts of higher botryococcenes and methylated squalenes (Figure 7) in B. braunii Showa requires a robust supply of S-adenosylmethionine (SAM) that is used as a donor for methylation reactions. Transfer of the methyl group of SAM to the substrate by these methyltransferases yields S-adenosylhomocysteine that is hydrolysed by S-adenosyl-L-homocysteine hydrolase (AhcY, E.C. 22.214.171.124) to homocysteine and adenosine (Additional file 1: Table S13). Homocysteine is methylated by 5-homocysteine S-methyltransferases, including the S-methylmethionine-dependent MmuM (E.C. 126.96.36.199), the 5-methyltetrahydrofolate- and cobalamin-dependent MetH (E.C. 188.8.131.52), and the 5-methyltetrahydropteroyl-triglutamate-utilizing but cobalamin-independent MetE (E.C. 184.108.40.206), yielding L-methionine. Finally, S-adenosylmethionine synthase (MetK, E.C. 220.127.116.11) transfers the adenosyl moiety of ATP to methionine, yielding SAM and releasing phosphate and pyrophosphate. Machine-assembled contigs for all these deduced enzymes have been identified in the Showa transcriptome, with low (MetK, MmuM), low to moderate (MetH) and high to extremely high sequence coverage (AhcY, MetE, Additional file 1: Table S13).
Competing storage compounds: biosynthesis of triacylglycerols
Photosynthetic carbon and energy intended for storage is partitioned in B. braunii Showa amongst terpenoid hydrocarbons, triacylglycerols (TAGs), and carbohydrates. We have generated inventories for machine-assembled contigs predicted to encode crucial enzymes in hydrocarbon-competing storage compound biosynthetic pathways.
Long chain fatty acids may undergo further desaturations and chain elongations to afford long chain and very long chain polyunsaturated fatty acids (LC-PUFA and VLC-PUFA, respectively, Figure 11). Oxygen-dependent Ω6 (Δ12) fatty acid desaturases (FAD6, E.C. 1.14.19.-) and Ω3 (Δ15) fatty acid desaturases (FAD8, E.C. 1.14.19.-) in the chloroplast and the endoplasmic reticulum generate linoleic acid (C18:2 n-6) and α-linolenic acid (C18:3 n-3) and their longer chain equivalents, respectively, using acyl-glycerolipid substrates [23, 99, 100]. Both of these predicted enzymes are encoded in the Showa transcriptome at low to moderate sequence coverage. An apparent front-end desaturase, FADS2 (Δ6 fatty acid desaturase, E.C. 1.14.19.-) that may yield stearidonic acid (C18:4 n-3) or its longer chain equivalents, perhaps by using acyl-CoA esters , is also present in the transcriptome with low sequence coverage. For VLCFA and VLC-PUFA, recursive cycles of 2-carbon additions from malonyl-CoA and the following three desaturation steps outside of the chloroplast (probably in the microsomes) parallel that of de novo long-chain fatty acid biosynthesis in the chloroplast, but without the involvement of acyl carrier proteins . Thus, multiple contigs with moderate to high sequence coverage that encode putative very long chain fatty acid elongase (ELOVL5, E.C. 2.3.1.-), β-ketoacyl-CoA reductase (KAR, E.C. 1.1.1.-), 3-hydroxyacyl-CoA dehydratase (PHS1, E.C. 4.2.1.-), and enoyl-CoA reductase (TER, E.C. 1.3.1.-) enzymes have been found in the Showa transcriptome. However, a very long chain fatty acyl-CoA hydrolase (E.C. 3.1.2-) that may act as a thioesterase has not been identified. VLCFA, VLC-PUFA and hydrocarbon pyrolysis products of these have been observed in B. braunii, albeit primarily in the A race strains [80, 101–103].
Triacylglycerol storage lipids are assembled in the endoplasmic reticulum by two sequential acylations of sn-glycerol-3-phosphate, followed by dephosphorylation and a final acyl transfer (Figure 11). Exchange of acyl chains amongst TAG, glycerolipids and the acyl-CoA pool (acyl editing) provides flexibility to channel carbons for storage or for functional lipid biosynthesis . Following the ATP-dependent phosphorylation of glycerol by glycerol kinase (GlpK, E.C. 18.104.22.168), glycerol-3-phosphate O-acyltransferase (GPAT, E.C. 22.214.171.124) and lysophosphatidic acid acyltransferase (LPAAT, E.C. 126.96.36.199) generate phosphatidic acid using the acyl-CoA pool (Additional file 1: Table S14). Phosphatidic acid phosphatase (PAP, E.C. 188.8.131.52) affords sn- 1,2-diacylglycerol, to be acylated by diacylglycerol acyltransferase (DGAT, E.C. 184.108.40.206) to produce TAG. Diacylglycerols are also the precursors for the various polar lipids (glycosyl-glycerolipids and phosphoglycerolipids), while phospholipid:diacylglycerol acyltransferase (PDAT, E.C. 220.127.116.11) shuttles acyl groups amongst phosphoglycerolipids, betaine lipids, and TAG in an acyl-CoA independent process [19, 97, 103]. All these key TAG biosynthetic enzymes are predicted to be present in the Showa transcriptome with low sequence coverage ( Additional file 1: Table S14) which may be reflective the low level of TAGs predicted to be present in the B race of B. braunii.
In land plants, TAGs accumulate in oil bodies whose lipid/water interface features the structural proteins oleosins. Hydrophobic proteins (MLDP, major lipid droplet protein) that may be functionally equivalent to, but not structurally similar to plant oleosins have recently been identified by proteomic approaches in the oil bodies of Chlorophyta algae [20, 104]. Multiple contigs encoding presumed MLDPs are also featured in the Showa transcriptome, some with very high sequence coverage (contigs 0772, 35177 and 42893, with 237–295 reads/kb) that may indicate active transcription of the corresponding genes.
Competing storage compounds: biosynthesis of starch and other carbohydrates
A major sink for photosynthetic carbon intended for storage in algae are polysaccharides including starch. The biosynthesis of these compounds competes with those of hydrocarbon oils and TAG lipids in B. braunii, thereby reducing biofuel yield. On the other hand, starch and cellulosic biomass, after hydrolysis, may be utilized as a feedstock for the fermentative production of biofuel using non-photosynthetic microorganisms [105, 106].
Starch biosynthesis in the chloroplast is initiated by the phosphorylation of α-D-glucose at the C6 position by hexokinase (HK, E.C. 18.104.22.168) and glucokinase (Glk, E.C. 22.214.171.124) in ATP-dependent reactions. D-glucose 6-phosphate is then converted by phosphoglucomutase (PGM, E.C. 126.96.36.199) to α-D-glucose 1-phosphate that serves as a substrate for the ATP-dependent glucose-1-phosphate adenylyltransferase (GlgC, E.C. 188.8.131.52), a major rate-controlling enzyme of the pathway in plants and bacteria . The resulting ADP-glucose is then polymerized by starch synthase (GlgA, E.C. 184.108.40.206) to generate amylose with linear α-1,4-glycosidic linkages. Branched α-1,6 glycosidic linkages between α-1,4-glucan chains is generated by the 1,4-α-glucan branching enzyme (GlgB, E.C. 220.127.116.11) to yield water-insoluble amylopectin. Conversely, amylo-α-1,6-glucosidase (AGL, E.C. 18.104.22.168) hydrolyses α-1,6 glycosidic linkages to limit branching and to mobilise glucose from starch. Multiple contigs with low to moderate sequence coverage encoding these predicted enzymes are present in the Showa transcriptome, with the deduced glycogen debranching enzyme AGL encoded by a single contig that might have originated from a fungal cohabitant (Additional file 1: Table S15).
Considering that cellulose is a major constituent of the cell wall in green algae (up to 80% in Chlorella sp., ), the biosynthesis of this 1,4-β-D-glucan presents a major demand from the available photosynthetic carbon. Since Chlorophyta cell walls do not contain lignin [109, 110], cellulose from the biomass of these algae may provide a relatively more easily accessible source of sugar for the production of biofuels by fermentation. For the biosynthesis of cellulose, α-D-glucose 1-phosphate is activated by uridylation by UTP:α-D-glucose-1-phosphate uridylyltransferase (UGP, E.C. 22.214.171.124). The resulting UDP-glucose is the substrate for cellulose synthase (BcsA, E.C. 126.96.36.199). Orthologous genes are present in single copies in green algae . Contigs with low to moderate sequence coverage have been identified for these two inferred enzymes in the B. braunii Showa transcriptome (Additional file 1: Table S15).
Localization of liquid terpenoid compounds into the extracellular matrix
Apart from the biosynthesis of hydrocarbons with fuel characteristics and utility resembling that of fossil crude oil, B. braunii is also remarkable in depositing these compounds into a communal extracellular matrix that holds the colony together. Forming a fibrous network, this matrix consists of highly complex cross-linked polymers, originating from methylated squalenes and/or very long chain fatty acids . The hydrocarbons initially accumulate within the cells as oil bodies, but the large majority (95%) of the extractable liquid oils are found in the matrix . Traffic of these hydrocarbons seems to coincide with maturation of botryococcenes to C34 and higher homologues in B. braunii, race B [13, 75]. While the exact mechanism of the excretion of hydrocarbons into the extracellular space is unknown, the characterization and subsequent engineering of this trait into other (micro)organisms holds great promise for relieving product toxicity and simplifying extraction of advanced biofuels from biomass.
Efflux pumps  are one of the obvious candidates for the cellular export system of hydrocarbons. The Showa transcriptome contains numerous contigs encoding potential ABC (ATP-binding cassette) transporters (also known as multidrug-resistance-related proteins, MRPs) that mediate ATP-dependent transport of a bewildering array of molecules across organellar and cellular membranes (Additional file 1: Table S16). Plant ABC transporters fall into several subfamilies with varied (and often uncharacterized) functions but with a substantial functional redundancy of the family members within a single organism . Subfamily A members take part in the transport of various lipids including sterols and lipoproteins across membranes in animals. Although members of this subfamily have been identified in plants, their functions have not been characterized. Subfamily B members are multidrug resistance factors that take part in auxin, secondary metabolite, and xenobiotic traffic in plants, with mitochondrion-located members involved in iron-sulphur cluster trafficking. Subfamily C members play a role in detoxification and in vacuolar transport of glucuronides, chlorophyll degradation products and anthocyanins. Subfamily D members are essential for the import of VLCFA into the peroxisomes for β-oxidation. Subfamily G transporters are involved in the export of alkanes and other lipids that form the waxy cuticle and in the secretion of volatile compounds in flowers and roots of plants. Other group members convey resistance to various substances including terpenoids and herbicides [113, 114]. Subfamily G is generally expanded in plants compared to animals , with a large number of transcripts encoding presumed members of this subfamily also present in the Showa transcriptome. The proteins represented by these contigs, and to a lesser extent those of putative Subfamily D and A members (Additional file 1: Table S16), are prime candidates for the role of a liquid hydrocarbon exporter in B. braunii race B strains including Showa.
Conclusions and outlook
Functional annotation of the transcriptome of the B. braunii race B, strain Showa uncovered diverse biological processes operational in this hydrocarbon-producing green alga. Global comparisons with other algal transcriptomes reveal similar degrees of annotation coverage, with the number and the range of annotations supporting a strong degree of conservation between algal genomes and transcriptomes. The level of annotation coverage in B. braunii is supported by the identification of the majority of predicted enzymes within diverse biological pathways, for example within those related to the production of terpenes and lipids.
This study describes reconstituted metabolic pathways related to the biosynthesis of terpenoid hydrocarbons in the non-model green microalga B. braunii Showa, following the fate of photosynthetic carbon from 3-phosphoglycerate to the general terpenoid precursors IPP and DMAPP, the production of linear polyprenyl backbones, the biosynthesis of triterpenoids, the decoration/tailoring of botryococcene and squalene to yield liquid hydrocarbon compounds and matrix structural materials, and possible routes for the extracellular localization of these compounds. A recurrent theme is the expansion of particular gene families. This allows the adaptation of the paralogs to structurally orthogonal substrates (botryococcene methyltransferases), and permits neofunctionalization to support novel biochemical reactions (botryococcene synthase). Paralogs may yield increased metabolic flux, or may provide additional flexibility in terms of regulation, compertmentalization, and biochemical properties (deoxyxylulose phosphate synthase, and possibly ABC transporters for hydrocarbon export). Metabolic pathways leading to other terpenoids have also been reconstructed, and anabolic pathways for competing storage compounds (TAG and polysaccharides) were similarly mapped.
The reconstructed metabolic networks, their participating enzymes and the corresponding cDNA sequences provide a genetic and metabolic framework that should empower biosynthetic engineering approaches targeting the increased production of hydrocarbons in B. braunii, or the mobilization of these pathways into genetically tractable photosynthetic (algal or land plant) hosts or heterotrophic microbial strains [50, 120–122]. In particular, increasing the flux of photosynthetic carbon towards terpene precursors in the chloroplast or the cytosol of B. braunii by transplanting the cytoplasmic mevalonate pathway, or by fine tuning the expression of DXR (deoxyxylulose phosphate reductase) and Idi (isopentenyl-diphosphate delta-isomerase) may be interesting approaches [60, 120, 123]. Further experiments should investigate the importance of the pentose phosphate cycle as an anaplerotic pathway feeding into the MEP/DOXP pathway, and clarify the actual metabolite(s) and enzyme(s) connecting these two routes . Conversely, tuning down, or blocking the biosynthesis of competing storage carbons (like that of starch in Chlamydomonas[20, 124, 125]) may also increase the accumulation of liquid hydrocarbons. Varying the level and timing of the expression for the two-component botryococcene synthase [34, 38] and various storage hydrocarbon decorating enzymes (including the recently described methyltransferases for C32 triterpenoids , the still unknown methyltransferases for C34 triterpenoids, and the various oxidases that yield oxidized (methyl)squalenes, ether polymers and matrix materials) should help to tailor liquid hydrocarbon biosynthesis for particular biofuel applications. Investigations into the export of hydrocarbons into the extracellular matrix of B. braunii should uncover a particularly valuable trait that may have a widespread application to increase yield, reduce end product toxicity and simplify product recovery during biofuel manufacture by fermentation. Interesting questions are presented by the interconnection of the biosyntheses of terpenoid hydrocarbons and very long chain fatty acids with the formation of the extracellular matrix materials in B. braunii that provide the building blocks for an extracellular carbon storage organ, but also the physical basis for colony organization in this organism. Further transcriptomic, proteomic, metabolomic, and metabolic flux analyses that compare varied growth conditions influencing hydrocarbon accumulation would shed light on the regulatory networks and pathway interactions channeling carbon and energy flow in algal cells.
In addition to biofuel biosynthetic pathway discovery, the integrated data-mining environment offered by the publicly available web annotation tool described here allows researchers to query the B. braunii transcriptome for any specific transcript sets to rapidly and efficiently extract biologically relevant information related to different contexts. As the transcriptome dataset is revised and supplemented with addtional manually curated annotations, the most current functional data will be made available to the community via the public portal and annotation tool.
Algal cell cultures
A near axenic culture of B. braunii B race, Showa strain (a.k.a the Berkeley strain)  was obtained as previously described . Briefly, cultures were diluted with sterilized media, individual colonies were isolated under a microscope, and transferred to fresh, sterile media for growth. The resulting cultures were grown in modified Chu 13 media  at 22.5°C under a 12:12 light:dark (L:D) cycle using 13 W compact fluorescent 65 K light bulbs at a distance of 7.62 cm, which produced a light intensity of 280 μmol photons · m-2 · s-1. The cultures were continuously aerated with filter-sterilized, enriched air containing 2.5% CO2. Fifty milliliters of culture was used to inoculate 750 ml of subsequent subcultures every 4 weeks. Bacterial contamination of the cultures was routinely monitored by inoculating LB plates with 100 μl of undiluted media from a B. braunii culture that has had the algae removed by filtration (see below), and incubating the plates at 37°C for 4 days. This routinely yields an average of 5 bacterial colonies.
Sample collection and RNA isolation
B. braunii samples for total RNA isolation were collected at seven time points (days 0, 3, 5, 8, 14, 18, and 22) representing the four-week culture cycle. For days 0, 3, and 5, cells were collected from a full culture flask (750 ml) each, yielding a biomass of approximately 0.5 g (days 0 and 3) and 1 g (day 5) of wet cell weight. For all other time points, one half of a flask culture (375 ml) was harvested, yielding a biomass of approximately 1.5 g (day 8) and 3 g (day 22) of wet cell weight. Thus, a total of 5 separate culture flasks were used for the biomass collection (1 flask each for days 0, 3, and 5; 1 flask for days 8, and 14; one flask for days 18 and 22), and each of these 5 flasks originated from a single 4-week-old B. braunii culture. Samples were harvested by vacuum filtration using 35 μm nylon mesh (Aquatic Ecosystems Inc., Apopka, FL). The harvested colonies were rinsed with sterilized deionized water, frozen in liquid nitrogen, and stored at −80°C for RNA extraction.
Total RNA from each sample was extracted by pulverizing liquid nitrogen-frozen algae using a Tissuelyser II (Qiagen, Valencia, CA). Approximately 200 mg of the pulverized frozen tissue was added to 1 ml of TRIzol reagent (Life Technologies, Grand Island, NY) and total RNA was extracted following the protocol of the manufacturer. The resulting RNA pellet was washed with 75% EtOH, dried in a SpeedVac, and the RNA was selectively precipitated from contaminating polysaccharides by resuspending the RNA in 500 μl of 2M LiCl and incubating at room temperature for 5 min. The sample was centrifuged at 12,000 x g for 15 min to pellet the RNA and the supernatant (containing the polysaccharides) was discarded. The LiCl precipitation steps were repeated until the size of the RNA pellet appeared constant. The final RNA pellet was resuspended in 1x TE, extracted against an equal volume of phenol/CHCl3/isoamyl alcohol (25/24/1) and centrifuged at 10,000 x g. The aqueous supernatant containg the RNA was extracted with CHCl3 and centrifuged at 10,000 x g. The RNA in the aqueous supernatant was then precipitated by adding 0.1 volumes of 3M sodium acetate and 2.5 volumes 100% EtOH, incubating at −20°C for 20 min, and centrifuging at 10,000 x g for 15 min. The RNA pellet was washed with 70% EtOH, dried, and dissolved in 50 μl RNase free dH2O. The RNA was quantified by absorbance at 260 nm (A260), and the A260/A280 (protein contamination) and A260/A230 (polysaccharide contamination) ratios were determined. An A260/A280 ratio of 1.8 or higher and an A260/A230 ratio of 2.0 or higher were deemed as suitable. The later time points (days 18 and 22) contained higher amounts of polysaccharides that were not completely removable by repeated LiCl precipitation and the A260/A230 ratios for these samples were in the range of 1.5 to 1.6. The RNA from each sample was then pooled into one sample and treated at 37°C for 30 min with 1 unit of RQ1 RNAse-free DNAse (Promega, Fitchburg, WI) for each μg of RNA. After phenol/CHCl3/isoamyl alcohol and subseaquent CHCl3 extractions, the RNA was again precipitated with 3M sodium acetate/100% EtOH as above. The resulting RNA was resuspended in 100 μl RNAse-free dH2O and the RNA quantified and protein and polysaccharide contamination determined as above. Finally, the RNA was sent to the DOE Joint Genome Institute (JGI) for library construction and 454 sequencing.
RNA quality assessment
In addition to the A260/A280 and A260/A280 ratios, the quality of the isolated RNA was analyzed by several methods. The RNA from days 0, 3, and 5 were analyzed by RT-PCR for the presence of several cDNAs including squalene synthase (SS) , squalene synthase-like-1 (SSL-1) , β-actin , ferredoxin-1, and ferredoxin NADPH reductase ( Additional file 1: Figure S2). The pooled RNA was analyzed for the presence of active nucleases by incubating one 5-μg aliquot of RNA at 37°C and one 5-μg aliquot of RNA at 0°C for 1 hr, followed by visual inspection for degradation of the RNA on a denaturing agarose gel.
The data were assembled using an in house cleaning, clustering, and assembling pipeline developed for 454 transcriptome sequences . This pipeline uses a combination of assembly methods, which has been shown to produce assemblies that are superior to any single method, resulting in fewer contigs with shorter cumulative length . Raw reads were first cleaned using SnoWhite v1.1.3 (http://dlugoschlab.arizona.edu/software.html) with the TagDust cleaning step implemented . The cleaned reads were clustered using the wcd EST clustering tool with an error threshold value of 5 . Based on the wcd output, a perl script split the dataset into multiple fasta and associated quality files based on cluster sizes. These cluster sizes correspond to read depth (i.e. expression level) allowing the assembly parameters to be tailored for the expression level of transcripts. In this way, transcript contigs are initially assembled in separate bins according to their expression level, ameliorating the effects of large differences in sequence representation in transcriptome datasets. The multiple files were assembled separately using Mira v3.0.0 . The Mira outputs were assembled a second time using CAP3 with an overlap percent identity cutoff of 94 , creating an “assembly of assemblies”. Coverage and read depth in this final assembly were determined by mapping the reads to the assembly in CLC Genomics Workbench (v4.9). The quality of the assembly was benchmarked against the core set of eukaryotic genes using the Core Eukaryotic Genes Mapping Approach (CEGMA) algorithm . Transcript contigs were deposited into the GenBank Transcriptome Shotgun Assembly Sequence Database and were assigned accession numbers KA089548 - KA133805 and KA659919 - KA660048. The number of contigs that might have originated from transcripts of other organisms were estimated using the Metagenomics RAST server  from the taxonomical distribution of the top BlastX hits to proteins in GenBank with the transcripts of the B. braunii Showa database.
The Botryococcus braunii web-based annotation tool and data depository
Contigs were annotated using Blast2GO to assign Gene Ontology classifications . KEGG annotation was done with the KAAS annotation server (http://www.genome.jp/tools/kaas/) using the bi-directional best hit method . Custom Perl sripts (avalible from JDH by request) were used to generate the Venn diagram in Figure 3 to compare KEGG annotations for Chlorophyta genomes and the B. braunii transcriptome, using predicted proteins from the genomes of Chlamydomonas reinhardtii (v3.0, http://genome.jgi-psf.org/Chlre3/Chlre3.download.ftp.html), Chlorella variabilis NC64A (http://genome.jgi-psf.org/ChlNC64A_1/ChlNC64A_1.download.ftp.html) and Micromonas sp. RCC299 (v3.0, http://genome.jgi-psf.org/MicpuN3/MicpuN3.download.ftp.html).
To assign pathway annotations, unique transcript sequences were aligned using the BlastX alignment program to a target database containing KEGG proteins with pathway annotations as described in . Alignments were filtered using an E-value threshold of 10-5, and the top KEGG protein hit for each transcript meeting the threshold was kept. KEGG pathways assigned to the protein in the target database were subsequently assigned to the corresponding transcripts. MetaCyc, Panther, and Reactome pathways were assigned by performing BlastX alignments against UniProt entries with annotations in the respective databases and filtering using an E-value threshold of 10-5. Gene Ontology, MapMan Ontology, and KOG annotations were assigned by identifying orthologous proteins in the model alga C. reinhardtii and the plant Arabidopsis thaliana using reciprocal BlastX and tBlastN alignments and keeping alignments that were pairwise best hits. Annotations assigned to orthologous proteins in those organisms were subsequently transferred to Botryococcus transcripts. Pfam annotations were assigned using the publicly available web-based batch search feature of the Pfam database.
Manual curation of selected transcripts and pathways
Contigs with automated annotations related to terpenoid biosynthesis were extracted from the dataset and confirmed by reciprocal BlastX searches using the non-redundant protein sequences at NCBI as comparators. Further contigs representing target genes were collected by exhaustive Blast searches of the B. braunii Showa contig and singleton databases using BlastStation-Free v 1.1 (TM Software, Inc.). Machine-identified B. braunii Showa contigs and representative proteins from plants were used as sequence baits. Further gene candidates were identified by searching for appropriate conserved domains and enzyme name keywords using the B. braunii Showa web-based annotation tool (http://pathways-pellegrini.mcdb.ucla.edu/botryo1). Contigs and singletons that may overlap and extend contigs of interest were collected by BlastN searches using BlastStation-Free v1.1. Overlapping sequences were assembled into isotigs using Sequencher 5.0 (GeneCodes Corp.) implementing Tablet v1.11 (The James Hutton Institute). Protein models were built by identifying protein-coding regions in “noisy” sequences using FrameDP (http://iant.toulouse.inra.fr/FrameDP/cgi-bin/framedp.cgi) that integrates protein similarities and probabilistic models into a single prediction. Translational start sites were further evaluated using the NetStart 1.0 Server (http://www.cbs.dtu.dk/services/NetStart/). Predicted frameshifts were corrected by scrutinizing the assemblies for the questionable regions using Tablet. Contigs and curated contigs whose deduced protein products showed the highest similarities to fungal, animal or bacterial proteins with no or very low similarities to plant proteins were flagged, and their codon usage frequencies were plotted against those of B. braunii (http://www.kazusa.or.jp/codon/index.html) using the Graphical Codon Usage Analyser (http://gcua.schoedl.de/). Flagged contigs and isotigs whose mean differences of codon usage frequencies and relative adaptiveness values significantly exceeded (by 30% or more) those of non-flagged contigs and isotigs of the same pathway were treated as potentially originating from a different organism. Subcellular localization of deduced proteins with protein models covering the N-terminus were predicted using the TargetP 1.1 Server (http://www.cbs.dtu.dk/services/TargetP/).
Dedicated to the memory of Prof. Michael Cusanovich (1942–2010).
The work in the authors’ laboratories is supported by the Department of Energy (contract DE-EE0003046 to the NAABB consortium). We are thankful to Andrew Koppisch (Northern Arizona Unviersity), Joe Chappell and Tom D. Niehaus (University of Kentucky), and Shigeru Okada (University of Tokyo) for their expertise and their efforts of initiating the B. braunii sequencing project as part of JGI Program CSP 2009, and to the scientists at the JGI who carried out the sequencing and publicly released their data into the NCBI Sequence Reads Archive. We are also grateful to the three anonymous reviewers for their comprehensive and constructive comments that strenghthened this manuscript.
- Stephens E, Ross IL, Mussgnug JH, Wagner LD, Borowitzka MA, Posten C, Kruse O, Hankamer B: Future prospects of microalgal biofuel production systems. Trends Plant Sci. 2010, 15 (10): 554-564.PubMedGoogle Scholar
- Chisti Y: Biodiesel from microalgae. Biotechnol Adv. 2007, 25 (3): 294-306.PubMedGoogle Scholar
- Chisti Y: Biodiesel from microalgae beats bioethanol. Trends Biotechnol. 2008, 26 (3): 126-131.PubMedGoogle Scholar
- Scott SA, Davey MP, Dennis JS, Horst I, Howe CJ, Lea-Smith DJ, Smith AG: Biodiesel from algae: challenges and prospects. Curr Opin Biotechnol. 2010, 21 (3): 277-286.PubMedGoogle Scholar
- Dismukes GC, Carrieri D, Bennette N, Ananyev GM, Posewitz MC: Aquatic phototrophs: efficient alternatives to land-based crops for biofuels. Curr Opin Biotechnol. 2008, 19 (3): 235-240.PubMedGoogle Scholar
- Singh A, Nigam PS, Murphy JD: Mechanism and challenges in commercialisation of algal biofuels. Bioresour Technol. 2011, 102 (1): 26-34.PubMedGoogle Scholar
- Singh A, Nigam PS, Murphy JD: Renewable fuels from algae: An answer to debatable land based fuels. Bioresour Technol. 2011, 102 (1): 10-16.PubMedGoogle Scholar
- Greenwell HC, Laurens LM, Shields RJ, Lovitt RW, Flynn KJ: Placing microalgae on the biofuels priority list: a review of the technological challenges. J R Soc Interface. 2010, 7 (46): 703-726.PubMed CentralPubMedGoogle Scholar
- Banerjee A, Sharma R, Chisti Y, Banerjee UC: Botryococcus braunii: a renewable source of hydrocarbons and other chemicals. Crit Rev Biotechnol. 2002, 22 (3): 245-279.PubMedGoogle Scholar
- Metzger P, Largeau C: Chemicals of Botryococcus braunii. Chemicals from microalgae. Edited by: Cohen Z. 1999, Philadelphia, PA: Taylor & Francis, 205-260.Google Scholar
- Metzger P, Largeau C: Botryococcus braunii: a rich source for hydrocarbons and related ether lipids. Appl Microbiol Biotechnol. 2005, 66 (5): 486-496.PubMedGoogle Scholar
- Brown AC, Knights BA, Conway E: Hydrocarbon content and its relation to physiological state in the green alga Botryococcus braunii. Phytochem. 1969, 8 (3): 543-547.Google Scholar
- Weiss TL, Chun HJ, Okada S, Vitha S, Holzenburg A, Laane J, Devarenne TP: Raman spectroscopy analysis of botryococcene hydrocarbons from the green microalga Botryococcus braunii. J Biol Chem. 2010, 285 (42): 32458-32466.PubMed CentralPubMedGoogle Scholar
- Merchant SS, Prochnik SE, Vallon O, Harris EH, Karpowicz SJ, Witman GB, Terry A, Salamov A, Fritz-Laylin LK, Marechal-Drouard L, et al: The Chlamydomonas genome reveals the evolution of key animal and plant functions. Science. 2007, 318 (5848): 245-250.PubMed CentralPubMedGoogle Scholar
- Blanc G, Duncan G, Agarkova I, Borodovsky M, Gurnon J, Kuo A, Lindquist E, Lucas S, Pangilinan J, Polle J, et al: The Chlorella variabilis NC64A genome reveals adaptation to photosymbiosis, coevolution with viruses, and cryptic sex. Plant Cell. 2010, 22 (9): 2943-2955.PubMed CentralPubMedGoogle Scholar
- Prochnik SE, Umen J, Nedelcu AM, Hallmann A, Miller SM, Nishii I, Ferris P, Kuo A, Mitros T, Fritz-Laylin LK, et al: Genomic analysis of organismal complexity in the multicellular green alga Volvox carteri. Science. 2010, 329 (5988): 223-226.PubMed CentralPubMedGoogle Scholar
- Kroth PG, Chiovitti A, Gruber A, Martin-Jezequel V, Mock T, Parker MS, Stanley MS, Kaplan A, Caron L, Weber T, et al: A model for carbohydrate metabolism in the diatom Phaeodactylum tricornutum deduced from comparative whole genome analysis. PLoS One. 2008, 3 (1): e1426-PubMed CentralPubMedGoogle Scholar
- Chang RL, Ghamsari L, Manichaikul A, Hom EF, Balaji S, Fu W, Shen Y, Hao T, Palsson BO, Salehi-Ashtiani K, et al: Metabolic network reconstruction of Chlamydomonas offers insight into light-driven algal metabolism. Mol Syst Biol. 2011, 7: 518-PubMed CentralPubMedGoogle Scholar
- Khozin-Goldberg I, Cohen Z: Unraveling algal lipid metabolism: Recent advances in gene identification. Biochimie. 2011, 93 (1): 91-100.PubMedGoogle Scholar
- Merchant SS, Kropat J, Liu B, Shaw J, Warakanont J: TAG, You're it! Chlamydomonas as a reference organism for understanding algal triacylglycerol accumulation. Curr Opin Biotechnol. 2012, 23 (3): 352-363.PubMedGoogle Scholar
- Castruita M, Casero D, Karpowicz SJ, Kropat J, Vieler A, Hsieh SI, Yan W, Cokus S, Loo JA, Benning C, et al: Systems biology approach in Chlamydomonas reveals connections between copper nutrition and multiple metabolic steps. Plant Cell. 2011, 23 (4): 1273-1292.PubMed CentralPubMedGoogle Scholar
- Gonzalez-Ballester D, Casero D, Cokus S, Pellegrini M, Merchant SS, Grossman AR: RNA-seq analysis of sulfur-deprived Chlamydomonas cells reveals aspects of acclimation critical for cell survival. Plant Cell. 2010, 22 (6): 2058-2084.PubMed CentralPubMedGoogle Scholar
- 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-PubMed CentralPubMedGoogle Scholar
- Guarnieri MT, Nag A, Smolinski SL, Darzins A, Seibert M, Pienkos PT: Examination of triacylglycerol biosynthetic pathways via de novo transcriptomic and proteomic analyses in an unsequenced microalga. PLoS One. 2011, 6 (10): e25851-PubMed CentralPubMedGoogle Scholar
- 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 (Database issue): D355-D360.PubMed CentralPubMedGoogle Scholar
- Consortium GO: The Gene Ontology in 2010: extensions and refinements. Nucleic Acids Res. 2010, 38 (Database issue): D331-D335.Google Scholar
- Finn RD, Mistry J, Tate J, Coggill P, Heger A, Pollington JE, Gavin OL, Gunasekaran P, Ceric G, Forslund K, et al: The Pfam protein families database. Nucleic Acids Res. 2010, 38 (Database issue): D211-D222.PubMed CentralPubMedGoogle Scholar
- Lopez D, Casero D, Cokus SJ, Merchant SS, Pellegrini M: Algal Functional Annotation Tool: a web-based analysis suite to functionally interpret large gene lists using integrated annotation and expression data. BMC Bioinformatics. 2011, 12 (1): 282-PubMed CentralPubMedGoogle Scholar
- Ioki M, Baba M, Bidadi H, Suzuki I, Shiraiwa Y, Watanabe MM, Nakajima N: Modes of hydrocarbon oil biosynthesis by comparative gene expression analysis of race A and race B strains of Botryococcus braunii. Bioresour Technol. 2012, 109: 271-276.PubMedGoogle Scholar
- Baba M, Ioki M, Nakajima N, Shiraiwa Y, Watanabe MM: Transcriptome analysis of an oil-rich race A strain of Botryococcus braunii (BOT-88-2) by de novo assembly of pyrosequencing cDNA reads. Bioresour Technol. 2012, 109: 282-286.PubMedGoogle Scholar
- Ioki M, Baba M, Nakajima N, Shiraiwa Y, Watanabe MM: Transcriptome analysis of an oil-rich race B strain of Botryococcus braunii (BOT-70) by de novo assembly of 5'-end sequences of full-length cDNA clones. Bioresour Technol. 2012, 109: 277-281.PubMedGoogle Scholar
- Ioki M, Baba M, Nakajima N, Shiraiwa Y, Watanabe MM: Transcriptome analysis of an oil-rich race B strain of Botryococcus braunii (BOT-22) by de novo assembly of pyrosequencing cDNA reads. Bioresour Technol. 2012, 109: 292-296.PubMedGoogle Scholar
- Niehaus T, Kinison S, Okada S, Yeo YS, Bell SA, Cui P, Devarenne TP, Chappell J: Functional identification of triterpene methyltransferases from Botryococcus braunii race B. J Biol Chem. 2012, 287: 8163-8173.PubMed CentralPubMedGoogle Scholar
- Niehaus TD, Okada S, Devarenne TP, Watt DS, Sviripa V, Chappell J: Identification of unique mechanisms for triterpene biosynthesis in Botryococcus braunii. Proc Natl Acad Sci USA. 2011, 108 (30): 12260-12265.PubMed CentralPubMedGoogle Scholar
- Matsushima D, Jenke-Kodama H, Sato Y, Fukunaga Y, Sumimoto K, Kuzuyama T, Matsunaga S, Okada S: The single cellular green microalga Botryococcus braunii, race B possesses three distinct 1-deoxy-D-xylulose 5-phosphate synthases. Plant Sci. 2012, 185–186: 309-320.PubMedGoogle Scholar
- Weiss TL, Johnston JS, Fujisawa K, Sumimoto K, Okada S, Chappell J, Devarenne TP: Phylogenetic placement, genome size, and GC content of the liquid-hydrocarbon-producing green microalga Botryococcus braunii strain Berkeley (Showa) (Chlorophyta). J Phycol. 2010, 46: 534-540.Google Scholar
- Okada S, Devarenne TP, Chappell J: Molecular characterization of squalene synthase from the green microalga Botryococcus braunii, race B. Arch Biochem Biophys. 2000, 373 (2): 307-317.PubMedGoogle Scholar
- Okada S, Devarenne TP, Murakami M, Abe H, Chappell J: Characterization of botryococcene synthase enzyme activity, a squalene synthase-like activity from the green microalga Botryococcus braunii, Race B. Arch Biochem Biophys. 2004, 422 (1): 110-118.PubMedGoogle Scholar
- Wisecaver JH, Hackett JD: Transcriptome analysis reveals nuclear-encoded proteins for the maintenance of temporary plastids in the dinoflagellate Dinophysis acuminata. BMC Genomics. 2010, 11: 366-PubMed CentralPubMedGoogle Scholar
- Kumar S, Blaxter ML: Comparing de novo assemblers for 454 transcriptome data. BMC Genomics. 2010, 11: 571-PubMed CentralPubMedGoogle Scholar
- Parra G, Bradnam K, Korf I: CEGMA: a pipeline to accurately annotate core genes in eukaryotic genomes. Bioinformatics. 2007, 23 (9): 1061-1067.PubMedGoogle Scholar
- Hu H, Bandyopadhyay PK, Olivera BM, Yandell M: Characterization of the Conus bullatus genome and its venom-duct transcriptome. BMC Genomics. 2011, 12: 60-PubMed CentralPubMedGoogle Scholar
- Caspi R, Altman T, Dale JM, Dreher K, Fulcher CA, Gilham F, Kaipa P, Karthikeyan AS, Kothari A, Krummenacker M, et al: The MetaCyc database of metabolic pathways and enzymes and the BioCyc collection of pathway/genome databases. Nucleic Acids Res. 2010, 38 (Database issue): D473-D479.PubMed CentralPubMedGoogle Scholar
- Croft D, O'Kelly G, Wu G, Haw R, Gillespie M, Matthews L, Caudy M, Garapati P, Gopinath G, Jassal B, et al: Reactome: a database of reactions, pathways and biological processes. Nucleic Acids Res. 2011, 39 (Database issue): D691-D697.PubMed CentralPubMedGoogle Scholar
- Mi H, Dong Q, Muruganujan A, Gaudet P, Lewis S, Thomas PD: PANTHER version 7: improved phylogenetic trees, orthologs and collaboration with the Gene Ontology Consortium. Nucleic Acids Res. 2010, 38 (Database issue): D204-D210.PubMed CentralPubMedGoogle Scholar
- Punta M, Coggill PC, Eberhardt RY, Mistry J, Tate J, Boursnell C, Pang N, Forslund K, Ceric G, Clements J, et al: The Pfam protein families database. Nucleic Acids Res. 2012, 40 (Database issue): D290-D301.PubMed CentralPubMedGoogle Scholar
- Thimm O, Blasing O, Gibon Y, Nagel A, Meyer S, Kruger P, Selbig J, Muller LA, Rhee SY, Stitt M: MAPMAN: a user-driven tool to display genomics data sets onto diagrams of metabolic pathways and other biological processes. Plant J. 2004, 37 (6): 914-939.PubMedGoogle Scholar
- Meyer F, Paarmann D, D'Souza M, Olson R, Glass EM, Kubal M, Paczian T, Rodriguez A, Stevens R, Wilke A, et al: The metagenomics RAST server - a public resource for the automatic phylogenetic and functional analysis of metagenomes. BMC Bioinformatics. 2008, 9: 386-PubMed CentralPubMedGoogle Scholar
- Lohr M, Schwender J, Polle JEW: Isoprenoid biosynthesis in eukaryotic phototrophs: A spotlight on algae. Plant Sci. 2011, 185–186: 9-22.PubMedGoogle Scholar
- Kirby J, Keasling JD: Biosynthesis of plant isoprenoids: perspectives for microbial engineering. Annu Rev Plant Biol. 2009, 60: 335-355.PubMedGoogle Scholar
- Metzger P, Rager M-N, Largeau C: Polyacetals based on polymethylsqualene diols, precursors of algaenan in Botryococcus braunii race B. Organic Geochemistry. 2007, 38 (4): 566-581.Google Scholar
- Weiss TL, Spencer JJ, Fujisawa K, Okada S, Devarenne TP: Genome size and phylogenetic analysis of the A and L races of Botryococcus braunii. J Appl Phycol. 2011, 23: 833-839.Google Scholar
- Sato Y, Ito Y, Okada S, Murakami M, Abe H: Biosynthesis of the triterpenoids, botryococcenes and tetramethylsqualene in the B race of Botryococcus braunii via the non-mevalonate pathway. Tetrahedron Lett. 2003, 44 (37): 7035-7037.Google Scholar
- Miziorko HM: Enzymes of the mevalonate pathway of isoprenoid biosynthesis. Arch Biochem Biophys. 2011, 505 (2): 131-143.PubMed CentralPubMedGoogle Scholar
- Schwender J, Seemann M, Lichtenthaler HK, Rohmer M: Biosynthesis of isoprenoids (carotenoids, sterols, prenyl side-chains of chlorophylls and plastoquinone) via a novel pyruvate/glyceraldehyde 3-phosphate non-mevalonate pathway in the green alga Scenedesmus obliquus. Biochem J. 1996, 316 (Pt 1): 73-80.PubMed CentralPubMedGoogle Scholar
- Hunter WN: The non-mevalonate pathway of isoprenoid precursor biosynthesis. J Biol Chem. 2007, 282 (30): 21573-21577.PubMedGoogle Scholar
- Rohmer M: Methylerythritol phosphate pathway. Comprehensive Natural Products II: Chemistry and biology vol. 1. Edited by: Mander L, Liu H-W. 2010, New York, NY: Elsevier Inc, 517-555.Google Scholar
- Cordoba E, Porta H, Arroyo A, San Roman C, Medina L, Rodriguez-Concepcion M, Leon P: Functional characterization of the three genes encoding 1-deoxy-D-xylulose 5-phosphate synthase in maize. J Exp Bot. 2011, 62 (6): 2023-2038.PubMedGoogle Scholar
- Rohdich F, Hecht S, Gartner K, Adam P, Krieger C, Amslinger S, Arigoni D, Bacher A, Eisenreich W: Studies on the nonmevalonate terpene biosynthetic pathway: metabolic role of IspH (LytB) protein. Proc Natl Acad Sci U S A. 2002, 99 (3): 1158-1163.PubMed CentralPubMedGoogle Scholar
- Poliquin K, Ershov YV, Cunningham FX, Woreta TT, Gantt RR, Gantt E: Inactivation of sll1556 in Synechocystis strain PCC 6803 impairs isoprenoid biosynthesis from pentose phosphate cycle substrates in vitro. J Bacteriol. 2004, 186 (14): 4685-4693.PubMed CentralPubMedGoogle Scholar
- Kaneda K, Kuzuyama T, Takagi M, Hayakawa Y, Seto H: An unusual isopentenyl diphosphate isomerase found in the mevalonate pathway gene cluster from Streptomyces sp. strain CL190. Proc Natl Acad Sci U S A. 2001, 98 (3): 932-937.PubMed CentralPubMedGoogle Scholar
- Ershov YV, Gantt RR, Cunningham FX, Gantt E: Isoprenoid biosynthesis in Synechocystis sp. strain PCC6803 is stimulated by compounds of the pentose phosphate cycle but not by pyruvate or deoxyxylulose-5-phosphate. J Bacteriol. 2002, 184 (18): 5045-5051.PubMed CentralPubMedGoogle Scholar
- Chen F, Tholl D, Bohlmann J, Pichersky E: The family of terpene synthases in plants: a mid-size family of genes for specialized metabolism that is highly diversified throughout the kingdom. Plant J. 2011, 66 (1): 212-229.PubMedGoogle Scholar
- Degenhardt J, Kollner TG, Gershenzon J: Monoterpene and sesquiterpene synthases and the origin of terpene skeletal diversity in plants. Phytochemistry. 2009, 70 (15–16): 1621-1637.PubMedGoogle Scholar
- Dewick PM: The biosynthesis of C5-C25 terpenoid compounds. Nat Prod Rep. 2002, 19 (2): 181-222.PubMedGoogle Scholar
- Grung M, Metzger P, Liaaen-Jensen S: Algal carotenoids. Part 42. Primary and secondary carotenoids in two races of the green algaBotryococcus braunii. Biochem Syst Ecol. 1989, 17 (4): 263-269.Google Scholar
- Blagg BS, Jarstfer MB, Rogers DH, Poulter CD: Recombinant squalene synthase. A mechanism for the rearrangement of presqualene diphosphate to squalene. J Am Chem Soc. 2002, 124 (30): 8846-8853.PubMedGoogle Scholar
- Polle JEW, Tran D: Regulating the production of long chain hydrocarbons. 2009, The City University of New York, USA: World Intellectual Property Organization, 158433-Google Scholar
- Metzger P, Casadevall E, Pouet MJ, Pouet Y: Structures of some botryococcenes: branched hydrocarbons from the B-race of the green alga Botryococcus braunii. Phytochemistry. 1985, 24 (12): 2995-3002.Google Scholar
- Metzger P, Casadevall E, Coute A: Botryococcene distribution in strains of the green alga Botryococcus braunii. Phytochemistry. 1988, 27 (5): 1383-1388.Google Scholar
- Wolf FR, Nonomura AM, Bassham JA: Growth and branched hydrocarbon production in a strain of Botryococcus braunii (Chlorophyta). J Phycol. 1985, 21 (3): 388-396.Google Scholar
- Achitouv E, Metzger P, Rager M-N, Largeau C: C31-C34 methylated squalenes from a Bolivian strain of Botryococcus braunii. Phytochem. 2004, 65 (23): 3159-3165.Google Scholar
- Diener AC, Li H, Zhou W, Whoriskey WJ, Nes WD, Fink GR: Sterol methyltransferase 1 controls the level of cholesterol in plants. Plant Cell. 2000, 12 (6): 853-870.PubMed CentralPubMedGoogle Scholar
- Nes WD: Enzyme mechanisms for sterol C-methylations. Phytochemistry. 2003, 64 (1): 75-95.PubMedGoogle Scholar
- Metzger P, David M, Casadevall E: Biosynthesis of triterpenoid hydrocarbons in the B-race of the green algaBotryococcus braunii. Sites on production and nature of the methylating agent. Phytochemistry. 1986, 26 (1): 129-134.Google Scholar
- Bouvier-Nave P, Husselstein T, Benveniste P: Two families of sterol methyltransferases are involved in the first and the second methylation steps of plant sterol biosynthesis. Eur J Biochem. 1998, 256 (1): 88-96.PubMedGoogle Scholar
- Okada S, Tonegawa I, Matsuda H, Murakami M, Yamaguchi K: Braunixanthins 1 and 2, new carotenoids from the green microalga Botryococcus braunii. Tetrahedron. 1997, 53 (33): 11307-11316.Google Scholar
- Behrman EJ, Gopalan V: Cholesterol and plants. J Chem Educ. 2005, 82: 1791-1793.Google Scholar
- Metzger P, Allard B, Casadevall E, Berkaloff C, Coute A: Structure and chemistry of a new chemical race of Botryococcus braunii (Chlorophyceae) that produces lycopadiene, a tetraterpenoid hydrocarbon. J Phycol. 1990, 26 (2): 258-266.Google Scholar
- Barupal DK, Kind T, Kothari SL, Lee Do Y, Fiehn O: Hydrocarbon phenotyping of algal species using pyrolysis-gas chromatography mass spectrometry. BMC Biotechnol. 2010, 10: 40-PubMed CentralPubMedGoogle Scholar
- Brumfield KM, Moroney JV, Moore TS, Simms TA, Donze D: Functional characterization of the Chlamydomonas reinhardtii ERG3 ortholog, a gene involved in the biosynthesis of ergosterol. PLoS One. 2010, 5 (1): e8659-PubMed CentralPubMedGoogle Scholar
- Metzger P, Casadevall E: Aldehydes, very long chain alkenylphenols, epoxides and other lipids from an alkadiene-producing strain of Botryococcus braunii. Phytochemistry. 1989, 28 (8): 2097-2104.Google Scholar
- Ruban AV, Johnson MP: Xanthophylls as modulators of membrane protein function. Arch Biochem Biophys. 2010, 504 (1): 78-85.PubMedGoogle Scholar
- Xiao FG, Shen L, Ji HF: On photoprotective mechanisms of carotenoids in light harvesting complex. Biochem Biophys Res Commun. 2011, 414 (1): 1-4.PubMedGoogle Scholar
- Takaichi S: Carotenoids in algae: distributions, biosyntheses and functions. Mar Drugs. 2011, 9 (6): 1101-1118.PubMed CentralPubMedGoogle Scholar
- Ranga Rao A, Raghunath Reddy RL, Baskaran V, Sarada R, Ravishankar GA: Characterization of microalgal carotenoids by mass spectrometry and their bioavailability and antioxidant properties elucidated in rat model. J Agric Food Chem. 2010, 58 (15): 8553-8559.PubMedGoogle Scholar
- Tonegawa I, Okada S, Murakami M, Yamaguchi K: Pigment composition of the green microalga Botryococcus braunii Kawaguchi-1. Fisheries Science. 1998, 64 (2): 305-308.Google Scholar
- Okada S, Tonegawa I, Matsuda H, Murakami M, Yamaguchi K: Botryoxanthin B and α-botryoxanthin A from the green microalga Botryococcus braunii Kawaguchi-1. Phytochem. 1998, 47 (6): 1111-1115.Google Scholar
- Okada S, Matsuda H, Murakami M, Yamaguchi K: Botryoxanthin A, a member of a new class of carotenoids from the green microalga Botryococcus braunii Berkeley. Tetrahedron Lett. 1996, 37 (7): 1065-1068.Google Scholar
- Zhang Z, Metzger P, Sachs JP: Co-occurrence of long chain diols, keto-ols, hydroxy acids and keto acids in recent sediments of Lake El Junco, Galápagos Islands. Organic Geochemistry. 2011, 42: 823-837.Google Scholar
- Falk J, Munne-Bosch S: Tocochromanol functions in plants: antioxidation and beyond. J Exp Bot. 2010, 61 (6): 1549-1566.PubMedGoogle Scholar
- Lochner K, Doring O, Bottger M: Phylloquinone, what can we learn from plants?. Biofactors. 2003, 18 (1–4): 73-78.PubMedGoogle Scholar
- Swiezewska E: Ubiquinone and plastoquinone metabolism in plants. Methods Enzymol. 2004, 378: 124-131.PubMedGoogle Scholar
- Depuydt S, Hardtke CS: Hormone signalling crosstalk in plant growth regulation. Curr Biol. 2011, 21 (9): R365-R373.PubMedGoogle Scholar
- Tarakhovskaya ER, Maslov YI, Shishova MF: Phytohormones in algae. Russ J Plant Physisol. 2007, 54 (2): 163-170.Google Scholar
- Yasuno R, Wada H: The biosynthetic pathway for lipoic acid is present in plastids and mitochondria in Arabidopsis thaliana. FEBS Lett. 2002, 517 (1–3): 110-114.PubMedGoogle Scholar
- Riekhof WR, Benning C: Glycerolipid biosynthesis. The Chlamydomonas sourcebook: Organellar and metabolic processes vol. 2. Edited by: Stern DB. 2009, San Diego, CA: Academic Press, 41-68.Google Scholar
- Sasaki Y, Nagano Y: Plant acetyl-CoA carboxylase: structure, biosynthesis, regulation, and gene manipulation for plant breeding. Biosci Biotechnol Biochem. 2004, 68 (6): 1175-1184.PubMedGoogle Scholar
- Harwood JL, Guschina IA: The versatility of algae and their lipid metabolism. Biochimie. 2008, 91 (6): 679-684.PubMedGoogle Scholar
- Guschina IA, Harwood JL: Lipids and lipid metabolism in eukaryotic algae. Prog Lipid Res. 2006, 45 (2): 160-186.PubMedGoogle Scholar
- Ranga Rao A, Sarada R, Ravishankar GA: Influence of CO2 on growth and hydrocarbon production in Botryococcus braunii. J Microbiol Biotechnol. 2007, 17 (3): 414-419.PubMedGoogle Scholar
- Ben-Amotz A, Tornabene TG: Chemical profile of selected species of microalgae with emphasis on lipids. J Phycol. 1985, 21: 72-81.Google Scholar
- MacDougall KM, McNichol J, McGinn PJ, O'Leary SJ, Melanson JE: Triacylglycerol profiling of microalgae strains for biofuel feedstock by liquid chromatography-high-resolution mass spectrometry. Anal Bioanal Chem. 2011, 401 (8): 2609-2616.PubMed CentralPubMedGoogle Scholar
- Nguyen HM, Baudet M, Cuine S, Adriano JM, Barthe D, Billon E, Bruley C, Beisson F, Peltier G, Ferro M, et al: Proteomic profiling of oil bodies isolated from the unicellular green microalga Chlamydomonas reinhardtii: with focus on proteins involved in lipid metabolism. Proteomics. 2011, 11 (21): 4266-4273.PubMedGoogle Scholar
- John RP, Anisha GS, Nampoothiri KM, Pandey A: Micro and macroalgal biomass: a renewable source for bioethanol. Bioresour Technol. 2011, 102 (1): 186-193.PubMedGoogle Scholar
- Rude MA, Schirmer A: New microbial fuels: a biotech perspective. Curr Opin Microbiol. 2009, 12 (3): 274-281.PubMedGoogle Scholar
- Ball SG, Morell MK: From bacterial glycogen to starch: understanding the biogenesis of the plant starch granule. Annu Rev Plant Biol. 2003, 54: 207-233.PubMedGoogle Scholar
- Rodrigues MA, da Silva Bon EP: Evaluation of Chlorella (Chlorophyta) as source of fermentable sugars via cell wall enzymatic hydrolysis. Enzyme Res. 2011, 2011: 405603-PubMed CentralPubMedGoogle Scholar
- Popper ZA, Tuohy MG: Beyond the green: understanding the evolutionary puzzle of plant and algal cell walls. Plant Physiol. 2010, 153 (2): 373-383.PubMed CentralPubMedGoogle Scholar
- Popper ZA, Michel G, Herve C, Domozych DS, Willats WGT, Tuohy MG, Kloareg B, Stengel DB: Evolution and diversity of plant cell walls: from algae to flowering plants. Annu Rev Plant Biol. 2011, 62: 567-590.PubMedGoogle Scholar
- Yin Y, Huang J, Xu Y: The cellulose synthase superfamily in fully sequenced plants and algae. BMC plant biology. 2009, 9: 99-PubMed CentralPubMedGoogle Scholar
- Dunlop MJ, Dossani ZY, Szmidt HL, Chu HC, Lee TS, Keasling JD, Hadi MZ, Mukhopadhyay A: Engineering microbial biofuel tolerance and export using efflux pumps. Mol Syst Biol. 2011, 7: 487-PubMed CentralPubMedGoogle Scholar
- Yazaki K, Shitan N, Sugiyama A, Takanashi K: Cell and molecular biology of ATP-binding cassette proteins in plants. Int Rev Cell Mol Biol. 2009, 276: 263-299.PubMedGoogle Scholar
- Verrier PJ, Bird D, Burla B, Dassa E, Forestier C, Geisler M, Klein M, Kolukisaoglu U, Lee Y, Martinoia E, et al: Plant ABC proteins - a unified nomenclature and updated inventory. Trends Plant Sci. 2008, 13 (4): 151-159.PubMedGoogle Scholar
- Berkaloff C, Rousseau B, Coute A, Casadevall E, Metzger P, Chirac C: Variability of cell wall structure and hydrocarbon type in different strains of Botryococcus braunii. J Phycol. 1984, 20 (3): 377-389.Google Scholar
- Glick D, Barth S, Macleod KF: Autophagy: cellular and molecular mechanisms. J Pathol. 2010, 221 (1): 3-12.PubMed CentralPubMedGoogle Scholar
- Lockshin RA, Zakeri Z: Apoptosis, autophagy, and more. Int J Biochem Cell Biol. 2004, 36 (12): 2405-2419.PubMedGoogle Scholar
- Perez-Perez ME, Crespo JL: Autophagy in the model alga Chlamydomonas reinhardtii. Autophagy. 2010, 6 (4): 562-563.PubMedGoogle Scholar
- Perez-Perez ME, Florencio FJ, Crespo JL: Inhibition of target of rapamycin signaling and stress activate autophagy in Chlamydomonas reinhardtii. Plant Physiol. 2010, 152 (4): 1874-1888.PubMed CentralPubMedGoogle Scholar
- Withers ST, Keasling JD: Biosynthesis and engineering of isoprenoid small molecules. Appl Microbiol Biotechnol. 2007, 73 (5): 980-990.PubMedGoogle Scholar
- Klein-Marcuschamer D, Ajikumar PK, Stephanopoulos G: Engineering microbial cell factories for biosynthesis of isoprenoid molecules: beyond lycopene. Trends Biotechnol. 2007, 25 (9): 417-424.PubMedGoogle Scholar
- Zhang F, Rodriguez S, Keasling JD: Metabolic engineering of microbial pathways for advanced biofuels production. Curr Opin Biotechnol. 2011, 22 (6): 775-783.PubMedGoogle Scholar
- Kumar S, Hahn FM, Baidoo E, Kahlon TS, Wood DF, McMahan CM, Cornish K, Keasling JD, Daniell H, Whalen MC: Remodeling the isoprenoid pathway in tobacco by expressing the cytoplasmic mevalonate pathway in chloroplasts. Metab Eng. 2011, 14 (1): 19-28.PubMedGoogle Scholar
- Goodson C, Roth R, Wang ZT, Goodenough U: Structural correlates of cytoplasmic and chloroplast lipid body synthesis in Chlamydomonas reinhardtii and stimulation of lipid body production with acetate boost. Eukaryot Cell. 2011, 10 (12): 1592-1606.PubMed CentralPubMedGoogle Scholar
- Work VH, Radakovits R, Jinkerson RE, Meuser JE, Elliott LG, Vinyard DJ, Laurens LM, Dismukes GC, Posewitz MC: Increased lipid accumulation in the Chlamydomonas reinhardtii sta7-10 starchless isoamylase mutant and increased carbohydrate synthesis in complemented strains. Eukaryot Cell. 2010, 9 (8): 1251-1261.PubMed CentralPubMedGoogle Scholar
- Nonomura AM: Botryococcus braunii var. Showa (Chlorophyceae) from Berkeley, California, United States of America. Jpn J Phycol. 1988, 36: 285-291.Google Scholar
- Lassmann T, Hayashizaki Y, Daub CO: TagDust - a program to eliminate artifacts from next generation sequencing data. Bioinformatics. 2009, 25 (21): 2839-2840.PubMed CentralPubMedGoogle Scholar
- Hazelhurst S, Hide W, Liptak Z, Nogueira R, Starfield R: An overview of the wcd EST clustering tool. Bioinformatics. 2008, 24 (13): 1542-1546.PubMed CentralPubMedGoogle Scholar
- Chevreux B, Pfisterer T, Drescher B, Driesel AJ, Müller WE, Wetter T, Suhai S: Using the miraEST assembler for reliable and automated mRNA transcript assembly and SNP detection in sequenced ESTs. Genome Res. 2004, 14 (6): 1147-1159.PubMed CentralPubMedGoogle Scholar
- Huang X, Madan A: CAP3: A DNA sequence assembly program. Genome Res. 1999, 9 (9): 868-877.PubMed CentralPubMedGoogle Scholar
- Conesa A, Götz S, Garcia-Gómez JM, Terol J, Talón M, Robles M: Blast2GO: a universal tool for annotation, visualization and analysis in functional genomics research. Bioinformatics. 2005, 21 (18): 3674-3676.PubMedGoogle Scholar
- Moriya Y, Itoh M, Okuda S, Yoshizawa AC, Kanehisa M: KAAS: an automatic genome annotation and pathway reconstruction server. Nucleic Acids Res. 2007, 35 (Web Server issue): W182-W185.PubMed CentralPubMedGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.