Thermal regime and host clade, rather than geography, drive Symbiodinium and bacterial assemblages in the scleractinian coral Pocillopora damicornis sensu lato
© The Author(s). 2018
Received: 19 July 2017
Accepted: 11 February 2018
Published: 20 February 2018
Although the term holobiont has been popularized in corals with the advent of the hologenome theory of evolution, the underlying concepts are still a matter of debate. Indeed, the relative contribution of host and environment and especially thermal regime in shaping the microbial communities should be examined carefully to evaluate the potential role of symbionts for holobiont adaptation in the context of global changes. We used the sessile, long-lived, symbiotic and environmentally sensitive reef-building coral Pocillopora damicornis to address these issues.
We sampled Pocillopora damicornis colonies corresponding to two different mitochondrial lineages in different geographic areas displaying different thermal regimes: Djibouti, French Polynesia, New Caledonia, and Taiwan. The community composition of bacteria and the algal endosymbiont Symbiodinium were characterized using high-throughput sequencing of 16S rRNA gene and internal transcribed spacer, ITS2, respectively. Bacterial microbiota was very diverse with high prevalence of Endozoicomonas, Arcobacter, and Acinetobacter in all samples. While Symbiodinium sub-clade C1 was dominant in Taiwan and New Caledonia, D1 was dominant in Djibouti and French Polynesia. Moreover, we also identified a high background diversity (i.e., with proportions < 1%) of A1, C3, C15, and G Symbiodinum sub-clades. Using redundancy analyses, we found that the effect of geography was very low for both communities and that host genotypes and temperatures differently influenced Symbiodinium and bacterial microbiota. Indeed, while the constraint of host haplotype was higher than temperatures on bacterial composition, we showed for the first time a strong relationship between the composition of Symbiodinium communities and minimal sea surface temperatures.
Because Symbiodinium assemblages are more constrained by the thermal regime than bacterial communities, we propose that their contribution to adaptive capacities of the holobiont to temperature changes might be higher than the influence of bacterial microbiota. Moreover, the link between Symbiodinium community composition and minimal temperatures suggests low relative fitness of clade D at lower temperatures. This observation is particularly relevant in the context of climate change, since corals will face increasing temperatures as well as much frequent abnormal cold episodes in some areas of the world.
All partners (bionts) involved in a stable symbiosis, and thus being part of the entire organism, constitute the holobiont . A decade after this term has been defined; it has been popularized in corals  and subsequently led to the hologenome theory of evolution [3–6]. The hologenome is defıned as the sum of the genetic information of the host and its symbiotic microorganisms. In this context, phenotypes are the product of the collective genomes of the holobiont partners, being the true unit of biological organization and thus the object of natural selection [7–9]. This concept has gained increased attention for many issues on the functioning, homeostasis, or evolution of living organisms, extending our knowledge of microbial community associated to them (see  for a review on metaorganisms).
Scleractinian corals (the major reef-building organisms of coral reef ecosystems) are considered as the most diverse symbiotic ecosystem studied to date , forming a complex consortium composed by the cnidarian host, as well as microbial eukaryotes (including the dinoflagellate endosymbiont Symbiodinium), prokaryotes (bacteria and archae), and viruses. The symbiosis between corals and dinoflagellate algae of the genus Symbiodinium provides the foundation for the ecological success of coral reefs over millions of years . In this phototrophic and potentially mutualistic association, the coral host provides inorganic nutrients in exchange of photosynthetically fixed carbon (photosynthates) and amino acids from the algal symbiont [12–14]. Algae from the Symbiodinium genus are classified into nine clades (from A to I) [15, 16], and the physiology of the coral holobiont is affected by the clade of the symbiont [17, 18]. Although the establishment of a specific symbiosis occurs during early stages of host larvae colonization , some coral species can switch endosymbiotic algae during their lifetime in response to major environmental stressors such as thermally induced bleaching events [20, 21], although the long-term persistence of such changes is matter of debate . Bacterial communities associated to corals have also been extensively studied comparing different species [2, 23], disease states [24–26], or environmental conditions [27, 28] (see also  for a review).
Despite this large corpus of studies, none addressed the effect of natural thermal regimes on microbial assemblages in scleractinian corals. However, sea surface temperature increase is the main factor of ongoing climate changes affecting reef-building corals (32.8% of species being considered at risk of extinction)  with mass mortalities following severe and recurrent bleaching events .
In this study, we investigated the effect of thermal regimes, as well as host clade and geographical distribution, on the bacterial and Symbiodinium assemblages in the complex Pocillopora damicornis sensus lato (Veron and Pichon 1976), a functional group of environmentally sensitive scleractinian corals  that was recently split into five clades [32–34]. High-throughput metabarcoding allowed us to access for the first time the unculturable diversity of both bacterial communities and Symbiodinium assemblages in coral populations from different locations with contrasting thermal regimes.
Sampling sites and study design
Details of the sampling of Pocillopora damicornis sensu lato conducted in this study
11° 43′ 49.27″ N
43° 13′ 26.69″ E
11° 46′ 46.15″ N
42° 55′ 27.65″ E
11° 34′ 57.41″ N
42° 47′ 45.86″ E
17° 29′ 31.72″ S
149° 52′ 8.34″ O
17° 29′ 22.90″ S
149° 53′ 48.65″ O
16° 47′ 22.00″ S
151° 23′ 30.63″ O
16° 36′ 51.03″ S
151° 32′ 33.91″ O
16° 40′ 36.26″ S
151° 27′ 18.92″ O
17° 34′ 27.93″ S
149° 37′ 11.07″ O
Baie de Ste. Marie
22° 17′ 35.86″ S
166° 28′ 14.44″ E
Baie de Ma’a
22° 11′ 45.17″ S
166° 20′ 23.67″ E
22° 14′ 38.63″ S
166° 25′ 7.58″ E
21° 55′ 48.4″ N
120° 44′ 41.9″ E
21° 56′ 18.6″ N
120° 44′ 45.7″ E
21° 59′ 37.7″ N
120° 42′ 22.7″ E
Thermal regime descriptors were compiled from weekly mean sea surface temperature data collected from IGOSS, Integrated Global Ocean Services System Products Bulletin (http://iridl.ldeo.columbia.edu/SOURCES/.IGOSS/) for quadrats of 1° longitude × 1° latitude and from 1982 to the year of sampling (2008–2009).
DNA extraction was performed using CTAB (cetyl trimethylammonium bromide)-based extraction method . Briefly, coral tips were lysed 2 h in 600 μL CTAB buffer (2% CTAB, 0.2% β-mercaptoethanol, 1.4 M NaCI, 20 mm EDTA pH 8, I00 mm Tris-HCl pH 8, 100 μg/mL proteinase K) at 60 °C. An equal volume of chloroform:isoamyl-alcohol (24:1) was then added. After centrifugation for 10 min at 14,000g, the aqueous phase was transferred to a new tube and DNA was precipitated by adding 400 μL of ice-cold isopropanol and incubating 1 h at − 20 °C. After centrifugation for 15 min at 14,000g, the supernatant was discarded and the pellet washed with 70% ethanol. The pellet was air-dried for 5 min and resuspended in water.
Host haplotype identification
The mitochondrial variable open reading frame (ORF) was amplified with FATP6.1 (5′-TTTGGGSATTCGTTTAGCAG-3′) and RORF (5′-SCCAATATGTTAAACASCATGTCA-3′) primers  and submitted to sanger sequencing (GenBank submission ID 2033564). Protein-coding sequences were analyzed using MEGA version 6 . Sequence alignment was performed using MUSCLE. The best model (Kimura-2 parameters with uniform substitution rates) was selected for the lowest BIC (Bayesian Information Criterion). Maximum-likelihood tree was computed with the best model, and the robustness of the tree was tested with 1000 bootstrap replicates.
16S rRNA gene and internal transcribed spacer ITS2 metabarcoding
For all 94 individual samples, amplicon sequencing was performed for the two markers separately. For bacterial communities, 16S rRNA gene libraries were generated using the 341F (5′-CCTACGGGNGGCWGCAG-3′) and 805R (5′-GACTACHVGGGTATCTAATCC-3′) primers targeting the variable V3V4 loops . For each sample, a first PCR amplification of 35 cycles was performed for the generation of amplicons (96 °C for 30 s, 55 °C for 30 s, 72 °C for 1 min) followed by 15 cycles with Illumina adaptors and sequencing index for thelibrary (95 °C for 15 s, 58 °C for 30s, 72 °C for 1 min).
Symbiodinium assemblages were analyzed using ITS2 (internal transcribed spacer of the ribosomal RNA gene) amplicon libraries with specific primers targeting a ~ 350 bp sequence: ITS2-F (5′-GTGAATTGCAGAACTCCGTG-3′) and ITS2-R (5′-CCTCCGCTTACTTATATGCTT-3′) [41, 42]. For each sample, a first PCR amplification of 35 cycles was performed for the generation of amplicons (96 °C for 30 s, 54 °C for 30 s, 72 °C for 1 min) followed by 15 cycles with Illumina adaptors and sequencing index for the library (95 °C for 15 s, 58 °C for 30 s, 72 °C for 1 min).
For both 16S rRNA gene and ITS2 markers, paired-end sequencing with 250 bp read length was performed on two different flow cells using the MiSeq system (Illumina) using the v2 chemistry according to the manufacturer’s protocol. Sequencing was performed at the McGill University in the Génome Québec Innovation Centre, Montréal, Canada. Raw sequence data are available in the SRA database (BioProject ID PRJNA393088, BioSample ID SUB2829485).
Sequence analysis of metabarcoding datasets
The FROGS pipeline (Find Rapidly OTU with Galaxy Solution) implemented on a galaxy instance (http://sigenae-workbench.toulouse.inra.fr/galaxy/) was used for data processing . In brief, paired reads were merged using FLASH . After denoising and primer/adapters removal with cutadapt , de novo clustering was done using SWARM that uses a local clustering threshold , with aggregation distance d = 3 after denoising. Chimera were removed using VSEARCH . We filtered the dataset for singletons and performed affiliation using Blast+ against the Silva database (release 128, September 2016) for 16S rRNA gene amplicons. For ITS2 metabarcoding, the Symbiodinium clade was assessed using blastn best hit against the nr/nt database of the NCBI . To confirm blast identification, we performed phylogenetic analysis using MAFFT to produce sequence alignment and FastTree (GTR + CAT model) to compute tree with the approximately maximum-likelihood method. Finally, OTU tables were produced in a standard BIOM format for subsequent analyses.
All statistical analyses were done using R version 3.3.1 (, http://www.R-project.org). We used the phyloseq R package for community composition analysis  to infer alpha diversity metrics at the OTU level, as well as beta diversity (between sample distance) from the OTU table. Community similarity was assessed by principal coordinate analysis (PCoA) using the Bray-Curtis dissimilarity index.
We performed non-parametric Kruskal-Wallis tests (since normality of residuals and homogeneity of variances were rejected (Shapiro and Bartlett test, respectively)) to compare alpha diversity metrics (Chao1 and Shannon). When Kruskal-Wallis tests were significant, we then computed pairwise comparisons between group levels (post hoc analyses) with Bonferroni corrections for multiple testing using Dunn tests.
Redundancy analysis  (hereafter named RDA) was used to investigate the variations of the different bacterial OTU/Symbiodinium clades under the constraint of the environmental variables. Bacterial OTUs and Symbiodinium clades were Hellinger-transformed before performing RDAs  on datasets of geography, temperature, and host genotypes (Additional file 1: Table S1). Then significant variables (i.e., variables that significantly explained changes in the distribution of OTUs/clades) were identified using a forward-selection procedure (999 permutations), implemented in the R “vegan” package and in the “rda” and “ordiR2step” functions. For all analyses, the threshold significance level was set at 0.05.
Sampling sites display diverse and contrasted thermal regimes
Thermal regime descriptors of the four sampled regions calculated from weekly SST data (from 1982 to the year of sampling, 2008 or 2009, according to regions)
Mean temperature (°C)
Max temperature (°C)
Min temperature (°C)
Mean of the warmer months (°C)
Mean of the colder months (°C)
Pocillopora damicornis sensu lato morphotypes corresponded to two different clades
Analysis of mitochondrial ORF for all samples revealed two clades corresponding to Pocillopora types 3c, 3e, 3g, 3h, and 7a (clade 2) for all samples from Djibouti and Pocillopora type 5 (clade 1) for samples from New Caledonia, Taiwan, and French Polynesia  (Additional file 2: Figure S1 and Additional file 1: Table S1). This is consistent with the known geographical distribution of Pocillopora lineages with a wide distribution of clade 1, i.e., types 4 and 5, namely Pocillopora damicornis and Pocillopora acuta, respectively, throughout the Pacific Ocean, and the presence of types 3 and 7 (P. damicornis sensu lato morphotype) in the Western Indian Ocean [32, 34].
Symbiodinium assemblages revealed high specificity as well as high background diversity
The relative proportions of Symbiodinium types in each coral sample were analyzed using high-throughput sequencing method. ITS2 amplicon sequencing with MiSeq yielded a total of 4,076,779 informative clusters. After clustering and filtering for OTUs containing less than 100 sequence tags, 53 clusters of ITS2 were obtained, corresponding to 3,990,373 sequences (Additional file 3: Table S2). Taxonomic affiliation was performed using blastn in comparison to NT followed by phylogenetic analysis. Maximum-likelihood trees confirmed the presence of clades A1, C1, C3, C15, D1, and G with largely unresolved polytomies in the clade C (Additional file 4: Figure S2), such as observed previously . Nevertheless, the topology is congruent with previous works on the molecular taxonomy of Symbiodinium .
Bacterial microbiota was highly diverse with a small number of core taxa
Phylum-level assignment of bacterial OTUs indicated the dominance of Proteobacteria in all samples (Fig. 4b). Within this phylum, the majority of sequences were affiliated to the genus Endozoicomonas (family Hahellaceae) representing 66.9% of the overall sequences (Additional file 7: Table S5). Three other genera, Arcobacter, Acinetobacter, and Sphingomonas were present at appreciable relative proportions (6.6, 1.9, and 1%, respectively, of total phyla tags). Arcobacter (Campylobacteraceae) was particularly abundant in samples from Djibouti (25% of sequence tags), whereas Acinetobacter (Moraxellaceae) represented 9.4% of the sequence tags in New Caledonia.
Number of samples for each region for core taxa (shared by at least 50% of samples within each region and 50% of the overall samples) at the family and genus levels
Considering each region independently, strikingly much more taxa were shared between at least 50% of samples from Djibouti (27 families and 46 genera) (Additional file 7: Table S5), which could be linked to higher equitability. Twelve families and 5 genera were core to French Polynesia samples, 12 families and 4 genera for New Caledonia, whereas only 7 families and 3 genera were shared by at least 50% of samples from Taiwan (Additional file 7: Table S5).
Symbiodinium communities were mostly influenced by minimal temperatures, whereas bacterial community structure was associated to host haplotype and mean temperatures
Diversity of Pocillopora damicornis sensu lato samples
In this study, we sampled Pocillopora damicornis sensu lato colonies in four regions and phylogenetic analyses revealed a complex pattern. Samples from French Polynesia, New Caledonia, and Taiwan corresponded to Pocillopora clade 1 (more precisely type 5, newly assigned to P. acuta sensu) [33, 34]. Samples from Djibouti corresponded to Pocillopora clade 2 (types 3 and 7) later assigned to SSH12 and SSH13 and all consistent with P. damicornis morphotype . Morphological characters have been shown to be insufficient to discriminate species in the P. damicornis complex, and genetic markers are necessary to discriminate the species . This species complex displays in addition high phenotypic plasticity and cryptic lineages . Nevertheless, as we used a single mitochondrial marker for Pocillopora clade identification, we cannot completely rule out the possibility of mitochondrial introgression or hybridization that seem to be common within this genus [55, 56]. Finally, although it was proposed that Symbiodinium assemblages may be useful for integrative taxonomy of Pocillopora species , we found in this study that Symbiodinium communities did not discriminate Pocillopora haplotypes.
Diversity of coral microbiota
The core microbiome can be broadly defined as the stable and consistent taxonomic groups associated to a particular habitat . These commonly associated microbial communities are likely of ecological and functional importance in the holobiont fitness. The definition of the core microbiome in corals is variable among authors, ranging from 30 to 100% of shared bacterial phylotypes at all different taxonomic levels from kingdom to strain (reviewed in ). These studies already highlighted the high diversity and variability of microbial communities in corals. We defined here the core microbiome as taxa (OTUs/genera/families) present in at least 50% of samples within each population and thus reflecting stable associations with P. damicornis sensu lato.
All core taxa of bacteria were among the most abundant in the dataset, which is in contradiction with previous study identifying rare taxa as ubiquitous endosymbionts . Bacteria of the Endozoicomonas genus were present in all samples and were the most abundant in the whole dataset, making up from 33.3% of the total number of sequence tags in Djibouti (consistent with a higher evenness) to 87.4% in French Polynesia where rare OTUs were the most diverse. This group has been described for the first time only recently in marine slugs . So far, they have been identified in numerous coral species and are now considered as ubiquitous endosymbionts of many marine host species (see  for a review). Notably, comparative genomic analysis of different strains of Endozoicomonas from different hosts suggests a common role in carbohydrate cycling with potential specificities in amino-acid synthesis . In addition, Acinetobacter and Arcobacter were respectively present in 86 and 84 samples over 94 and were the third and second most abundant genus after Endozoicomonas. Although considered as a terrestrial bacteria, Acinetobacter (gamma-proteobacteria) can be dominant in bleached corals . It has also been identified in healthy corals where it can be abundant [64–66], but its function remains elusive. Bacteria of the genus Arcobacter (epsilon-proteobacteria) are associated to a wide range of habitats, as free-living or pathogenic, and are especially abundant in marine environments  and have been identified specifically in diseased corals in some studies , sometimes also in healthy samples , and has been show to increase during pathogenesis [70, 71]. Arcobacter has also been found in high abundances in necrotic sponges  as well as moribund oysters  where it is associated with decreased bacterial diversity and may act as an opportunistic pathogen. Although detectable at background abundances in the majority of samples, Arcobacter was highly abundant in all samples from Djibouti (25% of sequence tags vs. less than 1% in the other regions), indicating that this genus can also be abundant in healthy hosts and moreover associated with high bacterial diversity.
Metabarcoding targeting ITS2 provides highly sensitive measurements of Symbiodinium diversity and relative proportions . It was used with 454-pyrosequencing [74, 75] and more recently with Illumina MiSeq sequencing . Although we did not find ubiquitous clades for the whole samples, we showed using MiSeq sequencing striking differences in terms of Symbiodinium associations between populations. Clade D1 was dominant in French Polynesia. This is consistent with a recent study in Moorea that showed association with Symbiodinium clade D and few clade C 
Colonies from Djibouti were also mainly associated with clade D1, while clade C1 was dominant in Taiwan and New Caledonia. To our knowledge, this is the first study of symbiotic communities associated with Pocillopora in Djibouti, New Caledonia, and Taiwan. A comparative study of Symbiodinium clades associated with scleractinian corals from the Persian Gulf with annual temperatures from 16 up to 36 °C revealed the dominance of clade C3 associated with a lower proportion of A1 and minor quantities of C15, whereas clade D was not detected .
Even if one single clade was dominant in almost all samples, we detected background abundances of the other clades in all regions constituting the so called “Symbiodinium rare biosphere” . The cryptic diversity that we reported here is actually very high (H′ Shannon index up to 0.6 for French Polynesia samples, Fig. 2a) and falls into a similar range such as that reported in the latter work for three Pocilloporid species . Such high diversity can potentially confer higher adaptive potential to their hosts through redundant or complementary symbiotic functions and/or potential for symbiont switching or shuffling . High-throughput methods (ITS metabarcoding using 454) applied to the generalist and environmentally sensitive Acropora species also revealed the same pattern across Western Australia with clade C being dominant in all samples with background abundances of other clades .
Finally, bacterial communities in corals may also be shaped by the presence and type of photosymbionts  through the use of DMSP released by the algae , but at a finer scale, we could not find any relationship for the same host species between Symbiodinium type and bacterial communities as was also the case in Seriatopora hystrix . In addition, we observed contrasting patterns in terms of alpha diversity between the Symbiodinium and bacteria communities, with notably an opposite pattern on the Shannon index for samples from Djibouti (highest Shannon index for 16S, Fig. 2a) and French Polynesia (highest Shannon index for ITS2, Fig. 4a). Previous work in Isopora palifera at the seasonal scale also revealed important variations in Shannon diversity that are uncoupled from changes in Symbiodinium community composition .
Influence of geography on coral microbiota is very low
Because we found different microbiota composition between samples, we used redundancy analyses in order to identify variables that significantly constrained their composition. We first found that the influence of geography was low, since longitude had no significant effect and the effect of latitude in our dataset was null for Symbiodinium assemblages and was marginal on bacterial communities (less than 1% when added to the variance explained by host haplotype and mean temperature). In contradiction with this result, geographic clustering of symbionts is expected to be high in brooding species with vertical transmission, compared to those with horizontal transmission that would in contrary favor propagation of locally adapted symbionts. Pocillopora damicornis is a hermaphrodite brooder [83, 84] thus transmitting its symbionts vertically  and is able to release sexual as well as parthenogenetic larvae [86, 87]. However, these results are similar to others studies showing that bacterial communities associated with different coral species are stable through space and time . Thus, we hypothesize that low parental effect on symbiont communities may be explained by the relatively high connectivity and long-range dispersal for this species compared to other brooders [35, 88], which would mitigate the effect of geography and allow higher symbiont mixing with the local microbial environment.
Influence of host haplotype is higher for bacterial communities than for Symbiodinium
Although we did not find a significant link between microbiota and geography, we highlighted links with host haplotypes, particularly for bacterial microbiota. As the two host clades sampled in this study were from potentially different microbial environments, we cannot rule out a confounding effect between host haplotype and environmental communities present in the seawater at the time of larval settlement and bacterial microbiota acquisition. Nevertheless, others studies observed that host genotypes might influence bacterial microbiota, since (i) different coral species harbored different microbiota composition  and (ii) differential specificity between closely related corals and abundant Endozoicomonas endosymbionts were found across global scales .
In contrast, we found a significant but less important link between Symbiodinium and host haplotype. This observation might also be linked to the adaptive bleaching hypothesis, which proposes that coral bleaching is an adaptive mechanism through the modification of Symbiodinium composition, with acquisition of clades with higher thermal tolerances . For example, endosymbiotic flexibility (the ability of a coral species to associate with multiple Symbiodinium clades) was observed to be linked with resistance . In addition, Pocillopora is a generalist host, displaying a high intraspecific and interspecific flexibility in terms of Symbiodinium assemblages (mainly among clades A, C, and D) .
Influence of temperature is higher for Symbiodinium than for bacterial microbiota
Although host genotypes mostly constrained bacterial assemblages, we also found a significant link with mean temperatures. Interestingly, reciprocal transplants of Acropora corals at a reduced geographic scale revealed as well that their microbiome were different across thermally variable habitats and changed over time after transplantation . As a consequence, bacterial communities might be linked to heat tolerance of their coral hosts, but further studies of their dynamics during a thermal stress are needed to understand their contribution to the holobiont response.
Strikingly, redundancy analysis revealed that most variation in Symbiodinium communities was actually explained by annual minimal temperatures (70%), whereas we observed no significant correlation with mean, maximal, or annual fluctuations of temperatures. In Acropora, it has been shown experimentally that Symbiodinium type alters larvae settlement in a temperature-dependent fashion . In the latter work, lower temperatures and not higher temperatures adversely affected recruitment by reducing larval survival and settlement. Accordingly, we propose that host/Symbiodinium associations are stable in space at a regional scale and are actually more sensitive to a minimal, threshold temperature, than to annual variations or to maximal temperatures.
In particular, clade D is common in all areas colonized by scleractinian corals , and many studies have shown that it confers higher thermotolerance to its host through photoprotection [94, 95] and that it increases in proportion after bleaching events [20, 95, 96]. Clade D is thus commonly considered as stress resistant and potentially opportunistic for compromised corals facing stressful conditions . Nevertheless and in accordance to the link we observed with minimal temperatures, some studies have revealed that cold stress may be detrimental for coral harboring this symbiont [98, 99]. More recently, it has been shown experimentally that clade D1a suffered more photodamages at low temperatures than clade C3 . Thus, it seems that clade D may not display higher plasticity in terms of thermal range, but rather that it is only adapted to hot thermal environments. Because climate changes will not only be associated with global increase of sea surface temperatures, but also to extreme thermal events including cold episodes , further studies are required to understand the potential response of coral reef in the future, particularly for regions with highly fluctuating thermal regimes and low temperature thresholds.
The hologenome theory of evolution reboots elements of the Lamarckian evolution and has been thus a matter of much debate (reviewed in [7, 102]). Although the host genome follows a Mendelian framework, potential changes in microbial assemblages may be inherited to the next generation. Conversely, for the selection to operate on holobiont phenotypes at evolutionary scales, specificity of the interaction through coevolution between partners would be expected.
Using combined high-throughput barcoding approaches for both bacterial and Symbiodinium communities, we showed that variation in Symbiodinium composition is mostly explained by thermal regime, especially minimal temperatures, whereas bacterial communities are much less related to temperature modifications. In this context, we propose that Symbiodinium types might confer more enhanced adaptive capacities to temperature modifications than the bacterial microbiota. However, Symbiodinium clade D was known to be adapted to high temperature, and we found a negative relationship with low temperatures, suggesting low plasticity for this clade. Such low plasticity might limit the adaptive capabilities of a coral associated to clade D and living in highly variable thermal regime. However, a high background diversity of Symbiodinium was also observed, providing the potential for coral colonies to adapt or acclimatize to future environmental changes via symbiont shuffling (i.e., changes in the relative proportion of Symbiodinium types constituting the within host community).
Our study may thus contribute to new insights into the importance of microbial (i.e., bacterial and Symbiodinium) communities for holobiont functioning as well as the relative importance of host and environmental factors in shaping the interaction.
We would like to thank Lorenzo Bramanti, Pauline Bosserelle, and Tung-Yung Fan for their help in the coral sampling. We are grateful to the genotoul bioinformatics platform Toulouse Midi-Pyrenees and Sigenae group for providing help and computing resources through Galaxy instance http://sigenae-workbench.toulouse.inra.fr.
This project has been funded by the ADACNI programme of the French National Research Agency (ANR) (project no. ANR-12-ADAP-0016; http://adacni.imbe.fr).
Availability of data and materials
The datasets generated during the current study are available in the Sequence Read Archive repository under BioProject ID PRJNA393088, BioSample ID SUB2829485 (to be released upon publication).
ET, JVD, MA, DA, FB, and GM were involved in the study concept and design. JVD and MA were involved in the collection of samples. KBR, CC, FB, DA, MP, and ET were involved in the data acquisition and analysis. ET drafted the manuscript, and all authors contributed to critical revisions and approved the final manuscript.
Ethics approval and consent to participate
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
- Margulis L, Fester R. Symbiosis as a source of evolutionary innovation. Cambridge: MIT Press; 1991. pp 2–14.Google Scholar
- Rohwer F, Seguritan V, Azam F. Diversity and distribution of coral-associated bacteria. Marine Ecology Progress. 2002;243:1–10.View ArticleGoogle Scholar
- Hernandez-Agreda A, Gates RD, Ainsworth TD. Defining the core microbiome in corals’ microbial soup. Trends Microbiol. 2017;25(2):125–40.View ArticlePubMedGoogle Scholar
- Jefferson R. The Hologenome. Agriculture, environment and the developing world: a future of PCR. New York: Cold Spring Harbor; 1994.Google Scholar
- Rosenberg E, Koren O, Reshef L, Efrony R, Zilber-Rosenberg I. The role of microorganisms in coral health, disease and evolution. Nat Rev Micro. 2007;5:355–62.View ArticleGoogle Scholar
- Zilber-Rosenberg I, Rosenberg E. Role of microorganisms in the evolution of animals and plants: the hologenome theory of evolution. FEMS Microbiol Rev. 2008;32:723–35.View ArticlePubMedGoogle Scholar
- Bordenstein SR, Theis KR. Host biology in light of the microbiome: ten principles of holobionts and hologenomes. PLoS Biol. 2015;13:e1002226.View ArticlePubMedPubMed CentralGoogle Scholar
- Guerrero R, Margulis L, Berlanga M. Symbiogenesis: the holobiont as a unit of evolution. Int Microbiol. 2013;16:133–43.PubMedGoogle Scholar
- McFall-Ngai M, Hadfield MG, Bosch TC, Carey HV, Domazet-Lošo T, Douglas AE, Dubilier N, Eberl G, Fukami T, Gilbert SF. Animals in a bacterial world, a new imperative for the life sciences. Proc Natl Acad Sci U S A. 2013;110:3229–36.View ArticlePubMedPubMed CentralGoogle Scholar
- Bosch TCG, McFall-Ngai MJ. Metaorganisms as the new frontier. Zoology. 2011;114:185–90.View ArticlePubMedPubMed CentralGoogle Scholar
- Blackall LL, Wilson B, van Oppen MJH. Coral-the world's most diverse symbiotic ecosystem. Mol Ecol. 2015;24:5330–47.View ArticlePubMedGoogle Scholar
- Muscatine L, Porter JW. Reef corals: mutualistic symbioses adapted to nutrient-poor environments. Bioscience. 1977;27:454–60.View ArticleGoogle Scholar
- Davy SK, Allemand D, Weis VM. Cell biology of cnidarian-dinoflagellate symbiosis. Microbiol Mol Biol Rev. 2012;76:229–61.View ArticlePubMedPubMed CentralGoogle Scholar
- DeSalvo MK, Sunagawa S, Fisher PL, Voolstra CR, Iglesias-Prieto R, Medina M. Coral host transcriptomic states are correlated with Symbiodiniumgenotypes. Mol Ecol. 2010;19:1174–86.View ArticlePubMedGoogle Scholar
- Baker AC. Flexibility and specificity in coral-algal symbiosis: diversity, ecology, and biogeography of Symbiodinium. Annual Review of Ecology. 2003;34:661–89.View ArticleGoogle Scholar
- Coffroth MA, Santos SR. Genetic diversity of symbiotic dinoflagellates in the genus Symbiodinium. Protist. 2005;156:19–34.View ArticlePubMedGoogle Scholar
- Mieog JC, van Oppen M, Cantin NE, Stam WT. Real-time PCR reveals a high incidence of Symbiodinium clade D at low levels in four scleractinian corals across the Great Barrier Reef: implications for symbiont shuffling. Coral Reefs. 2007;26:449–57.View ArticleGoogle Scholar
- Rouzé H, Lecellier G, Saulnier D, Berteaux-Lecellier V. Symbiodinium clades A and D differentially predispose Acropora cythereato disease and vibrio spp. colonization. Ecol Evol. 2016;6:560–72.View ArticlePubMedPubMed CentralGoogle Scholar
- Weis V, Reynolds W, deBoer M, Krupp D. Host-symbiont specificity during onset of symbiosis between the dinoflagellates Symbiodinium spp. and planula larvae of the scleractinian coral Fungia scutaria. Coral Reefs. 2001;20:301–8.View ArticleGoogle Scholar
- Baker AC, Starger CJ, McClanahan TR, Glynn PW. Coral reefs: corals’ adaptive response to climate change. Nature. 2004;430:741.View ArticlePubMedGoogle Scholar
- Jones AM, Berkelmans R, van Oppen MJH, Mieog JC, Sinclair W. A community change in the algal endosymbionts of a scleractinian coral following a natural bleaching event: field evidence of acclimatization. Proc R Soc B-Biol Sci. 2008;275:1359–65.View ArticleGoogle Scholar
- Thornhill DJ, Lajeunesse TC, Kemp DW, Fitt WK. Multi-year, seasonal genotypic surveys of coral-algal symbioses reveal prevalent stability or post-bleaching reversion. Mar Biol. 2006;148:711–22.View ArticleGoogle Scholar
- Littman R, Singh Y, Willis BL, Ahmad J, Pfeffer C, Musarrat J, Bourne DG, Ehtesham NZ, Hasnain SE. 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:152–63.View ArticlePubMedGoogle Scholar
- Cróquer A, Bastidas C, Elliott A, Sweet M. Bacterial assemblages shifts from healthy to yellow band disease states in the dominant reef coral Montastraea faveolata. Environ Microbiol Rep. 2013;5:90–6.View ArticlePubMedGoogle Scholar
- Meyer JL, Rodgers JM, Dillard BA, Paul VJ, Teplitski M. Epimicrobiota associated with the decay and recovery of Orbicella corals exhibiting dark spot syndrome. Front Microbiol. 2016;7:893.View ArticlePubMedPubMed CentralGoogle Scholar
- Sato Y, Willis BL, Bourne DG. Successional changes in bacterial communities during the development of black band disease on the reef coral, Montipora hispida. ISME J. 2009;4:203–14.View ArticlePubMedGoogle Scholar
- Bourne D, Bourne D, Iida Y, Iida Y, Uthicke S, Uthicke S, Smith-Keune C, Smith-Keune C. Changes in coral-associated microbial communities during a bleaching event. ISME J. 2007;2:350–63.View ArticlePubMedGoogle Scholar
- Pantos O, Bongaerts P, Dennis PG, Tyson GW, Hoegh-Guldberg O. Habitat-specific environmental conditions primarily control the microbiomes of the coral Seriatopora hystrix. ISME J. 2015;9:1916–27.View ArticlePubMedPubMed CentralGoogle Scholar
- Carpenter KE, Abrar M, Aeby G, Aronson RB. One-third of reef-building corals face elevated extinction risk from climate change and local impacts. Science. 2008;321:560–3.View ArticlePubMedGoogle Scholar
- Hughes TP, Kerry JT, Álvarez-Noriega M, Álvarez-Romero JG, Anderson KD, Baird AH, Babcock RC, Beger M, Bellwood DR, Berkelmans R, Bridge TC, Butler IR, Byrne M, Cantin NE, Comeau S, Connolly SR, Cumming GS, Dalton SJ, Diaz-Pulido G, Eakin CM, Figueira WF, Gilmour JP, Harrison HB, Heron SF, Hoey AS, J-PA H, Hoogenboom MO, Kennedy EV, Kuo C-Y, Lough JM, et al. Global warming and recurrent mass bleaching of corals. Nature. 2017;543:373–7.View ArticlePubMedGoogle Scholar
- van Woesik R, Sakai K, Ganase A, Loya Y. Revisiting the winners and the losers a decade after coral bleaching. Mar Ecol Prog Ser. 2011;434:67–76.View ArticleGoogle Scholar
- Gélin P, Postaire B, Fauvelot C, Magalon H. Reevaluating species number, distribution and endemism of the coral genus Pocillopora Lamarck, 1816 using species delimitation methods and microsatellites. Mol Phylogenet Evol. 2017;109:430–46.View ArticlePubMedGoogle Scholar
- Pinzón JH, Sampayo E, Cox E, Chauka LJ, Chen CA, Voolstra CR, LaJeunesse TC. Blind to morphology: genetics identifies several widespread ecologically common species and few endemics among Indo-Pacific cauliflower corals (Pocillopora, Scleractinia). J Biogeogr. 2013;40:1595–608.View ArticleGoogle Scholar
- Schmidt-Roach S, Miller KJ, Lundgren P, Andreakis N. With eyes wide open: a revision of species within and closely related to the Pocillopora damicornisspecies complex (Scleractinia; Pocilloporidae) using morphology and genetics. Zool J Linnean Soc. 2014;170:1–33.View ArticleGoogle Scholar
- Adjeroud M, Guérécheau A, Vidal-Dupiol J, Flot J-F, Arnaud-Haond S, Bonhomme F. Genetic diversity, clonality and connectivity in the scleractinian coral Pocillopora damicornis: a multi-scale analysis in an insular, fragmented reef system. Mar Biol. 2013;161:531–41.View ArticleGoogle Scholar
- Flot JF. Vers une taxonomie moléculaire des coraux du genre Pocillopora. 2007.Google Scholar
- Winnepenninckx B, Backeljau T, De Wachter R. Extraction of high molecular weight DNA from molluscs. Trends Genet. 1993;9:407.View ArticlePubMedGoogle Scholar
- Flot J-F, Tillier S. The mitochondrial genome of Pocillopora (Cnidaria: Scleractinia) contains two variable regions: the putative D-loop and a novel ORF of unknown function. Gene. 2007;401:80–7.View ArticlePubMedGoogle Scholar
- Tamura K, Stecher G, Peterson D, Filipski A, Kumar S. MEGA6: molecular evolutionary genetics analysis version 6.0. Mol Biol Evol. 2013;30:2725–9.View ArticlePubMedPubMed CentralGoogle Scholar
- Klindworth A, Pruesse E, Schweer T, Peplies J, Quast C, Horn M, Glockner FO. Evaluation of general 16S ribosomal RNA gene PCR primers for classical and next-generation sequencing-based diversity studies. Nucleic Acids Res. 2012;41:e1.View ArticlePubMedPubMed CentralGoogle Scholar
- Lajeunesse TC, Trench RK. Biogeography of two species of Symbiodinium (Freudenthal) inhabiting the intertidal sea anemone Anthopleura elegantissima (Brandt). Biol Bull. 2000;199:126–34.View ArticlePubMedGoogle Scholar
- Quigley KM, Quigley KM, Davies SW, Davies SW, Kenkel CD, Kenkel CD, Willis BL, Willis BL, Matz MV, Matz MV, Bay LK, Bay LK. Deep-sequencing method for quantifying background abundances of Symbiodinium types: exploring the rare Symbiodinium biosphere in reef-building corals. PLoS One. 2014;9:e94297.View ArticlePubMedPubMed CentralGoogle Scholar
- Escudié F, Auer L, Bernard M, Mariadassou M, Cauquil L, Vidal K, Maman S, Hernandez-Raquet G, Combes S, Pascal G. FROGS: find, rapidly, OTUs with galaxy solution. Bioinformatics. 2017; https://doi.org/10.1093/bioinformatics/btx791.
- Magoč T, Salzberg SL. FLASH: fast length adjustment of short reads to improve genome assemblies. Bioinformatics. 2011;27:2957–63.View ArticlePubMedPubMed CentralGoogle Scholar
- Martin M. Cutadapt removes adapter sequences from high-throughput sequencing reads. EMBnetjournal. 2011;17:10.Google Scholar
- Mahé F, Rognes T, Quince C, De Vargas C, Dunthorn M. Swarm: robust and fast clustering method for amplicon-based studies. PeerJ. 2014;2:e593.View ArticlePubMedPubMed CentralGoogle Scholar
- Rognes T, Flouri T, Nichols B, Quince C, Mahé F. VSEARCH: a versatile open source tool for metagenomics. PeerJ. 2016;4:e2584.View ArticlePubMedPubMed CentralGoogle Scholar
- Altschul SF, Gish W, Miller W, Myers EW, Lipman DJ. Basic local alignment search tool. J Mol Biol. 1990;215:403–10.View ArticlePubMedGoogle Scholar
- Team RC. R Core team. 2013. R: a language and environment for statistical computing. R Foundation for Statistical Computing, Vienna. 2013.Google Scholar
- McMurdie PJ, Holmes S. Phyloseq: an R package for reproducible interactive analysis and graphics of microbiome census data. PLoS One. 2013;8:e61217.View ArticlePubMedPubMed CentralGoogle Scholar
- Ramette A. Multivariate analyses in microbial ecology. FEMS Microbiol Ecol. 2007;62:142–60.View ArticlePubMedPubMed CentralGoogle Scholar
- Legendre P, Gallagher ED. Ecologically meaningful transformations for ordination of species data. Oecologia. 2001;129:271–80.View ArticlePubMedGoogle Scholar
- Lajeunesse TC. “Species” radiations of symbiotic dinoflagellates in the Atlantic and Indo-Pacific since the Miocene–Pliocene transition. Mol Biol Evol. 2004;22:570–81.View ArticlePubMedGoogle Scholar
- Pochon X, Montoya-Burgos JI, Stadelmann B, Pawlowski J. Molecular phylogeny, evolutionary rates, and divergence timing of the symbiotic dinoflagellate genus Symbiodinium. Mol Phylogenet Evol. 2006;38:20–30.View ArticlePubMedGoogle Scholar
- Combosch DJ, Vollmer SV. Trans-Pacific RAD-Seq population genomics confirms introgressive hybridization in Eastern Pacific Pocillopora corals. Mol Phylogenet Evol. 2015;88:154–62.View ArticlePubMedGoogle Scholar
- Combosch DJ, Guzman HM, Schuhmacher H, Vollmer SV. Interspecific hybridization and restricted trans-Pacific gene flow in the tropical Eastern Pacific Pocillopora. Mol Ecol. 2008;17:1304–12.View ArticlePubMedGoogle Scholar
- Pinzón JH, LaJeunesse TC. Species delimitation of common reef corals in the genus Pocillopora using nucleotide sequence phylogenies, population genetics and symbiosis ecology. Mol Ecol. 2011;20:311–25.View ArticlePubMedGoogle Scholar
- Shade A, Handelsman J. Beyond the Venn diagram: the hunt for a core microbiome. Environ Microbiol. 2011;14:4–12.View ArticlePubMedGoogle Scholar
- D Ainsworth T, Krause L, Bridge T, Torda G, Raina J-B, Zakrzewski M, Gates RD, Padilla-Gamiño JL, Spalding HL, Smith C, Woolsey ES, Bourne DG, Bongaerts P, Hoegh-Guldberg O, Leggat W. The coral core microbiome identifies rare bacterial taxa as ubiquitous endosymbionts. ISME J. 2015;9:2261–74.View ArticlePubMedPubMed CentralGoogle Scholar
- Kurahashi M, Yokota A. Endozoicomonas elysicola gen. nov., sp. nov., a gamma-proteobacterium isolated from the sea slug Elysia ornata. Syst Appl Microbiol. 2007;30:202–6.View ArticlePubMedGoogle Scholar
- Neave MJ, Apprill A, Ferrier-Pages C, Voolstra CR. Diversity and function of prevalent symbiotic marine bacteria in the genus Endozoicomonas. Appl Microbiol Biotechnol. 2016;100:8315–24.View ArticlePubMedPubMed CentralGoogle Scholar
- Neave MJ, Michell CT, Apprill A, Voolstra CR. Endozoicomonas genomes reveal functional adaptation and plasticity in bacterial strains symbiotically associated with diverse marine hosts. Sci Rep. 2017;7:40579.View ArticlePubMedPubMed CentralGoogle Scholar
- Koren O, Rosenberg E. Bacteria associated with the bleached and cave coral Oculina patagonica. Microb Ecol. 2008;55:523–9.View ArticlePubMedGoogle Scholar
- Chen C-P, Tseng C-H, Chen CA, Tang S-L. The dynamics of microbial partnerships in the coral Isopora palifera. ISME J. 2011;5:728–40.View ArticlePubMedGoogle Scholar
- Li J, Chen Q, Long L-J, Dong J-D, Yang J, Zhang S. Bacterial dynamics within the mucus, tissue and skeleton of the coral Porites Lutea during different seasons. Sci Rep. 2014;4:7320.View ArticlePubMedPubMed CentralGoogle Scholar
- Morrow KM, Moss AG, Chadwick NE, Liles MR. Bacterial associates of two Caribbean coral species reveal species-specific distribution and geographic variability. Appl Environ Microbiol. 2012;78:6438–49.View ArticlePubMedPubMed CentralGoogle Scholar
- Miller WG, Parker CT, Rubenfield M, Mendz GL, Wösten MMSM, Ussery DW, Stolz JF, Binnewies TT, Hallin PF, Wang G, Malek JA, Rogosin A, Stanker LH, Mandrell RE. The complete genome sequence and analysis of the epsilonproteobacterium Arcobacter butzleri. PLoS One. 2006;2:e1358.View ArticleGoogle Scholar
- Frias-Lopez J, Zerkle AL, Bonheyo GT, Fouke BW. Partitioning of bacterial communities between seawater and healthy, black band diseased, and dead coral surfaces. Appl Environ Microbiol. 2002;68:2214–28.View ArticlePubMedPubMed CentralGoogle Scholar
- de Castro AP, Araújo SD, Reis AMM, Moura RL, Francini-Filho RB, Pappas G, Rodrigues TB, Thompson FL, Krüger RH. Bacterial community associated with healthy and diseased reef coral Mussismilia hispida from eastern Brazil. Microb Ecol. 2010;59:658–67.View ArticlePubMedGoogle Scholar
- Séré M, Wilkinson DA, Schleyer MH, Chabanet P, Quod J-P, Tortosa P. Characterisation of an atypical manifestation of black band disease on Porites lutea in the Western Indian Ocean. PeerJ. 2015;4:e2073.View ArticleGoogle Scholar
- Sweet M, Bythell J. Ciliate and bacterial communities associated with White syndrome and Brown Band disease in reef-building corals. Environ Microbiol. 2012;14:2184–99.View ArticlePubMedPubMed CentralGoogle Scholar
- Fan L, Liu M, Simister R, Webster NS, Thomas T. Marine microbial symbiosis heats up: the phylogenetic and functional response of a sponge holobiont to thermal stress. ISME J. 2013;7:991–1002.View ArticlePubMedPubMed CentralGoogle Scholar
- Lokmer A, Mathias Wegner K. Hemolymph microbiome of Pacific oysters in response to temperature, temperature stress and infection. ISME J. 2014;9:670–82.View ArticlePubMedPubMed CentralGoogle Scholar
- Boulotte NM, Dalton SJ, Carroll AG, Harrison PL, Putnam HM, Peplow LM, van Oppen MJ. Exploring the Symbiodinium rare biosphere provides evidence for symbiont switching in reef-building corals. ISME J. 2016;10:2693–701.View ArticlePubMedPubMed CentralGoogle Scholar
- Thomas L, Kendrick GA, Kennington WJ, Richards ZT, Stat M. Exploring Symbiodiniumdiversity and host specificity in Acropora corals from geographical extremes of Western Australia with 454 amplicon pyrosequencing. Mol Ecol. 2014;23:3113–26.View ArticlePubMedGoogle Scholar
- Ziegler M, Seneca FO, Yum LK, Palumbi SR, Voolstra CR. Bacterial community dynamics are linked to patterns of coral heat tolerance. Nat Comms. 2017;8:14213.View ArticleGoogle Scholar
- Rouzé H, Lecellier GJ, Saulnier D, Planes S, Gueguen Y, Wirshing HH, Berteaux-Lecellier V. Data S2: Symbiodinium 28S DNA sequences mentioned in the paper. PeerJ. 2017;5:e2856.View ArticlePubMedPubMed CentralGoogle Scholar
- Hume B, D'Angelo C, Burt J, Baker AC, Riegl B, Wiedenmann J. Corals from the Persian/Arabian Gulf as models for thermotolerant reef-builders: prevalence of clade C3 Symbiodinium, host fluorescence and ex situ temperature tolerance. Mar Pollut Bull. 2013;72:313–22.View ArticlePubMedGoogle Scholar
- Weber L, DeForce E, Apprill A. Optimization of DNA extraction for advancing coral microbiota investigations. Microbiome. 2017;5:18.View ArticlePubMedPubMed CentralGoogle Scholar
- Fabina NS, Putnam HM, Franklin EC, Stat M, Gates RD. Symbiotic specificity, association patterns, and function determine community responses to global changes: defining critical research areas for coral-Symbiodiniumsymbioses. Glob Change Biol. 2013;19:3306–16.Google Scholar
- Bourne DG, Bourne DG, Dennis PG, Dennis PG, Uthicke S, Uthicke S, Soo RM, Soo RM, Tyson GW, Tyson GW, Webster N, Webster N. Coral reef invertebrate microbiomes correlate with the presence of photosymbionts. ISME J. 2013;7:1452–8.View ArticlePubMedPubMed CentralGoogle Scholar
- Raina J-B, Dinsdale EA, Willis BL, Bourne DG. Do the organic sulfur compounds DMSP and DMS drive coral microbial associations? Trends Microbiol. 2010;18:101–8.View ArticlePubMedGoogle Scholar
- Fan TY, Lin KH, Kuo FW, Soong K, Liu LL. Diel patterns of larval release by five brooding scleractinian corals. Mar Ecol Prog Ser. 2006;321:133–42.View ArticleGoogle Scholar
- Harrison PL, Wallace CC. Reproduction, dispersal and recruitment of scleractinian corals: Ecosystems of the world. Amsterdam: Ed. Z Dubinsky Elsevier; 1990. pp. 133-207.Google Scholar
- Baird AH, Guest JR, Willis BL. Systematic and biogeographical patterns in the reproductive biology of scleractinian corals. Annu Rev Ecol Evol Syst. 2009;40:551–71.View ArticleGoogle Scholar
- Whitaker K. Genetic evidence for mixed modes of reproduction in the coral Pocillopora damicornis and its effect on population structure. Mar Ecol Prog Ser. 2006;306:115–24.View ArticleGoogle Scholar
- Yeoh S-R, Dai C-F. The production of sexual and asexual larvae within single broods of the scleractinian coral, Pocillopora damicornis. Mar Biol. 2009;157:351–9.View ArticleGoogle Scholar
- Gélin P, Fauvelot C, Mehn V, Bureau S, Rouzé H, Magalon H. Superclone expansion, long-distance clonal dispersal and local genetic structuring in the coral Pocillopora damicornis type β in Reunion Island, South Western Indian Ocean. PLoS One. 2017;12:e0169692.View ArticlePubMedPubMed CentralGoogle Scholar
- Neave MJ, Rachmawati R, Xun L, Michell CT, Bourne DG, Apprill A, Voolstra CR. Differential specificity between closely related corals and abundant Endozoicomonas endosymbionts across global scales. ISME J. 2016;11:186–200.View ArticlePubMedPubMed CentralGoogle Scholar
- Buddemeier RW, Fautin DG. Coral bleaching as an adaptive mechanism. Bioscience. 1993;43:320–6.View ArticleGoogle Scholar
- Putnam HM, Stat M, Pochon X, Gates RD. Endosymbiotic flexibility associates with environmental sensitivity in scleractinian corals. Proc R Soc B-Biol Sci. 2012;279:4352–61.View ArticleGoogle Scholar
- Winkler NS, Pandolfi JM, Sampayo EM. Symbiodinium identity alters the temperature-dependent settlement behaviour of Acropora millepora coral larvae before the onset of symbiosis. Proc Biol Sci. 2015;282:–20142260.Google Scholar
- Knowlton N, Rohwer F. Multispecies microbial mutualisms on coral reefs: the host as a habitat. Am Nat. 2003;162:S51–62.View ArticlePubMedGoogle Scholar
- Jones A, Berkelmans R. The photokinetics of thermo-tolerance in Symbiodinium. Mar Ecol. 2012;33:490–8.View ArticleGoogle Scholar
- Rowan R. Coral bleaching: thermal adaptation in reef coral symbionts. Nature. 2004;430:742.View ArticlePubMedGoogle Scholar
- Fabricius KE, Mieog JC, Colin PL, Idip D, van Oppen MJH. Identity and diversity of coral endosymbionts (zooxanthellae) from three Palauan reefs with contrasting bleaching, temperature and shading histories. Mol Ecol. 2004;13:2445–58.View ArticlePubMedGoogle Scholar
- Stat M, Gates RD. Clade D Symbiodiniumin scleractinian corals: A “Nugget” of hope, a selfish opportunist, an ominous sign, or all of the above? Journal of Marine Biology. 2011;2011:1–9.View ArticleGoogle Scholar
- Roth MS, Goericke R, Deheyn DD. Cold induces acute stress but heat is ultimately more deleterious for the reef-building coral Acropora yongei. Sci Rep. 2012;2:240.View ArticlePubMedPubMed CentralGoogle Scholar
- Tsang R, Ang P Jr. Cold temperature stress and predation effects on corals: their possible roles in structuring a nonreefal coral community. Coral Reefs. 2015;34:97–108.View ArticleGoogle Scholar
- Silverstein RN, Cunning R, Baker AC, Tenacious D. Symbiodiniumin clade D remain in reef corals at both high and low temperature extremes despite impairment. J Exp Biol. 2017;220:1192–6.View ArticlePubMedGoogle Scholar
- Kim B-M, Son S-W, Min S-K, Jeong J-H, Kim S-J, Zhang X, Shim T, Yoon J-H. Weakening of the stratospheric polar vortex by Arctic sea-ice loss. Nat Comms. 2013;5:4646.View ArticleGoogle Scholar
- Theis KR, Dheilly NM, Klassen JL, Brucker RM, Baines JF, Bosch TCG, Cryan JF, Gilbert SF, Goodnight CJ, Lloyd EA, Sapp J, Vandenkoornhuyse P, Zilber-Rosenberg I, Rosenberg E, Bordenstein SR. Getting the hologenome concept right: an eco-evolutionary framework for hosts and their microbiomes. mSystems. 2016;1:e00028–16.View ArticlePubMedPubMed CentralGoogle Scholar