Dinoflagellate phylogeny revisited: using ribosomal proteins to resolve deep branching dinoflagellate clades

Characterisation and Cultivation of New Lineages of Colponemids, a Critical Assemblage for Inferring Alveolate Evolution

There are several alveolate groups outside the well-studied trio - ciliates, dinoflagellates, and apicomplexans - that are crucial for understanding the evolution of this major taxon. One such assemblage is the "colponemids", which are eukaryotrophic biflagellates, usually with a ventral groove associated with the posterior flagellum. Previous phylogenetic studies show colponemids forming up to three distinct deep branches within alveolates (e.g. sister groups to Myzozoa or all other alveolates). We have developed dieukaryotic (predator-prey) cultures of four colponemid isolates. One represents the first stable culture of the halophile Palustrimonas (feeding on Pharyngomonas), while SSU rDNA phylogenies show the other isolates as two distinct new lineages. Neocolponema saponarium gen. et sp. nov. is a swimming alkaliphile with a large groove, which feeds on a kinetoplastid. Loeffela hirca gen. et sp. nov. is halophilic, has a subtle groove, usually moves along surfaces, and feeds on Pharyngomonas and Percolomonas. Prey capture in both new genera is raptorial, involves a specialized structure/region to the right of the proximal posterior flagellum, and presumed extrusomes. The relationships amongst Myzozoa, ciliates, and the (now) five described colponemid clades are unresolved, signaling that colponemid diversity represents both a challenge and important resource for tracing deep alveolate evolution.

A Comparison of Dinoflagellate Thiolation Domain Binding Proteins Using In Vitro and Molecular Methods

Dinoflagellates play important roles in ecosystems as primary producers and consumers making natural products that can benefit or harm environmental and human health but are also potential therapeutics with unique chemistries. Annotations of dinoflagellate genes have been hampered by large genomes with many gene copies that reduce the reliability of transcriptomics, quantitative PCR, and targeted knockouts. This study aimed to functionally characterize dinoflagellate proteins by testing their interactions through in vitro assays. Specifically, nine Amphidinium carterae thiolation domains that scaffold natural product synthesis were substituted into an indigoidine synthesizing gene from the bacterium Streptomyces lavendulae and exposed to three A. carterae phosphopantetheinyl transferases that activate synthesis. Unsurprisingly, several of the dinoflagellate versions inhibited the ability to synthesize indigoidine despite being successfully phosphopantetheinated. However, all the transferases were able to phosphopantetheinate all the thiolation domains nearly equally, defying the canon that transferases participate in segregated processes via binding specificity. Moreover, two of the transferases were expressed during growth in alternating patterns while the final transferase was only observed as a breakdown product common to all three. The broad substrate recognition and compensatory expression shown here help explain why phosphopantetheinyl transferases are lost throughout dinoflagellate evolution without a loss in a biochemical process.

Vitamin B12-auxotrophy in dinoflagellates caused by incomplete or absent cobalamin-independent methionine synthase genes (metE)

Dinoflagellates are responsible for most marine harmful algal blooms (HABs) and play vital roles in many ocean processes. More than 90% of dinoflagellates are vitamin B12 auxotrophs and that B12 availability can control dinoflagellate HABs, yet the genetic basis of B12 auxotrophy in dinoflagellates in the framework of the ecology of dinoflagellates and particularly HABs, which was the objective of this work. Here, we investigated the presence, phylogeny, and transcription of two methionine synthase genes (B12-dependent metH and B12-independent metE) via searching and assembling transcripts and genes from transcriptomic and genomic databases, cloning 38 cDNA isoforms of the two genes from 14 strains of dinoflagellates, measuring the expression at different scenarios of B12, and comprehensive phylogenetic analyses of more than 100 organisms. We found that 1) metH was present in all 58 dinoflagellates accessible and metE was present in 40 of 58 species, 2) all metE genes lacked N-terminal domains, 3) metE of dinoflagellates were phylogenetically distinct from other known metE genes, and 4) expression of metH in dinoflagellates was responsive to exogenous B12 levels while expression of metE was not responding as that of genuine metE genes. We conclude that most, hypothetically all, dinoflagellates have either non-functional metE genes lacking N-terminal domain for most species, or do not possess metE for other species, which provides the genetic basis for the widespread nature of B12 auxotrophy in dinoflagellates. The work elucidated a fundamental aspect of the nutritional ecology of dinoflagellates.

Integrative omics analysis highlights the immunomodulatory effects of the parasitic dinoflagellate hhematodinium on crustacean hemocytes

Parasitic dinoflagellates in genus Hematodinium have caused substantial economic losses to multiple commercially valuable marine crustaceans around the world. Recent efforts to better understand the life cycle and biology of the parasite have improved our understanding of the disease ecology. However, studies on the host-parasite interaction, especially how Hematodinium parasites evade the host immune response are lacking. To address this shortfall, we used the comprehensive omics approaches (miRNA transcriptomics, iTRAQ-based proteomics) to get insights into the host-parasite interaction between hemocytes from Portunus trituberculatus and Hematodinium perezi in the present study. The parasitic dinoflagellate H. perezi remodeled the miRNome and proteome of hemocytes from challenged hosts, modulated the host immune response at both post-transcriptional and translational levels and caused post-transcriptional regulation to the host immune response. Multiple important cellular and humoral immune-related pathways (ex. Apoptosis, Endocytosis, ECM-receptor interaction, proPO activation pathway, Toll-like signaling pathway, Jak-STAT signaling pathway) were significantly affected by Hematodinium parasites. Through modulation of the host miRNome, the host immune responses of nodulation, proPO activation and antimicrobial peptides were significantly suppressed. Cellular homeostasis was imbalanced via post-transcriptional dysregulation of the phagosome and peroxisome pathways. Cellular structure and communication was seriously impacted by post-transcriptional downregulation of ECM-receptor interaction and focal adhesion pathways. In conclusion, H. perezi parasites could trigger striking changes in the miRNome and proteome of crustacean hemocytes, and this parasite exhibited multifaceted immunomodulatory effects and potential immune-suppressive mechanisms in crustacean hosts.

A Novel Parasitoid of Marine Dinoflagellates, Pararosarium dinoexitiosum gen. et sp. nov. (Perkinsozoa, Alveolata), Showing Characteristic Beaded Sporocytes

The phylum Perkinsozoa is known as an exclusively parasitic group within alveolates and is widely distributed in various aquatic environments from marine to freshwater environments. Nonetheless, their morphology, life cycle, the identity of the host, and physiological characteristics remain still poorly understood. During intensive sampling along the west coast of Korea in October and November 2017, a new parasitoid, which shares several characteristics with the extant families Perkinsidae and Parviluciferaceae, was discovered and three strains of the new parasitoid were successfully established in cultures. Cross-infection experiments showed that among the examined planktonic groups, only dinoflagellates were susceptible to the new parasitoid, with infections observed in species belonging to eight genera. Even though the new parasitoid shared many morphological and developmental characteristics with other Perkinsozoan parasites, it differed from them by its densely packed trophocyte structure without a large vacuole or hyaline material during the growth stage. These characteristics are common among Parviluciferaceae members. Furthermore, through palintomic extracellular sporogenesis, it produced characteristic interconnected sporocytes resembling a string of beads. Phylogenetic analyses based on the small subunit and large subunit ribosomal DNA sequences revealed that the new parasitoid was distantly related to the family Parviluciferaceae and was more closely related to the families Perkinsidae and Xcellidae. Morphological, ultrastructural, and molecular data on the new parasitoid raised the need to erect a new family, i.e., Pararosariidae, within the phylum Perkinsozoa with Pararosarium dinoexitiosum gen. et sp. nov. as the type species. The isolation and establishment in culture of the new parasitoid outside the family Parviluciferaceae in the present study would contribute to the better understanding of the diversity of Perkinsozoan parasites and provide useful material for comparisons to other parasite species in the further study.

A Global Approach to Estimating the Abundance and Duplication of Polyketide Synthase Domains in Dinoflagellates

Many dinoflagellate species make toxins in a myriad of different molecular configurations but the underlying chemistry in all cases is presumably via modular synthases, primarily polyketide synthases. In many organisms modular synthases occur as discrete synthetic genes or domains within a gene that act in coordination thus forming a module that produces a particular fragment of a natural product. The modules usually occur in tandem as gene clusters with a syntenic arrangement that is often predictive of the resultant structure. Dinoflagellate genomes however are notoriously complex with individual genes present in many tandem repeats and very few synthetic modules occurring as gene clusters, unlike what has been seen in bacteria and fungi. However, modular synthesis in all organisms requires a free thiol group that acts as a carrier for sequential synthesis called a thiolation domain. We scanned 47 dinoflagellate transcriptomes for 23 modular synthase domain models and compared their abundance among 10 orders of dinoflagellates as well as their co-occurrence with thiolation domains. The total count of domain types was quite large with over thirty-thousand identified, 29 000 of which were in the core dinoflagellates. Although there were no specific trends in domain abundance associated with types of toxins, there were readily observable lineage specific differences. The Gymnodiniales, makers of long polyketide toxins such as brevetoxin and karlotoxin had a high relative abundance of thiolation domains as well as multiple thiolation domains within a single transcript. Orders such as the Gonyaulacales, makers of small polyketides such as spirolides, had fewer thiolation domains but a relative increase in the number of acyl transferases. Unique to the core dinoflagellates, however, were thiolation domains occurring alongside tetratricopeptide repeats that facilitate protein-protein interactions, especially hexa and hepta-repeats, that may explain the scaffolding required for synthetic complexes capable of making large toxins. Clustering analysis for each type of domain was also used to discern possible origins of duplication for the multitude of single domain transcripts. Single domain transcripts frequently clustered with synonymous domains from multi-domain transcripts such as the BurA and ZmaK like genes as well as the multi-ketosynthase genes, sometimes with a large degree of apparent gene duplication, while fatty acid synthesis genes formed distinct clusters. Surprisingly the acyl-transferases and ketoreductases involved in fatty acid synthesis (FabD and FabG, respectively) were found in very large clusters indicating an unprecedented degree of gene duplication for these genes. These results demonstrate a complex evolutionary history of core dinoflagellate modular synthases with domain specific duplications throughout the lineage as well as clues to how large protein complexes can be assembled to synthesize the largest natural products known.

Dawn of the dinophytes: A first attempt to date origin and diversification of harmful algae

Unicellular dinophytes include a considerable number of harmful algae and potent producers of toxins. The dinophyte fossil record is one of the richest among protists and indicates a geologically old origin of the group dating back to the Triassic. Besides of these records, very few molecular dating studies of dinophytes have been published to date, precluding an understanding of the diversification history of these organisms. In the present study, we used first appearances in the fossil record for the best-represented dinophyte lineages, namely Gonyaulacales and Peridiniales, to calibrate a molecular phylogeny. It is inferred from ribosomal RNA sequence data covering a representative taxon sampling of all currently recognised lineages. Dinophytes may have started diversifying during main tectonic events of the supercontinent Pangaea, witnessing and surviving some of the biggest mass extinction events on Earth. Groups including harmful dinophytes originated at different points in time, but they all predate the Cretaceous-Paleogene boundary. Our chronogram provides a first time frame and may stimulate studies in future bringing molecular phylogenetics of dinophytes and their impressive fossil record together in more detail.

New Perspectives Related to the Bioluminescent System in Dinoflagellates: Pyrocystis lunula, a Case Study

Pyrocystis lunula is considered a model organism due to its bioluminescence capacity linked to circadian rhythms. The mechanisms underlying the bioluminescent phenomenon have been well characterized in dinoflagellates; however, there are still some aspects that remain an enigma. Such is the case of the presence and diversity of the luciferin-binding protein (LBP), as well as the synthesis process of luciferin. Here we carry out a review of the literature in relation to the molecular players responsible for bioluminescence in dinoflagellates, with particular interest in P. lunula. We also carried out a phylogenetic analysis of the conservation of protein sequence, structure and evolutionary pattern of these key players. The basic structure of the luciferase (LCF) is quite conserved among the sequences reported to date for dinoflagellate species, but not in the case of the LBP, which has proven to be more variable in terms of sequence and structure. In the case of luciferin, its synthesis has been shown to be complex process with more than one metabolic pathway involved. The glutathione S-transferase (GST) and the P630 or blue compound, seem to be involved in this process. In the same way, various hypotheses regarding the role of bioluminescence in dinoflagellates are exposed.

Nuclear Gene Transformation in the Dinoflagellate Oxyrrhis marina

The lack of a robust gene transformation tool that allows proper expression of foreign genes and functional testing for the vast number of nuclear genes in dinoflagellates has greatly hampered our understanding of the fundamental biology in this ecologically important and evolutionarily unique lineage of microeukaryotes. Here, we report the development of a dinoflagellate expression vector containing various DNA elements from phylogenetically separate dinoflagellate lineages, an electroporation protocol, and successful expression of introduced genes in an early branching dinoflagellate, Oxyrrhis marina. This protocol, involving the use of Lonza's Nucleofector and a codon-optimized antibiotic resistance gene, has been successfully used to produce consistent results in several independent experiments for O. marina. It is anticipated that this protocol will be adaptable for other dinoflagellates and will allow characterization of many novel dinoflagellate genes.

The Biochemistry and Evolution of the Dinoflagellate Nucleus

Dinoflagellates are known to possess a highly aberrant nucleus-the so-called dinokaryon-that exhibits a multitude of exceptional biological features. These include: (1) Permanently condensed chromosomes; (2) DNA in a cholesteric liquid crystalline state, (3) extremely large DNA content (up to 200 pg); and, perhaps most strikingly, (4) a deficit of histones-the canonical building blocks of all eukaryotic chromatin. Dinoflagellates belong to the Alveolata clade (dinoflagellates, apicomplexans, and ciliates) and, therefore, the biological oddities observed in dinoflagellate nuclei are derived character states. Understanding the sequence of changes that led to the dinokaryon has been difficult in the past with poor resolution of dinoflagellate phylogeny. Moreover, lack of knowledge of their molecular composition has constrained our understanding of the molecular properties of these derived nuclei. However, recent advances in the resolution of the phylogeny of dinoflagellates, particularly of the early branching taxa; the realization that divergent histone genes are present; and the discovery of dinoflagellate-specific nuclear proteins that were acquired early in dinoflagellate evolution have all thrown new light nature and evolution of the dinokaryon.

Lacking catalase, a protistan parasite draws on its photosynthetic ancestry to complete an antioxidant repertoire with ascorbate peroxidase

Background

Antioxidative enzymes contribute to a parasite’s ability to counteract the host’s intracellular killing mechanisms. The facultative intracellular oyster parasite, Perkinsus marinus, a sister taxon to dinoflagellates and apicomplexans, is responsible for mortalities of oysters along the Atlantic coast of North America. Parasite trophozoites enter molluscan hemocytes by subverting the phagocytic response while inhibiting the typical respiratory burst. Because P. marinus lacks catalase, the mechanism(s) by which the parasite evade the toxic effects of hydrogen peroxide had remained unclear. We previously found that P. marinus displays an ascorbate-dependent peroxidase (APX) activity typical of photosynthetic eukaryotes. Like other alveolates, the evolutionary history of P. marinus includes multiple endosymbiotic events. The discovery of APX in P. marinus raised the questions: From which ancestral lineage is this APX derived, and what role does it play in the parasite’s life history?

Results

Purification of P. marinus cytosolic APX activity identified a 32 kDa protein. Amplification of parasite cDNA with oligonucleotides corresponding to peptides of the purified protein revealed two putative APX-encoding genes, designated PmAPX1 and PmAPX2. The predicted proteins are 93% identical, and PmAPX2 carries a 30 amino acid N-terminal extension relative to PmAPX1. The P. marinus APX proteins are similar to predicted APX proteins of dinoflagellates, and they more closely resemble chloroplastic than cytosolic APX enzymes of plants. Immunofluorescence for PmAPX1 and PmAPX2 shows that PmAPX1 is cytoplasmic, while PmAPX2 is localized to the periphery of the central vacuole. Three-dimensional modeling of the predicted proteins shows pronounced differences in surface charge of PmAPX1 and PmAPX2 in the vicinity of the aperture that provides access to the heme and active site.

Conclusions

PmAPX1 and PmAPX2 phylogenetic analysis suggests that they are derived from a plant ancestor. Plant ancestry is further supported by the presence of ascorbate synthesis genes in the P. marinus genome that are similar to those in plants. The localizations and 3D structures of the two APX isoforms suggest that APX fulfills multiple functions in P. marinus within two compartments. The possible role of APX in free-living and parasitic stages of the life history of P. marinus is discussed.

Electronic supplementary material

The online version of this article (10.1186/s12862-019-1465-5) contains supplementary material, which is available to authorized users.

DNA Damage Response Pathways in Dinoflagellates

Dinoflagellates are a general group of phytoplankton, ubiquitous in aquatic environments. Most dinoflagellates are non-obligate autotrophs, subjected to potential physical and chemical DNA-damaging agents, including UV irradiation, in the euphotic zone. Delay of cell cycles by irradiation, as part of DNA damage responses (DDRs), could potentially lead to growth inhibition, contributing to major errors in the estimation of primary productivity and interpretations of photo-inhibition. Their liquid crystalline chromosomes (LCCs) have large amount of abnormal bases, restricted placement of coding sequences at the chromosomes periphery, and tandem repeat-encoded genes. These chromosome characteristics, their large genome sizes, as well as the lack of architectural nucleosomes, likely contribute to possible differential responses to DNA damage agents. In this study, we sought potential dinoflagellate orthologues of eukaryotic DNA damage repair pathways, and the linking pathway with cell-cycle control in three dinoflagellate species. It appeared that major orthologues in photoreactivation, base excision repair, nucleotide excision repair, mismatch repair, double-strand break repair and homologous recombination repair are well represented in dinoflagellate genomes. Future studies should address possible differential DNA damage responses of dinoflagellates over other planktonic groups, especially in relation to possible shift of life-cycle transitions in responses to UV irradiation. This may have a potential role in the persistence of dinoflagellate red tides with the advent of climatic change.

An aerobic eukaryotic parasite with functional mitochondria that likely lacks a mitochondrial genome

Dinoflagellates are microbial eukaryotes that have exceptionally large nuclear genomes; however, their organelle genomes are small and fragmented and contain fewer genes than those of other eukaryotes. The genus Amoebophrya (Syndiniales) comprises endoparasites with high genetic diversity that can infect other dinoflagellates, such as those forming harmful algal blooms (e.g., Alexandrium). We sequenced the genome (~100 Mb) of Amoebophrya ceratii to investigate the early evolution of genomic characters in dinoflagellates. The A. ceratii genome encodes almost all essential biosynthetic pathways for self-sustaining cellular metabolism, suggesting a limited dependency on its host. Although dinoflagellates are thought to have descended from a photosynthetic ancestor, A. ceratii appears to have completely lost its plastid and nearly all genes of plastid origin. Functional mitochondria persist in all life stages of A. ceratii, but we found no evidence for the presence of a mitochondrial genome. Instead, all mitochondrial proteins appear to be lost or encoded in the A. ceratii nucleus.

A precedented nuclear genetic code with all three termination codons reassigned as sense codons in the syndinean Amoebophrya sp. ex Karlodinium veneficum

Amoebophrya is part of an enigmatic, diverse, and ubiquitous marine alveolate lineage known almost entirely from anonymous environmental sequencing. Two cultured Amoebophrya strains grown on core dinoflagellate hosts were used for transcriptome sequencing. BLASTx using different genetic codes suggests that Amoebophyra sp. ex Karlodinium veneficum uses the three typical stop codons (UAA, UAG, and UGA) to encode amino acids. When UAA and UAG are translated as glutamine about half of the alignments have better BLASTx scores, and when UGA is translated as tryptophan one fifth have better scores. However, the sole stop codon appears to be UGA based on conserved genes, suggesting contingent translation of UGA. Neither host sequences, nor sequences from the second strain, Amoebophrya sp. ex Akashiwo sanguinea had similar results in BLASTx searches. A genome survey of Amoebophyra sp. ex K. veneficum showed no evidence for transcript editing aside from mitochondrial transcripts. The dynein heavy chain (DHC) gene family was surveyed and of 14 transcripts only two did not use UAA, UAG, or UGA in a coding context. Overall the transcriptome displayed strong bias for A or U in third codon positions, while the tRNA genome survey showed bias against codons ending in U, particularly for amino acids with two codons ending in either C or U. Together these clues suggest contingent translation mechanisms in Amoebophyra sp. ex K. veneficum and a phylogenetically distinct instance of genetic code modification.

Distribution and Evolution of Peroxisomes in Alveolates (Apicomplexa, Dinoflagellates, Ciliates)

The peroxisome was the last organelle to be discovered and five decades later it is still the Cinderella of eukaryotic compartments. Peroxisomes have a crucial role in the detoxification of reactive oxygen species, the beta-oxidation of fatty acids, and the biosynthesis of etherphospholipids, and they are assumed to be present in virtually all aerobic eukaryotes. Apicomplexan parasites including the malaria and toxoplasmosis agents were described as the first group of mitochondriate protists devoid of peroxisomes. This study was initiated to reassess the distribution and evolution of peroxisomes in the superensemble Alveolata (apicomplexans, dinoflagellates, ciliates). We established transcriptome data from two chromerid algae (Chromera velia, Vitrella brassicaformis), and two dinoflagellates (Prorocentrum minimum, Perkinsus olseni) and identified the complete set of essential peroxins in all four reference species. Our comparative genome analysis provides unequivocal evidence for the presence of peroxisomes in Toxoplasma gondii and related genera. Our working hypothesis of a common peroxisomal origin of all alveolates is supported by phylogenetic analyses of essential markers such as the import receptor Pex5. Vitrella harbors the most comprehensive set of peroxisomal proteins including the catalase and the glyoxylate cycle and it is thus a promising model organism to investigate the functional role of this organelle in Apicomplexa.

Delineating closely related dinoflagellate lineages using phylotranscriptomics

Recently radiated dinoflagellates Apocalathium aciculiferum (collected in Lake Erken, Sweden), Apocalathium malmogiense (Baltic Sea) and Apocalathium aff. malmogiense (Highway Lake, Antarctica) represent a lineage with an unresolved phylogeny. We determined their phylogenetic relationships using phylotranscriptomics based on 792 amino acid sequences. Our results showed that A. aciculiferum diverged from the other two closely related lineages, consistent with their different morphologies in cell size, relative cell length and presence of spines. We hypothesized that A. aff. malmogiense and A. malmogiense, which inhabit different hemispheres, are evolutionarily more closely related because they diverged from a marine common ancestor, adapting to a wide salinity range, while A. aciculiferum colonized a freshwater habitat, by acquiring adaptations to this environment, in particular, salinity intolerance. We show that phylotranscriptomics can resolve the phylogeny of recently diverged protists. This has broad relevance, given that many phytoplankton species are morphologically very similar, and single genes sometimes lack the information to determine species' relationships.

Hv 1 Proton Channels in Dinoflagellates: Not Just for Bioluminescence?

Bioluminescence in dinoflagellates is controlled by H_V1 proton channels. Database searches of dinoflagellate transcriptomes and genomes yielded hits with sequence features diagnostic of all confirmed H_V1, and show that H_V1 is widely distributed in the dinoflagellate phylogeny including the basal species Oxyrrhis marina. Multiple sequence alignments followed by phylogenetic analysis revealed three major subfamilies of H_V1 that do not correlate with presence of theca, autotrophy, geographic location, or bioluminescence. These data suggest that most dinoflagellates express a H_V1 which has a function separate from bioluminescence. Sequence evidence also suggests that dinoflagellates can contain more than one H_V1 gene.

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.

Translation and Translational Control in Dinoflagellates

Dinoflagellates are unicellular protists that feature a multitude of unusual nuclear features, including large genomes, packaging of DNA without histones, and multiple gene copies organized as tandem gene arrays. Furthermore, all dinoflagellate mRNAs experience trans-splicing with a common 22-nucleotide splice leader (SL) sequence. These features challenge some of the concepts and assumptions about the regulation of gene expression derived from work on model eukaryotes such as yeasts and mammals. Translational control in the dinoflagellates, based on extensive study of circadian bioluminescence and by more recent microarray and transcriptome analyses, is now understood to be a crucial element in regulating gene expression. A picture of the translation machinery of dinoflagellates is emerging from the recent availability of transcriptomes of multiple dinoflagellate species and the first complete genome sequences. The components comprising the translational control toolkit of dinoflagellates are beginning to take shape and are outlined here.

dinoref: A curated dinoflagellate (Dinophyceae) reference database for the 18S rRNA gene

Dinoflagellates are a heterogeneous group of protists present in all aquatic ecosystems where they occupy various ecological niches. They play a major role as primary producers, but many species are mixotrophic or heterotrophic. Environmental metabarcoding based on high-throughput sequencing is increasingly applied to assess diversity and abundance of planktonic organisms, and reference databases are definitely needed to taxonomically assign the huge number of sequences. We provide an updated 18S rRNA reference database of dinoflagellates: dinoref. Sequences were downloaded from genbank and filtered based on stringent quality criteria. All sequences were taxonomically curated, classified taking into account classical morphotaxonomic studies and molecular phylogenies, and linked to a series of metadata. dinoref includes 1,671 sequences representing 149 genera and 422 species. The taxonomic assignation of 468 sequences was revised. The largest number of sequences belongs to Gonyaulacales and Suessiales that include toxic and symbiotic species. dinoref provides an opportunity to test the level of taxonomic resolution of different 18S barcode markers based on a large number of sequences and species. As an example, when only the V4 region is considered, 374 of the 422 species included in dinoref can still be unambiguously identified. Clustering the V4 sequences at 98% similarity, a threshold that is commonly applied in metabarcoding studies, resulted in a considerable underestimation of species diversity.

Kingdom Chromista and its eight phyla: a new synthesis emphasising periplastid protein targeting, cytoskeletal and periplastid evolution, and ancient divergences

In 1981 I established kingdom Chromista, distinguished from Plantae because of its more complex chloroplast-associated membrane topology and rigid tubular multipartite ciliary hairs. Plantae originated by converting a cyanobacterium to chloroplasts with Toc/Tic translocons; most evolved cell walls early, thereby losing phagotrophy. Chromists originated by enslaving a phagocytosed red alga, surrounding plastids by two extra membranes, placing them within the endomembrane system, necessitating novel protein import machineries. Early chromists retained phagotrophy, remaining naked and repeatedly reverted to heterotrophy by losing chloroplasts. Therefore, Chromista include secondary phagoheterotrophs (notably ciliates, many dinoflagellates, Opalozoa, Rhizaria, heliozoans) or walled osmotrophs (Pseudofungi, Labyrinthulea), formerly considered protozoa or fungi respectively, plus endoparasites (e.g. Sporozoa) and all chromophyte algae (other dinoflagellates, chromeroids, ochrophytes, haptophytes, cryptophytes). I discuss their origin, evolutionary diversification, and reasons for making chromists one kingdom despite highly divergent cytoskeletons and trophic modes, including improved explanations for periplastid/chloroplast protein targeting, derlin evolution, and ciliary/cytoskeletal diversification. I conjecture that transit-peptide-receptor-mediated 'endocytosis' from periplastid membranes generates periplastid vesicles that fuse with the arguably derlin-translocon-containing periplastid reticulum (putative red algal trans-Golgi network homologue; present in all chromophytes except dinoflagellates). I explain chromist origin from ancestral corticates and neokaryotes, reappraising tertiary symbiogenesis; a chromist cytoskeletal synapomorphy, a bypassing microtubule band dextral to both centrioles, favoured multiple axopodial origins. I revise chromist higher classification by transferring rhizarian subphylum Endomyxa from Cercozoa to Retaria; establishing retarian subphylum Ectoreta for Foraminifera plus Radiozoa, apicomonad subclasses, new dinozoan classes Myzodinea (grouping Colpovora gen. n., Psammosa), Endodinea, Sulcodinea, and subclass Karlodinia; and ranking heterokont Gyrista as phylum not superphylum.

Chromosome structures: reduction of certain problems with unequal gene content and gene paralogs to integer linear programming

Background

Chromosome structure is a very limited model of the genome including the information about its chromosomes such as their linear or circular organization, the order of genes on them, and the DNA strand encoding a gene. Gene lengths, nucleotide composition, and intergenic regions are ignored. Although highly incomplete, such structure can be used in many cases, e.g., to reconstruct phylogeny and evolutionary events, to identify gene synteny, regulatory elements and promoters (considering highly conserved elements), etc. Three problems are considered; all assume unequal gene content and the presence of gene paralogs. The distance problem is to determine the minimum number of operations required to transform one chromosome structure into another and the corresponding transformation itself including the identification of paralogs in two structures. We use the DCJ model which is one of the most studied combinatorial rearrangement models. Double-, sesqui-, and single-operations as well as deletion and insertion of a chromosome region are considered in the model; the single ones comprise cut and join. In the reconstruction problem, a phylogenetic tree with chromosome structures in the leaves is given. It is necessary to assign the structures to inner nodes of the tree to minimize the sum of distances between terminal structures of each edge and to identify the mutual paralogs in a fairly large set of structures. A linear algorithm is known for the distance problem without paralogs, while the presence of paralogs makes it NP-hard. If paralogs are allowed but the insertion and deletion operations are missing (and special constraints are imposed), the reduction of the distance problem to integer linear programming is known. Apparently, the reconstruction problem is NP-hard even in the absence of paralogs. The problem of contigs is to find the optimal arrangements for each given set of contigs, which also includes the mutual identification of paralogs.

Results

We proved that these problems can be reduced to integer linear programming formulations, which allows an algorithm to redefine the problems to implement a very special case of the integer linear programming tool. The results were tested on synthetic and biological samples.

Conclusions

Three well-known problems were reduced to a very special case of integer linear programming, which is a new method of their solutions. Integer linear programming is clearly among the main computational methods and, as generally accepted, is fast on average; in particular, computation systems specifically targeted at it are available. The challenges are to reduce the size of the corresponding integer linear programming formulations and to incorporate a more detailed biological concept in our model of the reconstruction.

Evolutionary Trends of Perkinsozoa (Alveolata) Characters Based on Observations of Two New Genera of Parasitoids of dinoflagellates, Dinovorax gen. nov. and Snorkelia gen. nov

Parasites are one of the ecologically most relevant groups of marine food webs, but their taxonomic and biological complexity hampers the assessment of their diversity and evolutionary trends. Moreover, the within-host processes that govern parasitoid infection, development and reproduction are often unknown. In this study, we describe a new species of a perkinsozoan endoparasitoid that infects the toxic dinoflagellate Dinophysis sacculus, by including observations of its morphology, ultrastructure, life-cycle development and phylogeny. The SSU rDNA sequence and main morphological features were also obtained for a second parasitoid species infecting the bloom-forming dinoflagellate Levanderina fissa. Phylogenetic analyses including the sequences obtained show that all known Perkinsozoa species infecting dinoflagellates cluster together. However, sequences of Parvilucifera prorocentri and those obtained in this study cluster at the base of the clade, while the rest of Parvilucifera representatives form a separated highly-supported cluster. These results, together with differing morphological characters like the formation of a germ-tube, the presence of trichocysts, or the heterochromatin presence in zoospores nucleus justify the erection of Dinovorax pyriformis gen. nov. et sp. nov., and Snorkelia prorocentri gen. nov. et comb. nov. (=Parvilucifera prorocentri). The morphological features and phylogenetic position of these parasitoids represent ancestral characters for the Perkinsozoa phylum, and also for Dinozoa clade, allowing the inference of the evolutionary framework of these Alveolata members.

Use of Antibiotics for Maintenance of Axenic Cultures of Amphidinium carterae for the Analysis of Translation

Most dinoflagellates in culture are bacterized, complicating the quantification of protein synthesis, as well as the analysis of its regulation. In bacterized cultures of Amphidinium carterae Hulbert, up to 80% of protein synthetic activity appears to be predominantly bacterial based on responses to inhibitors of protein synthesis. To circumvent this, axenic cultures of A. carterae were obtained and shown to respond to inhibitors of protein synthesis in a manner characteristic of eukaryotes. However, these responses changed with time in culture correlating with the reappearance of bacteria. Here we show that culture with kanamycin (50 μg/mL), carbenicillin (100 μg/mL), and streptomycin sulfate (50 μg/mL) (KCS), but not 100 units/mL of penicillin and streptomycin (PS), prevents the reappearance of bacteria and allows A. carterae protein synthesis to be quantified without the contribution of an associated bacterial community. We demonstrate that A. carterae can grow in the absence of a bacterial community. Furthermore, maintenance in KCS does not inhibit the growth of A. carterae cultures but slightly extends the growth phase and allows accumulation to somewhat higher saturation densities. We also show that cultures of A. carterae maintained in KCS respond to the eukaryotic protein synthesis inhibitors cycloheximide, emetine, and harringtonine. Establishment of these culture conditions will facilitate our ability to use polysome fractionation and ribosome profiling to study mRNA recruitment. Furthermore, this study shows that a simple and fast appraisal of the presence of a bacterial community in A. carterae cultures can be made by comparing responses to cycloheximide and chloramphenicol rather than depending on lengthier culture-based assessments.

Robust Dinoflagellata phylogeny inferred from public transcriptome databases

Dinoflagellates are dominant members of the plankton and play key roles in ocean ecosystems as primary producers, predators, parasites, coral photobionts, and causative agents of algal blooms that produce toxins harmful to humans and commercial fisheries. These unicellular protists exhibit remarkable trophic and morphological diversity and include species with some of the largest reported nuclear genomes. Despite their high ecological and economic importance, comprehensive genome (or transcriptome) based dinoflagellate trees of life are few in number. To address this issue, we used recently generated public sequencing data, including from the Moore Microbial Eukaryote Transcriptome Sequencing Project, to identify dinoflagellate-specific ortholog groups. These orthologs were combined to create a broadly sampled and highly resolved phylogeny of dinoflagellates. Our results emphasize the scope and utility of public sequencing databases in creating broad and robust phylogenies for large and complex taxonomic lineages, while also providing unique insights into the evolution of thecate dinoflagellates.

Characterization of Acetyl-CoA Carboxylases in the Basal Dinoflagellate Amphidinium carterae

Dinoflagellates make up a diverse array of fatty acids and polyketides. A necessary precursor for their synthesis is malonyl-CoA formed by carboxylating acetyl CoA using the enzyme acetyl-CoA carboxylase (ACC). To date, information on dinoflagellate ACC is limited. Through transcriptome analysis in Amphidinium carterae, we found three full-length homomeric type ACC sequences; no heteromeric type ACC sequences were found. We assigned the putative cellular location for these ACCs based on transit peptide predictions. Using streptavidin Western blotting along with mass spectrometry proteomics, we validated the presence of ACC proteins. Additional bands showing other biotinylated proteins were also observed. Transcript abundance for these ACCs follow the global pattern of expression for dinoflagellate mRNA messages over a diel cycle. This is one of the few descriptions at the transcriptomic and protein level of ACCs in dinoflagellates. This work provides insight into the enzymes which make the CoA precursors needed for fatty acid and toxin synthesis in dinoflagellates.

Transcriptome Analysis of Core Dinoflagellates Reveals a Universal Bias towards "GC" Rich Codons

Although dinoflagellates are a potential source of pharmaceuticals and natural products, the mechanisms for regulating and producing these compounds are largely unknown because of extensive post-transcriptional control of gene expression. One well-documented mechanism for controlling gene expression during translation is codon bias, whereby specific codons slow or even terminate protein synthesis. Approximately 10,000 annotatable genes from fifteen “core” dinoflagellate transcriptomes along a range of overall guanine and cytosine (GC) content were used for codonW analysis to determine the relative synonymous codon usage (RSCU) and the GC content at each codon position. GC bias in the analyzed dataset and at the third codon position varied from 51% and 54% to 66% and 88%, respectively. Codons poor in GC were observed to be universally absent, but bias was most pronounced for codons ending in uracil followed by adenine (UA). GC bias at the third codon position was able to explain low abundance codons as well as the low effective number of codons. Thus, we propose that a bias towards codons rich in GC bases is a universal feature of core dinoflagellates, possibly relating to their unique chromosome structure, and not likely a major mechanism for controlling gene expression.

Identification of a vacuolar proton channel that triggers the bioluminescent flash in dinoflagellates

In 1972, J. Woodland Hastings and colleagues predicted the existence of a proton selective channel (HV1) that opens in response to depolarizing voltage across the vacuole membrane of bioluminescent dinoflagellates and conducts protons into specialized luminescence compartments (scintillons), thereby causing a pH drop that triggers light emission. HV1 channels were subsequently identified and demonstrated to have important functions in a multitude of eukaryotic cells. Here we report a predicted protein from Lingulodinium polyedrum that displays hallmark properties of bona fide HV1, including time-dependent opening with depolarization, perfect proton selectivity, and characteristic ΔpH dependent gating. Western blotting and fluorescence confocal microscopy of isolated L. polyedrum scintillons immunostained with antibody to LpHV1 confirm LpHV1's predicted organellar location. Proteomics analysis demonstrates that isolated scintillon preparations contain peptides that map to LpHV1. Finally, Zn2+ inhibits both LpHV1 proton current and the acid-induced flash in isolated scintillons. These results implicate LpHV1 as the voltage gated proton channel that triggers bioluminescence in L. polyedrum, confirming Hastings' hypothesis. The same channel likely mediates the action potential that communicates the signal along the tonoplast to the scintillon.

Major transitions in dinoflagellate evolution unveiled by phylotranscriptomics

Dinoflagellates are key species in marine environments, but they remain poorly understood in part because of their large, complex genomes, unique molecular biology, and unresolved in-group relationships. We created a taxonomically representative dataset of dinoflagellate transcriptomes and used this to infer a strongly supported phylogeny to map major morphological and molecular transitions in dinoflagellate evolution. Our results show an early-branching position of Noctiluca, monophyly of thecate (plate-bearing) dinoflagellates, and paraphyly of athecate ones. This represents unambiguous phylogenetic evidence for a single origin of the group's cellulosic theca, which we show coincided with a radiation of cellulases implicated in cell division. By integrating dinoflagellate molecular, fossil, and biogeochemical evidence, we propose a revised model for the evolution of thecal tabulations and suggest that the late acquisition of dinosterol in the group is inconsistent with dinoflagellates being the source of this biomarker in pre-Mesozoic strata. Three distantly related, fundamentally nonphotosynthetic dinoflagellates, Noctiluca, Oxyrrhis, and Dinophysis, contain cryptic plastidial metabolisms and lack alternative cytosolic pathways, suggesting that all free-living dinoflagellates are metabolically dependent on plastids. This finding led us to propose general mechanisms of dependency on plastid organelles in eukaryotes that have lost photosynthesis; it also suggests that the evolutionary origin of bioluminescence in nonphotosynthetic dinoflagellates may be linked to plastidic tetrapyrrole biosynthesis. Finally, we use our phylogenetic framework to show that dinoflagellate nuclei have recruited DNA-binding proteins in three distinct evolutionary waves, which included two independent acquisitions of bacterial histone-like proteins.

Investigation into the Physiological Significance of the Phytohormone Abscisic Acid in Perkinsus marinus, an Oyster Parasite Harboring a Nonphotosynthetic Plastid

Some organisms have retained plastids even after they have lost the ability to photosynthesize. Several studies of nonphotosynthetic plastids in apicomplexan parasites have shown that the isopentenyl pyrophosphate biosynthesis pathway in the organelle is essential for their survival. A phytohormone, abscisic acid, one of several compounds biosynthesized from isopentenyl pyrophosphate, regulates the parasite cell cycle. Thus, it is possible that the phytohormone is universally crucial, even in nonphotosynthetic plastids. Here, we examined this possibility using the oyster parasite Perkinsus marinus, which is a plastid‐harboring cousin of apicomplexan parasites and has independently lost photosynthetic ability. Fluridone, an inhibitor of abscisic acid biosynthesis, blocked parasite growth and induced cell clustering. Nevertheless, abscisic acid and its intermediate carotenoids did not affect parasite growth or rescue the parasite from inhibition. Moreover, abscisic acid was not detected from the parasite using liquid chromatography mass spectrometry. Our findings show that abscisic acid does not play any significant roles in P. marinus.

Evolutionary cell biology: functional insight from "endless forms most beautiful"

In animal and fungal model organisms, the complexities of cell biology have been analyzed in exquisite detail and much is known about how these organisms function at the cellular level. However, the model organisms cell biologists generally use include only a tiny fraction of the true diversity of eukaryotic cellular forms. The divergent cellular processes observed in these more distant lineages are still largely unknown in the general scientific community. Despite the relative obscurity of these organisms, comparative studies of them across eukaryotic diversity have had profound implications for our understanding of fundamental cell biology in all species and have revealed the evolution and origins of previously observed cellular processes. In this Perspective, we will discuss the complexity of cell biology found across the eukaryotic tree, and three specific examples of where studies of divergent cell biology have altered our understanding of key functional aspects of mitochondria, plastids, and membrane trafficking.

A method for identification of highly conserved elements and evolutionary analysis of superphylum Alveolata

Background

Perfectly or highly conserved DNA elements were found in vertebrates, invertebrates, and plants by various methods. However, little is known about such elements in protists. The evolutionary distance between apicomplexans can be very high, in particular, due to the positive selection pressure on them. This complicates the identification of highly conserved elements in alveolates, which is overcome by the proposed algorithm.

Results

A novel algorithm is developed to identify highly conserved DNA elements. It is based on the identification of dense subgraphs in a specially built multipartite graph (whose parts correspond to genomes). Specifically, the algorithm does not rely on genome alignments, nor pre-identified perfectly conserved elements; instead, it performs a fast search for pairs of words (in different genomes) of maximum length with the difference below the specified edit distance. Such pair defines an edge whose weight equals the maximum (or total) length of words assigned to its ends. The graph composed of these edges is then compacted by merging some of its edges and vertices. The dense subgraphs are identified by a cellular automaton-like algorithm; each subgraph defines a cluster composed of similar inextensible words from different genomes. Almost all clusters are considered as predicted highly conserved elements. The algorithm is applied to the nuclear genomes of the superphylum Alveolata, and the corresponding phylogenetic tree is built and discussed.

Conclusion

We proposed an algorithm for the identification of highly conserved elements. The multitude of identified elements was used to infer the phylogeny of Alveolata.

Electronic supplementary material

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

Nuclear Architecture and Patterns of Molecular Evolution Are Correlated in the Ciliate Chilodonella uncinata

The relationship between nuclear architecture and patterns of molecular evolution in lineages across the eukaryotic tree of life is not well understood, partly because molecular evolution is traditionally explored as changes in base pairs along a linear sequence without considering the context of nuclear position of chromosomes. The ciliate Chilodonella uncinata is an ideal system to address the relationship between nuclear architecture and patterns of molecular evolution as the somatic macronucleus of this ciliate is composed of a peripheral DNA-rich area (orthomere) and a DNA-poor central region (paramere) to form a “heteromeric” macronucleus. Moreover, because the somatic chromosomes of C. uncinata are highly processed into “gene-sized” chromosomes (i.e., nanochromosomes), we can assess fine-scale relationships between location and sequence evolution. By combining fluorescence microscopy and analyses of transcriptome data from C. uncinata, we find that highly expressed genes have the greatest codon usage bias and are enriched in DNA-poor regions. In contrast, genes with less biased sequences tend to be concentrated in DNA abundant areas, at least during vegetative growth. Our analyses are consistent with recent work in plants and animals where nuclear architecture plays a role in gene expression. At the same time, the unusual localization of nanochromosomes suggests that the highly structured nucleus in C. uncinata may create a “gene bank” that facilitates rapid changes in expression of genes required only in specific life history stages. By using “nonmodel” organisms like C. uncinata, we can explore the universality of eukaryotic features while also providing examples of novel properties (i.e., the presence of a gene bank) that build from these features.

Zoosporic parasites infecting marine diatoms - A black box that needs to be opened

Living organisms in aquatic ecosystems are almost constantly confronted by pathogens. Nevertheless, very little is known about diseases of marine diatoms, the main primary producers of the oceans. Only a few examples of marine diatoms infected by zoosporic parasites are published, yet these studies suggest that diseases may have significant impacts on the ecology of individual diatom hosts and the composition of communities at both the producer and consumer trophic levels of food webs. Here we summarize available ecological and morphological data on chytrids, aphelids, stramenopiles (including oomycetes, labyrinthuloids, and hyphochytrids), parasitic dinoflagellates, cercozoans and phytomyxids, all of which are known zoosporic parasites of marine diatoms. Difficulties in identification of host and pathogen species and possible effects of environmental parameters on the prevalence of zoosporic parasites are discussed. Based on published data, we conclude that zoosporic parasites are much more abundant in marine ecosystems than the available literature reports, and that, at present, both the diversity and the prevalence of such pathogens are underestimated.

Do All Dinoflagellates have an Extranuclear Spindle?

The syndinean dinoflagellates are a diverse assemblage of alveolate endoparasites that branch basal to the core dinoflagellates. Because of their phylogenetic position, the syndineans are considered key model microorganisms in understanding early evolution in the dinoflagellates. Closed mitosis with an extranuclear spindle that traverses the nucleus in cytoplasmic grooves or tunnels is viewed as one of the morphological features shared by syndinean and core dinoflagellates. Here we describe nuclear morphology and mitosis in the syndinean dinoflagellate Amoebophrya sp. from Akashiwo sanguinea, a member of the A. ceratii complex, as revealed by protargol silver impregnation, DNA specific fluorochromes, and transmission electron microscopy. Our observations show that not all species classified as dinoflagellates have an extranuclear spindle. In Amoebophrya sp. from A. sanguinea, an extranuclear microtubule cylinder located in a depression in the nuclear surface during interphase moves into the nucleoplasm via sequential membrane fusion events and develops into an entirely intranuclear spindle. Results suggest that the intranuclear spindle of Amoebophrya spp. may have evolved from an ancestral extranuclear spindle and indicate the need for taxonomic revision of the Amoebophryidae.

Complex Ancestries of Isoprenoid Synthesis in Dinoflagellates

Isoprenoid metabolism occupies a central position in the anabolic metabolism of all living cells. In plastid-bearing organisms, two pathways may be present for de novo isoprenoid synthesis, the cytosolic mevalonate pathway (MVA) and nuclear-encoded, plastid-targeted nonmevalonate pathway (DOXP). Using transcriptomic data we find that dinoflagellates apparently make exclusive use of the DOXP pathway. Using phylogenetic analyses of all DOXP genes we inferred the evolutionary origins of DOXP genes in dinoflagellates. Plastid replacements led to a DOXP pathway of multiple evolutionary origins. Dinoflagellates commonly referred to as dinotoms due to their relatively recent acquisition of a diatom plastid, express two completely redundant DOXP pathways. Dinoflagellates with a tertiary plastid of haptophyte origin, by contrast, express a hybrid pathway of dual evolutionary origin. Here, changes in the targeting motif of signal/transit peptide likely allow for targeting the new plastid by the proteins of core isoprenoid metabolism proteins. Parasitic dinoflagellates of the Amoebophyra species complex appear to have lost the DOXP pathway, suggesting that they may rely on their host for sterol synthesis.

Integration of plastids with their hosts: Lessons learned from dinoflagellates

After their endosymbiotic acquisition, plastids become intimately connected with the biology of their host. For example, genes essential for plastid function may be relocated from the genomes of plastids to the host nucleus, and pathways may evolve within the host to support the plastid. In this review, we consider the different degrees of integration observed in dinoflagellates and their associated plastids, which have been acquired through multiple different endosymbiotic events. Most dinoflagellate species possess plastids that contain the pigment peridinin and show extreme reduction and integration with the host biology. In some species, these plastids have been replaced through serial endosymbiosis with plastids derived from a different phylogenetic derivation, of which some have become intimately connected with the biology of the host whereas others have not. We discuss in particular the evolution of the fucoxanthin-containing dinoflagellates, which have adapted pathways retained from the ancestral peridinin plastid symbiosis for transcript processing in their current, serially acquired plastids. Finally, we consider why such a diversity of different degrees of integration between host and plastid is observed in different dinoflagellates and how dinoflagellates may thus inform our broader understanding of plastid evolution and function.

Cryptic infection of a broad taxonomic and geographic diversity of tadpoles by Perkinsea protists

The decline of amphibian populations, particularly frogs, is often cited as an example in support of the claim that Earth is undergoing its sixth mass extinction event. Amphibians seem to be particularly sensitive to emerging diseases (e.g., fungal and viral pathogens), yet the diversity and geographic distribution of infectious agents are only starting to be investigated. Recent work has linked a previously undescribed protist with mass-mortality events in the United States, in which infected frog tadpoles have an abnormally enlarged yellowish liver filled with protist cells of a presumed parasite. Phylogenetic analyses revealed that this infectious agent was affiliated with the Perkinsea: a parasitic group within the alveolates exemplified by Perkinsus sp., a "marine" protist responsible for mass-mortality events in commercial shellfish populations. Using small subunit (SSU) ribosomal DNA (rDNA) sequencing, we developed a targeted PCR protocol for preferentially sampling a clade of the Perkinsea. We tested this protocol on freshwater environmental DNA, revealing a wide diversity of Perkinsea lineages in these environments. Then, we used the same protocol to test for Perkinsea-like lineages in livers of 182 tadpoles from multiple families of frogs. We identified a distinct Perkinsea clade, encompassing a low level of SSU rDNA variation different from the lineage previously associated with tadpole mass-mortality events. Members of this clade were present in 38 tadpoles sampled from 14 distinct genera/phylogroups, from five countries across three continents. These data provide, to our knowledge, the first evidence that Perkinsea-like protists infect tadpoles across a wide taxonomic range of frogs in tropical and temperate environments, including oceanic islands.

The use of -omic tools in the study of disease processes in marine bivalve mollusks

Our understanding of disease processes and host-pathogen interactions in model species has benefited greatly from the application of medium and high-throughput genomic, metagenomic, epigenomic, transcriptomic, and proteomic analyses. The rate at which new, low-cost, high-throughput -omic technologies are being developed has also led to an expansion in the number of studies aimed at gaining a better understanding of disease processes in bivalves. This review provides a catalogue of the genetic and -omic tools available for bivalve species and examples of how -omics has contributed to the advancement of marine bivalve disease research, with a special focus in the areas of immunity, bivalve-pathogen interactions, mechanisms of disease resistance and pathogen virulence, and disease diagnosis. The analysis of bivalve genomes and transcriptomes has revealed that many immune and stress-related gene families are expanded in the bivalve taxa examined thus far. In addition, the analysis of proteomes confirms that responses to infection are influenced by epigenetic, post-transcriptional, and post-translational modifications. The few studies performed in bivalves show that epigenetic modifications are non-random, suggesting a role for epigenetics in regulating the interactions between bivalves and their environments. Despite the progress -omic tools have enabled in the field of marine bivalve disease processes, there is much more work to be done. To date, only three bivalve genomes have been sequenced completely, with assembly status at different levels of completion. Transcriptome datasets are relatively easy and inexpensive to generate, but their interpretation will benefit greatly from high quality genome assemblies and improved data analysis pipelines. Finally, metagenomic, epigenomic, proteomic, and metabolomic studies focused on bivalve disease processes are currently limited but their expansion should be facilitated as more transcriptome datasets and complete genome sequences become available for marine bivalve species.

Endosymbiosis undone by stepwise elimination of the plastid in a parasitic dinoflagellate

Organelle gain through endosymbiosis has been integral to the origin and diversification of eukaryotes, and, once gained, plastids and mitochondria seem seldom lost. Indeed, discovery of nonphotosynthetic plastids in many eukaryotes--notably, the apicoplast in apicomplexan parasites such as the malaria pathogen Plasmodium--highlights the essential metabolic functions performed by plastids beyond photosynthesis. Once a cell becomes reliant on these ancillary functions, organelle dependence is apparently difficult to overcome. Previous examples of endosymbiotic organelle loss (either mitochondria or plastids), which have been invoked to explain the origin of eukaryotic diversity, have subsequently been recognized as organelle reduction to cryptic forms, such as mitosomes and apicoplasts. Integration of these ancient symbionts with their hosts has been too well developed to reverse. Here, we provide evidence that the dinoflagellate Hematodinium sp., a marine parasite of crustaceans, represents a rare case of endosymbiotic organelle loss by the elimination of the plastid. Extensive RNA and genomic sequencing data provide no evidence for a plastid organelle, but, rather, reveal a metabolic decoupling from known plastid functions that typically impede organelle loss. This independence has been achieved through retention of ancestral anabolic pathways, enzyme relocation from the plastid to the cytosol, and metabolic scavenging from the parasite's host. Hematodinium sp. thus represents a further dimension of endosymbiosis--life after the organelle.

The alveolate translation initiation factor 4E family reveals a custom toolkit for translational control in core dinoflagellates

Background

Dinoflagellates are eukaryotes with unusual cell biology and appear to rely on translational rather than transcriptional control of gene expression. The eukaryotic translation initiation factor 4E (eIF4E) plays an important role in regulating gene expression because eIF4E binding to the mRNA cap is a control point for translation. eIF4E is part of an extended, eukaryote-specific family with different members having specific functions, based on studies of model organisms. Dinoflagellate eIF4E diversity could provide a mechanism for dinoflagellates to regulate gene expression in a post-transcriptional manner. Accordingly, eIF4E family members from eleven core dinoflagellate transcriptomes were surveyed to determine the diversity and phylogeny of the eIF4E family in dinoflagellates and related lineages including apicomplexans, ciliates and heterokonts.

Results

The survey uncovered eight to fifteen (on average eleven) different eIF4E family members in each core dinoflagellate species. The eIF4E family members from heterokonts and dinoflagellates segregated into three clades, suggesting at least three eIF4E cognates were present in their common ancestor. However, these three clades are distinct from the three previously described eIF4E classes, reflecting diverse approaches to a central eukaryotic function. Heterokonts contain four clades, ciliates two and apicomplexans only a single recognizable eIF4E clade. In the core dinoflagellates, the three clades were further divided into nine sub-clades based on the phylogenetic analysis and species representation. Six of the sub-clades included at least one member from all eleven core dinoflagellate species, suggesting duplication in their shared ancestor. Conservation within sub-clades varied, suggesting different selection pressures.

Conclusions

Phylogenetic analysis of eIF4E in core dinoflagellates revealed complex layering of duplication and conservation when compared to other eukaryotes. Our results suggest that the diverse eIF4E family in core dinoflagellates may provide a toolkit to enable selective translation as a strategy for controlling gene expression in these enigmatic eukaryotes.

Electronic supplementary material

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

De novo assembly and characterization of the transcriptome of the toxic dinoflagellate Karenia brevis

Background

Karenia brevis is a harmful algal species that blooms in the Gulf of Mexico and produces brevetoxins that cause neurotoxic shellfish poisoning. Elevated brevetoxin levels in K. brevis cells have been measured during laboratory hypo-osmotic stress treatments. To investigate mechanisms underlying K. brevis osmoacclimation and osmoregulation and establish a valuable resource for gene discovery, we assembled reference transcriptomes for three clones: Wilson-CCFWC268, SP3, and SP1 (a low-toxin producing variant). K. brevis transcriptomes were annotated with gene ontology terms and searched for putative transmembrane proteins that may elucidate cellular responses to hypo-osmotic stress. An analysis of single nucleotide polymorphisms among clones was used to characterize genetic divergence.

Results

K. brevis reference transcriptomes were assembled with 58.5 (Wilson), 78.0 (SP1), and 51.4 million (SP3) paired reads. Transcriptomes contained 86,580 (Wilson), 93,668 (SP1), and 84,309 (SP3) predicted transcripts. Approximately 40% of the transcripts were homologous to proteins in the BLAST nr database with an E value ≤ 1.00E-6. Greater than 80% of the highly conserved CEGMA core eukaryotic genes were identified in each transcriptome, which supports assembly completeness. Seven putative voltage-gated Na⁺ or Ca²⁺ channels, two aquaporin-like proteins, and twelve putative VATPase subunits were discovered in all clones using multiple bioinformatics approaches. Furthermore, 45% (Wilson) and 43% (SP1 and SP3) of the K. brevis putative peptides > 100 amino acids long produced significant hits to a sequence in the NCBI nr protein database. Of these, 77% (Wilson and SP1) and 73% (SP3) were successfully assigned gene ontology functional terms. The predicted single nucleotide polymorphism (SNP) frequencies between clones were 0.0028 (Wilson to SP1), 0.0030 (Wilson to SP3), and 0.0028 (SP1 to SP3).

Conclusions

The K. brevis transcriptomes assembled here provide a foundational resource for gene discovery and future RNA-seq experiments. The identification of ion channels, VATPases, and aquaporins in all three transcriptomes indicates that K. brevis regulates cellular ion and water concentrations via transmembrane proteins. Additionally, > 40,000 unannotated loci may include potentially novel K. brevis genes. Ultimately, the SNPs identified among the three ecologically diverse clones with different toxin profiles may help to elucidate variations in K. brevis brevetoxin production.

Electronic supplementary material

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

Genomic insights into processes driving the infection of Alexandrium tamarense by the Parasitoid Amoebophrya sp

The regulatory circuits during infection of dinoflagellates by their parasites are largely unknown on the molecular level. Here we provide molecular insights into these infection dynamics. Alexandrium tamarense is one of the most prominent harmful algal bloom dinoflagellates. Its pathogen, the dinoflagellate parasitoid Amoebophrya sp., has been observed to infect and control the blooms of this species. We generated a data set of transcripts from three time points (0, 6, and 96 h) during the infection of this parasite-host system. Assembly of all transcript data from the parasitoid (>900,000 reads/313 Mbp with 454/Roche next-generation sequencing [NGS]) yielded 14,455 contigs, to which we mapped the raw transcript reads of each time point of the infection cycle. We show that particular surface lectins are expressed at the beginning of the infection cycle which likely mediate the attachment to the host cell. In a later phase, signal transduction-related genes together with transmembrane transport and cytoskeleton proteins point to a high integration of processes involved in host recognition, adhesion, and invasion. At the final maturation stage, cell division- and proliferation-related genes were highly expressed, reflecting the fast cell growth and nuclear division of the parasitoid. Our molecular insights into dinoflagellate parasitoid interactions point to general mechanisms also known from other eukaryotic parasites, especially from the Alveolata. These similarities indicate the presence of fundamental processes of parasitoid infection that have remained stable throughout evolution within different phyla.

Metatranscriptome profiling of a harmful algal bloom

Metagenomic methods provide a powerful means to investigate complex ecological phenomena. Developed originally for study of Bacteria and Archaea, the application of these methods to eukaryotic microorganisms is yet to be fully realized. Most prior environmental molecular studies of eukaryotes have relied heavily on PCR amplification with eukaryote-specific primers. Here we apply high throughput short-read sequencing of poly-A selected RNA to capture the metatranscriptome of an estuarine dinoflagellate bloom. To validate the metatranscriptome assembly process we simulated metatranscriptomic datasets using short-read sequencing data from clonal cultures of four algae of varying phylogenetic distance. We find that the proportion of chimeric transcripts reconstructed from community transcriptome sequencing is low, suggesting that metatranscriptomic sequencing can be used to accurately reconstruct the transcripts expressed by bloom-forming communities of eukaryotes. To further validate the bloom metatransciptome assembly we compared it to a transcriptomic assembly from a cultured, clonal isolate of the dominant bloom-causing alga and found that the two assemblies are highly similar. Eukaryote-wide phylogenetic analyses reveal the taxonomic composition of the bloom community, which is comprised of several dinoflagellates, ciliates, animals, and fungi. The assembled metatranscriptome reveals the functional genomic composition of a metabolically active community. Highlighting the potential power of these methods, we found that relative transcript abundance patterns suggest that the dominant dinoflagellate might be expressing toxin biosynthesis related genes at a higher level in the presence of competitors, predators and prey compared to it growing in monoculture.

Phylogeny of a relapsing fever Borrelia species transmitted by the hard tick Ixodes scapularis

The discovery of Borrelia species that were related to the agents of relapsing fever but were transmitted by hard ticks rather than soft ticks challenged previous taxonomies based largely on microbe-host specificities and geographic considerations. One of these newly-identified organisms is the Borrelia miyamotoi sensu lato strain LB-2001 from North America and transmitted by Ixodes scapularis. This or related strains have been identified as the cause of human disease, but comparatively little is known about their biology or genetics. Using recently acquired chromosome sequence of LB-2001 together with database sequences and additional sequences determined here, I carried out comparisons of the several species of Borrelia, including those in the two major clades: the relapsing fever group of species and the Lyme disease group of species. Phylogenetic inference at the species level was based on four data sets: whole chromosomes of ∼1Mb each, and concatenated sequences of 19 ribosomal protein genes, 3 conserved nucleic acid enzymes (rpoC, recC, and dnaE), and 4 contiguous genes for nucleotide salvage on a large plasmid. Analyses using neighbor-joining, maximum likelihood, and Bayesian methods were largely concordant for each of the trees. They showed that LB-2001 and related hard tick-associated organisms, like Borrelia lonestari, are deeply positioned within the RF group of species and that these organisms did not, as some earlier estimations had suggested, constitute a paraphyletic group. The analyses also provided further evidence that major changes in host ranges and life cycles, such as hard to soft ticks or vice versa, may not correlate well with overall sequence differences. The genetic differences between LB-2001 and B. miyamotoi sensu stricto justify provisional use of the "sensu lato" designation for LB-2001.