Multifunctional polyketide synthase genes identified by genomic survey of the symbiotic dinoflagellate, Symbiodinium minutum

Transcriptomic analysis of polyketide synthesis in dinoflagellate, Prorocentrum lima

The benthic dinoflagellate Prorocentrum lima is among the most common toxic morphospecies with a cosmopolitan distribution. P. lima can produce polyketide compounds, such as okadaic acid (OA), dinophysistoxin (DTX) and their analogues, which are responsible for diarrhetic shellfish poisoning (DSP). Studying the molecular mechanism of DSP toxin biosynthesis is crucial for understanding the environmental driver influencing toxin biosynthesis as well as for better monitoring of marine ecosystems. Commonly, polyketides are produced by polyketide synthases (PKS). However, no gene has been confirmatively assigned to DSP toxin production. Here, we assembled a transcriptome from 94,730,858 Illumina RNAseq reads using Trinity, resulting in 147,527 unigenes with average sequence length of 1035 nt. Using bioinformatics analysis methods, we found 210 unigenes encoding single-domain PKS with sequence similarity to type I PKSs, as reported in other dinoflagellates. In addition, 15 transcripts encoding multi-domain PKS (forming typical type I PKSs modules) and 5 transcripts encoding hybrid nonribosomal peptide synthetase (NRPS)/PKS were found. Using comparative transcriptome and differential expression analysis, a total of 16 PKS genes were identified to be up-regulated in phosphorus-limited cultures, which was related to the up regulation of toxin expression. In concert with other recent transcriptome analyses, this study contributes to the building consensus that dinoflagellates may utilize a combination of Type I multi-domain and single-domain PKS proteins, in an as yet undefined manner, to synthesize polyketides. Our study provides valuable genomic resource for future research in order to understand the complex mechanism of toxin production in this dinoflagellate.

The Influence of Symbiosis on the Proteome of the Exaiptasia Endosymbiont Breviolum minutum

The cellular mechanisms responsible for the regulation of nutrient exchange, immune response, and symbiont population growth in the cnidarian-dinoflagellate symbiosis are poorly resolved. Here, we employed liquid chromatography-mass spectrometry to elucidate proteomic changes associated with symbiosis in Breviolum minutum, a native symbiont of the sea anemone Exaiptasia diaphana ('Aiptasia'). We manipulated nutrients available to the algae in culture and to the holobiont in hospite (i.e., in symbiosis) and then monitored the impacts of our treatments on host-endosymbiont interactions. Both the symbiotic and nutritional states had significant impacts on the B. minutum proteome. B. minutum in hospite showed an increased abundance of proteins involved in phosphoinositol metabolism (e.g., glycerophosphoinositol permease 1 and phosphatidylinositol phosphatase) relative to the free-living alga, potentially reflecting inter-partner signalling that promotes the stability of the symbiosis. Proteins potentially involved in concentrating and fixing inorganic carbon (e.g., carbonic anhydrase, V-type ATPase) and in the assimilation of nitrogen (e.g., glutamine synthase) were more abundant in free-living B. minutum than in hospite, possibly due to host-facilitated access to inorganic carbon and nitrogen limitation by the host when in hospite. Photosystem proteins increased in abundance at high nutrient levels irrespective of the symbiotic state, as did proteins involved in antioxidant defences (e.g., superoxide dismutase, glutathione s-transferase). Proteins involved in iron metabolism were also affected by the nutritional state, with an increased iron demand and uptake under low nutrient treatments. These results detail the changes in symbiont physiology in response to the host microenvironment and nutrient availability and indicate potential symbiont-driven mechanisms that regulate the cnidarian-dinoflagellate symbiosis.

The Natural Product Domain Seeker version 2 (NaPDoS2) webtool relates ketosynthase phylogeny to biosynthetic function

The Natural Product Domain Seeker (NaPDoS) webtool detects and classifies ketosynthase (KS) and condensation domains from genomic, metagenomic, and amplicon sequence data. Unlike other tools, a phylogeny-based classification scheme is used to make broader predictions about the polyketide synthase (PKS) and nonribosomal peptide synthetase (NRPS) genes in which these domains are found. NaPDoS is particularly useful for the analysis of incomplete biosynthetic genes or gene clusters, as are often observed in poorly assembled genomes and metagenomes, or when loci are not clustered, as in eukaryotic genomes. To help support the growing interest in sequence-based analyses of natural product biosynthetic diversity, here we introduce version 2 of the webtool, NaPDoS2, available at http://napdos.ucsd.edu/napdos2. This update includes the addition of 1417 KS sequences, representing a major expansion of the taxonomic and functional diversity represented in the webtool database. The phylogeny-based KS classification scheme now recognizes 41 class and subclass assignments, including new type II PKS subclasses. Workflow modifications accelerate run times, allowing larger datasets to be analyzed. In addition, default parameters were established using statistical validation tests to maximize KS detection and classification accuracy while minimizing false positives. We further demonstrate the applications of NaPDoS2 to assess PKS biosynthetic potential using genomic, metagenomic, and PCR amplicon datasets. These examples illustrate how NaPDoS2 can be used to predict biosynthetic potential and detect genes involved in the biosynthesis of specific structure classes or new biosynthetic mechanisms.

Expression of fatty acid synthase genes and their role in development and arboviral infection of Aedes aegypti

BACKGROUND

Fatty acids are the building blocks of complex lipids essential for living organisms. In mosquitoes, fatty acids are involved in cell membrane production, energy conservation and expenditure, innate immunity, development and reproduction. Fatty acids are synthesized by a multifunctional enzyme complex called fatty acid synthase (FAS). Several paralogues of FAS were found in the Aedes aegypti mosquito. However, the molecular characteristics and expression of some of these paralogues have not been investigated.

METHODS

Genome assemblies of Ae. aegypti were analyzed, and orthologues of human FAS was identified. Phylogenetic analysis and in silico molecular characterization were performed to identify the functional domains of the Ae. aegypti FAS (AaFAS). Quantitative analysis and loss-of-function experiments were performed to determine the significance of different AaFAS transcripts in various stages of development, expression following different diets and the impact of AaFAS on dengue virus, serotype 2 (DENV2) infection and transmission.

RESULTS

We identified seven putative FAS genes in the Ae. aegypti genome assembly, based on nucleotide similarity to the FAS proteins (tBLASTn) of humans, other mosquitoes and invertebrates. Bioinformatics and molecular analyses suggested that only five of the AaFAS genes produce mRNA and therefore represent complete gene models. Expression levels of AaFAS varied among developmental stages and between male and female Ae. aegypti. Quantitative analyses revealed that expression of AaFAS1, the putative orthologue of the human FAS, was highest in adult females. Transient knockdown (KD) of AaFAS1 did not induce a complete compensation by other AaFAS genes but limited DENV2 infection of Aag2 cells in culture and the midgut of the mosquito.

CONCLUSION

AaFAS1 is the predominant AaFAS in adult mosquitoes. It has the highest amino acid similarity to human FAS and contains all enzymatic domains typical of human FAS. AaFAS1 also facilitated DENV2 replication in both cell culture and in mosquito midguts. Our data suggest that AaFAS1 may play a role in transmission of dengue viruses and could represent a target for intervention strategies.

Full-Length Transcriptome Sequencing of the Scleractinian Coral Montipora foliosa Reveals the Gene Expression Profile of Coral-Zooxanthellae Holobiont

Coral-zooxanthellae holobionts are one of the most productive ecosystems in the ocean. With global warming and ocean acidification, coral ecosystems are facing unprecedented challenges. To save the coral ecosystems, we need to understand the symbiosis of coral-zooxanthellae. Although some Scleractinia (stony corals) transcriptomes have been sequenced, the reliable full-length transcriptome is still lacking due to the short-read length of second-generation sequencing and the uncertainty of the assembly results. Herein, PacBio Sequel II sequencing technology polished with the Illumina RNA-seq platform was used to obtain relatively complete scleractinian coral M. foliosa transcriptome data and to quantify M. foliosa gene expression. A total of 38,365 consensus sequences and 20,751 unique genes were identified. Seven databases were used for the gene function annotation, and 19,972 genes were annotated in at least one database. We found 131 zooxanthellae transcripts and 18,829 M. foliosa transcripts. A total of 6328 lncRNAs, 847 M. foliosa transcription factors (TFs), and 2 zooxanthellae TF were identified. In zooxanthellae we found pathways related to symbiosis, such as photosynthesis and nitrogen metabolism. Pathways related to symbiosis in M. foliosa include oxidative phosphorylation and nitrogen metabolism, etc. We summarized the isoforms and expression level of the symbiont recognition genes. Among the membrane proteins, we found three pathways of glycan biosynthesis, which may be involved in the organic matter storage and monosaccharide stabilization in M. foliosa. Our results provide better material for studying coral symbiosis.

Isolation from a fish kill and transcriptomic characterization of Gyrodinium jinhaense off Long Island Sound

First found in Korean coastal water, the dinoflagellate Gyrodinium jinhaense is a recently established species with unclear global distribution and unexplored genomic characteristics. From a laboratory fish mortality event off Long Island Sound, USA, we isolated a dinoflagellate, and by microscopic and molecular (18S rRNA gene; >99% identical) analyses found that it resembles G. jinhaense, hence named G. jinhaense strain AP17. Towards developing a genetic database for this dinoflagellate, a transcriptome of this species was sequenced using RNA-seq, producing 6 Gbp of data that was assembled into over 70,000 unigenes. The assembled transcriptome GC content was approximately 56% and the total Benchmarking Universal Single-Copy Orthologs for Eukaryota and Alveolata databases were approximately 50% and 57%, respectfully. Genes involved in grazing, energy generation, genome architecture, and protein synthesis, processing, and degradation were highly represented in the transcriptome. Moreover, fragments of polyketide synthase and saxitoxin genes were found but saxitoxins were not detected in high performance liquid chromatography measurements. With the first reported transcriptome for the Gyrodinium genus, this study will serve as a baseline for future Gyrodinium genomics and toxicological studies.

A New Dinoflagellate Genome Illuminates a Conserved Gene Cluster Involved in Sunscreen Biosynthesis

Photosynthetic dinoflagellates of the Family Symbiodiniaceae live symbiotically with many organisms that inhabit coral reefs and are currently classified into fifteen groups, including seven genera. Draft genomes from four genera, Symbiodinium, Breviolum, Fugacium, and Cladocopium, which have been isolated from corals, have been reported. However, no genome is available from the genus Durusdinium, which occupies an intermediate phylogenetic position in the Family Symbiodiniaceae and is well known for thermal tolerance (resistance to bleaching). We sequenced, assembled, and annotated the genome of Durusdinium trenchii, isolated from the coral, Favia speciosa, in Okinawa, Japan. Assembled short reads amounted to 670 Mb with ∼47% GC content. This GC content was intermediate among taxa belonging to the Symbiodiniaceae. Approximately 30,000 protein-coding genes were predicted in the D. trenchii genome, fewer than in other genomes from the Symbiodiniaceae. However, annotations revealed that the D. trenchii genome encodes a cluster of genes for synthesis of mycosporine-like amino acids, which absorb UV radiation. Interestingly, a neighboring gene in the cluster encodes a glucose-methanol-choline oxidoreductase with a flavin adenine dinucleotide domain that is also found in Symbiodinium tridacnidorum. This conservation seems to partially clarify an ancestral genomic structure in the Symbiodiniaceae and its loss in late-branching lineages, including Breviolum and Cladocopium, after splitting from the Durusdinium lineage. Our analysis suggests that approximately half of the taxa in the Symbiodiniaceae may maintain the ability to synthesize mycosporine-like amino acids. Thus, this work provides a significant genomic resource for understanding the genomic diversity of Symbiodiniaceae in corals.

Bioluminescence and toxicity as driving factors in harmful algal blooms: Ecological functions and genetic variability

Dinoflagellates are an ecologically important group of marine microbial eukaryotes with a remarkable array of adaptive strategies. It is ironic that two of the traits for which dinoflagellates are best known, toxin production and bioluminescence, are rarely linked when considering the ecological significance of either. Although dinoflagellate species that form some of the most widespread and frequent harmful algal blooms (HABs) are bioluminescent, the molecular and eco-evolutionary associations between these two traits has received little attention. Here, the major themes of biochemistry and genetics, ecological functions, signaling mechanisms, and evolution are addressed, with parallels and connections drawn between the two. Of the 17 major classes of dinoflagellate toxins, only two are produced by bioluminescent species: saxitoxin (STX) and yessotoxin. Of these, STX has been extensively studied, including the identification of the STX biosynthetic genes. While numerous theories have been put forward as to the eco-evolutionary roles of both bioluminescence and toxicity, a general consensus is that both function as grazing deterrents. Thus, both bioluminescence and toxicity may aid in HAB initiation as they alleviate grazing pressure on the HAB species. A large gap in our understanding is the genetic variability among natural bloom populations, as both toxic and non-toxic strains have been isolated from the same geographic location. The same applies to bioluminescence, as there exist both bioluminescent and non-bioluminescent strains of the same species. Recent evidence demonstrating that blooms are not monoclonal events necessitates a greater level of understanding as to the genetic variability of these traits among sub-populations as well as the mechanisms by which cells acquire or lose the trait, as sequence analysis of STX+ and STX- species indicate the key gene required for toxicity is lost rather than gained. While the extent of genetic variability for both bioluminescence and toxicity among natural HAB sub-populations remains unknown, it is an area that needs to be explored in order to gain greater insights into the molecular mechanisms and environmental parameters driving HAB evolution.

Transcriptome survey and toxin measurements reveal evolutionary modification and loss of saxitoxin biosynthesis genes in the dinoflagellates Amphidinium carterae and Prorocentrum micans

In the present study, we characterized the potential toxin genes for polyketide synthase (PKS) and saxitoxin (STX) biosynthesis using the transcriptomes of two non-STX producing dinoflagellates Amphidinium carterae and Prorocentrum micans. RNA sequencing revealed 94 and 166 PKS contigs in A. carterae and P. micans, respectively. We first detected type III PKS, which was closely related to bacteria. In addition, dozens of homologs of 20 STX biosynthesis genes were identified. Interestingly, the core STX-synthesizing genes sxtA and sxtB were only found in P. micans, whereas sxtD was detected in A. carterae alone. Bioinformatic analysis showed that the first two core genes (sxtA and sxtG) had a low sequence similarity (37.0-67.6%) and different domain organization compared to those of other toxigenic dinoflagellates, such as Alexandrium pacificum. These might result in the breakdown of the initial reactions in STX production and ultimately the loss of the ability to synthesize the toxins in both dinoflagellates. Our findings suggest that toxin-related PKS and sxt genes are commonly found in non-STX producing dinoflagellates. In addition to their involvement in the synthesis of toxins, our result indicates that genes may also have other molecular metabolic functions.

Transcriptomic analysis of polyketide synthases in a highly ciguatoxic dinoflagellate, Gambierdiscus polynesiensis and low toxicity Gambierdiscus pacificus, from French Polynesia

Marine dinoflagellates produce a diversity of polyketide toxins that are accumulated in marine food webs and are responsible for a variety of seafood poisonings. Reef-associated dinoflagellates of the genus Gambierdiscus produce toxins responsible for ciguatera poisoning (CP), which causes over 50,000 cases of illness annually worldwide. The biosynthetic machinery for dinoflagellate polyketides remains poorly understood. Recent transcriptomic and genomic sequencing projects have revealed the presence of Type I modular polyketide synthases in dinoflagellates, as well as a plethora of single domain transcripts with Type I sequence homology. The current transcriptome analysis compares polyketide synthase (PKS) gene transcripts expressed in two species of Gambierdiscus from French Polynesia: a highly toxic ciguatoxin producer, G. polynesiensis, versus a non-ciguatoxic species G. pacificus, each assembled from approximately 180 million Illumina 125 nt reads using Trinity, and compares their PKS content with previously published data from other Gambierdiscus species and more distantly related dinoflagellates. Both modular and single-domain PKS transcripts were present. Single domain β-ketoacyl synthase (KS) transcripts were highly amplified in both species (98 in G. polynesiensis, 99 in G. pacificus), with smaller numbers of standalone acyl transferase (AT), ketoacyl reductase (KR), dehydratase (DH), enoyl reductase (ER), and thioesterase (TE) domains. G. polynesiensis expressed both a larger number of multidomain PKSs, and larger numbers of modules per transcript, than the non-ciguatoxic G. pacificus. The largest PKS transcript in G. polynesiensis encoded a 10,516 aa, 7 module protein, predicted to synthesize part of the polyether backbone. Transcripts and gene models representing portions of this PKS are present in other species, suggesting that its function may be performed in those species by multiple interacting proteins. This study contributes to the building consensus that dinoflagellates utilize a combination of Type I modular and single domain PKS proteins, in an as yet undefined manner, to synthesize polyketides.

Comparative transcriptomics of toxin synthesis genes between the non-toxin producing dinoflagellate Cochlodinium polykrikoides and toxigenic Alexandrium pacificum

In the present study, we extensively characterized potential toxin-related genes, including polyketide synthase (PKS), saxitoxin (STX) and fatty acid synthase (FAS) from the non-toxin producing marine dinoflagellate Cochlodinium polykrikoides, comparing to those of a toxigenic dinoflagellate Alexandrium pacificum. RNA sequencing revealed 50 and 271 PKS contigs from C. polykrikoides and A. pacificum, respectively. According to domain constitute and amino acid alteration, we further classified the dinoflagellate type I PKS genes into 4 sub-groups. Type III PKS was first identified in C. polykrikoides. Interestingly, we detected a large number (242 and 288) of homologs of 18 sxt genes from two studied dinoflagellates. Most of the eight key genes (sxtA, sxtB, sxtD, sxtG, sxtH/T, sxtI, sxtS and sxtU) for STX synthesis were detected in both dinoflatellates, whereas a core STX biosynthesis gene sxtG was not detected in C. polykrikoides. This may partially explain the absence of saxitoxin production in C. polykrikoides. In addition, we identified several type I and type II FAS genes, including FabD, FabF, FabG, FabH, FabI, and FabZ, whereas FabB was not found in C. polykrikoides. Overall, the numbers of the toxin-related genes in C. polykrikoides were less than that of A. pacificum. Phylogenetic analyses showed that type I PKS/FASs of dinoflagellates had close relationships with apicomplexans and bacteria. These suggest that the toxin-related PKS and sxt genes are commonly present in toxigenic and non-toxin producing dinoflagellates, and may be involved not only in the toxin synthesis, but also in other related molecular metabolic functions.

Evolution and Diversity of Assembly-Line Polyketide Synthases

Assembly-line polyketide synthases (PKSs) are among the most complex protein machineries known in nature, responsible for the biosynthesis of numerous compounds used in the clinic. Their present-day diversity is the result of an evolutionary path that has involved the emergence of a multimodular architecture and further diversification of assembly-line PKSs. In this review, we provide an overview of previous studies that investigated PKS evolution and propose a model that challenges the currently prevailing view that gene duplication has played a major role in the emergence of multimodularity. We also analyze the ensemble of orphan PKS clusters sequenced so far to evaluate how large the entire diversity of assembly-line PKS clusters and their chemical products could be. Finally, we examine the existing techniques to access the natural PKS diversity in natural and heterologous hosts and describe approaches to further expand this diversity through engineering.

Research Progress in the Biosynthetic Mechanisms of Marine Polyether Toxins

Marine polyether toxins, mainly produced by marine dinoflagellates, are novel, complex, and diverse natural products with extensive toxicological and pharmacological effects. Owing to their harmful effects during outbreaks of marine red tides, as well as their potential value for the development of new drugs, marine polyether toxins have been extensively studied, in terms of toxicology, pharmacology, detection, and analysis, structural identification, as well as their biosynthetic mechanisms. Although the biosynthetic mechanisms of marine polyether toxins are still unclear, certain progress has been made. In this review, research progress and current knowledge on the biosynthetic mechanisms of polyether toxins are summarized, including the mechanisms of carbon skeleton deletion, pendant alkylation, and polyether ring formation, along with providing a summary of mined biosynthesis-related genes. Finally, future research directions and applications of marine polyether toxins are discussed.

Transcriptomic investigation into polyketide toxin synthesis in Ostreopsis (Dinophyceae) species

In marine ecosystems, dinoflagellates can become highly abundant and even dominant at times, despite their comparatively slow growth. Their ecological success may be related to their production of complex toxic polyketide compounds. Ostreopsis species produce potent palytoxin-like compounds (PLTX), which are associated with human skin and eye irritations, and illnesses through the consumption of contaminated seafood. To investigate the genetic basis of PLTX-like compounds, we sequenced and annotated transcriptomes from two PLTX-producing Ostreopsis species; O. cf. ovata, O. cf. siamensis, one non-PLTX producing species, O. rhodesae and compared them to a close phylogenetic relative and non-PLTX producer, Coolia malayensis. We found no clear differences in the presence or diversity of ketosynthase and ketoreductase transcripts between PLTX producing and non-producing Ostreopsis and Coolia species, as both groups contained >90 and > 10 phylogenetically diverse ketosynthase and ketoreductase transcripts, respectively. We report for the first-time type I single-, multi-domain polyketide synthases (PKSs) and hybrid non-ribosomal peptide synthase/PKS transcripts from all species. The long multi-modular PKSs were insufficient by themselves to synthesize the large complex polyether backbone of PLTX-like compounds. This implies that numerous PKS domains, including both single and multi-, work together on the biosynthesis of PLTX-like and other related polyketide compounds.

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.

Diversified secondary metabolite biosynthesis gene repertoire revealed in symbiotic dinoflagellates

Symbiodiniaceae dinoflagellates possess smaller nuclear genomes than other dinoflagellates and produce structurally specialized, biologically active, secondary metabolites. Till date, little is known about the evolution of secondary metabolism in dinoflagellates as comparative genomic approaches have been hampered by their large genome sizes. Here, we overcome this challenge by combining genomic and metabolomics approaches to investigate how chemical diversity arises in three decoded Symbiodiniaceae genomes (clades A3, B1 and C). Our analyses identify extensive diversification of polyketide synthase and non-ribosomal peptide synthetase genes from two newly decoded genomes of Symbiodinium tridacnidorum (A3) and Cladocopium sp. (C). Phylogenetic analyses indicate that almost all the gene families are derived from lineage-specific gene duplications in all three clades, suggesting divergence for environmental adaptation. Few metabolic pathways are conserved among the three clades and we detect metabolic similarity only in the recently diverged clades, B1 and C. We establish that secondary metabolism protein architecture guides substrate specificity and that gene duplication and domain shuffling have resulted in diversification of secondary metabolism genes.

Two divergent Symbiodinium genomes reveal conservation of a gene cluster for sunscreen biosynthesis and recently lost genes

Background

The marine dinoflagellate, Symbiodinium, is a well-known photosynthetic partner for coral and other diverse, non-photosynthetic hosts in subtropical and tropical shallows, where it comprises an essential component of marine ecosystems. Using molecular phylogenetics, the genus Symbiodinium has been classified into nine major clades, A-I, and one of the reported differences among phenotypes is their capacity to synthesize mycosporine-like amino acids (MAAs), which absorb UV radiation. However, the genetic basis for this difference in synthetic capacity is unknown. To understand genetics underlying Symbiodinium diversity, we report two draft genomes, one from clade A, presumed to have been the earliest branching clade, and the other from clade C, in the terminal branch.

Results

The nuclear genome of Symbiodinium clade A (SymA) has more gene families than that of clade C, with larger numbers of organelle-related genes, including mitochondrial transcription terminal factor (mTERF) and Rubisco. While clade C (SymC) has fewer gene families, it displays specific expansions of repeat domain-containing genes, such as leucine-rich repeats (LRRs) and retrovirus-related dUTPases. Interestingly, the SymA genome encodes a gene cluster for MAA biosynthesis, potentially transferred from an endosymbiotic red alga (probably of bacterial origin), while SymC has completely lost these genes.

Conclusions

Our analysis demonstrates that SymC appears to have evolved by losing gene families, such as the MAA biosynthesis gene cluster. In contrast to the conservation of genes related to photosynthetic ability, the terminal clade has suffered more gene family losses than other clades, suggesting a possible adaptation to symbiosis. Overall, this study implies that Symbiodinium ecology drives acquisition and loss of gene families.

Electronic supplementary material

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

Analysis of the genomic basis of functional diversity in dinoflagellates using a transcriptome-based sequence similarity network

Dinoflagellates are one of the most abundant and functionally diverse groups of eukaryotes. Despite an overall scarcity of genomic information for dinoflagellates, constantly emerging high-throughput sequencing resources can be used to characterize and compare these organisms. We assembled de novo and processed 46 dinoflagellate transcriptomes and used a sequence similarity network (SSN) to compare the underlying genomic basis of functional features within the group. This approach constitutes the most comprehensive picture to date of the genomic potential of dinoflagellates. A core-predicted proteome composed of 252 connected components (CCs) of putative conserved protein domains (pCDs) was identified. Of these, 206 were novel and 16 lacked any functional annotation in public databases. Integration of functional information in our network analyses allowed investigation of pCDs specifically associated with functional traits. With respect to toxicity, sequences homologous to those of proteins found in species with toxicity potential (e.g., sxtA4 and sxtG) were not specific to known toxin-producing species. Although not fully specific to symbiosis, the most represented functions associated with proteins involved in the symbiotic trait were related to membrane processes and ion transport. Overall, our SSN approach led to identification of 45,207 and 90,794 specific and constitutive pCDs of, respectively, the toxic and symbiotic species represented in our analyses. Of these, 56% and 57%, respectively (i.e., 25,393 and 52,193 pCDs), completely lacked annotation in public databases. This stresses the extent of our lack of knowledge, while emphasizing the potential of SSNs to identify candidate pCDs for further functional genomic characterization.

Both modular and single-domain Type I polyketide synthases are expressed in the brevetoxin-producing dinoflagellate, Karenia brevis (Dinophyceae)

Dinoflagellates are prolific producers of polyketide compounds, many of which are potent toxins with adverse impacts on human and marine animal health. To identify polyketide synthase (PKS) genes in the brevetoxin-producing dinoflagellate, Karenia brevis, we assembled a transcriptome from 595 million Illumina reads, sampled under different growth conditions. The assembly included 125,687 transcripts greater than 300 nt in length, with over half having >100× coverage. We found 121 transcripts encoding Type I ketosynthase (KS) domains, of which 99 encoded single KS domains, while 22 contained multiple KS domains arranged in 1-3 protein modules. Phylogenetic analysis placed all single domain and a majority of multidomain KSs within a monophyletic clade of protist PKSs. In contrast with the highly amplified single-domain KSs, only eight single-domain ketoreductase transcripts were found in the assembly, suggesting that they are more evolutionarily conserved. The multidomain PKSs were dominated by trans-acyltransferase architectures, which were recently shown to be prevalent in other algal protists. Karenia brevis also expressed several hybrid nonribosomal peptide synthetase (NRPS)/PKS sequences, including a burA-like sequence previously reported in a wide variety of dinoflagellates. This contrasts with a similarly deep transcriptome of Gambierdiscus polynesiensis, which lacked NRPS/PKS other than the burA-like transcript, and may reflect the presence of amide-containing polyketides in K. brevis and their absence from G. polynesiensis. In concert with other recent transcriptome analyses, this study provides evidence for both single domain and multidomain PKSs in the synthesis of polyketide compounds in dinoflagellates.

Marine microorganisms as a promising and sustainable source of bioactive molecules

There is an urgent need to discover new drug entities due to the increased incidence of severe diseases as cancer and neurodegenerative pathologies, and reducing efficacy of existing antibiotics. Recently, there is a renewed interest in exploring the marine habitat for new pharmaceuticals also thanks to the advancement in cultivation technologies and in molecular biology techniques. Microorganisms represent a still poorly explored resource for drug discovery. The possibility of obtaining a continuous source of bioactives from marine microorganisms, more amenable to culturing compared to macro-organisms, may be able to meet the challenging demands of pharmaceutical industries. This would enable a more environmentally-friendly approach to drug discovery and overcome the over-utilization of marine resources and the use of destructive collection practices. The importance of the topic is underlined by the number of EU projects funded aimed at improving the exploitation of marine organisms for drug discovery.

Role of Modular Polyketide Synthases in the Production of Polyether Ladder Compounds in Ciguatoxin-Producing Gambierdiscus polynesiensis and G. excentricus (Dinophyceae)

Gambierdiscus, a benthic dinoflagellate, produces ciguatoxins that cause the human illness Ciguatera. Ciguatoxins are polyether ladder compounds that have a polyketide origin, indicating that polyketide synthases (PKS) are involved in their production. We sequenced transcriptomes of Gambierdiscus excentricus and Gambierdiscus polynesiensis and found 264 contigs encoding single domain ketoacyl synthases (KS; G. excentricus: 106, G. polynesiensis: 143) and ketoreductases (KR; G. excentricus: 7, G. polynesiensis: 8) with sequence similarity to type I PKSs, as reported in other dinoflagellates. In addition, 24 contigs (G. excentricus: 3, G. polynesiensis: 21) encoding multiple PKS domains (forming typical type I PKSs modules) were found. The proposed structure produced by one of these megasynthases resembles a partial carbon backbone of a polyether ladder compound. Seventeen contigs encoding single domain KS, KR, s-malonyltransacylase, dehydratase and enoyl reductase with sequence similarity to type II fatty acid synthases (FAS) in plants were found. Type I PKS and type II FAS genes were distinguished based on the arrangement of domains on the contigs and their sequence similarity and phylogenetic clustering with known PKS/FAS genes in other organisms. This differentiation of PKS and FAS pathways in Gambierdiscus is important, as it will facilitate approaches to investigating toxin biosynthesis pathways in dinoflagellates.