Deciphering the nature of the coral–Chromera association

Since the discovery of Chromera velia as a novel coral-associated microalga, this organism has attracted interest because of its unique evolutionary position between the photosynthetic dinoflagellates and the parasitic apicomplexans. The nature of the relationship between Chromera and its coral host is controversial. Is it a mutualism, from which both participants benefit, a parasitic relationship, or a chance association? To better understand the interaction, larvae of the common Indo-Pacific reef-building coral Acropora digitifera were experimentally infected with Chromera, and the impact on the host transcriptome was assessed at 4, 12, and 48 h post-infection using Illumina RNA-Seq technology. The transcriptomic response of the coral to Chromera was complex and implies that host immunity is strongly suppressed, and both phagosome maturation and the apoptotic machinery is modified. These responses differ markedly from those described for infection with a competent strain of the coral mutualist Symbiodinium, instead resembling those of vertebrate hosts to parasites and/or pathogens such as Mycobacterium tuberculosis. Consistent with ecological studies suggesting that the association may be accidental, the transcriptional response of A. digitifera larvae leads us to conclude that Chromera could be a coral parasite, commensal, or accidental bystander, but certainly not a beneficial mutualist.


Introduction
Although the association between photosynthetic algae of the genus Symbiodinium and corals has been known for many years, recent work has shown that several lineages of apicomplexans are also associated with corals (e.g., [1][2][3][4]). Ribosomal small subunit RNA (rRNA) sequences previously misidentified as of bacterial origin in fact originate from eight distinct novel apicomplexan-related lineages (ARLs), a subset of which are associated with corals [4]. Two of these ARLs, Chromera velia and Vitrella brassicaformis, are alveolates of the phylum Chromerida [2]. Chromera is the closest known photosynthetic relative of the apicomplexan parasites, but is also related to the photosynthetic dinoflagellates including the coral mutualist Symbiodinium [2]. Chromera was initially isolated from the scleractinian corals Plesiastrea versipora from Sydney harbor and Leptastrea purpurea from One Tree Island on the Great Barrier Reef (GBR), Australia [2]. More recently, Cumbo et al. [5] isolated Chromera from another scleractinian coral, Montipora digitata (Acroporidae), from Nelly Bay, Magnetic Island, in the inner central region of the GBR, and established that this strain of Chromera can infect larvae of both Acropora digitifera and A. tenuis. However, Chromera is clearly not limited to Australian corals, as it has also been isolated from Agaricia agaricites in the Caribbean [6], and identified in sequence data from Curaçao [4,7]. While the photosynthetic capacity of Chromera raises the possibility of it having a beneficial relationship with corals, it is also possible that the interaction may not be mutualistic, given that the majority of apicomplexans are parasites. Little is known about the coral response to pathogens, but, in general, host transcriptional responses to potential mutualists differ markedly from those to parasites [8]. While the innate immune repertoire of corals is surprisingly complex and vertebrate-like [9], to date few studies have addressed immune responses at the transcriptomic level [10][11][12][13][14][15]. However, the response of Acropora to a competent Symbiodinium strain (a strain able to infect coral larvae and establish a mutualistic relationship) is both subtle [16] and transient [17], whereas non-competent strains (that cannot establish a mutualistic relationship) trigger immune responses [16]. The host transcriptional response may, therefore, provide clues as to the nature of the interaction between coral and Chromera.
To investigate the coral response to Chromera infection, gene expression levels in Chromera-infected A. digitifera larvae were compared to uninfected larval controls at 4, 12, and 48 h post-infection by mapping Illumina RNA-Seq reads onto the A. digitifera transcriptome [18]. In contrast to the response to a competent Symbiodinium strain, the response to Chromera resembled that to a non-competent strain in taking place on a longer timescale and involving differential expression of larger numbers of genes. The coral response to Chromera infection was complex, involving modulation of the endocytic pathway and suppression of both immunity and apoptosis, and resembles the vertebrate responses to pathogens and some parasites. The possibility that Chromera has a beneficial relationship with coral hosts requires reevaluation in light of the data presented here.

Materials and methods
Chromera culture Chromera velia CCMP2878 was used in this experiment and was originally isolated from Plesiastrea versipora (Faviidae) from Sydney harbor [2]. Cultures were grown in axenic f/2 medium (G0154, Sigma-Aldrich) [19] and maintained at 25°C under a 12/12-h light/dark cycle before they were used for infection of coral larvae.

Coral larvae and Chromera infection experiment
Approximately 700 planula larvae (6 day post-fertilization) were distributed into each of six 1 L plastic containers containing 700 ml of 0.2 μm FSW, giving three replicates each of uninfected and Chromera-infected larvae. Cultured Chromera were washed three times in 0.2 μm FSW and added at a density of 5 × 10 3 cells/ml. Containers were held at 26°C under fluorescent lamps that provided light (86 ± 2 μmol photon/ m 2 /s measured at the surface) on a 12/12-h light/dark cycle. At 4, 12 and 48 h post-infection,~150 larvae from each replicate were washed in 0.2 μm FSW (ensuring as little liquid carry-over as possible), snap frozen, and stored at −80°C until further treatment. Total RNA from these larvae was isolated as described by Mohamed et al. [17] using TRIzol® reagent (Ambion Life Technologies, Austin, TX, USA) according to manufacturer's instructions. RNA was dissolved in 40 μl of RNase-free water and stored at −80°C. RNA was quantitated and its integrity assessed using a NanoDrop ND-1000 spectrometer (Wilmington, DE, USA) and Agilent 2100 Bioanalyzer (Santa Clara, CA, USA).

Sequencing and RNA-Seq data analysis
RNA-Seq data were obtained and analyzed as previously described by Mohamed et al. [17]. In brief, Illumina reads were mapped onto the A. digitifera transcriptome assembly v1.0 [18] using BOWTIE [20] with default parameters, then transcript counts in each sample were obtained using RSEM software version 1.1.17 [21]. Differential gene expression was inferred based on the mapping counts using edgeR [22]. Chromera-infected samples were compared to control samples at each of the three time points; 4, 12, and 48 h post-infection. P-values for differential gene expression were corrected for multiple testing using the Benjamini and Hochberg algorithm [23], with Pvalue ≤ 0.05 being considered significant. Clustering analysis was conducted on the DEGs detected at 4-h and the 1086 most-highly expressed genes detected at 48-h. Expression values (FPKM) were normalized by the TMMnormalization method (genes were log 2 -transformed and median centered by transcript; Roninson & Oshlack [24]). For DEGs with absolute log 2 (fold-change) > 1, functions were inferred in a two-stage analysis, as previously described [17].
The RNA-Seq data on which this study was based have been submitted to the NCBI Gene Expression Omnibus (GEO) database under Accession number GSE102664.

Microscopic observation of the infection process
When Acropora digitifera larvae were exposed to the CCMP strain of Chromera, intense red fluorescence was observed in the ectoderm of larvae at 4 h post-infection ( Fig. 1). As shown in Fig. 1b, the fluorescence is particulate and quite different from the fluorescence of coral RFP [25,26]. Moreover, no obvious candidates (such as GFP-related or RFP-related proteins) emerged in analyses of early response transcriptomic data (see below). The intense ectodermal fluorescence observed early in the infection process was, however, transient; by 48 h post-infection, only limited fluorescence was visible. Some larvae died during the experiment, but most survived and these contained only small numbers of Chromera (if any) after a few days.

Transcriptome analyses
To understand the nature of the coral-Chromera relationship, host transcriptome-wide gene expression changes were determined following exposure of A. digitifera larvae to Chromera (Fig. 1). At 4, 12, and 48 h post-infection, whole transcriptome expression profiles were compared between Chromera-infected and uninfected coral larvae. Sequencing yielded an average of 20 million paired-end (PE) reads per library, and on average 34% of reads were successfully mapped onto the A. digitifera transcriptome (Supplementary Tables S1 and S2; Supplementary Fig. S1). Hierarchical clustering of expression data revealed a high level of agreement between the biological replicates, as shown in the pairwise Spearman correlations between samples (Supplementary Figs. S2 and S3). Moreover, multidimensional scaling (based on the 500 genes that best differentiate the samples) resolved the Chromera-infected and control samples along the BCV distance 1 (Supplementary Fig. S4).

Differential gene expression analysis
Distinct transcriptomic responses were detected in the coral host at 4 and 48 h post Chromera infection, whereas no significant (using adjusted P ≤ 0.05) differential expression was detected at 12 h. At 4 h, relatively few genes (n = 48; 0.13% of the coral transcriptome) were downregulated ( Fig. 2; Supplementary Fig. S5). However, at 48 h postinfection, the response involved many more genes; 5748 DEGs (about 16% of the coral transcriptome) were differentially expressed, of which 3594 and 2154 genes were downregulated or upregulated, respectively, (

The early response profile
Forty-eight Acropora genes responded to infection at 4 h, whereas many more (5748) responded at 48 h (Supplementary Fig. S9). Twenty of the 48 had reliable Swiss-Prot annotations (BLASTX E-values ≤ 10 −10 ), and amongst these were genes encoding several glycoproteins and members of the immunoglobulin superfamily (Supplementary Table S3). Comparison of the 4 h responses of A. digitifera larvae to Chromera and Symbiodinium [17] revealed only four transcripts in common ( Supplementary  Fig. S10). Amongst these was adi_EST_assem_1403, which encodes a homolog of the major pancreatic secretory granule membrane glycoprotein GP2. adi_ES-T_assem_1403 was downregulated to a similar extent in both Chromera-infected and Symbiodinium-infected larvae (2.72-fold and 2.1-fold, respectively; Supplementary Table S4), suggesting that it may play a central role during coral-algae interactions.

The late (48 h) response profile-GO and KEGG enrichment
Of those genes differentially expressed at 48 h in response to Chromera infection, 4004 had reliable Swiss-Prot annotation (BLASTX E ≤ 10 −10 ) (Fig. 2a). Three pathways were enriched amongst the downregulated genes (Benjamini-corrected P ≤ 0.05)-regulation of actin cytoskeleton, ECM-receptor interaction, and focal adhesion (Supplementary Tables S5-S8). Mapping against the KEGG database provides more detailed views of the enriched pathways (Supplementary Figs. S11-S13).
Downregulated genes at 48 h showed significant GO enrichment with respect to 27 GO-BP (Biological Process), 43 GO-CC (Cellular Component) and 45 GO-MF (Molecular Function) terms (Supplementary Tables S9-S11). Of the genes downregulated at 48 h, GTPase regulator activity was the most overrepresented GO-MF term; also amongst the 122 transcripts in this category were homologs of RAB GTPase activating and binding proteins and members of the TBC1 domain family, all of which play roles in early endosome formation in mammals (Supplementary Table S12, RAB GTPase regulator activity).
Amongst the genes upregulated at 48-h, 14 GO-BP, 63 GO-CC, and 7 GO-MF terms showed significant enrichment (Supplementary Tables S12-S15), the GO-CC term mitochondrion and the GO-MF term structural constituent of ribosome being the most highly overrepresented terms. Several genes encoding proteins involved in ribosome functions and translation were upregulated (Supplementary Table S16, Ribosome). Many genes encoding mitochondrial ribosomal proteins and components of the electron transport chain including many ATP synthase subunits and proteins of the mitochondrial inner membrane complexes were also upregulated (Supplementary Table S17, Mitochondrion).

Focus on genes implicated in host-microbe interactions
The GO and KEGG analyses provided a general overview of pathways enriched during the infection process, after Fig. 2 a Summary of the differential gene expression profile in Chromera-infected larvae compared to controls at 4, 12, and 48 h. An adjusted P ≤ 0.05 and E-value cut off ≤10 −10 were used to filter differentially expressed genes and for BLASTX searches against the Swiss-Prot database. b, c Volcano plots showing the coral genes differentially expressed at 4 h b and 48 h c after Chromera infection compared to the uninfected control condition. The red dots represent the significantly differentially expressed transcripts at adjusted P ≤ 0.05 which another phase of data analysis was undertaken by grouping the annotated DEGs into categories based on literature searches highlighting functions likely to be involved in host-microbe interactions. During the late response to Chromera infection, strong responses were detected in categories of genes associated with immunity, the endocytic pathway, and apoptosis (Fig. 3). The composition of each of these categories of genes is explored in more detail below.

Suppression of the host immune response at 48 h post-infection
Analyses of the late response of the coral to Chromera are consistent with downregulation of host immunity late in the infection process. Sixty A. digitifera genes with potential roles in immunity were downregulated at the late time point, including homologs of pattern recognition receptors (PRRs), components of the complement system, genes involved in the TLR/NLR signaling pathways, antimicrobial activities, and genes involved in ROS and inflammatory responses ( Fig. 4; Table 1 and Supplementary  Table S18).
Candidate PRRs whose expression was downregulated included homologs of mannose receptor 2 (MRC2; a c-type lectin) and lectin domain family 4-member G (CLC4G), as well as two scavenger receptors-the cysteine-rich type 1 protein M160 and CD163 molecule-like 1. Expression of three genes encoding homologs of complement system components C3 and C6 precursors was strongly suppressed. The coral homolog of mammalian GP2, which was downregulated at 4 h post-infection (see above), was further downregulated at the 48 h time point, and genes encoding other glycoproteins, including dystroglycan 1, were also downregulated at this time.
Downregulated genes implicated in TLR and NLR signaling pathways included those encoding homologs of two distinct toll-like receptors and NOD2, as well as TRAF4, NF-kappa-B, and TNF alpha-induced protein 3 (TNFAIP3), which are components of the corresponding downstream signaling pathways.   Fig. 3 Integrative model of the genes and pathways of larvae of Acropora digitifera involved in the interaction with Chromera.Upregulated and downregulated genes/functions are in blue and red text, respectively. The initial contact phase (left panel) involves upregulation of ribosomal and mitochondrial functions (based on enrichment of GO terms), and suppression of host immune responses including the downregulation of the pattern recognition receptors (PRRs) that can recognize signature microbial compounds (i.e., microbe-associated molecular patterns or MAMPs). The PRRs that were detected included toll-like receptors (TLRs), a nucleotide-binding oligomerization domain (Nod) protein, scavenger receptors (SRs), lectins, and the complement protein C3. Moreover, the pancreatic secretory granule membrane major glycoprotein GP2, which serves as an uptake receptor for pathogenic bacteria in man, was strongly downregulated. Suppression of PRR-MAMP interactions leads to inactivation of nuclear factor kappa-B (NF-kappa-B), which is a master regulator of immunity. The phagosome formation phase (central panel) involves downregulation of genes involved in phagocytosis and actin remodeling as well as differential expression of genes implicated in endosomal trafficking that enhance the maturation of the phagosome and lysosome fusion (see Fig. 4 for more details about those genes). The Chromera tolerance phase (right panel) involves complex changes in the apoptotic network. During this phase, genes likely to have both anti-apoptotic and pro-apoptotic functions were differentially expressed Genes encoding proteins implicated in antimicrobial activity including RNA polymerase III polypeptides A and B (which are key players in sensing and limiting bacterial and DNA viral infection), guanylate-binding protein 7 (which promotes oxidative killing and delivers antimicrobial peptides to phagolysosomes), and lactoperoxidase (which uses hydrogen peroxide and thiocyanate to generate the antimicrobial compound hypothiocyanous acid) were all downregulated. Genes involved in ROS and inflammatory responses including those encoding catalase, superoxide dismutase, and arachidonate 5-lipoxygenase were also downregulated ( Fig. 4; Table 1; Supplementary Table S18).   I1  I2  I3  C1  C2  C3   adi_EST_assem_638  adi_EST_assem_3715  adi_EST_assem_35905  adi_EST_assem_8256  adi_EST_assem_14742  adi_EST_assem_10380  adi_EST_assem_1432  adi_EST_assem_5779  adi_EST_assem_7004  adi_EST_assem_2688  adi_EST_assem_3531  adi_EST_assem_8762  adi_EST_assem_3532  adi_EST_assem_4981  adi_EST_assem_461  adi_EST_assem_22790  adi_EST_assem_5267  adi_EST_assem_3877  adi_EST_assem_13491  adi_EST_assem_2273  adi_EST_assem_1695  adi_EST_assem_3655  adi_EST_assem_1190  adi_EST_assem_5314  adi_EST_assem_4678  adi_EST_assem_14075  adi_EST_assem_2070  adi_EST_assem_3633  adi_EST_assem_28915  adi_EST_assem_15402  adi_EST_assem_31512  adi_EST_assem_10795  adi_EST_assem_2961  adi_EST_assem_9676  adi_EST_assem_858  adi_EST_assem_23391  adi_EST_assem_7676  adi_EST_assem_8717  adi_EST_assem_1463  adi_EST_assem_580  adi_EST_assem_17218  adi_EST_assem_4182  adi_EST_assem_5555  adi_EST_assem_26435  adi_EST_assem_1987  adi_EST_assem_10633  adi_EST_assem_3259  adi_EST_assem_8683  adi_EST_assem_3427  adi_EST_assem_5384  adi_EST_assem_9934  adi_EST_assem_18500  adi_EST_assem_26531  adi_EST_assem_6149  adi_EST_assem_9748  adi_EST_assem_4061  adi_EST_assem_4797  adi_EST_assem_683 Table S19), all but one of these being downregulated. Twelve genes implicated in the early stages of phagosome maturation were also downregulated, including those encoding the early endosome markers EEA1 (early endosome antigen1) and ALS2 (amyotrophic lateral sclerosis 2) ( Fig. 5a; Table 2;  Supplementary Table S20); ALS2 is a Rab5 effector protein, which localizes with Rab5 on early endosomal compartments. Conversely, ten genes implicated in the later phase of phagosome maturation were upregulated at the same time point, including those encoding the late phagosome markers Rab7A, CD63, LAMP1 (lysosomal-associated membrane protein 1), TBC1 domain family member 5 (which can displace Rab7A), and LAMTOR2-B, a late endosome/lysosome adapter protein ( Table 2; Fig. 6). Differential expression of genes implicated in autophagy and lysosome function was also observed,   Table S21). In addition, 37 genes encoding proteins involved in endosomal trafficking were differentially expressed including Rab proteins and vacuolar protein sorting proteins (Supplementary Table S22).  Tables S23-S25), but the diversity of roles (likely pro-apoptotic and/or anti-apoptotic functions) complicates interpretation of the biological significance of these data. Eleven genes predicted to have pro-apoptotic functions were downregulated at 48 h ( Fig. 5b; Supplementary Table S23), including genes encoding homologs of the death-domaincontaining protein CRADD (which acts as an apoptotic adapter for caspase-2 by recruiting it to the TNFR-1 signaling complex), FEM-1 (which acts as a death receptor-associated protein, thus mediating apoptosis), and the apoptosome scaffold protein APAF1 ( Fig. 5b; Supplementary Table S23). Consistent with suppression of apoptosis, genes encoding homologs of six mammalian genes with known anti-apoptotic functions ( Fig. 5c; Supplementary Table S24) were upregulated, including a Bcl-2 family member (this A. digitifera gene encodes a clear ortholog of the A. millepora Mcl1-like protein, for which anti-apoptotic activity has been demonstrated; [27]) and the FAS apoptotic inhibitory molecule.

Coral-Chromera initial contact
Conversely, some genes whose products are likely to have anti-apoptotic roles, such as IAP proteins, were downregulated, and others with pro-apoptotic roles upregulated at the same time point (Fig. 5c; Supplementary  Tables S23 and S24). Additionally, 28 genes whose products are implicated in modulation of apoptotic pathway activity, including programmed cell death proteins, were differentially expressed (Supplementary Table S25).

Discussion
Chromera velia was initially isolated from a coral and has been considered a possible mutualist, but its mode of entry into coral larvae is very different to that of Symbiodinium. Although Symbiodinium infection occurs directly into endodermal cells of Acropora larvae [28] and does not occur until after the oral pore has opened, Chromera is able to infect via the ectoderm. Very rapid entry into host cells is also characteristic of parasitic apicomplexans such as Toxoplasma and Plasmodium, where the active process can be effected in 20-30 s [29]. In contrast to these other apicomplexans, little is known about the infective stage of Chromera. However, Oborník et al. [30] have described three distinct morphological stages (coccoid, cystic, and flagellated), of which the flagellated stage seems the best candidate for infection. The flagellated stage develops a unique structure known as the "chromerosome", which Oborník et al. [30] hypothesized might be involved in penetration of coral cells.
Massive invasion of coral planulae by Chromera (i.e., as shown in Fig. 1b) has not previously been documented despite similar experiments having been conducted with Chromera at the same location and with the same host species. However, the initial sampling point used by Cumbo et al. [5] was one day after exposure, by which time

Phagosome maturation
Cathepsin D CD63 CD63 Fig. 6 The coral host modulates the endocytic pathway and enhances phagosome maturation and lysosome fusion in response to Chromera-Upon phagocytosis, the phagosome acquires the GTPase Rab5 via fusion with early endosomes. The Rab5 effector ALS2 was also upregulated. Rab5 acts to recruit phosphatidylinositol-3-kinase (PI3K), which generates phosphatidylinositol-3-phosphate (PI3P) and recruits early endosomal antigen (EEA1) from endosomes. EEA1 is a Rab5 effector protein and its upregulation triggers fusion of the phagosome with a late endosome. During the phagosome maturation process, Rab7 replaces Rab5 and the intermediate phagosome fuses with late endosomal vesicles, acquiring a suite of proteins that includes the proton-ATPase pump (V-ATPase), lysosome-associated membrane glycoprotein 1 (LAMP1), CD63, and lysosomal hydrolases. Vacuoles containing Chromera fuse with late endosomes/lysosomes as indicated by the upregulation of genes encoding Rab7, LAMP1, and CD63 (plus late endosomal/ lysosomal adapter, SNAP-associated protein and vesicleassociated membrane protein 7 (VAMP7) "not shown"). The observed upregulation of the lysosomal V-ATPase signals the formation of phagolysosomes, which are typically rich in hydrolytic enzymes and have an extremely low pH, presumably in order to eliminate Chromera. Genes in red text are downregulated, while those in blue are upregulated fluorescence was greatly reduced in the present study. Another possibility is that the coral response to the CCMP strain of Chromera (used here) differs fundamentally from the response to stains derived from tropical corals (as used in the [5] study). The timing and scale of the transcriptomic response of A. digitifera larvae to Chromera differed markedly from that to a competent Symbiodinium strain [17]. In the latter case, significant changes were detected only at the 4 h time point (although note that a recent study documented subtle transcriptomic changes at 6 d post-infection [31] and involved <3% of transcripts, whereas in the case of Chromera infection, the response occurred predominantly on a longer timescale, but involved many more host genes (>16% of transcripts). While no strictly comparable data are yet available, the coral response to Chromera was more similar to the response to a non-competent strain of Symbiodinium than to a competent strain in that the former provoked more extensive responses on longer time scale [16]. Direct comparison of the data presented here with the literature is complicated by differences in timing -in the present case 4, 12, and 48 h post-infection, whereas Voolstra et al. [16] used 30 min and 6 d sampling times. Some of the same genes were differentially expressed at the late time point in both cases, but the direction of change sometimes differed between the two studies (Supplementary Table S26). Genes implicated in signal transduction, transport, cell adhesion, and recognition were downregulated in response to Chromera, but were upregulated in response to non-competent Symbiodinium, implying that while the coral reacts by attempting to limit interaction with Chromera, in the case of non-competent Symbiodinium a negotiation phase may occur, involving (for example) upregulation of c-type lectin.
Despite the distinct responses of coral larvae to Chromera and Symbiodinium [17], one common characteristic was downregulation of the gene encoding the coral homolog of pancreatic secretory granule membrane major glycoprotein GP2 at 4 h post-infection. In the case of infection with a competent Symbiodinium strain, GP2 expression had returned to baseline levels by 12 h postinfection [17], whereas the gene was further downregulated 48 h after Chromera infection. In man, GP2 is specifically expressed on membranous cells (M cells) where it acts as a bacterial uptake receptor and is involved in initiating an innate immune response [32][33][34]. While the role of coral GP2 protein is unknown, its downregulation in response to the presence of potential symbionts suggests that it might function in host-algal recognition in an analogous manner to its mammalian homologs.

The coral late (48 h) response to Chromera infection
In contrast to the 4 h response, the 48 h response of the coral host to Chromera was extensive and complex. As indicated above, the GO molecular function GTPase regulator activity was highly overrepresented amongst downregulated genes, indicating modulation of the host endocytic pathway, which plays critical roles not only in host-pathogen interactions, such as the mammalian response to Mycobacterium tuberculosis [35], but also in the establishment of symbiosis in corals [17] and sea anemones [36,37]. Amongst upregulated genes, the GO terms mitochondrion and structural constituent of ribosome (ribosomal proteins) were highly overrepresented, implying that Chromera infection resulted in the activation of protein synthesis and metabolism. While responses like this have been reported in the sea fan Gorgonia ventalina after exposure to the parasite Aplanochytrium [38] and in other invertebrates [39,40], these results are in complete contrast to the responses of coral larvae to a competent strain of Symbiodinium, where protein synthesis and metabolism were temporarily suppressed [17].

Modulation of phagosome maturation in the coral host in order to eliminate Chromera
Phagocytosis is well understood in vertebrate phagocytes, such as macrophages, which engulf and destroy non-self cells or microbial invaders. Immediately after phagocytosis, a specific marker, the small GTPase Rab5, is recruited to the phagosome resulting in fusion with an early endosome [41]. Rab5 recruits phosphatidylinositol-3-phosphate (PI3K), which generates phosphoinositol 3-phosphate PI(3)P, which in turn recruits early endosome antigen1 (EEA1) from the endosome. EEA1 is a Rab5 effector protein that triggers fusion of a phagosome with a late endosome [35]. The phago-lysosome is a microbicidal environment that includes antimicrobial peptides and hydrolases. Some pathogenic microbes have evolved strategies enabling them to survive inside host phagocytes (see below), whereas Leishmania mexicana and Coxiella burnetii are well adapted to the highly acidic environment of phagolysosomes and Trypanosoma cruzi has evolved mechanisms enabling escape from the phagosome [35,42,43].
Host endocytic pathway gene expression in response to Chromera infection is very different to that observed with a competent strain of Symbiodinium [17]. During Chromera infection, many genes involved in phagocytosis (e.g., those encoding PRRs) were downregulated, resulting in reduced recognition of invading Chromera and suppression of host immune responses. Downregulation of genes involved in actin remodeling and enrichment of the KEGG pathways regulation of actin cytoskeleton and focal adhesion amongst downregulated genes suggest that actin remodeling is required at the site of phagocytosis in order to prevent the formation and extension of pseudopods. Genes encoding a range of early endosome markers (Rab5 effector proteins (EAA1 and ALS2)) and phosphoinositide-3 and 4-kinases (PI3K and PI4K) were downregulated during Chromera infection ( Table 2). In addition, genes involved in phagosome maturation were upregulated, including markers of the late endosome (Rab7a, LAMP1, VAMPs, and CD63) and of the vesicle fusion process (late endosomal / lysosomal adapter, SNAP-associated protein and VAMP7). Genes encoding two subunits of the lysosomal proton-ATPase were also upregulated, and others encoding vacuolar sorting proteins, TBC1 domain members, other Rab proteins, and lysosomal hydrolases had altered expression compared to the control, indicating the dynamic nature of the endosomal pathway. In contrast to the Chromera data, during infection of coral larvae with a competent strain of Symbiodinium, early endosome markers were upregulated, whereas molecular data for both the coral and a symbiotic sea anemone are consistent with the symbiosome being essentially an arrested early phagosome [17,36,37], these results indicate that the coral host responds to Chromera by upregulating phagosome maturation, presumably to promote fusion with late endosomes and/or lysosomes in order to eliminate Chromera.
In addition to its role in the establishment of a stable cnidarian-Symbiodinium interaction, the endocytic pathway is manipulated in various ways by bacterial pathogens and parasites to gain entry to the cell [44]. For example, like Symbiodinium [17], pathogenic strains of M. tuberculosis [35] arrest phagosome maturation at an early stage, retaining Rab5, but inhibiting both Rab7 acquisition [45] and recruitment of the H + -ATPase to the phagosomal membrane [46]. Non-pathogenic strains of M. tuberculosis do not have the ability to arrest phagosome maturation [35].

Suppression of the host immune response and apoptosis after Chromera infection
Cnidarian genomes encode many of the innate immune components known from vertebrates, including those involved in pattern recognition, signaling cascades, and effector responses [9,47]. At 48 h post-infection, expression of a range of PRRs and other genes, including those encoding C-type lectins, scavenger receptors, complement components, TLRs, NLRs, and GP2 was strongly suppressed, implying that host immunity may be downregulated or compromised in response to Chromera infection. Downregulation of some PRRs could reflect the host attempting to limit interactions with non-beneficial organisms; for example, complement C3 and the C-type lectin, mannose receptor 2 (MRC2), have been implicated in symbiont recognition. Kvennefors et al. [48] localized a C3 homolog near the resident symbionts in the coral A. millepora, suggesting that C3 might act as an opsonin with respect to Symbiodinium and play a role in host-symbiont communication. MRC2 was highly upregulated in A. digitifera larvae infected with a competent Symbiodinium strain [17], suggesting a function in host-symbiont recognition during establishment of coral-algal symbioses. Thus, downregulating the expression of these genes might reflect the host attempting to limit uptake of Chromera. Downregulation of a scavenger receptor might also reflect a response of this type, as a scavenger receptor has been implicated in communication between the symbiotic sea anemone Anthopleura elegantissima and its resident Symbiodinium [49].
On the other hand, some of the observed differences in immune gene expression may reflect suppression of host immunity by Chromera. Downregulation of genes encoding antimicrobial activities (such as guanylate-binding protein 7 and lactoperoxidase) in response to Chromera infection clearly suggests the ability of the invading Chromera to deactivate the host-killing mechanisms, resulting in survival inside the host. Scavenger receptors have been implicated in the recognition of pathogenic Mycobacterium [50], and TLRs and NLRs are involved in sensing the malaria parasite Plasmodium spp. in different vertebrate hosts [51]. Antimicrobial activities are often induced in response to TLR activation [52]. Pathogenic strains of Mycobacterium use the mannose-capped lipoglycan, mannosylated lipoarabinomannan, in order to inhibit TLR signaling and impair T-cell activation.
However, the observed suppression of the immune response by Chromera contrasts with the response of the Caribbean Sea fan, Gorgonia ventalina, to infection with the parasite Aplanochytrium, where genes encoding PRRs and likely antimicrobial activities (including guanylatebinding protein) were upregulated in response to the parasite [38]. The Chromera result also stands in contrast with data from corals infected with white syndrome and white band diseases [11,53]. Moreover, activation of coral immune responses has been reported in response to viral and bacterial mimics [15] and challenge with pathogenic bacteria [54].
The observed downregulation of a suite of pro-apoptotic genes (those encoding TNF receptors and corresponding pathway components, APAF1 and cytochrome c-mediated cell death pathway genes) and upregulation of antiapoptotic genes (Bcl-2, TRAF2, and TRAF5), implies that apoptosis is inhibited in coral larvae following Chromera infection. It is unclear, however, whether these changes reflect host adaptation to tolerate Chromera infection, or manipulation of the coral response by Chromera. Pathogenic microbes including M. tuberculosis have adapted to their host cells by developing anti-apoptosis mechanisms that may reduce the effects of apoptosis-inducing components [55], and a similar strategy appears to be in place in Helicobacter pylori [56].
The relationship between Chromera and coralsparasite or innocent bystander?
The association between apicomplexans and corals has been known for some time (e.g., [1,57]), but the description of Chromera in 2008 [2] and the availability of abundant sequence data from corals have spurred interest in their apicomplexan associates. Kirk et al. [3] used PCR to investigate the presence and seasonal abundance of apicomplexans in four coral species in Florida and the Bahamas for periods of up to nine years and found the associations to be ubiquitous, although seasonably variable. However, rather than Chromera, the most abundant apicomplexan was an undescribed organism presently known only as "genotype N". In an effort to clarify the diversity of apicomplexans associated with coral and coral reefs, Janouškovec et al. [4] examined prokaryotic sequence surveys for DNA from eukaryotic plastids. From these surveys they recognized eight groups of apicomplexan-related lineages, only two of which-Chromera and Vitrellahave been described and characterized. They also suggested that Chromera was associated with the surface of corals rather than with coral tissue, whereas the presently uncharacterized apicomplexan ARL-V was associated specifically with coral tissue rather than with the surface of the coral. That ARL-V was closely associated with the coral while both Chromera and Vitrella were more-closely associated with macroalgae growing near the reef was later confirmed by analysis of data from a fine-scale 16S profiling survey along transects that included a reef and associated macroalgae [7,58]. On the basis of this survey, the authors state "This suggests, contrary to common belief, that Vitrella and Chromera may not be obligate coral symbionts, and possibly interact indirectly with coral."-a conclusion with which the transcriptomic data reported here are consistent. Also, it should be stressed that the absolute abundance of Chromera sequences in these surveys is very low, again suggesting that the association with corals may be somewhat accidental, and is certainly not obligate. However, abundance is not necessarily correlated with functional importance in coral-associated microbes [59].

Conclusion
Despite the level of interest in Chromera, very little is known about the relationship of this novel alga with corals. In the present case, the response of A. digitifera larvae to the CCMP strain of Chromera differed markedly from that to a competent strain of the coral mutualist Symbiodinium, instead being more similar to the response to non-competent ("incompetent") strains of Symbiodinium [16] in extent and timescale (although note that strictly comparable data are not available). Although infection with a competent Symbiodinium strain resulted in arrest of the host endocytic pathway at the early phagosome stage, in the case of Chromera infection, phagosome maturation appeared to be stimulated.
Together with the observed modulation of both the immune and apoptotic networks, these data point to a hostile response to Chromera by coral larvae, and macroscopic observation suggests that Chromera cells are eventually largely cleared despite their initial efficient uptake (Fig. 1); one apparent paradox is that clearing appears to occur more rapidly than can be accounted for by the transcriptional responses documented here. The most likely explanations are that either gene expression changes occurring between the 12 and 48 h sampling time points, or post-transcriptional processes (such as the activation of apoptotic pathways)-or both of these-are responsible for Chromera clearance. Regardless, the abundance of Chromera associated with corals in nature appears to be lowour own attempts to survey their distribution support thispresumably as a consequence of hostile response of the coral to Chromera infection.
While the work described here, and the sequence survey literature, imply that Chromera is not a coral mutualist, the extent to which the observations made here with the CCMP (Sydney Harbor) strain of Chromera hold more generally is not clear. Anecdotally, Chromera is a diverse genus and, in the case of Symbiodinium, high genetic diversity plays out in a spectrum of host responses (for example, [16]). A better understanding of the nature of the coral-Chromera interaction will require estimation of the genetic diversity of Chromera.
AM, EB, and DM interpreted the data and wrote the manuscript. All authors read the article and approved the final version.

Compliance with ethical standards
Conflict of interest The authors declare no conflict of interest.
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as 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 images or other third party material in this article are included in the article's Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons. org/licenses/by/4.0/.