The genome of the polar eukaryotic microalga Coccomyxa subellipsoidea reveals traits of cold adaptation

A Walkthrough to the Use of GreeNC: The Plant lncRNA Database

Experimentally validated plant lncRNAs have been shown to regulate important agronomic traits such as phosphate starvation response, flowering time, and interaction with symbiotic organisms, making them of great interest in plant biology and in breeding. We developed a pipeline to annotate lncRNAs and applied it to 37 plant species and 6 algae, resulting in the annotation of more than 120,000 lncRNAs. To facilitate the study of lncRNAs for the plant research community, the information gathered is organized in the Green Non-Coding Database (GreeNC, http://greenc.sciencedesigners.com/) . This chapter contains a detailed explanation of the content of GreeNC and how to access both programmatically and with a web browser.

Evolutionary, structural and functional analysis of the caleosin/peroxygenase gene family in the Fungi

BACKGROUND

Caleosin/peroxygenases, CLO/PXG, (designated PF05042 in Pfam) are a group of genes/proteins with anomalous distributions in eukaryotic taxa. We have previously characterised CLO/PXGs in the Viridiplantae. The aim of this study was to investigate the evolution and functions of the CLO/PXGs in the Fungi and other non-plant clades and to elucidate the overall origin of this gene family.

RESULTS

CLO/PXG-like genes are distributed across the full range of fungal groups from the basal clades, Cryptomycota and Microsporidia, to the largest and most complex Dikarya species. However, the genes were only present in 243 out of 844 analysed fungal genomes. CLO/PXG-like genes have been retained in many pathogenic or parasitic fungi that have undergone considerable genomic and structural simplification, indicating that they have important functions in these species. Structural and functional analyses demonstrate that CLO/PXGs are multifunctional proteins closely related to similar proteins found in all major taxa of the Chlorophyte Division of the Viridiplantae. Transcriptome and physiological data show that fungal CLO/PXG-like genes have complex patterns of developmental and tissue-specific expression and are upregulated in response to a range of biotic and abiotic stresses as well as participating in key metabolic and developmental processes such as lipid metabolism, signalling, reproduction and pathogenesis. Biochemical data also reveal that the Aspergillus flavus CLO/PXG has specific functions in sporulation and aflatoxin production as well as playing roles in lipid droplet function.

CONCLUSIONS

In contrast to plants, CLO/PXGs only occur in about 30% of sequenced fungal genomes but are present in all major taxa. Fungal CLO/PXGs have similar but not identical roles to those in plants, including stress-related oxylipin signalling, lipid metabolism, reproduction and pathogenesis. While the presence of CLO/PXG orthologs in all plant genomes sequenced to date would suggest that they have core housekeeping functions in plants, the selective loss of CLO/PXGs in many fungal genomes suggests more restricted functions in fungi as accessory genes useful in particular environments or niches. We suggest an ancient origin of CLO/PXG-like genes in the 'last eukaryotic common ancestor' (LECA) and their subsequent loss in ancestors of the Metazoa, after the latter had diverged from the ancestral fungal lineage.

Out in the Cold: Identification of Genomic Regions Associated With Cold Tolerance in the Biocontrol Fungus Clonostachys rosea Through Genome-Wide Association Mapping

There is an increasing importance for using biocontrol agents in combating plant diseases sustainably and in the long term. As large scale genomic sequencing becomes economically viable, the impact of single nucleotide polymorphisms (SNPs) on biocontrol-associated phenotypes can be easily studied across entire genomes of fungal populations. Here, we improved a previously reported genome assembly of the biocontrol fungus Clonostachys rosea strain IK726 using the PacBio sequencing platform, which resulted in a total genome size of 70.7 Mbp and 21,246 predicted genes. We further performed whole-genome re-sequencing of 52 additional C. rosea strains isolated globally using Illumina sequencing technology, in order to perform genome-wide association studies in conditions relevant for biocontrol activity. One such condition is the ability to grow at lower temperatures commonly encountered in cryic or frigid soils in temperate regions, as these will be prevalent for protecting growing crops in temperate climates. Growth rates at 10°C on potato dextrose agar of the 53 sequenced strains of C. rosea were measured and ranged between 0.066 and 0.413 mm/day. Performing a genome wide association study, a total of 1,478 SNP markers were significantly associated with the trait and located in 227 scaffolds, within or close to (< 1000 bp distance) 265 different genes. The predicted gene products included several chaperone proteins, membrane transporters, lipases, and proteins involved in chitin metabolism with possible roles in cold tolerance. The data reported in this study provides a foundation for future investigations into the genetic basis for cold tolerance in fungi, with important implications for biocontrol.

Contrasting Symbiotic Patterns in Two Closely Related Lineages of Trimembered Lichens of the Genus Peltigera

Species circumscription is key to the characterization of patterns of specificity in symbiotic systems at a macroevolutionary scale. Here, a worldwide phylogenetic framework was used to assess the biodiversity and symbiotic patterns of association among partners in trimembered lichens from the genus Peltigera, section Chloropeltigera. We sequenced six loci of the main fungal partner and performed species discovery and validation analyses to establish putative species boundaries. Single locus phylogenies were used to establish the identity of both photobionts, Nostoc (cyanobacterium) and Coccomyxa (green alga). Distribution and specificity patterns were compared to the closely related clade, section Peltidea, which includes mainly Peltigera species with trimembered thalli. For section Chloropeltigera, eight fungal species (including five newly delimited putative species) were found in association with nine Nostoc phylogroups and two Coccomyxa species. In contrast, eight fungal species (including three newly delimited putative species) in section Peltidea were found in association with only four Nostoc phylogroups and the same two Coccomyxa species as for section Chloropeltigera. This difference in cyanobiont biodiversity between these two sections can potentially be explained by a significantly higher frequency of sexual reproductive structures in species from section Chloropeltigera compared to section Peltidea. Therefore, horizontal transmission of the cyanobiont might be more prevalent in Chloropeltigera species, while vertical transmission might be more common in Peltidea species. All Peltigera species in section Chloropeltigera are generalists in their association with Nostoc compared to more specialized Peltigera species in section Peltidea. Constrained distributions of Peltigera species that associate strictly with one species of green algae (Coccomyxa subellipsoidea) indicate that the availability of the green alga and the specificity of the interaction might be important factors limiting geographic ranges of trimembered Peltigera, in addition to constraints imposed by their interaction with Nostoc partners and by climatic factors.

Coccomyxagreatwallensis sp. nov. (Trebouxiophyceae, Chlorophyta), a lichen epiphytic alga from Fildes Peninsula, Antarctica

A single-celled green alga Coccomyxagreatwallensis Shunan Cao & Qiming Zhou, sp. nov., isolated from a specimen of Antarctic lichen Psoromahypnorum (Vahl) Gray, is described and illustrated based on a comprehensive investigation of morphology, ultrastructure, ecology and phylogeny. The cells of C.greatwallensis are ovoid to long ellipsoidal and measured 3-5 µm × 6-12 µm. The new species has distinct ITS rDNA and SSU rDNA sequences and differs from the phylogenetic closely related species C.antarctica, C.arvernensis and C.viridis in cell size, distribution and habitat.

Predicting microalgae growth and phosphorus removal in cold region waste stabilization ponds using a stochastic modelling approach

A stochastic ecological model with an integrated equilibrium temperature model was developed to predict microalgae growth and phosphorus removal in cold region waste stabilization ponds (WSPs). The model utilized a Monte Carlo simulation to account for parameter uncertainty. The equilibrium temperature model was parameterized using field data collected from two WSPs in Nunavut, Canada, from 2012 to 2014. The equilibrium temperature model provided good agreement with field data on a daily time step. The full model was run using historic (1956-2005) temperature and solar radiation data from five communities (Baker Lake, Cambridge Bay, Coral Harbour, Hall Beach, Resolute) in Nunavut, Canada. The communities represented a range of geographical locations and environmental conditions. Logistic regression on pooled model outputs showed that mean July temperature and mean treatment season temperature (June 1-September 15, ice-free period) provided the best predictors for microalgae growth. They had a predictive success rate of 93 and 88%, respectively. The modelled threshold (50% probability from the Monte Carlo simulation) for microalgae growth was 8.7 and 5.6 °C for the July temperature and mean treatment season temperature, respectively. The logistic regression was applied to each community (except Sanikiluaq) in Nunavut using historic climate data and a probability of microalgae growth was calculated. Based on the model results, soluble phosphorus concentrations consistent with secondary treatment could be achieved if WSP depth is less than 2 m. The model demonstrated a robust method to predict whether a microalgae bloom will occur under a range of model parameters.

Induction of oil accumulation by heat stress is metabolically distinct from N stress in the green microalgae Coccomyxa subellipsoidea C169

Algae are often promoted as feedstock organisms to produce a sustainable petroleum fossil fuel alternative. However, to induce lipid accumulation most often requires a severe stress that is difficult to induce in large batch cultures. The objective of this study is to analyze and mathematically model heat stress on growth, chlorophyll content, triacylglyceride, and starch synthesis in algae. We initially screened 30 algal species for the most pronounced induction of lipid droplets from heat stress using confocal microscopy and mass spectroscopy techniques. One species, Coccomyxa subellipsoidea C169, was selected and subjected to further biochemical analyses using a jacketed bioreactor amended with 1% CO2 at 25°C, 30°C, 32°C, 33°C, 34°C, 35°C, and 36°C. Lipid and starch accumulation was less extreme than N stress. Growth was reduced above 25°C, but heat stress induced lipid droplet synthesis was negatively correlated with growth only past a demonstrated threshold temperature above 32°C. The optimal temperature for lipid accumulation was 35°C, which led to 6% of dry weight triglyceride content and a 72% reduction from optimal growth after 5 days. Fatty acid influx rates into triglycerides and 15N labeling of amino acids and proteins indicate that heat stress is mechanistically distinct from N stress. Thus, this study lends support to a novel hypothesis that lipid droplet triglycerides result from a redistribution of carbon flux as fatty acids to neutral storage lipids over membrane or other lipids.

The potential effect of low cell osmolarity on cell function through decreased concentration of enzyme substrates

Some freshwater algae have lower (<130 osmol m-3) intracellular osmolarities than most others (>180 osmol m-3). Low osmolarities are related to the presence of flagella and the low energy cost of active water efflux following downhill water influx unconstrained by cell walls covering the plasmalemma, and the low resource cost of cell wall synthesis with the same mechanical degree of safety. One consequence of low intracellular osmolarity is limitation on the concentration of metabolites, that is, substrates and products of enzyme activity. Models of the flux through metabolic pathways, and hence the specific growth rate, using steady-state concentrations of enzymes and metabolites have involved organisms with intracellular metabolite osmolarities >280 osmol m-3, where the metabolite concentrations are much greater than the total osmolarity of some freshwater algae. Since the protein concentration (mol m-3) in the cells and the specific growth rates of freshwater cells with low and with higher intracellular osmolarity are highly similar, the models of trade-offs between enzyme and metabolite concentrations for cells with high intracellular osmolarity need modification for cells with low intracellular osmolarity. The soluble free-radical scavenger ascorbate can constitute as little as 0.2% of the low intracellular metabolite concentration (mol m-3) of low-intracellular-osmolarity cells.

Isolation and characterization of 4-hydroxy-3-methylbut-2-enyl diphosphate reductase gene from Botryococcus braunii, race B

The B race of a green microalga Botryococcus braunii Kützing produces triterpene hydrocarbons that is a promising source for biofuel. In this algal race, precursors of triterpene hydrocarbons are provided from the 2-C-methyl-D-erythritol 4-phosphate (MEP) pathway. The terminal enzyme of this pathway, 4-hydroxy-3-methylbut-2-enyl diphosphate reductase (HDR) is regarded as one of the key enzymes that affect yields of products in terpene biosynthesis. In order to better understand the MEP pathway of the alga, cDNA and genomic clones of HDR were obtained from B. braunii Showa strain. B. braunii HDR (BbHDR) is encoded on a single copy gene including a 1509-bp open reading frame that was intervened by 6 introns. The exon-intron structure of BbHDR genes did not show clear relation to phylogeny, while its amino acid sequence reflected phyla and classes well. BbHDR sequence was distinctive from that of the HDR protein from Escherichia coli in the residues involved in hydrogen-bond network that surrounds substrate. Introduction of BbHDR cDNA into an E. coli HDR deficient mutant resulted in recovery of its auxotrophy. BbHDR expression level was upregulated from the onset of liquid culture to the 24th day after inoculation with a 2.5-fold increase and retained its level in the subsequent period.

Metabolic co-dependence drives the evolutionarily ancient Hydra-Chlorella symbiosis

Many multicellular organisms rely on symbiotic associations for support of metabolic activity, protection, or energy. Understanding the mechanisms involved in controlling such interactions remains a major challenge. In an unbiased approach we identified key players that control the symbiosis between Hydra viridissima and its photosynthetic symbiont Chlorella sp. A99. We discovered significant up-regulation of Hydra genes encoding a phosphate transporter and glutamine synthetase suggesting regulated nutrition supply between host and symbionts. Interestingly, supplementing the medium with glutamine temporarily supports in vitro growth of the otherwise obligate symbiotic Chlorella, indicating loss of autonomy and dependence on the host. Genome sequencing of Chlorella sp. A99 revealed a large number of amino acid transporters and a degenerated nitrate assimilation pathway, presumably as consequence of the adaptation to the host environment. Our observations portray ancient symbiotic interactions as a codependent partnership in which exchange of nutrients appears to be the primary driving force.

A transposable element annotation pipeline and expression analysis reveal potentially active elements in the microalga Tisochrysis lutea

BACKGROUND

Transposable elements (TEs) are mobile DNA sequences known as drivers of genome evolution. Their impacts have been widely studied in animals, plants and insects, but little is known about them in microalgae. In a previous study, we compared the genetic polymorphisms between strains of the haptophyte microalga Tisochrysis lutea and suggested the involvement of active autonomous TEs in their genome evolution.

RESULTS

To identify potentially autonomous TEs, we designed a pipeline named PiRATE (Pipeline to Retrieve and Annotate Transposable Elements, download: https://doi.org/10.17882/51795 ), and conducted an accurate TE annotation on a new genome assembly of T. lutea. PiRATE is composed of detection, classification and annotation steps. Its detection step combines multiple, existing analysis packages representing all major approaches for TE detection and its classification step was optimized for microalgal genomes. The efficiency of the detection and classification steps was evaluated with data on the model species Arabidopsis thaliana. PiRATE detected 81% of the TE families of A. thaliana and correctly classified 75% of them. We applied PiRATE to T. lutea genomic data and established that its genome contains 15.89% Class I and 4.95% Class II TEs. In these, 3.79 and 17.05% correspond to potentially autonomous and non-autonomous TEs, respectively. Annotation data was combined with transcriptomic and proteomic data to identify potentially active autonomous TEs. We identified 17 expressed TE families and, among these, a TIR/Mariner and a TIR/hAT family were able to synthesize their transposase. Both these TE families were among the three highest expressed genes in a previous transcriptomic study and are composed of highly similar copies throughout the genome of T. lutea. This sum of evidence reveals that both these TE families could be capable of transposing or triggering the transposition of potential related MITE elements.

CONCLUSION

This manuscript provides an example of a de novo transposable element annotation of a non-model organism characterized by a fragmented genome assembly and belonging to a poorly studied phylum at genomic level. Integration of multi-omics data enabled the discovery of potential mobile TEs and opens the way for new discoveries on the role of these repeated elements in genomic evolution of microalgae.

A transposable element annotation pipeline and expression analysis reveal potentially active elements in the microalga Tisochrysis lutea

Background

Transposable elements (TEs) are mobile DNA sequences known as drivers of genome evolution. Their impacts have been widely studied in animals, plants and insects, but little is known about them in microalgae. In a previous study, we compared the genetic polymorphisms between strains of the haptophyte microalga Tisochrysis lutea and suggested the involvement of active autonomous TEs in their genome evolution.

Results

To identify potentially autonomous TEs, we designed a pipeline named PiRATE (Pipeline to Retrieve and Annotate Transposable Elements, download: 10.17882/51795), and conducted an accurate TE annotation on a new genome assembly of T. lutea. PiRATE is composed of detection, classification and annotation steps. Its detection step combines multiple, existing analysis packages representing all major approaches for TE detection and its classification step was optimized for microalgal genomes. The efficiency of the detection and classification steps was evaluated with data on the model species Arabidopsis thaliana. PiRATE detected 81% of the TE families of A. thaliana and correctly classified 75% of them. We applied PiRATE to T. lutea genomic data and established that its genome contains 15.89% Class I and 4.95% Class II TEs. In these, 3.79 and 17.05% correspond to potentially autonomous and non-autonomous TEs, respectively. Annotation data was combined with transcriptomic and proteomic data to identify potentially active autonomous TEs. We identified 17 expressed TE families and, among these, a TIR/Mariner and a TIR/hAT family were able to synthesize their transposase. Both these TE families were among the three highest expressed genes in a previous transcriptomic study and are composed of highly similar copies throughout the genome of T. lutea. This sum of evidence reveals that both these TE families could be capable of transposing or triggering the transposition of potential related MITE elements.

Conclusion

This manuscript provides an example of a de novo transposable element annotation of a non-model organism characterized by a fragmented genome assembly and belonging to a poorly studied phylum at genomic level. Integration of multi-omics data enabled the discovery of potential mobile TEs and opens the way for new discoveries on the role of these repeated elements in genomic evolution of microalgae.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-4763-1) contains supplementary material, which is available to authorized users.

Evolutionary and genomic analysis of the caleosin/peroxygenase (CLO/PXG) gene/protein families in the Viridiplantae

Bioinformatics analyses of caleosin/peroxygenases (CLO/PXG) demonstrated that these genes are present in the vast majority of Viridiplantae taxa for which sequence data are available. Functionally active CLO/PXG proteins with roles in abiotic stress tolerance and lipid droplet storage are present in some Trebouxiophycean and Chlorophycean green algae but are absent from the small number of sequenced Prasinophyceaen genomes. CLO/PXG-like genes are expressed during dehydration stress in Charophyte algae, a sister clade of the land plants (Embryophyta). CLO/PXG-like sequences are also present in all of the >300 sequenced Embryophyte genomes, where some species contain as many as 10–12 genes that have arisen via selective gene duplication. Angiosperm genomes harbour at least one copy each of two distinct CLO/PX isoforms, termed H (high) and L (low), where H-forms contain an additional C-terminal motif of about 30–50 residues that is absent from L-forms. In contrast, species in other Viridiplantae taxa, including green algae, non-vascular plants, ferns and gymnosperms, contain only one (or occasionally both) of these isoforms per genome. Transcriptome and biochemical data show that CLO/PXG-like genes have complex patterns of developmental and tissue-specific expression. CLO/PXG proteins can associate with cytosolic lipid droplets and/or bilayer membranes. Many of the analysed isoforms also have peroxygenase activity and are involved in oxylipin metabolism. The distribution of CLO/PXG-like genes is consistent with an origin >1 billion years ago in at least two of the earliest diverging groups of the Viridiplantae, namely the Chlorophyta and the Streptophyta, after the Viridiplantae had already diverged from other Archaeplastidal groups such as the Rhodophyta and Glaucophyta. While algal CLO/PXGs have roles in lipid packaging and stress responses, the Embryophyte proteins have a much wider spectrum of roles and may have been instrumental in the colonisation of terrestrial habitats and the subsequent diversification as the major land flora.

The Antarctic sea ice alga Chlamydomonas sp. ICE-L provides insights into adaptive patterns of chloroplast evolution

BACKGROUND

The ice alga Chlamydomonas sp. ICE-L is the main contributor to primary productivity in Antarctic sea ice ecosystems and is well adapted to the extremely harsh environment. However, the adaptive mechanism of Chlamydomonas sp. ICE-L to sea-ice environment remains unclear. To study the adaptive strategies in Chlamydomonas sp. ICE-L, we investigated the molecular evolution of chloroplast photosynthetic genes that are essential for the accumulation of carbohydrate and energy living in Antarctic sea ice.

RESULTS

The 60 chloroplast protein-coding genes of Chlamydomonas sp. ICE-L were obtained, and the branch-site test detected significant signatures of positive selection on atpB, psaB, and rbcL genes in Chlamydomonas sp. ICE-L associated with the photosynthetic machinery. These positively selected genes were further identified as being under convergent evolution between Chlamydomonas sp. ICE-L and the halotolerant alga Dunaliella salina.

CONCLUSIONS

Our study provides evidence that the phototrophic component of Chlamydomonas sp. ICE-L exhibits adaptive evolution under extreme environment. The positive Darwinian selection operates on the chloroplast protein-coding genes of Antarctic ice algae adapted to extreme environment following functional-specific and lineages-specific patterns. In addition, three positively selected genes with convergent substitutions in Chlamydomonas sp. ICE-L were identified, and the adaptive modifications in these genes were in functionally important regions of the proteins. Our study provides a foundation for future experiments on the biochemical and physiological impacts of photosynthetic genes in green algae.

The Physcomitrella patens chromosome-scale assembly reveals moss genome structure and evolution

The draft genome of the moss model, Physcomitrella patens, comprised approximately 2000 unordered scaffolds. In order to enable analyses of genome structure and evolution we generated a chromosome-scale genome assembly using genetic linkage as well as (end) sequencing of long DNA fragments. We find that 57% of the genome comprises transposable elements (TEs), some of which may be actively transposing during the life cycle. Unlike in flowering plant genomes, gene- and TE-rich regions show an overall even distribution along the chromosomes. However, the chromosomes are mono-centric with peaks of a class of Copia elements potentially coinciding with centromeres. Gene body methylation is evident in 5.7% of the protein-coding genes, typically coinciding with low GC and low expression. Some giant virus insertions are transcriptionally active and might protect gametes from viral infection via siRNA mediated silencing. Structure-based detection methods show that the genome evolved via two rounds of whole genome duplications (WGDs), apparently common in mosses but not in liverworts and hornworts. Several hundred genes are present in colinear regions conserved since the last common ancestor of plants. These syntenic regions are enriched for functions related to plant-specific cell growth and tissue organization. The P. patens genome lacks the TE-rich pericentromeric and gene-rich distal regions typical for most flowering plant genomes. More non-seed plant genomes are needed to unravel how plant genomes evolve, and to understand whether the P. patens genome structure is typical for mosses or bryophytes.

The Glaciozyma antarctica genome reveals an array of systems that provide sustained responses towards temperature variations in a persistently cold habitat

Extremely low temperatures present various challenges to life that include ice formation and effects on metabolic capacity. Psyhcrophilic microorganisms typically have an array of mechanisms to enable survival in cold temperatures. In this study, we sequenced and analysed the genome of a psychrophilic yeast isolated in the Antarctic region, Glaciozyma antarctica. The genome annotation identified 7857 protein coding sequences. From the genome sequence analysis we were able to identify genes that encoded for proteins known to be associated with cold survival, in addition to annotating genes that are unique to G. antarctica. For genes that are known to be involved in cold adaptation such as anti-freeze proteins (AFPs), our gene expression analysis revealed that they were differentially transcribed over time and in response to different temperatures. This indicated the presence of an array of adaptation systems that can respond to a changing but persistent cold environment. We were also able to validate the activity of all the AFPs annotated where the recombinant AFPs demonstrated anti-freeze capacity. This work is an important foundation for further collective exploration into psychrophilic microbiology where among other potential, the genes unique to this species may represent a pool of novel mechanisms for cold survival.

Localization and dimer stability of a newly identified microbial rhodopsin from a polar, non-motile green algae

OBJECTIVE

The eukaryotic plasma membrane localized light-gated proton-pumping rhodopsins possesses great optogenetic applications for repolarization (silencing) of the neuronal activity simply by light illumination. Very few plasma membrane localized proton-pumping rhodopsins of a eukaryotic origin are known that have optogenetic potential. Our objective was to identify and characterize microbial rhodopsin of an eukaryotic origin that expresses on plasma membrane. The plasma membrane localized light-gated proton pump of an eukaryotic origin hold great promise to be used as an optogenetic tools for the neurobiology.

RESULTS

Here, we had characterized the cellular expression and membrane localization of a new rhodopsin in Antarctican algae Coccomyxa subellipsoidea. It is the first algal ion pumping rhodopsin that localizes to the plasma membrane of the eukaryotic cells. Coccomyxa subellipsoidea rhodopsin exists in the monomeric and dimeric state both the in vivo and in vitro. The dimeric form of the Coccomyxa subellipsoidea rhodopsin is resistant to heat and detergent denaturants.

Genome-wide analysis of carbohydrate-active enzymes in Pyramimonas parkeae (Prasinophyceae)

The wall-less green flagellate Pyramimonas parkeae is classified in clade I of the prasinophytes, a paraphyletic assemblage representing the last common ancestor of Viridiplantae, a monophyletic group composed of the green algae and land plants. Consequently, P. parkeae and other prasinophytes illuminate early-evolved Viridiplantae traits likely fundamental in the systems biology of green algae and land plants. Cellular structure and organellar genomes of P. parkeae are now well understood, and transcriptomic sequence data are also publically available for one strain of this species, but corresponding nuclear genomic sequence data are lacking. For this reason, we obtained shotgun genomic sequence and assembled a draft nuclear genome for P. parkeaeNIES254 to use along with existing transcriptomic sequence to focus on carbohydrate-active enzymes. We found that the P. parkeae nuclear genome encodes carbohydrate-active protein families similar to those previously observed for other prasinophytes, green algae, and early-diverging embryophytes for which full nuclear genomic sequence is publically available. Sequences homologous to genes related to biosynthesis of starch and cell wall carbohydrates were identified in the P. parkeae genome, indicating molecular traits common to Viridiplantae. For example, the P. parkeae genome includes sequences clustering with bacterial genes that encode cellulose synthases (Bcs), including regions coding for domains common to bacterial and plant cellulose synthases; these new sequences were incorporated into phylogenies aimed at illuminating the evolutionary history of cellulose production by Viridiplantae. Genomic sequences related to biosynthesis of xyloglucans, pectin, and starch likewise shed light on the origin of key Viridiplantae traits.

Acidophilic green algal genome provides insights into adaptation to an acidic environment

Some microalgae are adapted to extremely acidic environments in which toxic metals are present at high levels. However, little is known about how acidophilic algae evolved from their respective neutrophilic ancestors by adapting to particular acidic environments. To gain insights into this issue, we determined the draft genome sequence of the acidophilic green alga Chlamydomonas eustigma and performed comparative genome and transcriptome analyses between Ceustigma and its neutrophilic relative Chlamydomonas reinhardtii The results revealed the following features in Ceustigma that probably contributed to the adaptation to an acidic environment. Genes encoding heat-shock proteins and plasma membrane H⁺-ATPase are highly expressed in Ceustigma This species has also lost fermentation pathways that acidify the cytosol and has acquired an energy shuttle and buffering system and arsenic detoxification genes through horizontal gene transfer. Moreover, the arsenic detoxification genes have been multiplied in the genome. These features have also been found in other acidophilic green and red algae, suggesting the existence of common mechanisms in the adaptation to acidic environments.

Transcriptome sequencing and characterization of Symbiodinium muscatinei and Elliptochloris marina, symbionts found within the aggregating sea anemone Anthopleura elegantissima

There is a growing body of literature using transcriptomic data to study how tropical cnidarians and their photosynthetic endosymbionts respond to environmental stressors and participate in metabolic exchange. Despite these efforts, our understanding of how essential genes function to facilitate symbiosis establishment and maintenance remains limited. The inclusion of taxonomically and ecologically diverse endosymbionts will enhance our understanding of these interactions. Here we characterize the transcriptomes of two very different symbionts found within the temperate sea anemone Anthopleura elegantissima: the chlorophyte Elliptochloris marina and the dinoflagellate Symbiodinium muscatinei. We use a multi-level approach to assess the diversity of genes found across S. muscatinei and E. marina transcriptomes, and compare their overall protein domains with other dinoflagellates and chlorophytes. Our analysis identified several genes that are potentially involved in mitigating stress response (e.g., heat shock proteins pathways for mediating reactive oxygen species) and metabolic exchange (e.g., ion transporters). Finally, we show that S. muscatinei and other Symbiodinium strains are equipped with a high salt peridinin-chl-protein (HSPCP) gene previously identified only in free-living dinoflagellates. The addition of these transcriptomes to the cnidarian-symbiont molecular toolkit will aid in understanding how these vitally important symbiotic relationships are established and maintained across a variety of environmental conditions.

Multiple origins of endosymbionts in Chlorellaceae with no reductive effects on the plastid or mitochondrial genomes

Ancient endosymbiotic relationships have led to extreme genomic reduction in many bacterial and eukaryotic algal endosymbionts. Endosymbionts in more recent and/or facultative relationships can also experience genomic reduction to a lesser extent, but little is known about the effects of the endosymbiotic transition on the organellar genomes of eukaryotes. To understand how the endosymbiotic lifestyle has affected the organellar genomes of photosynthetic green algae, we generated the complete plastid genome (plastome) and mitochondrial genome (mitogenome) sequences from three green algal endosymbionts (Chlorella heliozoae, Chlorella variabilis and Micractinium conductrix). The mitogenomes and plastomes of the three newly sequenced endosymbionts have a standard set of genes compared with free-living trebouxiophytes, providing no evidence for functional genomic reduction. Instead, their organellar genomes are generally larger and more intron rich. Intron content is highly variable among the members of Chlorella, suggesting very high rates of gain and/or loss of introns during evolution. Phylogenetic analysis of plastid and mitochondrial genes demonstrated that the three endosymbionts do not form a monophyletic group, indicating that the endosymbiotic lifestyle has evolved multiple times in Chlorellaceae. In addition, M. conductrix is deeply nested within the Chlorella clade, suggesting that taxonomic revision is needed for one or both genera.

Molecular cloning and functional characterization of NADPH-dependent cytochrome P450 reductase from the green microalga Botryococcus braunii, B race

The green microalga Botryococcus braunii of the B race accumulates various lipophilic compounds containing a 10,11-oxidosqualene epoxide moiety in addition to large amounts of triterpene hydrocarbons. While 2,3-squalene epoxidases have already been isolated and characterized from the alga, the enzyme that catalyzes the 10,11-epoxidation of squalene has remained elusive. In order to obtain a molecular tool to explore a 10,11-squalene epoxidase, cDNA cloning of an NADPH-dependent cytochrome P450 reductase (CPR) that is required by both squalene epoxidases and cytochrome P450 enzymes was carried out. The isolated cDNA contained an open reading frame (1998 bp) that encoded for a protein with 665 amino acid residues with a predicted molecular weight of 71.46 kDa and a theoretical pI of 5.49. Analysis of the deduced amino acid sequence revealed the presence of conserved motifs, including FMN, FAD, and NADPH binding domains, which are typical of other CPRs and necessary for enzyme activity. By truncation of the N-terminal transmembrane anchor and addition of a 6× His-tag, BbCPR was heterologously produced in Escherichia coli and purified by Ni-NTA affinity chromatography. The purified recombinant enzyme showed optimal reducing activity of cytochrome c at around a neutral pH at a temperature range of 30-37°C. For steady state kinetic parameters, the recombinant enzyme had a k_m for cytochrome c and NADPH of 11.7±1.6 and 9.4±1.4 μM, and a k_cat for cytochrome c and NADPH of 2.78±0.09 and 3.66±0.11 μmol/min/mg protein, respectively. This is the first study to perform the functional characterization of a CPR from eukaryotic microalgae.

The genome and phenome of the green alga Chloroidium sp. UTEX 3007 reveal adaptive traits for desert acclimatization

To investigate the phenomic and genomic traits that allow green algae to survive in deserts, we characterized a ubiquitous species, Chloroidium sp. UTEX 3007, which we isolated from multiple locations in the United Arab Emirates (UAE). Metabolomic analyses of Chloroidium sp. UTEX 3007 indicated that the alga accumulates a broad range of carbon sources, including several desiccation tolerance-promoting sugars and unusually large stores of palmitate. Growth assays revealed capacities to grow in salinities from zero to 60 g/L and to grow heterotrophically on >40 distinct carbon sources. Assembly and annotation of genomic reads yielded a 52.5 Mbp genome with 8153 functionally annotated genes. Comparison with other sequenced green algae revealed unique protein families involved in osmotic stress tolerance and saccharide metabolism that support phenomic studies. Our results reveal the robust and flexible biology utilized by a green alga to successfully inhabit a desert coastline.

DOI:http://dx.doi.org/10.7554/eLife.25783.001

Single-celled green algae, also known as green microalgae, play an important role for the world’s ecosystems, in part, because they can harness energy from sunlight to produce carbon-rich compounds. Microalgae are also important for biotechnology and people have harnessed them to make food, fuel and medicines. Green microalgae live in many types of habitats from streams to oceans, and they can also be found on the land, including in deserts. Like plants that live in the desert, these microalgae have likely evolved specific traits that allow them to live in these hot and dry regions. Yet, fewer scientists have studied microalgae compared to land plants, and until now it was not well understood how microalgae could survive in the desert.

Nelson et al. analyzed green microalgae from different locations around the United Arab Emirates and found that one microalga, known as Chloroidium, is one of the most dominant algae in this area. This included samples from beaches, mangroves, desert oases, buildings and public fresh water sources. Chloroidium has a unique set of genes and proteins and grew particularly well in freshwater and saltwater. Rather than just harnessing sunlight, the microalgae were able to consume over 40 different varieties of carbon sources to produce energy. The microalgae also accumulated oily molecules with a similar composition to palm oil, which may help this species to survive in desert regions.

A next step will be to develop biotechnological assets based on the information obtained. In the future, microalgae could be used to make an oil that represents an alternative to palm oil; this would reduce the demand for palm tree plantations, which pose a major threat to the natural environment. Moreover, understanding how microalgae can colonize a desert region will help us to understand the effects of climate change in the region.

DOI:http://dx.doi.org/10.7554/eLife.25783.002

Chromosome-level genome assembly and transcriptome of the green alga Chromochloris zofingiensis illuminates astaxanthin production

Microalgae have potential to help meet energy and food demands without exacerbating environmental problems. There is interest in the unicellular green alga Chromochloris zofingiensis, because it produces lipids for biofuels and a highly valuable carotenoid nutraceutical, astaxanthin. To advance understanding of its biology and facilitate commercial development, we present a C. zofingiensis chromosome-level nuclear genome, organelle genomes, and transcriptome from diverse growth conditions. The assembly, derived from a combination of short- and long-read sequencing in conjunction with optical mapping, revealed a compact genome of ∼58 Mbp distributed over 19 chromosomes containing 15,274 predicted protein-coding genes. The genome has uniform gene density over chromosomes, low repetitive sequence content (∼6%), and a high fraction of protein-coding sequence (∼39%) with relatively long coding exons and few coding introns. Functional annotation of gene models identified orthologous families for the majority (∼73%) of genes. Synteny analysis uncovered localized but scrambled blocks of genes in putative orthologous relationships with other green algae. Two genes encoding beta-ketolase (BKT), the key enzyme synthesizing astaxanthin, were found in the genome, and both were up-regulated by high light. Isolation and molecular analysis of astaxanthin-deficient mutants showed that BKT1 is required for the production of astaxanthin. Moreover, the transcriptome under high light exposure revealed candidate genes that could be involved in critical yet missing steps of astaxanthin biosynthesis, including ABC transporters, cytochrome P450 enzymes, and an acyltransferase. The high-quality genome and transcriptome provide insight into the green algal lineage and carotenoid production.

γ-Aminobutyric acid (GABA) signalling in plants

The role of γ-aminobutyric acid (GABA) as a signal in animals has been documented for over 60 years. In contrast, evidence that GABA is a signal in plants has only emerged in the last 15 years, and it was not until last year that a mechanism by which this could occur was identified-a plant 'GABA receptor' that inhibits anion passage through the aluminium-activated malate transporter family of proteins (ALMTs). ALMTs are multigenic, expressed in different organs and present on different membranes. We propose GABA regulation of ALMT activity could function as a signal that modulates plant growth, development, and stress response. In this review, we compare and contrast the plant 'GABA receptor' with mammalian GABAA receptors in terms of their molecular identity, predicted topology, mode of action, and signalling roles. We also explore the implications of the discovery that GABA modulates anion flux in plants, its role in signal transduction for the regulation of plant physiology, and predict the possibility that there are other GABA interaction sites in the N termini of ALMT proteins through in silico evolutionary coupling analysis; we also explore the potential interactions between GABA and other signalling molecules.

RNAi-based targeted gene knockdown in the model oleaginous microalgae Nannochloropsis oceanica

Microalgae are promising feedstock for renewable fuels such as biodiesel, yet development of industrial oleaginous strains has been hindered by the paucity and inefficiency of reverse genetics tools. Here we established an efficient RNAi-based targeted gene-knockdown method for Nannochloropsis spp., which are emerging model organisms for industrial microalgal oil production. The method achieved a 40-80% success rate in Nannochloropsis oceanica strain IMET1. When transcript level of one carbonic anhydrase (CA) was inhibited by 62-83% via RNAi, mutant cells exhibited photosynthetic oxygen evolution (POE) rates that were 68-100% higher than wild-type (WT) at pH 6.0, equivalent to WT at pH 8.2, yet 39-45% lower than WT at pH 9.0. Moreover, the mutant POE rates were negatively correlated with the increase of culture pH, an exact opposite of WT. Thus, a dynamic carbon concentration mechanism (CCM) that is highly sensitive to pH homeostasis was revealed, where the CA inhibition likely partially abrogated the mechanism that normally deactivates CCM under a high level of dissolved CO₂ . Extension of the method to another sequenced N. oceanica strain of CCMP 1779 demonstrated comparable performance. Finally, McrBC-PCR followed by bisulfite sequencing revealed that the gene knockdown is mediated by the CG, CHG and CHH types of DNA methylation at the coding region of the targeted gene. The efficiency, robustness and general applicability of this reverse genetics approach suggested the possibility of large-scale RNAi-based gene function screening in industrial microalgae.

Exploring the transcriptome of non-model oleaginous microalga Dunaliella tertiolecta through high-throughput sequencing and high performance computing

Background

RNA-Seq technology has received a lot of attention in recent years for microalgal global transcriptomic profiling. It is widely used in transcriptome-wide analysis of gene expression., particularly for microalgal strains with potential as biofuel sources. However, insufficient genomic or transcriptomic information of non-model microalgae has limited the understanding of their regulatory mechanisms and hampered genetic manipulation to enhance biofuel production. As such, an optimal microalgal transcriptomic database construction is a subject of urgent investigation.

Results

Dunaliella tertiolecta, a non-model oleaginous microalgal species, was sequenced via Illumina MISEQ and HISEQ 4000 in RNA-Seq studies. The high quality high-throughout sequencing data were explored using high performance computing (HPC) in a petascale data center and subjected to de novo assembly and parallelized mpiBLASTX search with multiple species. As a result, a transcriptome database of 17,845 was constructed (~95% completeness). This enlarged database constructed fueled the RNA-Seq data analysis, which was validated by a nitrogen deprivation (ND) study that induces triacylglycerol (TAG) production.

Conclusions

The new paralleled assembly and annotation method under HPC presented here allows the solution of large-scale data processing problems in acceptable computation time. There is significant increase in the number of transcriptomic data achieved and observable heterogeneity in the performance to identify differentially expressed genes in the ND treatment paradigm. The results provide new insights as to how response to ND treatment in microalgae is regulated. ND analyses highlight the advantages of this database generated in this study that could also serve as a useful resource for future gene manipulation and transcriptome-wide analysis. We thus demonstrate the usefulness of exploring the transcriptome as an informative platform for functional studies and genetic manipulations in similar species.

Electronic supplementary material

The online version of this article (doi:10.1186/s12859-017-1551-x) contains supplementary material, which is available to authorized users.

Horizontal Gene Transfer of Phytochelatin Synthases from Bacteria to Extremophilic Green Algae

Transcriptomic sequencing together with bioinformatic analyses and an automated annotation process led us to identify novel phytochelatin synthase (PCS) genes from two extremophilic green algae (Chlamydomonas acidophila and Dunaliella acidophila). These genes are of intermediate length compared to known PCS genes from eukaryotes and PCS-like genes from prokaryotes. A detailed phylogenetic analysis gives new insight into the complicated evolutionary history of PCS genes and provides evidence for multiple horizontal gene transfer events from bacteria to eukaryotes within the gene family. A separate subgroup containing PCS-like genes within the PCS gene family is not supported since the PCS genes are monophyletic only when the PCS-like genes are included. The presence and functionality of the novel genes in the organisms were verified by genomic sequencing and qRT-PCR. Furthermore, the novel PCS gene in Chlamydomonas acidophila showed very strong induction by cadmium. Cloning and expression of the gene in Escherichia coli clearly improves its cadmium resistance. The gene in Dunaliella was not induced, most likely due to gene duplication.

Phosphoribosyl Diphosphate (PRPP): Biosynthesis, Enzymology, Utilization, and Metabolic Significance

Phosphoribosyl diphosphate (PRPP) is an important intermediate in cellular metabolism. PRPP is synthesized by PRPP synthase, as follows: ribose 5-phosphate + ATP → PRPP + AMP. PRPP is ubiquitously found in living organisms and is used in substitution reactions with the formation of glycosidic bonds. PRPP is utilized in the biosynthesis of purine and pyrimidine nucleotides, the amino acids histidine and tryptophan, the cofactors NAD and tetrahydromethanopterin, arabinosyl monophosphodecaprenol, and certain aminoglycoside antibiotics. The participation of PRPP in each of these metabolic pathways is reviewed. Central to the metabolism of PRPP is PRPP synthase, which has been studied from all kingdoms of life by classical mechanistic procedures. The results of these analyses are unified with recent progress in molecular enzymology and the elucidation of the three-dimensional structures of PRPP synthases from eubacteria, archaea, and humans. The structures and mechanisms of catalysis of the five diphosphoryltransferases are compared, as are those of selected enzymes of diphosphoryl transfer, phosphoryl transfer, and nucleotidyl transfer reactions. PRPP is used as a substrate by a large number phosphoribosyltransferases. The protein structures and reaction mechanisms of these phosphoribosyltransferases vary and demonstrate the versatility of PRPP as an intermediate in cellular physiology. PRPP synthases appear to have originated from a phosphoribosyltransferase during evolution, as demonstrated by phylogenetic analysis. PRPP, furthermore, is an effector molecule of purine and pyrimidine nucleotide biosynthesis, either by binding to PurR or PyrR regulatory proteins or as an allosteric activator of carbamoylphosphate synthetase. Genetic analyses have disclosed a number of mutants altered in the PRPP synthase-specifying genes in humans as well as bacterial species.

Carbon and Acyl Chain Flux during Stress-induced Triglyceride Accumulation by Stable Isotopic Labeling of the Polar Microalga Coccomyxa subellipsoidea C169

Deriving biofuels and other lipoid products from algae is a promising future technology directly addressing global issues of atmospheric CO₂ balance. To better understand the metabolism of triglyceride synthesis in algae, we examined their metabolic origins in the model species, Coccomyxa subellipsoidea C169, using stable isotopic labeling. Labeling patterns arising from [U-¹³C]glucose, ¹³CO₂, or D₂O supplementation were analyzed by GC-MS and/or LC-MS over time courses during nitrogen starvation to address the roles of catabolic carbon recycling, acyl chain redistribution, and de novo fatty acid (FA) synthesis during the expansion of the lipid bodies. The metabolic origin of stress-induced triglyceride was found to be a continuous 8:2 ratio between de novo synthesized FA and acyl chain transfer from pre-stressed membrane lipids with little input from lipid remodeling. Membrane lipids were continually synthesized with associated acyl chain editing during nitrogen stress, in contrast to an overall decrease in total membrane lipid. The incorporation rates of de novo synthesized FA into lipid classes were measured over a time course of nitrogen starvation. The synthesis of triglycerides, phospholipids, and galactolipids followed a two-stage pattern where nitrogen starvation resulted in a 2.5-fold increase followed by a gradual decline. Acyl chain flux into membrane lipids was dominant in the first stage followed by triglycerides. These data indicate that the level of metabolic control that determines acyl chain flux between membrane lipids and triglycerides during nitrogen stress relies primarily on the Kennedy pathway and de novo FA synthesis with limited, defined input from acyl editing reactions.

Conserved Roles of CrRLK1L Receptor-Like Kinases in Cell Expansion and Reproduction from Algae to Angiosperms

Receptor-like kinases (RLKs) are regulators of plant development through allowing cells to sense their extracellular environment. They facilitate detection of local endogenous signals, in addition to external biotic and abiotic stimuli. The Catharanthus roseus RLK1-like (CrRLK1L) protein kinase subfamily, which contains FERONIA, plays a central role in regulating fertilization and in cell expansion mechanisms such as cell elongation and tip growth, as well as having indirect links to plant-pathogen interactions. Several components of CrRLK1L signaling pathways have been identified, including an extracellular ligand, coreceptors, and downstream signaling elements. The presence and abundance of the CrRLK1L proteins in the plant kingdom suggest an origin within the Streptophyta lineage, with a notable increase in prevalence in the seeded land plants. Given the function of the sole CrRLK1L protein in a charophycean alga, the possibility of a conserved role in detection and/or regulation of cell wall integrity throughout the Strephtophytes is discussed. Orthologs of signaling pathway components are also present in extant representatives of non-vascular land plants and early vascular land plants including the liverwort Marchantia polymorpha, the moss Physcomitrella patens and the lycophyte Selaginella moellendorffii. Deciphering the roles in development of the CrRLK1L protein kinases in early diverging land plants will provide insights into their ancestral function, furthering our understanding of this diversified subfamily of receptors in higher plants.

The genome and transcriptome of Trichormus sp. NMC-1: insights into adaptation to extreme environments on the Qinghai-Tibet Plateau

The Qinghai-Tibet Plateau (QTP) has the highest biodiversity for an extreme environment worldwide, and provides an ideal natural laboratory to study adaptive evolution. In this study, we generated a draft genome sequence of cyanobacteria Trichormus sp. NMC-1 in the QTP and performed whole transcriptome sequencing under low temperature to investigate the genetic mechanism by which T. sp. NMC-1 adapted to the specific environment. Its genome sequence was 5.9 Mb with a G+C content of 39.2% and encompassed a total of 5362 CDS. A phylogenomic tree indicated that this strain belongs to the Trichormus and Anabaena cluster. Genome comparison between T. sp. NMC-1 and six relatives showed that functionally unknown genes occupied a much higher proportion (28.12%) of the T. sp. NMC-1 genome. In addition, functions of specific, significant positively selected, expanded orthogroups, and differentially expressed genes involved in signal transduction, cell wall/membrane biogenesis, secondary metabolite biosynthesis, and energy production and conversion were analyzed to elucidate specific adaptation traits. Further analyses showed that the CheY-like genes, extracellular polysaccharide and mycosporine-like amino acids might play major roles in adaptation to harsh environments. Our findings indicate that sophisticated genetic mechanisms are involved in cyanobacterial adaptation to the extreme environment of the QTP.

New features of desiccation tolerance in the lichen photobiont Trebouxia gelatinosa are revealed by a transcriptomic approach

Trebouxia is the most common lichen-forming genus of aero-terrestrial green algae and all its species are desiccation tolerant (DT). The molecular bases of this remarkable adaptation are, however, still largely unknown. We applied a transcriptomic approach to a common member of the genus, T. gelatinosa, to investigate the alteration of gene expression occurring after dehydration and subsequent rehydration in comparison to cells kept constantly hydrated. We sequenced, de novo assembled and annotated the transcriptome of axenically cultured T. gelatinosa by using Illumina sequencing technology. We tracked the expression profiles of over 13,000 protein-coding transcripts. During the dehydration/rehydration cycle c. 92 % of the total protein-coding transcripts displayed a stable expression, suggesting that the desiccation tolerance of T. gelatinosa mostly relies on constitutive mechanisms. Dehydration and rehydration affected mainly the gene expression for components of the photosynthetic apparatus, the ROS-scavenging system, Heat Shock Proteins, aquaporins, expansins, and desiccation related proteins (DRPs), which are highly diversified in T. gelatinosa, whereas Late Embryogenesis Abundant Proteins were not affected. Only some of these phenomena were previously observed in other DT green algae, bryophytes and resurrection plants, other traits being distinctive of T. gelatinosa, and perhaps related to its symbiotic lifestyle. Finally, the phylogenetic inference extended to DRPs of other chlorophytes, embryophytes and bacteria clearly pointed out that DRPs of chlorophytes are not orthologous to those of embryophytes: some of them were likely acquired through horizontal gene transfer from extremophile bacteria which live in symbiosis within the lichen thallus.

Comparative genomics explains the evolutionary success of reef-forming corals

Transcriptome and genome data from twenty stony coral species and a selection of reference bilaterians were studied to elucidate coral evolutionary history. We identified genes that encode the proteins responsible for the precipitation and aggregation of the aragonite skeleton on which the organisms live, and revealed a network of environmental sensors that coordinate responses of the host animals to temperature, light, and pH. Furthermore, we describe a variety of stress-related pathways, including apoptotic pathways that allow the host animals to detoxify reactive oxygen and nitrogen species that are generated by their intracellular photosynthetic symbionts, and determine the fate of corals under environmental stress. Some of these genes arose through horizontal gene transfer and comprise at least 0.2% of the animal gene inventory. Our analysis elucidates the evolutionary strategies that have allowed symbiotic corals to adapt and thrive for hundreds of millions of years.

DOI:http://dx.doi.org/10.7554/eLife.13288.001

For millions of years, reef-building stony corals have created extensive habitats for numerous marine plants and animals in shallow tropical seas. Stony corals consist of many small, tentacled animals called polyps. These polyps secrete a mineral called aragonite to create the reef – an external ‘skeleton’ that supports and protects the corals.

Photosynthesizing algae live inside the cells of stony corals, and each species depends on the other to survive. The algae produce the coral’s main source of food, although they also produce some waste products that can harm the coral if they build up inside cells. If the oceans become warmer and more acidic, the coral are more likely to become stressed and expel the algae from their cells in a process known as coral bleaching. This makes the coral more likely to die or become diseased. Corals have survived previous periods of ocean warming, although it is not known how they evolved to do so.

The evolutionary history of an organism can be traced by studying its genome – its complete set of DNA – and the RNA molecules encoded by these genes. Bhattacharya et al. performed this analysis for twenty stony coral species, and compared the resulting genome and RNA sequences with the genomes of other related marine organisms, such as sea anemones and sponges. In particular, Bhattacharya et al. examined “ortholog” groups of genes, which are present in different species and evolved from a common ancestral gene. This analysis identified the genes in the corals that encode the proteins responsible for constructing the aragonite skeleton. The coral genome also encodes a network of environmental sensors that coordinate how the polyps respond to temperature, light and acidity.

Bhattacharya et al. also uncovered a variety of stress-related pathways, including those that detoxify the polyps of the damaging molecules generated by algae, and the pathways that enable the polyps to adapt to environmental stress. Many of these genes were recruited from other species in a process known as horizontal gene transfer.

The oceans are expected to become warmer and more acidic in the coming centuries. Provided that humans do not physically destroy the corals’ habitats, the evidence found by Bhattacharya et al. suggests that the genome of the corals contains the diversity that will allow them to adapt to these new conditions.

DOI:http://dx.doi.org/10.7554/eLife.13288.002

The Gonium pectorale genome demonstrates co-option of cell cycle regulation during the evolution of multicellularity

The transition to multicellularity has occurred numerous times in all domains of life, yet its initial steps are poorly understood. The volvocine green algae are a tractable system for understanding the genetic basis of multicellularity including the initial formation of cooperative cell groups. Here we report the genome sequence of the undifferentiated colonial alga, Gonium pectorale, where group formation evolved by co-option of the retinoblastoma cell cycle regulatory pathway. Significantly, expression of the Gonium retinoblastoma cell cycle regulator in unicellular Chlamydomonas causes it to become colonial. The presence of these changes in undifferentiated Gonium indicates extensive group-level adaptation during the initial step in the evolution of multicellularity. These results emphasize an early and formative step in the evolution of multicellularity, the evolution of cell cycle regulation, one that may shed light on the evolutionary history of other multicellular innovations and evolutionary transitions.

Elucidating the composition and conservation of the autophagy pathway in photosynthetic eukaryotes

Aquatic photosynthetic eukaryotes represent highly diverse groups (green, red, and chromalveolate algae) derived from multiple endosymbiosis events, covering a wide spectrum of the tree of life. They are responsible for about 50% of the global photosynthesis and serve as the foundation for oceanic and fresh water food webs. Although the ecophysiology and molecular ecology of some algal species are extensively studied, some basic aspects of algal cell biology are still underexplored. The recent wealth of genomic resources from algae has opened new frontiers to decipher the role of cell signaling pathways and their function in an ecological and biotechnological context. Here, we took a bioinformatic approach to explore the distribution and conservation of TOR and autophagy-related (ATG) proteins (Atg in yeast) in diverse algal groups. Our genomic analysis demonstrates conservation of TOR and ATG proteins in green algae. In contrast, in all 5 available red algal genomes, we could not detect the sequences that encode for any of the 17 core ATG proteins examined, albeit TOR and its interacting proteins are conserved. This intriguing data suggests that the autophagy pathway is not conserved in red algae as it is in the entire eukaryote domain. In contrast, chromalveolates, despite being derived from the red-plastid lineage, retain and express ATG genes, which raises a fundamental question regarding the acquisition of ATG genes during algal evolution. Among chromalveolates, Emiliania huxleyi (Haptophyta), a bloom-forming coccolithophore, possesses the most complete set of ATG genes, and may serve as a model organism to study autophagy in marine protists with great ecological significance.

Characterization of an epoxide hydrolase from the Florida red tide dinoflagellate, Karenia brevis

Epoxide hydrolases (EH, EC 3.3.2.3) have been proposed to be key enzymes in the biosynthesis of polyether (PE) ladder compounds such as the brevetoxins which are produced by the dinoflagellate Karenia brevis. These enzymes have the potential to catalyze kinetically disfavored endo-tet cyclization reactions. Data mining of K. brevis transcriptome libraries revealed two classes of epoxide hydrolases: microsomal and leukotriene A4 (LTA4) hydrolases. A microsomal EH was cloned and expressed for characterization. The enzyme is a monomeric protein with molecular weight 44kDa. Kinetic parameters were evaluated using a variety of epoxide substrates to assess substrate selectivity and enantioselectivity, as well as its potential to catalyze the critical endo-tet cyclization of epoxy alcohols. Monitoring of EH activity in high and low toxin producing cultures of K. brevis over a three week period showed consistently higher activity in the high toxin producing culture implicating the involvement of one or more EH in brevetoxin biosynthesis.

Microalgae as a Source for VLC-PUFA Production

Microalgae present a huge and still insufficiently tapped resource of very long-chain omega-3 and omega-6 polyunsaturated fatty acids (VLC-PUFA) for human nutrition and medicinal applications. This chapter describes the diversity of unicellular eukaryotic microalgae in respect to VLC-PUFA biosynthesis. Then, we outline the major biosynthetic pathways mediating the formation of VLC-PUFA by sequential desaturation and elongation of C18-PUFA acyl groups. We address the aspects of spatial localization of those pathways and elaborate on the role for VLC-PUFA in microalgal cells. Recent progress in microalgal genetic transformation and molecular engineering has opened the way to increased production efficiencies for VLC-PUFA. The perspectives of photobiotechnology and metabolic engineering of microalgae for altered or enhanced VLC-PUFA production are also discussed.

Transcriptome analysis reveals global regulation in response to CO2 supplementation in oleaginous microalga Coccomyxa subellipsoidea C-169

BACKGROUND

Microalgae are emerging as suitable feedstock for renewable biofuel production and providing a promising way to alleviate green house gas CO2. Characterizing the metabolic pathways involved in the biosynthesis of energy-rich compounds and their global regulation upon elevated CO2 is necessary to explore the mechanism underlying rapid growth and lipid accumulation, so as to realize the full potential of these organisms as energy resources.

RESULTS

In the present study, 2 and 5 % CO2 increased growth rate and lipid accumulation in autotrophically cultured green alga Coccomyxa subellipsoidea C-169. Overall biomass productivity as 222 mg L(-1) day(-1) and fatty acid content as 48.5 % dry cell weight were attained in 2 % CO2, suggesting C-169 as a great candidate for lipid production via CO2 supplementation. Transcriptomic analysis of 2 % against 0.04 % CO2-cultured C-169 unveiled the global regulation of important metabolic processes. Other than enhancing gene expression in the Calvin cycle, C-169 upregulated the expression of phosphoenolpyruvate carboxylase, pyruvate carboxylase and carbamoyl-phosphate synthetase II to enhance the anaplerotic carbon assimilation reactions upon elevated CO2. Upregulation of ferredoxin and ferredoxin-NADP(+) reductase implied that plentiful energy captured through photosynthesis was transferred through ferredoxin to sustain rapid growth and lipid accumulation. Genes involved in the glycolysis, TCA cycle and oxidative phosphorylation were predominantly upregulated presumably to provide abundant intermediates and metabolic energy for anabolism. Coordinated upregulation of nitrogen acquisition and assimilation genes, together with activation of specific carbamoyl-phosphate synthetase and ornithine pathway genes, might help C-169 to maintain carbon/nitrogen balance upon elevated CO2. Significant downregulation of fatty acid degradation genes, as well as the upregulation of fatty acid synthesis genes at the later stage might contribute to the tremendous lipid accumulation.

CONCLUSION

Global and collaborative regulation was employed by C-169 to assimilate more carbon and maintain carbon/nitrogen balance upon elevated CO2, which provide abundant carbon skeleton and affluent metabolic energy to sustain rapid growth and lipid accumulation. Data here for the first time bring significant insights into the regulatory profile of metabolism and acclimation to elevated CO2 in C-169, which provide important information for future metabolic engineering in the development of sustainable microalgae-based biofuels.

Plastoquinone and Ubiquinone in Plants: Biosynthesis, Physiological Function and Metabolic Engineering

Plastoquinone (PQ) and ubiquinone (UQ) are two important prenylquinones, functioning as electron transporters in the electron transport chain of oxygenic photosynthesis and the aerobic respiratory chain, respectively, and play indispensable roles in plant growth and development through participating in the biosynthesis and metabolism of important chemical compounds, acting as antioxidants, being involved in plant response to stress, and regulating gene expression and cell signal transduction. UQ, particularly UQ₁₀, has also been widely used in people's life. It is effective in treating cardiovascular diseases, chronic gingivitis and periodontitis, and shows favorable impact on cancer treatment and human reproductive health. PQ and UQ are made up of an active benzoquinone ring attached to a polyisoprenoid side chain. Biosynthesis of PQ and UQ is very complicated with more than thirty five enzymes involved. Their synthetic pathways can be generally divided into two stages. The first stage leads to the biosynthesis of precursors of benzene quinone ring and prenyl side chain. The benzene quinone ring for UQ is synthesized from tyrosine or phenylalanine, whereas the ring for PQ is derived from tyrosine. The prenyl side chains of PQ and UQ are derived from glyceraldehyde 3-phosphate and pyruvate through the 2-C-methyl-D-erythritol 4-phosphate pathway and/or acetyl-CoA and acetoacetyl-CoA through the mevalonate pathway. The second stage includes the condensation of ring and side chain and subsequent modification. Homogentisate solanesyltransferase, 4-hydroxybenzoate polyprenyl diphosphate transferase and a series of benzene quinone ring modification enzymes are involved in this stage. PQ exists in plants, while UQ widely presents in plants, animals and microbes. Many enzymes and their encoding genes involved in PQ and UQ biosynthesis have been intensively studied recently. Metabolic engineering of UQ₁₀ in plants, such as rice and tobacco, has also been tested. In this review, we summarize and discuss recent research progresses in the biosynthetic pathways of PQ and UQ and enzymes and their encoding genes involved in side chain elongation and in the second stage of PQ and UQ biosynthesis. Physiological functions of PQ and UQ played in plants as well as the practical application and metabolic engineering of PQ and UQ are also included.

Genomic analysis of Luteimonas abyssi XH031(T): insights into its adaption to the subseafloor environment of South Pacific Gyre and ecological role in biogeochemical cycle

Background

Luteimonas abyssi XH031^T, which was previously isolated from subseafloor environment of the South Pacific Gyre (SPG), was an aerobic, gram-negative bacterium, and was identified to be a novel species of the genus Luteimonas in the family of Xanthomonadaceae. The nutrients utilization and metabolic mechanisms of XH031^T indicate its plasticity. In view of the above characteristics, its genome was sequenced, and an in-depth analysis of the XH031^T genome was performed to elucidate its adaption to extreme ecological environment.

Results

Various macromolecules including polysaccharide, protein, lipid and DNA could be degraded at low temperature by XH031^T under laboratory conditions, and its degradation abilities to starch, gelatin and casein were considerably strong. Genome sequence analysis indicated that XH031^T possesses extensive enzyme-encoding genes compared with four other Luteimonas strains. In addition, intricate systems (such as two-component regulatory systems, secretion systems, etc.), which are often used by bacteria to modulate the interactions of bacteria with their environments, were predicted in the genome of XH031^T. Genes encoding a choline-glycine betaine transporter and 99 extracellular peptidases featured with halophilicity were predicted in the genome, which might help the bacterium to adapt to the salty marine environment. Moreover, there were many gene clusters in the genome encoding ATP-binding cassette superfamily transporters, major facilitator superfamily transporters and cytochrome P450s that might function in the process of various substrate transportation and metabolisms. Furthermore, drug resistance genes harbored in the genome might signify that XH031^T has evolved hereditary adaptation to toxic environment. Finally, the annotation of metabolic pathways of the elements (such as carbon, nitrogen, sulfur, phosphor and iron) in the genome elucidated the degradation of organic matter in the deep sediment of the SPG.

Conclusions

The genome analysis showed that XH031^T had genetic advantages to adapt to subseafloor environment. The material metabolism manifests that the strain may play an important ecological role in the biogeochemical cycle of the SPG, and various cold-adapted extracelluar enzymes produced by the strain may have significant value in application.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-2326-2) contains supplementary material, which is available to authorized users.

Reverse transcriptase genes are highly abundant and transcriptionally active in marine plankton assemblages

Genes encoding reverse transcriptases (RTs) are found in most eukaryotes, often as a component of retrotransposons, as well as in retroviruses and in prokaryotic retroelements. We investigated the abundance, classification and transcriptional status of RTs based on Tara Oceans marine metagenomes and metatranscriptomes encompassing a wide organism size range. Our analyses revealed that RTs predominate large-size fraction metagenomes (>5 μm), where they reached a maximum of 13.5% of the total gene abundance. Metagenomic RTs were widely distributed across the phylogeny of known RTs, but many belonged to previously uncharacterized clades. Metatranscriptomic RTs showed distinct abundance patterns across samples compared with metagenomic RTs. The relative abundances of viral and bacterial RTs among identified RT sequences were higher in metatranscriptomes than in metagenomes and these sequences were detected in all metatranscriptome size fractions. Overall, these observations suggest an active proliferation of various RT-assisted elements, which could be involved in genome evolution or adaptive processes of plankton assemblage.

Eukaryotic Components Remodeled Chloroplast Nucleoid Organization during the Green Plant Evolution

Chloroplast (cp) DNA is thought to originate from the ancestral endosymbiont genome and is compacted to form nucleoprotein complexes, cp nucleoids. The structure of cp nucleoids is ubiquitously observed in diverse plants from unicellular algae to flowering plants and is believed to be a multifunctional platform for various processes, including cpDNA replication, repair/recombination, transcription, and inheritance. Despite its fundamental functions, the protein composition for cp nucleoids in flowering plants was suggested to be divergent from those of bacteria and algae, but the evolutionary process remains elusive. In this research, we aimed to reveal the evolutionary history of cp nucleoid organization by analyzing the key organisms representing the three evolutionary stages of eukaryotic phototrophs: the chlorophyte alga Chlamydomonas reinhardtii, the charophyte alga Klebsormidium flaccidum, and the most basal land plant Marchantia polymorpha. To clarify the core cp nucleoid proteins in C. reinhardtii, we performed an LC-MS/MS analysis using highly purified cp nucleoid fractions and identified a novel SAP domain-containing protein with a eukaryotic origin as a constitutive core component. Then, homologous genes for cp nucleoid proteins were searched for in C. reinhardtii, K. flaccidum, and M. polymorpha using the genome databases, and their intracellular localizations and DNA binding activities were investigated by cell biological/biochemical analyses. Based on these results, we propose a model that recurrent modification of cp nucleoid organization by eukaryotic factors originally related to chromatin organization might have been the driving force for the diversification of cp nucleoids since the early stage of green plant evolution.

RNA-Seq transcriptomic analysis with Bag2D software identifies key pathways enhancing lipid yield in a high lipid-producing mutant of the non-model green alga Dunaliella tertiolecta

BACKGROUND

For many years, increasing demands for fossil fuels have met with limited supply. As a potential substitute and renewable source of biofuel feedstock, microalgae have received significant attention. However, few of the current algal species produce high lipid yields to be commercially viable. To discover more high yielding strains, next-generation sequencing technology is used to elucidate lipid synthetic pathways and energy metabolism involved in lipid yield. When subjected to manipulation by genetic and metabolic engineering, enhancement of such pathways may further enhance lipid yield.

RESULTS

In this study, transcriptome profiling of a random insertional mutant with enhanced lipid production generated from a non-model marine microalga Dunaliella tertiolecta is presented. D9 mutant has a lipid yield that is 2- to 4-fold higher than that of wild type. Using novel Bag2D-workflow scripts developed and reported here, the non-redundant transcripts from de novo assembly were annotated based on the best hits in five model microalgae, namely Chlamydomonas reinhardtii, Coccomyxa subellipsoidea C-169, Ostreococcus lucimarinus, Volvox carteri, Chlorella variabilis NC64A and a high plant species Arabidopsis thaliana. The assembled contigs (~181 Mb) includes 481,381 contigs, covering 10,185 genes. Pathway analysis showed that a pathway from inositol phosphate metabolism to fatty acid biosynthesis is the most significantly correlated with higher lipid yield in this mutant.

CONCLUSIONS

Herein, we described a pipeline to analyze RNA-Seq data without pre-existing transcriptomic information. The draft transcriptome of D. tertiolecta was constructed and annotated, which offered useful information for characterizing high lipid-producing mutants. D. tertiolecta mutant was generated with an enhanced photosynthetic efficiency and lipid production. RNA-Seq data of the mutant and wild type were compared, providing biological insights into the expression patterns of contigs associated with energy metabolism and carbon flow pathways. Comparison of D. tertiolecta genes with homologs of five other green algae and a model high plant species can facilitate the annotation of D. tertiolecta and lead to a more complete annotation of its sequence database, thus laying the groundwork for optimization of lipid production pathways based on genetic manipulation.

Conversion of a light-driven proton pump into a light-gated ion channel

Interest in microbial rhodopsins with ion pumping activity has been revitalized in the context of optogenetics, where light-driven ion pumps are used for cell hyperpolarization and voltage sensing. We identified an opsin-encoding gene (CsR) in the genome of the arctic alga Coccomyxa subellipsoidea C-169 that can produce large photocurrents in Xenopus oocytes. We used this property to analyze the function of individual residues in proton pumping. Modification of the highly conserved proton shuttling residue R83 or its interaction partner Y57 strongly reduced pumping power. Moreover, this mutation converted CsR at moderate electrochemical load into an operational proton channel with inward or outward rectification depending on the amino acid substitution. Together with molecular dynamics simulations, these data demonstrate that CsR-R83 and its interacting partner Y57 in conjunction with water molecules forms a proton shuttle that blocks passive proton flux during the dark-state but promotes proton movement uphill upon illumination.

Transitions between the Arabidopsis-type and the human-type telomere sequence in green algae (clade Caudivolvoxa, Chlamydomonadales)

Telomeres are nucleoprotein structures that distinguish native chromosomal ends from double-stranded breaks. They are maintained by telomerase that adds short G-rich telomeric repeats at chromosomal ends in most eukaryotes and determines the TnAmGo sequence of canonical telomeres. We employed an experimental approach that was based on detection of repeats added by telomerase to identify the telomere sequence type forming the very ends of chromosomes. Our previous studies that focused on the algal order Chlamydomonadales revealed several changes in telomere motifs that were consistent with the phylogeny and supported the concept of the Arabidopsis-type sequence being the ancestral telomeric motif for green algae. In addition to previously described independent transitions to the Chlamydomonas-type sequence, we report that the ancestral telomeric motif was replaced by the human-type sequence in the majority of algal species grouped within a higher order clade, Caudivolvoxa. The Arabidopsis-type sequence was apparently retained in the Polytominia clade. Regarding the telomere sequence, the Chlorogonia clade within Caudivolvoxa bifurcates into two groups, one with the human-type sequence and the other group with the Arabidopsis-type sequence that is solely formed by the Chlorogonium species. This suggests that reversion to the Arabidopsis-type telomeric motif occurred in the common ancestral Chlorogonium species. The human-type sequence is also synthesized by telomerases of algal strains from Arenicolinia, Dunaliellinia and Stephanosphaerinia, except a distinct subclade within Stephanosphaerinia, where telomerase activity was not detected and a change to an unidentified telomeric motif might arise. We discuss plausible reasons why changes in telomeric motifs were tolerated during evolution of green algae.