Figures
Abstract
Global change and its associated temperature increase has directly or indirectly changed the distributions of hosts and pathogens, and has affected host immunity, pathogen virulence and growth rates. This has resulted in increased disease in natural plant and animal populations worldwide, including scleractinian corals. While the effects of temperature increase on immunity and pathogen virulence have been clearly identified, their interaction, synergy and relative weight during pathogenesis remain poorly documented. We investigated these phenomena in the interaction between the coral Pocillopora damicornis and the bacterium Vibrio coralliilyticus, for which the infection process is temperature-dependent. We developed an experimental model that enabled unraveling the effects of thermal stress, and virulence vs. non-virulence of the bacterium. The physiological impacts of various treatments were quantified at the transcriptome level using a combination of RNA sequencing and targeted approaches. The results showed that thermal stress triggered a general weakening of the coral, making it more prone to infection, non-virulent bacterium induced an ‘efficient’ immune response, whereas virulent bacterium caused immuno-suppression in its host.
Citation: Vidal-Dupiol J, Dheilly NM, Rondon R, Grunau C, Cosseau C, Smith KM, et al. (2014) Thermal Stress Triggers Broad Pocillopora damicornis Transcriptomic Remodeling, while Vibrio coralliilyticus Infection Induces a More Targeted Immuno-Suppression Response. PLoS ONE 9(9): e107672. https://doi.org/10.1371/journal.pone.0107672
Editor: Pikul Jiravanichpaisal, Fish Vet Group, Thailand
Received: May 16, 2014; Accepted: August 13, 2014; Published: September 26, 2014
This is an open-access article, free of all copyright, and may be freely reproduced, distributed, transmitted, modified, built upon, or otherwise used by anyone for any lawful purpose. The work is made available under the Creative Commons CC0 public domain dedication.
Data Availability: The authors confirm that all data underlying the findings are fully available without restriction. The reference transcriptome used is accessible at http://2ei.univ-perp.fr/telechargement/transcriptomes/blast2go_fasta_Pdamv2.zip. The raw data (untreated reads) for all treatments are publicly available at http://www.ncbi.nlm.nih.gov/sra with the accession number SRP029998.
Funding: This study was supported by the Agence Nationale de la Recherche through the Program BIOADAPT (ADACNI ANR-12-ADAP-0016-03) and the French-Israeli High Council for Science and Technology (P2R n °29702YG). Work in the Freitag laboratory was supported by start-up funds from the OSU Computational and Genome Biology Initiative and Oregon State University. The funders had no role in the study design, data collection and analysis, decision to publish, or preparation of the manuscript.
Competing interests: The co-author MF is a PLOS ONE editorial board member and this does not alter the authors’ adherence to PLOS ONE editorial policies and criteria. In addition, the authors have declared that no competing interests exist.
Introduction
Greenhouse gas emissions have increased since the industrial revolution, and the exceptionally high concentrations now reached have caused global climate change [1]. One of the major consequences of global climate change is that the overall seawater temperature has increased by a mean value of 0.6°C over the last 100 years [2]. The direct consequences of these changes are mass mortality of species having limited ability to adapt to temperature increase [2], [3], but indirect consequences can include mortalities triggered by temperature-dependent epizootics [4].
Marine disease risk is clearly enhanced under global warming conditions [5], [6], but the causes of the increase in the frequency and severity of epizootics are unclear and probably multi-factorial. It is known that global warming can change host and pathogen repartition to enhance the probability of host/pathogen encounters, or to facilitate new interactions. If some of these new host/pathogen interactions are compatible (susceptible host and/or virulent pathogen), a rapid spread of disease can occur [4]–[7], leading to massive terrestrial and marine vertebrate or invertebrate host die-offs [8]–[13]. In addition, an increase in temperature can diminish host immune abilities, leading to increased susceptibility to disease [14], [15] or increased pathogen virulence [16]–[18]. Taken together, these data show that the emergence of epizootics under climate change conditions is probably a result of multiple factors, the relative contribution and synergy of which must be studied to better understand and predict the epizootic risk in the current context of climate change [5], [19].
Among marine host/pathogen models, the interaction between the coral Pocillopora damicornis and the bacterium Vibrio coralliilyticus is useful for studying and unraveling the effects of temperature and bacterial virulence in pathogenesis. P. damicornis is widely distributed in the Indo-Pacific region [20], [21] and is highly susceptible to a wide range of disturbances [22]–[25], including disease [26]–[29]. V. coralliilyticus is a common pathogen [30] whose virulence is temperature-dependent [26]. The genomes of two V. coralliilyticus strains were recently sequenced, and a proteomic study performed under a range of temperature conditions has revealed the thermo-dependent expression of several putative virulence factors [31], [32], but their effects on coral physiology remain to be characterized. Recent data obtained using this infection model showed that the bacterium modulates the expression of several immune genes in P. damicornis [33], and that in its virulent state the Vibrio strongly decreases the expression of a gene encoding the antimicrobial peptide damicornin [34]. In this interaction the symptoms differ depending on the populations or strains involved. When V. coralliilyticus YB1 (isolated from Zanzibar) contacts its sympatric host, it triggers coral bleaching (i.e., substantial or partial loss of endosymbiotic dinoflagellate microalgae - commonly referred to as zooxanthellae - from coral tissues, and/or the loss or reduction of photosynthetic pigment concentrations within zooxanthellae) when the temperature reaches 29.5°C [22]. YB1 caused no pathogenicity in an allopatric coral from the Gulf of Eilat (Red Sea) at temperatures less than 24°C, triggered bleaching between 24 and 25°C, and caused tissue lysis at temperatures exceeding 25°C [35]. In another allopatric interaction involving a P. damicornis isolate from Lombok (Indonesia), no signs of infection were evident at temperatures up to 25°C, but at 28°C the bacteria penetrated the host tissues inducing rapid tissue lysis but no bleaching [33]. These results suggest that this bacterium has a high epizootic potential because of its ability to infect allopatric host populations and the various factors modulating its pathogenesis [32].
In this context, we used the infection model described above to investigate the relationship between temperature and bacteria (non-virulent and virulent), and their effects on coral physiology and pathogenesis. To this end, coral nubbins were exposed to stable or increasing temperature with or without bacterial addition, and the physiological changes caused at key stages of the interaction by the various treatments were assessed at the genome-wide scale using RNA sequencing (RNA-seq) approaches. We also used the quantitative reverse transcriptase polymerase chain reaction (q-RT-PCR) to study the response of 44 candidate immune genes during the stress response.
Methods
Biological material
The P. damicornis (Linnaeus, 1758) isolate used in this study was obtained from Lombok, Indonesia (Indonesian CITES Management Authority, CITES number 06832/VI/SATS/LN/2001-E; France Direction de l’Environnement, CITES number 06832/VI/SATS/LN/2001-I) and has been maintained in aquaria since the year 2001. Analysis of the first 600 bp of the mitochondrial ORF marker [36] showed that morphologically and genetically this P. damicornis isolate corresponds to the recently re-characterized P. damicornis clade [20], [37], [38]. For use in the experiments, coral explants (7 cm high, 6 cm diameter) were detached from the parent colony, and left to recover for a period of 1 month prior to use.
The coral pathogen Vibrio coralliilyticus strain YB1 [39]; CIP 107925, Institut Pasteur, Paris, France) was used to challenge or infect P. damicornis [22]. V. coralliilyticus was cultured as previously described [33].
Stress protocol
As the objective of the study was to detect stress effects and to avoid the potential influence of inter-individual variability and genetic background on these effects, the experiments were performed using clones of the same P. damicornis isolate and a unique bacterial strain (YB1). Coral nubbins were randomly placed in 120 L tanks (n = 27 per tank) and acclimatized at 25°C for a period of two weeks prior to initiating the treatments. To distinguish the effects of bacterial stress (exposure to V. coralliilyticus under non-virulence conditions), thermal stress and bacterial infection (V. coralliilyticus virulence activated by temperature increase) on the coral host, three independent treatments and a control were established as previously described [33]. For the non-virulent treatment, V. coralliilyticus was regularly added to coral nubbins held at 25°C, which is a sub-virulent temperature in this host/pathogen interaction. For the virulent treatment, V. coralliilyticus was regularly added to coral nubbins while the temperature was gradually increased from 25°C to 32.5°C, triggering activation of bacterial virulence and infection. For the control, the coral nubbins were held at 25°C without added bacteria. For the thermal stress treatment, the coral nubbins were subjected to a gradual increase in temperature without the addition of bacteria, triggering thermal stress that induces bleaching. For all treatments and the control, three nubbins were randomly sampled every three days and immediately frozen and stored in liquid nitrogen. Coral health and the stability/breakdown of the coral symbiosis were monitored over time in each experimental group by assessing the density of zooxanthellae and visual monitoring, as previously described [33].
cDNA library construction and high-throughput sequencing
Four cDNA libraries, corresponding to the control and treatment groups, were sequenced using Illumina GAIIx technology. One lane per library was used to generate an RNA-seq dataset of four 80-nucleotide paired-end short sequence reads. The samples for this analysis were collected at day 12, which was the last sampling occasion prior to the appearance of symptoms of bleaching or bacterial infection in the thermal stress and virulent treatments, respectively.
The total RNA from the 3 coral nubbins of each of the control and treatment groups was extracted using TRIzol reagent (Invitrogen), as previously described [40]. The cDNA libraries were then constructed following an established protocol [41]. Briefly, to obtain high quality mRNA, two cycles of oligo (dT) purification were performed using the Ambion Poly(A) Purification kit (Ambion). For the second round of purification, the mRNA purified in the first round was used in excess (4 µg mRNA per library). Subsequently, first strand cDNA synthesis was performed by combining 1 µg of the purified mRNA, 300 µg of random hexamer primers and the superscript III reverse transcription kit (Invitrogen) reagent. The excess of random primers enabled generation of first strand cDNA with an average length of 400 bp, and maximized coverage at the 5′ and 3′ ends of the mRNA. The second strand synthesis was performed using a combination of RNase H enzyme and DNA polymerase I (New England Biolabs). The resulting cDNA was purified using the QIAquick PCR purification kit (Qiagen). The amount and quality of nucleic acid obtained using this protocol was determined using a bioanalyzer and a nanodrop apparatus. The purified cDNA (1 µg) was used to generate each of the paired-end Illumina sequencing libraries. The libraries were prepared using Illumina adapter and PCR primers, according to previously published protocols [42], [43]. Libraries with an average insert size of 400–500 bp were isolated, and the concentration was adjusted to 10 nM. The samples (7 pM per sample) were loaded into separate channels of an Illumina GAIIx sequencer. The sequencing was performed at the Oregon State University Center for Genome Research and Biocomputing. The raw data (untreated reads) for all treatments are publicly available (http://www.ncbi.nlm.nih.gov/sra study accession number SRP029998).
Differential gene expression analysis
To assess changes in gene expression induced by the various treatments, reads from each library were mapped against a previously assembled reference transcriptome [44]. This transcriptome was assembled de novo from 80-nucleotide paired-end short sequence reads based on 6 lanes of Illumina sequencing, and contained 72,890 contigs. Among these contigs, 27.7% and 69.8% were predicted to belong to the symbiont and the host transcriptomes, respectively. Each sequencing lane contained cDNA prepared from coral nubbins of the same genotype. The sequencing lanes were loaded with cDNA from corals exposed to: 1) thermal stress (lane 1); 2) V. coralliilyticus in a non-virulent state (lane 2); 3) V. coralliilyticus in a virulent state (lane 3); 4) a constant temperature of 25°C (lane 4); 5) a pH of 7.4 for 3 weeks (lane 5); and 6) a pH of 8.1 for 3 weeks. This reference transcriptome is accessible at http://2ei.univ-perp.fr/telechargement/transcriptomes/blast2go_fasta_Pdamv2.zip.
The mapping was conducted using Burrows-Wheeler Aligner (BWA), using the default options: aln -n = 0.04; aln -o = 1; aln -e = -1; aln -d = 16; aln -i = 5; aln -l = -1; aln -k = 2; aln -M = 3; aln -O = 11; aln -E = 4 [45]. As RNA-seq data are a function of both the molar concentration and the transcript length, the results of the mapping step were corrected and expressed as reads per kilobase per million mapped reads (RPKM; [46]). To identify significantly different gene expressions among the control and treatment groups, the MARS method (MA plot-based methods using a random sampling model) of the R package in DEGseq was used [47].
q-RT-PCR
Quantitative real-time PCR (q-RT-PCR) was used to validate the expression profiles obtained from the MARS DEGseq analysis of the RNA-seq data. It was also used to measure expression levels of selected candidate genes during the non-virulent, thermal stress and virulent treatments. Total RNA was extracted and treated with DNase, and the poly(A) RNA was purified as described above. Approximately 50 ng of purified poly(A) RNA was reverse transcribed with hexamer random primers using ReverTAid H Minus Reverse Transcriptase (Fermentas). The q-RT-PCR experiments were performed using cDNA obtained from three coral nubbins per control and treatment group, as described previously [33]. For each candidate gene the level of transcription was normalized using the mean geometric transcription rate of three reference sequences encoding ribosomal protein genes from P. damicornis (60S ribosomal protein L22, GenBank accession number HO112261; 60S ribosomal protein L40A, accession number HO112283; and 60S acidic ribosomal phosphoprotein P0, accession number HO112666). The stable expression status of these three genes during biotic and abiotic stress has been demonstrated previously [33].
Statistical analysis
Differences in gene expression data obtained using the RNA-seq analyses were analyzed for statistical significance using the MARS method, as described above [47]. Differences in transcript levels among experimental groups were considered statistically significant at p<0.0001. The statistical analyses used to identify biological processes (Gene Ontology term) that were significantly up- or down-regulated [48] were performed using Blast2GO software (version 2.6.4). Increases in processes were considered statistically significant at p<0.05. The hierarchical clustering of q-RT-PCR data was performed using Multiple Array Viewer software (version 4.8.1), with average linkage clustering using the Pearson correlation as a default distance metric. The normality of the RPKM distribution was assessed using the Kolmogorov-Smirnov test. As the data were not normally distributed, we used non-parametric statistical procedures. The Mann-Whitney U-test was used to compare the expression level between the core set of up-regulated genes and the core set of down-regulated genes. These statistical analyses were conducted using SPSS 10.0 (Kolmogorov-Smirnov and Mann-Whitney), and were considered statistically significant at p<0.05.
Results
Gene expression analysis and validation
The global transcriptomic approach was conducted using RNA-seq methodology that was applied to samples collected 12 days following initiation of the various treatments. To distinguish the bacterial virulence effect from the thermal stress effect in the virulent treatment, in which corals were submitted to a combination of added bacteria and increased temperature, we used the same methods to compare the transcript content in the thermal stress and the virulent treatments. This comparison was called the virulence effect. For each treatments and comparison the results of the: i) Illumina sequencing, ii) quality control and filtering, iii) mapping to the reference transcriptome, and the iv) statistical approach (DEGseq) are summarized in the Table 1. To validate results from the RNA-seq analysis we used an alternative method (q-RT-PCR) for transcript quantification. For this analysis, 22 host transcripts were selected along the gradient of expression from highly up-regulated to highly down-regulated genes, based on RNA-seq data. A significant correlation between the log2-fold change in expression of the RNA-seq and the qRT-PCR data was obtained (log2-fold change qRT-PCR vs. log2-fold change RNA-seq: control vs. non-virulent, r2 = 0.87 and p<0.0001; control vs. thermal stress, r2 = 0.89 and p<0.0001; control vs. virulent, r2 = 0.875 and p<0.0001; thermal stress vs. virulent, r2 = 0.64 and p<0.0001; Fig. 1) and validated the results obtained through the RNA-seq analytical pipeline.
Twenty-two genes were arbitrary selected, from highly up-regulated to highly down-regulated contigs. Their levels of expression were quantified by q-RT-PCR, and the results were compared with those obtained using the RNA-seq approach. The log2 changes in expression based on q-RT-PCR and RNA-seq analyses were closely correlated for all treatments, indicating the accuracy of the RNA-seq approach for quantification.
Specific biological processes modulated in the various experimental conditions
A Gene Ontology (GO) term enrichment analysis [48] was performed on the sets of transcripts showing a significant change in expression. The results showed that in response to the non-virulent treatment, 12 biological processes were significantly enriched in the set of up-regulated genes (p<0.05; Table 2), and 11 in the set of down-regulated genes (p<0.05; Table 3). In response to the thermal stress treatment 20 biological processes were significantly more represented in the set of up-regulated transcripts relative to the control (p<0.05; Table 2) and 20 from the set of down-regulated genes (p<0.05; Table 3). In response to the virulent treatment six biological processes were significantly enriched (p<0.05) from the set of up-regulated genes (p<0.05; Table 2) and 20 from the set of down-regulated genes (p<0.05; Table 3).
We then investigated whether there were conserved up- or down-regulated genes in response to all treatments. This analysis showed that 229 genes were significantly up-regulated in response to the non-virulent, thermal stress and virulent treatments, while 1372 were significantly down-regulated in all experimental groups. All the genes showing conserved regulation among groups are referred to as core response genes. Four biological processes were significantly enriched in the set of core up-regulated genes (p<0.05; Table 4), and 15 from the core down-regulated genes (p<0.05; Table 4). To distinguish the bacterial virulence effect in the virulent treatment, in which corals were submitted to a combination of added bacteria and increased temperature, the same enrichment analysis was performed comparing transcriptomic data obtained from the thermal stress and the virulent treatments. This enabled identification of six enriched biological processes from the set of genes up-regulated in the virulence treatment (p<0.05; Table 4) and nine in the set of genes down-regulated (p<0.05; Table 4).
The annotation and expression levels of the genes belonging to each enriched GO category are shown in Table S1 & S2.
Selection of immune candidate genes and expression analysis
As the enrichment analysis revealed a clear effect of the various treatments on the immune function, we investigated and compared the expression of putative P. damicornis immune genes, amongst the various treatments. These genes, annotated manually or using Blast2GO, corresponded to the immune toolbox of P. damicornis. They included genes encoding proteins involved in recognition (e.g. TLR and lectins), signaling pathways (e.g. NF-kB, AP1/ATF, c-Jun), complement system (e.g. C3, C-type lectin, and membrane attack complex), prophenoloxidase cascade (e.g. prophenoloxidase activating enzyme, laccase), leukotriene cascade (e.g. 5-lypoxigenase, leukotriene A4 hydrolase), antimicrobial molecules (e.g. LPBPI and antimicrobial peptides) and ROS scavengers (e.g. peroxidases, catalases, GFP-like molecules). The results supporting the annotation of these genes are presented in the Table S3.
To assess the impact of each treatment on the regulation of these genes, their expression was assessed every 3 days during the experimental period, using q-RT-PCR. To confirm that the observed regulation of gene expression was not because of physiological collapse of the coral or an experimental artifact, four housekeeping genes were also monitored during the experimental period. These genes were selected from amongst host genes that were not regulated in any treatment. As expected, the results showed that their expression remained stable in all treatments (Fig. 2; Table S4).
The data included q-RT-PCR results for nubbins sampled during the non-virulent, thermal stress, virulent and bacterial virulence (thermal stress vs. virulent treatment) treatments (days 3, 6, 9, 12, 15 and 18). Quantification was normalized to the control conditions for the non-virulent, thermal stress and virulent treatment, and with results for colonies sampled at the same temperature as that for the thermal stress vs. virulent comparison (bacterial virulence effect only). The results are presented as a log2-fold change in expression. The hierarchical clustering of the q-RT-PCR data was done using Multiple Array Viewer software (version 4.8.1), with average linkage clustering based on the Pearson correlation as a default distance metric. Cluster C1 represents the response to thermal stress, cluster C2 represents the response to bacteria (non-virulent and virulent), and cluster C3 represents the response to the virulence of the bacteria. The numbers at the bottom of the figure correspond to the following genes: 1, epsilon isoform 1 (housekeeping control); 2, MASP3; 3, peroxidase2; 4, peroxidase1; 5, cyclin d2 (housekeeping control); 6, preprotein translocase SecY subunit (housekeeping control); 7, laccase; 8, prophenol oxidase activating enzyme; 9, ubiquitin-conjugating enzyme E2 (housekeeping control); 10, 5-lypoxigenase; 11, catalase1; 12, nucleoredoxin; 13, TAK1; 14, MKK7; 15, MKK4; 16, JNK; 17, TRAF6; 18, IKBa; 19, NF-kB; 20, TIR2; 21, AP1; 22, ATF; 23, Bf; 24, C3; 25, MASP1; 26, catalase2; 27, MKK3/6; 28, IKK; 29, p38; 30, LRR2; 31, apextrin; 32, SOD1; 33, leukotriene A4 hydrolase; 34,Leukotriene C4 synthase; 35, MyD88; 36, GFP-Like2; 37, TIR3; 38, SOD2; 39, MEKK1; 40, LPBPI; 41, LRR-TIR-IGG; 42, GFP-Like1; 43, Tx60A2; 44, Tx60A1; 45, PdC-Lectin; 46,damicornin; 47, phospholipase A2; 48, mytimacin-like.
A hierarchical clustering approach was used to create sample and gene trees for the 24 samples (six for each treatment) and the 48 genes (44 immune genes and four housekeeping genes; Fig. 2). Three distinct clusters were identified. The first cluster mainly represented the response to thermal stress (cluster C1; Fig. 2), and contained two subgroups of genes that were highly modulated. The first subgroup comprised genes displaying a high level of down-regulation, including those encoding the antimicrobial peptide damicornin (average 40.0-fold decrease; Fig. 2), the mannose binding lectin PdC-Lectin (average 6.9-fold decrease; Figs 2 and 3), and two members of the membrane attack complex (Tx60A1 and A2; average 91.2-fold decrease; Fig. 2). The second subgroup comprised genes displaying a high level of up-regulation, including those encoding a putative MASP1, putative C3 and Bf components of the complement pathway (average 26.4-fold increase; Figs 2 and 3), and two immune-related transcription factors ATF and AP1 (average 19.3-fold increase; Figs 2 and 3). The second cluster (cluster C2; Figs 2 and 3) corresponded to genes involved in the response to bacteria (non-virulent and virulent treatments). It included genes subject to down-regulation (e.g. a gene encoding a putative apextrin protein; average decrease 3.2-fold; Figs 2 and 3), and some showing up-regulation (e.g. genes encoding prophenoloxidase activating enzyme and laccase, which are two key elements of the prophenoloxidase pathway; average 8.1-fold increase; Figs 2 and 3). The third cluster (cluster C3; Fig. 2) corresponded to genes involved in the response to the virulence of the bacteria. This cluster was mainly characterized by general down-regulation of the innate immune recognition and signaling pathways; the genes encoding the LRR and TIR domains containing proteins, TAK1, IKK, NF-kB, MKK3/6, P38, IkBa, AP1/ATF and MKK4/7 were all down-regulated by an average factor of 4.1 (Figs 2 and 3).
Reconstitution of the immune pathways identified in previous studies and from the present study (see Table 5 for references). Arrows highlight the average response (if any) of each pathway to each treatment or comparison. Green arrow: response to the non-virulent treatment; yellow arrow: response to the thermal stress treatment; red arrow: response to the virulent treatment; violet arrow: response to the virulence effect (the comparison between the thermal stress and the virulent treatment).
In relation to specific biological functions, the q-RT-PCR experiments provided noteworthy results for several immune genes and the main ones are highlighted below. For antimicrobial effectors, we found that expression of the damicornin, mytimacin-like and LBP–BPI genes was decreased in the thermal stress and virulent treatments (Figs 2 and 3). The two genes of the prophenoloxidase pathway (prophenoloxidase activating enzyme and laccase) were co-up-regulated in the presence of bacteria (non-virulent and virulent treatments) whereas in response to bacterial virulence, there were strongly down-regulated at day 18. The complement pathway (PdC-lectin, MASP1 and 2, C3 and Bf) was mainly up-regulated by the presence of non-virulent bacteria and during the virulent treatment. Members of the membrane attack complex Tx60 A1 and A2, and apextrin (a downstream component of the complement pathway; Fig. 3) were strongly down-regulated under temperature stress conditions (thermal stress and virulent treatment). The expression of some key components of the backbone of the recognition and signaling pathways (LRR, TIR, Myd88 and TRAF6) was disturbed during the virulent, virulence and thermal stress treatments. The downstream pathways of this backbone (NF-kB, ATF/AP1 and JNK) showed a similar trend of regulation in response to the non-virulent bacteria, thermal stress and virulent bacteria. In response to bacterial virulence the NF-kB and ATF/AP1 pathways were completely down-regulated (Figs 2 and 3). The leukotriene cascade responded mainly to temperature stress by general up-regulation of its component. Among the genes encoding ROS scavengers, two groups could be clearly distinguished. The first group comprised ‘classical’ ROS scavengers (peroxidase, catalase, superoxide dismutase and nucleoredoxin), which showed strong up-regulation in response to thermal stress, non-virulent bacteria and virulent bacteria. The second group (encoding GFP-like protein); was initially up-regulated by thermal stress but then down-regulated at higher temperature (Figs 2 and 3). In response to the non-virulent bacterium these two genes were up-regulated, but were down-regulated immediately following initiation of virulence in the bacterium (virulent treatment, days 6–18; Figs 2 and 3).
Discussion
By combining the analysis of natural coral/Vibrio interaction, experimental exposures of corals to the bacteria, global transcriptomic studies (RNA-seq) confirmed by q-RT-PCR, we showed that thermal stress induces a general decrease in coral gene expression; this included decreased expression of immune genes, hence reducing the immune capacities of the coral. We also found that virulent bacteria triggered a marked suppression of host immunity. The RNA-seq approach unveiled a large range of transcripts expressed either by the coral host or by its dinoflagellate symbionts (entire holobiont response). A clear illustration of the response of the symbiont was the down-regulation of photosynthesis genes and for the host, the up-regulation of immune genes in response to the non-virulent treatment). In general it was possible to assign transcripts to either the host or the symbiont partner (as explained in the methods section), although for some genes the discrimination was not always obvious. Assembled genomes for the two partners in the symbiosis are lacking, and dinoflagellate genomes are known to contain numerous copies of orthologous genes, which can be divergent. The latter factor increases the complexity of the transcriptome [49]–[52] and can make assignment difficult. In addition, our transcriptome probably contained xeno-contaminant sequences from bacteria, eukaryotic species and RNA viruses [53], which are part of the coral holobiont [54].
The sophisticated ancestral immune toolkit responds to experimental bacterial exposure
The immunity of scleractinian corals is an expanding field in basic research because of the basal position of cnidarians in the eumetazoan tree of life, and its relevance to understanding the evolution of defense systems against pathogens. It is also crucial to understand the immune capabilities of corals, because of the increase of coral epizootics. In this context, most past studies have focused on the identification of immune pathways and effectors, mainly through searches for orthologous genes in the expanding genome and transcriptome databases. These studies have revealed that the common ancestor of animal species contained several immune genes and pathways that are also present in higher invertebrate and vertebrate species (Table 5). This discovery led to the hypothesis that early eumetazoan species had an ancestral immune toolkit [61], [81], which evolved to the more sophisticated systems found in higher lineage. This hypothesis remains to be proven because the functions of the putative immune genes and their involvement in immunity remain to be investigated [82]. Demonstration for this hypothesis was not trivial in our non-model organisms. Indeed, several immune genes are involved in a number of cellular processes (e.g. NF-kB transcription factor), but demonstrating their involvement in immunity will require knock-out or knock-down approaches that will have to be developed for non-model organisms, including corals. In this context, a targeted and organized response to a pathogen could be considered initial evidence that these genes and pathways are directly involved in the immune response. Some of the genes and proteins have been shown to respond to experimental infection or to pathogen elicitors (Table 5), but studies at the entire transcriptome level have not occurred. The present study partly addressed this shortcoming by showing that most of the known immune factors in the coral were modulated following a natural bacterial challenge. Indeed, we observed co-regulation of genes belonging to the same pathways, including genes for: (i) MASP1, C3 and Bf in the complement pathway; ii) the prophenoloxidase activating enzyme and laccase in the melanization pathway; and (iii) proteins involved in signaling pathways, including the transcription factors NF-kB and JNK (Fig. 3). In summary, the results of this study, and recent studies of A. millepora [83] and Gorgonia ventalina [84] indicate that corals are able to mount an immune response using an ancestral immune toolkit.
The effect of non-virulent bacteria
In the presence of the non-virulent bacteria, we found that innate immune pathways (Table 2 and Fig. 2) were activated in the coral, which is consistent with previous data showing that non-virulent V. coralliilyticus triggers an immune response [33], [34]. The absence of bacteria in host tissues during the non-virulent treatment [33] suggests that the coral either: (i) detected the presence of bacteria, and the coral cells developed an immune response directed against non-internalized bacteria; or (ii) was able to develop an efficient immune response that killed all bacteria on entry.
The mechanisms underlying the recognition of V. coralliilyticus by P. damicornis remain to be identified. However, in Hydra spp. (also a cnidarian) the unconventional pathogen sensors HyLRR, which are expressed at the surface of epithelial cells, interact with an intracellular HyTRR (a TIR domain-containing protein) during immune challenge [55], and this leads to the expression of an AMP (periculin-1). These LRR and TRR molecules were identified in the P. damicornis transcriptome, suggesting that P. damicornis can detect V. coralliilyticus using similar mechanisms. The ability to mount an immune response has been reported in numerous cnidarian species exposed to lipopolysaccharide (LPS) including Montastraea faveolata, Stephanocoenia intersepta, Porites astreoides and Acropora millepora [70], [77]. However, our coral/pathogen model is the first to have enabled an association to be made between the absence of infection and an organized up-regulation of multiple and interlinked pathways that reflect the three key steps in the immune response: i) pathogen recognition; ii) signal transduction; and iii) effector responses (Figs 2 and 3).
The thermal stress effect
Consistent with several previous studies (for review, see; [85], we found that the response to thermal stress is characterized by up-regulation of several heat shock proteins (HSPs) and genes encoding ROS scavengers. However, we showed that this was accompanied by significant down-regulation of genes involved in innate immunity and apoptosis (Table 3). These observations support the hypothesis that high temperature favors infection because it impacts coral immune function [86]–[90]. However, this phenomenon may not occur with all immune genes. Indeed, our candidate gene approach showed that several immune genes were up-regulated in the thermal stress treatment (Fig. 2). Similar results have been reported in previous studies investigating the immune abilities of corals in response to thermal stress. Prophenoloxidase activity was shown to be higher in thermally stressed M. faveolata corals relative to healthy or diseased corals, whereas the opposite occurred for antibacterial and lysozyme-like activities [71]. Such results raise questions about the uniform nature of thermal stress effects on the coral immune response, and challenge the hypothesis that high temperature negatively impacts all coral immune functions and always favors infection [86]–[90]. Nevertheless, we showed that several encoding pattern recognition receptors (PRR; Toll-like, LBP–BP and mannose-binding lectin) followed a general trend of down-regulation, as did key factors in the downstream signaling pathway, including Myd88 (Figs 2 and 3). Other essential components of the immune response, and more particularly some effectors, were also down-regulated (damicornin, TX60A-1 and TX60A-2) during thermal stress, which may significantly alter the immune capabilities of the coral. Thus, the combination of our results and those reported previously suggest that, even if the immune abilities of P. damicornis are not completely inhibited by thermal stress, the changes induced by increased temperature reduce the capacity of the coral to mount an efficient immune response, thus facilitating pathogen colonization of host tissues.
The virulence +/− thermal stress effect
Analysis of the transcriptomic changes between the virulent treatment and the control (thermal stress+virulence effect) or the thermal stress (virulence effect) treatment indicated profound remodeling of the transcriptome. Coral immune function was highly affected by the infection process, and many genes were down-regulated, particularly following internalization of the bacteria. This finding is consistent with evidence from a previous study of P. damicornis, which showed down-regulation of the gene encoding the AMP damicornin during infection by V. coralliilyticus [34]. Nevertheless, this down-regulation response is not absolute, as revealed using the candidate gene approach, which demonstrated the up-regulation of some genes. This approach also highlights that the coral response to the virulent and non-virulent bacteria was very similar (Fig. 2; cluster C2). This suggests that the response is mediated following interaction with the bacteria, and is sufficient to neutralize the bacteria under non-virulent conditions, but not under virulent conditions. This shows that under virulent conditions the bacteria were able to circumvent the coral response and to establish the pathogenic intracellular and intra-vesicular form, which involves localization that facilitates protection of the bacteria from systemic immune attack [26], [33]. This strategy is also used by the oyster (Crassostrea gigas) pathogen V. splendidus: following internalization in immune cells this bacterium manipulates and evades host defenses by preventing acidic vacuole formation and limiting ROS production [91]. At high temperature and at the same stage of infection, V. coralliilyticus has also been shown to be able to express a series of putative virulence factors involved in host degradation, secretion, antimicrobial resistance and transcriptional regulation, which may help this pathogen to survive and spread in the intracellular environment [31], [32].
The greater the stress, the greater the transcriptomic disturbance: a trade-off mechanism?
Based on the symptoms evident in each treatment, exposure to the bacteria under high temperature conditions triggering virulence resulted in greater stress to the coral than exposure to thermal stress alone, which in turn was more stressful than exposure to the non-virulent bacteria. In the virulent treatment tissue lysis and death of the coral colonies occurred within 18 days, in the high temperature treatment the symbiosis broke down and coral bleaching occurred within 18 days, and no symptoms were observed in the treatment involving non-virulent bacteria [33]. Our RNA-seq results showed that 12 days following initiation of the treatments there was a clear correlation between the intensity of the stress and the number of down-regulated genes: the greater the stress, the greater the down-regulation (Table 1). Extensive down-regulation in response to environmental stressors is often reported in transcriptomic studies, but most research has focused on biological functions that are expected to be affected, rather than using the global analysis of the physiological status of the impacted organism. Indeed, it is easier and more straightforward to focus on down-regulated genes clearly associated with specific functions than to interpret the modulation of numerous genes associated with diverse, general biological processes. Nonetheless, the phenomenon of widespread down-regulation during stress responses has been verified through genome-wide transcriptomic studies in corals and other organisms exposed to various environmental stressors [44], [92]–[97]. This suggests that down-regulation is a conserved phenomenon under stress conditions.
In this context, the coral stress response was recently compared with the environmental stress response (ESR) of a budding yeast [98], which involves a conserved, extreme, rapid, and genome-wide response to a broad set of environmental stressors that trigger a common response among a large set (approximately 900) of genes [99]. In corals the core stress response involves the up-regulation of HSP and ROS scavenger encoding genes, and disruption of the expression of genes involved in Ca2+ homeostasis, cytoskeleton organization, cell signaling and transcriptional regulation [98]. However, as was found in the yeast, in which approximately 66% of the genes of the ESR were down-regulated [99], the broad down-regulation observed in corals exposed to various stressors may also be a fundamental part of the coral stress response. This hypothesis is supported by parallels between the coral response in our study, and the main biological processes down-regulated in the yeast ESR. As in the yeast ESR [99], our GO enrichment analysis (Table 4) suggested a general down-regulation of what are usually considered to be housekeeping biological processes, including translation, nucleic acid metabolic processes, secretion, regulation of gene expression and others biosynthetic processes. These housekeeping functions are generally very costly, and monopolize most energy production and the cell’s transcriptional and translational machinery [100]. Under stress conditions, down-regulation may help to minimize energetic expenditure, enabling rapid and efficient responses to the new environmental conditions [101]. In yeast, this energetic economy was shown to affect functions that are not essential for survival, but also involved the selective down-regulation of genes encoding high molecular weight and/or highly expressed proteins [101]. Demonstrating such phenomena in corals is difficult because of the lack of tools enabling complete transcriptome analysis. However, the results obtained in our study for the host core genes were consistent with this hypothesis. Indeed, we found that under the control conditions the core coral genes were more represented among the down-regulated genes (RPKM average, 789.40 reads) than the up-regulated genes (RPKM average, 546.74 reads; Mann-Whitney U test, p<0.01).
Conclusion
Taken together, our results show that scleractinian corals have an immune system that is able to respond to pathogenic agents, and support the “sophisticated ancestral immune toolkit” hypothesis [61], [81], [102]. The results also have important implications for understanding the immune system of ancestral metazoans, and its evolution and conservation through the eumetazoan lineage. Our results demonstrate that thermal stress alters immune capacities, especially during recognition and antibacterial processes, and that this probably facilitates pathogenesis. However, bacterial virulence and the intracellular localization of the pathogen seem to be the major factors responsible for pathogenesis, especially through an immuno-suppressive effect that may decrease the efficiency of the immune response, and lead to bacterial proliferation and physiological collapse of the coral. In addition to effects on immunity, thermal stress induces a strong and energetically costly response that may weaken the coral and favor infection by specific or opportunistic pathogens, through a trade-off mechanism. In summary, the absolute effect of thermal stress on the coral is less than that of the virulent bacteria during pathogenesis, but is a clear facilitator of infection.
Supporting Information
Table S1.
Detailed results of the GO term enrichment analysis for the up-regulated set of genes.
https://doi.org/10.1371/journal.pone.0107672.s001
(XLS)
Table S2.
Detailed results of the GO term enrichment analysis for the up-regulated set of genes.
https://doi.org/10.1371/journal.pone.0107672.s002
(XLSX)
Table S3.
Candidate innate immune genes, Annotation, Contig name, Sequence, Primer, BlastX results and protein domain.
https://doi.org/10.1371/journal.pone.0107672.s003
(XLSX)
Table S4.
Expression of innate immune candidate genes, numerical values.
https://doi.org/10.1371/journal.pone.0107672.s004
(XLSX)
Acknowledgments
The facilities (bioinformatics and molecular biology) of the Tecnoviv platform were used for the study. The authors are indebted to Jérôme Bossier for his help in statistical procedures, and to Richard Galinier and Jean-François Allienne for technical assistance on the Tecnoviv platform.
Author Contributions
Conceived and designed the experiments: JVD CG CC MF GM. Performed the experiments: JVD RR CG CC KMS MF. Analyzed the data: JVD NMD RR CG KMS MA GM. Contributed reagents/materials/analysis tools: NMD RR MF MA. Contributed to the writing of the manuscript: JVD NMD RR CC CG KMS MA MF GM.
References
- 1.
Pachauri RK (2007) Climate change 2007: Synthesis report. Geneva: IPCC Secreteria.
- 2. Hoegh-Guldberg O, Bruno JF (2010) The impact of climate change on the world’s marine ecosystems. Science 328(5985): 1523–1528.
- 3. Walther GR, Post E, Convey P, Menzel A, Parmesan C, et al. (2002) Ecological responses to recent climate change. Nature 416(6879): 389–395.
- 4. Harvell CD, Mitchell CE, Ward JR, Altizer S, Dobson AP, et al. (2002) Climate warming and disease risks for terrestrial and marine biota. Science 296(5576): 2158–2162.
- 5. Harvell D, Altizer S, Cattadori IM, Harrington L, Weil E (2009) Climate change and wildlife diseases: When does the host matter the most? Ecology 90(4): 912–920.
- 6. Ward JR, Lafferty KD (2004) The elusive baseline of marine disease: Are diseases in ocean ecosystems increasing? PLoS Biol 2(4): e120.
- 7. Harvell CD, Kim K, Burkholder JM, Colwell RR, Epstein PR, et al. (1999) Emerging marine diseases–climate links and anthropogenic factors. Science 285(5433): 1505–1510.
- 8. Burreson EM, Ragone Calvo LM (1996) Epizootiology of Perkinsus marinus disease of oysters in Chesapeake Bay, with emphasis on data since 1985. J Shellfish Res 15: 17–34.
- 9. Cook T, Folli M, Klinck J, Ford S, Miller J (2004) The relationship between increasing sea-surface temperature and the northward spread of Perkinsus marinus (Dermo) disease epizootics in oysters. Estuar Coast Shelf Sci 46(4): 587–597.
- 10. Devaux CA (2012) Emerging and re-emerging viruses: A global challenge illustrated by Chikungunya virus outbreaks. World J Virol 1(1): 11–22.
- 11. Pounds AJ, Bustamante MR, Coloma LA, Consuegra JA, Fogden MPL, et al. (2006) Widespread amphibian extinctions from epidemic disease driven by global warming. Nature 439(7073): 161–167.
- 12. Scheibling RE, Hennigar AW (1997) Recurrent outbreaks of disease in sea urchins Strongylocentrotus droebachiensis in Nova Scotia: evidence for a link with large-scale meteorologic and oceanographic events. Mar Ecol-Prog Ser 152(1): 155–165.
- 13. Skerratt L, Berger L, Speare R, Cashins S, McDonald K, et al. (2007) Spread of chytridiomycosis has caused the rapid global decline and extinction of frogs. EcoHealth 4(2): 125–134.
- 14. Campbell AH, Harder T, Nielsen S, Kjelleberg S, Steinberg PD (2011) Climate change and disease: bleaching of a chemically defended seaweed. Glob chang biol 17(9): 2958–2970.
- 15. Vargas-Albores F, Hinojosa-Baltazar P, Portillo-Clark G, Magallon-Barajas F (1998) Influence of temperature and salinity on the yellowleg shrimp, Penaeus californieinsis Holmes, prophenoloxidase system. Aquac Res 29(8): 549–553.
- 16. Case RJ, Longford SR, Campbell AH, Low A, Tujula N, et al. (2011) Temperature induced bacterial virulence and bleaching disease in a chemically defended marine macroalga. Environ Microbiol 13(2): 529–537.
- 17. Decostere A, Haesebrouck F, Turnbull JF, Charlier G (1999) Influence of water quality and temperature on adhesion of high and low virulence Flavobacterium columnare strains to isolated gill arches. J Fish Dis 22(1): 1–11.
- 18. Kushmaro A, Rosenberg E, Fine M, Loya Y (1997) Bleaching of the coral Oculina patagonica by Vibrio AK-1. Mar Ecol-Prog Ser 147: 159–165.
- 19. Ellner SP, Jones LE, Mydlarz LD, Harvell CD (2007) Within-host disease ecology in the sea fan Gorgonia ventalina: Modeling the spatial immunodynamics of a coral-pathogen interaction. Am Nat 170(6): E143–E161.
- 20. Pinzon JH, LaJeunesse T (2011) Species delimitation of common reef corals in the genus Pocillopora using nucleotide sequence phylogenies, population genetics and symbiosis ecology. Mol Ecol 20(2): 311–325.
- 21.
Veron JEN (2000) Corals of the World; Stafford-Smith M, editor. Townsville: Australian Institute of Marine Science. 463 p.
- 22. Ben-Haim Y, Rosenberg E (2002) A novel Vibrio sp. pathogen of the coral Pocillopora damicornis. Mar Biol 141: 47–55.
- 23. Hashimoto K, Shibuno T, Murayama-Kayano E, Tanaka H, Kayano T (2004) Isolation and characterization of stress-responsive genes from the scleractinian coral Pocillopora damicornis. Coral Reefs 23: 485–491.
- 24. Loya Y, Sakai K, Yamazato K, Nakano Y, Sambali R, et al. (2001) Coral bleaching: the winners and the losers. Ecol Lett 4(2): 122–131.
- 25. Stimson J, Sakai K, Sembali H (2002) Interspecific comparison of the symbiotic relationship in corals with high and low rates of bleaching-induced mortality. Coral Reefs V21(4): 409–421.
- 26.
Ben-Haim Rozenblat Y, Rosenberg E (2004) Temperature-regulated bleaching and tissue lysis of Pocillopora damicornis by the novel pathogen Vibrio coralliilyticus. In: Rosenberg E, Loya Y, editors. Coral health and disease. New-York: Spinger-Verlag. pp. 301–324.
- 27.
Dinsdale EA. Abundance of black-band disease on corals from one location on the Great Barrier Reef: a comparison with abundance in the Caribbean region; 2002; Bali. pp. 1239–1243.
- 28. Luna GM, Biavasco F, Danovaro R (2007) Bacteria associated with the rapid tissue necrosis of stony corals. Environ Microbiol 9: 1851–1857.
- 29. Luna GM, Bongiorni L, Gili C, Biavasco F, Danovaro R (2010) Vibrio harveyi as a causative agent of the White Syndrome in tropical stony corals. Environ Microbiol Rep 2(1): 120–127.
- 30. Pollock FJ, Wilson B, Johnson WR, Morris PJ, Willis BL, et al. (2010) Phylogeny of the coral pathogen Vibrio coralliilyticus. Environ Microbiol Rep 2(1): 172–178.
- 31. De O Santos E, Alves N Jr, Dias GM, Mazotto AM, Vermelho A, et al. (2011) Genomic and proteomic analyses of the coral pathogen Vibrio coralliilyticus reveal a diverse virulence repertoire. ISME J 5(9): 1471–1483.
- 32. Kimes NE, Grim CJ, Johnson WR, Hasan NA, Tall BD, et al. (2012) Temperature regulation of virulence factors in the pathogen Vibrio coralliilyticus. ISME J 6: 835–846.
- 33. Vidal-Dupiol J, Ladrière O, Meistertzheim AL, Fouré L, Adjeroud M, et al. (2011) Physiological responses of the scleractinian coral Pocillopora damicornis to bacterial stress from Vibrio coralliilyticus. J Exp Biol 214: 1533–1545.
- 34. Vidal-Dupiol J, Ladrière O, Destoumieux-Garzon D, Sautière P-E, Meistertzheim AL, et al. (2011) Innate immune responses of a scleractinian coral to vibriosis. J Biol Chem 286(25): 22688–22698.
- 35. Ben-Haim Y, Zicherman-Keren M, Rosenberg E (2003) Temperature-regulated bleaching and lysis of the coral Pocillopora damicornis by the novel pathogen Vibrio coralliilyticus. Appl Environ Microbiol 69(7): 4236–4242.
- 36. Flot JF, Tillier S (2007) The mitochondrial genome of Pocillopora (Cnidaria: Scleractinia) contains two variable regions: the putative D-loop and a novel ORF of unknown function. Gene 401(1–2): 80–87.
- 37. Marti-Puig P, Forsman ZH, Haverkort-Yeh RD, Knapp ISS, Maragos JE, et al. (2013) Extreme phenotypic polymorphism in the coral genus Pocillopora; micro-morphology corresponds to mitochondrial groups, while colony morphology does not. Bull Mar Sci 90(1): 1–21.
- 38. Pinzon JH, Sampayo E, Cox E, Chauka LJ, Chen CA, et al. (2013) Blind to morphology: genetics identifies several widespread ecologically common species and few endemics among Indo-Pacific cauliflower corals (Pocillopora, Scleractinia). Journal of Biogeography 40(8): 1595–1608.
- 39. Ben-Haim Y, Thompson FL, Thompson CC, Cnockaert MC, Hoste B, et al. (2003) Vibrio coralliilyticus sp. nov., a temperature-dependent pathogen of the coral Pocillopora damicornis. Int J Syst Evol Microbiol 53(1): 309–315.
- 40. Vidal-Dupiol J, Adjeroud M, Roger E, Foure L, Duval D, et al. (2009) Coral bleaching under thermal stress: putative involvement of host/symbiont recognition mechanisms. BMC Physiol 9: 14.
- 41.
Fox S, Sergei F, Mockler TC (2010) Applications of ultra-hignt-throughput sequencing. In: Belostotky DA, editor. Plant systems Biology. New York: Humana Press. pp. 79–108.
- 42.
Pomraning KR, Smith KM, Bredeweg EL, Phatale PA, L.R C et al.. (2012) Paired-end library preparation for rapid genome sequencing. Fungal Secondary Metabolism. New York: Humana Press.
- 43. Pomraning KR, Smith KM, Freitag M (2009) Genome-wide high throughput analysis of DNA methylation in eukaryotes. Methods 47: 142–150.
- 44. Vidal-Dupiol J, Zoccola D, Tambutté E, Grunau C, Cosseau C, et al. (2013) Genes related to ion-transport and energy production are upregulated in response to CO2-driven pH decrease in corals: New insights from transcriptome analysis. PloS ONE 8(3): e58652.
- 45. Li H, Durbin R (2009) Fast and accurate short read alignment with Burrows-Wheeler transform. Bioinformatics 25(14): 1754–1760.
- 46. Mortazavi A, Williams BA, McCue K, Schaeffer L, Wold B (2008) Mapping and quantifying mammalian transcriptomes by RNA-Seq. Nat Methods 5(7): 621–628.
- 47. Wang L, Feng Z, Wang X, Wang X, Zhang X (2010) DEGseq: an R package for identifying differentially expressed genes from RNA-seq data. Bioinformatics 26(1): 136–138.
- 48.
Bluthgen N, Brand K, Cajavec B, Swat M, Herzel H et al.. (2004) Biological profiling of gene groups utilizing Gene Ontology. Arxiv preprint q-bio/0407034.
- 49. Le QH, Markovic P, Hastings JW, Jovine RVM, Morse D (1997) Structure and organization of the peridinin-chlorophyll a-binding protein gene in Gonyaulax polyedra. Mol Gen Genet 255(6): 595–604.
- 50. Lee DH, Mittag M, Sczekan S, Morse D, Hastings JW (1993) Molecular cloning and genomic organization of a gene for luciferin-binding protein from the dinoflagellate Gonyaulax polyedra. J Biol Chem 268(12): 8842–8850.
- 51. Shoguchi E, Shinzato C, Kawashima T, Gyoja F, Mungpakdee S, et al. (2013) Draft assembly of the Symbiodinium minutum nuclear genome reveals dinoflagellate gene structure. Curr Biol 23(15): 1399–1408.
- 52. Zhang H, Hou Y, Lin S (2006) Isolation and characterization of proliferating cell nuclear antigen from the Dinoflagellate Pfiesteria piscicida. J Eukaryot Microbiol 53(2): 142–150.
- 53. Rudd S (2003) Expressed sequence tags: alternative or complement to whole genome sequences? Trends in Plant Science 8(7): 321–329.
- 54. Rosenberg E, Koren O, Reshef L, Efrony R, Zilber-Rosenberg I (2007) The role of microorganisms in coral health, disease and evolution. Nature Review Microbiology 5(5): 355–362.
- 55. Bosch TCG, Augustin R, Anton-Erxleben F, Fraune S, Hemmrich G, et al. (2009) Uncovering the evolutionary history of innate immunity: The simple metazoan Hydra uses epithelial cells for host defence. Dev Comp Immunol 33(4): 559–569.
- 56. Brower DL, Brower SM, Hayward DC, Ball EE (1997) Molecular evolution of integrins: Genes encoding integrin beta subunits from a coral and a sponge. Proc Natl Acad Sci U S A 94(17): 9182–9187.
- 57. Hayes ML, Eytan RI, Hellberg ME (2010) High amino acid diversity and positive selection at a putative coral immunity gene (tachylectin-2). BMC Evol Biol 10(1): 150.
- 58. Kazumitsu H, Masa-aki K (1983) The effects of lectins on the feeding response in Hydra japonica. Comp Biochem Physiol A-Mol Integr Physiol 76(2): 283–287.
- 59. Knack BA, Iguchi A, Shinzato C, Hayward DC, Ball EE, et al. (2008) Unexpected diversity of cnidarian integrins: expression during coral gastrulation. BMC Evol Biol 8(1): 136.
- 60. Kvennefors ECE, Leggat W, Hoegh-Guldberg O, Degnan BM, Barnes AC (2008) An ancient and variable mannose-binding lectin from the coral Acropora millepora binds both pathogens and symbionts. Dev Comp Immunol 32(12): 1582–1592.
- 61. Miller D, Hemmrich G, Ball E, Hayward D, Khalturin K, et al. (2007) The innate immune repertoire in Cnidaria - ancestral complexity and stochastic gene loss. Genome Biol 8(4): R59.
- 62. Reitzel AM, Sullivan JC, Traylor-Knowles N, Finnerty JR (2008) Genomic survey of candidate stress-response genes in the estuarine anemone Nematostella vectensis. Biol Bull 214(3): 233–254.
- 63. Schwarz J, Brokstein P, Voolstra C, Terry A, Miller D, et al. (2008) Coral life history and symbiosis: Functional genomic resources for two reef building Caribbean corals, Acropora palmata and Montastraea faveolata. BMC Genomics 9(1): 97.
- 64. Shinzato C, Shoguchi E, Kawashima T, Hamada M, Hisata K, et al. (2011) Using the Acropora digitifera genome to understand coral responses to environmental change. Nature 476: 320–323.
- 65. Wood-Charlson EM, Weis VM (2009) The diversity of C-type lectins in the genome of a basal metazoan, Nematostella vectensis. Dev Comp Immunol 33(8): 881–889.
- 66. Putnam NH, Srivastava M, Hellsten U, Dirks B, Chapman J, et al. (2007) Sea anemone genome reveals ancestral eumetazoan gene repertoire and genomic organization. Science 317(5834): 86–94.
- 67. Souter P, Bay LK, Andreakis N, Csaszar N, Seneca FO, et al. (2011) A multilocus, temperature stress-related gene expression profile assay in Acropora millepora, a dominant reef-building coral. Mol Ecol Resour 11(2): 328–334.
- 68. Dishaw L, Smith S, Bigger C (2005) Characterization of a C3-like cDNA in a coral: phylogenetic implications. Immunogenetics 57(7): 535–548.
- 69. Kimura A, Sakaguchi E, Nonaka M (2009) Multi-component complement system of Cnidaria: C3, Bf, and MASP genes expressed in the endodermal tissues of a sea anemone, Nematostella vectensis. Immunobiology 214(3): 165–178.
- 70. Kvennefors ECE, Leggat W, Kerr CC, Ainsworth TD, Hoegh-Guldberg O, et al. (2010) Analysis of evolutionarily conserved innate immune components in coral links immunity and symbiosis. Dev Comp Immunol 34(11): 1219–1229.
- 71. Mydlarz LD, Couch CS, Weil E, Smith G, Harvell CD (2009) Immune defenses of healthy, bleached and diseased Montastraea faveolata during a natural bleaching event. Dis Aquat Org 87(1–2): 67–78.
- 72. Mydlarz LD, Holthouse SF, Peters EC, Harvell CD (2008) Cellular responses in sea fan corals: Granular amoebocytes react to pathogen and climate stressors. PloS ONE 3(3): e1811.
- 73. Mydlarz LD, Palmer CV (2011) The presence of multiple phenoloxidases in Caribbean reef-building corals. Comp Biochem Physiol A-Mol Integr Physiol 159(4): 372–378.
- 74. Palmer CV, Bythell JC, Willis BL (2010) Levels of immunity parameters underpin bleaching and disease susceptibility of reef corals. Faseb J 24(6): 1935–1946.
- 75. Palmer CV, Bythell JC, Willis BL (2011) A comparative study of phenoloxidase activity in diseased and bleached colonies of the coral Acropora millepora. Dev Comp Immunol 35(10): 1098–1101.
- 76. Palmer CV, Bythell JC, Willis BL (2012) Enzyme activity demonstrates multiple pathways of innate immunity in Indo-Pacific anthozoans. Proc R Soc B-Biol Sci 279(1743): 3879–3887.
- 77. Palmer CV, McGinty ES, Cummings DJ, Smith SM, Bartels E, et al. (2011) Patterns of coral ecological immunology: variation in the responses of Caribbean corals to elevated temperature and a pathogen elicitor. J Exp Biol 214(24): 4240–4249.
- 78. Palmer CV, Mydlarz LD, Willis BL (2008) Evidence of an inflammatory-like response in non-normally pigmented tissues of two scleractinian corals. Proc R Soc B-Biol Sci 275(1652): 2687–2693.
- 79. Jung S, Dingley AJ, Augustin R, Anton-Erxleben F, Stanisak M, et al. (2009) Hydramacin-1, Structure and Antibacterial Activity of a Protein from the Basal Metazoan Hydra. J Biol Chem 284(3): 1896–1905.
- 80. Ovchinnikova TV, Balandin SV, Aleshina GM, Tagaev AA, Leonova YF, et al. (2006) Aurelin, a novel antimicrobial peptide from jellyfish Aurelia aurita with structural features of defensins and channel-blocking toxins. Biochem Biophys Res Commun 348(2): 514–523.
- 81. Hemmrich G, Miller DJ, Bosch TCG (2007) The evolution of immunity: a low-life perspective. Trends Immunol 28(10): 449–454.
- 82. Palmer CV, Traylor-Knowles N (2012) Towards an integrated network of coral immune mechanisms. Proc R Soc B-Biol Sci 279(1745): 4106–4114.
- 83. Weiss Y, Foret S, Hayward D, Ainsworth T, King R, et al. (2013) The acute transcriptional response of the coral Acropora millepora to immune challenge: expression of GiMAP/IAN genes links the innate immune responses of corals with those of mammals and plants. BMC Genomics 14(1): 400.
- 84. Burge CA, Mouchka ME, Harvell CD, Roberts S (2013) Immune response of the Caribbean sea fan, Gorgonia ventalina, exposed to an Aplanochytrium parasite as revealed by transcriptome sequencing. Front Invert Physiol 4: 1–9.
- 85. Weis VM (2008) Cellular mechanisms of Cnidarian bleaching: stress causes the collapse of symbiosis. J Exp Biol 211: 3059–3066.
- 86. Bourne DG, Garren M, Work TM, Rosenberg E, Smith GW, et al. (2009) Microbial disease and the coral holobiont. Trends Microbiol 17(12): 554–562.
- 87. Bruno JF, Selig ER, Casey KS, Page CA, Willis BL, et al. (2007) Thermal stress and coral cover as drivers of coral disease outbreaks. PLoS Biol 5(6): e124.
- 88. Lesser MP, Bythell JC, Gates RD, Johnstone RW, Hoegh-Guldberg O (2007) Are infectious diseases really killing corals? Alternative interpretations of the experimental and ecological data. J Exp Mar Biol Ecol 346(1–2): 36–44.
- 89. Muller EM, van Woesik R (2012) Caribbean coral diseases: primary transmission or secondary infection? Glob chang biol 18(12): 3529–3535.
- 90. Mydlarz LD, McGinty ES, Harvell CD (2010) What are the physiological and immunological responses of coral to climate warming and disease? J Exp Biol 213(6): 934–945.
- 91. Duperthuy M, Schmitt P, Garzon E, Caro A, Rosa RD, et al. (2011) Use of OmpU porins for attachment and invasion of Crassostrea gigas immune cells by the oyster pathogen Vibrio splendidus. Proc Natl Acad Sci USA 108(7): 2993–2998.
- 92. Kaniewska P, Campbell PR, Kline DI, Rodriguez-Lanetty M, Miller DJ, et al. (2012) Major cellular and physiological impacts of ocean acidification on a reef-building coral. PloS ONE 7(4): e34659.
- 93. Meyer E, Aglyamova GV, Matz MV (2011) Profiling gene expression responses of coral larvae (Acropora millepora) to elevated temperature and settlement inducers using a novel RNA-Seq procedure. Mol Ecol 20(17): 3599–3616.
- 94. Moya A, Ganot P, Furla P, Sabourault C (2012) The transcriptomic response to thermal stress is immediate, transient and potentiated by ultraviolet radiation in the sea anemone Anemonia viridis. Mol Ecol 21: 1158–1174.
- 95. Moya A, Huisman L, Ball EE, Hayward DC, Grasso LC, et al. (2012) Whole transcriptome analysis of the coral Acropora millepora reveals complex responses to CO2-driven acidification during the initiation of calcification. Mol Ecol 21(10): 2440–2454.
- 96. Rodriguez-Lanetty M, Harii S, Hoegh-Guldberg O (2009) Early molecular responses of coral larvae to hyperthermal stress. Mol Ecol 18(24): 5101–5114.
- 97. Zhao X, Yu H, Kong L, Li Q (2012) Transcriptomic responses to salinity stress in the pacific oyster Crassostrea gigas. PloS ONE 7(9): e46244.
- 98. Barshis DJ, Ladner JT, Oliver TA, Seneca FoO, Traylor-Knowles N, et al. (2013) Genomic basis for coral resilience to climate change. Proc Natl Acad Sci USA 110(4): 1387–1392.
- 99. Gasch AP, Spellman PT, Kao CM, Carmel-Harel O, Eisen MB, et al. (2000) Genomic expression programs in the response of yeast cells to environmental changes. Mol Biol Cell 11(12): 4241–4257.
- 100. Warner JR (1999) The economics of ribosome biosynthesis in yeast. Trends Biochem Sci 24(11): 437–440.
- 101. Vilaprinyo E, Alves R, Sorribas A (2010) Minimization of biosynthetic costs in adaptive gene expression responses of yeast to environmental changes. PLoS Comput Biol 6(2): e1000674.
- 102. Irazoqui JE, Urbach JM, Ausubel FM (2010) Evolution of host innate defence: insights from Caenorhabditis elegans and primitive invertebrates. Nat Rev Immunol 10(1): 47–58.