Draft Assembly of the Symbiodinium minutum Nuclear Genome Reveals Dinoflagellate Gene Structure

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.

Three Novel Bacteria Associated with Two Centric Diatom Species from the Mediterranean Sea, Thalassiosira rotula and Skeletonema marinoi

Diatoms are a successful group of microalgae at the base of the marine food web. For hundreds of millions of years, they have shared common habitats with bacteria, which favored the onset of interactions at different levels, potentially driving the synthesis of biologically active molecules. To unveil their presence, we sequenced the genomes of bacteria associated with the centric diatom Thalassiosira rotula from the Gulf of Naples. Annotation of the metagenome and its analysis allowed the reconstruction of three bacterial genomes that belong to currently undescribed species. Their investigation showed the existence of novel gene clusters coding for new polyketide molecules, antibiotics, antibiotic-resistance genes and an ectoine production pathway. Real-time PCR was used to investigate the association of these bacteria with three different diatom clones and revealed their preference for T. rotula FE80 and Skeletonema marinoi FE7, but not S. marinoi FE60 from the North Adriatic Sea. Additionally, we demonstrate that although all three bacteria could be detected in the culture supernatant (free-living), their number is up to 45 times higher in the cell associated fraction, suggesting a close association between these bacteria and their host. We demonstrate that axenic cultures of T. rotula are unable to grow in medium with low salinity (<28 ppt NaCl) whereas xenic cultures can tolerate up to 40 ppt NaCl with concomitant ectoine production, likely by the associated bacteria.

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

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

Spatial organization of dinoflagellate genomes: Novel insights and remaining critical questions

As is true for many other aspects, genome architecture, evolution, and function in dinoflagellates are enigmatic and, in the meantime, continuous inspiration for scientific quests. Recent third-generation sequencing and Hi-C linkage analyses brought new insights into the spatial organization of symbiodiniacean genomes, revealing the topologically associated domains, discrete gene clusters and their cis and trans orientations, and relationships with transcription. Where do these new findings bring us in dinoflagellate genomics? Here, we aim to place these new results in the backdrop of the long history of research on this topic and in the context of what critical questions remain to be pursued in the future. The new data suggest, pending verification of other complete chromosome assemblies, a potential evolutionary trend in chromosome number decrease and length increase within the Symbiodiniaceae. While questions remain about the mechanics of the three-dimensional chromosome structure and cell cycle-related DNA replication, the mechanisms of gene transcription and genome size evolution, these latest findings set new starting points for further inquiries.

Dependence of genome size and copy number of rRNA gene on cell volume in dinoflagellates

Dinoflagellates are an ecologically important group of protists in aquatic environment and have evolved many unusual and enigmatic genomic features such as immense genome sizes, high repeated genes, and a large portion of hydroxymethyluracil in DNA. Although previous studies have observed positive correlations between the large subunit (LSU) rRNA gene copy number and genome size of a variety of eukaryotic organisms (e.g. higher plants and animals), or between cell volume and LSU rRNA gene copy number, and/or between genome size and cell size, which suggests a possible co-evolution among these three features in different lineages of life, it remains an open question regarding the relationships among these three parameters in dinoflagellates. For the first time, we estimated the copy numbers of the LSU rRNA gene, the genome sizes, and cell volumes within a broad range of dinoflagellates (covering 15 species of 11 genera) using single-cell qPCR-based assay (determining LSU rRNA gene copy number), FlowCAM (cell volume measurement), and ultraviolet spectrophotometry (genome size estimation). The measured copy number of LSU rRNA gene ranged from 398 ± 184 (Prorocentrum minimum) to 152,078 ± 33,555 copies•cell^-1 (Alexandrium pacificum), while the genome size and the cell volume ranged from 5.6 ± 0.2 (Karlodinium veneficum) to 853 ± 19.9 pg•cell^-1 (Pseliodinium pirum), and from 1,070 ± 225 (Kar. veneficum) to 168,474 ± 124,180 μm³ (Ps. pirum), respectively. Together with the three parameters measured in literature, there are significant positive linear correlations between LSU rRNA gene copy numbers and genome sizes, cell volumes and LSU rRNA gene copy numbers, and between genome sizes and cell volumes via comparisons of multi-model regression analyses, suggesting a dependence of genome size and rRNA gene copy number on the cell volumes of dinoflagellates. Validation of the measurement methods was conducted via comparisons between reported data in the literature and that predicted using the linear equations we obtained, and between genome size measured by flow cytometry (FCM) and ultraviolet spectrophotometry (Nanodrop). These results provide insightful understandings of dinoflagellate evolution in terms of the relationships among genomes, gene copy number, and cell volume, and of rRNA gene-based studies in intra-populational and intra-individual genetic diversity, taxonomy, and diversity assessment in the environment of dinoflagellates. The results also provide a dataset useful for reads calibration in environmental metabarcoding studies of dinoflagellates and selection of candidate species for whole genome sequencing.

A Rhizobium bacterium and its population dynamics under different culture conditions of its associated toxic dinoflagellate Gambierdiscus balechii

Rhizobium bacteria are known as symbionts of legumes for developing nodules on plant roots and fixing N₂ for the host plants but unknown for associations with dinoflagellates. Here, we detected, isolated, and characterized a Rhizobium species from the marine toxic dinoflagellate Gambierdiscus culture. Its 16S rRNA gene (rDNA) is 99% identical to that of Rhizobium rosettiformans, and the affiliation is supported by the phylogenetic placement of its cell wall hydrolase -encoding gene (cwh). Using quantitative PCR of 16S rDNA and cwh, we found that the abundance of this bacterium increased during the late exponential growth phase of Gambierdiscus and under nitrogen limitation, suggesting potential physiological interactions between the dinoflagellate and the bacterium. This is the first report of dinoflagellate-associated Rhizobium bacterium, and its prevalence and ecological roles in dinoflagellate-Rhizobium relationships remain to be investigated in the future.

Supplementary Information

The online version contains supplementary material available at 10.1007/s42995-021-00102-1.

Elucidating gene expression adaptation of phylogenetically divergent coral holobionts under heat stress

As coral reefs struggle to survive under climate change, it is crucial to know whether they have the capacity to withstand changing conditions, particularly increasing seawater temperatures. Thermal tolerance requires the integrative response of the different components of the coral holobiont (coral host, algal photosymbiont, and associated microbiome). Here, using a controlled thermal stress experiment across three divergent Caribbean coral species, we attempt to dissect holobiont member metatranscriptome responses from coral taxa with different sensitivities to heat stress and use phylogenetic ANOVA to study the evolution of gene expression adaptation. We show that coral response to heat stress is a complex trait derived from multiple interactions among holobiont members. We identify host and photosymbiont genes that exhibit lineage-specific expression level adaptation and uncover potential roles for bacterial associates in supplementing the metabolic needs of the coral-photosymbiont duo during heat stress. Our results stress the importance of integrative and comparative approaches across a wide range of species to better understand coral survival under the predicted rise in sea surface temperatures.

Genetic changes involving the coral gastrovascular system support the transition between colonies and bailed-out polyps: evidence from a Pocillopora acuta transcriptome

Background

A coral colony is composed of physiologically integrated polyps. In stony corals, coloniality adopts a wide diversity of forms and involves complex ontogenetic dynamics. However, molecular mechanisms underlying coloniality have been little studied. To understand the genetic basis of coloniality and its contribution to coral ecology, we induced polyp bail-out in a colonial coral, Pocillopora acuta, and compared transcription profiles of bailed-out polyps and polyps in normal colonies, and their responses to heat shock and hyposalinity.

Results

Consistent with morphological formation of a gastrovascular system and its neural transmission and molecular transport functions, we found genetic activation of neurogenesis and development of tube-like structures in normal colonies that is absent in bailed-out polyps. Moreover, relative to bailed-out polyps, colonies showed significant overexpression of genes for angiotensin-converting enzymes and endothelin-converting enzymes. In response to hyperthermal and hyposaline treatments, a high proportion of genetic regulation proved specific to either bailed-out polyps or colonies. Elevated temperatures even activated NF-κB signaling in colonies. On the other hand, colonies showed no discernible advantage over bailed-out polyps in regard to hyposalinity.

Conclusions

The present study provides a first look at the genetic basis of coloniality and documents different responses to environmental stimuli in P. acuta colonies versus those in bailed-out polyps. Overexpression of angiotensin-converting enzymes and endothelin-converting enzymes in colonies suggests possible involvement of these genes in development of the gastrovascular system in P. acuta. Functional characterization of these coral genes and further investigation of other forms of the transition to coloniality in stony corals should be fruitful areas for future research.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-08026-x.

Insights into Alexandrium minutum Nutrient Acquisition, Metabolism and Saxitoxin Biosynthesis through Comprehensive Transcriptome Survey

Simple Summary

Alexandrium minutum is one of the causing organisms for the occurrence of harmful algae bloom (HABs) in marine ecosystems. This species produces saxitoxin, one of the deadliest neurotoxins which can cause human mortality. However, molecular information such as genes and proteins catalog on this species is still lacking. Therefore, this study has successfully characterized several new molecular mechanisms regarding A. minutum environmental adaptation and saxitoxin biosynthesis. Ultimately, this study provides a valuable resource for facilitating future dinoflagellates’ molecular response to environmental changes.

Abstract

The toxin-producing dinoflagellate Alexandrium minutum is responsible for the outbreaks of harmful algae bloom (HABs). It is a widely distributed species and is responsible for producing paralytic shellfish poisoning toxins. However, the information associated with the environmental adaptation pathway and toxin biosynthesis in this species is still lacking. Therefore, this study focuses on the functional characterization of A. minutum unigenes obtained from transcriptome sequencing using the Illumina Hiseq 4000 sequencing platform. A total of 58,802 (47.05%) unigenes were successfully annotated using public databases such as NCBI-Nr, UniprotKB, EggNOG, KEGG, InterPRO and Gene Ontology (GO). This study has successfully identified key features that enable A. minutum to adapt to the marine environment, including several carbon metabolic pathways, assimilation of various sources of nitrogen and phosphorus. A. minutum was found to encode homologues for several proteins involved in saxitoxin biosynthesis, including the first three proteins in the pathway of saxitoxin biosynthesis, namely sxtA, sxtG and sxtB. The comprehensive transcriptome analysis presented in this study represents a valuable resource for understanding the dinoflagellates molecular metabolic model regarding nutrient acquisition and biosynthesis of saxitoxin.

Large-scale genome sequencing reveals the driving forces of viruses in microalgal evolution

Being integral primary producers in diverse ecosystems, microalgal genomes could be mined for ecological insights, but representative genome sequences are lacking for many phyla. We cultured and sequenced 107 microalgae species from 11 different phyla indigenous to varied geographies and climates. This collection was used to resolve genomic differences between saltwater and freshwater microalgae. Freshwater species showed domain-centric ontology enrichment for nuclear and nuclear membrane functions, while saltwater species were enriched in organellar and cellular membrane functions. Further, marine species contained significantly more viral families in their genomes (p = 8e-4). Sequences from Chlorovirus, Coccolithovirus, Pandoravirus, Marseillevirus, Tupanvirus, and other viruses were found integrated into the genomes of algal from marine environments. These viral-origin sequences were found to be expressed and code for a wide variety of functions. Together, this study comprehensively defines the expanse of protein-coding and viral elements in microalgal genomes and posits a unified adaptive strategy for algal halotolerance.

Effects of the herbicide Irgarol 1051 on the transcriptome of hermatypic coral Acropora tenuis and its symbiotic dinoflagellates

Coral reefs face multiple threats, including climate change, agricultural runoff, shipping activities, coastal development, and chemical pollutants. Irgarol 1051, a PSII herbicide, has been used as an antifouling booster since the previously used antibiofouling agent tributyltin (TBT) was banned worldwide. Although the mechanisms through which elevated temperatures cause coral bleaching have been reported, it remains unclear how PSII herbicides cause bleaching. Thus, in this study, we investigated the transcriptomes of Acropora tenuis and its symbiotic dinoflagellates by RNA-sequencing (RNA-Seq) to elucidate the molecular mechanisms underlying Irgarol-induced bleaching. Coral exposure to 10 μg/L Irgarol for 7 d affected coral body colour, specifically by an increase in their red, green, and blue (RGB) values; however, no such effect was observed in corals exposed to 1 μg/L Irgarol. RNA-Seq revealed the differentially expressed genes (DEGs) in corals and symbiotic dinoflagellates following Irgarol exposure. Coral DEGs encoded green fluorescent protein, blue-light-sensing photoreceptor (cryptochrome), chromoprotein, caspase 8, and nuclear receptors; DEGs in symbiotic dinoflagellates encoded light-harvesting proteins, photosystem II proteins, and heat shock proteins (i.e. HSP70 and HSP90), and ubiquitin. Bioinformatic analyses revealed that both Irgarol treatments disrupted various gene ontology terms, pathways, and protein interaction networks; these are different in corals (e.g. oxidative phosphorylation, metabolic pathway, transforming growth factor-β signalling pathway, adherens junction, and apoptosis) and symbiotic dinoflagellates (e.g. protein processing in endoplasmic reticulum, carbon fixation in photosynthetic organisms, metabolic pathway, and photosynthesis). Our data suggest that Irgarol disrupts the expression of various coral genes, thereby affecting various gene ontology terms, pathways, and protein interaction networks. Our study provides new insights into the potential molecular mechanisms underlying the bleaching effect of PSII herbicides, such as Irgarol, on corals and symbiotic dinoflagellates.

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.

Cation and Anion Channelrhodopsins: Sequence Motifs and Taxonomic Distribution

Cation and anion channelrhodopsins (CCRs and ACRs, respectively) primarily from two algal species, Chlamydomonas reinhardtii and Guillardia theta, have become widely used as optogenetic tools to control cell membrane potential with light. We mined algal and other protist polynucleotide sequencing projects and metagenomic samples to identify 75 channelrhodopsin homologs from four channelrhodopsin families, including one revealed in dinoflagellates in this study. We carried out electrophysiological analysis of 33 natural channelrhodopsin variants from different phylogenetic lineages and 10 metagenomic homologs in search of sequence determinants of ion selectivity, photocurrent desensitization, and spectral tuning in channelrhodopsins. Our results show that association of a reduced number of glutamates near the conductance path with anion selectivity depends on a wider protein context, because prasinophyte homologs with a glutamate pattern identical to that in cryptophyte ACRs are cation selective. Desensitization is also broadly context dependent, as in one branch of stramenopile ACRs and their metagenomic homologs, its extent roughly correlates with phylogenetic relationship of their sequences. Regarding spectral tuning, we identified two prasinophyte CCRs with red-shifted spectra to 585 nm. They exhibit a third residue pattern in their retinal-binding pockets distinctly different from those of the only two types of red-shifted channelrhodopsins known (i.e., the CCR Chrimson and RubyACRs). In cryptophyte ACRs we identified three specific residue positions in the retinal-binding pocket that define the wavelength of their spectral maxima. Lastly, we found that dinoflagellate rhodopsins with a TCP motif in the third transmembrane helix and a metagenomic homolog exhibit channel activity.

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

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

Tentacle Morphological Variation Coincides with Differential Expression of Toxins in Sea Anemones

Phylum Cnidaria is an ancient venomous group defined by the presence of cnidae, specialised organelles that serve as venom delivery systems. The distribution of cnidae across the body plan is linked to regionalisation of venom production, with tissue-specific venom composition observed in multiple actiniarian species. In this study, we assess whether morphological variants of tentacles are associated with distinct toxin expression profiles and investigate the functional significance of specialised tentacular structures. Using five sea anemone species, we analysed differential expression of toxin-like transcripts and found that expression levels differ significantly across tentacular structures when substantial morphological variation is present. Therefore, the differential expression of toxin genes is associated with morphological variation of tentacular structures in a tissue-specific manner. Furthermore, the unique toxin profile of spherical tentacular structures in families Aliciidae and Thalassianthidae indicate that vesicles and nematospheres may function to protect branched structures that host a large number of photosynthetic symbionts. Thus, hosting zooxanthellae may account for the tentacle-specific toxin expression profiles observed in the current study. Overall, specialised tentacular structures serve unique ecological roles and, in order to fulfil their functions, they possess distinct venom cocktails.

Environmental specialization and cryptic genetic divergence in two massive coral species from the Florida Keys Reef Tract

Broadcast-spawning coral species have wide geographical ranges spanning strong environmental gradients, but it is unclear how much spatially varying selection these gradients actually impose. Strong divergent selection might present a considerable barrier for demographic exchange between disparate reef habitats. We investigated whether the cross-shelf gradient is associated with spatially varying selection in two common coral species, Montastraea cavernosa and Siderastrea siderea, in the Florida Keys. To this end, we generated a de novo genome assembly for M. cavernosa and used 2bRAD to genotype 20 juveniles and 20 adults of both species from each of the three reef zones to identify signatures of selection occurring within a single generation. Unexpectedly, each species was found to be composed of four genetically distinct lineages, with gene flow between them still ongoing but highly reduced in 13.0%-54.7% of the genome. Each species includes two sympatric lineages that are only found in the deep (20 m) habitat, while the other lineages are found almost exclusively on the shallower reefs (3-10 m). The two "shallow" lineages of M. cavernosa are also specialized for either nearshore or offshore: comparison between adult and juvenile cohorts indicates that cross-shelf migrants are more than twice as likely to die before reaching adulthood than local recruits. S. siderea and M. cavernosa are among the most ecologically successful species on the Florida Keys Reef Tract, and this work offers important insight into the genomic background of divergent selection and environmental specialization that may in part explain their resilience and broad environmental range.

Genomic variation of an endosymbiotic dinoflagellate (Symbiodinium 'fitti') among closely related coral hosts

Mutualisms where hosts are coupled metabolically to their symbionts often exhibit high partner fidelity. Most reef-building coral species form obligate symbioses with a specific species of photosymbionts, dinoflagellates in the family Symbiodiniaceae, despite needing to acquire symbionts early in their development from environmental sources. Three Caribbean acroporids (Acropora palmata, A. cervicornis and their F₁ hybrid) are sympatric across much of their range, but often occupy different depth and light habitats. Throughout this range, both species and their hybrid associate with the endosymbiotic dinoflagellate Symbiodinium 'fitti'. Because light (and therefore depth) influences the physiology of dinoflagellates, we investigated whether S. 'fitti' populations from each host taxon were differentiated genetically. Single nucleotide polymorphisms (SNPs) among S. 'fitti' strains were identified by aligning shallow metagenomic sequences of acroporid colonies sampled from across the Caribbean to a ~600-Mb draft assembly of the S. 'fitti' genome (from the CFL14120 A. cervicornis metagenome). Phylogenomic and multivariate analyses revealed that genomic variation among S. 'fitti' strains partitioned to each host taxon rather than by biogeographical origin. This is particularly noteworthy because the hybrid has a sparse fossil record and may be of relatively recent origin. A subset (37.6%) of the SNPs putatively under selection were nonsynonymous mutations predicted to alter protein efficiency. Differences in genomic variation of S. 'fitti' strains from each host taxon may reflect the unique selection pressures created by the microenvironments associated with each host. The nonrandom sorting among S. 'fitti' strains to different hosts could be the basis for lineage diversification via disruptive selection, leading to ecological specialization and ultimately speciation.

Transcription-dependent domain-scale three-dimensional genome organization in the dinoflagellate Breviolum minutum

Dinoflagellate chromosomes represent a unique evolutionary experiment, as they exist in a permanently condensed, liquid crystalline state; are not packaged by histones; and contain genes organized into tandem gene arrays, with minimal transcriptional regulation. We analyze the three-dimensional genome of Breviolum minutum, and find large topological domains (dinoflagellate topologically associating domains, which we term 'dinoTADs') without chromatin loops, which are demarcated by convergent gene array boundaries. Transcriptional inhibition disrupts dinoTADs, implicating transcription-induced supercoiling as the primary topological force in dinoflagellates.

Micro-size plankton abundance and assemblages in the western North Pacific Subtropical Gyre under microscopic observation

While primary productivity in the oligotrophic North Pacific Subtropical Gyre (NPSG) is changing, the micro-size plankton community has not been evaluated in the last 4 decades, prompting a re-evaluation. We collected samples over three years (2016–2018) from depths of 10 to 200 m (n = 127), and the micro-size plankton were identified and counted to understand the heterogeneity of micro-size plankton community structure. The assemblages were consistent to the those of 4 decades ago. Dinophyceae (dinoflagellates) were the most numerically abundant, followed by Cryptophyceae and Bacillariophyceae (diatoms). The other micro-size plankton classes (Cyanophyceae, Haptophyceae, Dictyochophyceae, Euglenophyceae, and Prasinophyceae) were not always detected, whereas only Trichodesmium spp. was counted in the Cyanophyceae. Other unidentified autotrophic and heterotrophic flagellates were also significantly present, and their numeric abundance was higher than or at the same level as was that of the Dinophyceae. In the Dinophyceae, Gymnodiniaceae and Peridiniales were abundant. The chlorophyll a concentration and these class-level assemblages suggested micro-size plankton is not a major primary producer in this area. We applied generalized additive models (GAMs) and principal coordination analyses (PCoAs) to evaluate the habitats of every plankton group and the heterogeneity of the assemblages. The GAMs suggested that every classified plankton abundance showed a similar response to salinity, and we observed differences in habitats in terms of temperature and nitrate concentrations. Based on the PCoAs, we observed unique communities at the 200 m depth layer compared with those at the other sampling layers. The site scores of PCoAs indicated that the micro-size plankton assemblages are most heterogeneous at the 10 m depth layer. At such depth, diazotrophic Cyanophyceae (Trichodesmium spp.) are abundant, particularly in less-saline water. Therefore, nitrogen fixation may contribute to the heterogeneity in the abundance and assemblages in the western NPSG.

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.

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.

Photosynthesis acclimation under severely fluctuating light conditions allows faster growth of diatoms compared with dinoflagellates

BACKGROUND

Diatoms contribute 20% of the global primary production and are adaptable in dynamic environments. Diatoms always bloom earlier in the annual phytoplankton succession instead of dinoflagellates. However, how diatoms acclimate to a dynamic environment, especially under changing light conditions, remains unclear.

RESULTS

We compared the growth and photosynthesis under fluctuating light conditions of red tide diatom Skeletonema costatum, red tide dinoflagellate Amphidinium carterae, Prorocentrum donghaiense, Karenia mikimotoi, model diatom Phaeodactylum tricornutum, Thalassiosira pseudonana and model dinoflagellate Dinophycae Symbiodinium. Diatoms grew faster and maintained a consistently higher level of photosynthesis. Diatoms were sensitive to the specific inhibitor of Proton Gradient Regulation 5 (PGR5) depending photosynthetic electron flow, which is a crucial mechanism to protect their photosynthetic apparatus under fluctuating light. In contrast, the dinoflagellates were not sensitive to this inhibitor. Therefore, we investigate how PGR5 functions under light fluctuations in the model diatom P. tricornutum by knocking down and overexpressing PGR5. Overexpression of PGR5 reduced the photosystem I acceptor side limitation (Y (NA)) and increased growth rate under severely fluctuating light in contrast to the knockdown of PGR5.

CONCLUSION

Diatoms acclimatize to fluctuating light conditions better than dinoflagellates. PGR5 in diatoms can regulate their photosynthetic electron flow and accelerate their growth under severe light fluctuation, supporting fast biomass accumulation under dynamic environments in pioneer blooms.

Chiral Systems Made from DNA

The very chemical structure of DNA that enables biological heredity and evolution has non-trivial implications for the self-organization of DNA molecules into larger assemblies and provides limitless opportunities for building functional nanostructures. This progress report discusses the natural organization of DNA into chiral structures and recent advances in creating synthetic chiral systems using DNA as a building material. How nucleic acid chirality naturally comes into play in a diverse array of situations is considered first, at length scales ranging from an individual nucleotide to entire chromosomes. Thereafter, chiral liquid crystal phases formed by dense DNA mixtures are discussed, including the ongoing efforts to understand their origins. The report then summarizes recent efforts directed toward building chiral structures, and other structures of complex topology, using the principle of DNA self-assembly. Discussed last are existing and proposed functional man-made nanostructures designed to either probe or harness DNA's chirality, from plasmonics and spintronics to biosensing.

Conservation and turnover of miRNAs and their highly complementary targets in early branching animals

MicroRNAs (miRNAs) are crucial post-transcriptional regulators that have been extensively studied in Bilateria, a group comprising the majority of extant animals, where more than 30 conserved miRNA families have been identified. By contrast, bilaterian miRNA targets are largely not conserved. Cnidaria is the sister group to Bilateria and thus provides a unique opportunity for comparative studies. Strikingly, like their plant counterparts, cnidarian miRNAs have been shown to predominantly have highly complementary targets leading to transcript cleavage by Argonaute proteins. Here, we assess the conservation of miRNAs and their targets by small RNA sequencing followed by miRNA target prediction in eight species of Anthozoa (sea anemones and corals), the earliest-branching cnidarian class. We uncover dozens of novel miRNAs but only a few conserved ones. Further, given their high complementarity, we were able to computationally identify miRNA targets in each species. Besides evidence for conservation of specific miRNA target sites, which are maintained between sea anemones and stony corals across 500 Myr of evolution, we also find indications for convergent evolution of target regulation by different miRNAs. Our data indicate that cnidarians have only few conserved miRNAs and corresponding targets, despite their high complementarity, suggesting a high evolutionary turnover.

Correlation between Organelle Genetic Variation and RNA Editing in Dinoflagellates Associated with the Coral Acropora digitifera

In order to develop successful strategies for coral reef preservation, it is critical that the biology of both host corals and symbiotic algae are investigated. In the Ryukyu Archipelago, which encompasses many islands spread over ∼500 km of the Pacific Ocean, four major populations of the coral Acropora digitifera have been studied using whole-genome shotgun (WGS) sequence analysis (Shinzato C, Mungpakdee S, Arakaki N, Satoh N. 2015. Genome-wide single-nucleotide polymorphism (SNP) analysis explains coral diversity and recovery in the Ryukyu Archipelago. Sci Rep. 5:18211.). In contrast, the diversity of the symbiotic dinoflagellates associated with these A. digitifera populations is unknown. It is therefore unclear if these two core components of the coral holobiont share a common evolutionary history. This issue can be addressed for the symbiotic algal populations by studying the organelle genomes of their mitochondria and plastids. Here, we analyzed WGS data from ∼150 adult A. digitifera, and by mapping reads to the available reference genome sequences, we extracted 2,250 sequences representing 15 organelle genes of Symbiodiniaceae. Molecular phylogenetic analyses of these mitochondrial and plastid gene sets revealed that A. digitifera from the southern Yaeyama islands harbor a different Symbiodiniaceae population than the islands of Okinawa and Kerama in the north, indicating that the distribution of symbiont populations partially matches that of the four host populations. Interestingly, we found that numerous SNPs correspond to known RNA-edited sites in 14 of the Symbiodiniaceae organelle genes, with mitochondrial genes showing a stronger correspondence than plastid genes. These results suggest a possible correlation between RNA editing and SNPs in the two organelle genomes of symbiotic dinoflagellates.

Intracellular bacteria are common and taxonomically diverse in cultured and in hospite algal endosymbionts of coral reefs

Corals house a variety of microorganisms which they depend on for their survival, including endosymbiotic dinoflagellates (Symbiodiniaceae) and bacteria. While cnidarian–microorganism interactions are widely studied, Symbiodiniaceae–bacteria interactions are only just beginning to receive attention. Here, we describe the localization and composition of the bacterial communities associated with cultures of 11 Symbiodiniaceae strains from nine species and six genera. Three-dimensional confocal laser scanning and electron microscopy revealed bacteria are present inside the Symbiodiniaceae cells as well as closely associated with their external cell surface. Bacterial pure cultures and 16S rRNA gene metabarcoding from Symbiodiniaceae cultures highlighted distinct and highly diverse bacterial communities occur intracellularly, closely associated with the Symbiodiniaceae outer cell surface and loosely associated (i.e., in the surrounding culture media). The intracellular bacteria are highly conserved across Symbiodiniaceae species, suggesting they may be involved in Symbiodiniaceae physiology. Our findings provide unique new insights into the biology of Symbiodiniaceae.

Genomic adaptations to an endolithic lifestyle in the coral-associated alga Ostreobium

The green alga Ostreobium is an important coral holobiont member, playing key roles in skeletal decalcification and providing photosynthate to bleached corals that have lost their dinoflagellate endosymbionts. Ostreobium lives in the coral's skeleton, a low-light environment with variable pH and O₂ availability. We present the Ostreobium nuclear genome and a metatranscriptomic analysis of healthy and bleached corals to improve our understanding of Ostreobium's adaptations to its extreme environment and its roles as a coral holobiont member. The Ostreobium genome has 10,663 predicted protein-coding genes and shows adaptations for life in low and variable light conditions and other stressors in the endolithic environment. This alga presents a rich repertoire of light-harvesting complex proteins but lacks many genes for photoprotection and photoreceptors. It also has a large arsenal of genes for oxidative stress response. An expansion of extracellular peptidases suggests that Ostreobium may supplement its energy needs by feeding on the organic skeletal matrix, and a diverse set of fermentation pathways allows it to live in the anoxic skeleton at night. Ostreobium depends on other holobiont members for vitamin B12, and our metatranscriptomes identify potential bacterial sources. Metatranscriptomes showed Ostreobium becoming a dominant agent of photosynthesis in bleached corals and provided evidence for variable responses among coral samples and different Ostreobium genotypes. Our work provides a comprehensive understanding of the adaptations of Ostreobium to its extreme environment and an important genomic resource to improve our comprehension of coral holobiont resilience, bleaching, and recovery.

Shimiella gen. nov. and Shimiella gracilenta sp. nov. (Dinophyceae, Kareniaceae), a Kleptoplastidic Dinoflagellate from Korean Waters and its Survival under Starvation

A small dinoflagellate, ~13 μm in cell length, was isolated from Jinhae Bay, Korea. Light microscopy showed that it was similar to the kleptoplastidic dinoflagellate Gymnodinium gracilentum nom. inval. rDNA sequences were obtained and its anatomy and morphology described using light and scanning and transmission electron microscopy. Phylogenetic analyses indicated that it belonged to the family Kareniaceae. However, its large subunit (LSU) rDNA sequences were 5.2-9.5% different from those of the other five genera in the family, and its clade was clearly divergent from that of each genus. Its overall morphology was different from those of the other five genera in the family and from Gymnodinium. Unlike Gymnodinium, this dinoflagellate did not have a horseshoe-shaped apical groove, nuclear envelope chambers, or a nuclear fibrous connective (NFC). It had an apical line of narrow amphiesmal vesicles and an elongated apical furrow crossing the apex. Cells were covered with polygonal amphiesmal vesicles arranged in 16 rows. Starved cells did not contain their own plastids, eyespots, pyrenoids, peridinin, or fucoxanthin. However, they could survive without added prey for approximately one month using chloroplasts from the cryptophyte prey Teleaulax amphioxeia, indicating kleptoplastidy. Because this taxon is genetically distinct at the generic rank from the other genera in Kareniaceae, it is placed in Shimiella gen. nov., and because G. gracilentum was invalid, the new bionomial S. gracilenta sp. nov. is proposed.

Rapid protein evolution, organellar reductions, and invasive intronic elements in the marine aerobic parasite dinoflagellate Amoebophrya spp

BACKGROUND

Dinoflagellates are aquatic protists particularly widespread in the oceans worldwide. Some are responsible for toxic blooms while others live in symbiotic relationships, either as mutualistic symbionts in corals or as parasites infecting other protists and animals. Dinoflagellates harbor atypically large genomes (~ 3 to 250 Gb), with gene organization and gene expression patterns very different from closely related apicomplexan parasites. Here we sequenced and analyzed the genomes of two early-diverging and co-occurring parasitic dinoflagellate Amoebophrya strains, to shed light on the emergence of such atypical genomic features, dinoflagellate evolution, and host specialization.

RESULTS

We sequenced, assembled, and annotated high-quality genomes for two Amoebophrya strains (A25 and A120), using a combination of Illumina paired-end short-read and Oxford Nanopore Technology (ONT) MinION long-read sequencing approaches. We found a small number of transposable elements, along with short introns and intergenic regions, and a limited number of gene families, together contribute to the compactness of the Amoebophrya genomes, a feature potentially linked with parasitism. While the majority of Amoebophrya proteins (63.7% of A25 and 59.3% of A120) had no functional assignment, we found many orthologs shared with Dinophyceae. Our analyses revealed a strong tendency for genes encoded by unidirectional clusters and high levels of synteny conservation between the two genomes despite low interspecific protein sequence similarity, suggesting rapid protein evolution. Most strikingly, we identified a large portion of non-canonical introns, including repeated introns, displaying a broad variability of associated splicing motifs never observed among eukaryotes. Those introner elements appear to have the capacity to spread over their respective genomes in a manner similar to transposable elements. Finally, we confirmed the reduction of organelles observed in Amoebophrya spp., i.e., loss of the plastid, potential loss of a mitochondrial genome and functions.

CONCLUSION

These results expand the range of atypical genome features found in basal dinoflagellates and raise questions regarding speciation and the evolutionary mechanisms at play while parastitism was selected for in this particular unicellular lineage.

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.

Speciation-by-depth on coral reefs: Sympatric divergence with gene flow or cryptic transient isolation?

The distributions of many sister species in the sea overlap geographically but are partitioned along depth gradients. The genetic changes leading to depth segregation may evolve in geographic isolation as a prerequisite to coexistence or may emerge during primary divergence leading to new species. These alternatives can now be distinguished via the power endowed by the thousands of scorable loci provided by second-generation sequence data. Here, we revisit the case of two depth-segregated, genetically isolated ecotypes of the nominal Caribbean candelabrum coral Eunicea flexuosa. Previous analyses based on a handful of markers could not distinguish between models of genetic exchange after a period of isolation (consistent with secondary contact) and divergence with gene flow (consistent with primary divergence). Analyses of the history of isolation, genetic exchange and population size based on 15,640 new SNP markers derived from RNAseq data best support models where divergence began 800K BP and include epochs of divergence with gene flow, but with an intermediate period of transient isolation. Results also supported the previous conclusion that recent exchange between the ecotypes occurs asymmetrically from the Shallow lineage to the Deep. Parallel analyses of data from two other corals with depth-segregated populations (Agaricia fragilis and Pocillopora damicornis) suggest divergence leading to depth-segregated populations may begin with a period of symmetric exchange, but that an epoch of population isolation precedes more complete isolation marked by asymmetric introgression. Thus, while divergence-with-gene flow may account for much of the differentiation that separates closely related, depth-segregated species, it remains to be seen whether any critical steps in the speciation process only occur when populations are isolated.

Bioactivity Potential of Marine Natural Products from Scleractinia-Associated Microbes and In Silico Anti-SARS-COV-2 Evaluation

Marine organisms and their associated microbes are rich in diverse chemical leads. With the development of marine biotechnology, a considerable number of research activities are focused on marine bacteria and fungi-derived bioactive compounds. Marine bacteria and fungi are ranked on the top of the hierarchy of all organisms, as they are responsible for producing a wide range of bioactive secondary metabolites with possible pharmaceutical applications. Thus, they have the potential to provide future drugs against challenging diseases, such as cancer, a range of viral diseases, malaria, and inflammation. This review aims at describing the literature on secondary metabolites that have been obtained from Scleractinian-associated organisms including bacteria, fungi, and zooxanthellae, with full coverage of the period from 1982 to 2020, as well as illustrating their biological activities and structure activity relationship (SAR). Moreover, all these compounds were filtered based on ADME analysis to determine their physicochemical properties, and 15 compounds were selected. The selected compounds were virtually investigated for potential inhibition for SARS-CoV-2 targets using molecular docking studies. Promising potential results against SARS-CoV-2 RNA dependent RNA polymerase (RdRp) and methyltransferase (nsp16) are presented.

Phylogenetic analysis of cell-cycle regulatory proteins within the Symbiodiniaceae

In oligotrophic waters, cnidarian hosts rely on symbiosis with their photosynthetic dinoflagellate partners (family Symbiodiniaceae) to obtain the nutrients they need to grow, reproduce and survive. For this symbiosis to persist, the host must regulate the growth and proliferation of its symbionts. One of the proposed regulatory mechanisms is arrest of the symbiont cell cycle in the G₁ phase, though the cellular mechanisms involved remain unknown. Cell-cycle progression in eukaryotes is controlled by the conserved family of cyclin-dependent kinases (CDKs) and their partner cyclins. We identified CDKs and cyclins in different Symbiodiniaceae species and examined their relationship to homologs in other eukaryotes. Cyclin proteins related to eumetazoan cell-cycle-related cyclins A, B, D, G/I and Y, and transcriptional cyclin L, were identified in the Symbiodiniaceae, alongside several alveolate-specific cyclin A/B proteins, and proteins related to protist P/U-type cyclins and apicomplexan cyclins. The largest expansion of Symbiodiniaceae cyclins was in the P/U-type cyclin groups. Proteins related to eumetazoan cell-cycle-related CDKs (CDK1) were identified as well as transcription-related CDKs. The largest expansion of CDK groups was, however, in alveolate-specific groups which comprised 11 distinct CDK groups (CDKA-J) with CDKB being the most widely distributed CDK protein. As a result of its phylogenetic position, conservation across Symbiodiniaceae species, and the presence of the canonical CDK motif, CDKB emerged as a likely candidate for a Saccharomyces cerevisiae Cdc28/Pho85-like homolog in Symbiodiniaceae. Similar to cyclins, two CDK-groups found in Symbiodiniaceae species were solely associated with apicomplexan taxa. A comparison of Breviolum minutum CDK and cyclin gene expression between free-living and symbiotic states showed that several alveolate-specific CDKs and two P/U-type cyclins exhibited altered expression in hospite, suggesting that symbiosis influences the cell cycle of symbionts on a molecular level. These results highlight the divergence of Symbiodiniaceae cell-cycle proteins across species. These results have important implications for host control of the symbiont cell cycle in novel cnidarian–dinoflagellate symbioses.

Transcriptome Reprogramming of Symbiodiniaceae Breviolum minutum in Response to Casein Amino Acids Supplementation

Dinoflagellates in the family Symbiodiniaceae can live freely in ocean waters or form a symbiosis with a variety of cnidarians including corals, sea anemones, and jellyfish. Trophic plasticity of Symbiodiniaceae is critical to its ecological success as it moves between environments. However, the molecular mechanisms underlying these trophic shifts in Symbiodiniaceae are still largely unknown. Using Breviolum minutum strain SSB01 (designated SSB01) as a model, we showed that Symbiodiniaceae go through a physiological and transcriptome reprogramming when the alga is grown with the organic nitrogen containing nutrients in hydrolyzed casein, but not with inorganic nutrients. SSB01 grows at a much faster rate and maintains stable photosynthetic efficiency when supplemented with casein amino acids compared to only inorganic nutrients or seawater. These physiological changes are driven by massive transcriptome changes in SSB01 supplemented with casein amino acids. The levels of transcripts encoding proteins involved in altering DNA conformation such as DNA topoisomerases, histones, and chromosome structural components were all significantly changed. Functional enrichment analysis also revealed processes involved in translation, ion transport, generation of second messengers, and phosphorylation. The physiological and molecular changes that underlie in vitro trophic transitions in Symbiodiniaceae can serve as an orthogonal platform to further understand the factors that impact the Symbiodiniaceae lifestyle.

Whole-Genome Transcriptome Analyses of Native Symbionts Reveal Host Coral Genomic Novelties for Establishing Coral-Algae Symbioses

Reef-building corals and photosynthetic, endosymbiotic algae of the family Symbiodiniaceae establish mutualistic relationships that are fundamental to coral biology, enabling coral reefs to support a vast diversity of marine species. Although numerous types of Symbiodiniaceae occur in coral reef environments, Acropora corals select specific types in early life stages. In order to study molecular mechanisms of coral–algal symbioses occurring in nature, we performed whole-genome transcriptomic analyses of Acropora tenuis larvae inoculated with Symbiodinium microadriaticum strains isolated from an Acropora recruit. In order to identify genes specifically involved in symbioses with native symbionts in early life stages, we also investigated transcriptomic responses of Acropora larvae exposed to closely related, nonsymbiotic, and occasionally symbiotic Symbiodinium strains. We found that the number of differentially expressed genes was largest when larvae acquired native symbionts. Repertoires of differentially expressed genes indicated that corals reduced amino acid, sugar, and lipid metabolism, such that metabolic enzymes performing these functions were derived primarily from S. microadriaticum rather than from A. tenuis. Upregulated gene expression of transporters for those metabolites occurred only when coral larvae acquired their natural symbionts, suggesting active utilization of native symbionts by host corals. We also discovered that in Acropora, genes for sugar and amino acid transporters, prosaposin-like, and Notch ligand-like, were upregulated only in response to native symbionts, and included tandemly duplicated genes. Gene duplications in coral genomes may have been essential to establish genomic novelties for coral–algae symbiosis.

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

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

Different levels of energetic coupling between photosynthesis and respiration do not determine the occurrence of adaptive responses of Symbiodiniaceae to global warming

Disentangling the metabolic functioning of corals' endosymbionts (Symbiodiniaceae) is relevant to understanding the response of coral reefs to warming oceans. In this work, we first question whether there is an energetic coupling between photosynthesis and respiration in Symbiodiniaceae (Symbiodinium, Durusdinium and Effrenium), and second, how different levels of energetic coupling will affect their adaptive responses to global warming. Coupling between photosynthesis and respiration was established by determining the variation of metabolic rates during thermal response curves, and how inhibition of respiration affects photosynthesis. Adaptive (irreversible) responses were studied by exposing two Symbiodinium species with different levels of photosynthesis-respiration interaction to high temperature conditions (32°C) for 1 yr. We found that some Symbiodiniaceae have a high level of energetic coupling; that is, photosynthesis and respiration have the same temperature dependency, and photosynthesis is negatively affected when respiration is inhibited. Conversely, photosynthesis and respiration are not coupled in other species. In any case, prolonged exposure to high temperature caused adjustments in both photosynthesis and respiration, but these changes were fully reversible. We conclude that energetic coupling between photosynthesis and respiration exhibits wide variation amongst Symbiodiniaceae and does not determine the occurrence of adaptive responses in Symbiodiniaceae to temperature increase.

De novo transcriptome assembly from the gonads of a scleractinian coral, Euphyllia ancora: molecular mechanisms underlying scleractinian gametogenesis

BACKGROUND

Sexual reproduction of scleractinians has captured the attention of researchers and the general public for decades. Although extensive ecological data has been acquired, underlying molecular and cellular mechanisms remain largely unknown. In this study, to better understand mechanisms underlying gametogenesis, we isolated ovaries and testes at different developmental phases from a gonochoric coral, Euphyllia ancora, and adopted a transcriptomic approach to reveal sex- and phase-specific gene expression profiles. In particular, we explored genes associated with oocyte development and maturation, spermiogenesis, sperm motility / capacitation, and fertilization.

RESULTS

1.6 billion raw reads were obtained from 24 gonadal samples. De novo assembly of trimmed reads, and elimination of contigs derived from symbiotic dinoflagellates (Symbiodiniaceae) and other organisms yielded a reference E. ancora gonadal transcriptome of 35,802 contigs. Analysis of 4 developmental phases identified 2023 genes that were differentially expressed during oogenesis and 678 during spermatogenesis. In premature/mature ovaries, 631 genes were specifically upregulated, with 538 in mature testes. Upregulated genes included those involved in gametogenesis, gamete maturation, sperm motility / capacitation, and fertilization in other metazoans, including humans. Meanwhile, a large number of genes without homology to sequences in the SWISS-PROT database were also observed among upregulated genes in premature / mature ovaries and mature testes.

CONCLUSIONS

Our findings show that scleractinian gametogenesis shares many molecular characteristics with that of other metazoans, but it also possesses unique characteristics developed during cnidarian and/or scleractinian evolution. To the best of our knowledge, this study is the first to create a gonadal transcriptome assembly from any scleractinian. This study and associated datasets provide a foundation for future studies regarding gametogenesis and differences between male and female colonies from molecular and cellular perspectives. Furthermore, our transcriptome assembly will be a useful reference for future development of sex-specific and/or stage-specific germ cell markers that can be used in coral aquaculture and ecological studies.

Integrated omics unveil the secondary metabolic landscape of a basal dinoflagellate

BACKGROUND

Some dinoflagellates cause harmful algal blooms, releasing toxic secondary metabolites, to the detriment of marine ecosystems and human health. Our understanding of dinoflagellate toxin biosynthesis has been hampered by their unusually large genomes. To overcome this challenge, for the first time, we sequenced the genome, microRNAs, and mRNA isoforms of a basal dinoflagellate, Amphidinium gibbosum, and employed an integrated omics approach to understand its secondary metabolite biosynthesis.

RESULTS

We assembled the ~ 6.4-Gb A. gibbosum genome, and by probing decoded dinoflagellate genomes and transcriptomes, we identified the non-ribosomal peptide synthetase adenylation domain as essential for generation of specialized metabolites. Upon starving the cells of phosphate and nitrogen, we observed pronounced shifts in metabolite biosynthesis, suggestive of post-transcriptional regulation by microRNAs. Using Iso-Seq and RNA-seq data, we found that alternative splicing and polycistronic expression generate different transcripts for secondary metabolism.

CONCLUSIONS

Our genomic findings suggest intricate integration of various metabolic enzymes that function iteratively to synthesize metabolites, providing mechanistic insights into how dinoflagellates synthesize secondary metabolites, depending upon nutrient availability. This study provides insights into toxin production associated with dinoflagellate blooms. The genome of this basal dinoflagellate provides important clues about dinoflagellate evolution and overcomes the large genome size, which has been a challenge previously.

Clade-Specific Sterol Metabolites in Dinoflagellate Endosymbionts Are Associated with Coral Bleaching in Response to Environmental Cues

Cnidarians cannot synthesize sterols (which play essential roles in growth and development) de novo but often use sterols acquired from endosymbiotic dinoflagellates. While sterol availability can impact the mutualistic interaction between coral host and algal symbiont, the biosynthetic pathways (in the dinoflagellate endosymbionts) and functional roles of sterols in these symbioses are poorly understood. In this study, we found that itraconazole, which perturbs sterol metabolism by inhibiting the sterol 14-demethylase CYP51 in dinoflagellates, induces bleaching of the anemone Heteractis crispa and that bleaching perturbs sterol metabolism of the dinoflagellate. While Symbiodiniaceae have clade-specific sterol metabolites, they share features of the common sterol biosynthetic pathway but with distinct architecture and substrate specificity features of participating enzymes. Tracking sterol profiles and transcripts of enzymes involved in sterol biosynthesis across time in response to different environmental cues revealed similarities and idiosyncratic features of sterol synthesis in the endosymbiont Breviolum minutum Exposure of algal cultures to high levels of light, heat, and acidification led to alterations in sterol synthesis, including blocks through downregulation of squalene synthase transcript levels accompanied by marked growth reductions.IMPORTANCE These results indicate that sterol metabolites in Symbiodiniaceae are clade specific, that their biosynthetic pathways share architectural and substrate specificity features with those of animals and plants, and that environmental stress-specific perturbation of sterol biosynthesis in dinoflagellates can impair a key mutualistic partnership for healthy reefs.

Dual RNA-sequencing analyses of a coral and its native symbiont during the establishment of symbiosis

Despite the ecological significance of the mutualistic relationship between Symbiodiniaceae and reef-building corals, the molecular interactions during establishment of this relationship are not well understood. This is particularly true of the transcriptional changes that occur in the symbiont. In the current study, a dual RNA-sequencing approach was used to better understand transcriptional changes on both sides of the coral-symbiont interaction during the colonization of Acropora tenuis by a compatible Symbiodiniaceae strain (Cladocopium goreaui; ITS2 type C1). Comparison of transcript levels of the in hospite symbiont 3, 12, 48 and 72 hr after exposure to those of the same strain in culture revealed that extensive and generalized down-regulation of symbiont gene expression occurred during the infection process. Included in this "symbiosis-derived transcriptional repression" were a range of stress response and immune-related genes. In contrast, a suite of symbiont genes implicated in metabolism was upregulated in the symbiotic state. The coral data support the hypothesis that immune-suppression and arrest of phagosome maturation play important roles during the establishment of compatible symbioses, and additionally imply the involvement of some SCRiP family members in the colonization process. Consistent with previous ecological studies, the transcriptomic data suggest that active translocation of metabolites to the host may begin early in the colonization process, and thus that the mutualistic relationship can be established at the larval stage. This dual RNA-sequencing study provides insights into the transcriptomic remodelling that occurs in C. goreaui during transition to a symbiotic lifestyle and the novel coral genes implicated in symbiosis.

The 5S rRNA genes in Alexandrium: their use as a FISH chromosomal marker in studies of the diversity, cell cycle and sexuality of dinoflagellates

Chromosomal markers of the diversity and evolution of dinoflagellates are scarce because the genomes of these organisms are unique among eukaryotes in terms of their base composition and chromosomal structure. Similarly, a lack of appropriate tools has hindered studies of the chromosomal localization of 5S ribosomal DNA (rDNA) in the nucleosome-less chromosomes of dinoflagellates. In this study, we isolated and cloned 5S rDNA sequences from various toxin-producing species of the genus Alexandrium and developed a fluorescence in situ hybridization (FISH) probe that allows their chromosomal localization. Our results can be summarized as follows: 1) The 5S rDNA unit is composed of a highly conserved 122-bp coding region and an intergenic spacer (IGS), the length and sequence of which are variable even within strains. 2) Three different IGS types, one containing the U6 small nuclear RNA (snRNA) gene, were found among four of the studied species (A. minutum, A. tamarense, A. catenella and A. pacificum). 3) In all strains investigated by FISH (A. minutum, A. tamarense, A. pacificum, A. catenella, A. andersonii and A. ostenfeldii), 5S rDNA gene arrays were separate from the nucleolar organizer region, which contains the genes for the large 45S pre-ribosomal RNA. 4) One to three 5S rDNA sites per haploid genome were detected, depending on the strains/species. Intraspecific variability in the number of 5S rDNA sites was determined among strains of A. minutum and A. pacificum. 5) 5S rDNA is a useful chromosomal marker of mitosis progression and can be employed to differentiate vegetative (haploid) vs. planozygotes (diploid) cells. Thus, the FISH probe (oligo-Dino5Smix5) developed in this study facilitates analyses of the diversity, cell cycle and life stages of the genus Alexandrium.

Microbial community and transcriptional responses to increased temperatures in coral Pocillopora damicornis holobiont

A few studies have holistically examined successive changes in coral holobionts in response to increased temperatures. Here, responses of the coral host Pocillopora damicornis, its Symbiodiniaceae symbionts, and associated bacteria to increased water temperatures were investigated. High temperatures induced bleaching, but no coral mortality was observed. Transcriptome analyses showed that P. damicornis responded more quickly to elevated temperatures than its algal symbionts. Numerous genes putatively associated with apoptosis, exocytosis, and autophagy were upregulated in P. damicornis, suggesting that Symbiodiniaceae can be eliminated or expelled through these mechanisms when P. damicornis experiences heat stress. Furthermore, apoptosis in P. damicornis is presumably induced through tumour necrosis factor and p53 signalling and caspase pathways. The relative abundances of several coral disease-associated bacteria increased at 32°C, which may affect immune responses in heat-stressed corals and potentially accelerates the loss of algal symbionts. Additionally, consistency of Symbiodiniaceae community structures under heat stress suggests non-selective loss of Symbiodiniaceae. We propose that heat stress elicits interrelated response mechanisms in all parts of the coral holobiont.

Population genetics of the coral Acropora millepora: Toward genomic prediction of bleaching

Although reef-building corals are declining worldwide, responses to bleaching vary within and across species and are partly heritable. Toward predicting bleaching response from genomic data, we generated a chromosome-scale genome assembly for the coral Acropora millepora We obtained whole-genome sequences for 237 phenotyped samples collected at 12 reefs along the Great Barrier Reef, among which we inferred little population structure. Scanning the genome for evidence of local adaptation, we detected signatures of long-term balancing selection in the heat-shock co-chaperone sacsin We conducted a genome-wide association study of visual bleaching score for 213 samples, incorporating the polygenic score derived from it into a predictive model for bleaching in the wild. These results set the stage for genomics-based approaches in conservation strategies.

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.

N-Linked Surface Glycan Biosynthesis, Composition, Inhibition, and Function in Cnidarian-Dinoflagellate Symbiosis

The success of symbioses between cnidarian hosts (e.g., corals and sea anemones) and micro-algal symbionts hinges on the molecular interactions that govern the establishment and maintenance of intracellular mutualisms. As a fundamental component of innate immunity, glycan-lectin interactions impact the onset of marine endosymbioses, but our understanding of the effects of cell surface glycome composition on symbiosis establishment remains limited. In this study, we examined the canonical N-glycan biosynthesis pathway in the genome of the dinoflagellate symbiont Breviolum minutum (family Symbiodiniaceae) and found it to be conserved with the exception of the transferase GlcNAc-TII (MGAT2). Using coupled liquid chromatography-mass spectrometry (LC-MS/MS), we characterized the cell surface N-glycan content of B. minutum, providing the first insight into the molecular composition of surface glycans in dinoflagellates. We then used the biosynthesis inhibitors kifunensine and swainsonine to alter the glycan composition of B. minutum. Successful high-mannose enrichment via kifunensine treatment resulted in a significant decrease in colonization of the model sea anemone Aiptasia (Exaiptasia pallida) by B. minutum. Hybrid glycan enrichment via swainsonine treatment, however, could not be confirmed and did not impact colonization. We conclude that functional Golgi processing of N-glycans is critical for maintaining appropriate cell surface glycan composition and for ensuring colonization success by B. minutum.

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.

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

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