Analysis of EST data of the marine protist Oxyrrhis marina, an emerging model for alveolate biology and evolution

Nuclear genome of dinoflagellates: Size variation and insights into evolutionary mechanisms

Recent progress in high-throughput sequencing technologies has dramatically increased availability of genome data for prokaryotes and eukaryotes. Dinoflagellates have distinct chromosomes and a huge genome size, which make their genomic analysis complicated. Here, we reviewed the nuclear genomes of core dinoflagellates, focusing on the genome and cell size. Till now, the genome sizes of several dinoflagellates (more than 25) have been measured by certain methods (e.g., flow cytometry), showing a range of 3-250 pg of genomic DNA per cell. In contrast to their relatively small cell size, their genomes are huge (about 1-80 times the human haploid genome). In the present study, we collected the genome and cell size data of dinoflagellates and compared their relationships. We found that dinoflagellate genome size exhibits a positive correlation with cell size. On the other hand, we recognized that the genome size is not correlated with phylogenetic relatedness. These may be caused by genome duplication, increased gene copy number, repetitive non-coding DNA, transposon expansion, horizontal gene transfer, organelle-to-nucleus gene transfer, and/or mRNA reintegration into the genome. Ultimate verification of these factors as potential causative mechanisms would require sequencing of more dinoflagellate genomes in the future.

Multiple waves of viral invasions in Symbiodiniaceae algal genomes

Dinoflagellates from the family Symbiodiniaceae are phototrophic marine protists that engage in symbiosis with diverse hosts. Their large and distinct genomes are characterized by pervasive gene duplication and large-scale retroposition events. However, little is known about the role and scale of horizontal gene transfer (HGT) in the evolution of this algal family. In other dinoflagellates, high levels of HGTs have been observed, linked to major genomic transitions, such as the appearance of a viral-acquired nucleoprotein that originated via HGT from a large DNA algal virus. Previous work showed that Symbiodiniaceae from different hosts are actively infected by viral groups, such as giant DNA viruses and ssRNA viruses, that may play an important role in coral health. Latent viral infections may also occur, whereby viruses could persist in the cytoplasm or integrate into the host genome as a provirus. This hypothesis received experimental support; however, the cellular localization of putative latent viruses and their taxonomic affiliation are still unknown. In addition, despite the finding of viral sequences in some genomes of Symbiodiniaceae, viral origin, taxonomic breadth, and metabolic potential have not been explored. To address these questions, we searched for putative viral-derived proteins in thirteen Symbiodiniaceae genomes. We found fifty-nine candidate viral-derived HGTs that gave rise to twelve phylogenies across ten genomes. We also describe the taxonomic affiliation of these virus-related sequences, their structure, and their genomic context. These results lead us to propose a model to explain the origin and fate of Symbiodiniaceae viral acquisitions.

DMSP synthesis genes distinguish two types of DMSP producer phenotypes

Dimethylsulfoniopropionate (DMSP) is an important organic carbon and sulfur source in the surface ocean that fuels microbial activity and significantly impacts Earth's climate. After three decades of research, the cellular role(s) of DMSP and environmental drivers of production remain enigmatic. Recent work suggests that cellular DMSP concentrations, and changes in these concentrations in response to environmental stressors, define two major groups of DMSP producers: high DMSP producers that contain ≥ 50 mM intracellular DMSP and low DMSP producers that contain < 50 mM. Here we show that two recently described DMSP synthesis genes (DSYB and TpMT2) may differentiate these two DMSP phenotypes. A survey of prokaryotic and eukaryotic isolates found a significant correlation between the presence of DSYB and TpMT2 genes and previous measurements of high and low DMSP concentrations, respectively. Phylogenetic analysis demonstrated that DSYB and TpMT2 form two distinct clades. DSYB and TpMT2 were also found to be globally abundant in in situ surface communities, and their taxonomic annotations were similar to those observed for isolates. The strong correlation of the DSYB and TpMT2 synthesis genes with high and low producer phenotypes establishes a foundation for direct quantification of DMSP producers, enabling significantly improved predictions of DMSP in situ.

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.

Fugacium Spliced Leader Genes Identified from Stranded RNA-Seq Datasets

Trans-splicing mechanisms have been documented in many lineages that are widely distributed phylogenetically, including dinoflagellates. The spliced leader (SL) sequence itself is conserved in dinoflagellates, although its gene sequences and arrangements have diversified within or across different species. In this study, we present 18 Fugacium kawagutii SL genes identified from stranded RNA-seq reads. These genes typically have a single SL but can contain several partial SLs with lengths ranging from 103 to 292 bp. Unexpectedly, we find the SL gene transcripts contain sequences upstream of the canonical SL, suggesting that generation of mature transcripts will require additional modifications following trans-splicing. We have also identified 13 SL-like genes whose expression levels and length are comparable to Dino-SL genes. Lastly, introns in these genes were identified and a new site for Sm-protein binding was proposed. Overall, this study provides a strategy for fast identification of SL genes and identifies new sequences of F. kawagutii SL genes to supplement our understanding of trans-splicing.

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.

Distinctive Nuclear Features of Dinoflagellates with A Particular Focus on Histone and Histone-Replacement Proteins

Dinoflagellates are important eukaryotic microorganisms that play critical roles as producers and grazers, and cause harmful algal blooms. The unusual nuclei of dinoflagellates "dinokaryon" have led researchers to investigate their enigmatic nuclear features. Their nuclei are unusual in terms of their permanently condensed nucleosome-less chromatin, immense genome, low protein to DNA ratio, guanine-cytosine rich methylated DNA, and unique mitosis process. Furthermore, dinoflagellates are the only known group of eukaryotes that apparently lack histone proteins. Over the course of evolution, dinoflagellates have recruited other proteins, e.g., histone-like proteins (HLPs), from bacteria and dinoflagellates/viral nucleoproteins (DVNPs) from viruses as histone substitutes. Expression diversity of these nucleoproteins has greatly influenced the chromatin structure and gene expression regulation in dinoflagellates. Histone replacement proteins (HLPs and DVNPs) are hypothesized to perform a few similar roles as histone proteins do in other eukaryotes, i.e., gene expression regulation and repairing DNA. However, their role in bulk packaging of DNA is not significant as low amounts of proteins are associated with the gigantic genome. This review intends to summarize the discoveries encompassing unique nuclear features of dinoflagellates, particularly focusing on histone and histone replacement proteins. In addition, a comprehensive view of the evolution of dinoflagellate nuclei is presented.

Dinoflagellates, a Unique Lineage for Retrogene Research

The birth and evolution of retrogenes have played crucial roles in genome evolution. Dinoflagellates represent a unique lineage for retrogene research because the retrogenes can be reliably identified by the presence of a 22 nucleotide splice leader called DinoSL, which is post-transcriptionally added to the 5' terminus of all mRNAs. Compared to studies of retrogenes conducted in other model genomes, dinoflagellate retrogenes can potentially be more comprehensively characterized because intron-containing retrogenes have already been detected. Unfortunately, dinoflagellate retrogene research has long been neglected. Here, we review the work on dinoflagellate retrogenes and show their distinct character. Like the dinoflagellate genome itself, dinoflagellate retrogenes are also characterized by many unusual features, including a high survival rate and large numbers in the genome. These data are critical complements to what we know about retrogenes, and will further frame our understanding of retroposition and its roles in genome evolution, as well as providing new insights into retrogene studies in other genomes.

Comparative Genomics Reveals Two Major Bouts of Gene Retroposition Coinciding with Crucial Periods of Symbiodinium Evolution

Gene retroposition is an important mechanism of genome evolution but the role it plays in dinoflagellates, a critical player in marine ecosystems, is not known. Until recently, when the genomes of two coral-symbiotic dinoflagellate genomes, Symbiodinium kawagutii and S. minutum, were released, it has not been possible to systematically study these retrogenes. Here we examine the abundant retrogenes (∼23% of the total genes) in these species. The hallmark of retrogenes in the genome is the presence of DCCGTAGCCATTTTGGCTCAAG, a spliced leader (DinoSL) constitutively trans-spliced to the 5′-end of all nucleus-encoded mRNAs. Although the retrogenes have often lost part of the 22-nt DinoSL, the putative promoter motif from the DinoSL, TTT(T/G), is consistently retained in the upstream region of these genes, providing an explanation for the high survival rate of retrogenes in dinoflagellates. Our analysis of DinoSL sequence divergence revealed two major bursts of retroposition in the evolutionary history of Symbiodinium, occurring at ∼60 and ∼6 Ma. Reconstruction of the evolutionary trajectory of the Symbiodinium genomes mapped these 2 times to the origin and rapid radiation of this dinoflagellate lineage, respectively. GO analysis revealed differential functional enrichment of the retrogenes between the two episodes, with a broad impact on transport in the first bout and more localized influence on symbiosis-related processes such as cell adhesion in the second bout. This study provides the first evidence of large-scale retroposition as a major mechanism of genome evolution for any organism and sheds light on evolution of coral symbiosis.

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.

Cell Cycle-Dependent Expression Dynamics of G1/S Specific Cyclin, Cellulose Synthase and Cellulase in the Dinoflagellate Prorocentrum donghaiense

Dinoflagellates undergo a typical eukaryotic cell cycle consisting of G1, S, G2, and M phases and some of the typical cell cycle related genes have been computationally identified. However, very few of these genes have been experimentally linked to the cell cycle phases. Besides, although thecate dinoflagellates are known to possess theca composed of cellulose, information on cellulose synthesis and degradation associated with the cell cycle is also limited. In this study, we isolated G1/S cyclin, cellulose synthase and cellulase encoding genes in dinoflagellate Prorocentrum donghaiense. Further, using reverse transcription quantitative PCR (RT-qPCR), we characterized the expression profiles of the three genes throughout the cell cycle. All three showed clear expression dynamics throughout the cell cycle, with fold changes of 26, 2.4 and 9.3 for G1/S cyclin, cellulose synthase and cellulase gene, respectively. The transcript abundance of G1/S cyclin increased in late G1 phase and dropped in early S phase, indicating that this protein is involved in the G1/S transition. Throughout the cell cycle, the average transcript level of cellulose synthase was 4.5-fold higher than that of cellulase. Cellulose synthase and cellulase gene expressions showed peak transcript abundances at middle G1 phase and G2M phase, respectively, indicating the respective roles of these enzymes in the growth of newly divided cells and in cytokinesis. Our results suggest that G1/S cyclin, cellulase, and cellulose synthase genes associated with G1/S transition, G2M, and G1 phases of the cell cycle and are candidates of biomarkers for assessing growth status of P. donghaiense.

Evidence for the retention of two evolutionary distinct plastids in dinoflagellates with diatom endosymbionts

Dinoflagellates harboring diatom endosymbionts (termed “dinotoms”) have undergone a process often referred to as “tertiary endosymbiosis”—the uptake of algae containing secondary plastids and integration of those plastids into the new host. In contrast to other tertiary plastids, and most secondary plastids, the endosymbiont of dinotoms is distinctly less reduced, retaining a number of cellular features, such as their nucleus and mitochondria and others, in addition to their plastid. This has resulted in redundancy between host and endosymbiont, at least between some mitochondrial and cytosolic metabolism, where this has been investigated. The question of plastidial redundancy is particularly interesting as the fate of the host dinoflagellate plastid is unclear. The host cytosol possesses an eyespot that has been postulated to be a remnant of the ancestral peridinin plastid, but this has not been tested, nor has its possible retention of plastid functions. To investigate this possibility, we searched for plastid-associated pathways and functions in transcriptomic data sets from three dinotom species. We show that the dinoflagellate host has indeed retained genes for plastid-associated pathways and that these genes encode targeting peptides similar to those of other dinoflagellate plastid-targeted proteins. Moreover, we also identified one gene encoding an essential component of the dinoflagellate plastid protein import machinery, altogether suggesting the presence of a functioning plastid import system in the host, and by extension a relict plastid. The presence of the same plastid-associated pathways in the endosymbiont also extends the known functional redundancy in dinotoms, further confirming the unusual state of plastid integration in this group of dinoflagellates.

Paralytic shellfish toxin content is related to genomic sxtA4 copy number in Alexandrium minutum strains

Dinoflagellates are microscopic aquatic eukaryotes with huge genomes and an unusual cell regulation. For example, most genes are present in numerous copies and all copies seem to be obligatorily transcribed. The consequence of the gene copy number (CPN) for final protein synthesis is, however, not clear. One such gene is sxtA, the starting gene of paralytic shellfish toxin (PST) synthesis. PSTs are small neurotoxic compounds that can accumulate in the food chain and cause serious poisoning incidences when ingested. They are produced by dinoflagellates of the genera Alexandrium, Gymnodium, and Pyrodinium. Here we investigated if the genomic CPN of sxtA4 is related to PST content in Alexandrium minutum cells. SxtA4 is the 4th domain of the sxtA gene and its presence is essential for PST synthesis in dinoflagellates. We used PST and genome size measurements as well as quantitative PCR to analyze sxtA4 CPN and toxin content in 15 A. minutum strains. Our results show a strong positive correlation between the sxtA4 CPN and the total amount of PST produced in actively growing A. minutum cells. This correlation was independent of the toxin profile produced, as long as the strain contained the genomic domains sxtA1 and sxtA4.

Diversification of the light-harvesting complex gene family via intra- and intergenic duplications in the coral symbiotic alga Symbiodinium

The light-harvesting complex (LHC) is an essential component in light energy capture and transduction to facilitate downstream photosynthetic reactions in plant and algal chloroplasts. The unicellular dinoflagellate alga Symbiodinium is an endosymbiont of cnidarian animals, including corals and sea anemones, and provides carbohydrates generated through photosynthesis to host animals. Although Symbiodinium possesses a unique LHC gene family, called chlorophyll a-chlorophyll c₂-peridinin protein complex (acpPC), its genome-level diversity and evolutionary trajectories have not been investigated. Here, we describe a phylogenetic analysis revealing that many of the LHCs are encoded by highly duplicated genes with multi-subunit polyprotein structures in the nuclear genome of Symbiodinium minutum. This analysis provides an extended list of the LHC gene family in a single organism, including 80 loci encoding polyproteins composed of 145 LHC subunits recovered in the phylogenetic tree. In S. minutum, 5 phylogenetic groups of the Lhcf-type gene family, which is exclusively conserved in algae harboring secondary plastids of red algal origin, were identified. Moreover, 5 groups of the Lhcr-type gene family, of which members are known to be associated with PSI in red algal plastids and secondary plastids of red algal origin, were identified. Notably, members classified within a phylogenetic group of the Lhcf-type (group F1) are highly duplicated, which may explain the presence of an unusually large number of LHC genes in this species. Some gene units were homologous to other units within single loci of the polyprotein genes, whereas intergenic homologies between separate loci were conspicuous in other cases, implying that gene unit ‘shuffling’ by gene conversion and/or genome rearrangement might have been a driving force for diversification. These results suggest that vigorous intra- and intergenic gene duplication events have resulted in the genomic framework of photosynthesis in coral symbiont dinoflagellate algae.