Transcriptome analysis of functional differentiation between haploid and diploid cells of Emiliania huxleyi, a globally significant photosynthetic calcifying cell

Phagocytosis in Marine Coccolithophore Gephyrocapsa huxleyi: Comparison between Calcified and Non-Calcified Strains

Coccolithophores play a significant role in marine calcium carbonate production and carbon cycles, attributing to their unique feature of producing calcareous plates, coccoliths. Coccolithophores also possess a haplo-diplontic life cycle, presenting distinct morphology types and calcification states. However, differences in nutrient acquisition strategies and mixotrophic behaviors of the two life phases remain unclear. In this study, we conducted a series of phagocytosis experiments of calcified diploid and non-calcified haploid strains of coccolithophore Gephyrocapsa huxleyi under light and dark conditions. The phagocytosis capability of each strain was examined based on characteristic fluorescent signals from ingested beads using flow cytometry and fluorescence microscopy. The results show a significantly higher phagocytosis percentage on fluorescent beads in the bacterial prey surrogates of the non-calcified haploid Gephyrocapsa huxleyi strain, than the calcified diploid strain with or without light. In addition, the non-calcified diploid cells seemingly to presented a much higher phagocytosis percentage in darkness than under light. The differential phagocytosis capacities between the calcified diploid and non-calcified haploid Gephyrocapsa huxleyi strains indicate potential distinct nutritional strategies at different coccolithophore life and calcifying stages, which may further shed light on the potential strategies that coccolithophore possesses in unfavorable environments such as twilight zones and the expanding coccolithophore niches in the natural marine environment under the climate change scenario.

Gephyrocapsa huxleyi (Emiliania huxleyi) as a model system for coccolithophore biology

Coccolithophores are the most abundant calcifying organisms in modern oceans and are important primary producers in many marine ecosystems. Their ability to generate a cellular covering of calcium carbonate plates (coccoliths) plays a major role in marine biogeochemistry and the global carbon cycle. Coccolithophores also play an important role in sulfur cycling through the production of the climate-active gas dimethyl sulfide. The primary model organism for coccolithophore research is Emiliania huxleyi, now named Gephyrocapsa huxleyi. G. huxleyi has a cosmopolitan distribution, occupying coastal and oceanic environments across the globe, and is the most abundant coccolithophore in modern oceans. Research in G. huxleyi has identified many aspects of coccolithophore biology, from cell biology to ecological interactions. In this perspective, we summarize the key advances made using G. huxleyi and examine the emerging tools for research in this model organism. We discuss the key steps that need to be taken by the research community to advance G. huxleyi as a model organism and the suitability of other species as models for specific aspects of coccolithophore biology.

Maintaining beneficial alga-associated bacterial communities under heat stress: insights from controlled co-culture experiments using antibiotic-resistant bacterial strains

Brown algae, like many eukaryotes, possess diverse microbial communities. Ectocarpus-a model brown alga-relies on these communities for essential processes, such as growth development. Controlled laboratory systems are needed for functional studies of these algal-bacterial interactions. We selected bacterial strains based on their metabolic networks to provide optimal completion of the algal metabolism, rendered them resistant to two antibiotics, and inoculate them to establish controlled co-cultures with Ectocarpus under continuous antibiotic treatment. We then monitored the stability of the resulting associations under control conditions and heat stress using 16S metabarcoding. Antibiotics strongly reduced bacterial diversity both in terms of taxonomy and predicted metabolic functions. In the inoculated sample, 63%-69% of reads corresponded to the inoculated strains, and the communities remained stable during temperature stress. They also partially restored the predicted metabolic functions of the natural community. Overall, the development of antibiotic-resistant helper cultures offers a promising route to fully controlled laboratory experiments with algae and microbiota and thus represents an important step towards generating experimental evidence for specific host-microbe interactions in the systems studied. Further work will be required to achieve full control and progressively expand our repertoire of helper strains including those currently 'unculturable'.

Genome-Wide Identification and Analysis of MYB Transcription Factors in Pyropia yezoensis

MYB transcription factors are one of the largest transcription factor families in plants, and they regulate numerous biological processes. Red algae are an important taxonomic group and have important roles in economics and research. However, no comprehensive analysis of the MYB gene family in any red algae, including Pyropia yezoensis, has been conducted. To identify the MYB gene members of Py. yezoensis, and to investigate their family structural features and expression profile characteristics, a study was conducted. In this study, 3 R2R3-MYBs and 13 MYB-related members were identified in Py. yezoensis. Phylogenetic analysis indicated that most red algae MYB genes could be clustered with green plants or Glaucophyta MYB genes, inferring their ancient origins. Synteny analysis indicated that 13 and 5 PyMYB genes were orthologous to Pyropia haitanensis and Porphyra umbilicalis, respectively. Most Bangiaceae MYB genes contain several Gly-rich motifs, which may be the result of an adaptation to carbon limitations and maintenance of important regulatory functions. An expression profile analysis showed that PyMYB genes exhibited diverse expression profiles. However, the expression patterns of different members appeared to be diverse, and PyMYB5 was upregulated in response to dehydration, low temperature, and Pythium porphyrae infection. This is the first comprehensive study of the MYB gene family in Py. Yezoensis and it provides vital insights into the functional divergence of MYB genes.

A joint proteomic and genomic investigation provides insights into the mechanism of calcification in coccolithophores

Coccolithophores are globally abundant, calcifying microalgae that have profound effects on marine biogeochemical cycles, the climate, and life in the oceans. They are characterized by a cell wall of CaCO3 scales called coccoliths, which may contribute to their ecological success. The intricate morphologies of coccoliths are of interest for biomimetic materials synthesis. Despite the global impact of coccolithophore calcification, we know little about the molecular machinery underpinning coccolithophore biology. Working on the model Emiliania huxleyi, a globally distributed bloom-former, we deploy a range of proteomic strategies to identify coccolithogenesis-related proteins. These analyses are supported by a new genome, with gene models derived from long-read transcriptome sequencing, which revealed many novel proteins specific to the calcifying haptophytes. Our experiments provide insights into proteins involved in various aspects of coccolithogenesis. Our improved genome, complemented with transcriptomic and proteomic data, constitutes a new resource for investigating fundamental aspects of coccolithophore biology.

High-Depth Transcriptome Reveals Differences in Natural Haploid Ginkgo biloba L. Due to the Effect of Reduced Gene Dosage

As a representative of gymnosperms, the discovery of natural haploids of Ginkgo biloba L. has opened a new door for its research. Haploid germplasm has always been a research material of interest to researchers because of its special characteristics. However, we do not yet know the special features and mechanisms of haploid ginkgo following this significant discovery. In this study, we conducted a homogenous garden experiment on haploid and diploid ginkgo to explore the differences in growth, physiology and biochemistry between the two. Additionally, a high-depth transcriptome database of both was established to reveal their transcriptional differences. The results showed that haploid ginkgo exhibited weaker growth potential, lower photosynthesis and flavonoid accumulation capacity. Although the up-regulated expression of DEGs in haploid ginkgo reached 46.7% of the total DEGs in the whole transcriptome data, the gene sets of photosynthesis metabolic, glycolysis/gluconeogenesis and flavonoid biosynthesis pathways, which were significantly related to these differences, were found to show a significant down-regulated expression trend by gene set enrichment analysis (GSEA). We further found that the major metabolic pathways in the haploid ginkgo transcriptional database were down-regulated in expression compared to the diploid. This study reveals for the first time the phenotypic, growth and physiological differences in haploid ginkgos, and demonstrates their transcriptional patterns based on high-depth transcriptomic data, laying the foundation for subsequent in-depth studies of haploid ginkgos.

Rampant nuclear-mitochondrial-plastid phylogenomic discordance in globally distributed calcifying microalgae

Incongruent phylogenies have been widely observed between nuclear and plastid or mitochondrial genomes in terrestrial plants and animals. However, few studies have examined these patterns in microalgae or the discordance between the two organelles. Here we investigated the nuclear-mitochondrial-plastid phylogenomic incongruence in Emiliania-Gephyrocapsa, a group of cosmopolitan calcifying phytoplankton with enormous populations and recent speciations. We assembled mitochondrial and plastid genomes of 27 strains from across global oceans and temperature regimes, and analyzed the phylogenomic histories of the three compartments using concatenation and coalescence methods. Six major clades with varying morphology and distribution are well recognized in the nuclear phylogeny, but such relationships are absent in the mitochondrial and plastid phylogenies, which also differ substantially from each other. The rampant phylogenomic discordance is due to a combination of organellar capture (introgression), organellar genome recombination, and incomplete lineage sorting of ancient polymorphic organellar genomes. Hybridization can lead to replacements of whole organellar genomes without introgression of nuclear genes and the two organelles are not inherited as a single cytoplasmic unit. This study illustrates the convoluted evolution and inheritance of organellar genomes in isogamous haplodiplontic microalgae and provides a window into the phylogenomic complexity of marine unicellular eukaryotes.

Evolutionary Rates in the Haptophyta: Exploring Molecular and Phenotypic Diversity

Haptophytes are photosynthetic protists found in both freshwater and marine environments with an origin possibly dating back to the Neoproterozoic era. The most recent molecular phylogeny reveals several haptophyte “mystery clades” that await morphological verification, but it is otherwise highly consistent with morphology-based phylogenies, including that of the coccolithophores (calcifying haptophytes). The fossil coccolith record offers unique insights into extinct lineages, including the adaptive radiations that produced extant descendant species. By combining molecular data of extant coccolithophores and phenotype-based studies of their ancestral lineages, it has become possible to probe the modes and rates of speciation in more detail, although this approach is still limited to only few taxa because of the lack of whole-genome datasets. The evolution of calcification likely involved several steps, but its origin can be traced back to an early association with organic scales typical for all haptophytes. Other key haptophyte traits, including the haplo-diplontic life cycle, are herein mapped upon the coccolithophorid phylogeny to help navigate a discussion of their ecological benefits and trade-offs in a rapidly changing ocean.

Trade-off between sex and growth in diatoms: Molecular mechanisms and demographic implications

Diatoms are fast-growing and winning competitors in aquatic environments, possibly due to optimized growth performance. However, their life cycles are complex, heteromorphic, and not fully understood. Here, we report on the fine control of cell growth and physiology during the sexual phase of the marine diatom Pseudo-nitzschia multistriata. We found that mating, under nutrient replete conditions, induces a prolonged growth arrest in parental cells. Transcriptomic analyses revealed down-regulation of genes related to major metabolic functions from the early phases of mating. Single-cell photophysiology also pinpointed an inhibition of photosynthesis and storage lipids accumulated in the arrested population, especially in gametes and zygotes. Numerical simulations revealed that growth arrest affects the balance between parental cells and their siblings, possibly favoring the new generation. Thus, in addition to resources availability, life cycle traits contribute to shaping the species ecological niches and must be considered to describe and understand the structure of plankton communities.

Resolving the microalgal gene landscape at the strain level: A novel hybrid transcriptome of Emiliania huxleyi CCMP3266

Microalgae are key ecological players with a complex evolutionary history. Genomic diversity, in addition to limited availability of high-quality genomes, challenge studies that aim to elucidate molecular mechanisms underlying microalgal ecophysiology. Here, we present a novel and comprehensive transcriptomic hybrid approach to generate a reference for genetic analyses, and resolve the microalgal gene landscape at the strain level. The approach is demonstrated for a strain of the coccolithophore microalga Emiliania huxleyi, which is a species complex with considerable genome variability. The investigated strain is commonly studied as a model for algal-bacterial interactions, and was therefore sequenced in the presence of bacteria to elicit the expression of interaction-relevant genes. We applied complementary PacBio Iso-Seq full-length cDNA, and poly(A)-independent Illumina total RNA sequencing, which resulted in a de novo assembled, near complete hybrid transcriptome. In particular, hybrid sequencing improved the reconstruction of long transcripts and increased the recovery of full-length transcript isoforms. To use the resulting hybrid transcriptome as a reference for genetic analyses, we demonstrate a method that collapses the transcriptome into a genome-like dataset, termed "synthetic genome" (sGenome). We used the sGenome as a reference to visually confirm the robustness of the CCMP3266 gene assembly, to conduct differential gene expression analysis, and to characterize novel E. huxleyi genes. The newly-identified genes contribute to our understanding of E. huxleyi genome diversification, and are predicted to play a role in microbial interactions. Our transcriptomic toolkit can be implemented in various microalgae to facilitate mechanistic studies on microalgal diversity and ecology. Importance Microalgae are key players in the ecology and biogeochemistry of our oceans. Efforts to implement genomic and transcriptomic tools in laboratory studies involving microalgae suffer from the lack of published genomes. In the case of coccolithophore microalgae, the problem has long been recognized; the model species Emiliania huxleyi is a species complex with genomes composed of a core, and a large variable portion. To study the role of the variable portion in niche adaptation, and specifically in microbial interactions, strain-specific genetic information is required. Here we present a novel transcriptomic hybrid approach, and generated strain-specific genome-like information. We demonstrate our approach on an E. huxleyi strain that is co-cultivated with bacteria. By constructing a "synthetic genome", we generated comprehensive gene annotations that enabled accurate analyses of gene expression patterns. Importantly, we unveiled novel genes in the variable portion of E. huxleyi that play putative roles in microbial interactions.

Dimethyl sulfide mediates microbial predator-prey interactions between zooplankton and algae in the ocean

Phytoplankton are key components of the oceanic carbon and sulfur cycles¹. During bloom events, some species can emit large amounts of the organosulfur volatile dimethyl sulfide (DMS) into the ocean and consequently the atmosphere, where it can modulate aerosol formation and affect climate^2,3. In aquatic environments, DMS plays an important role as a chemical signal mediating diverse trophic interactions. Yet, its role in microbial predator-prey interactions remains elusive with contradicting evidence for its role in either algal chemical defence or in the chemo-attraction of grazers to prey cells^4,5. Here we investigated the signalling role of DMS during zooplankton-algae interactions by genetic and biochemical manipulation of the algal DMS-generating enzyme dimethylsulfoniopropionate lyase (DL) in the bloom-forming alga Emiliania huxleyi⁶. We inhibited DL activity in E. huxleyi cells in vivo using the selective DL-inhibitor 2-bromo-3-(dimethylsulfonio)-propionate⁷ and overexpressed the DL-encoding gene in the model diatom Thalassiosira pseudonana. We showed that algal DL activity did not serve as an anti-grazing chemical defence but paradoxically enhanced predation by the grazer Oxyrrhis marina and other microzooplankton and mesozooplankton, including ciliates and copepods. Consumption of algal prey with induced DL activity also promoted O. marina growth. Overall, our results demonstrate that DMS-mediated grazing may be ecologically important and prevalent during prey-predator dynamics in aquatic ecosystems. The role of algal DMS revealed here, acting as an eat-me signal for grazers, raises fundamental questions regarding the retention of its biosynthetic enzyme through the evolution of dominant bloom-forming phytoplankton in the ocean.

An Emiliania huxleyi pan-transcriptome reveals basal strain specificity in gene expression patterns

Emiliania huxleyi is a cosmopolitan coccolithophore widespread in temperate oceans. This unicellular photoautotroph forms massive recurring blooms that play an important role in large biogeochemical cycles of carbon and sulfur, which play a role in climate change. The mechanism of bloom formation and demise, controlled by giant viruses that routinely infect these blooms, is poorly understood. We generated a pan-transcriptome of E. huxleyi, derived from three strains with different susceptibility to viral infection. Expression profiling of E. huxleyi sensitive and resistant strains showed major basal differences, including many genes that are induced upon viral infection. This suggests that basal gene expression can affect the host metabolic state and the susceptibility of E. huxleyi to viruses. Due to its ecological importance, the pan-transcriptome and its protein translation, applicable to many E. huxleyi strains, is a powerful resource for investigation of eukaryotic microbial communities.

Virus-induced spore formation as a defense mechanism in marine diatoms

Algal viruses are important contributors to carbon cycling, recycling nutrients and organic material through host lysis. Although viral infection has been described as a primary mechanism of phytoplankton mortality, little is known about host defense responses. We show that viral infection of the bloom-forming, planktonic diatom Chaetoceros socialis induces the mass formation of resting spores, a heavily silicified life cycle stage associated with carbon export due to rapid sinking. Although viral RNA was detected within spores, mature virions were not observed. 'Infected' spores were capable of germinating, but did not propagate or transmit infectious viruses. These results demonstrate that diatom spore formation is an effective defense strategy against viral-mediated mortality. They provide a possible mechanistic link between viral infection, bloom termination, and mass carbon export events and highlight an unappreciated role of viruses in regulating diatom life cycle transitions and ecological success.

Association of Phosphatidylinositol-Specific Phospholipase C with Calcium-Induced Biomineralization in the Coccolithophore Emiliania huxleyi

Biomineralization by calcifying microalgae is a precisely controlled intracellular calcification process that produces delicate calcite scales (or coccoliths) in the coccolithophore Emiliania huxleyi (Haptophycea). Despite its importance in biogeochemical cycles and the marine environment globally, the underlying molecular mechanism of intracellular coccolith formation, which requires calcium, bicarbonate, and coccolith-polysaccharides, remains unclear. In E. huxleyi CCMP 371, we demonstrated that reducing the calcium concentration from 10 (ambient seawater) to 0.1 mM strongly restricted coccolith production, which was then recovered by adding 10 mM calcium, irrespective of inorganic phosphate conditions, indicating that coccolith production could be finely controlled by the calcium supply. Using this strain, we investigated the expression of differentially expressed genes (DEGs) to observe the cellular events induced by changes in calcium concentrations. Intriguingly, DEG analysis revealed that the phosphatidylinositol-specific phospholipase C (PI-PLC) gene was upregulated and coccolith production by cells was blocked by the PI-PLC inhibitor U73122 under conditions closely associated with calcium-induced calcification. These findings imply that PI-PLC plays an important role in the biomineralization process of the coccolithophore E. huxleyi.

De novo transcriptome profile of coccolithophorid alga Emiliania huxleyi CCMP371 at different calcium concentrations with proteome analysis

The haptophyte alga Emiliania huxleyi is the most abundant coccolithophore in the modern ocean and produces elaborate calcite crystals, called coccolith, in a separate intracellular compartment known as the coccolith vesicle. Despite the importance of biomineralization in coccolithophores, the molecular mechanism underlying it remains unclear. Understanding this precise machinery at the molecular level will provide the knowledge needed to enable further manipulation of biomineralization. In our previous study, altering the calcium concentration modified the calcifying ability of E. huxleyi CCMP371. Therefore in this study, we tested E. huxleyi cells acclimated to three different calcium concentrations (0, 0.1, and 10 mM). To understand the whole transcript profile at different calcium concentrations, RNA-sequencing was performed and used for de novo assembly and annotation. The differentially expressed genes (DEGs) among the three different calcium concentrations were analyzed. The functional classification by gene ontology (GO) revealed that 'intrinsic component of membrane' was the most enriched of the GO terms at the ambient calcium concentration (10 mM) compared with the limited calcium concentrations (0 and 0.1 mM). Moreover, the DEGs in those comparisons were enriched mainly in 'secondary metabolites biosynthesis, transport and catabolism' and 'signal transduction mechanisms' in the KOG clusters and 'processing in endoplasmic reticulum', and 'ABC transporters' in the KEGG pathways. Furthermore, metabolic pathways involved in protein synthesis were enriched among the differentially expressed proteins. The results of this study provide a molecular profile for understanding the expression of transcripts and proteins in E. huxleyi at different calcium concentrations, which will help to identify the detailed mechanism of its calcification.

Ocean acidification has little effect on the biochemical composition of the coccolithophore Emiliania huxleyi

Owing to the hierarchical organization of biology, from genomes over transcriptomes and proteomes down to metabolomes, there is continuous debate about the extent to which data and interpretations derived from one level, e.g. the transcriptome, are in agreement with other levels, e.g. the metabolome. Here, we tested the effect of ocean acidification (OA; 400 vs. 1000 μatm CO₂) and its modulation by light intensity (50 vs. 300 μmol photons m^-2 s^-1) on the biomass composition (represented by 75 key metabolites) of diploid and haploid life-cycle stages of the coccolithophore Emiliania huxleyi (RCC1216 and RCC1217) and compared these data with interpretations from previous physiological and gene expression screenings. The metabolite patterns showed minor responses to OA in both life-cycle stages. Whereas previous gene expression analyses suggested that the observed increased biomass buildup derived from lipid and carbohydrate storage, this dataset suggests that OA slightly increases overall biomass of cells, but does not significantly alter their metabolite composition. Generally, light was shown to be a more dominant driver of metabolite composition than OA, increasing the relative abundances of amino acids, mannitol and storage lipids, and shifting pigment contents to accommodate increased irradiance levels. The diploid stage was shown to contain vastly more osmolytes and mannitol than the haploid stage, which in turn had a higher relative content of amino acids, especially aromatic ones. Besides the differences between the investigated cell types and the general effects on biomass buildup, our analyses indicate that OA imposes only negligible effects on E. huxleyi´s biomass composition.

Exploring Molecular Signs of Sex in the Marine Diatom Skeletonema marinoi

Sexual reproduction plays a fundamental role in diatom life cycles. It contributes to increasing genetic diversity through meiotic recombination and also represents the phase where large-sized cells are produced to counteract the cell size reduction process that characterizes these microalgae. With the aim to identify genes linked to the sexual phase of the centric planktonic diatom Skeletonema marinoi, we carried out an RNA-seq experiment comparing the expression level of transcripts in sexualized cells with that of large cells not competent for sex. A set of genes involved in meiosis were found upregulated. Despite the fact that flagellate gametes were observed in the sample, we did not detect the expression of genes involved in the synthesis of flagella that were upregulated during sexual reproduction in another centric diatom. A comparison with the set of genes changing during the first phases of sexual reproduction of the pennate diatom Pseudo-nitzschia multistriata revealed the existence of commonalities, including the strong upregulation of genes with an unknown function that we named Sex Induced Genes (SIG). Our results further broadened the panel of genes that can be used as a marker for sexual reproduction of diatoms, crucial for the interpretation of metatranscriptomic datasets.

Estimating the Fitness Effect of Deleterious Mutations During the Two Phases of the Life Cycle: A New Method Applied to the Root-Rot Fungus Heterobasidion parviporum

Many eukaryote species, including taxa such as fungi or algae, have a lifecycle with substantial haploid and diploid phases. A recent theoretical model predicts that such haploid-diploid lifecycles are stable over long evolutionary time scales when segregating deleterious mutations have stronger effects in homozygous diploids than in haploids and when they are partially recessive in heterozygous diploids. The model predicts that effective dominance-a measure that accounts for these two effects-should be close to 0.5 in these species. It also predicts that diploids should have higher fitness than haploids on average. However, an appropriate statistical framework to conjointly investigate these predictions is currently lacking. In this study, we derive a new quantitative genetic model to test these predictions using fitness data of two haploid parents and their diploid offspring, and genome-wide genetic distance between haploid parents. We apply this model to the root-rot basidiomycete fungus Heterobasidion parviporum-a species where the heterokaryotic (equivalent to the diploid) phase is longer than the homokaryotic (haploid) phase. We measured two fitness-related traits (mycelium growth rate and the ability to degrade wood) in both homokaryons and heterokaryons, and we used whole-genome sequencing to estimate nuclear genetic distance between parents. Possibly due to a lack of power, we did not find that deleterious mutations were recessive or more deleterious when expressed during the heterokaryotic phase. Using this model to compare effective dominance among haploid-diploid species where the relative importance of the two phases varies should help better understand the evolution of haploid-diploid life cycles.

Blueprints for the Next Generation of Bioinspired and Biomimetic Mineralised Composites for Bone Regeneration

Coccolithophores are unicellular marine phytoplankton, which produce intricate, tightly regulated, exoskeleton calcite structures. The formation of biogenic calcite occurs either intracellularly, forming ‘wheel-like’ calcite plates, or extracellularly, forming ‘tiled-like’ plates known as coccoliths. Secreted coccoliths then self-assemble into multiple layers to form the coccosphere, creating a protective wall around the organism. The cell wall hosts a variety of unique species-specific inorganic morphologies that cannot be replicated synthetically. Although biomineralisation has been extensively studied, it is still not fully understood. It is becoming more apparent that biologically controlled mineralisation is still an elusive goal. A key question to address is how nature goes from basic building blocks to the ultrafine, highly organised structures found in coccolithophores. A better understanding of coccolithophore biomineralisation will offer new insight into biomimetic and bioinspired synthesis of advanced, functionalised materials for bone tissue regeneration. The purpose of this review is to spark new interest in biomineralisation and gain new insight into coccolithophores from a material science perspective, drawing on existing knowledge from taxonomists, geologists, palaeontologists and phycologists.

Morphological switch to a resistant subpopulation in response to viral infection in the bloom-forming coccolithophore Emiliania huxleyi

Recognizing the life cycle of an organism is key to understanding its biology and ecological impact. Emiliania huxleyi is a cosmopolitan marine microalga, which displays a poorly understood biphasic sexual life cycle comprised of a calcified diploid phase and a morphologically distinct biflagellate haploid phase. Diploid cells (2N) form large-scale blooms in the oceans, which are routinely terminated by specific lytic viruses (EhV). In contrast, haploid cells (1N) are resistant to EhV. Further evidence indicates that 1N cells may be produced during viral infection. A shift in morphology, driven by meiosis, could therefore constitute a mechanism for E. huxleyi cells to escape from EhV during blooms. This process has been metaphorically coined the 'Cheshire Cat' (CC) strategy. We tested this model in two E. huxleyi strains using a detailed assessment of morphological and ploidy-level variations as well as expression of gene markers for meiosis and the flagellate phenotype. We showed that following the CC model, production of resistant cells was triggered during infection. This led to the rise of a new subpopulation of cells in the two strains that morphologically resembled haploid cells and were resistant to EhV. However, ploidy-level analyses indicated that the new resistant cells were diploid or aneuploid. Thus, the CC strategy in E. huxleyi appears to be a life-phase switch mechanism involving morphological remodeling that is decoupled from meiosis. Our results highlight the adaptive significance of morphological plasticity mediating complex host-virus interactions in marine phytoplankton.

Schrödinger's Cheshire Cat: Are Haploid Emiliania huxleyi Cells Resistant to Viral Infection or Not?

Emiliania huxleyi is the main calcite producer on Earth and is routinely infected by a virus (EhV); a double stranded DNA (dsDNA) virus belonging to the family Phycodnaviridae. E. huxleyi exhibits a haplodiploid life cycle; the calcified diploid stage is non-motile and forms extensive blooms. The haploid phase is a non-calcified biflagellated cell bearing organic scales. Haploid cells are thought to resist infection, through a process deemed the “Cheshire Cat” escape strategy; however, a recent study detected the presence of viral lipids in the same haploid strain. Here we report on the application of an E. huxleyi CCMP1516 EhV-86 combined tiling array (TA) that further confirms an EhV infection in the RCC1217 haploid strain, which grew without any signs of cell lysis. Reverse transcription polymerase chain reaction (RT-PCR) and PCR verified the presence of viral RNA in the haploid cells, yet indicated an absence of viral DNA, respectively. These infected cells are an alternative stage of the virus life cycle deemed the haplococcolithovirocell. In this instance, the host is both resistant to and infected by EhV, i.e., the viral transcriptome is present in haploid cells whilst there is no evidence of viral lysis. This superimposed state is reminiscent of Schrödinger’s cat; of being simultaneously both dead and alive.

The origin of carbon isotope vital effects in coccolith calcite

Calcite microfossils are widely used to study climate and oceanography in Earth's geological past. Coccoliths, readily preserved calcite plates produced by a group of single-celled surface-ocean dwelling algae called coccolithophores, have formed a significant fraction of marine sediments since the Late Triassic. However, unlike the shells of foraminifera, their zooplankton counterparts, coccoliths remain underused in palaeo-reconstructions. Precipitated in an intracellular chemical and isotopic microenvironment, coccolith calcite exhibits large and enigmatic departures from the isotopic composition of abiogenic calcite, known as vital effects. Here we show that the calcification to carbon fixation ratio determines whether coccolith calcite is isotopically heavier or lighter than abiogenic calcite, and that the size of the deviation is determined by the degree of carbon utilization. We discuss the theoretical potential for, and current limitations of, coccolith-based CO2 paleobarometry, that may eventually facilitate use of the ubiquitous and geologically extensive sedimentary archive.

Coccolithophore Cell Biology: Chalking Up Progress

Coccolithophores occupy a special position within the marine phytoplankton because of their production of intricate calcite scales, or coccoliths. Coccolithophores are major contributors to global ocean calcification and long-term carbon fluxes. The intracellular production of coccoliths requires modifications to cellular ultrastructure and metabolism that are surveyed here. In addition to calcification, which appears to have evolved with a diverse range of functions, several other remarkable features that likely underpin the ecological and evolutionary success of coccolithophores have recently been uncovered. These include complex and varied life cycle strategies related to abiotic and biotic interactions as well as a range of novel metabolic pathways and nutritional strategies. Together with knowledge of coccolithophore genetic and physiological variability, these findings are beginning to shed new light on species diversity, distribution, and ecological adaptation. Further advances in genetics and functional characterization at the cellular level will likely to lead to a rapid increase in this understanding.

Interactions between Genetic and Ecological Effects on the Evolution of Life Cycles

Sexual reproduction leads to an alternation between haploid and diploid phases, whose relative length varies widely across taxa. Previous genetical models showed that diploid or haploid life cycles may be favored, depending on dominance interactions and on effective recombination rates. By contrast, niche differentiation between haploids and diploids may favor biphasic life cycles, in which development occurs in both phases. In this article, we explore the interplay between genetical and ecological factors, assuming that deleterious mutations affect the competitivity of individuals within their ecological niche and allowing different effects of mutations in haploids and diploids (including antagonistic selection). We show that selection on a modifier gene affecting the relative length of both phases can be decomposed into a direct selection term favoring the phase with the highest mean fitness (due to either ecological differences or differential effects of mutations) and an indirect selection term favoring the phase in which selection is more efficient. When deleterious alleles occur at many loci and in the presence of ecological differentiation between haploids and diploids, evolutionary branching often occurs and leads to the stable coexistence of alleles coding for haploid and diploid cycles, while temporal variations in niche sizes may stabilize biphasic cycles.

Acidification, not carbonation, is the major regulator of carbon fluxes in the coccolithophore Emiliania huxleyi

A combined increase in seawater [CO2 ] and [H(+) ] was recently shown to induce a shift from photosynthetic HCO3 (-) to CO2 uptake in Emiliania huxleyi. This shift occurred within minutes, whereas acclimation to ocean acidification (OA) did not affect the carbon source. To identify the driver of this shift, we exposed low- and high-light acclimated E. huxleyi to a matrix of two levels of dissolved inorganic carbon (1400, 2800 μmol kg(-1) ) and pH (8.15, 7.85) and directly measured cellular O2 , CO2 and HCO3 (-) fluxes under these conditions. Exposure to increased [CO2 ] had little effect on the photosynthetic fluxes, whereas increased [H(+) ] led to a significant decline in HCO3 (-) uptake. Low-light acclimated cells overcompensated for the inhibition of HCO3 (-) uptake by increasing CO2 uptake. High-light acclimated cells, relying on higher proportions of HCO3 (-) uptake, could not increase CO2 uptake and photosynthetic O2 evolution consequently became carbon-limited. These regulations indicate that OA responses in photosynthesis are caused by [H(+) ] rather than by [CO2 ]. The impaired HCO3 (-) uptake also provides a mechanistic explanation for lowered calcification under OA. Moreover, it explains the OA-dependent decrease in photosynthesis observed in high-light grown phytoplankton.

A Bacterial Pathogen Displaying Temperature-Enhanced Virulence of the Microalga Emiliania huxleyi

Emiliania huxleyi is a globally abundant microalga that plays a significant role in biogeochemical cycles. Over the next century, sea surface temperatures are predicted to increase drastically, which will likely have significant effects on the survival and ecology of E. huxleyi. In a warming ocean, this microalga may become increasingly vulnerable to pathogens, particularly those with temperature-dependent virulence. Ruegeria is a genus of Rhodobacteraceae whose population size tracks that of E. huxleyi throughout the alga’s bloom–bust lifecycle. A representative of this genus, Ruegeria sp. R11, is known to cause bleaching disease in a red macroalga at elevated temperatures. To investigate if the pathogenicity of R11 extends to microalgae, it was co-cultured with several cell types of E. huxleyi near the alga’s optimum (18°C), and at an elevated temperature (25°C) known to induce virulence in R11. The algal populations were monitored using flow cytometry and pulse-amplitude modulated fluorometry. Cultures of algae without bacteria remained healthy at 18°C, but lower cell counts in control cultures at 25°C indicated some stress at the elevated temperature. Both the C (coccolith-bearing) and S (scale-bearing swarming) cell types of E. huxleyi experienced a rapid decline resulting in apparent death when co-cultured with R11 at 25°C, but had no effect on N (naked) cell type at either temperature. R11 had no initial negative impact on C and S type E. huxleyi population size or health at 18°C, but caused death in older co-cultures. This differential effect of R11 on its host at 18 and 25°C suggest it is a temperature-enhanced opportunistic pathogen of E. huxleyi. We also detected caspase-like activity in dying C type cells co-cultured with R11, which suggests that programmed cell death plays a role in the death of E. huxleyi triggered by R11 – a mechanism induced by viruses (EhVs) and implicated in E. huxleyi bloom collapse. Given that E. huxleyi has recently been shown to have acquired resistance against EhVs at elevated temperature, bacterial pathogens with temperature-dependent virulence, such as R11, may become much more important in the ecology of E. huxleyi in a warming climate.

Recent Reticulate Evolution in the Ecologically Dominant Lineage of Coccolithophores

The coccolithophore family Noëlaerhabdaceae contains a number of taxa that are very abundant in modern oceans, including the cosmopolitan bloom-forming Emiliania huxleyi. Introgressive hybridization has been suggested to account for incongruences between nuclear, mitochondrial and plastidial phylogenies of morphospecies within this lineage, but the number of species cultured to date remains rather limited. Here, we present the characterization of 5 new Noëlaerhabdaceae culture strains isolated from samples collected in the south-east Pacific Ocean. These were analyzed morphologically using scanning electron microscopy and phylogenetically by sequencing 5 marker genes (nuclear 18S and 28S rDNA, plastidial tufA, and mitochondrial cox1 and cox3 genes). Morphologically, one of these strains corresponded to Gephyrocapsa ericsonii and the four others to Reticulofenestra parvula. Ribosomal gene sequences were near identical between these new strains, but divergent from G. oceanica, G. muellerae, and E. huxleyi. In contrast to the clear distinction in ribosomal phylogenies, sequences from other genomic compartments clustered with those of E. huxleyi strains with which they share an ecological range (i.e., warm temperate to tropical waters). These data provide strong support for the hypothesis of past (and potentially ongoing) introgressive hybridization within this ecologically important lineage and for the transfer of R. parvula to Gephyrocapsa. These results have important implications for understanding the role of hybridization in speciation in vast ocean meta-populations of phytoplankton.

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.

A paneukaryotic genomic analysis of the small GTPase RABL2 underscores the significance of recurrent gene loss in eukaryote evolution

Background

The cilium (flagellum) is a complex cellular structure inherited from the last eukaryotic common ancestor (LECA). A large number of ciliary proteins have been characterized in a few model organisms, but their evolutionary history often remains unexplored. One such protein is the small GTPase RABL2, recently implicated in the assembly of the sperm tail in mammals.

Results

Using the wealth of currently available genome and transcriptome sequences, including data from our on-going sequencing projects, we systematically analyzed the phylogenetic distribution and evolutionary history of RABL2 orthologs. Our dense taxonomic sampling revealed the presence of RABL2 genes in nearly all major eukaryotic lineages, including small “obscure” taxa such as breviates, ancyromonads, malawimonads, jakobids, picozoans, or palpitomonads. The phyletic pattern of RABL2 genes indicates that it was present already in the LECA. However, some organisms lack RABL2 as a result of secondary loss and our present sampling predicts well over 30 such independent events during the eukaryote evolution. The distribution of RABL2 genes correlates with the presence/absence of cilia: not a single well-established cilium-lacking species has retained a RABL2 ortholog. However, several ciliated taxa, most notably nematodes, some arthropods and platyhelminths, diplomonads, and ciliated subgroups of apicomplexans and embryophytes, lack RABL2 as well, suggesting some simplification in their cilium-associated functions. On the other hand, several algae currently unknown to form cilia, e.g., the “prasinophytes” of the genus Prasinoderma or the ochrophytes Pelagococcus subviridis and Pinguiococcus pyrenoidosus, turned out to encode not only RABL2, but also homologs of some hallmark ciliary proteins, suggesting the existence of a cryptic flagellated stage in their life cycles. We additionally obtained insights into the evolution of the RABL2 gene architecture, which seems to have ancestrally consisted of eight exons subsequently modified not only by lineage-specific intron loss and gain, but also by recurrent loss of the terminal exon encoding a poorly conserved C-terminal extension.

Conclusions

Our comparative analysis supports the notion that RABL2 is an ancestral component of the eukaryotic cilium and underscores the still underappreciated magnitude of recurrent gene loss, or reductive evolution in general, in the history of eukaryotic genomes and cells.

Reviewers

This article was reviewed by Berend Snel and James O. McInerney.

Electronic supplementary material

The online version of this article (doi:10.1186/s13062-016-0107-8) contains supplementary material, which is available to authorized users.

MARINE SULFUR CYCLE. Identification of the algal dimethyl sulfide-releasing enzyme: A missing link in the marine sulfur cycle

Algal blooms produce large amounts of dimethyl sulfide (DMS), a volatile with a diverse signaling role in marine food webs that is emitted to the atmosphere, where it can affect cloud formation. The algal enzymes responsible for forming DMS from dimethylsulfoniopropionate (DMSP) remain unidentified despite their critical role in the global sulfur cycle. We identified and characterized Alma1, a DMSP lyase from the bloom-forming algae Emiliania huxleyi. Alma1 is a tetrameric, redox-sensitive enzyme of the aspartate racemase superfamily. Recombinant Alma1 exhibits biochemical features identical to the DMSP lyase in E. huxleyi, and DMS released by various E. huxleyi isolates correlates with their Alma1 levels. Sequence homology searches suggest that Alma1 represents a gene family present in major, globally distributed phytoplankton taxa and in other marine organisms.

Changes in gene expression of Prymnesium parvum induced by nitrogen and phosphorus limitation

Prymnesium parvum is a globally distributed prymnesiophyte alga commonly found in brackish water marine ecosystems and lakes. It possesses a suite of toxins with ichthyotoxic, cytotoxic and hemolytic effects which, along with its mixotrophic nutritional capabilities, allows it to form massive Ecosystem Disruptive Algal Blooms (EDABs). While blooms of high abundance coincide with high levels of nitrogen (N) and phosphorus (P), reports of field and laboratory studies have noted that P. parvum toxicity appears to be augmented at high N:P ratios or P-limiting conditions. Here we present the results of a comparative analysis of P. parvum RNA-Seq transcriptomes under nutrient replete conditions, and N or P deficiency to understand how this organism responds at the transcriptional level to varying nutrient conditions. In nutrient limited conditions we found diverse transcriptional responses for genes involved in nutrient uptake, protein synthesis and degradation, photosynthesis, and toxin production. As anticipated, when either N or P was limiting, transcription levels of genes encoding transporters for the respective nutrient were higher than those under replete condition. Ribosomal and lysosomal protein genes were expressed at higher levels under either nutrient-limited condition compared to the replete condition. Photosynthesis genes and polyketide synthase genes were more highly expressed under P-limitation but not under N-limitation. These results highlight the ability of P. parvum to mount a coordinated and varied cellular and physiological response to nutrient limitation. Results also provide potential marker genes for further evaluating the physiological response and toxin production of P. parvum populations during bloom formation or to changing environmental conditions.

Morphological and Phylogenetic Characterization of New Gephyrocapsa Isolates Suggests Introgressive Hybridization in the Emiliania/Gephyrocapsa Complex (Haptophyta)

The coccolithophore genus Gephyrocapsa contains a cosmopolitan assemblage of pelagic species, including the bloom-forming Gephyrocapsa oceanica, and is closely related to the emblematic coccolithophore Emiliania huxleyi within the Noëlaerhabdaceae. These two species have been extensively studied and are well represented in culture collections, whereas cultures of other species of this family are lacking. We report on three new strains of Gephyrocapsa isolated into culture from samples from the Chilean coastal upwelling zone using a novel flow cytometric single-cell sorting technique. The strains were characterized by morphological analysis using scanning electron microscopy and phylogenetic analysis of 6 genes (nuclear 18S and 28S rDNA, plastidial 16S and tufA, and mitochondrial cox1 and cox3 genes). Morphometric features of the coccoliths indicate that these isolates are distinct from G. oceanica and best correspond to G. muellerae. Surprisingly, both plastidial and mitochondrial gene phylogenies placed these strains within the E. huxleyi clade and well separated from G. oceanica isolates, making Emiliania appear polyphyletic. The only nuclear sequence difference, 1bp in the 28S rDNA region, also grouped E. huxleyi with the new Gephyrocapsa isolates and apart from G. oceanica. Specifically, the G. muellerae morphotype strains clustered with the mitochondrial β clade of E. huxleyi, which, like G. muellerae, has been associated with cold (temperate and sub-polar) waters. Among putative evolutionary scenarios that could explain these results we discuss the possibility that E. huxleyi is not a valid taxonomic unit, or, alternatively the possibility of past hybridization and introgression between each E. huxleyi clade and older Gephyrocapsa clades. In either case, the results support the transfer of Emiliania to Gephyrocapsa. These results have important implications for relating morphological species concepts to ecological and evolutionary units of diversity.

Phenotypic diversity of diploid and haploid Emiliania huxleyi cells and of cells in different growth phases revealed by comparative metabolomics

In phytoplankton a high species diversity of microalgae co-exists at a given time. But diversity is not only reflected by the species composition. Within these species different life phases as well as different metabolic states can cause additional diversity. One important example is the coccolithophore Emiliania huxleyi. Diploid cells play an important role in marine ecosystems since they can form massively abundant algal blooms but in addition the less abundant haploid life phase of E. huxleyi occurs in lower quantities. Both life phases may fulfill different functions in the plankton. We hypothesize that in addition to the functional diversity caused by this life phase transition the growth stage of cells can also influence the metabolic composition and thus the ecological impact of E. huxleyi. Here we introduce a metabolomic survey in dependence of life phases as well as different growth phases to reveal such changes. The comparative metabolomic approach is based on the extraction of intracellular metabolites from intact microalgae, derivatization and analysis by gas chromatography coupled to mass spectrometry (GC-MS). Automated data processing and statistical analysis using canonical analysis of principal coordinates (CAP) revealed unique metabolic profiles for each life phase. Concerning the correlations of metabolites to growth phases, complex patterns were observed. As for example the saccharide mannitol showed its highest concentration in the exponential phase, whereas fatty acids were correlated to stationary and sterols to declining phase. These results are indicative for specific ecological roles of these stages of E. huxleyi and are discussed in the context of previous physiological and ecological studies.

Life-cycle modification in open oceans accounts for genome variability in a cosmopolitan phytoplankton

Emiliania huxleyi is the most abundant calcifying plankton in modern oceans with substantial intraspecific genome variability and a biphasic life cycle involving sexual alternation between calcified 2N and flagellated 1N cells. We show that high genome content variability in Emiliania relates to erosion of 1N-specific genes and loss of the ability to form flagellated cells. Analysis of 185 E. huxleyi strains isolated from world oceans suggests that loss of flagella occurred independently in lineages inhabiting oligotrophic open oceans over short evolutionary timescales. This environmentally linked physiogenomic change suggests life cycling is not advantageous in very large/diluted populations experiencing low biotic pressure and low ecological variability. Gene loss did not appear to reflect pressure for genome streamlining in oligotrophic oceans as previously observed in picoplankton. Life-cycle modifications might be common in plankton and cause major functional variability to be hidden from traditional taxonomic or molecular markers.

Emiliania huxleyi endures N-limitation with an efficient metabolic budgeting and effective ATP synthesis

Background

Global change will affect patterns of nutrient upwelling in marine environments, potentially becoming even stricter regulators of phytoplankton primary productivity. To better understand phytoplankton nutrient utilization on the subcellular basis, we assessed the transcriptomic responses of the life-cycle stages of the biogeochemically important microalgae Emiliania huxleyi to nitrogen-limitation. Cells grown in batch cultures were harvested at ‘early’ and ‘full’ nitrogen-limitation and were compared with non-limited cells. We applied microarray-based transcriptome profilings, covering ~10.000 known E. huxleyi gene models, and screened for expression patterns that indicate the subcellular responses.

Results

The diploid life-cycle stage scavenges nitrogen from external organic sources and -like diatoms- uses the ornithine-urea cycle to rapidly turn over cellular nitrogen. The haploid stage reacts similarly, although nitrogen scavenging is less pronounced and lipid oxidation is more prominent. Generally, polyamines and proline appear to constitute major organic pools that back up cellular nitrogen. Both stages induce a malate:quinone-oxidoreductase that efficiently feeds electrons into the respiratory chain and drives ATP generation with reduced respiratory carbon throughput.

Conclusions

The use of the ornithine-urea cycle to budget the cellular nitrogen in situations of limitation resembles the responses observed earlier in diatoms. This suggests that underlying biochemical mechanisms are conserved among distant clades of marine phototrophic protists. The ornithine-urea cycle and proline oxidation appear to constitute a sensory-regulatory system that monitors and controls cellular nitrogen budgets under limitation. The similarity between the responses of the life-cycle stages, despite the usage of different genes, also indicates a strong functional consistency in the responses to nitrogen-limitation that appears to be owed to biochemical requirements. The malate:quinone-oxidoreductase is a genomic feature that appears to be absent from diatom genomes, and it is likely to strongly contribute to the uniquely high endurance of E. huxleyi under nutrient limitation.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-1051) contains supplementary material, which is available to authorized users.

Transcriptomic Insights into the Life History of Bolidophytes, the Sister Lineage to Diatoms

Diatoms are perhaps the most diverse lineage of eukaryotic algae, with their siliceous cell wall and diplontic life history often considered to have played important roles in their extraordinary diversification. The characteristic diminution of the diatom cell wall over the course of vegetative growth provides a reliable, intrinsic trigger for sexual reproduction, establishing a direct link between the evolution of their cell-wall and life-history features. It is unclear, however, whether the diplontic life cycle of diatoms represents an ancestral or derived trait. This uncertainty is based in part on our lack of understanding of the life cycle of the sister lineage to diatoms, which includes a mix of two free-living and separately classified forms: naked biflagellate unicells in the genus Bolidomonas and silicified forms in the order Parmales. These two forms might represent different life-history stages, although directly establishing such links can be difficult. We sequenced transcriptomes for Bolidomonas and two diatoms and found that ~0.1% of the coding regions in the two diploid diatoms are heterozygous, whereas Bolidomonas is virtually devoid of heterozygous alleles, consistent with expectations for a haploid genome. These results suggest that Bolidomonas is haploid and predict that parmaleans represent the diploid phase of a haplodiplontic life cycle. These data fill an important gap in our understanding of the origin of the diplontic life history of diatoms, which may represent an evolutionarily derived, adaptive feature.

De novo transcriptome analysis of an aerial microalga Trentepohlia jolithus: pathway description and gene discovery for carbon fixation and carotenoid biosynthesis

Background

Algae in the order Trentepohliales have a broad geographic distribution and are generally characterized by the presence of abundant β-carotene. The many monographs published to date have mainly focused on their morphology, taxonomy, phylogeny, distribution and reproduction; molecular studies of this order are still rare. High-throughput RNA sequencing (RNA-Seq) technology provides a powerful and efficient method for transcript analysis and gene discovery in Trentepohlia jolithus.

Methods/Principal Findings

Illumina HiSeq 2000 sequencing generated 55,007,830 Illumina PE raw reads, which were assembled into 41,328 assembled unigenes. Based on NR annotation, 53.28% of the unigenes (22,018) could be assigned to gene ontology classes with 54 subcategories and 161,451 functional terms. A total of 26,217 (63.44%) assembled unigenes were mapped to 128 KEGG pathways. Furthermore, a set of 5,798 SSRs in 5,206 unigenes and 131,478 putative SNPs were identified. Moreover, the fact that all of the C4 photosynthesis genes exist in T. jolithus suggests a complex carbon acquisition and fixation system. Similarities and differences between T. jolithus and other algae in carotenoid biosynthesis are also described in depth.

Conclusions/Significance

This is the first broad transcriptome survey for T. jolithus, increasing the amount of molecular data available for the class Ulvophyceae. As well as providing resources for functional genomics studies, the functional genes and putative pathways identified here will contribute to a better understanding of carbon fixation and fatty acid and carotenoid biosynthesis in T. jolithus.

Genome variations associated with viral susceptibility and calcification in Emiliania huxleyi

Emiliania huxleyi, a key player in the global carbon cycle is one of the best studied coccolithophores with respect to biogeochemical cycles, climatology, and host-virus interactions. Strains of E. huxleyi show phenotypic plasticity regarding growth behaviour, light-response, calcification, acidification, and virus susceptibility. This phenomenon is likely a consequence of genomic differences, or transcriptomic responses, to environmental conditions or threats such as viral infections. We used an E. huxleyi genome microarray based on the sequenced strain CCMP1516 (reference strain) to perform comparative genomic hybridizations (CGH) of 16 E. huxleyi strains of different geographic origin. We investigated the genomic diversity and plasticity and focused on the identification of genes related to virus susceptibility and coccolith production (calcification). Among the tested 31940 gene models a core genome of 14628 genes was identified by hybridization among 16 E. huxleyi strains. 224 probes were characterized as specific for the reference strain CCMP1516. Compared to the sequenced E. huxleyi strain CCMP1516 variation in gene content of up to 30 percent among strains was observed. Comparison of core and non-core transcripts sets in terms of annotated functions reveals a broad, almost equal functional coverage over all KOG-categories of both transcript sets within the whole annotated genome. Within the variable (non-core) genome we identified genes associated with virus susceptibility and calcification. Genes associated with virus susceptibility include a Bax inhibitor-1 protein, three LRR receptor-like protein kinases, and mitogen-activated protein kinase. Our list of transcripts associated with coccolith production will stimulate further research, e.g. by genetic manipulation. In particular, the V-type proton ATPase 16 kDa proteolipid subunit is proposed to be a plausible target gene for further calcification studies.

Emiliania huxleyi increases calcification but not expression of calcification-related genes in long-term exposure to elevated temperature and pCO2

Increased atmospheric pCO2 is expected to render future oceans warmer and more acidic than they are at present. Calcifying organisms such as coccolithophores that fix and export carbon into the deep sea provide feedbacks to increasing atmospheric pCO2. Acclimation experiments suggest negative effects of warming and acidification on coccolithophore calcification, but the ability of these organisms to adapt to future environmental conditions is not well understood. Here, we tested the combined effect of pCO2 and temperature on the coccolithophore Emiliania huxleyi over more than 700 generations. Cells increased inorganic carbon content and calcification rate under warm and acidified conditions compared with ambient conditions, whereas organic carbon content and primary production did not show any change. In contrast to findings from short-term experiments, our results suggest that long-term acclimation or adaptation could change, or even reverse, negative calcification responses in E. huxleyi and its feedback to the global carbon cycle. Genome-wide profiles of gene expression using RNA-seq revealed that genes thought to be essential for calcification are not those that are most strongly differentially expressed under long-term exposure to future ocean conditions. Rather, differentially expressed genes observed here represent new targets to study responses to ocean acidification and warming.

Dissecting the impact of CO2 and pH on the mechanisms of photosynthesis and calcification in the coccolithophore Emiliania huxleyi

Coccolithophores are important calcifying phytoplankton predicted to be impacted by changes in ocean carbonate chemistry caused by the absorption of anthropogenic CO2 . However, it is difficult to disentangle the effects of the simultaneously changing carbonate system parameters (CO2 , bicarbonate, carbonate and protons) on the physiological responses to elevated CO2 . Here, we adopted a multifactorial approach at constant pH or CO2 whilst varying dissolved inorganic carbon (DIC) to determine physiological and transcriptional responses to individual carbonate system parameters. We show that Emiliania huxleyi is sensitive to low CO2 (growth and photosynthesis) and low bicarbonate (calcification) as well as low pH beyond a limited tolerance range, but is much less sensitive to elevated CO2 and bicarbonate. Multiple up-regulated genes at low DIC bear the hallmarks of a carbon-concentrating mechanism (CCM) that is responsive to CO2 and bicarbonate but not to pH. Emiliania huxleyi appears to have evolved mechanisms to respond to limiting rather than elevated CO2 . Calcification does not function as a CCM, but is inhibited at low DIC to allow the redistribution of DIC from calcification to photosynthesis. The presented data provides a significant step in understanding how E. huxleyi will respond to changing carbonate chemistry at a cellular level.

Ocean acidification reduces growth and calcification in a marine dinoflagellate

Ocean acidification is considered a major threat to marine ecosystems and may particularly affect calcifying organisms such as corals, foraminifera and coccolithophores. Here we investigate the impact of elevated pCO₂ and lowered pH on growth and calcification in the common calcareous dinoflagellate Thoracosphaera heimii. We observe a substantial reduction in growth rate, calcification and cyst stability of T. heimii under elevated pCO₂. Furthermore, transcriptomic analyses reveal CO₂ sensitive regulation of many genes, particularly those being associated to inorganic carbon acquisition and calcification. Stable carbon isotope fractionation for organic carbon production increased with increasing pCO₂ whereas it decreased for calcification, which suggests interdependence between both processes. We also found a strong effect of pCO₂ on the stable oxygen isotopic composition of calcite, in line with earlier observations concerning another T. heimii strain. The observed changes in stable oxygen and carbon isotope composition of T. heimii cysts may provide an ideal tool for reconstructing past seawater carbonate chemistry, and ultimately past pCO₂. Although the function of calcification in T. heimii remains unresolved, this trait likely plays an important role in the ecological and evolutionary success of this species. Acting on calcification as well as growth, ocean acidification may therefore impose a great threat for T. heimii.

Responses of the Emiliania huxleyi proteome to ocean acidification

Ocean acidification due to rising atmospheric CO₂ is expected to affect the physiology of important calcifying marine organisms, but the nature and magnitude of change is yet to be established. In coccolithophores, different species and strains display varying calcification responses to ocean acidification, but the underlying biochemical properties remain unknown. We employed an approach combining tandem mass-spectrometry with isobaric tagging (iTRAQ) and multiple database searching to identify proteins that were differentially expressed in cells of the marine coccolithophore species Emiliania huxleyi (strain NZEH) between two CO₂ conditions: 395 (∼current day) and ∼1340 p.p.m.v. CO₂. Cells exposed to the higher CO₂ condition contained more cellular particulate inorganic carbon (CaCO₃) and particulate organic nitrogen and carbon than those maintained in present-day conditions. These results are linked with the observation that cells grew slower under elevated CO₂, indicating cell cycle disruption. Under high CO₂ conditions, coccospheres were larger and cells possessed bigger coccoliths that did not show any signs of malformation compared to those from cells grown under present-day CO₂ levels. No differences in calcification rate, particulate organic carbon production or cellular organic carbon: nitrogen ratios were observed. Results were not related to nutrient limitation or acclimation status of cells. At least 46 homologous protein groups from a variety of functional processes were quantified in these experiments, of which four (histones H2A, H3, H4 and a chloroplastic 30S ribosomal protein S7) showed down-regulation in all replicates exposed to high CO₂, perhaps reflecting the decrease in growth rate. We present evidence of cellular stress responses but proteins associated with many key metabolic processes remained unaltered. Our results therefore suggest that this E. huxleyi strain possesses some acclimation mechanisms to tolerate future CO₂ scenarios, although the observed decline in growth rate may be an overriding factor affecting the success of this ecotype in future oceans.

Ocean acidification affects redox-balance and ion-homeostasis in the life-cycle stages of Emiliania huxleyi

Ocean Acidification (OA) has been shown to affect photosynthesis and calcification in the coccolithophore Emiliania huxleyi, a cosmopolitan calcifier that significantly contributes to the regulation of the biological carbon pumps. Its non-calcifying, haploid life-cycle stage was found to be relatively unaffected by OA with respect to biomass production. Deeper insights into physiological key processes and their dependence on environmental factors are lacking, but are required to understand and possibly estimate the dynamics of carbon cycling in present and future oceans. Therefore, calcifying diploid and non-calcifying haploid cells were acclimated to present and future CO(2) partial pressures (pCO(2); 38.5 Pa vs. 101.3 Pa CO(2)) under low and high light (50 vs. 300 µmol photons m(-2) s(-1)). Comparative microarray-based transcriptome profiling was used to screen for the underlying cellular processes and allowed to follow up interpretations derived from physiological data. In the diplont, the observed increases in biomass production under OA are likely caused by stimulated production of glycoconjugates and lipids. The observed lowered calcification under OA can be attributed to impaired signal-transduction and ion-transport. The haplont utilizes distinct genes and metabolic pathways, reflecting the stage-specific usage of certain portions of the genome. With respect to functionality and energy-dependence, however, the transcriptomic OA-responses resemble those of the diplont. In both life-cycle stages, OA affects the cellular redox-state as a master regulator and thereby causes a metabolic shift from oxidative towards reductive pathways, which involves a reconstellation of carbon flux networks within and across compartments. Whereas signal transduction and ion-homeostasis appear equally OA-sensitive under both light intensities, the effects on carbon metabolism and light physiology are clearly modulated by light availability. These interactive effects can be attributed to the influence of OA and light on the redox equilibria of NAD and NADP, which function as major sensors for energization and stress. This generic mode of action of OA may therefore provoke similar cell-physiological responses in other protists.

Coccolithogenesis In Scyphosphaera apsteinii (Prymnesiophyceae)

Coccolithophores are the most significant producers of marine biogenic calcite, although the intracellular calcification process is poorly understood. In the case of Scyphosphaera apsteinii Lohmann 1902, flat ovoid muroliths and bulky, vase-shaped lopadoliths with a range of intermediate morphologies may be produced by a single cell. This polymorphic species is within the Zygodiscales, a group that remains understudied with respect to ultrastructure and coccolith ontogeny. We therefore undertook an analysis of cell ultrastructure, morphology, and coccolithogenesis. The cell ultrastructure showed many typical haptophyte features, with calcification following a similar pattern to that described for other heterococcolith bearing species including Emiliania huxleyi. Of particular significance was the reticular body role in governing fine-scale morphology, specifically the central pore formation of the coccolith. Our observations also highlighted the essential role of the inter- and intracrystalline organic matrix in growth and arrangement of the coccolith calcite. S. apsteinii secreted mature coccoliths that attached to the plasma membrane via fibrillar material. Time-lapse light microscopy demonstrated secretion of lopadoliths occurred base first before being actively repositioned at the cell surface. Significantly, growth irradiance influenced the coccosphere composition with fewer lopadoliths being formed relative to muroliths at higher light intensities. Overall, our observations support dynamic metabolic (i.e., in response to growth irradiance), sensory and cytoskeletal control over the morphology and secretion of polymorphic heterococcoliths. With a basic understanding of calcification established, S. apsteinii could be a valuable model to further study coccolithophore calcification and cell physiological responses to ocean acidification.