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

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.

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.

Multi-omics analysis reveals the molecular response to heat stress in a "red tide" dinoflagellate

BACKGROUND

"Red tides" are harmful algal blooms caused by dinoflagellate microalgae that accumulate toxins lethal to other organisms, including humans via consumption of contaminated seafood. These algal blooms are driven by a combination of environmental factors including nutrient enrichment, particularly in warm waters, and are increasingly frequent. The molecular, regulatory, and evolutionary mechanisms that underlie the heat stress response in these harmful bloom-forming algal species remain little understood, due in part to the limited genomic resources from dinoflagellates, complicated by the large sizes of genomes, exhibiting features atypical of eukaryotes.

RESULTS

We present the de novo assembled genome (~ 4.75 Gbp with 85,849 protein-coding genes), transcriptome, proteome, and metabolome from Prorocentrum cordatum, a globally abundant, bloom-forming dinoflagellate. Using axenic algal cultures, we study the molecular mechanisms that underpin the algal response to heat stress, which is relevant to current ocean warming trends. We present the first evidence of a complementary interplay between RNA editing and exon usage that regulates the expression and functional diversity of biomolecules, reflected by reduction in photosynthesis, central metabolism, and protein synthesis. These results reveal genomic signatures and post-transcriptional regulation for the first time in a pelagic dinoflagellate.

CONCLUSIONS

Our multi-omics analyses uncover the molecular response to heat stress in an important bloom-forming algal species, which is driven by complex gene structures in a large, high-G+C genome, combined with multi-level transcriptional regulation. The dynamics and interplay of molecular regulatory mechanisms may explain in part how dinoflagellates diversified to become some of the most ecologically successful organisms on Earth.

Genes possibly related to symbiosis in early life stages of Acropora tenuis inoculated with Symbiodinium microadriaticum

Due to the ecological importance of mutualism between reef-building corals and symbiotic algae (Family Symbiodiniaceae), various transcriptomic studies on coral-algal symbiosis have been performed; however, molecular mechanisms, especially genes essential to initiate and maintain these symbioses remain unknown. We investigated transcriptomic responses of Acropora tenuis to inoculation with the native algal symbiont, Symbiodinium microadriaticum, during early life stages, and identified possible symbiosis-related genes. Genes involved in immune regulation, protection against oxidative stress, and metabolic interactions between partners are particularly important for symbiosis during Acropora early life stages. In addition, molecular phylogenetic analysis revealed that some possible symbiosis-related genes originated by gene duplication in the Acropora lineage, suggesting that gene duplication may have been the driving force to establish stable mutualism in Acropora, and that symbiotic molecular mechanisms may vary among coral lineages.

Stony coral tissue loss disease intervention with amoxicillin leads to a reversal of disease-modulated gene expression pathways

Stony coral tissue loss disease (SCTLD) remains an unprecedented disease outbreak due to its high mortality rate and rapid spread throughout Florida's Coral Reef and wider Caribbean. A collaborative effort is underway to evaluate strategies that mitigate the spread of SCTLD across coral colonies and reefs, including restoration of disease-resistant genotypes, genetic rescue, and disease intervention with therapeutics. We conducted an in-situ experiment in Southeast Florida to assess molecular responses among SCTLD-affected Montastraea cavernosa pre- and post-application of the most widely used intervention method, CoreRx Base 2B with amoxicillin. Through Tag-Seq gene expression profiling of apparently healthy, diseased, and treated corals, we identified modulation of metabolomic and immune gene pathways following antibiotic treatment. In a complementary ex-situ disease challenge experiment, we exposed nursery-cultured M. cavernosa and Orbicella faveolata fragments to SCTLD-affected donor corals to compare transcriptomic profiles among clonal individuals from unexposed controls, those exposed and displaying disease signs, and corals exposed and not displaying disease signs. Suppression of metabolic functional groups and activation of stress gene pathways as a result of SCTLD exposure were apparent in both species. Amoxicillin treatment led to a 'reversal' of the majority of gene pathways implicated in disease response, suggesting potential recovery of corals following antibiotic application. In addition to increasing our understanding of molecular responses to SCTLD, we provide resource managers with transcriptomic evidence that disease intervention with antibiotics appears to be successful and may help to modulate coral immune responses to SCTLD. These results contribute to feasibility assessments of intervention efforts following disease outbreaks and improved predictions of coral reef health across the wider Caribbean.

Coral-algal endosymbiosis characterized using RNAi and single-cell RNA-seq

Corals form an endosymbiotic relationship with the dinoflagellate algae Symbiodiniaceae, but ocean warming can trigger algal loss, coral bleaching and death, and the degradation of ecosystems. Mitigation of coral death requires a mechanistic understanding of coral-algal endosymbiosis. Here we report an RNA interference (RNAi) method and its application to study genes involved in early steps of endosymbiosis in the soft coral Xenia sp. We show that a host endosymbiotic cell marker called LePin (lectin and kazal protease inhibitor domains) is a secreted Xenia lectin that binds to algae to initiate phagocytosis of the algae and coral immune response modulation. The evolutionary conservation of domains in LePin among marine anthozoans performing endosymbiosis suggests a general role in coral-algal recognition. Our work sheds light on the phagocytic machinery and posits a mechanism for symbiosome formation, helping in efforts to understand and preserve coral-algal relationships in the face of climate change.

The coral microbiome: towards an understanding of the molecular mechanisms of coral-microbiota interactions

Corals live in a complex, multipartite symbiosis with diverse microbes across kingdoms, some of which are implicated in vital functions, such as those related to resilience against climate change. However, knowledge gaps and technical challenges limit our understanding of the nature and functional significance of complex symbiotic relationships within corals. Here, we provide an overview of the complexity of the coral microbiome focusing on taxonomic diversity and functions of well-studied and cryptic microbes. Mining the coral literature indicate that while corals collectively harbour a third of all marine bacterial phyla, known bacterial symbionts and antagonists of corals represent a minute fraction of this diversity and that these taxa cluster into select genera, suggesting selective evolutionary mechanisms enabled these bacteria to gain a niche within the holobiont. Recent advances in coral microbiome research aimed at leveraging microbiome manipulation to increase coral's fitness to help mitigate heat stress-related mortality are discussed. Then, insights into the potential mechanisms through which microbiota can communicate with and modify host responses are examined by describing known recognition patterns, potential microbially derived coral epigenome effector proteins and coral gene regulation. Finally, the power of omics tools used to study corals are highlighted with emphasis on an integrated host-microbiota multiomics framework to understand the underlying mechanisms during symbiosis and climate change-driven dysbiosis.

Rapid shifts in thermal reaction norms and tolerance of brooded coral larvae following parental heat acclimation

Thermal priming of reef corals can enhance their heat tolerance; however, the legacy effects of heat stress during parental brooding on larval resilience remain understudied. This study investigated whether preconditioning adult coral Pocillopora damicornis to high temperatures (29°C and 32°C) could better prepare their larvae for heat stress. Results showed that heat-acclimated adults brooded larvae with reduced symbiont density and shifted thermal performance curves. Reciprocal transplant experiments demonstrated higher bleaching resistance and better photosynthetic and autotrophic performance in heat-exposed larvae from acclimated adults compared to unacclimated adults. RNA-seq revealed strong cellular stress responses in larvae from heat-acclimated adults that could have been effective in rescuing host cells from stress, as evidenced by the widespread upregulation of genes involved in cell cycle and mitosis. For symbionts, a molecular coordination between light harvesting, photoprotection and carbon fixation was detected in larvae from heat-acclimated adults, which may help optimize photosynthetic activity and yield under high temperature. Furthermore, heat acclimation led to opposing regulations of symbiont catabolic and anabolic pathways and favoured nutrient translocation to the host and thus a functional symbiosis. Notwithstanding, the improved heat tolerance was paralleled by reduced light-enhanced dark respiration, indicating metabolic depression for energy saving. Our findings suggest that adult heat acclimation can rapidly shift thermal tolerance of brooded coral larvae and provide integrated physiological and molecular evidence for this adaptive plasticity, which could increase climate resilience. However, the metabolic depression may be maladaptive for long-term organismal performance, highlighting the importance of curbing carbon emissions to better protect corals.

Coral larvae suppress heat stress response during the onset of symbiosis decreasing their odds of survival

The endosymbiosis between most corals and their photosynthetic dinoflagellate partners begins early in the host life history, when corals are larvae or juvenile polyps. The capacity of coral larvae to buffer climate-induced stress while in the process of symbiont acquisition could come with physiological trade-offs that alter behaviour, development, settlement and survivorship. Here we examined the joint effects of thermal stress and symbiosis onset on colonization dynamics, survival, metamorphosis and host gene expression of Acropora digitifera larvae. We found that thermal stress decreased symbiont colonization of hosts by 50% and symbiont density by 98.5% over 2 weeks. Temperature and colonization also influenced larval survival and metamorphosis in an additive manner, where colonized larvae fared worse or prematurely metamorphosed more often than noncolonized larvae under thermal stress. Transcriptomic responses to colonization and thermal stress treatments were largely independent, while the interaction of these treatments revealed contrasting expression profiles of genes that function in the stress response, immunity, inflammation and cell cycle regulation. The combined treatment either cancelled or lowered the magnitude of expression of heat-stress responsive genes in the presence of symbionts, revealing a physiological cost to acquiring symbionts at the larval stage with elevated temperatures. In addition, host immune suppression, a hallmark of symbiosis onset under ambient temperature, turned to immune activation under heat stress. Thus, by integrating the physical environment and biotic pressures that mediate presettlement event in corals, our results suggest that colonization may hinder larval survival and recruitment under projected climate scenarios.

Larval transcriptomic responses of a stony coral, Acropora tenuis, during initial contact with the native symbiont, Symbiodinium microadriaticum

Although numerous dinoflagellate species (Family Symbiodiniaceae) are present in coral reef environments, Acropora corals tend to select a single species, Symbiodinium microadriaticum, in early life stages, even though this species is rarely found in mature colonies. In order to identify molecular mechanisms involved in initial contact with native symbionts, we analyzed transcriptomic responses of Acropora tenuis larvae at 1, 3, 6, 12, and 24 h after their first contact with S. microadriaticum, as well as with non-native symbionts, including the non-symbiotic S. natans and the occasional symbiont, S. tridacnidorum. Some gene expression changes were detected in larvae inoculated with non-native symbionts at 1 h post-inoculation, but those returned to baseline levels afterward. In contrast, when larvae were exposed to native symbionts, we found that the number of differentially expressed genes gradually increased in relation to inoculation time. As a specific response to native symbionts, upregulation of pattern recognition receptor-like and transporter genes, and suppression of cellular function genes related to immunity and apoptosis, were exclusively observed. These findings indicate that coral larvae recognize differences between symbionts, and when the appropriate symbionts infect, they coordinate gene expression to establish stable mutualism.

Hypoxia as a physiological cue and pathological stress for coral larvae

Ocean deoxygenation events are intensifying worldwide and can rapidly drive adult corals into a state of metabolic crisis and bleaching-induced mortality, but whether coral larvae are subject to similar stress remains untested. We experimentally exposed apo-symbiotic coral larvae of Acropora selago to deoxygenation stress with subsequent reoxygenation aligned to their night-day light cycle, and followed their gene expression using RNA-Seq. After 12 h of deoxygenation stress (~2 mg O₂ /L), coral planulae demonstrated a low expression of HIF-targeted hypoxia response genes concomitant with a significantly high expression of PHD2 (a promoter of HIFα proteasomal degradation), similar to corresponding adult corals. Despite exhibiting a consistent swimming phenotype compared to control samples, the differential gene expression observed in planulae exposed to deoxygenation-reoxygenation suggests a disruption of pathways involved in developmental regulation, mitochondrial activity, lipid metabolism, and O₂ -sensitive epigenetic regulators. Importantly, we found that treated larvae exhibited a disruption in the expression of conserved HIF-targeted developmental regulators, for example, Homeobox (HOX) genes, corroborating how changes in external oxygen levels can affect animal development. We discuss how the observed deoxygenation responses may be indicative of a possible acclimation response or alternatively may imply negative latent impacts for coral larval fitness.

Transcriptome Analysis of Durusdinium Associated with the Transition from Free-Living to Symbiotic

To detect the change during coral-dinoflagellate endosymbiosis establishment, we compared transcriptome data derived from free-living and symbiotic Durusdinium, a coral symbiont genus. We detected differentially expressed genes (DEGs) using two statistical methods (edgeR using raw read data and the Student's t-test using bootstrap resampling read data) and detected 1214 DEGs between the symbiotic and free-living states, which we subjected to gene ontology (GO) analysis. Based on the representative GO terms and 50 DEGs with low false discovery rates, changes in Durusdinium during endosymbiosis were predicted. The expression of genes related to heat-shock proteins and microtubule-related proteins tended to decrease, and those of photosynthesis genes tended to increase. In addition, a phylogenetic analysis of dapdiamide A (antibiotics) synthase, which was upregulated among the 50 DEGs, confirmed that two genera in the Symbiodiniaceae family, Durusdinium and Symbiodinium, retain dapdiamide A synthase. This antibiotic synthase-related gene may contribute to the high stress tolerance documented in Durusdinium species, and its increased expression during endosymbiosis suggests increased antibacterial activity within the symbiotic complex.

Cnidarian Pattern Recognition Receptor Repertoires Reflect Both Phylogeny and Life History Traits

Pattern recognition receptors (PRRs) are evolutionarily ancient and crucial components of innate immunity, recognizing danger-associated molecular patterns (DAMPs) and activating host defenses. Basal non-bilaterian animals such as cnidarians must rely solely on innate immunity to defend themselves from pathogens. By investigating cnidarian PRR repertoires we can gain insight into the evolution of innate immunity in these basal animals. Here we utilize the increasing amount of available genomic resources within Cnidaria to survey the PRR repertoires and downstream immune pathway completeness within 15 cnidarian species spanning two major cnidarian clades, Anthozoa and Medusozoa. Overall, we find that anthozoans possess prototypical PRRs, while medusozoans appear to lack these immune proteins. Additionally, anthozoans consistently had higher numbers of PRRs across all four classes relative to medusozoans, a trend largely driven by expansions in NOD-like receptors and C-type lectins. Symbiotic, sessile, and colonial cnidarians also have expanded PRR repertoires relative to their non-symbiotic, mobile, and solitary counterparts. Interestingly, cnidarians seem to lack key components of mammalian innate immune pathways, though similar to PRR numbers, anthozoans possess more complete immune pathways than medusozoans. Together, our data indicate that anthozoans have greater immune specificity than medusozoans, which we hypothesize to be due to life history traits common within Anthozoa. Overall, this investigation reveals important insights into the evolution of innate immune proteins within these basal animals.

Host-Parasite Interaction of Atlantic salmon (Salmo salar) and the Ectoparasite Neoparamoeba perurans in Amoebic Gill Disease

Marine farmed Atlantic salmon (Salmo salar) are susceptible to recurrent amoebic gill disease (AGD) caused by the ectoparasite Neoparamoeba perurans over the growout production cycle. The parasite elicits a highly localized response within the gill epithelium resulting in multifocal mucoid patches at the site of parasite attachment. This host-parasite response drives a complex immune reaction, which remains poorly understood. To generate a model for host-parasite interaction during pathogenesis of AGD in Atlantic salmon the local (gill) and systemic transcriptomic response in the host, and the parasite during AGD pathogenesis was explored. A dual RNA-seq approach together with differential gene expression and system-wide statistical analyses of gene and transcription factor networks was employed. A multi-tissue transcriptomic data set was generated from the gill (including both lesioned and non-lesioned tissue), head kidney and spleen tissues naïve and AGD-affected Atlantic salmon sourced from an in vivo AGD challenge trial. Differential gene expression of the salmon host indicates local and systemic upregulation of defense and immune responses. Two transcription factors, znfOZF-like and znf70-like, and their associated gene networks significantly altered with disease state. The majority of genes in these networks are candidates for mediators of the immune response, cellular proliferation and invasion. These include Aurora kinase B-like, rho guanine nucleotide exchange factor 25-like and protein NDNF-like inhibited. Analysis of the N. perurans transcriptome during AGD pathology compared to in vitro cultured N. perurans trophozoites, as a proxy for wild type trophozoites, identified multiple gene candidates for virulence and indicates a potential master regulatory gene system analogous to the two-component PhoP/Q system. Candidate genes identified are associated with invasion of host tissue, evasion of host defense mechanisms and formation of the mucoid lesion. We generated a novel model for host-parasite interaction during AGD pathogenesis through integration of host and parasite functional profiles. Collectively, this dual transcriptomic study provides novel molecular insights into the pathology of AGD and provides alternative theories for future research in a step towards improved management of AGD.

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.