RNA-Seq as an Emerging Tool for Marine Dinoflagellate Transcriptome Analysis: Process and Challenges

Modification of the Trizol Method for the Extraction of RNA from Prorocentrum triestinum ACIZ_LEM2

In samples of harmful algal blooms (HABs), seawater can contain a high abundance of microorganisms and elemental ions. Along with the hardness of the walls of key HAB dinoflagellates such as Prorocentrum triestinum, this makes RNA extraction very difficult. These components interfere with RNA isolation, causing its degradation, in addition to the complex seawater properties of HABs that could hinder RNA isolation for effective RNA sequencing and transcriptome profiling. In this study, an RNA isolation technique was established through the modification of the Trizol method by applying the Micropestle System on cell pellets of P. triestinum frozen at -20 °C, obtained from 400 mL of culture with a total of 107 cells/mL. The results of the modified Trizol protocol generated quality RNA samples for transcriptomics sequencing, as determined by their measurement in Analyzer Agilent 4150.

Insights into Alexandrium minutum Nutrient Acquisition, Metabolism and Saxitoxin Biosynthesis through Comprehensive Transcriptome Survey

Simple Summary

Alexandrium minutum is one of the causing organisms for the occurrence of harmful algae bloom (HABs) in marine ecosystems. This species produces saxitoxin, one of the deadliest neurotoxins which can cause human mortality. However, molecular information such as genes and proteins catalog on this species is still lacking. Therefore, this study has successfully characterized several new molecular mechanisms regarding A. minutum environmental adaptation and saxitoxin biosynthesis. Ultimately, this study provides a valuable resource for facilitating future dinoflagellates’ molecular response to environmental changes.

Abstract

The toxin-producing dinoflagellate Alexandrium minutum is responsible for the outbreaks of harmful algae bloom (HABs). It is a widely distributed species and is responsible for producing paralytic shellfish poisoning toxins. However, the information associated with the environmental adaptation pathway and toxin biosynthesis in this species is still lacking. Therefore, this study focuses on the functional characterization of A. minutum unigenes obtained from transcriptome sequencing using the Illumina Hiseq 4000 sequencing platform. A total of 58,802 (47.05%) unigenes were successfully annotated using public databases such as NCBI-Nr, UniprotKB, EggNOG, KEGG, InterPRO and Gene Ontology (GO). This study has successfully identified key features that enable A. minutum to adapt to the marine environment, including several carbon metabolic pathways, assimilation of various sources of nitrogen and phosphorus. A. minutum was found to encode homologues for several proteins involved in saxitoxin biosynthesis, including the first three proteins in the pathway of saxitoxin biosynthesis, namely sxtA, sxtG and sxtB. The comprehensive transcriptome analysis presented in this study represents a valuable resource for understanding the dinoflagellates molecular metabolic model regarding nutrient acquisition and biosynthesis of saxitoxin.

Recent developments in quantitative PCR for monitoring harmful marine microalgae

Marine microalgae produce a variety of specialised metabolites that have toxic effects on humans, farmed fish, and marine wildlife. Alarmingly, many of these compounds bioaccumulate in the tissues of shellfish and higher trophic organisms, including species consumed by humans. Molecular methods are emerging as a potential alternative and complement to the conventional microscopic diagnosis of toxic or otherwise harmful microalgal species. Quantitative PCR (qPCR) in particular, has gained popularity over the past decade as a sensitive, rapid, and cost-effective method for monitoring harmful microalgae. Assays targeting taxonomic marker genes provide the opportunity to identify and quantify (or semi-quantify) microalgal species and importantly to pre-empt bloom events. Moreover, the discovery of paralytic shellfish toxin biosynthesis genes in dinoflagellates has enabled researchers to directly monitor toxigenic species in coastal waters and fisheries. This review summarises the recent developments in qPCR detection methods for harmful microalgae, with emphasis on emerging toxin gene monitoring technologies.

From the sxtA4 Gene to Saxitoxin Production: What Controls the Variability Among Alexandrium minutum and Alexandrium pacificum Strains?

Paralytic shellfish poisoning (PSP) is a human foodborne syndrome caused by the consumption of shellfish that accumulate paralytic shellfish toxins (PSTs, saxitoxin group). In PST-producing dinoflagellates such as Alexandrium spp., toxin synthesis is encoded in the nuclear genome via a gene cluster (sxt). Toxin production is supposedly associated with the presence of a 4th domain in the sxtA gene (sxtA4), one of the core genes of the PST gene cluster. It is postulated that gene expression in dinoflagellates is partially constitutive, with both transcriptional and post-transcriptional processes potentially co-occurring. Therefore, gene structure and expression mode are two important features to explore in order to fully understand toxin production processes in dinoflagellates. In this study, we determined the intracellular toxin contents of twenty European Alexandrium minutum and Alexandrium pacificum strains that we compared with their genome size and sxtA4 gene copy numbers. We observed a significant correlation between the sxtA4 gene copy number and toxin content, as well as a moderate positive correlation between the sxtA4 gene copy number and genome size. The 18 toxic strains had several sxtA4 gene copies (9-187), whereas only one copy was found in the two observed non-toxin producing strains. Exploration of allelic frequencies and expression of sxtA4 mRNA in 11 A. minutum strains showed both a differential expression and specific allelic forms in the non-toxic strains compared with the toxic ones. Also, the toxic strains exhibited a polymorphic sxtA4 mRNA sequence between strains and between gene copies within strains. Finally, our study supported the hypothesis of a genetic determinism of toxin synthesis (i.e., the existence of several genetic isoforms of the sxtA4 gene and their copy numbers), and was also consistent with the hypothesis that constitutive gene expression and moderation by transcriptional and post-transcriptional regulation mechanisms are the cause of the observed variability in the production of toxins by A. minutum.

Biosynthesis of Saxitoxin in Marine Dinoflagellates: An Omics Perspective

Saxitoxin is an alkaloid neurotoxin originally isolated from the clam Saxidomus giganteus in 1957. This group of neurotoxins is produced by several species of freshwater cyanobacteria and marine dinoflagellates. The saxitoxin biosynthesis pathway was described for the first time in the 1980s and, since then, it was studied in more than seven cyanobacterial genera, comprising 26 genes that form a cluster ranging from 25.7 kb to 35 kb in sequence length. Due to the complexity of the genomic landscape, saxitoxin biosynthesis in dinoflagellates remains unknown. In order to reveal and understand the dynamics of the activity in such impressive unicellular organisms with a complex genome, a strategy that can carefully engage them in a systems view is necessary. Advances in omics technology (the collective tools of biological sciences) facilitated high-throughput studies of the genome, transcriptome, proteome, and metabolome of dinoflagellates. The omics approach was utilized to address saxitoxin-producing dinoflagellates in response to environmental stresses to improve understanding of dinoflagellates gene–environment interactions. Therefore, in this review, the progress in understanding dinoflagellate saxitoxin biosynthesis using an omics approach is emphasized. Further potential applications of metabolomics and genomics to unravel novel insights into saxitoxin biosynthesis in dinoflagellates are also reviewed.

The Genetic Basis of Toxin Biosynthesis in Dinoflagellates

In marine ecosystems, dinoflagellates can become highly abundant and even dominant at times, despite their comparatively slow growth rates. One factor that may play a role in their ecological success is the production of complex secondary metabolite compounds that can have anti-predator, allelopathic, or other toxic effects on marine organisms, and also cause seafood poisoning in humans. Our knowledge about the genes involved in toxin biosynthesis in dinoflagellates is currently limited due to the complex genomic features of these organisms. Most recently, the sequencing of dinoflagellate transcriptomes has provided us with valuable insights into the biosynthesis of polyketide and alkaloid-based toxin molecules in dinoflagellate species. This review synthesizes the recent progress that has been made in understanding the evolution, biosynthetic pathways, and gene regulation in dinoflagellates with the aid of transcriptomic and other molecular genetic tools, and provides a pathway for future studies of dinoflagellates in this exciting omics era.

Comparative Time-Scale Gene Expression Analysis Highlights the Infection Processes of Two Amoebophrya Strains

Understanding factors that generate, maintain, and constrain host-parasite associations is of major interest to biologists. Although little studied, many extremely virulent micro-eukaryotic parasites infecting microalgae have been reported in the marine plankton. This is the case for Amoebophrya, a diverse and highly widespread group of Syndiniales infecting and potentially controlling dinoflagellate populations. Here, we analyzed the time-scale gene expression of a complete infection cycle of two Amoebophrya strains infecting the same host (the dinoflagellate Scrippsiella acuminata), but diverging by their host range (one infecting a single host, the other infecting more than one species). Over two-thirds of genes showed two-fold differences in expression between at least two sampled stages of the Amoebophrya life cycle. Genes related to carbohydrate metabolism as well as signaling pathways involving proteases and transporters were overexpressed during the free-living stage of the parasitoid. Once inside the host, all genes related to transcription and translation pathways were actively expressed, suggesting the rapid and extensive protein translation needed following host-cell invasion. Finally, genes related to cellular division and components of the flagellum organization were overexpressed during the sporont stage. In order to gain a deeper understanding of the biological basis of the host-parasitoid interaction, we screened proteins involved in host-cell recognition, invasion, and protection against host-defense identified in model apicomplexan parasites. Very few of the genes encoding critical components of the parasitic lifestyle of apicomplexans could be unambiguously identified as highly expressed in Amoebophrya. Genes related to the oxidative stress response were identified as highly expressed in both parasitoid strains. Among them, the correlated expression of superoxide dismutase/ascorbate peroxidase in the specialist parasite was consistent with previous studies on Perkinsus marinus defense. However, this defense process could not be identified in the generalist Amoebophrya strain, suggesting the establishment of different strategies for parasite protection related to host specificity.