Community shifts in the coral Oculina patagonica holobiont in response to confinement, temperature and vibrio infections

doi:10.21203/rs.3.rs-3893459/v1

Download PDF

Research Article

Community shifts in the coral Oculina patagonica holobiont in response to confinement, temperature and vibrio infections

https://doi.org/10.21203/rs.3.rs-3893459/v1

This work is licensed under a CC BY 4.0 License

Journal Publication

published 29 Oct, 2024

Read the published version in Microbiome →

You are reading this latest preprint version

Background

Extensive research on the diversity and functional roles of the microorganisms associated with reef-building corals has been promoted as a consequence of the rapid global decline of coral reefs attributed to climate change. Several studies have highlighted the importance of coral-associated algae (Symbiodinium) and bacteria and their potential roles in promoting coral host fitness and survival. However, the complex coral holobiont extends beyond these components to encompass other entities such as protists, fungi, and viruses. While each constituent has been individually investigated in corals, a comprehensive understanding of their collective roles is imperative for a holistic comprehension of coral health and resilience.

Results

A metagenomic analysis was performed to investigate the entire microbiome of the coral Oculina patagonica under various stress conditions, including aquaria confinement, heat stress, and Vibrio infections. Results revealed significant shifts in the abundances and richness of dominant bacterial phyla within the microbiome. In general, aquaria confinement led to an increase in Rhodobacteraceae, Flavobacteraceae, and Vibrionaceae families, accompanied by a decrease in Streptomycetaceae. Concurrently, it was observed that both the abundance and richness of microeukaryotic species significantly declined. Some fungi such as Aspergillus, Fusarium and Rhizofagus genera decrease (or increased, in the case of Saccharomyces). Furthermore, our analysis also identified a reduction in genes associated with antimicrobial compound production by the coral itself, as well as by Symbiodinium and fungi.

Conclusion

Our findings suggest that the interplay between microeukaryotic and prokaryotic components of the coral holobiont may be disrupted by stress conditions, such as confinement, increase of seawater temperature or Vibrio infection, leading to a dysbiosis in the global microbial community that may increase coral susceptibility to diseases. Further, microeukaryotic community seems to exert influence on the prokaryotic community dynamics, possibly through predation or the production of secondary metabolites with anti-bacterial activity.

Oculina patagonica

coral

holobiont

metagenome

Vibrio

Fungi

Corals form a dynamic meta-organism known as the coral holobiont, which involves a multipartite relationship between the cnidarian host, its endosymbiotic dinoflagellate algae (family Symbiodiniaceae; [1]), a diverse array of prokaryotes (Archaea and Bacteria), viruses, and eukaryotes (fungi and non-Symbiodiniaceae protists) [2–4]. Symbiodiniaceae, the primary photosymbiont in corals, has long been recognized for their crucial role in fulfilling most of the host energy needs, facilitating effective calcification and the formation of modern reefs [5]. Prokaryotes also play a pivotal role in the health and fitness of coral holobionts and, although most of their specific roles and functions remain unknown, they are involved in nutrient acquisition, production of bioactive secondary metabolites, and probiotic mechanisms such as competitive exclusion of pathogenic bacteria [6].

All these microorganisms coexist and interact with their host, maintaining the coral´s health and facilitating its ability to adapt to climate change. When environmental stressors, such as rising temperatures, occur, coral species often lose Symbiodiniaceae algae, leading to coral bleaching due to an energy deficit [7–9]. Additional changes can range from the replacement of mutualistic species by commensalism or parasitic species, as well as an increase in potential pathogens such as Vibrio species [10–12]. The simultaneous occurrence of these changes makes it difficult to determine whether the variations in microbial composition in thermally stressed corals are a consequence or cause of bleaching.

The variability observed in coral microbial composition across seasons and spatial scales poses a significant challenge in unraveling the specific roles of microbial species within the coral holobiont under natural conditions. Consequently, many studies aiming to reduce complexity conduct experiments under controlled conditions; e.g., responses to pH gradients, nutrient fluxes, or thermal stress [13–16]. However, the applicability of data obtained from these controlled settings to natural conditions is still not well understood. For example, in coral maintained in aquaria some taxa rapidly decline, or even disappear, while other microorganisms remain unaffected for extended periods [17, 18]. This phenomenon has also been observed in incubations of seawater samples, where a reduction in microbial diversity occurs, accompanied by the replacement of autotrophic microbes with heterotrophs [19–21]. Therefore, within confined corals, additional factors, not well-understood and beyond nutrient imbalances or temperature variations, may contribute to the observed variations.

Omics tools have revolutionized coral research, allowing for a deeper understanding of the various partnerships within the coral holobiont. However, metagenomic approaches aiming to describe the entire holobiont face technical challenges, particularly due to the limited availability of genomes for individual holobiont components, especially microeukaryotes. In fact, only a few coral metagenomes have been sequenced, and most of these studies have analyzed specific partnerships separately, such as focusing on only the prokaryotic [22, 23] or the viral community [24].

In this work, we conducted a metagenomic analysis of the entire Oculina patagonica microbiome. This coral, known as an opportunistic colonizer, is currently expanding throughout the Mediterranean Sea coast [25–27]. The expansion of this coral is attributed in part to its broad tolerance for varying light levels and trophic conditions [28, 29], as well as its heightened thermal resilience compared to other scleractinian corals [27, 30]. However, O. patagonica is still susceptible to high temperature, which can lead to changes in the coral microbiome, coral bleaching, and, in some cases, irreversible mortality [28, 31–33]. In the Western Mediterranean Sea, these coral bleaching episodes have been linked to the presence of two Vibrio species (Vibrio coralliilyticus and Vibrio mediterranei), whose pathogenicity increases when both bacteria infect the coral simultaneously [28].

Here, we provide the first characterization of the O. patagonica microbiome and investigate its changes in response to confinement, thermal stress and the presence of two Vibrio pathogens. To achieve this, we conducted a comparative analysis between the metagenome of an O. patagonica sample recovered in its natural environment and those from the same sampled coral maintained under controlled aquarium experiments. Our study sheds light on how alterations in the microeukaryotic community affect the prokaryotic community and whether these microbiome shifts result in a dysbiosis that promotes greater disease susceptibility in corals.

O. patagonica samples. Twenty-five fragments of apparently healthy specimens, each approximately 5 cm in diameter, of the coral O. patagonica were collected in June 2016 from the Marine Reserve of Tabarca (Mediterranean Sea, 38°09'60.00" N 0°27'59.99" E, Spain) and transported to the laboratory. One coral fragment was used to assess the natural holobiont and is referred to as the Ocu-2016 dataset. The remaining 24 coral fragments were utilized in the aquaria experiment, as described in Rubio-Portillo et al. (2020) [34]. In brief, corals were distributed in two sets of four tanks, with one set incubated at 20°C and another set at 28°C. Each set of tanks contained the following treatments: (I) control corals with no Vibrio added; (II) corals inoculated with the Vibrio coculture; (III) corals with a Vibrio coculture growing inside a dialysis membrane, and (IV) corals in the presence of monoculture from each of the two Vibrio species growing independently inside a dialysis membrane. Dialysis membranes prevented contact between the coral and the Vibrio cells but allowed the free dilution of molecules with a molecular weight lower than 100 kDa. Vibrios used in this experiment were Vibrio mediterranei (strain Vic-Oc-097) and Vibrio corallilyticus (strain Vic-Oc-068).

DNA extraction and metagenomes sequencing. Coral fragments were gently washed three times with 50 ml of sterile filtered seawater (SFSW) to remove any adventitious microbes. Coral pieces were broken into small pieces and centrifuged for 3 min at 2,900 g (Labofuge 400R, Heraeus instruments) in order to remove the mucus from the coral. After centrifugation, the coral pieces were crushed in SFSW using a mortar and pestle, the CaCO₃ skeleton was allowed to settle for 15 min and the supernatant (i.e., the crushed tissue) was collected and used for DNA extraction. DNA from the crushed coral tissue was extracted using the UltraClean Soil DNA Kit (Mo Bio; Carlsbad, CA, United States) following the manufacturer’s instructions for maximum yield.

The O. patagonica sample (Ocu-2016) was sequenced using paired-end reads on an Illumina MySeq (Fisabio, Spain) and also a HiSeq4000 platform (Macrogen, South Corea), with a total of 69 Gb of sequence. Reads from metagenomes obtained from the aquaria experiment were published previously in [34] to assess Vibrio abundance in the corals but metagenomes assemblages reported herein have not been published.

Metagenomic analysis. Raw reads were processed by trimming adaptor sequences and low-quality ends (quality score < 30) using Prinseq [35]. Trimmed reads were used to assess the coverage of metagenomic sequencing with Nonpareil [36]. Subsets of ten million reads from each dataset were randomly separated to perform a taxonomic affiliation using Kaiju [37].

Reads were assembled using MEGAHIT [38] for Ocu-2016 metagenome and meta-SPADES [39] in the case of the aquaria experiments. Contigs larger than 5 Kb were selected for metagenome binning using two programs, MetaBAT [40] and MaxBin 2.0 [41]. The final optimized bins were obtained with DAS Tool 1.0 [42]. Genome completeness and putative contamination within the MAGs were determined with Anvi’o v.2.1.0 [43]. According to the standards suggested by Konstantinidis et al. [44], the quality of the MAGs was considered acceptable when completeness was higher than 80% and contamination lower than 5%. MAGs were manually checked for consistent coverage and taxonomy across contigs. The post-curated MAGs were subjected to fast pairwise comparison via a Mash algorithm [45] and grouped into primary clusters at 95% average nucleotide identity (ANI). After this de-replication step, taxonomic affiliation of MAGs was obtained using the GTDB-Tk pipeline, which uses a marker gene set to place genomes in the GTDB reference tree [46]. Gene predictions on the individual reads and also from the contigs of the MAGs were carried out using Prodigal [47]. CDS from the contigs bigger than 1 Kb were annotated with EggNOG v5.0 [48]. Functional annotation of the predicted ORFs from contigs and MAGs was also performed by BLASTp comparisons against the NCBI nr database, Pfam [49], COG [50], and TIGRFAM [51]. To reconstruct the metabolic potential for the MAGs, GhostKOALA was used [52]. The search for genes encoding biosynthesis of secondary metabolites was performed with antiSMASH v. 3.0 [53]. For the identification of putative enzymes involved in the breakdown, biosynthesis or modification of carbohydrates, the CDSs from the contigs bigger than 1 Kb were compared against the CAZy database (dbCAN) [54].

Recruitment analyses. To track the relative abundance of each MAG in the different samples, BBMap was used to map metagenomic short reads to MAG contigs [55]. Matches were filtered for single best alignments, using a minimum of 90% query cover alignment length and 95% nucleotide identity of reads mapping against the reference genome (ANIr). In order to remove biases from highly conserved regions and contig edges, the 80% central truncated average of sequencing depth of all bases (TAD80) was used as described previously [56]. MAG abundance in each metagenomic dataset (as percentage of total community) was calculated as the quotient of the MAG’s TAD80 value and the genome equivalents (GE) from MicrobeCensus [57].

In the case of the recruitment analysis of the dinoflagellate genomes against each of the metagenomes, a database containing 20 dinoflagellate assemblies was constructed (Breviolum minutum Mf 1.05b.01 (GCA_000507305.1), Symbiodinium microadriaticum (GCA_001939145.1), Symbiodinium microadriaticum (GCA_905231925.1), Symbiodinium microadriaticum AJIS2-C2 (GCA_018327485.1), Symbiodinium sp. clade A Y106 (GCA_003297005.1), Symbiodinium sp. clade C Y103 (GCA_003297045.1), Symbiodinium kawagutii (GCA_009767595.1), Symbiodinium natans (GCA_905221605.1), Symbiodinium tridacnidorum CCMP2592 (GCA_905221615.1), Symbiodinium sp. KB8 (GCA_905221625.1), Symbiodinium linuchaea CCMP2456 (GCA_905221635.1), Symbiodinium pilosum (GCA_905231905.1), Symbiodinium necroappetens (GCA_905231915.1), Amoebophrya sp. AT5.2 (GCA_005223375.1), Amoebophrya sp. A120 (GCA_905178155.14), Amoebophrya sp. A25 (GCA_905178165.1), Polarella glacialis (GCA_905237085.1), Polarella glacialis (GCA_905237095.1), Cladocopium goreaui (GCA_947184155.1) and Amphidinium carterae (GCA_019702695.1). BlastN comparisons of subsets of 10 million reads were performed and all the reads that hit over 90% identity with a coverage over 70% were considered. When a read matched several dinoflagellate genomes equally well, it was considered as a match for all genomes using a modification of the BlastTab.best_hit_sorted of Enveomics Toolbox. To calculate the frequency of the dinoflagellates in the metagenomes, the number of hits were normalized against the dinoflagellate size-dataset.

Comparison with other environmental datasets.

For the MetaFast comparison [58], a subset of ten million reads randomly extracted from each collection was first obtained. The metagenome collections used came from Ocu-2016, samples from the aquaria experiments, together with other metagenomes obtained from different corals (Acropora palmata (SRR8789034), Cyphastrea (Australia, PRJNA364879), Diploira (Mexico, PRJNA364874), Siderastrea (Mexico, PRJNA364873), a sponge (Ircinia ramosa, PRJNA555144), a metagenome from soil (Soil Forest Indian, JGI-Ga0247671), several seawater metagenomes (two from the TARA collections: TARA-009, 5 m, ERR594288 and TARA-009, 55 m, ERR594315) [59] [59] [59] [59] [59] [59]; 4 metagenomes from the MED collection [60]: Med-OCT2015-15m (SRR5007106); Med-OCT2015-45m (SRR5007115) and two seawater metagenomes affected by a “bottle-effect”, Med-OCT2015-15m-7h (SRR8503606) and Med-OCT2015-15m-14h (SRR8503605).

Viral sequence Analyses Viral contigs were included if their length was over ≥ 10 kb and they were identified as a viral sequence of category 1 or 2 by VIRSorter v1.03 [61]. Also, contigs which were identified as prophages by VIRSorter v1.03 were also annotated with PHASTER [62] and only those contigs identified as intact or questionable prophages were used. To cluster viral sequences into viral operational taxonomic units (vOTUs), CD-HIT was used with a cut-off of 95% average nucleotide identity and 80% alignment fraction threshold for the smallest contig [63]. vConTACT2 0.9.09 was used to compare the viral contigs using the Prokaryotic Viral RefSeq version 201 (with ICTV and NCBI taxonomies). Only those connections with a score > 1 were considered. The nets were visualized with Cytoscape [64].

To characterize the coral holobiont community, encompassing eukaryotes, prokaryotes, and viruses, a metagenome of approximately 70 Gb of sequence from a sample of O. patagonica collected during the summer of 2016 in Tabarca (Alicante, Spain) was sequenced. This is referred to as the Ocu-2016 dataset. Additionally, to investigate alterations in these communities in response to heat stress or the presence of Vibrio coral pathogens, we conducted complementary analyses using published metagenomes from Rubio-Portillo et al. [34]. The published metagenomes came from the same coral sample that had been split into several experimental regimes, though the published analysis primarily focused on characterizing the Vibrio assemblages and did not encompass the complete holobiont. For all nine metagenomes, non-pareil analysis indicated a coverage of approximately 94% for Ocu-2016 and between 50 and 68% for the metagenomes from the aquaria (Table 1), indicating a relatively good representation of the biodiversity contained within them.

Oculina patagonica holobiont

First, a comparison of all O. patagonica metagenomes was conducted against other published coral and Mediterranean seawater datasets. As expected, samples from the seawater clustered separately from the coral samples, which formed a cluster comprising metagenomes from O. patagonica and other corals metagenomes (Supplementary Fig. 1).

The reads annotated within the Ocu-2016 coral metagenome were distributed as follows: 10.86% belonged to Bacteria, 3.85% to Eukarya, 0.43% to Archaea and 0.07% to Viruses, with 84.78% remaining unclassified (Supplementary Table 1). To assess whether this pattern was specific to our sample, a similar analysis was conducted using a dataset from the coral Porites lutea (SRR9182857) [23]. The results were comparable, with only 1.4% of the reads identified as eukaryotes and 7.7% classified as prokaryotes (Supplementary Table 1). The low percentage of eukaryotes found suggests that, due to the relatively lower availability of genomes from eukaryotes compared to prokaryotes, it is possible that some of the unclassified reads may belong to the coral and other microeukaryotic organisms.

O. patagonica -associated eukaryotes

Among coral symbionts, Symbiodoniaceae were the first and most important to be recognized, and mutual transport of nutrients between both taxa has been well described [1]. O. patagonica harbours Symbiodinium species belonging to three different clades (A, B and C) from which representative genomes are available (Fig. 1). Among them, the genome of Breviolum minutum Mf1.05b (also known as Symbiodinium minutum), which belongs to clade B, exhibited the highest read recruitment rates. This result was consistent with previous findings that identified Symbiodinium type B1 as the primary clade in O. patagonica [25]. It is also consistent with prior investigations in other coral species, demonstrating the concurrent association of corals with multiple Symbiodiniaceae, with a prevailing clone [65].

Here it was found that protista other than Symbiodinium were also abundant in the O. patagonica holobiont, with the Evosea, Euglenozoa, and Apicomplexa phyla being particularly prevalent (Fig. 2). Previous studies have already documented the presence of other protists in the coral holobiont, suggesting that they may play a role in assisting coral hosts in obtaining sufficient nutrients or serve as an additional food source during recovery from stress, such as tissue loss and bleaching events [66]. Among the detected protists were two Choanoflagellate species, Monosiga brevicollis and Salpingoeca rosetta, which were not previously known to be associated with corals. The unexpected occurrence of reads matching the parasitic protist Plasmodium may have been attributed to the presence of the "apicoplast", an organelle in Plasmodium and other apicomplexans [67] that is derived from endosymbiotic cyanobacteria. In other coral reef samples, sequences from apicomplexan-related plastids from Chromera and Vitrella have been consistently detected [68, 69], but given the low recruitment of these genera in Ocu-2016, it is possible that many of the reads matching Plasmodium belong to a yet unknown microalgae plastid(s) present in O. patagonica.

Numerous non-photosynthetic microeukaryotes associated with corals have been identified, exhibiting diverse roles within the holobiont, ranging from beneficial symbiosis to parasitic relationships, and even acting as primary pathogens [70–72]. In this study, most of the non-photosynthetic microeukaryotic reads were assigned as Fungi. Among them, the most prominently represented phyla were Ascomycota, Basidiomycota, and Mucoromycota (Fig. 2), in good agreement with previous results reporting Ascomycetes and Basidiomycetes as the two major groups of fungi associated with corals [71, 73]. Within the Ascomycota, the genera Aspergillus and Fusarium were the most abundant in O. patagonica. Also detected was a significant abundance of reads assigned to the genus Rhizophagus (Glomeraceae family), a beneficial mycorrhizal fungus commonly used as a soil inoculant in agriculture and forest ecosystems to enhance phosphorus uptake [74]. Recruitment analysis of Ocu-2016 reads against the reference genome of Rhizophagus irregularis (DAOM 181602) revealed the presence of reads with similarity to their ribosomal operons and housekeeping genes. These findings suggest that although Rhizophagus species have not been previously described in aquatic environments, related fungi may possess beneficial traits in corals, as similar results were also obtained using the microbial metagenomes of the corals P. lutea and Acropora palmata (results not shown).

It has been suggested that microeukaryotes play a role in regulating microbial communities within corals through processes like phagocytosis and the production of antimicrobial compounds [2, 69, 71]. Searching for the presence of genes that encode the synthesis of secondary metabolites was done to investigate the potential antimicrobial activities of O. patagonica microbiome. Detected genes included those that encode polyketide synthases (PKSs) and multienzymatic nonribosomal peptide synthetases (NRPSs). When the PKSs genes in Ocu-2016 were analyzed according to their taxonomy, a high number of these sequences belonged almost exclusively to Eukarya (Symbiodinium, Scleractinia, Blastocladiomycota and Dictyosteliales), whereas no bacterial PKSs were detected (Fig. 3).

O. patagonica associated prokaryotes

Among the prokaryotic community, bacterial reads greatly outnumbered the archaea, maintaining a ratio of approximately 15:1 (Fig. 2 and Supplementary Table 1). The archaeal community was composed of 0.2% Euryarchaeota and 0.1% Thaumarchaeota, with the most abundant family being Nitrosopumilaceae (Fig. 2). The genus Nitrosopumilus exhibited the highest abundance within this family (Fig. 4). This specific genus has been previously linked to the crucial ammonium oxidation process occurring within the coral's mucus layer [75]. Notably, Nitrosopumilus has been previously detected in O. patagonica, and its presence has been exclusively associated with healthy colonies [76]. Also in accordance with previous findings obtained by 16S rRNA [76], the Bacteria were primarily composed of Proteobacteria (4.3%), Actinobacteria (1.6%), Firmicutes (1.2%), and Bacteroidetes (1.1%) (Fig. 2), with Streptomyces, Pseudomonas and Bacillus as the most abundant genera (Fig. 4).

Binning of the assembled contigs from the Ocu-2016 metagenome resulted in five metagenome assembled genomes (MAGs), which belonged to Desulfobacterales (MAG1-Ocu2016), Flavobacteriales (MAG2-Ocu2016), Holosporales (MAG3-Ocu2016), Parvularculales (MAG4-Ocu2016) and Rhizobiales (MAG5-Ocu2016) (Supplementary Table 2). These MAGs showed recruitment rates below 0.001% (Fig. 5), and ANIr values over 99% (Supplementary Table 3), suggesting that only MAGs from low abundance bacteria displaying a very low intra-population diversity could be retrieved. Excluding MAG1, the average genome size for the other four MAGs from the Ocu-2016 dataset was 1.7 ± 0.5 Mb, a size small enough to be consistent with a host-associated lifestyle, reflecting the loss of non-essential genes [77]. For example, MAG3-Ocu2016 (Candidatus Hepatobacter penai; 87.3% complete) lacks key biosynthetic pathways for essential amino acids and shows incomplete synthesis of purines and pyrimidines (Supplementary Table 4). However, it harbors genes for biotin (vitamin B7) synthesis. This could be crucial for the coral holobiont, as both corals and Symbiodinium, are suggested to be auxotrophs for different B vitamins, potentially obtaining them from associated bacteria [78, 79].

Furthermore, within the MAGs, genes encoding mechanisms for stable symbiosis with the host, such as ankyrin repeats proteins (ARPs), were identified (Supplementary Table 2). ARPs are common protein interaction motifs that modulate intracellular processes, promoting stable symbiotic or pathogenic associations [80]. Previous analysis of ARP distribution in microbial genomes showed that species dedicating more than 0.2% of their protein-coding genes to ARPs are typically obligate intracellular or facultative host-associated species [81]. In this case, MAG3-Ocu2016 and MAG4-Ocu2016, both from the Alphaproteobacteria class, surpass this percentage. Moreover, in the case of MAG4-Ocu2016 (Unknown Parvularculales), these proteins were affiliated with other corals and sponges, suggesting its potential as a symbiont for marine invertebrates.

MAG1-Ocu2016 represents a dissimilatory sulfate-reducing bacterium (SRB) belonging to the Desulfobacter genus; its genome harbored genes responsible for nitrogen fixation (Supplementary Table 4). The identification of a diazotrophic SRB, not previously described in corals, carries substantial implications for coral communities that commonly inhabit nutrient-depleted environments. These bacteria may enhance the coral's ability to efficiently convert gaseous nitrogen into a more usable form such as ammonia, a vital process for sustaining a steady nitrogen supply for Symbiodinium-based primary production within corals [82].

O. patagonica associated virus

A total of 275 different viral operational taxonomic units (vOTUs) were retrieved from Ocu-2016 metagenome, which accounted for 0.08% of the metagenomics reads. Only 14% of these vOTUs exhibited similarity to known reference genomes and most of them (86%) were classified as dsDNA phages belonging to the class Caudoviricetes (Supplementary Table 5). A phage phylogenetic network using phage reference genomes showed that most of the phages identified in the O. patagonica holobiont were new. Others grouped mainly with phages infecting Pseudoalteromonas, Pseudomonas and Roseobacter, bacterial genera consistently associated with corals. Three vOTUs presented a partial similarity to Suoliviridae sequences (Crassvirales), with 30 to 63% similarity to a major head protein, DNA polymerase or hypothetical proteins. Despite marine crassviruses having been described [83], O. patagonica assembled viruses clustered separately from those (Fig. 6).

Five of the identified vOTUs were assigned as eukaryotic viruses belonging to the Phycodnaviridae and Mimiviridae families within the Nucleocytoplasmic Large DNA Viruses (NCLDVs) (Supplementary Table 5). The family Phycodnaviridae infects phytoplankton and has been previously detected in both heat-stressed corals and cultures of Symbiodinium spp., suggesting that these viruses may have a role in the destruction of algal symbionts or the dysfunction of symbiont–host mutualism, although the extent of such infections is unknown [84–86].

Oculina patagonica holobiont changes under experimental conditions

When compared with collection of other metagenomes, the Ocu-2016 dataset clustered alongside the coral samples maintained at 20ºC, separate from the 28ºC samples (Supplementary Fig. 1). Read-level analysis further supported these findings, with Ocu-2016 sharing 46.9% of reads with C20, only 15% with C28, and lower similarities were observed when compared to the Vibrio infection datasets (Fig. 7). These results highlight significant shifts in the microbial community of O. patagonica under heat stress, particularly in the presence of Vibrio species. These community shifts in the aquaria confinements were mainly enrichments of heterotrophic microbes. This phenomenon, known as the "bottle effect," has been previously observed when seawater samples or corals are kept under prolonged confinement conditions [18, 87].

Changes due to confinement

It is important to note that the proportion of reads assigned to Eukarya within the metagenome from the coral maintained at 20ºC decreased to 2.2% compared to Ocu-2016, while the Bacteria reads increased to 17% (Supplementary Table 1). In aquaria, the analysis indicated a decline in Symbiodinium, as well as other microeukaryotes and prokaryotes compared to Ocu-2016. Meanwhile an expansion of the rare biosphere was observed under confinement conditions (Fig. 4). Microeukaryotes that proliferated in the aquarium experiments were mainly members of the Saccharomycodaceae family, which were not detected in Ocu-2016 (Fig. 2). The majority of these reads showed similarity to the yeast Saccharomycodes ludwigii, with some matching the Hanseniaspora genus. The presence of these sugar-consuming microorganisms may be attributed to the increased carbohydrate content in the mucus of stressed corals [88].

Also under confinement conditions, there was a notable increase in certain prokaryotic genera, including potential human pathogens like Acinetobacter, Bordetella, Neisseria, Klebsiella, and Salmonella, which were not detected in the seawater recovered O. patagonica (Fig. 4). Although speculative, this shift may be linked to changes in the microeukaryotic community associated with the coral. As previously commented, in Ocu-2016, most of the detected polyketide synthases (PKSs) primarily belonged to Eukarya (Symbiodinium, Scleractinia, Blastocladiomycota, and Dictyosteliales); abundances of these PKSs decreased under confinement conditions (Fig. 3). The absence or reduction of these natural eukaryotic PKS could potentially lead to an increase in the proliferation of fast-growing organisms, as discussed below. This would be in agreement with previous results indicating that the coral-associated fungi that decreased under confinement conditions, such as species from Aspergillus or Fusarium, displayed antimicrobial activities against human pathogenic bacteria [89, 90].

Changes due to thermal stress under experimental conditions

The proportion of annotated eukaryal reads in the metagenome from the coral maintained at 28ºC remained similar to that at 20ºC. However, there was a substantial increase in the bacterial reads, up to 40%. The number of detected bacterial families also increased, from 112 at 20ºC to 171 at 28ºC (Supplementary Table 1).

Under heat stress, various genera from the Rhodobacteraceae family, including Marivitia, Ruegeria, Loktanella, and Yoonia, were identified as bacterial bloomers (Fig. 4). In line with our results, the Rhodobacteraceae family has been recently suggested as an indicator species for thermal stress in corals [91]. Furthermore, a significant number of Rhodobacteraceae genera detected in our study are known to be involved in the breakdown of organic sulfur compounds like dimethylsulfoniopropionate (DMSP) and dimethylsulfide (DMS). The production of these compounds increases in corals under thermal stress [92–94]. In good agreement, genes responsible for DMSP catabolism (DmdA, DmdB, DmdC, DddD, DddP, and DddL) increased in abundance in corals kept in aquaria, particularly in corals maintained at 28ºC (from 0.00001% in Ocu-2016 to 0.0001–0.001% of reads). Furthermore, two MAGs recovered from the aquaria metagenomes (MAG9-MS28 and MAG10-MX28, see Supplementary Table 2) corresponded to two Rhodobacteraceae genera potentially involved in DMSP metabolism, Pelagibaca and Yoonia, respectively. In fact, an orthologue gene encoding the DddL enzyme, responsible for DMSP catabolism, was detected in a Yoonia-related MAG, with a more pronounced increase in coral samples maintained at 28ºC (Fig. 5). This suggests that sulfur compounds produced in response to stressors, like confinement or thermal stress, may contribute to shaping coral-associated bacterial communities. This supports the hypothesis proposed by [95] that DMSP and DMS play a pivotal role in structuring coral-associated bacterial communities.

Other genera that increase under confinement conditions and particularly under heat stress were Marinifilum, a member of the Bacteroidetes genus, and Halodesulfovibrio, originally affiliated with Desulfovibrio genus (Fig. 4). These genera have previously been observed in coral samples, and their proliferation under aquarium conditions has been documented by [96]. Marinifilum may play a role in the sulfur and carbon cycle, as well as lipid catabolism. On the other hand, Halodesulfovibrio has been identified as a secondary pathogen responsible for initiating and progressing black band disease in coral hosts, producing sulfide as a product of dissimilatory sulfate reduction [97].

The increase of bacterial bloomers under heat stress may be partially influenced by the recycling of metabolic waste products within the holobiont [98] and subsequent increase of available nutrients. Additionally, corals regularly release mucus into the surrounding seawater, which carries elevated concentrations of nutrients like organic carbon [99]. In heat-stressed corals, an increase in carbohydrate content in mucus is detected [88], thus serving as a carbon source for heterotrophic microbes. To explore this idea, variation in the abundances of genes encoding carbohydrate-active enzymes (CAZymes) was examined. Results showed an increase in the proportion of glycosylhydrolases, glycosyltransferases, and carbohydrate binding modules in corals maintained under experimental conditions, particularly in the thermally stressed ones, compared to Ocu-2016 sample (Supplementary Fig. 2). This suggests that mucus released under stressful conditions could provide newly accessible nutrients, serving as sustenance for fast-growers microbes. This process could potentially be enhanced by the decrease of other bacterial taxa and microeukaryotes under stress conditions, which may act as microbial regulators within the coral holobiont, as explained below.

Regarding viruses, it is well established that VLPs increase in corals during stressful conditions or bleaching events [84, 100]. In this case, recruitment analysis showed that one of the viruses, identified as Phycodnaviridae in the Ocu-2016 metagenome, clearly increased in abundance under thermal stress (Supplementary Fig. 3) when Symbiodinium abundances decreased. This fact suggests that this virus could be involved in in the destruction of algal symbionts or the dysfunction of symbiont–host mutualism as mentioned above.

Changes due to Vibrio infection under experimental conditions

All the above-described changes observed in confined or heat stress were enhanced in corals exposed to Vibrio cells (V. coralliilyticus and V. mediterranea), directly or indirectly. Specifically, corals at 28ºC in the presence of Vibrio pathogens experienced a more significant increase in Vibrio species abundance compared to controls. Accordingly, a previous transcriptomic study carried out with these samples [34] revealed that pathogenic Vibrio release quorum sensing molecules, triggering alterations in coral-associated bacteria and an increase in other potential pathogens already present in the coral sample, thus explaining the increase of Vibrio spp. detected here. Furthermore, the increase in species related to DMSP metabolism was also more pronounced under the presence of Vibrio pathogens, which also suggests vibrios play key role in coral microbiome modulation.

Analysis of the viral contigs in these coral metagenomes revealed a significant decrease in the proportion of genomes classified as prophages in corals experiencing heat stress and in those exposed to Vibrio coral pathogens compared to the natural sample (Fig. 8A). To delve into the lysogenic dynamics within the coral holobiont, recruitment analyses were performed using only prophage genomes identified within bacterial contigs; this approach confirmed their integration into their respective bacterial hosts. Our data unveiled the induction of three distinct prophages in corals subjected to Vibrio infection (I20) or heat-induced stress (C28), and two of these prophages were found to be integrated within contigs belonging to the Rhodobacterales order (see Fig. 8B). These findings align well with recent research suggesting that coral pathogens, such as V. coralliilyticus, produce hydrogen peroxide to initiate the lytic cycle of prophages in their competitors, thereby providing the coral pathogen with an advantage by reducing competition during coral colonization [101]. This mechanism further contributes to host dysbiosis by shifting the balance from symbionts to pathobionts as observed in the human gut [102].

Despite the potential bias due to the absence of eukaryal reference genomes, our findings suggest the existence of a highly diverse microeukaryotic community within O. patagonica. Importantly, many of these microeukaryotic symbionts suffered a significant decline (below detectable levels) in O. patagonica aquaria experiments, particularly at 28ºC. Concurrently, a notable increase in abundance of the prokaryotic community was observed, shedding light on the complex interplay between these two communities. Our data suggest that the microeukaryotic community potentially exerts influence on the prokaryotic community dynamics, possibly through predation or the production of secondary metabolites with anti-bacterial activity.

The increase in the contribution of the prokaryotic community to the holobiont appears to coincide with an increase of CAZymes and genes responsible for DMSP catabolism. This correlation suggests the existence of a potential recycling mechanism for organic products generated by the holobiont under stress conditions. Thus, the alteration of conditions (such as temperature, confinement and pathogen presence) within the coral holobiont appears to enhance the production of enzymes dedicated to carbohydrate degradation. This discovery underscores the dynamic and adaptive nature of microbial communities within coral ecosystems and offers exciting prospects for further research into the exploitation of these enzymatic resources in various industrial and environmental contexts.

Availability of data and materials

The raw sequences of the wild type O. patagonica metagenome dataset (Ocu-2016) have been deposited in the SRA repository under the BioProject number PRJNA661426. Metagenomes from the infection experiments and their respective controls can be found in the [34] (PRJNA612159).

Competing interests

The authors declare that they have no competing interests

Funding

This research was supported in part by the EU-H2020 MetaFluidics and BlueTools projects with grant agreement numbers 685474 and 101081957 (to J.A.) E.R.-P. was funded by the postdoctoral program Valid (GVA) (grant number APOSTD-2016-091).

Authors' contributions

E.R.-P. and J.A. conceived and designed the study. E.R.-P conducted experiments. A.B.M.-C and E.R.-P. analyzed the metagenomic data. A.B.M.-C, E.R.-P and J.A wrote the manuscript, and all of us contributed substantially with discussion.

Acknowledgements

We are grateful to Forest Rohwer for providing funding for aquaria metagenomes sequencing (NSF-PIRE grant number OISE1243541). Special thanks to Heather Maughan for the professional English editing and critical reading of the manuscript.

LaJeunesse TC, Parkinson JE, Gabrielson PW, Jeong HJ, Reimer JD, Voolstra CR, et al. Systematic Revision of Symbiodiniaceae Highlights the Antiquity and Diversity of Coral Endosymbionts. Curr Biol. 2018;28(16):2570-80 e6; doi: 10.1016/j.cub.2018.07.008.
Knowlton N, Rohwer F. Multispecies microbial mutualisms on coral reefs: the host as a habitat. Am Nat. 2003;162(4 Suppl):S51-62.
Rohwer F SV, Azam F, Knowlton N. Diversity and distribution of coral-associated bacteria. Marine Ecology Progress Series. 2002;243:1–10; doi: 10.3354/meps243001.
van Oppen MJH, Blackall LL. Coral microbiome dynamics, functions and design in a changing world. Nat Rev Microbiol. 2019;17(9):557–67; doi: 10.1038/s41579-019-0223-4.
Muscatine LF, P.G.; Porter, J.W.; Dubinsky, Z. Fate of Photosynthetic Fixed Carbon in Light- and Shade-Adapted Colonies of the Symbiotic Coral Stylophora pistillata. Proceedings of the Royal Society of London Series B, Biological Sciences. 1984;222(1227):181–202.
Voolstra CRS, D.J.; Peixoto; R. S.; Parkinson, J.E.; Quigley, K.M.; Silveira, C.B.; Sweet, M.; Muller, E.M.; Barshis, D.J.; Bourne D. G. & Aranda M. Extending the natural adaptive capacity of coral holobionts. Nature Reviews Earth & Environment 2021;2:747–62.
Fine M, Loya Y. Endolithic algae: an alternative source of photoassimilates during coral bleaching. Proc Biol Sci. 2002;269(1497):1205–10; doi: 10.1098/rspb.2002.1983.
Schlichter D, Zscharnack B, Krisch H. Transfer of photoassimilates from endolithic algae to coral tissue. Naturwissenschaften. 1995;82(12):561–4; doi: Doi 10.1007/Bf01140246.
Sangsawang L, Casareto BE, Ohba H, Vu HM, Meekaew A, Suzuki T, et al. (13)C and (15)N assimilation and organic matter translocation by the endolithic community in the massive coral Porites lutea. R Soc Open Sci. 2017;4(12):171201; doi: 10.1098/rsos.171201.
Littman RA, Willis BL, Pfeffer C, Bourne DG. Diversities of coral-associated bacteria differ with location, but not species, for three acroporid corals on the Great Barrier Reef. FEMS Microbiol Ecol. 2009;68(2):152–63.
Bourne D, Iida Y, Uthicke S, Smith-Keune C. Changes in coral-associated microbial communities during a bleaching event. ISME J. 2008;2(4):350–63; doi: 10.1038/ismej.2007.112.
Zhou Z, Zhao S, Tang J, Liu Z, Wu Y, Wang Y, et al. Altered Immune Landscape and Disrupted Coral-Symbiodinium Symbiosis in the Scleractinian Coral Pocillopora damicornis by Vibrio coralliilyticus Challenge. Front Physiol. 2019;10:366; doi: 10.3389/fphys.2019.00366.
Meron D, Rodolfo-Metalpa R, Cunning R, Baker AC, Fine M, Banin E. Changes in coral microbial communities in response to a natural pH gradient. ISME J. 2012;6(9):1775–85; doi: 10.1038/ismej.2012.19.
Lee ST, Davy SK, Tang SL, Fan TY, Kench PS. Successive shifts in the microbial community of the surface mucus layer and tissues of the coral Acropora muricata under thermal stress. FEMS Microbiol Ecol. 2015;91(12); doi: 10.1093/femsec/fiv142.
Li J, Zhou Y, Qin Y, Wei J, Shigong P, Ma H, et al. Assessment of the juvenile vulnerability of symbiont-bearing giant clams to ocean acidification. Sci Total Environ. 2022;812:152265; doi: 10.1016/j.scitotenv.2021.152265.
Klinges JG, Patel SH, Duke WC, Muller EM, Vega Thurber RL. Phosphate enrichment induces increased dominance of the parasite Aquarickettsia in the coral Acropora cervicornis. FEMS Microbiol Ecol. 2022;98(2); doi: 10.1093/femsec/fiac013.
Galand PE, Chapron L, Meistertzheim AL, Peru E, Lartaud F. The Effect of Captivity on the Dynamics of Active Bacterial Communities Differs Between Two Deep-Sea Coral Species. Front Microbiol. 2018;9:2565; doi: 10.3389/fmicb.2018.02565.
Kooperman N, Ben-Dov E, Kramarsky-Winter E, Barak Z, Kushmaro A. Coral mucus-associated bacterial communities from natural and aquarium environments. FEMS Microbiol Lett. 2007;276(1):106–13; doi: 10.1111/j.1574-6968.2007.00921.x.
Calvo-Diaz A, Diaz-Perez L, Suarez LA, Moran XA, Teira E, Maranon E. Decrease in the autotrophic-to-heterotrophic biomass ratio of picoplankton in oligotrophic marine waters due to bottle enclosure. Appl Environ Microbiol. 2011;77(16):5739–46; doi: 10.1128/AEM.00066-11.
Dinasquet J, Kragh T, Schroter ML, Sondergaard M, Riemann L. Functional and compositional succession of bacterioplankton in response to a gradient in bioavailable dissolved organic carbon. Environ Microbiol. 2013;15(9):2616–28; doi: 10.1111/1462-2920.12178.
Massana RP-ACC, E.O.; Gasol, J.M. Changes in marine bacterioplankton phylogenetic composition during incubations designed to measure biogeochemically significant parameters. Limnology and Oceanography. 2001;46:1181–8; doi: doi.org/10.4319/lo.2001.46.5.1181.
Cai L, Tian RM, Zhou G, Tong H, Wong YH, Zhang W, et al. Exploring coral microbiome assemblages in the South China Sea. Sci Rep. 2018;8(1):2428.
Robbins SJ, Singleton CM, Chan CX, Messer LF, Geers AU, Ying H, et al. A genomic view of the reef-building coral Porites lutea and its microbial symbionts. Nat Microbiol. 2019;4(12):2090–100.
Weynberg KD, Wood-Charlson EM, Suttle CA, van Oppen MJ. Generating viral metagenomes from the coral holobiont. Front Microbiol. 2014;5:206; doi: 10.3389/fmicb.2014.00206.
Rubio-Portillo E, Souza-Egipsy V, Ascaso C, de Los Rios Murillo A, Ramos-Espla AA, Anton J. Eukarya associated with the stony coral Oculina patagonica from the Mediterranean Sea. Mar Genomics. 2014;17:17–23; doi: 10.1016/j.margen.2014.06.002.
Leydet KP, Grupstra CGB, Coma R, Ribes M, Hellberg ME. Host-targeted RAD-Seq reveals genetic changes in the coral Oculina patagonica associated with range expansion along the Spanish Mediterranean coast. Mol Ecol. 2018;27(11):2529–43; doi: 10.1111/mec.14702.
Serrano E, Ribes M, Coma R. Demographics of the zooxanthellate coral Oculina patagonica along the Mediterranean Iberian coast in relation to environmental parameters. Sci Total Environ. 2018;634:1580–92; doi: 10.1016/j.scitotenv.2018.04.032.
Rubio-Portillo E, Yarza P, Penalver C, Ramos-Espla AA, Anton J. New insights into Oculina patagonica coral diseases and their associated Vibrio spp. communities. ISME J. 2014;8(9):1794–807; doi: 10.1038/ismej.2014.33.
Zaquin TZ, P.; Pinkas, I.; Mass, T. Simulating Bleaching: Long-Term Adaptation to the Dark Reveals Phenotypic Plasticity of the Mediterranean Sea Coral Oculina patagonica. Frontiers in Marine Science. 2019;6; doi: doi.org/10.3389/fmars.2019.00662.
Rodolfo-Metalpa R, Hoogenboom MO, Rottier C, Ramos-Espla A, Baker AC, Fine M, et al. Thermally tolerant corals have limited capacity to acclimatize to future warming. Glob Chang Biol. 2014;20(10):3036–49; doi: 10.1111/gcb.12571.
Rubio-Portillo E, Kersting DK, Linares C, Ramos-Espla AA, Anton J. Biogeographic Differences in the Microbiome and Pathobiome of the Coral Cladocora caespitosa in the Western Mediterranean Sea. Front Microbiol. 2018;9:22; doi: 10.3389/fmicb.2018.00022.
Armoza-Zvuloni RS, R; Kramarsky-Winter, E; Loya, Y Repeated bleaching events may result in high tolerance and notable gametogenesis in stony corals: Oculina patagonica as a model. Marine Ecology Progress Series. 2011;426:149–59; doi: 10.3354/meps09018
Shenkar NF, M.; Loya Y. Size matters: bleaching dynamics of the coral Oculina patagonica. Marine Ecology Progress Series. 2005;294:181–8.
Rubio-Portillo E, Martin-Cuadrado AB, Caraballo-Rodriguez AM, Rohwer F, Dorrestein PC, Anton J. Virulence as a Side Effect of Interspecies Interaction in Vibrio Coral Pathogens. mBio. 2020;11(4); doi: 10.1128/mBio.00201-20.
Schmieder R, Edwards R. Quality control and preprocessing of metagenomic datasets. Bioinformatics. 2011;27(6):863–4; doi: 10.1093/bioinformatics/btr026.
Rodriguez RL, Konstantinidis KT. Nonpareil: a redundancy-based approach to assess the level of coverage in metagenomic datasets. Bioinformatics. 2014;30(5):629–35; doi: 10.1093/bioinformatics/btt584.
Menzel P, Ng KL, Krogh A. Fast and sensitive taxonomic classification for metagenomics with Kaiju. Nat Commun. 2016;7:11257; doi: 10.1038/ncomms11257.
Li D, Luo R, Liu CM, Leung CM, Ting HF, Sadakane K, et al. MEGAHIT v1.0: A fast and scalable metagenome assembler driven by advanced methodologies and community practices. Methods. 2016;102:3–11; doi: 10.1016/j.ymeth.2016.02.020.
Nurk S, Meleshko D, Korobeynikov A, Pevzner PA. metaSPAdes: a new versatile metagenomic assembler. Genome Res. 2017;27(5):824–34; doi: 10.1101/gr.213959.116.
Kang DD, Froula J, Egan R, Wang Z. MetaBAT, an efficient tool for accurately reconstructing single genomes from complex microbial communities. PeerJ. 2015;3:e1165; doi: 10.7717/peerj.1165.
Wu YW, Tang YH, Tringe SG, Simmons BA, Singer SW. MaxBin: an automated binning method to recover individual genomes from metagenomes using an expectation-maximization algorithm. Microbiome. 2014;2:26; doi: 10.1186/2049-2618-2-26.
Sieber CMK, Probst AJ, Sharrar A, Thomas BC, Hess M, Tringe SG, et al. Recovery of genomes from metagenomes via a dereplication, aggregation and scoring strategy. Nat Microbiol. 2018;3(7):836–43; doi: 10.1038/s41564-018-0171-1.
Eren AM, Esen OC, Quince C, Vineis JH, Morrison HG, Sogin ML, et al. Anvi'o: an advanced analysis and visualization platform for 'omics data. PeerJ. 2015;3:e1319; doi: 10.7717/peerj.1319.
Konstantinidis KT, Rossello-Mora R, Amann R. Uncultivated microbes in need of their own taxonomy. ISME J. 2017;11(11):2399–406; doi: 10.1038/ismej.2017.113.
Ondov BD, Treangen TJ, Melsted P, Mallonee AB, Bergman NH, Koren S, et al. Mash: fast genome and metagenome distance estimation using MinHash. Genome Biol. 2016;17(1):132; doi: 10.1186/s13059-016-0997-x.
Chaumeil PA, Mussig AJ, Hugenholtz P, Parks DH. GTDB-Tk: a toolkit to classify genomes with the Genome Taxonomy Database. Bioinformatics. 2019;36(6):1925–7; doi: 10.1093/bioinformatics/btz848.
Hyatt D, Chen GL, Locascio PF, Land ML, Larimer FW, Hauser LJ. Prodigal: prokaryotic gene recognition and translation initiation site identification. BMC Bioinformatics. 2010;11:119; doi: 10.1186/1471-2105-11-119.
Huerta-Cepas J, Szklarczyk D, Heller D, Hernandez-Plaza A, Forslund SK, Cook H, et al. eggNOG 5.0: a hierarchical, functionally and phylogenetically annotated orthology resource based on 5090 organisms and 2502 viruses. Nucleic Acids Res. 2019;47(D1):D309-D14; doi: 10.1093/nar/gky1085.
Bateman A, Coin L, Durbin R, Finn RD, Hollich V, Griffiths-Jones S, et al. The Pfam protein families database. Nucleic Acids Res. 2004;32(Database issue):D138-41; doi: 10.1093/nar/gkh121.
Tatusov RL, Natale DA, Garkavtsev IV, Tatusova TA, Shankavaram UT, Rao BS, et al. The COG database: new developments in phylogenetic classification of proteins from complete genomes. Nucleic Acids Res. 2001;29(1):22–8; doi: 10.1093/nar/29.1.22.
Haft DH, Selengut JD, Richter RA, Harkins D, Basu MK, Beck E. TIGRFAMs and Genome Properties in 2013. Nucleic Acids Res. 2013;41(Database issue):D387-95; doi: 10.1093/nar/gks1234.
Kanehisa M, Sato Y, Morishima K. BlastKOALA and GhostKOALA: KEGG Tools for Functional Characterization of Genome and Metagenome Sequences. J Mol Biol. 2016;428(4):726–31; doi: 10.1016/j.jmb.2015.11.006.
Medema MH, Blin K, Cimermancic P, de Jager V, Zakrzewski P, Fischbach MA, et al. antiSMASH: rapid identification, annotation and analysis of secondary metabolite biosynthesis gene clusters in bacterial and fungal genome sequences. Nucleic Acids Res. 2011;39(Web Server issue):W339-46; doi: 10.1093/nar/gkr466.
Zhang H, Yohe T, Huang L, Entwistle S, Wu P, Yang Z, et al. dbCAN2: a meta server for automated carbohydrate-active enzyme annotation. Nucleic Acids Res. 2018;46(W1):W95-W101; doi: 10.1093/nar/gky418.
Bushnell B. BBMap: A Fast, Accurate, Splice-Aware Aligner. United States: N p, Web. 2014.
Rodriguez RL, Tsementzi D, Luo C, Konstantinidis KT. Iterative subtractive binning of freshwater chronoseries metagenomes identifies over 400 novel species and their ecologic preferences. Environ Microbiol. 2020;22(8):3394–412; doi: 10.1111/1462-2920.15112.
Nayfach S, Pollard KS. Average genome size estimation improves comparative metagenomics and sheds light on the functional ecology of the human microbiome. Genome Biol. 2015;16(1):51; doi: 10.1186/s13059-015-0611-7.
Ulyantsev VI, Kazakov SV, Dubinkina VB, Tyakht AV, Alexeev DG. MetaFast: fast reference-free graph-based comparison of shotgun metagenomic data. Bioinformatics. 2016;32(18):2760–7.
Sunagawa S, Coelho LP, Chaffron S, Kultima JR, Labadie K, Salazar G, et al. Ocean plankton. Structure and function of the global ocean microbiome. Science. 2015;348(6237):1261359; doi: 10.1126/science.1261359.
Haro-Moreno JM, Lopez-Perez M, de la Torre JR, Picazo A, Camacho A, Rodriguez-Valera F. Fine metagenomic profile of the Mediterranean stratified and mixed water columns revealed by assembly and recruitment. Microbiome. 2018;6(1):128.
Roux S, Enault F, Hurwitz BL, Sullivan MB. VirSorter: mining viral signal from microbial genomic data. PeerJ. 2015;3:e985; doi: 10.7717/peerj.985.
Arndt D, Grant JR, Marcu A, Sajed T, Pon A, Liang Y, et al. PHASTER: a better, faster version of the PHAST phage search tool. Nucleic Acids Res. 2016;44(W1):W16-21; doi: 10.1093/nar/gkw387.
Fu L, Niu B, Zhu Z, Wu S, Li W. CD-HIT: accelerated for clustering the next-generation sequencing data. Bioinformatics. 2012;28(23):3150–2; doi: 10.1093/bioinformatics/bts565.
Shannon P, Markiel A, Ozier O, Baliga NS, Wang JT, Ramage D, et al. Cytoscape: a software environment for integrated models of biomolecular interaction networks. Genome Res. 2003;13(11):2498–504.
Thornhill DJ, Howells EJ, Wham DC, Steury TD, Santos SR. Population genetics of reef coral endosymbionts (Symbiodinium, Dinophyceae). Mol Ecol. 2017;26(10):2640–59; doi: 10.1111/mec.14055.
Kramarsky-Winter EH, M.; Sibon, N. ;Ben Dov, E. ; Brickner, I.; Loya, Y.; Kushmaro, A. Identification of a protist-coral association and its possible ecological role. Marine Ecology Progress Series. 2006;317:67–73; doi: 10.3354/meps317067.
Kohler S, Delwiche CF, Denny PW, Tilney LG, Webster P, Wilson RJ, et al. A plastid of probable green algal origin in Apicomplexan parasites. Science. 1997;275(5305):1485–9; doi: 10.1126/science.275.5305.1485.
Janouskovec J, Horak A, Barott KL, Rohwer FL, Keeling PJ. Global analysis of plastid diversity reveals apicomplexan-related lineages in coral reefs. Curr Biol. 2012;22(13):R518-9.
Janouskovec J, Horak A, Barott KL, Rohwer FL, Keeling PJ. Environmental distribution of coral-associated relatives of apicomplexan parasites. Isme J. 2013;7(2):444–7.
Bentis CJ, Kaufman L, Golubic S. Endolithic fungi in reef-building corals (Order: Scleractinia) are common, cosmopolitan, and potentially pathogenic. Biol Bull. 2000;198(2):254–60.
Raghukumar C, Ravindran J. Fungi and their role in corals and coral reef ecosystems. Prog Mol Subcell Biol. 2012;53:89–113.
Litchman E. Invisible invaders: non-pathogenic invasive microbes in aquatic and terrestrial ecosystems. Ecol Lett. 2010;13(12):1560–72; doi: 10.1111/j.1461-0248.2010.01544.x.
Yarden O, Ainsworth TD, Roff G, Leggat W, Fine M, Hoegh-Guldberg O. Increased prevalence of ubiquitous ascomycetes in an acropoid coral (Acropora formosa) exhibiting symptoms of Brown Band syndrome and skeletal eroding band disease. Appl Environ Microbiol. 2007;73(8):2755–7; doi: 10.1128/AEM.02738-06.
Savary R, Masclaux FG, Wyss T, Droh G, Cruz Corella J, Machado AP, et al. A population genomics approach shows widespread geographical distribution of cryptic genomic forms of the symbiotic fungus Rhizophagus irregularis. Isme J. 2018;12(1):17–30.
Siboni N, Ben-Dov E, Sivan A, Kushmaro A. Global distribution and diversity of coral-associated Archaea and their possible role in the coral holobiont nitrogen cycle. Environ Microbiol. 2008;10(11):2979–90; doi: 10.1111/j.1462-2920.2008.01718.x.
Rubio-Portillo E, Santos F, Martinez-Garcia M, de Los Rios A, Ascaso C, Souza-Egipsy V, et al. Structure and temporal dynamics of the bacterial communities associated to microhabitats of the coral Oculina patagonica. Environ Microbiol. 2016;18(12):4564–78; doi: 10.1111/1462-2920.13548.
McCutcheon JP, Moran NA. Extreme genome reduction in symbiotic bacteria. Nat Rev Microbiol. 2011;10(1):13–26; doi: 10.1038/nrmicro2670.
Hochart C, Paoli L, Ruscheweyh HJ, Salazar G, Boissin E, Romac S, et al. Ecology of Endozoicomonadaceae in three coral genera across the Pacific Ocean. Nat Commun. 2023;14(1):3037; doi: 10.1038/s41467-023-38502-9.
Tang YZ, Koch F, Gobler CJ. Most harmful algal bloom species are vitamin B1 and B12 auxotrophs. Proc Natl Acad Sci U S A. 2010;107(48):20756–61; doi: 10.1073/pnas.1009566107.
Al-Khodor S, Price CT, Kalia A, Abu Kwaik Y. Functional diversity of ankyrin repeats in microbial proteins. Trends Microbiol. 2010;18(3):132–9; doi: 10.1016/j.tim.2009.11.004.
Jernigan KK, Bordenstein SR. Ankyrin domains across the Tree of Life. PeerJ. 2014;2:e264; doi: 10.7717/peerj.264.
Lesser MP, Falcon, L. I., Rodriguez-Roman, A., Enriquez, S., Hoegh-Guldberg, O., and Iglesias-Prieto, R. Nitrogen fixation by symbiotic cyanobacteria provides a source of nitrogen for the scleractinian coral Montastraea cavernosa. Marine Ecology-progress Series. 2007;346:143–52; doi: doi.org/10.3354/meps07008.
Beaulaurier J, Luo E, Eppley JM, Uyl PD, Dai X, Burger A, et al. Assembly-free single-molecule sequencing recovers complete virus genomes from natural microbial communities. Genome Res. 2020;30(3):437–46; doi: 10.1101/gr.251686.119.
Correa AM, Ainsworth TD, Rosales SM, Thurber AR, Butler CR, Vega Thurber RL. Viral Outbreak in Corals Associated with an In Situ Bleaching Event: Atypical Herpes-Like Viruses and a New Megavirus Infecting Symbiodinium. Front Microbiol. 2021;7:127.
Lawrence SA, Wilson WH, Davy JE, Davy SK. Latent virus-like infections are present in a diverse range of Symbiodinium spp. (Dinophyta). J Phycol. 2014;50(6):984–97; doi: 10.1111/jpy.12242.
Correa AM, Welsh RM, Vega Thurber RL. Unique nucleocytoplasmic dsDNA and + ssRNA viruses are associated with the dinoflagellate endosymbionts of corals. ISME J. 2013;7(1):13–27; doi: 10.1038/ismej.2012.75.
Haro-Moreno JM, Rodriguez-Valera F, Lopez-Perez M. Prokaryotic Population Dynamics and Viral Predation in a Marine Succession Experiment Using Metagenomics. Front Microbiol. 2019;10:2926; doi: 10.3389/fmicb.2019.02926.
Wright RM, Strader ME, Genuise HM, Matz M. Effects of thermal stress on amount, composition, and antibacterial properties of coral mucus. PeerJ. 2019;7:e6849.
Xu JH, Lai KH, Su YD, Chang YC, Peng BR, Backlund A, et al. Briaviolides K-N, New Briarane-Type Diterpenoids from Cultured Octocoral Briareum violaceum. Mar Drugs. 2018;16(3); doi: 10.3390/md16030075.
Abd El-Rahman TMA, Tharwat NA, Abo El-Souad SMS, El-Beih AA, El-Diwany AI. Biological activities and variation of symbiotic fungi isolated from Coral reefs collected from Red Sea in Egypt. Mycology. 2020;11(3):243–55; doi: 10.1080/21501203.2020.1741470.
Pootakham W, Mhuantong W, Yoocha T, Putchim L, Jomchai N, Sonthirod C, et al. Heat-induced shift in coral microbiome reveals several members of the Rhodobacteraceae family as indicator species for thermal stress in Porites lutea. Microbiologyopen. 2019;8(12):e935; doi: 10.1002/mbo3.935.
Raina JB, Tapiolas DM, Foret S, Lutz A, Abrego D, Ceh J, et al. DMSP biosynthesis by an animal and its role in coral thermal stress response. Nature. 2013;502(7473):677–80; doi: 10.1038/nature12677.
Gardner SG, Raina JB, Ralph PJ, Petrou K. Reactive oxygen species (ROS) and dimethylated sulphur compounds in coral explants under acute thermal stress. J Exp Biol. 2017;220(Pt 10):1787–91; doi: 10.1242/jeb.153049.
Gardner SN, MRR; O'Brien, J.; Motti, CAA; Seymour, JRR; Ralph, PJJ; Petrou, K; Raina, J-B. IncreasedDMSPavailability duringthermalstressinfluences DMSP-degradingbacteriain coralmucus. Frontiers in Marine Science. 2022;9; doi: 10.3389/fmars.2022.912862.
Raina JB, Dinsdale EA, Willis BL, Bourne DG. Do the organic sulfur compounds DMSP and DMS drive coral microbial associations? Trends Microbiol. 2010;18(3):101–8; doi: 10.1016/j.tim.2009.12.002.
Zhang Y, Yang Q, Zhang Y, Ahmad M, Ling J, Tang X, et al. Shifts in abundance and network complexity of coral bacteria in response to elevated ammonium stress. Sci Total Environ. 2021;768:144631; doi: 10.1016/j.scitotenv.2020.144631.
Brownell A.C.; Richardson LL. Sulfate reducing bacteria as secondary and necessary pathogens in black band disease of corals. Rev biol trop. 2014;62.
Aranda M, Li Y, Liew YJ, Baumgarten S, Simakov O, Wilson MC, et al. Genomes of coral dinoflagellate symbionts highlight evolutionary adaptations conducive to a symbiotic lifestyle. Sci Rep. 2016;6:39734; doi: 10.1038/srep39734.
Nakajima Y, Nishikawa A, Isomura N, Iguchi A, Sakai K. Genetic connectivity in the broadcast-spawning coral Acropora digitifera analyzed by microsatellite markers on the Sekisei Reef, southwestern Japan. Zoolog Sci. 2009;26(3):209–15; doi: 10.2108/zsj.26.209.
Vega Thurber RL, Barott KL, Hall D, Liu H, Rodriguez-Mueller B, Desnues C, et al. Metagenomic analysis indicates that stressors induce production of herpes-like viruses in the coral Porites compressa. Proc Natl Acad Sci U S A. 2008;105(47):18413–8; doi: 10.1073/pnas.0808985105.
Wang W, Tang K, Wang P, Zeng Z, Xu T, Zhan W, et al. The coral pathogen Vibrio coralliilyticus kills non-pathogenic holobiont competitors by triggering prophage induction. Nat Ecol Evol. 2022;6(8):1132–44; doi: 10.1038/s41559-022-01795-y.
Mills S, Shanahan F, Stanton C, Hill C, Coffey A, Ross RP. Movers and shakers: influence of bacteriophages in shaping the mammalian gut microbiota. Gut Microbes. 2013;4(1):4–16; doi: 10.4161/gmic.22371.

No competing interests reported.

SupplementaryFigure1.pdf
Supplementary Figure 1. Whole metagenome comparison of O. patagonica samples compared to other coral metagenomes and seawater samples using MASH.
SupplementaryFigure2.pdf
Supplementary Figure 2. CAZymes detected (Carbon bidding modules, Glycosylhydrolases and Glycosyltransferases) using dbCAN in the O. patagonica samples. C: Control Corals; I: Infected Corals; MS: Corals in the Presence of Vibrio Monocultures; MX: Corals in the Presence of Vibrio Coculture. The numerical annotations 20 and 28 refer to corals maintained at 20°C and 28°C, respectively.
SupplementaryFigure3.pdf
Supplementary Figure 3. Number of reads recruited against the contig named prophages-C28-NODE_25_length_132664_cov_7.39709 and identified as Phycodnaviridae virus using subsets of 1 million reads from the O. patagonica metagenomes. C: Control Corals; I: Infected Corals; MS: Corals in the Presence of Vibrio Monocultures; MX: Corals in the Presence of Vibrio Coculture. The numerical annotations 20 and 28 refer to corals maintained at 20°C and 28°C, respectively.
SupplementaryTable1.Metagenomefeatures.xlsx
Supplementary Table 1. Characteristics of the viral metagenomes sequenced in this work.
SupplementaryTable2.ReadstaxonomyusingKaiju.xlsx
Supplementary Table 2. Number of phyla and orders in the metagenomes as determined using Kaiju.
SupplementaryTable3.MAGSgenomicfeatures.xlsx
Supplementary Table 3. Genomic features of the MAGs (completeness, 80%; contamination, 5%) assembled from the metagenome of O. patagonica samples.
SupplementaryTable4.ANIrfromBINS.xlsx
Supplementary Table 4 ANIr from the MAGs obtained in this work. Bolded text indicates the dataset from which the MAG was assembled.
SupplementaryTable5.KEGGmodulesMAGs.xlsx
Supplementary Table 5. KEGG modules from the MAGs assembled from the metagenomes obtained in this work. Number of genes in each module is indicated (RED= complete module, PINK= one block missing).
SupplementaryTable6.Kaijuviraltaxonomy.xlsx
Supplementary Table 6. Viral taxonomy of vOTUs as determined using Kaiju.

Download PDF

Journal Publication

published 29 Oct, 2024

Read the published version in Microbiome →

Editorial decision: Revision requested
06 Feb, 2024
Editor assigned by journal
06 Feb, 2024
Submission checks completed at journal
24 Jan, 2024
First submitted to journal
24 Jan, 2024

You are reading this latest preprint version

Community shifts in the coral Oculina patagonica holobiont in response to confinement, temperature and vibrio infections

Status:

Journal Publication

Version 1

Abstract

Figures

INTRODUCTION

MATERIAL AND METHODS

RESULTS AND DISCUSSION

Changes due to confinement

Changes due to thermal stress under experimental conditions

CONCLUSIONS

Declarations

References

Additional Declarations

Supplementary Files

Status:

Journal Publication

Version 1