Proof that Dinoflagellate Spliced Leader (DinoSL) is a Useful Hook for Fishing Dinoflagellate Transcripts from Mixed Microbial Samples: Symbiodinium kawagutii as a Case Study

Broad active metabolic pathways, autophagy, and antagonistic hormones regulate dinoflagellate cyst dormancy in marine sediments

This work aimed to reveal the molecular machinery regulating the dormancy of dinoflagellate resting cysts buried in marine sediments. Dinoflagellates play pivotal roles in marine ecosystems, particularly as major contributors of harmful algal blooms. Despite vital roles of cysts in blooming cycles and dinoflagellate ecology, the molecular processes controlling cyst dormancy have largely remained unexplored due to technological difficulties. Using DinoSL as a dinoflagellates-specific mRNA "hook" and SMRT sequencing, we analyzed metatranscriptomes of sediment-buried dinoflagellate cyst assemblages. The data show that most major metabolic and regulatory pathways, except photosynthesis, were transcriptionally active. This suggests the crucial importance of broad metabolic pathways in sustaining cyst viability and germination potential. Further expression analyses of 11 genes (relevant to autophagy and phytohormone gibberellin), lysosome/autolysosome staining, and germination experiments revealed vital roles of autophagy in energy generation, nutrient recycling, and of phytohormones abscisic acid/gibberellin in modulating dormancy/germination of resting cysts. Our findings lay a cornerstone for elucidating the molecular machinery regulating dinoflagellate cyst dormancy.

Seasonal dynamics and driving mechanisms of microbial biogenic elements cycling function, assembly process, and co-occurrence network in plateau lake sediments

Microbial community diversity significantly varies with seasonality. However, little is known about seasonal variation of microbial community functions in lake sediments and their associated environmental influences. In this study, metagenomic sequencing of sediments collected from winter, summer, and autumn from Caohai Lake, Guizhou Plateau, were used to evaluate the composition and function of sediment microbial communities, the potential interactions of functional genes, key genes associated with seasons, and community assembly mechanisms. The average concentrations of nitrogen (TN) and phosphorus (TP) in lake sediments were higher, which were 6.136 and 0.501 g/kg, respectively. TN and organic matter (OM) were the primary factors associated with sediment community composition and functional profiles. The diversity and structure of the microbial communities varied with seasons, and Proteobacteria relative abundances were significantly lower in summer than in other seasons (58.43-44.12 %). Seasons were also associated with the relative abundances of functional genes, and in particular korA, metF, narC, nrfA, pstC/S, and soxB genes. Network complexity was highest in the summer and key genes in the network also varied across seasons. Neutral community model analysis revealed that the assembly mechanisms related to carbon (C), nitrogen (N), phosphorus (P), and sulfur (S) cycle-related genes were primarily associated with random processes. In summary, diverse functional genes were identified in lake sediments and exhibited evidence for synergistic interactions (Positive proportion: 74.91-99.82 %), while seasonal factors influenced their distribution. The results of this study provide new insights into seasonal impacts on microbial-driven biogeochemical cycling in shallow lakes.

Endogenous viral elements reveal associations between a non-retroviral RNA virus and symbiotic dinoflagellate genomes

Endogenous viral elements (EVEs) offer insight into the evolutionary histories and hosts of contemporary viruses. This study leveraged DNA metagenomics and genomics to detect and infer the host of a non-retroviral dinoflagellate-infecting +ssRNA virus (dinoRNAV) common in coral reefs. As part of the Tara Pacific Expedition, this study surveyed 269 newly sequenced cnidarians and their resident symbiotic dinoflagellates (Symbiodiniaceae), associated metabarcodes, and publicly available metagenomes, revealing 178 dinoRNAV EVEs, predominantly among hydrocoral-dinoflagellate metagenomes. Putative associations between Symbiodiniaceae and dinoRNAV EVEs were corroborated by the characterization of dinoRNAV-like sequences in 17 of 18 scaffold-scale and one chromosome-scale dinoflagellate genome assembly, flanked by characteristically cellular sequences and in proximity to retroelements, suggesting potential mechanisms of integration. EVEs were not detected in dinoflagellate-free (aposymbiotic) cnidarian genome assemblies, including stony corals, hydrocorals, jellyfish, or seawater. The pervasive nature of dinoRNAV EVEs within dinoflagellate genomes (especially Symbiodinium), as well as their inconsistent within-genome distribution and fragmented nature, suggest ancestral or recurrent integration of this virus with variable conservation. Broadly, these findings illustrate how +ssRNA viruses may obscure their genomes as members of nested symbioses, with implications for host evolution, exaptation, and immunity in the context of reef health and disease.

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.

Possible functions of CobW domain-containing (CBWD) genes in dinoflagellates using Karlodinium veneficum as a representative

Since > 91% of dinoflagellates are proven auxotrophs of vitamin B₁₂ and the cobalamin synthetase W (CobW) is a key gene involved in vitamin B₁₂ synthesis pathway, a number of CobW domain-containing (CBWD) genes in dinoflagellates (DinoCBWDs) were surprisedly found from our transcriptomic and meta-transcriptomic studies. A total of 88 DinoCBWD genes were identified from the genomes and transcriptomes of four dinoflagellates, with five being cloned for full-lengths and characterized using the cosmopolitan and ecologically-important dinoflagellates Karlodinium veneficum and Scrippsiella trochoidea (synonym of Scrippsiella acuminata). DinoCBWDs were verified being irrelevant to vitamin B₁₂ biosynthesis due to their transcriptions irresponsive to vitamin B₁₂ levels and their phylogenetic positions. A comprehensive phylogenetic analysis demonstrated 75 out of the 88 DinoCBWD genes identified belong to three subfamilies of COG0523 protein family, of which most prokaryotic members are reported to be metallochaperones and the eukaryotic members are ubiquitously found but mostly unknown for their functions. Our results from K. veneficum demonstrated DinoCBWDs are associated with metal homeostasis and other divergent functions, with four KvCBWDs involving in zinc homeostasis and KvCBWD1 likely functioning as Fe-type nitrile hydratase activator. In addition, conserved motif analysis revealed the structural foundation of KvCBWD proteins that are consistent with previously described CBWD proteins with GTPase activity and metal binding. Our results provide a stepping-stone toward better understanding the functions of DinoCBWDs and the COG0523 family.

Hypothesis: Trans-splicing Generates Evolutionary Novelty in the Photosynthetic Amoeba Paulinella

Plastid primary endosymbiosis has occurred twice, once in the Archaeplastida ancestor and once in the Paulinella (Rhizaria) lineage. Both events precipitated massive evolutionary changes, including the recruitment and activation of genes that are horizontally acquired (HGT) and the redeployment of existing genes and pathways in novel contexts. Here we address the latter aspect in Paulinella micropora KR01 (hereafter, KR01) that has independently evolved spliced leader (SL) trans-splicing (SLTS) of nuclear-derived transcripts. We investigated the role of this process in gene regulation, novel gene origination, and endosymbiont integration. Our analysis shows that 20% of KR01 genes give rise to transcripts with at least one (but in some cases, multiple) sites of SL addition. This process, which often occurs at canonical cis-splicing acceptor sites (internal introns), results in shorter transcripts that may produce 5'-truncated proteins with novel functions. SL-truncated transcripts fall into four categories that may show: (i) altered protein localization, (ii) altered protein function, structure, or regulation, (iii) loss of valid alternative start codons, preventing translation, or (iv) multiple SL addition sites at the 5'-terminus. The SL RNA genes required for SLTS are putatively absent in the heterotrophic sister lineage of photosynthetic Paulinella species. Moreover, a high proportion of transcripts derived from genes of endosymbiotic gene transfer (EGT) and HGT origin contain SL sequences. We hypothesize that truncation of transcripts by SL addition may facilitate the generation and expression of novel gene variants and that SLTS may have enhanced the activation and fixation of foreign genes in the host genome of the photosynthetic lineages, playing a key role in primary endosymbiont integration.

The Implication Inferred from the Expression of Small Heat-Shock Protein Genes in Dinoflagellate Resting Cysts Buried in Marine Sediment

Dinoflagellates are unicellular eukaryotic microalgae, occupying pivotal niches in aquatic ecosystems with great ecological, biological, and economic significance. Small heat shock proteins (sHsps) are the most omnipresent, but the least conserved, family of molecular chaperones found in all domains of life. Although their common name (small Hsp) implies to exclusively stress their heat shock-responsive function, many sHsps in fact engage in a variety of physiological processes, from cell growth and proliferation to embryogenesis, development, differentiation, apoptosis, and even to human disease prevention. Recent years have greatly expanded our understanding of sHsps in higher plants; however, comprehensive study aiming to delineate the composition and expression pattern of dinoflagellate sHsp gene family has not yet been performed. In this study, we constructed dinoflagellate-specific environmental cDNA library from marine sediment and sequenced using the third-generation sequencing technique. Screening of sHsp genes from the library returned 13 entries with complete coding regions, which were considered to be transcriptionally activated in the natural community of dinoflagellate resting cysts. All the 13 dinoflagellate sHsps consisted of a solely characteristic α-crystallin domain, covering 88–123 amino acid residues with the typical A-X-X-X-N-G-V-L motif, flanked by variable N- and C-terminal extensions. Multiple alignment revealed considerable amino acid divergence (~26.7% average similarity) among them. An unexpected close relationship was revealed between dinoflagellate and green algal sHsps in the phylogenetic tree, seemingly reflecting a close evolutionary relationship of these sHsps themselves. We confirmed that sHsp mRNAs are expressed during dormancy of the resting cyst assemblages of dinoflagellates that were buried in marine sediment, which raised the possibility that the sHsp expression is part of the machinery of maintaining the dormancy or/and the adaptation to ambient conditions of dinoflagellate resting cysts. Our results, although preliminary, gained an important glance on the universal presence of sHsps in dinoflagellates and their active expressions in the assemblage of resting cysts that were buried in the marine sediment. The essentiality of sHsps functioning in resting cysts necessitate more intensive and extensive investigations on all possible functions of Hsps in dinoflagellates, a group of protists with vital ecological and biological importance.

Screening a Spliced Leader-Based Symbiodinium microadriaticum cDNA Library Using the Yeast-Two Hybrid System Reveals a Hemerythrin-Like Protein as a Putative SmicRACK1 Ligand

The dinoflagellate Symbiodiniaceae family plays a central role in the health of the coral reef ecosystem via the symbiosis that establishes with its inhabiting cnidarians and supports the host metabolism. In the last few decades, coral reefs have been threatened by pollution and rising temperatures which have led to coral loss. These events have raised interest in studying Symbiodiniaceae and their hosts; however, progress in understanding their metabolism, signal transduction pathways, and physiology in general, has been slow because dinoflagellates present peculiar characteristics. We took advantage of one of these peculiarities; namely, the post-transcriptional addition of a Dino Spliced Leader (Dino-SL) to the 5' end of the nuclear mRNAs, and used it to generate cDNA libraries from Symbiodinium microadriaticum. We compared sequences from two Yeast-Two Hybrid System cDNA Libraries, one based on the Dino-SL sequence, and the other based on the SMART technology (Switching Mechanism at 5' end of RNA Transcript) which exploits the template switching function of the reverse transcriptase. Upon comparison of the performance of both libraries, we obtained a significantly higher yield, number and length of sequences, number of transcripts, and better 5' representation from the Dino-SL based library than from the SMART library. In addition, we confirmed that the cDNAs from the Dino-SL library were adequately expressed in the yeast cells used for the Yeast-Two Hybrid System which resulted in successful screening for putative SmicRACK1 ligands, which yielded a putative hemerythrin-like protein.

Functional profiles of phycospheric microorganisms during a marine dinoflagellate bloom

Harmful algal blooms (HABs) are an ecological concern but relatively few studies have investigated the functional potential of bacterioplankton over a complete algal bloom cycle, which is critical for determining their contribution to the fate of algal blooms. To address this point, we carried out a time-series metagenomic analysis of the functional features of microbial communities at three different Gymnodinium catenatum bloom stages (pre-, peak-, and post-bloom). Different microbial composition were observed during the blooming stages. The environmental parameters and correlation networks co-contribute to microbial variability, and the former explained 38.4% of total variations of the bacterioplankton community composition. Functionally, a range of pathways involved in carbon, nitrogen, phosphorus and sulfur cycling were significantly different during the various HAB stages. Genes associated with carbohydrate-active enzymes, denitrification, and iron oxidation were enriched at the pre-bloom stage; genes involved in reductive citrate cycle for carbon fixation, carbon degradation, nitrification and phosphate transport were enhanced at the peak stage; and relative gene abundance related to sulfur oxidation, vitamin synthesis, and iron transport and storage was increased at the post-bloom stage. The ecological linkage analysis has shown that microbial functional potential especially the C/P/Fe metabolism were significantly linked to the fate of the algal blooms. Taken together, our results demonstrated that microorganisms displayed successional patterns not only at the community level, but also in the metabolic potential on HAB's progression. This work contributes to a growing understanding of microbial structural elasticity and functional plasticity and shed light on the potential mechanisms of microbial-mediated HAB trajectory.

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.

Conservation and Variability of the AUG Initiation Codon Context in Eukaryotes

Selection of the translation initiation site (TIS) is a crucial step during translation. In the 1980s Marylin Kozak performed key studies on vertebrate mRNAs to characterize the optimal TIS consensus sequence, the Kozak motif. Within this motif, conservation of nucleotides in crucial positions, namely a purine at -3 and a G at +4 (where the A of the AUG is numbered +1), is essential for TIS recognition. Ever since its characterization the Kozak motif has been regarded as the optimal sequence to initiate translation in all eukaryotes. We revisit here published in silico data on TIS consensus sequences, as well as experimental studies from diverse eukaryotic lineages, and propose that, while the -3A/G position is universally conserved, the remaining variability of the consensus sequences enables their classification as optimal, strong, and moderate TIS sequences.

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.

Transcriptomic Responses to Thermal Stress and Varied Phosphorus Conditions in Fugacium kawagutii

Coral reef-associated Symbiodiniaceae live in tropical and oligotrophic environments and are prone to heat and nutrient stress. How their metabolic pathways respond to pulses of warming and phosphorus (P) depletion is underexplored. Here, we conducted RNA-seq analysis to investigate transcriptomic responses to thermal stress, phosphate deprivation, and organic phosphorus (OP) replacement in Fugacium kawagutii. Using dual-algorithm (edgeR and NOIseq) to remedy the problem of no replicates, we conservatively found 357 differentially expressed genes (DEGs) under heat stress, potentially regulating cell wall modulation and the transport of iron, oxygen, and major nutrients. About 396 DEGs were detected under P deprivation and 671 under OP utilization, both mostly up-regulated and potentially involved in photosystem and defensome, despite different KEGG pathway enrichments. Additionally, we identified 221 genes that showed relatively stable expression levels across all conditions (likely core genes), mostly catalytic and binding proteins. This study reveals a wide range of, and in many cases previously unrecognized, molecular mechanisms in F. kawagutii to cope with heat stress and phosphorus-deficiency stress. Their quantitative expression dynamics, however, requires further verification with triplicated experiments, and the data reported here only provide clues for generating testable hypotheses about molecular mechanisms underpinning responses and adaptation in F. kawagutii to temperature and nutrient stresses.

The genetic intractability of Symbiodinium microadriaticum to standard algal transformation methods

Modern transformation and genome editing techniques have shown great success across a broad variety of organisms. However, no study of successfully applied genome editing has been reported in a dinoflagellate despite the first genetic transformation of Symbiodinium being published about 20 years ago. Using an array of different available transformation techniques, we attempted to transform Symbiodinium microadriaticum (CCMP2467), a dinoflagellate symbiont of reef-building corals, with the view to performing subsequent CRISPR-Cas9 mediated genome editing. Plasmid vectors designed for nuclear transformation containing the chloramphenicol resistance gene under the control of the CaMV p35S promoter as well as several putative endogenous promoters were used to test a variety of transformation techniques including biolistics, electroporation and agitation with silicon carbide whiskers. Chloroplast-targeted transformation was attempted using an engineered Symbiodinium chloroplast minicircle encoding a modified PsbA protein expected to confer atrazine resistance. We report that we have been unable to confer chloramphenicol or atrazine resistance on Symbiodinium microadriaticum strain CCMP2467.

Evidence for mitigation of coral bleaching by manganese

Unprecedented mass coral bleaching events due to global warming and overall seawater pollution have been observed worldwide over the last decades. Although metals are often considered as toxic substances for corals, some are essential at nanomolar concentrations for physiological processes such as photosynthesis and antioxidant defenses. This study was designed to elucidate, the individual and combined effects of nanomolar seawater enrichment in manganese (Mn) and iron (Fe), on the main physiological traits of Stylophora pistillata, maintained under normal growth and thermal stress conditions. We provide, for the first time, evidence that Mn is a key trace element for coral symbionts, enhancing cellular chlorophyll concentrations, photosynthetic efficiency and gross photosynthetic rates at ambient temperature. Our experiment also highlights the key role of Mn in increasing coral resistance to heat stress-induced bleaching. While Mn-enriched corals did not bleach and did not reduce their rates of photosynthesis and calcification, control corals experienced significant bleaching. On the contrary to Mn, Fe enrichment not only impaired calcification but induced significant bleaching. Such information is an important step towards a better understanding of the response of corals to seawater enrichment in metals. It can also explain, to some extent, species susceptibility to environmental stress.

Coral bleaching is linked to the capacity of the animal host to supply essential metals to the symbionts

Massive coral bleaching events result in extensive coral loss throughout the world. These events are mainly caused by seawater warming, but are exacerbated by the subsequent decrease in nutrient availability in surface waters. It has therefore been shown that nitrogen, phosphorus or iron limitation contribute to the underlying conditions by which thermal stress induces coral bleaching. Generally, information on the trophic ecology of trace elements (micronutrients) in corals, and on how they modulate the coral response to thermal stress is lacking. Here, we demonstrate for the first time that heterotrophic feeding (i.e. the capture of zooplankton prey by the coral host) and thermal stress induce significant changes in micro element concentrations and isotopic signatures of the scleractinian coral Stylophora pistillata. The results obtained first reveal that coral symbionts are the major sink for the heterotrophically acquired micronutrients and accumulate manganese, magnesium and iron from the food. These metals are involved in photosynthesis and antioxidant protection. In addition, we show that fed corals can maintain high micronutrient concentrations in the host tissue during thermal stress and do not bleach, whereas unfed corals experience a significant decrease in copper, zinc, boron, calcium and magnesium in the host tissue and bleach. In addition, the significant increase in δ⁶⁵ Cu and δ⁶⁶ Zn signature of symbionts and host tissue at high temperature suggests that these isotopic compositions are good proxy for stress in corals. Overall, present findings highlight a new way in which coral heterotrophy and micronutrient availability contribute to coral resistance to global warming and bleaching.

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.

Engineering Strategies to Decode and Enhance the Genomes of Coral Symbionts

Elevated sea surface temperatures from a severe and prolonged El Niño event (2014–2016) fueled by climate change have resulted in mass coral bleaching (loss of dinoflagellate photosymbionts, Symbiodinium spp., from coral tissues) and subsequent coral mortality, devastating reefs worldwide. Genetic variation within and between Symbiodinium species strongly influences the bleaching tolerance of corals, thus recent papers have called for genetic engineering of Symbiodinium to elucidate the genetic basis of bleaching-relevant Symbiodinium traits. However, while Symbiodinium has been intensively studied for over 50 years, genetic transformation of Symbiodinium has seen little success likely due to the large evolutionary divergence between Symbiodinium and other model eukaryotes rendering standard transformation systems incompatible. Here, we integrate the growing wealth of Symbiodinium next-generation sequencing data to design tailored genetic engineering strategies. Specifically, we develop a testable expression construct model that incorporates endogenous Symbiodinium promoters, terminators, and genes of interest, as well as an internal ribosomal entry site from a Symbiodinium virus. Furthermore, we assess the potential for CRISPR/Cas9 genome editing through new analyses of the three currently available Symbiodinium genomes. Finally, we discuss how genetic engineering could be applied to enhance the stress tolerance of Symbiodinium, and in turn, coral reefs.

Spliced leader-based analyses reveal the effects of polycyclic aromatic hydrocarbons on gene expression in the copepod Pseudodiaptomus poplesia

Polycyclic aromatic hydrocarbons (PAHs) are a group of toxic and carcinogenic pollutants that can adversely affect the development, growth and reproduction of marine organisms including copepods. However, knowledge on the molecular mechanisms regulating the response to PAH exposure in marine planktonic copepods is limited. In this study, we investigated the survival and gene expression of the calanoid copepod Pseudodiaptomus poplesia upon exposure to two PAHs, 1, 2-dimethylnaphthalene (1, 2-NAPH) and pyrene. Acute toxicity responses resulted in 96-h LC₅₀ of 788.98μgL^-1 and 54.68μgL^-1 for 1, 2-NAPH and pyrene, respectively. Using the recently discovered copepod spliced leader as a primer, we constructed full-length cDNA libraries from copepods exposed to sublethal concentrations and revealed 289 unique genes of diverse functions, including stress response genes and novel genes previously undocumented for this species. Eighty-three gene families were specifically expressed in PAH exposure libraries. We further analyzed the expression of seven target genes by reverse transcription-quantitative PCR in a time-course test with three sublethal concentrations. These target genes have primary roles in detoxification, oxidative defense, and signal transduction, and include different forms of glutathione S-transferase (GST), glutathione peroxidases (GPX), peroxiredoxin (PRDX), methylmalonate-semialdehyde dehydrogenase (MSDH) and ras-related C3 botulinum toxin substrate (RAC1). Expression stability of seven candidate reference genes were evaluated and the two most stable ones (RPL15 and RPS20 for 1, 2-NAPH exposure, RPL15 and EF1D for pyrene exposure) were used to normalize the expression levels of the target genes. Significant upregulation was detected in GST-T, GST-DE, GPX4, PRDX6 and RAC1 upon 1, 2-NAPH exposure, and GST-DE and MSDH upon pyrene exposure. These results indicated that the oxidative stress was induced and that signal transduction might be affected by PAH exposure in P. poplesia. However, gene upregulation was followed by a reduction in expression level towards 96h, indicating a threshold value of exposure time that leads to depressed gene expression. Prolonged exposure may cause dysfunction of detoxification and antioxidant machinery in P. poplesia. The transcriptional responses of GST-T, GPX2 and GPX4 upon pyrene exposure were minimal. Our results reveal the different sensitivity of P. poplesia to two PAHs at both the individual and transcriptional levels. As the first attempt, this study proved that copepod spliced leader is useful for obtaining full-length cDNA in P. poplesia exposed to PAHs and provided a valuable gene resource for this non-model species. This approach can be applied to other calanoid copepods exposed to various stressors, particularly under field conditions.

Molecular insights into a dinoflagellate bloom

In coastal waters worldwide, an increase in frequency and intensity of algal blooms has been attributed to eutrophication, with further increases predicted because of climate change. Yet, the cellular-level changes that occur in blooming algae remain largely unknown. Comparative metatranscriptomics was used to investigate the underlying molecular mechanisms associated with a dinoflagellate bloom in a eutrophied estuary. Here we show that under bloom conditions, there is increased expression of metabolic pathways indicative of rapidly growing cells, including energy production, carbon metabolism, transporters and synthesis of cellular membrane components. In addition, there is a prominence of highly expressed genes involved in the synthesis of membrane-associated molecules, including those for the production of glycosaminoglycans (GAGs), which may serve roles in nutrient acquisition and/or cell surface adhesion. Biotin and thiamine synthesis genes also increased expression along with several cobalamin biosynthesis-associated genes, suggesting processing of B₁₂ intermediates by dinoflagellates. The patterns in gene expression observed are consistent with bloom-forming dinoflagellates eliciting a cellular response to elevated nutrient demands and to promote interactions with their surrounding bacterial consortia, possibly in an effort to cultivate for enhancement of vitamin and nutrient exchanges and/or direct consumption. Our findings provide potential molecular targets for bloom characterization and management efforts.

Evidence for a role of viruses in the thermal sensitivity of coral photosymbionts

Symbiodinium, the dinoflagellate photosymbiont of corals, is posited to become more susceptible to viral infections when heat-stressed. To investigate this hypothesis, we mined transcriptome data of a thermosensitive and a thermotolerant type C1 Symbiodinium population at ambient (27 °C) and elevated (32°C) temperatures. We uncovered hundreds of transcripts from nucleocytoplasmic large double-stranded DNA viruses (NCLDVs) and the genome of a novel positive-sense single-stranded RNA virus (+ssRNAV). In the transcriptome of the thermosensitive population only, +ssRNAV transcripts had remarkable expression levels in the top 0.03% of all transcripts at 27 °C, but at 32 °C, expression levels of +ssRNAV transcripts decreased, while expression levels of anti-viral transcripts increased. In both transcriptomes, expression of NCLDV transcripts increased at 32 °C, but thermal induction of NCLDV transcripts involved in DNA manipulation was restricted to the thermosensitive population. Our findings reveal that viruses infecting Symbiodinium are affected by heat stress and may contribute to Symbiodinium thermal sensitivity.

Gene markers of fracture healing in early stage and the regulatory mechanism during the process using microarray analysis

BACKGROUND

The aim of this study was to explore crucial markers and uncover the regulatory mechanisms of fracture healing in the early stage.

METHODS

Gene expression profile of GSE45156 was downloaded, in which 3 fractured samples and 3 unfractured samples were used in our present study. Based on the threshold value, differentially expressed genes (DEGs) were selected between two kinds of samples using limma package in R. Enrichment analysis of these DEGs was performed by DAVID software. Furthermore, protein-protein interaction (PPI) network was established integrating information in STRING database, and visualized by Cytoscape software.

RESULTS

We identified a set of 960 DEGs including 509 up-regulated and 451 downregulated genes. Biological processes involving RNA splicing and cell cycle were significantly enriched for the up-regulated genes such as Snrpd2, Eftud2, Plk1 and Bub1b, whereas skeletal system development and bone development processes were predominant for down-regulated genes like Ubc. In the constructed PPI network, all the five genes were the predominant nodes, of which Snrpd2 was linked to Eftud2, while Bub1b was to interact with Plk1.

CONCLUSION

Five candidate genes crucial for indicating the process of fracture in early stage were identified. Eftud2, Snrpd2, Bub1b and Plk1 might function through the involvement of cell-cycle-related BP, while Ubc might influence the protein degradation during bone development. However, more experimental validations are needed to confirm these results.

Effects of Trace Metal Concentrations on the Growth of the Coral Endosymbiont Symbiodinium kawagutii

Symbiodinium is an indispensable endosymbiont in corals and the most important primary producer in coral reef ecosystems. During the past decades, coral bleaching attributed to the disruption of the symbiosis has frequently occurred resulting in reduction of coral reef coverage globally. Growth and proliferation of corals require some specific trace metals that are essential components of pertinent biochemical processes, such as in photosynthetic systems and electron transport chains. In addition, trace metals are vital in the survival of corals against oxidative stress because these metals serve as enzymatic cofactors in antioxidative defense mechanisms. The basic knowledge about trace metal requirements of Symbiodinium is lacking. Here we show that the requirement of Symbiodinium kawagutii for antioxidant-associated trace metals exhibits the following order: Fe >> Cu/Zn/Mn >> Ni. In growth media with Cu, Zn, Mn, and varying Fe concentrations, we observed that Cu, Zn, and Mn cellular quotas were inversely related to Fe concentrations. In the absence of Cu, Zn, and Mn, growth rates increased with increasing inorganic Fe concentrations up to 1250 pM, indicating the relatively high Fe requirement for Symbiodinium growth and potential functional complementarity of these metals. These results demonstrate the relative importance of trace metals to sustain Symbiodinium growth and a potential metal inter replacement strategy in Symbiodinium to ensure survival of coral reefs in an oligotrophic and stressful environment.

The Symbiodinium kawagutii genome illuminates dinoflagellate gene expression and coral symbiosis

Dinoflagellates are important components of marine ecosystems and essential coral symbionts, yet little is known about their genomes. We report here on the analysis of a high-quality assembly from the 1180-megabase genome of Symbiodinium kawagutii. We annotated protein-coding genes and identified Symbiodinium-specific gene families. No whole-genome duplication was observed, but instead we found active (retro)transposition and gene family expansion, especially in processes important for successful symbiosis with corals. We also documented genes potentially governing sexual reproduction and cyst formation, novel promoter elements, and a microRNA system potentially regulating gene expression in both symbiont and coral. We found biochemical complementarity between genomes of S. kawagutii and the anthozoan Acropora, indicative of host-symbiont coevolution, providing a resource for studying the molecular basis and evolution of coral symbiosis.

Spliced leader RNA trans-splicing discovered in copepods

Copepods are one of the most abundant metazoans in the marine ecosystem, constituting a critical link in aquatic food webs and contributing significantly to the global carbon budget, yet molecular mechanisms of their gene expression are not well understood. Here we report the detection of spliced leader (SL) trans-splicing in calanoid copepods. We have examined nine species of wild-caught copepods from Jiaozhou Bay, China that represent the major families of the calanoids. All these species contained a common 46-nt SL (CopepodSL). We further determined the size of CopepodSL precursor RNA (slRNA; 108-158 nt) through genomic analysis and 3′-RACE technique, which was confirmed by RNA blot analysis. Structure modeling showed that the copepod slRNA folded into typical slRNA secondary structures. Using a CopepodSL-based primer set, we selectively enriched and sequenced copepod full-length cDNAs, which led to the characterization of copepod transcripts and the cataloging of the complete set of 79 eukaryotic cytoplasmic ribosomal proteins (cRPs) for a single copepod species. We uncovered the SL trans-splicing in copepod natural populations, and demonstrated that CopepodSL was a sensitive and specific tool for copepod transcriptomic studies at both the individual and population levels and that it would be useful for metatranscriptomic analysis of copepods.

The diversity of algal phospholipase D homologs revealed by biocomputational analysis

Phospholipase D (PLD) participates in the formation of phosphatidic acid, a precursor in glycerolipid biosynthesis and a second messenger. PLDs are part of a superfamily of proteins that hydrolyze phosphodiesters and share a catalytic motif, HxKxxxxD, and hence a mechanism of action. Although HKD-PLDs have been thoroughly characterized in plants, animals and bacteria, very little is known about these enzymes in algae. To fill this gap in knowledge, we performed a biocomputational analysis by means of HMMER iterative profiling, using most eukaryotic algae genomes available. Phylogenetic analysis revealed that algae exhibit very few eukaryotic-type PLDs but possess, instead, many bacteria-like PLDs. Among algae eukaryotic-type PLDs, we identified C2-PLDs and PXPH-like PLDs. In addition, the dinoflagellate Alexandrium tamarense features several proteins phylogenetically related to oomycete PLDs. Our phylogenetic analysis also showed that algae bacteria-like PLDs (proteins with putative PLD activity) fall into five clades, three of which are novel lineages in eukaryotes, composed almost entirely of algae. Specifically, Clade II is almost exclusive to diatoms, whereas Clade I and IV are mainly represented by proteins from prasinophytes. The other two clades are composed of mitochondrial PLDs (Clade V or Mito-PLDs), previously found in mammals, and a subfamily of potentially secreted proteins (Clade III or SP-PLDs), which includes a homolog formerly characterized in rice. In addition, our phylogenetic analysis shows that algae have non-PLD members within the bacteria-like HKD superfamily with putative cardiolipin synthase and phosphatidylserine/phosphatidylglycerophosphate synthase activities. Altogether, our results show that eukaryotic algae possess a moderate number of PLDs that belong to very diverse phylogenetic groups.

Single-cell transcriptomics using spliced leader PCR: Evidence for multiple losses of photosynthesis in polykrikoid dinoflagellates

BACKGROUND

Most microbial eukaryotes are uncultivated and thus poorly suited to standard genomic techniques. This is the case for Polykrikos lebouriae, a dinoflagellate with ultrastructurally aberrant plastids. It has been suggested that these plastids stem from a novel symbiosis with either a diatom or haptophyte, but this hypothesis has been difficult to test as P. lebouriae dwells in marine sand rife with potential genetic contaminants.

RESULTS

We applied spliced-leader targeted PCR (SLPCR) to obtain dinoflagellate-specific transcriptomes on single-cell isolates of P. lebouriae from marine sediments. Polykrikos lebouriae expressed nuclear-encoded photosynthetic genes that were characteristic of the peridinin-plastids of dinoflagellates, rather than those from a diatom of haptophyte. We confirmed these findings at the genomic level using multiple displacement amplification (MDA) to obtain a partial plastome of P. lebouriae.

CONCLUSION

From these data, we infer that P. lebouriae has retained the peridinin plastids ancestral for dinoflagellates as a whole, while its closest relatives have lost photosynthesis multiple times independently. We discuss these losses with reference to mixotrophy in polykrikoid dinoflagellates. Our findings demonstrate new levels of variation associated with the peridinin plastids of dinoflagellates and the usefulness of SLPCR approaches on single cell isolates. Unlike other transcriptomic methods, SLPCR has taxonomic specificity, and can in principle be adapted to different splice-leader bearing groups.

Symbiodinium transcriptome and global responses of cells to immediate changes in light intensity when grown under autotrophic or mixotrophic conditions

Symbiosis between unicellular dinoflagellates (genus Symbiodinium) and their cnidarian hosts (e.g. corals, sea anemones) is the foundation of coral reef ecosystems. Dysfunction of this symbiosis under changing environmental conditions has led to global reef decline. Little information is known about Symbiodinium gene expression and mechanisms by which light impacts host-symbiont associations. To address these issues, we generated a transcriptome from axenic Symbiodinium strain SSB01. Here we report features of the transcriptome, including occurrence and length distribution of spliced leader sequences, the functional landscape of encoded proteins and the impact of light on gene expression. Expression of many Symbiodinium genes appears to be significantly impacted by light. Transcript encoding cryptochrome 2 declined in high light while some transcripts for Regulators of Chromatin Condensation (RCC1) declined in the dark. We also identified a transcript encoding a light harvesting AcpPC protein with homology to Chlamydomonas LHCSR2. The level of this transcript increased in high light autotrophic conditions, suggesting that it is involved in photo-protection and the dissipation of excess absorbed light energy. The most extensive changes in transcript abundances occurred when the algae were transferred from low light to darkness. Interestingly, transcripts encoding several cell adhesion proteins rapidly declined following movement of cultures to the dark, which correlated with a dramatic change in cell surface morphology, likely reflecting the complexity of the extracellular matrix. Thus, light-sensitive cell adhesion proteins may play a role in establishing surface architecture, which may in turn alter interactions between the endosymbiont and its host.

Transcriptomic characterization of the enzymatic antioxidants FeSOD, MnSOD, APX and KatG in the dinoflagellate genus Symbiodinium

Background

The diversity of the symbiotic dinoflagellate Symbiodinium sp., as assessed by genetic markers, is well established. To what extent this diversity is reflected on the amino acid level of functional genes such as enzymatic antioxidants that play an important role in thermal stress tolerance of the coral-Symbiodinium symbiosis is, however, unknown. Here we present a predicted structural analysis and phylogenetic characterization of the enzymatic antioxidant repertoire of the genus Symbiodinium. We also report gene expression and enzymatic activity under short-term thermal stress in Symbiodinium of the B1 genotype.

Results

Based on eight different ITS2 types, covering six clades, multiple protein isoforms for three of the four investigated antioxidants (ascorbate peroxidase [APX], catalase peroxidase [KatG], manganese superoxide dismutase [MnSOD]) are present in the genus Symbiodinium. Amino acid sequences of both SOD metalloforms (Fe/Mn), as well as KatG, exhibited a number of prokaryotic characteristics that were also supported by the protein phylogeny. In contrast to the bacterial form, KatG in Symbiodinium is characterized by extended functionally important loops and a shortened C-terminal domain. Intercladal sequence variations were found to be much higher in both peroxidases, compared to SODs. For APX, these variable residues involve binding sites for substrates and cofactors, and might therefore differentially affect the catalytic properties of this enzyme between clades. While expression of antioxidant genes was successfully measured in Symbiodinium B1, it was not possible to assess the link between gene expression and protein activity due to high variability in expression between replicates, and little response in their enzymatic activity over the three-day experimental period.

Conclusions

The genus Symbiodinium has a diverse enzymatic antioxidant repertoire that has similarities to prokaryotes, potentially as a result of horizontal gene transfer or events of secondary endosymbiosis. Different degrees of sequence evolution between SODs and peroxidases might be the result of potential selective pressure on the conserved molecular function of SODs as the first line of defence. In contrast, genetic redundancy of hydrogen peroxide scavenging enzymes might permit the observed variations in peroxidase sequences. Our data and successful measurement of antioxidant gene expression in Symbiodinium will serve as basis for further studies of coral health.

Electronic supplementary material

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

Light-promoted rhodopsin expression and starvation survival in the marine dinoflagellate Oxyrrhis marina

The discovery of microbial rhodopsins in marine proteobacteria changed the dogma that photosynthesis is the only pathway to use the solar energy for biological utilization in the marine environment. Although homologs of these rhodopsins have been identified in dinoflagellates, the diversity of the encoding genes and their physiological roles remain unexplored. As an initial step toward addressing the gap, we conducted high-throughput transcriptome sequencing on Oxyrrhis marina to retrieve rhodopsin transcripts, rapid amplification of cDNA ends to isolate full-length cDNAs of dominant representatives, and quantitative reverse-transcription PCR to investigate their expression under varying conditions. Our phylogenetic analyses showed that O. marina contained both the proton-pumping type (PR) and sensory type (SR) rhodopsins, and the transcriptome data showed that the PR type dominated over the SR type. We compared rhodopsin gene expression for cultures kept under light: dark cycle and continuous darkness in a time course of 24 days without feeding. Although both types of rhodopsin were expressed under the two conditions, the expression levels of PR were much higher than SR, consistent with the transcriptomic data. Furthermore, relative to cultures kept in the dark, rhodopsin expression levels and cell survival rate were both higher in cultures grown in the light. This is the first report of light-dependent promotion of starvation survival and concomitant promotion of PR expression in a eukaryote. While direct evidence needs to come from functional test on rhodopsins in vitro or gene knockout/knockdown experiments, our results suggest that the proton-pumping rhodopsin might be responsible for the light-enhanced survival of O. marina, as previously demonstrated in bacteria.

A genomic approach to coral-dinoflagellate symbiosis: studies of Acropora digitifera and Symbiodinium minutum

Far more intimate knowledge of scleractinian coral biology is essential in order to understand how diverse coral-symbiont endosymbioses have been established. In particular, molecular and cellular mechanisms enabling the establishment and maintenance of obligate endosymbiosis with photosynthetic dinoflagellates require further clarification. By extension, such understanding may also shed light upon environmental conditions that promote the collapse of this mutualism. Genomic data undergird studies of all symbiotic processes. Here we review recent genomic data derived from the scleractinian coral, Acropora digitifera, and the endosymbiotic dinoflagellate, Symbiodinium minutum. We discuss Acropora genes involved in calcification, embryonic development, innate immunity, apoptosis, autophagy, UV resistance, fluorescence, photoreceptors, circadian clocks, etc. We also detail gene loss in amino acid metabolism that may explain at least part of the Acropora stress-response. Characteristic features of the Symbiodinium genome are also reviewed, focusing on the expansion of certain gene families, the molecular basis for permanently condensed chromatin, unique spliceosomal splicing, and unusual gene arrangement. Salient features of the Symbiodinium plastid and mitochondrial genomes are also illuminated. Although many questions regarding these interdependent genomes remain, we summarize information necessary for future studies of coral-dinoflagellate endosymbiosis.

Biology of the Marine Heterotrophic Dinoflagellate Oxyrrhis marina: Current Status and Future Directions

Heterotrophic dinoflagellates are prevalent protists in marine environments, which play an important role in the carbon cycling and energy flow in the marine planktonic community. Oxyrrhismarina (Dinophyceae), a widespread heterotrophic dinoflagellate, is a model species used for a broad range of ecological, biogeographic, and evolutionary studies. Despite the increasing research effort on this species, there lacks a synthesis of the existing data and a coherent picture of this organism. Here we reviewed the literature to provide an overview of what is known regarding the biology of O. marina, and identify areas where further studies are needed. As an early branch of the dinoflagellate lineage, O. marina shares similarity with typical dinoflagellates in permanent condensed chromosomes, less abundant nucleosome proteins compared to other eukaryotes, multiple gene copies, the occurrence of trans-splicing in nucleus-encoded mRNAs, highly fragmented mitochondrial genome, and disuse of ATG as a start codon for mitochondrial genes. On the other hand, O. marina also exhibits some distinct cytological features (e.g., different flagellar structure, absence of girdle and sulcus or pustules, use of intranuclear spindle in mitosis, presence of nuclear plaque, and absence of birefringent periodic banded chromosomal structure) and genetic features (e.g., a single histone-like DNA-associated protein, cob-cox3 gene fusion, 5' oligo-U cap in the mitochondrial transcripts of protein-coding genes, the absence of mRNA editing, the presence of stop codon in the fused cob-cox3 mRNA produced by post-transcriptional oligoadenylation, and vestigial plastid genes). The best-studied biology of this dinoflagellate is probably the prey and predators types, which include a wide range of organisms. On the other hand, the abundance of this species in the natural waters and its controlling factors, genome organization and gene expression regulation that underlie the unusual cytological and ecological characteristics are among the areas that urgently need study.

Signal recognition particle RNA in dinoflagellates and the Perkinsid Perkinsus marinus

In dinoflagellates and perkinsids, the molecular structure of the protein translocating machinery is unclear. Here, we identified several types of full-length signal recognition particle (SRP) RNA genes from Karenia brevis (dinoflagellate) and Perkinsus marinus (perkinsid). We also identified the four SRP S-domain proteins, but not the two Alu domain proteins, from P. marinus and several dinoflagellates. We mapped both ends of SRP RNA transcripts from K. brevis and P. marinus, and obtained the 3' end from four other dinoflagellates. The lengths of SRP RNA are predicted to be ∼260-300 nt in dinoflagellates and 280-285 nt in P. marinus. Although these SRP RNA sequences are substantially variable, the predicted structures are similar. The genomic organization of the SRP RNA gene differs among species. In K. brevis, this gene is located downstream of the spliced leader (SL) RNA, either as SL RNA-SRP RNA-tRNA gene tandem repeats, or within a SL RNA-SRP RNA-tRNA-U6-5S rRNA gene cluster. In other dinoflagellates, SRP RNA does not cluster with SL RNA or 5S rRNA genes. The majority of P. marinus SRP RNA genes array as tandem repeats without the above-mentioned small RNA genes. Our results capture a snapshot of a potentially complex evolutionary history of SRP RNA in alveolates.

The PsbO homolog from Symbiodinium kawagutii (Dinophyceae) characterized using biochemical and molecular methods

A photosystem II component, the PsbO protein is essential for maximum rates of oxygen production during photosynthesis, and has been extensively characterized in plants and cyanobacteria but not in symbiotic dinoflagellates. Its close interaction with D1 protein has important environmental implications since D1 has been identified as the primary site of damage in endosymbiotic dinoflagellates after thermal stress. We identified and biochemically characterized the PsbO homolog from Symbiodinium kawagutii as a 28-kDa protein, and immunolocalized it to chloroplast membranes. Chloroplast association was further confirmed by western blot on photosynthetic membrane preparations. TX-114 phase partitioning, chromatography, and SDS-PAGE for single band separation and partial peptide sequencing yielded peptides identical or with high identity to PsbO from dinoflagellates. Analysis of a cDNA library revealed three genes differing by only one aminoacid residue in the in silico-translated ORFs despite greater differences at nucleotide level in the untranslated, putative regulatory sequences. The consensus full amino acid sequence displayed all the characteristic domains and features of PsbO from other sources, but changes in functionally critical, highly conserved motifs were detected. Our biochemical, molecular, and immunolocalization data led to the conclusion that the 28-kDa protein from S. kawagutii is the PsbO homolog, thereby named SkPsbO. We discuss the implications of critical amino acid substitutions for a putative regulatory role of this protein.