Symbiodinium transcriptomes: genome insights into the dinoflagellate symbionts of reef-building corals

Trade-off between photosymbiosis and innate immunity influences cnidarian's response to pathogenic bacteria

Mutualistic relationships with photosynthetic organisms are common in cnidarians, which form an intracellular symbiosis with dinoflagellates in the family Symbiodiniaceae. The establishment and maintenance of these symbionts are associated with the suppression of key host immune factors. Because of this, there are potential trade-offs between the nutrition that cnidarian hosts gain from their symbionts and their ability to successfully defend themselves from pathogens. To investigate these potential trade-offs, we utilized the facultatively symbiotic polyps of the upside-down jellyfish Cassiopea xamachana and exposed aposymbiotic and symbiotic polyps to the pathogen Serratia marcescens. Symbiotic polyps had a lower probability of survival following S. marcescens exposure. Gene expression analyses 24 hours following pathogen exposure indicate that symbiotic animals mounted a more damaging immune response, with higher levels of inflammation and oxidative stress likely resulting in more severe disruptions to cellular homeostasis. Underlying this more damaging immune response may be differences in constitutive and pathogen-induced expression of immune transcription factors between aposymbiotic and symbiotic polyps rather than broadscale immune suppression during symbiosis. Our findings indicate that in facultatively symbiotic polyps, hosting symbionts limits C. xamachana's ability to survive pathogen exposure, indicating a trade-off between symbiosis and immunity that has potential implications for coral disease research.

Heat stress in symbiotic dinoflagellates: Implications on oxidative stress and cellular changes

Global warming has been shown to harmfully affect symbiosis between Symbiodiniaceae and other marine invertebrates. When symbiotic dinoflagellates (the genus Breviolum) were in vitro exposed to acute heat stress of +7 °C for a period of 5 days, the results revealed the negative impact on all physiological and other cellular parameters measured. Elevated temperatures resulted in a severe reduction in algal density of up to 9.5-fold, as well as pigment concentrations, indicating the status of the physiological stress and early signs of photo-bleaching. Reactive oxygen species (ROS) were increased in all heated dinoflagellate cells, while the antioxidant-reduced glutathione levels initially dropped on day one but increased under prolonged temperature stress. The cell viability parameters were reduced by 97 % over the heating period, with an increased proportion of apoptotic and necrotic cells. Autofluorescence (AF) for Cy5-PE 660-20 was reduced from 1.7-fold at day 1 to up to 50-fold drop at the end of heating time, indicating that the AF changes were highly sensitive to heat stress and that it could be an extremely sensitive tool for assessing the functionality of algal photosynthetic machinery. The addition of the drug 5-AZA-2'-deoxycytidine (5-AZA), which inhibits DNA methylation processes, was assessed in parallel and contributed to some alterations in algal cellular stress response. The presence of drug 5-AZA combined with the temperature stress had an additional impact on Symbiodiniaceae density and cell complexity, including the AF levels. These variations in cellular stress response under heat stress and compromised DNA methylation conditions may indicate the importance of this epigenetic mechanism for symbiotic dinoflagellate thermal tolerance adaptability over a longer period, which needs further exploration. Consequently, the increased ROS levels and changes in AF signals reported during ongoing heat stress in dinoflagellate cells could be used as early stress biomarkers in these microalgae and potentially other photosynthetic species.

Multifaceted Dinoflagellates and the Marine Model Prorocentrum cordatum

BACKGROUND

Dinoflagellates are a monophyletic group within the taxon Alveolata, which comprises unicellular eukaryotes. Dinoflagellates have long been studied for their organismic and morphologic diversity as well as striking cellular features. They have a main size range of 10-100 µm, a complex "cell covering", exceptionally large genomes (∼1-250 Gbp with a mean of 50,000 protein-encoding genes) spread over a variable number of highly condensed chromosomes, and perform a closed mitosis with extranuclear spindles (dinomitosis). Photosynthetic, marine, and free-living Prorocentrum cordatum is a ubiquitously occurring, bloom-forming dinoflagellate, and an emerging model system, particularly with respect to systems biology.

SUMMARY

Focused ion beam/scanning electron microscopy (FIB/SEM) analysis of P. cordatum recently revealed (i) a flattened nucleus with unusual structural features and a total of 62 tightly packed chromosomes, (ii) a single, barrel-shaped chloroplast devoid of grana and harboring multiple starch granules, (iii) a single, highly reticular mitochondrion, and (iv) multiple phosphate and lipid storage bodies. Comprehensive proteomics of subcellular fractions suggested (i) major basic nuclear proteins to participate in chromosome condensation, (ii) composition of nuclear pores to differ from standard knowledge, (iii) photosystems I and II, chloroplast complex I, and chlorophyll a-b binding light-harvesting complex to form a large megacomplex (>1.5 MDa), and (iv) an extraordinary richness in pigment-binding proteins. Systems biology-level investigation of heat stress response demonstrated a concerted down-regulation of CO2-concentrating mechanisms, CO2-fixation, central metabolism, and monomer biosynthesis, which agrees with reduced growth yields.

KEY MESSAGES

FIB/SEM analysis revealed new insights into the remarkable subcellular architecture of P. cordatum, complemented by proteogenomic unraveling of novel nuclear structures and a photosynthetic megacomplex. These recent findings are put in the wider context of current understanding of dinoflagellates.

Genomic data reveals habitat partitioning in massive Porites on Guam, Micronesia

Corals in marginal reef habitats generally exhibit less bleaching and associated mortality compared to nearby corals in more pristine reef environments. It is unclear, however, if these differences are due to environmental differences, including turbidity, or genomic differences between the coral hosts in these different environments. One particularly interesting case is in the coral genus Porites, which contains numerous morphologically similar massive Porites species inhabiting a wide range of reef habitats, from turbid river deltas and stagnant back reefs to high-energy fore reefs. Here, we generate ddRAD data for 172 Porites corals from river delta and adjacent (<0.5 km) fore reef populations on Guam to assess the extent of genetic differentiation among massive Porites corals in these two contrasting environments and throughout the island. Phylogenetic and population genomic analyses consistently identify seven different clades of massive Porites, with the two largest clades predominantly inhabiting either river deltas or fore reefs, respectively. No population structure was detected in the two largest clades, and Cladocopium was the dominant symbiont genus in all clades and environments. The perceived bleaching resilience of corals in marginal reefs may therefore be attributed to interspecific differences between morphologically similar species, in addition to potentially mediating environmental differences. Marginal reef environments may therefore not provide a suitable refuge for many reef corals in a heating world, but instead host additional cryptic coral diversity.

Transcriptomic signatures across a critical sedimentation threshold in a major reef-building coral

Sedimentation is a major cause of global near-shore coral reef decline. Although the negative impacts of sedimentation on coral reef community composition have been well-documented, the effects of sedimentation on coral metabolism in situ have received comparatively little attention. Using transcriptomics, we identified gene expression patterns changing across a previously defined sedimentation threshold that was deemed critical due to changes in coral cover and community composition. We identified genes, pathways, and molecular processes associated with this transition that may allow corals, such as Porites lobata, to tolerate chronic, severe sedimentation and persist in turbid environments. Alternative energy generation pathways may help P. lobata maintain a persistent stress response to survive when the availability of light and oxygen is diminished. We found evidence for the expression of genes linked to increased environmental sensing and cellular communication that likely allow P. lobata to efficiently respond to sedimentation stress and associated pathogen challenges. Cell damage increases under stress; consequently, we found apoptosis pathways over-represented under severe sedimentation, a likely consequence of damaged cell removal to maintain colony integrity. The results presented here provide a framework for the response of P. lobata to sedimentation stress under field conditions. Testing this framework and its related hypotheses using multi-omics approaches can deepen our understanding of the metabolic plasticity and acclimation potential of corals to sedimentation and their resilience in turbid reef systems.

Multiomics data integration, limitations, and prospects to reveal the metabolic activity of the coral holobiont

Since their radiation in the Middle Triassic period ∼240 million years ago, stony corals have survived past climate fluctuations and five mass extinctions. Their long-term survival underscores the inherent resilience of corals, particularly when considering the nutrient-poor marine environments in which they have thrived. However, coral bleaching has emerged as a global threat to coral survival, requiring rapid advancements in coral research to understand holobiont stress responses and allow for interventions before extensive bleaching occurs. This review encompasses the potential, as well as the limits, of multiomics data applications when applied to the coral holobiont. Synopses for how different omics tools have been applied to date and their current restrictions are discussed, in addition to ways these restrictions may be overcome, such as recruiting new technology to studies, utilizing novel bioinformatics approaches, and generally integrating omics data. Lastly, this review presents considerations for the design of holobiont multiomics studies to support lab-to-field advancements of coral stress marker monitoring systems. Although much of the bleaching mechanism has eluded investigation to date, multiomic studies have already produced key findings regarding the holobiont's stress response, and have the potential to advance the field further.

Functional analyses of bacterial genomes found in Symbiodiniaceae genome assemblies

Bacterial-algal interactions strongly influence marine ecosystems. Bacterial communities in cultured dinoflagellates of the family Symbiodiniaceae have been characterized by metagenomics. However, little is known about whole-genome analysis of marine bacteria associated with these dinoflagellates. We performed in silico analysis of four bacterial genomes from cultures of four dinoflagellates of the genera Symbiodinium, Breviolum, Cladocopium and Durusdinium. Comparative analysis showed that the former three contain the alphaproteobacterial family Parvibaculaceae and that the Durusdinium culture includes the family Sphingomonadaceae. There were no large genomic reductions in the alphaproteobacteria with genome sizes of 2.9-3.9 Mb, implying they are not obligate intracellular bacteria. Genomic annotations of three Parvibaculaceae detected the gene for diacetylchitobiose deacetylase (Dac), which may be involved in the degradation of dinoflagellate cell surfaces. They also had metabolic genes for dissimilatory nitrate reduction to ammonium (DNRA) in the nitrogen (N) cycle and cobalamin (vitamin B12 ) biosynthetic genes in the salvage pathway. Those three characters were not found in the Sphingomonadaceae genome. Predicted biosynthetic gene clusters for secondary metabolites indicated that the Parvibaculaceae likely produce the same secondary metabolites. Our study suggests that the Parvibaculaceae is a major resident of Symbiodiniaceae cultures with antibiotics.

Symbiosis modulates gene expression of symbionts, but not coral hosts, under thermal challenge

Increasing ocean temperatures are causing dysbiosis between coral hosts and their symbionts. Previous work suggests that coral host gene expression responds more strongly to environmental stress compared to their intracellular symbionts; however, the causes and consequences of this phenomenon remain untested. We hypothesized that symbionts are less responsive because hosts modulate symbiont environments to buffer stress. To test this hypothesis, we leveraged the facultative symbiosis between the scleractinian coral Oculina arbuscula and its symbiont Breviolum psygmophilum to characterize gene expression responses of both symbiotic partners in and ex hospite under thermal challenges. To characterize host and in hospite symbiont responses, symbiotic and aposymbiotic O. arbuscula were exposed to three treatments: (1) control (18°C), (2) heat (32°C), and (3) cold (6°C). This experiment was replicated with B. psygmophilum cultured from O. arbuscula to characterize ex hospite symbiont responses. Both thermal challenges elicited classic environmental stress responses (ESRs) in O. arbuscula regardless of symbiotic state, with hosts responding more strongly to cold challenge. Hosts also exhibited stronger responses than in hospite symbionts. In and ex hospite B. psygmophilum both down-regulated gene ontology pathways associated with photosynthesis under thermal challenge; however, ex hospite symbionts exhibited greater gene expression plasticity and differential expression of genes associated with ESRs. Taken together, these findings suggest that O. arbuscula hosts may buffer environments of B. psygmophilum symbionts; however, we outline the future work needed to confirm this hypothesis.

Polyphyletic origin of saxitoxin biosynthesis genes in the marine dinoflagellate Alexandrium revealed by comparative transcriptomics

The marine dinoflagellate Alexandrium is known to form harmful algal blooms, and at least 14 species within the genus can produce saxitoxins (STXs). STX biosynthesis genes (sxt) are individually revealed in toxic dinoflagellates; however, the evolutionary history remains controversial. Herein, we determined the transcriptome sequences of toxic Alexandrium (A. catenella and A. pacificum) and non-toxic Alexandrium (A. fraterculus and A. fragae) and characterized their sxt by focusing on evolutionary events and STX production. Comparative transcriptome analysis revealed higher homology of the sxt in toxic Alexandrium than in non-toxic species. Notably, non-toxic Alexandrium spp. were found to have lost two sxt core genes, namely sxtA4 and sxtG. Expression levels of 28 transcripts related to eight sxt core genes showed that sxtA, sxtG, and sxtI were relatively high (>1.5) in the toxic group compared to the non-toxic group. In contrast, the non-toxic group showed high expression levels in sxtU (1.9) and sxtD (1.7). Phylogenetic tree comparisons revealed distinct evolutionary patterns between 28S rDNA and sxtA, sxtB, sxtI, sxtD, and sxtU. However, similar topology was observed between 28S rDNA, sxtS, and sxtH/T. In the sxtB and sxtI phylogeny trees, toxic Alexandrium and cyanobacteria were clustered together, separating from non-toxic species. These suggest that Alexandrium may acquire sxt genes independently via horizontal gene transfer from toxic cyanobacteria and other multiple sources, demonstrating monocistronic transcripts of sxt in dinoflagellates.

Intersection of coral molecular responses to a localized mortality event and ex situ deoxygenation

In July 2016, East Bank of Flower Garden Banks (FGB) National Marine Sanctuary experienced a localized mortality event (LME) of multiple invertebrate species that ultimately led to reductions in coral cover. Abiotic data taken directly after the event suggested that acute deoxygenation contributed to the mortality. Despite the large impact of this event on the coral community, there was no direct evidence that this LME was driven by acute deoxygenation, and thus we explored whether gene expression responses of corals to the LME would indicate what abiotic factors may have contributed to the LME. Gene expression of affected and unaffected corals sampled during the mortality event revealed evidence of the physiological consequences of the LME on coral hosts and their algal symbionts from two congeneric species (Orbicella franksi and Orbicella faveolata). Affected colonies of both species differentially regulated genes involved in mitochondrial regulation and oxidative stress. To further test the hypothesis that deoxygenation led to the LME, we measured coral host and algal symbiont gene expression in response to ex situ experimental deoxygenation (control = 6.9 ± 0.08 mg L-1, anoxic = 0.083 ± 0.017 mg L-1) in healthy O. faveolata colonies from the FGB. However, this deoxygenation experiment revealed divergent gene expression patterns compared to the corals sampled during the LME and was more similar to a generalized coral environmental stress response. It is therefore likely that while the LME was connected to low oxygen, it was a series of interconnected stressors that elicited the unique gene expression responses observed here. These in situ and ex situ data highlight how field responses to stressors are unique from those in controlled laboratory conditions, and that the complexities of deoxygenation events in the field likely arise from interactions between multiple environmental factors simultaneously.

Synergistic response to climate stressors in coral is associated with genotypic variation in baseline expression

As environments are rapidly reshaped due to climate change, phenotypic plasticity plays an important role in the ability of organisms to persist and is considered an especially important acclimatization mechanism for long-lived sessile organisms such as reef-building corals. Often, this ability of a single genotype to display multiple phenotypes depending on the environment is modulated by changes in gene expression, which can vary in response to environmental changes via two mechanisms: baseline expression and expression plasticity. We used transcriptome-wide expression profiling of eleven genotypes of common-gardened Acropora cervicornis to explore genotypic variation in the expression response to thermal and acidification stress, both individually and in combination. We show that the combination of these two stressors elicits a synergistic gene expression response, and that both baseline expression and expression plasticity in response to stress show genotypic variation. Additionally, we demonstrate that frontloading of a large module of coexpressed genes is associated with greater retention of algal symbionts under combined stress. These results illustrate that variation in the gene expression response of individuals to climate change stressors can persist even when individuals have shared environmental histories, affecting their performance under future climate change scenarios.

Unique high-temperature tolerance mechanisms of zoochlorellae Symbiochlorum hainanensis derived from scleractinian coral Porites lutea

Global warming is a key issue that causes coral bleaching mainly because of the thermosensitivity of zooxanthellae. Compared with the well-studied zooxanthellae Symbiodiniaceae in coral holobionts, we rarely know about other coral symbiotic algae, let alone their thermal tolerance. In this study, a zoochlorellae, Symbiochlorum hainanensis, isolated from the coral Porites lutea, was proven to have a threshold temperature of 38°C. Meanwhile, unique high-temperature tolerance mechanisms were suggested by integrated transcriptomics and real-time quantitative PCR, physiological and biochemical analyses, and electron microscopy observation. Under heat stress, S. hainanensis shared some similar response strategies with zooxanthellae Effrenium sp., such as increased ascorbate peroxidase, glutathione peroxidase, superoxide dismutase activities and chlorophyll a, thiamine, and thiamine phosphate contents. In particular, more chloroplast internal layered structure, increased CAT activity, enhanced selenate reduction, and thylakoid assembly pathways were highlighted for S. hainanensis's high-temperature tolerance. Notably, it is the first time to reveal a whole selenate reduction pathway from SeO42- to Se2- and its contribution to the high-temperature tolerance of S. hainanensis. These unique mechanisms, including antioxidation and maintaining photosynthesis homeostasis, efficiently ensure the high-temperature tolerance of S. hainanensis than Effrenium sp. Compared with the thermosensitivity of coral symbiotic zooxanthellae Symbiodiniaceae, this study provides novel insights into the high-temperature tolerance mechanisms of coral symbiotic zoochlorellae S. hainanensis, which will contribute to corals' survival in the warming oceans caused by global climate change.

IMPORTANCE

The increasing ocean temperature above 31°C-32°C might trigger a breakdown of the coral-Symbiodiniaceae symbioses or coral bleaching because of the thermosensitivity of Symbiodiniaceae; therefore, the exploration of alternative coral symbiotic algae with high-temperature tolerance is important for the corals' protection under warming oceans. This study proves that zoochlorellae Symbiochlorum hainanensis can tolerate 38°C, which is the highest temperature tolerance known for coral symbiotic algae to date, with unique high-temperature tolerance mechanisms. Particularly, for the first time, an internal selenium antioxidant mechanism of coral symbiotic S. hainanensis to high temperature was suggested.

Two novel glutathione S-transferase (GST) genes in the toxic marine dinoflagellate Alexandrium pacificum and their transcriptional responses to environmental contaminants

The present study identified two novel glutathione S-transferase (GST) genes from the toxic dinoflagellate Alexandrium pacificum and examined their molecular characteristics and transcriptional responses to algicides and environmental contaminants. Bioinformatic analysis revealed that both ApGSTs are cytosolic, belonging to the chi-like class (ApGST1) and an undefined class (ApGST2). The overall expression of ApGSTs showed similar patterns depending on the exposed contaminants, while they were differently regulated by polychlorinated biphenyl (PCB). Copper treatments (CuCl2 and CuSO4) did not significantly induce the expression of ApGSTs. The highest up-regulations of ApGST1 and ApGST2 were under 6-h treatments of 0.10 and 0.50 mg L-1 NaOCl. Interestingly, only ApGST1 increased significantly after 0.10, 0.50, and 1.00 mg L-1 of PCB exposure (6 h). Intracellular reactive oxygen species (ROS) increased considerably under NaOCl; however, it was not significantly higher in the PCB-treated cells. GST activity was increased by NaOCl and PCB treatments, but only PCB caused apoptosis. These results suggest that GSTs are involved in the first line of phase II detoxification, protecting dinoflagellate cells against oxidative damage.

Genome-wide transcriptome analysis reveals the diversity and function of long non-coding RNAs in dinoflagellates

Dinoflagellates are a diverse group of phytoplankton, ranging from harmful bloom-forming microalgae to photosymbionts of coral reefs. Genome-scale data from dinoflagellates reveal atypical genomic features, extensive genomic divergence, and lineage-specific innovation of gene functions. Long non-coding RNAs (lncRNAs), known to regulate gene expression in eukaryotes, are largely unexplored in dinoflagellates. Here, using high-quality genome and transcriptome data, we identified 48039 polyadenylated lncRNAs in three dinoflagellate species: the coral symbionts Cladocopium proliferum and Durusdinium trenchii, and the bloom-forming species, Prorocentrum cordatum. These lncRNAs have fewer introns and lower G+C content than protein-coding sequences; 37 768 (78.6%) are unique with respect to sequence similarity. We classified all lncRNAs based on conserved motifs (k-mers) into distinct clusters, following properties of protein-binding and/or subcellular localisation. Interestingly, 3708 (7.7%) lncRNAs are differentially expressed under heat stress, algal lifestyle, and/or growth phase, and share co-expression patterns with protein-coding genes. Based on inferred triplex interactions between lncRNA and putative promoter regions, we identified 19 460 putative gene targets for 3721 lncRNAs; 907 genes exhibit differential expression under heat stress. These results reveal, for the first time, the diversity of lncRNAs in dinoflagellates and how lncRNAs may regulate gene expression as a heat-stress response in these ecologically important microbes.

Multi-Chemical Omics Analysis of the Symbiodiniaceae Durusdinium trenchii under Heat Stress

The urgency of responding to climate change for corals necessitates the exploration of innovative methods to swiftly enhance our understanding of crucial processes. In this study, we employ an integrated chemical omics approach, combining elementomics, metabolomics, and volatilomics methodologies to unravel the biochemical pathways associated with the thermal response of the coral symbiont, Symbiodiniaceae Durusdinium trenchii. We outline the complimentary sampling approaches and discuss the standardised data corrections used to allow data integration and comparability. Our findings highlight the efficacy of individual methods in discerning differences in the biochemical response of D. trenchii under both control and stress-inducing temperatures. However, a deeper insight emerges when these methods are integrated, offering a more comprehensive understanding, particularly regarding oxidative stress pathways. Employing correlation network analysis enhanced the interpretation of volatile data, shedding light on the potential metabolic origins of volatiles with undescribed functions and presenting promising candidates for further exploration. Elementomics proves to be less straightforward to integrate, likely due to no net change in elements but rather elements being repurposed across compounds. The independent and integrated data from this study informs future omic profiling studies and recommends candidates for targeted research beyond Symbiodiniaceae biology. This study highlights the pivotal role of omic integration in advancing our knowledge, addressing critical gaps, and guiding future research directions in the context of climate change and coral reef preservation.

Protocol for the generation of Symbiodiniaceae mutants using UV mutagenesis

Genetic approaches are limited in the dinoflagellate family, Symbiodiniaceae, causing a bottleneck in the discovery of useful mutants toward the goal of preventing future coral bleaching events. In this protocol, we demonstrate the application of UV exposure, coupled with downstream phenotypic screening and mutant isolation, to form a UV mutagenesis pipeline. This pipeline provides an avenue to generate Symbiodiniaceae mutants to help link genotype to phenotype, as well as address previously unanswered questions surrounding relationships with host organisms, like coral. For complete details on the use and execution of this protocol, please refer to Jinkerson et al. (2022).1.

The enigmatic nucleus of the marine dinoflagellate Prorocentrum cordatum

The marine, bloom-forming dinoflagellate Prorocentrum cordatum CCMP 1329 (formerly P. minimum) has a genome atypical of eukaryotes, with a large size of ~4.15 Gbp, organized in plentiful, highly condensed chromosomes and packed in a dinoflagellate-specific nucleus (dinokaryon). Here, we apply microscopic and proteogenomic approaches to obtain new insights into this enigmatic nucleus of axenic P. cordatum. High-resolution focused ion beam/scanning electron microscopy analysis of the flattened nucleus revealed highest density of nuclear pores in the vicinity of the nucleolus, a total of 62 tightly packed chromosomes (~0.4-6.7 µm3), and interaction of several chromosomes with the nucleolus and other nuclear structures. A specific procedure for enriching intact nuclei was developed to enable proteomic analyses of soluble and membrane protein-enriched fractions. These were analyzed with geLC and shotgun approaches employing ion-trap and timsTOF (trapped-ion-mobility-spectrometry time-of-flight) mass spectrometers, respectively. This allowed identification of 4,052 proteins (39% of unknown function), out of which 418 were predicted to serve specific nuclear functions; additional 531 proteins of unknown function could be allocated to the nucleus. Compaction of DNA despite very low histone abundance could be accomplished by highly abundant major basic nuclear proteins (HCc2-like). Several nuclear processes including DNA replication/repair and RNA processing/splicing can be fairly well explained on the proteogenomic level. By contrast, transcription and composition of the nuclear pore complex remain largely elusive. One may speculate that the large group of potential nuclear proteins with currently unknown functions may serve yet to be explored functions in nuclear processes differing from those of typical eukaryotic cells. IMPORTANCE Dinoflagellates form a highly diverse group of unicellular microalgae. They provide keystone species for the marine ecosystem and stand out among others by their very large, unusually organized genomes embedded in the nuclei markedly different from other eukaryotic cells. Functional insights into nuclear and other cell biological structures and processes of dinoflagellates have long been hampered by the paucity of available genomic sequences. The here studied cosmopolitan P. cordatum belongs to the harmful algal bloom-forming, marine dinoflagellates and has a recently de novo assembled genome. We present a detailed 3D reconstruction of the P. cordatum nucleus together with comprehensive proteogenomic insights into the protein equipment mastering the broad spectrum of nuclear processes. This study significantly advances our understanding of mechanisms and evolution of the conspicuous dinoflagellate cell biology.

Stony coral tissue loss disease induces transcriptional signatures of in situ degradation of dysfunctional Symbiodiniaceae

Stony coral tissue loss disease (SCTLD), one of the most pervasive and virulent coral diseases on record, affects over 22 species of reef-building coral and is decimating reefs throughout the Caribbean. To understand how different coral species and their algal symbionts (family Symbiodiniaceae) respond to this disease, we examine the gene expression profiles of colonies of five species of coral from a SCTLD transmission experiment. The included species vary in their purported susceptibilities to SCTLD, and we use this to inform gene expression analyses of both the coral animal and their Symbiodiniaceae. We identify orthologous coral genes exhibiting lineage-specific differences in expression that correlate to disease susceptibility, as well as genes that are differentially expressed in all coral species in response to SCTLD infection. We find that SCTLD infection induces increased expression of rab7, an established marker of in situ degradation of dysfunctional Symbiodiniaceae, in all coral species accompanied by genus-level shifts in Symbiodiniaceae photosystem and metabolism gene expression. Overall, our results indicate that SCTLD infection induces symbiophagy across coral species and that the severity of disease is influenced by Symbiodiniaceae identity.

Transcriptomic analysis of polyketide synthesis in dinoflagellate, Prorocentrum lima

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

Genomic conservation and putative downstream functionality of the phosphatidylinositol signalling pathway in the cnidarian-dinoflagellate symbiosis

The mutualistic cnidarian-dinoflagellate symbiosis underpins the evolutionary success of stony corals and the persistence of coral reefs. However, a molecular understanding of the signalling events that lead to the successful establishment and maintenance of this symbiosis remains unresolved. For example, the phosphatidylinositol (PI) signalling pathway has been implicated during the establishment of multiple mutualistic and parasitic interactions across the kingdoms of life, yet its role within the cnidarian-dinoflagellate symbiosis remains unexplored. Here, we aimed to confirm the presence and assess the specific enzymatic composition of the PI signalling pathway across cnidaria and dinoflagellates by compiling 21 symbiotic anthozoan (corals and sea anemones) and 28 symbiotic dinoflagellate (Symbiodiniaceae) transcriptomic and genomic datasets and querying genes related to this pathway. Presence or absence of PI-kinase and PI-phosphatase orthologs were also compared between a broad sampling of taxonomically related symbiotic and non-symbiotic species. Across the symbiotic anthozoans analysed, there was a complete and highly conserved PI pathway, analogous to the pathway found in model eukaryotes. The Symbiodiniaceae pathway showed similarities to its sister taxon, the Apicomplexa, with the absence of PI 4-phosphatases. However, conversely to Apicomplexa, there was also an expansion of homologs present in the PI5-phosphatase and PI5-kinase groups, with unique Symbiodiniaceae proteins identified that are unknown from non-symbiotic unicellular organisms. Additionally, we aimed to unravel the putative functionalities of the PI signalling pathway in this symbiosis by analysing phosphoinositide (PIP)-binding proteins. Analysis of phosphoinositide (PIP)-binding proteins showed that, on average, 2.23 and 1.29% of the total assemblies of anthozoan and Symbiodiniaceae, respectively, have the potential to bind to PIPs. Enrichment of Gene Ontology (GO) terms associated with predicted PIP-binding proteins within each taxon revealed a broad range of functions, including compelling links to processes putatively involved in symbiosis regulation. This analysis establishes a baseline for current understanding of the PI pathway across anthozoans and Symbiodiniaceae, and thus a framework to target future research.

Building consensus around the assessment and interpretation of Symbiodiniaceae diversity

Within microeukaryotes, genetic variation and functional variation sometimes accumulate more quickly than morphological differences. To understand the evolutionary history and ecology of such lineages, it is key to examine diversity at multiple levels of organization. In the dinoflagellate family Symbiodiniaceae, which can form endosymbioses with cnidarians (e.g., corals, octocorals, sea anemones, jellyfish), other marine invertebrates (e.g., sponges, molluscs, flatworms), and protists (e.g., foraminifera), molecular data have been used extensively over the past three decades to describe phenotypes and to make evolutionary and ecological inferences. Despite advances in Symbiodiniaceae genomics, a lack of consensus among researchers with respect to interpreting genetic data has slowed progress in the field and acted as a barrier to reconciling observations. Here, we identify key challenges regarding the assessment and interpretation of Symbiodiniaceae genetic diversity across three levels: species, populations, and communities. We summarize areas of agreement and highlight techniques and approaches that are broadly accepted. In areas where debate remains, we identify unresolved issues and discuss technologies and approaches that can help to fill knowledge gaps related to genetic and phenotypic diversity. We also discuss ways to stimulate progress, in particular by fostering a more inclusive and collaborative research community. We hope that this perspective will inspire and accelerate coral reef science by serving as a resource to those designing experiments, publishing research, and applying for funding related to Symbiodiniaceae and their symbiotic partnerships.

The enigmatic clock of dinoflagellates, is it unique?

Dinoflagellate clocks are unique as they show no resemblance to any known model eukaryotic or prokaryotic clock architecture. Dinoflagellates are unicellular, photosynthetic, primarily marine eukaryotes are known for their unique biology and rhythmic physiology. Their physiological rhythms are driven by an internal oscillator whose molecular underpinnings are yet unknown. One of the primary reasons that slowed the progression of their molecular studies is their extremely large and repetitive genomes. Dinoflagellates are primary contributors to the global carbon cycle and oxygen levels, therefore, comprehending their internal clock architecture and its interaction with their physiology becomes a subject of utmost importance. The advent of high throughput Omics technology provided the momentum to understand the molecular architecture and functioning of the dinoflagellate clocks. We use these extensive databases to perform meta-analysis to reveal the status of clock components in dinoflagellates. In this article, we will delve deep into the various “Omics” studies that catered to various breakthroughs in the field of circadian biology in these organisms that were not possible earlier. The overall inference from these omics studies points toward an uncommon eukaryotic clock model, which can provide promising leads to understand the evolution of molecular clocks.

Temporal gene expression patterns in the coral Euphyllia paradivisa reveal the complexity of biological clocks in the cnidarian-algal symbiosis

Studying chronobiology in reef-building corals is challenging due to the tightly coupled symbiosis with their photosynthetic algae, Symbiodiniaceae. Although symbiosis requires metabolic synchronization and coordination of cellular processes in the holobiont, the cross-talk between the host and symbiont's clocks is still puzzling. Here, we use the mesophotic coral Euphyllia paradivisa to examine temporal gene expression patterns in symbiotic and aposymbiotic morphs exposed to natural light/dark cycles and constant darkness. Our comparative transcriptomic analyses revealed circadian and circatidal cycles of gene expression with a predominant diel pattern in both coral morphs. We found a substantial number of transcripts consistently rhythmic under both light conditions, including genes likely involved in the cnidarians' circadian clock, thus indicating that an endogenous clock, which can oscillate independently from the Symbiodiniaceae clock, exists in E. paradivisa. The analysis further manifests the remarkable impacts of symbiosis on transcriptional rhythms and implies that the algae's presence influences the host's biorhythm.

An overview of transcription in dinoflagellates

Dinoflagellates are a vital diverse family of unicellular algae widespread in various aquatic environments. Typically large genomes and permanently condensed chromosomes without histones make these organisms unique among eukaryotes in terms of chromatin structure and gene expression. Genomic and transcriptomic sequencing projects have provided new insight into the genetic foundation of dinoflagellate behaviors. Genes in tandem arrays, trans-splicing of mRNAs and lower levels of transcriptional regulation compared to other eukaryotes all contribute to the differences seen. Here we present a general overview of transcription in dinoflagellates based on previously described work.

Mechanisms and potential immune tradeoffs of accelerated coral growth induced by microfragmentation

Microfragmentation is the act of cutting corals into small pieces (~1 cm2) to accelerate the growth rates of corals relative to growth rates observed when maintaining larger-sized fragments. This rapid tissue and skeletal expansion technique offers great potential for supporting reef restoration, yet the biological processes and tradeoffs involved in microfragmentation-mediated accelerated growth are not well understood. Here we compared growth rates across a range of successively smaller fragment sizes in multiple genets of reef-building corals, Orbicella faveolata and Montastraea cavernosa. Our results confirm prior findings that smaller initial sizes confer accelerated growth after four months of recovery in a raceway. O. faveolata transcript levels associated with growth rate include genes encoding carbonic anhydrase and glutamic acid-rich proteins, which have been previously implicated in coral biomineralization, as well as a number of unannotated transcripts that warrant further characterization. Innate immunity enzyme activity assays and gene expression results suggest a potential tradeoff between growth rate after microfragmentation and immune investment. Microfragmentation-based restoration practices have had great success on Caribbean reefs, despite widespread mortality among wild corals due to infectious diseases. Future studies should continue to examine potential immune tradeoffs throughout the microfragmentation recovery period that may affect growout survival and disease transmission after outplanting.

Alphaflexivirus Genomes in Stony Coral Tissue Loss Disease-Affected, Disease-Exposed, and Disease-Unexposed Coral Colonies in the U.S. Virgin Islands

Stony coral tissue loss disease (SCTLD) is decimating Caribbean corals. Here, through the metatranscriptomic assembly and annotation of two alphaflexivirus-like strains, we provide genomic evidence of filamentous viruses in SCTLD-affected, -exposed, and -unexposed coral colonies. These data will assist in clarifying the roles of viruses in SCTLD.

Gene clusters for biosynthesis of mycosporine-like amino acids in dinoflagellate nuclear genomes: Possible recent horizontal gene transfer between species of Symbiodiniaceae (Dinophyceae)

Global warming increases the temperature of the ocean surface, which can disrupt dinoflagellate-coral symbioses and result in coral bleaching. Photosynthetic dinoflagellates of the family Symbiodiniaceae include bleaching-tolerant and bleaching-sensitive coral symbionts. Therefore, understanding the molecular mechanisms for changing symbiont diversity is potentially useful to assist recovery of coral holobionts (corals and their associated microbes, including multiple species of Symbiodiniaceae), although sexual reproduction has not been observed in the Symbiodiniaceae. Recent molecular phylogenetic analyses estimate that the Symbiodiniaceae appeared 160 million years ago and diversified into 15 groups, five genera of which now have available draft genomes (i.e., Symbiodinium, Durusdinium, Breviolum, Fugacium, and Cladocopium). Comparative genomic analyses have suggested that crown groups have fewer gene families than early-diverging groups, although many genes that were probably acquired via gene duplications and horizontal gene transfers (HGTs) have been found in each decoded genome. Because UV stress is likely a contributor to coral bleaching, and because the highly conserved gene cluster for mycosporine-like amino acid (MAA) biosynthesis has been found in thermal-tolerant symbiont genomes, I reviewed genomic features of the Symbiodiniaceae, focusing on possible acquisition of a biosynthetic gene cluster for MAAs, which absorb UV radiation. On the basis of highly conserved noncoding sequences, I hypothesized that HGTs have occurred among members of the Symbiodiniaceae and have contributed to the diversification of Symbiodiniaceae-host relationships. Finally, I proposed that bleaching tolerance may be strengthened by multiple MAAs from both symbiotic dinoflagellates and corals.

The tropical coral Pocillopora acuta displays an unusual chromatin structure and shows histone H3 clipping plasticity upon bleaching

Background: Pocillopora acuta is a hermatypic coral with strong ecological importance. Anthropogenic disturbances and global warming are major threats that can induce coral bleaching, the disruption of the mutualistic symbiosis between the coral host and its endosymbiotic algae. Previous works have shown that somaclonal colonies display different levels of survival depending on the environmental conditions they previously faced. Epigenetic mechanisms are good candidates to explain this phenomenon. However, almost no work had been published on the P. acuta epigenome, especially on histone modifications. In this study, we aim at providing the first insight into chromatin structure of this species. Methods: We aligned the amino acid sequence of P. acuta core histones with histone sequences from various phyla. We developed a centri-filtration on sucrose gradient to separate chromatin from the host and the symbiont. The presence of histone H3 protein and specific histone modifications were then detected by western blot performed on histone extraction done from bleached and healthy corals. Finally, micrococcal nuclease (MNase) digestions were undertaken to study nucleosomal organization. Results: The centri-filtration enabled coral chromatin isolation with less than 2% of contamination by endosymbiont material. Histone sequences alignments with other species show that P. acuta displays on average ~90% of sequence similarities with mice and ~96% with other corals. H3 detection by western blot showed that H3 is clipped in healthy corals while it appeared to be intact in bleached corals. MNase treatment failed to provide the usual mononucleosomal digestion, a feature shared with some cnidarian, but not all; suggesting an unusual chromatin structure. Conclusions: These results provide a first insight into the chromatin, nucleosome and histone structure of P. acuta. The unusual patterns highlighted in this study and partly shared with other cnidarian will need to be further studied to better understand its role in corals.

The specific inhibition of glycerol synthesis and the phosphorylation of a putative MAPK give insight into the mechanism of osmotic sensing in a dinoflagellate symbiont

Signaling pathways are fundamental for the establishment and maintenance of diverse symbioses. The symbiosis of cnidarians and dinoflagellate algae is the foundation for the ecological success of coral reefs, involving the transfer of photosynthetic products from symbiont to host. However, signal transduction pathways for this symbiosis remain uncharacterized. Cultured and natural cnidarian symbionts can produce glycerol, one of the main translocated photosynthates. Here, we investigate whether a signal transduction pathway may be involved in inducing glycerol synthesis in cultured symbionts under an osmotic stress model. We evaluated the effect of specific inhibitors of the main transduction pathways, p38, JNK, and ERK 1/2 in Brevolium minutum, the symbiont of the Aiptasia model system. We found that glycerol production and the specific activity of the enzyme Gpdh were selectively inhibited by a p38 MAPK inhibitor. Additionally, the phosphorylation of a putative p38-like protein was rapidly detected. Finally, we studied the presence of each of the components of the p38 MAPK pathway in silico, in genomes and transcriptomes reported up to date for different symbiont types. We propose a model for the arrangement of this pathway in the family of dinoflagellate symbionts known as Symbiodiniaceae.

Contrasting heat stress response patterns of coral holobionts across the Red Sea suggest distinct mechanisms of thermal tolerance

Corals from the northern Red Sea, in particular the Gulf of Aqaba (GoA), have exceptionally high bleaching thresholds approaching >5℃ above their maximum monthly mean (MMM) temperatures. These elevated thresholds are thought to be due to historical selection, as corals passed through the warmer Southern Red Sea during recolonization from the Arabian Sea. To test this hypothesis, we determined thermal tolerance thresholds of GoA versus central Red Sea (CRS) Stylophora pistillata corals using multi-temperature acute thermal stress assays to determine thermal thresholds. Relative thermal thresholds of GoA and CRS corals were indeed similar and exceptionally high (~7℃ above MMM). However, absolute thermal thresholds of CRS corals were on average 3℃ above those of GoA corals. To explore the molecular underpinnings, we determined gene expression and microbiome response of the coral holobiont. Transcriptomic responses differed markedly, with a strong response to the thermal stress in GoA corals and their symbiotic algae versus a remarkably muted response in CRS colonies. Concomitant to this, coral and algal genes showed temperature-induced expression in GoA corals, while exhibiting fixed high expression (front-loading) in CRS corals. Bacterial community composition of GoA corals changed dramatically under heat stress, whereas CRS corals displayed stable assemblages. We interpret the response of GoA corals as that of a resilient population approaching a tipping point in contrast to a pattern of consistently elevated thermal resistance in CRS corals that cannot further attune. Such response differences suggest distinct thermal tolerance mechanisms that may affect the response of coral populations to ocean warming.

The tropical coral Pocillopora acuta displays an unusual chromatin structure and shows histone H3 clipping plasticity upon bleaching

Background: Pocillopora acuta is a hermatypic coral with strong ecological importance. Anthropogenic disturbances and global warming are major threats that can induce coral bleaching, the disruption of the mutualistic symbiosis between the coral host and its endosymbiotic algae. Previous works have shown that somaclonal colonies display different levels of survival depending on the environmental conditions they previously faced. Epigenetic mechanisms are good candidates to explain this phenomenon. However, almost no work had been published on the P. acuta epigenome, especially on histone modifications. In this study, we aim at providing the first insight into chromatin structure of this species. Methods: We aligned the amino acid sequence of P. acuta core histones with histone sequences from various phyla. We developed a centri-filtration on sucrose gradient to separate chromatin from the host and the symbiont. The presence of histone H3 protein and specific histone modifications were then detected by western blot performed on histone extraction done from bleached and healthy corals. Finally, micrococcal nuclease (MNase) digestions were undertaken to study nucleosomal organization. Results: The centri-filtration enabled coral chromatin isolation with less than 2% of contamination by endosymbiont material. Histone sequences alignments with other species show that P. acuta displays on average ~90% of sequence similarities with mice and ~96% with other corals. H3 detection by western blot showed that H3 is clipped in healthy corals while it appeared to be intact in bleached corals. MNase treatment failed to provide the usual mononucleosomal digestion, a feature shared with some cnidarian, but not all; suggesting an unusual chromatin structure. Conclusions: These results provide a first insight into the chromatin, nucleosome and histone structure of P. acuta. The unusual patterns highlighted in this study and partly shared with other cnidarian will need to be further studied to better understand its role in corals.

Evaluation of Filter, Paramagnetic, and STAGETips Aided Workflows for Proteome Profiling of Symbiodiniaceae Dinoflagellate

The integrity of coral reef ecosystems worldwide rests on a fine-tuned symbiotic interaction between an invertebrate and a dinoflagellate microalga from the family Symbiodiniaceae. Recent advances in bottom-up shotgun proteomic approaches and the availability of vast amounts of genetic information about Symbiodiniaceae have provided a unique opportunity to better understand the molecular mechanisms underpinning the interactions of coral-Symbiodiniaceae. However, the resilience of this dinoflagellate cell wall, as well as the presence of polyanionic and phenolics cell wall components, requires the optimization of sample preparation techniques for successful implementation of bottom-up proteomics. Therefore, in this study we compare three different workflows—filter-aided sample preparation (FASP), single-pot solid-phase-enhanced sample preparation (SP3), and stop-and-go-extraction tips (STAGETips, ST)—to develop a high-throughput proteotyping protocol for Symbiodiniaceae algal research. We used the model isolate Symbiodinium tridacnidorum. We show that SP3 outperformed ST and FASP with regard to robustness, digestion efficiency, and contaminant removal, which led to the highest number of total (3799) and unique proteins detected from 23,593 peptides. Most of these proteins were detected with ≥2 unique peptides (73%), zero missed tryptic peptide cleavages (91%), and hydrophilic peptides (>70%). To demonstrate the functionality of this optimized SP3 sample preparation workflow, we examined the proteome of S. tridacnidorum to better understand the molecular mechanism of peridinin-chlorophyll-protein complex (PCP, light harvesting protein) accumulation under low light (LL, 30 μmol photon m−2 s−1). Cells exposed to LL for 7 days upregulated various light harvesting complex (LHCs) proteins through the mevalonate-independent pathway; proteins of this pathway were at 2- to 6-fold higher levels than the control of 120 μmol photon m−2 s−1. Potentially, LHCs which were maintained in an active phosphorylated state by serine/threonine-protein kinase were also upregulated to 10-fold over control. Collectively, our results show that the SP3 method is an efficient high-throughput proteotyping tool for Symbiodiniaceae algal research.

A stony coral cell atlas illuminates the molecular and cellular basis of coral symbiosis, calcification, and immunity

Stony corals are colonial cnidarians that sustain the most biodiverse marine ecosystems on Earth: coral reefs. Despite their ecological importance, little is known about the cell types and molecular pathways that underpin the biology of reef-building corals. Using single-cell RNA sequencing, we define over 40 cell types across the life cycle of Stylophora pistillata. We discover specialized immune cells, and we uncover the developmental gene expression dynamics of calcium-carbonate skeleton formation. By simultaneously measuring the transcriptomes of coral cells and the algae within them, we characterize the metabolic programs involved in symbiosis in both partners. We also trace the evolution of these coral cell specializations by phylogenetic integration of multiple cnidarian cell type atlases. Overall, this study reveals the molecular and cellular basis of stony coral biology.

Assessing nucleic acid binding activity of four dinoflagellate cold shock domain proteins from Symbiodinium kawagutii and Lingulodinium polyedra

Background

Dinoflagellates have a generally large number of genes but only a small percentage of these are annotated as transcription factors. Cold shock domain (CSD) containing proteins (CSPs) account for roughly 60% of these. CSDs are not prevalent in other eukaryotic lineages, perhaps suggesting a lineage-specific expansion of this type of transcription factors in dinoflagellates, but there is little experimental data to support a role for dinoflagellate CSPs as transcription factors. Here we evaluate the hypothesis that dinoflagellate CSPs can act as transcription factors by binding double-stranded DNA in a sequence dependent manner.

Results

We find that both electrophoretic mobility shift assay (EMSA) competition experiments and selection and amplification binding (SAAB) assays indicate binding is not sequence specific for four different CSPs from two dinoflagellate species. Competition experiments indicate all four CSPs bind to RNA better than double-stranded DNA.

Conclusion

Dinoflagellate CSPs do not share the nucleic acid binding properties expected for them to function as bone fide transcription factors. We conclude the transcription factor complement of dinoflagellates is even smaller than previously thought suggesting that dinoflagellates have a reduced dependance on transcriptional control compared to other eukaryotes.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12860-021-00368-4.

Gene expression associated with disease resistance and long-term growth in a reef-building coral

Rampant coral disease, exacerbated by climate change and other anthropogenic stressors, threatens reefs worldwide, especially in the Caribbean. Physically isolated yet genetically connected reefs such as Flower Garden Banks National Marine Sanctuary (FGBNMS) may serve as potential refugia for degraded Caribbean reefs. However, little is known about the mechanisms and trade-offs of pathogen resistance in reef-building corals. Here, we measure pathogen resistance in Montastraea cavernosa from FGBNMS. We identified individual colonies that demonstrated resistance or susceptibility to Vibrio spp. in a controlled laboratory environment. Long-term growth patterns suggest no trade-off between disease resistance and calcification. Predictive (pre-exposure) gene expression highlights subtle differences between resistant and susceptible genets, encouraging future coral disease studies to investigate associations between resistance and replicative age and immune cell populations. Predictive gene expression associated with long-term growth underscores the role of transmembrane proteins involved in cell adhesion and cell-cell interactions, contributing to the growing body of knowledge surrounding genes that influence calcification in reef-building corals. Together these results demonstrate that coral genets from isolated sanctuaries such as FGBNMS can withstand pathogen challenges and potentially aid restoration efforts in degraded reefs. Furthermore, gene expression signatures associated with resistance and long-term growth help inform strategic assessment of coral health parameters.

Comparison of 15 dinoflagellate genomes reveals extensive sequence and structural divergence in family Symbiodiniaceae and genus Symbiodinium

Background

Dinoflagellates in the family Symbiodiniaceae are important photosynthetic symbionts in cnidarians (such as corals) and other coral reef organisms. Breakdown of the coral-dinoflagellate symbiosis due to environmental stress (i.e. coral bleaching) can lead to coral death and the potential collapse of reef ecosystems. However, evolution of Symbiodiniaceae genomes, and its implications for the coral, is little understood. Genome sequences of Symbiodiniaceae remain scarce due in part to their large genome sizes (1–5 Gbp) and idiosyncratic genome features.

Results

Here, we present de novo genome assemblies of seven members of the genus Symbiodinium, of which two are free-living, one is an opportunistic symbiont, and the remainder are mutualistic symbionts. Integrating other available data, we compare 15 dinoflagellate genomes revealing high sequence and structural divergence. Divergence among some Symbiodinium isolates is comparable to that among distinct genera of Symbiodiniaceae. We also recovered hundreds of gene families specific to each lineage, many of which encode unknown functions. An in-depth comparison between the genomes of the symbiotic Symbiodinium tridacnidorum (isolated from a coral) and the free-living Symbiodinium natans reveals a greater prevalence of transposable elements, genetic duplication, structural rearrangements, and pseudogenisation in the symbiotic species.

Conclusions

Our results underscore the potential impact of lifestyle on lineage-specific gene-function innovation, genome divergence, and the diversification of Symbiodinium and Symbiodiniaceae. The divergent features we report, and their putative causes, may also apply to other microbial eukaryotes that have undergone symbiotic phases in their evolutionary history.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12915-021-00994-6.

Genotype by environment interactions in coral bleaching

Climate-driven reef decline has prompted the development of next-generation coral conservation strategies, many of which hinge on the movement of adaptive variation across genetic and environmental gradients. This process is limited by our understanding of how genetic and genotypic drivers of coral bleaching will manifest in different environmental conditions. We reciprocally transplanted 10 genotypes of Acropora cervicornis across eight sites along a 60 km span of the Florida Reef Tract and documented significant genotype × environment interactions in bleaching response during the severe 2015 bleaching event. Performance relative to site mean was significantly different between genotypes and can be mostly explained by ensemble models of correlations with genetic markers. The high explanatory power was driven by significant enrichment of loci associated DNA repair, cell signalling and apoptosis. No genotypes performed above (or below) bleaching average at all sites, so genomic predictors can provide practitioners with 'confidence intervals' about the chance of success in novel habitats. These data have important implications for assisted gene flow and managed relocation, and their integration with traditional active restoration.

Using glutamine synthetase 1 to evaluate the symbionts' potential of ammonia assimilation and their responses to illumination in five organs of the giant clam, Tridacna squamosa

Nitrogen-deficient symbiotic dinoflagellates (zooxanthellae) living inside the fluted giant clam, Tridacna squamosa, need to obtain nitrogen from the host. Glutamine synthetase 1 (GS1) is a cytosolic enzyme that assimilates ammonia into glutamine. We determined the transcript levels of zooxanthellal GS1 (Zoox-GS1), which represented comprehensively GS1 transcripts of Symbiodinium, Cladocopium and Durusdinium, in five organs of T. squamosa. The outer mantle had significantly higher transcript level of Zoox-GS1 than the inner mantle, foot muscle, hepatopancreas and ctenidium, but the transcript ratios of Zoox-GS1 to zooxanthellal form II ribulose-1,5-bisphosphate carboxylase/oxygenase (Zoox-rbcII), which represented the potential of ammonia assimilation relative to the phototrophic potential, were comparable among these five organs. Based on transcript ratios of Zoox-GS1 to zooxanthellal Urease (Zoox-URE), the outer mantle had the highest potential of urea degradation relative to ammonia assimilation among the five organs, probably because urea degradation could furnish CO₂ and NH₃ for photosynthesis and amino acid synthesis, respectively, in the symbionts therein. The protein abundance of Zoox-GS1 was upregulated in the outer mantle and the inner mantle during illumination. Zoox-GS1 could catalyze light-enhanced glutamine formation using ammonia absorbed from the host or ammonia released through urea degradation in the cytoplasm. The glutamine produced could be used to synthesize other nitrogenous compounds, including amino acids in the cytoplasm or in the plastid of the dinoflagellates. Some of the amino acids synthesized by the symbionts in the inner mantle and foot muscle could be donated to the host to support shell organic matrix formation and muscle production, respectively.

Genetic and spatial organization of the unusual chromosomes of the dinoflagellate Symbiodinium microadriaticum

Dinoflagellates are main primary producers in the oceans, the cause of algal blooms and endosymbionts of marine invertebrates. Much remains to be understood about their biology, including their peculiar crystalline chromosomes. We assembled 94 chromosome-scale scaffolds of the genome of the coral endosymbiont Symbiodinium microadriaticum and analyzed their organization. Genes are enriched towards the ends of chromosomes and are arranged in alternating unidirectional blocks. Some chromosomes are enriched for genes involved in specific biological processes. The chromosomes fold as linear rods and each is composed of a series of structural domains separated by boundaries. Domain boundaries are positioned at sites where transcription of two gene blocks converges and disappear when cells are treated with chemicals that block transcription, indicating correlations between gene orientation, transcription and chromosome folding. The description of the genetic and spatial organization of the S. microadriaticum genome provides a foundation for deeper exploration of the extraordinary biology of dinoflagellates and their chromosomes.

Symbiodiniaceae Dinoflagellates Express Urease in Three Subcellular Compartments and Upregulate its Expression Levels in situ in Three Organs of a Giant Clam (Tridacna squamosa) During Illumination

Giant clams harbor three genera of symbiotic dinoflagellates (Symbiodinium, Cladocopium, and Durusdinium) as extracellular symbionts (zooxanthellae). While symbiotic dinoflagellates can synthesize amino acids to benefit the host, they are nitrogen-deficient. Hence, the host must supply them with nitrogen including urea, which can be degraded to ammonia and carbon dioxide by urease (URE). Here, we report three complete coding cDNA sequences of URE, one for each genus of dinoflagellate, obtained from the colorful outer mantle of the giant clam, Tridacna squamosa. The outer mantle had higher transcript level of Tridacna squamosa zooxanthellae URE (TSZURE) than the whitish inner mantle, foot muscle, hepatopancreas, and ctenidium. TSZURE was immunolocalized strongly and atypically in the plastid, moderately in the cytoplasm, and weakly in the cell wall and plasma membrane of symbiotic dinoflagellates. In the outer mantle, illumination upregulated the protein abundance of TSZURE, which could enhance urea degradation in photosynthesizing dinoflagellates. The urea-nitrogen released could then augment synthesis of amino acids to be shared with the host for its general needs. Illumination also enhanced gene and protein expression levels of TSZURE/TSZURE in the inner mantle and foot muscle, which contain only small quantities of symbiotic dinoflagellate, have no iridocyte, and lack direct exposure to light. With low phototrophic potential, dinoflagellates in the inner mantle and foot muscle might need to absorb carbohydrates in order to assimilate the urea-nitrogen into amino acids. Amino acids donated by dinoflagellates to the inner mantle and the foot muscle could be used especially for synthesis of organic matrix needed for light-enhanced shell formation and muscle protein, respectively.

Innate immune gene expression in Acropora palmata is consistent despite variance in yearly disease events

Coral disease outbreaks are expected to increase in prevalence, frequency and severity due to climate change and other anthropogenic stressors. This is especially worrying for the Caribbean branching coral Acropora palmata which has already seen an 80% decrease in cover primarily due to disease. Despite the importance of this keystone species, there has yet to be a characterization of its transcriptomic response to disease exposure. In this study we provide the first transcriptomic analysis of 12 A. palmata genotypes and their symbiont Symbiodiniaceae exposed to disease in 2016 and 2017. Year was the primary driver of gene expression variance for A. palmata and the Symbiodiniaceae. We hypothesize that lower expression of ribosomal genes in the coral, and higher expression of transmembrane ion transport genes in the Symbiodiniaceae indicate that a compensation or dysbiosis may be occurring between host and symbiont. Disease response was the second driver of gene expression variance for A. palmata and included a core set of 422 genes that were significantly differentially expressed. Of these, 2 genes (a predicted cyclin-dependent kinase 11b and aspartate 1-decarboxylase) showed negative Log2 fold changes in corals showing transmission of disease, and positive Log2 fold changes in corals showing no transmission of disease, indicating that these may be important in disease resistance. Co-expression analysis identified two modules positively correlated to disease exposure, one enriched for lipid biosynthesis genes, and the other enriched in innate immune genes. The hub gene in the immune module was identified as D-amino acid oxidase, a gene implicated in phagocytosis and microbiome homeostasis. The role of D-amino acid oxidase in coral immunity has not been characterized but could be an important enzyme for responding to disease. Our results indicate that A. palmata mounts a core immune response to disease exposure despite differences in the disease type and virulence between 2016 and 2017. These identified genes may be important for future biomarker development in this Caribbean keystone species.

Population genetic structure of the great star coral, Montastraea cavernosa, across the Cuban archipelago with comparisons between microsatellite and SNP markers

Coral reef habitats surrounding Cuba include relatively healthy, well-developed shallow and mesophotic (30–150 m) scleractinian communities at the cross-currents of the Tropical Western Atlantic (TWA). However, Cuba’s coral communities are not immune to the declines observed throughout the TWA, and there is limited information available regarding genetic connectivity, diversity, and structure among these populations. This represents an immense gap in our understanding of coral ecology and population dynamics at both local and regional scales. To address this gap, we evaluated the population genetic structure of the coral Montastraea cavernosa across eight reef sites surrounding Cuba. Colonies were genotyped using nine microsatellite markers and > 9,000 single nucleotide polymorphism (SNP) markers generated using the 2bRAD approach to assess fine-scale genetic structure across these sites. Both the microsatellite and SNP analyses identified patterns of genetic differentiation among sample populations. While the microsatellite analyses did not identify significant genetic structure across the seven shallow M. cavernosa sampling sites, the SNP analyses revealed significant pairwise population differentiation, suggesting that differentiation is greater between eastern and western sites. This study provides insight into methodological differences between microsatellite and SNP markers including potential trade-offs between marker-specific biases, sample size, sequencing costs, and the ability to resolve subtle patterns of population genetic structure. Furthermore, this study suggests that locations in western Cuba may play important roles in this species’ regional metapopulation dynamics and therefore may merit incorporation into developing international management efforts in addition to the local management the sites receive.

Comparative transcriptomic analyses of Chromera and Symbiodiniaceae

Reef-building corals live in a mutualistic relationship with photosynthetic algae (family Symbiodiniaceae) that usually provide most of the energy required by the coral host. This relationship is sensitive to temperature stress; as little as a 1°C increase often leads to the collapse of the association. This sensitivity has led to an interest in the potential of more stress-tolerant algae to supplement or substitute for the normal Symbiodiniaceae mutualists. In this respect, the apicomplexan-like microalga Chromera is of particular interest due to its greater temperature tolerance. We generated a de novo transcriptome for a Chromera strain isolated from a GBR coral ('GBR Chromera') and compared with those of the reference strain of Chromera ('Sydney Chromera'), and to those of Symbiodiniaceae (Fugacium kawagutii, Cladocopium goreaui and Breviolum minutum), as well as the apicomplexan parasite, Plasmodium falciparum. In contrast to the high sequence divergence amongst representatives of different genera within the family Symbiodiniaceae, the two Chromera strains featured low sequence divergence at orthologous genes, implying that they are likely to be conspecifics. Although KEGG categories provide few criteria by which true coral mutualists might be identified, they do supply a molecular rationalization that explains the ecological dominance of Cladocopium spp. amongst Indo-Pacific reef corals. The presence of HSP20 genes may contribute to the high thermal tolerance of Chromera.

Sex in Symbiodiniaceae dinoflagellates: genomic evidence for independent loss of the canonical synaptonemal complex

Dinoflagellates of the Symbiodiniaceae family encompass diverse symbionts that are critical to corals and other species living in coral reefs. It is well known that sexual reproduction enhances adaptive evolution in changing environments. Although genes related to meiotic functions were reported in Symbiodiniaceae, cytological evidence of meiosis and fertilisation are however yet to be observed in these taxa. Using transcriptome and genome data from 21 Symbiodiniaceae isolates, we studied genes that encode proteins associated with distinct stages of meiosis and syngamy. We report the absence of genes that encode main components of the synaptonemal complex (SC), a protein structure that mediates homologous chromosomal pairing and class I crossovers. This result suggests an independent loss of canonical SCs in the alveolates, that also includes the SC-lacking ciliates. We hypothesise that this loss was due in part to permanently condensed chromosomes and repeat-rich sequences in Symbiodiniaceae (and other dinoflagellates) which favoured the SC-independent class II crossover pathway. Our results reveal novel insights into evolution of the meiotic molecular machinery in the ecologically important Symbiodiniaceae and in other eukaryotes.

Fine-scale structure among mesophotic populations of the great star coral Montastraea cavernosa revealed by SNP genotyping

Mesophotic reefs (30-150 m) have been proposed as potential refugia that facilitate the recovery of degraded shallow reefs following acute disturbances such as coral bleaching and disease. However, because of the technical difficulty of collecting samples, the connectivity of adjacent mesophotic reefs is relatively unknown compared with shallower counterparts. We used genotyping by sequencing to assess fine-scale genetic structure of Montastraea cavernosa at two sites at Pulley Ridge, a mesophotic coral reef ecosystem in the Gulf of Mexico, and downstream sites along the Florida Reef Tract. We found differentiation between reefs at Pulley Ridge (~68 m) and corals at downstream upper mesophotic depths in the Dry Tortugas (28-36 m) and shallow reefs in the northern Florida Keys (Key Biscayne, ~5 m). The spatial endpoints of our study were distinct, with the Dry Tortugas as a genetic intermediate. Most striking were differences in population structure among northern and southern sites at Pulley Ridge that were separated by just 12km. Unique patterns of clonality and outlier loci allele frequency support these sites as different populations and suggest that the long-distance horizontal connectivity typical of shallow-water corals may not be typical for mesophotic systems in Florida and the Gulf of Mexico. We hypothesize that this may be due to the spawning of buoyant gametes, which commits propagules to the surface, resulting in greater dispersal and lower connectivity than typically found between nearby shallow sites. Differences in population structure over small spatial scales suggest that demographic constraints and/or environmental disturbances may be more variable in space and time on mesophotic reefs compared with their shallow-water counterparts.

Adaptation to Bleaching: Are Thermotolerant Symbiodiniaceae Strains More Successful Than Other Strains Under Elevated Temperatures in a Model Symbiotic Cnidarian?

The ability of some symbiotic cnidarians to resist and better withstand stress factors that cause bleaching is a trait that is receiving increased attention. The adaptive bleaching hypothesis postulates that cnidarians that can form a stable symbiosis with thermotolerant Symbiodiniaceae strains may cope better with increasing seawater temperatures. We used polyps of the scyphozoan, Cassiopea xamachana, as a model system to test symbiosis success under heat stress. We sought to determine: (1) if aposymbiotic C. xamachana polyps could establish and maintain a symbiosis with both native and non-native strains of Symbiodiniaceae that all exhibit different tolerances to heat, (2) whether polyps with these newly acquired Symbiodiniaceae strains would strobilate (produce ephyra), and (3) if thermally tolerant Symbiodiniaceae strains that established and maintained a symbiosis exhibited greater success in response to heat stress (even if they are not naturally occurring in Cassiopea). Following recolonization of aposymbiotic C. xamachana polyps with different strains, we found that: (1) strains Smic, Stri, Slin, and Spil all established a stable symbiosis that promoted strobilation and (2) strains Bmin1 and Bmin2 did not establish a stable symbiosis and strobilation did not occur. Strains Smic, Stri, Slin, and Spil were used in a subsequent bleaching experiment; each of the strains was introduced to a subset of aposymbiotic polyps and once polyp tissues were saturated with symbionts they were subjected to elevated temperatures - 32°C and 34°C - for 2 weeks. Our findings indicate that, in general, pairings of polyps with Symbiodiniaceae strains that are native to Cassiopea (Stri and Smic) performed better than a non-native strain (Slin) even though this strain has a high thermotolerance. This suggests a degree of partner specificity that may limit the adaptive potential of certain cnidarians to increased ocean warming. We also observed that the free-living, non-native thermotolerant strain Spil was relatively successful in resisting bleaching during experimental trials. This suggests that free-living Symbiodiniaceae may provide a supply of potentially "new" thermotolerant strains to cnidarians following a bleaching event.

Shedding light: a phylotranscriptomic perspective illuminates the origin of photosymbiosis in marine bivalves

BACKGROUND

Photosymbiotic associations between metazoan hosts and photosynthetic dinoflagellates are crucial to the trophic and structural integrity of many marine ecosystems, including coral reefs. Although extensive efforts have been devoted to study the short-term ecological interactions between coral hosts and their symbionts, long-term evolutionary dynamics of photosymbiosis in many marine animals are not well understood. Within Bivalvia, the second largest class of mollusks, obligate photosymbiosis is found in two marine lineages: the giant clams (subfamily Tridacninae) and the heart cockles (subfamily Fraginae), both in the family Cardiidae. Morphologically, giant clams show relatively conservative shell forms whereas photosymbiotic fragines exhibit a diverse suite of anatomical adaptations including flattened shells, leafy mantle extensions, and lens-like microstructural structures. To date, the phylogenetic relationships between these two subfamilies remain poorly resolved, and it is unclear whether photosymbiosis in cardiids originated once or twice.

RESULTS

In this study, we establish a backbone phylogeny for Cardiidae utilizing RNASeq-based transcriptomic data from Tridacninae, Fraginae and other cardiids. A variety of phylogenomic approaches were used to infer the relationship between the two groups. Our analyses found conflicting gene signals and potential rapid divergence among the lineages. Overall, results support a sister group relationship between Tridacninae and Fraginae, which diverged during the Cretaceous. Although a sister group relationship is recovered, ancestral state reconstruction using maximum likelihood-based methods reveals two independent origins of photosymbiosis, one at the base of Tridacninae and the other within a symbiotic Fraginae clade.

CONCLUSIONS

The newly revealed common ancestry between Tridacninae and Fraginae brings a possibility that certain genetic, metabolic, and/or anatomical exaptations existed in their last common ancestor, which promoted both lineages to independently establish photosymbiosis, possibly in response to the modern expansion of reef habitats.

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

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

Gene regulation underpinning increased thermal tolerance in a laboratory-evolved coral photosymbiont

Small increases in ocean temperature can disrupt the obligate symbiosis between corals and dinoflagellate microalgae, resulting in coral bleaching. Little is known about the genes that drive the physiological and bleaching response of algal symbionts to elevated temperature. Moreover, many studies to-date have compared highly divergent strains, making it challenging to accredit specific genes to contrasting traits. Here, we compare transcriptional responses at ambient (27°C) and bleaching-relevant (31°C) temperatures in a monoclonal, wild-type (WT) strain of Symbiodiniaceae to those of a selected-strain (SS), derived from the same monoclonal culture and experimentally evolved to elevated temperature over 80 generations (2.5 years). Thousands of genes were differentially expressed at a log fold-change of >8 between the WT and SS over a 35 days temperature treatment period. At 31°C, WT cells exhibited a temporally unstable transcriptomic response upregulating genes involved in the universal stress response such as molecular chaperoning, protein repair, protein degradation and DNA repair. Comparatively, SS cells exhibited a temporally stable transcriptomic response and downregulated many stress response genes that were upregulated by the WT. Among the most highly upregulated genes in the SS at 31°C were algal transcription factors and a gene probably of bacterial origin that encodes a type II secretion system protein, suggesting interactions with bacteria may contribute to the increased thermal tolerance of the SS. Genes and functional pathways conferring thermal tolerance in the SS could be targeted in future genetic engineering experiments designed to develop thermally resilient algal symbionts for use in coral restoration and conservation.