Draft genomes and phenotypic characterization of Tisochrysis lutea strains. Toward the production of domesticated strains with high added value

Chromosome-scale genome assembly reveals insights into the evolution and ecology of the harmful algal bloom species Phaeocystis globosa Scherffel

The phytoplankton Phaeocystis globosa plays an important role in sulfur cycling and climate control, and can develop harmful algal blooms (HABs). Here we report a chromosome-scale reference genome assembly of P. globosa, which enable in-depth analysis of molecular underpinnings of important ecological characteristics. Comparative genomic analyses detected two-rounds of genome duplications that may have fueled evolutionary innovations. The genome duplication may have resulted in the formation of dual HiDP and LoDP dimethylsulphoniopropionate (DMSP) biosynthesis pathways in P. globosa. Selective gene family expansions may have strengthened biological pathways critical for colonial formation that is often associated with the development of algal blooms. The copy numbers of rhodopsin genes are variable in different strains, suggesting that rhodopsin genes may play a role in strain-specific adaptation to ecological factors. The successful reconstruction of the P. globosa genome sets up an excellent platform that facilitates in-depth research on bloom development and DMSP metabolism.

Long-read genome sequencing provides novel insights into the harmful algal bloom species Prymnesium parvum

Prymnesium parvum is a toxin-producing haptophyte that causes harmful algal blooms worldwide, which are often associated with massive fish-kills and subsequent economic losses. In here, we present nuclear and plastid genome assemblies using PacBio HiFi long reads and DNBseq short reads for the two P. parvum strains UTEX 2797 and CCMP 3037, representing producers of type A prymnesins. Our results show that the P. parvum strains have a moderate haptophyte genome size of 97.56 and 107.32 Mb. The genome assemblies present one of highest contiguous assembled contig sequences to date consisting of 463 and 362 contigs with a contig N50 of 596.99 kb and 968.39 kb for strain UTEX 2797 and CCMP 3037, respectively. The assembled contigs of UTEX 2797 and CCMP 3037 were anchored to 34 scaffolds, with a scaffold N50 of 5.35 Mb and 3.61 Mb, respectively, accounting for 93.2 % and 97.9 % of the total length. Each plastid genome comprises a circular contig. A total of 20,578 and 19,426 protein-coding genes were annotated for UTEX 2797 and CCMP 3037. The expanded gene family analysis showed that starch and sucrose metabolism, sulfur metabolism, energy metabolism and ABC transporters are involved in the evolution of P. parvum. Polyketide synthase (PKS) genes responsible for the production of secondary metabolites such as prymnesins displayed different expression patterns under nutrient limitation. Overlap with repeats and horizontal gene transfer may be two contributing factors to the high number of PKS genes found in this species. The two high quality P. parvum genomes will serve as valuable resources for ecological, genetic, and toxicological studies of haptophytes that can be used to monitor and potentially manage harmful blooms of ichthyotoxic P. parvum in the future.

Phylogeny and genetic variations of the three genome compartments in haptophytes shed light on the rapid evolution of coccolithophores

Haptophyte algae, including coccolithophores, play key roles in global carbon cycling and ecosystem. They exhibit exceptional morphological and functional diversity. However, their phylogeny is mostly based on short markers and genome researches are always limited to few species, hindering a better understanding about their evolution and diversification. In this study, by assembling 69 new plastid genomes, 65 new mitochondrial genomes, and 55 nuclear drafts, we systematically analyzed their genome variations and built the most comprehensive phylogenies in haptophytes and Noelaerhabdaceae, with the latter is the family of the model coccolithophore Emiliania huxleyi. The haptophyte genomes vary significantly in size, gene content, and structure. We detected phylogenetic incongruence of Prymnesiales between genome compartments. In Noelaerhabdaceae, by including Reticulofenestra sessilis and a proper outgroup, we found R. sessilis was not the basal taxon of this family. Noelaerhabdaceae strains have very similar genomic features and conserved sequences, but different gene content and dynamic structure. We speculate that was caused by DNA double-strand break repairs. Our results provide valuable genetic resources and new insights into the evolution of haptophytes, especially coccolithophores.

Evolution of a chimeric mitochondrial carbonic anhydrase through gene fusion in a haptophyte alga

Carbonic anhydrases (CAs) are a universal enzyme family that catalyses the interconversion of carbon dioxide and bicarbonate, and they are localized in most compartments including mitochondria and plastids. Thus far, eight classes of CAs (α-, β-, γ-, δ-, ζ-, η-, θ- and ι-CA) have been characterized. This study reports an interesting gene encoding a fusion protein of β-CA and ι-CA found in the haptophyte Isochrysis galbana. Recombinant protein assays demonstrated that the C-terminal ι-CA region catalyses CO₂ hydration, whereas the N-terminal β-CA region no longer exhibits enzymatic activity. Considering that haptophytes generally have mitochondrion-localized β-CAs and plastid-localized ι-CAs, the fusion CA would show an intermediate stage in which mitochondrial β-CA is replaced by ι-CA in a haptophyte species.

The Genome of the Haptophyte Diacronema lutheri (Pavlova lutheri, Pavlovales): A Model for Lipid Biosynthesis in Eukaryotic Algae

Haptophytes are biogeochemically and industrially important protists with underexplored genomic diversity. We present a nuclear genome assembly for the class Pavlovales, which was assembled with PacBio long-read data into highly contiguous sequences. We sequenced strain Diacronema lutheri NIVA-4/92, formerly known as Pavlova lutheri, because it has established roles in aquaculture and has been a key organism for studying microalgal lipid biosynthesis. Our data show that D. lutheri has the smallest and most streamlined haptophycean genome assembled to date, with an assembly size of 43.503 Mb and 14,446 protein-coding genes. Together with its high nuclear GC content, Diacronema is an important genus for investigating selective pressures on haptophyte genome evolution, contrasting with the much larger and more repetitive genome of the coccolithophore Emiliania huxleyi. The D. lutheri genome will be a valuable resource for resolving the genetic basis of algal lipid biosynthesis and metabolic remodeling that takes place during adaptation and stress response in natural and engineered environments.

Isolation of Four Microalgal Strains From the Lake Massaciuccoli: Screening of Common Pollutants Tolerance Pattern and Perspectives for Their Use in Biotechnological Applications

Aquatic ecosystems represent one of the largest reservoirs of phytoplankton accounting for most of the primary production of the Earth. The Lake Massaciuccoli located in Tuscany (Italy) is one of the largest swamps that in ancient times entirely covered the Versilia coastal plain. Despite its peculiar features, especially the eutrophic characteristics, its native microalgal consortia have never been explored up to now. In this work, we isolated and described four autochthonous microalgal strains from different sites in the lake (FB, Idr, CL_Sc, and CL_Ch); the four microalgal strains were identified within the Chlorella sorokiniana clade. We exposed them to ten of the most common or emerging environmental contaminants in order to describe their preliminary response to the tested substances: five metals (As, Fe, Ni, Cu, and Zn), two herbicides (Metolachlor and Sethoxydim), two antibiotics (Ciprofloxacin and Benzylpenicillin) and a non-steroidal anti-inflammatory drug (Ibuprofen). Physiological response of the strains highlighted intraspecific differences; strain CL_Sc was the most tolerant in presence of metals while strain Idr was the most sensitive. All strains were sensitive to sethoxydim and tolerant to metolachlor at all the tested concentrations. Strains FB and Idr were the most sensitive in presence of Ibuprofen while strain CL_Ch was the most sensitive to the highest Benzylpenicillin concentration. Resistance pattern of strain Idr somehow reflects both the phylogenetic and the geographic "isolation" from all other three strains. Finally, optical microscope observation confirmed some differences also in the microalgae morphological aspect. Overall, all the strains showed interesting responses in presence of high concentrations of the tested substances, representing putative interesting candidates for water remediation in wastewater treatment plants.

Enhancing PUFA-rich polar lipids in Tisochrysis lutea using adaptive laboratory evolution (ALE) with oscillating thermal stress

Adaptive laboratory evolution is a powerful tool for microorganism improvement likely to produce enhanced microalgae better tailored to their industrial uses. In this work, 12 wild-type strains of Tisochrysis lutea were co-cultivated under increasing thermal stress for 6 months. Indeed, temperature was oscillating daily between a high and a low temperature, with increasing amplitude along the experiment. The goal was to enhance the polyunsaturated fatty acid content of the polar lipids. Samples were taken throughout the evolution experiment and cultivated in standardized conditions to analyze the evolution of the lipid profile. Genomic analysis of the final population shows that two strains survived. The lipid content doubled, impacting all lipid classes. The fatty acid analyses show a decrease in SFAs correlated with an increase in monounsaturated fatty acids (MUFAs), while changes in polyunsaturated fatty acid (PUFAs) vary between both photobioreactors. Hence, the proportion of C18-MUFAs (18:1 n-9) and most C18-PUFAs (18:2 n-6, 18:3 n-3, and 18:4 n-3) increased, suggesting their potential role in adjusting membrane fluidity to temperature shifts. Of particular interest, DHA in polar lipids tripled in the final population while the growth rate was not affected. KEY POINTS: • Adaptive laboratory evolution on a mix of 12 T. lutea strains led to survival of 2 • Thermal stress impacted cell size, total lipid cell content, and all lipid classes • DHA cell content partitioned to polar lipids tripled throughout the experiment.

Repeated species radiations in the recent evolution of the key marine phytoplankton lineage Gephyrocapsa

Phytoplankton account for nearly half of global primary productivity and strongly affect the global carbon cycle, yet little is known about the forces that drive the evolution of these keystone microscopic organisms. Here we combine morphometric data from the fossil record of the ubiquitous coccolithophore genus Gephyrocapsa with genomic analyses of extant species to assess the genetic processes underlying Pleistocene palaeontological patterns. We demonstrate that all modern diversity in Gephyrocapsa (including Emiliania huxleyi) originated in a rapid species radiation during the last 0.6 Ma, coincident with the latest of the three pulses of Gephyrocapsa diversification and extinction documented in the fossil record. Our evolutionary genetic analyses indicate that new species in this genus have formed in sympatry or parapatry, with occasional hybridisation between species. This sheds light on the mode of speciation during evolutionary radiation of marine phytoplankton and provides a model of how new plankton species form.

The phytoplankton Gephyrocapsa have gone through repeated macroevolutionary shifts in size. Here, Bendif et al. combine fossil and genomic data to show the latest shift was coincident with a species radiation and suggest that previous shifts have also resulted from cycles of radiation and extinction.

How haptophytes microalgae mitigate vitamin B12 limitation

Vitamin B₁₂ (cobalamin) can control phytoplankton development and community composition, with around half of microalgal species requiring this vitamin for growth. B₁₂ dependency is determined by the absence of cobalamin-independent methionine synthase and is unrelated across lineages. Despite their important role in carbon and sulphur biogeochemistry, little is known about haptophytes utilization of vitamin B₁₂ and their ability to cope with its limitation. Here we report the first evaluation of B₁₂ auxotrophy among this lineage based on molecular data of 19 species from 9 families. We assume that all species encode only a B₁₂-dependent methionine synthase, suggesting ubiquitous B₁₂ auxotrophy in this phylum. We further address the effect of different B₁₂ limitations on the molecular physiology of the model haptophyte Tisochrysis lutea. By coupling growth assays in batch and chemostat to cobalamin quantification and expression analyses, we propose that haptophytes use three strategies to cope with B₁₂ limitation. Haptophytes may assimilate dissolved methionine, finely regulate genes involved in methionine cycle and B₁₂ transport and/or limit B₁₂ transport to the mitochondrion. Taken together, these results provide better understanding of B₁₂ metabolism in haptophytes and represent valuable data for deciphering how B₁₂-producing bacteria shape the structure and dynamics of this important phytoplankton community.

Comparative and Functional Algal Genomics

Over 100 whole-genome sequences from algae are published or soon to be published. The rapidly increasing availability of these fundamental resources is changing how we understand one of the most diverse, complex, and understudied groups of photosynthetic eukaryotes. Genome sequences provide a window into the functional potential of individual algae, with phylogenomics and functional genomics as tools for contextualizing and transferring knowledge from reference organisms into less well-characterized systems. Remarkably, over half of the proteins encoded by algal genomes are of unknown function, highlighting the volume of functional capabilities yet to be discovered. In this review, we provide an overview of publicly available algal genomes, their associated protein inventories, and their quality, with a summary of the statuses of protein function understanding and predictions.

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.