Impact of Ocean Warming and Acidification on Symbiosis Establishment and Gene Expression Profiles in Recruits of Reef Coral Acropora intermedia

The onset of symbiosis and the early development of most broadcast spawning corals play pivotal roles in recruitment success, yet these critical early stages are threatened by multiple stressors. However, molecular mechanisms governing these critical processes under ocean warming and acidification are still poorly understood. The present study investigated the interactive impact of elevated temperature (∼28.0°C and ∼30.5°C) and partial pressure of carbon dioxide (pCO2) (∼600 and ∼1,200 μatm) on early development and the gene expression patterns in juvenile Acropora intermedia over 33 days. The results showed that coral survival was >89% and was unaffected by high temperature, pCO2, or the combined treatment. Notably, high temperature completely arrested successful symbiosis establishment and the budding process, whereas acidification had a negligible effect. Moreover, there was a positive exponential relationship between symbiosis establishment and budding rates (y = 0.0004e6.43x, R = 0.72, P < 0.0001), which indicated the importance of symbiosis in fueling asexual budding. Compared with corals at the control temperature (28°C), those under elevated temperature preferentially harbored Durusdinium spp., despite unsuccessful symbiosis establishment. In addition, compared to the control, 351 and 153 differentially expressed genes were detected in the symbiont and coral host in response to experimental conditions, respectively. In coral host, some genes involved in nutrient transportation and tissue fluorescence were affected by high temperature. In the symbionts, a suite of genes related to cell growth, ribosomal proteins, photosynthesis, and energy production was downregulated under high temperatures, which may have severely hampered successful cell proliferation of the endosymbionts and explains the failure of symbiosis establishment. Therefore, our results suggest that the responses of symbionts to future ocean conditions could play a vital role in shaping successful symbiosis in juvenile coral.


INTRODUCTION
The obligate symbiosis between scleractinian corals and symbionts (family Symbiodiniaceae) is fundamental to the high-level production of reef corals and coral reefs in the oligotrophic tropical oceans. Over the past few decades, however, symbiosis has been affected by many environmental stressors, which can cause coral bleaching and even mortality (Hughes et al., 2017;Van Oppen and Lough, 2018). The main triggers for bleaching events are elevated sea surface temperature (SST) and ocean acidification (Anthony et al., 2008;Hughes et al., 2017). While ocean acidification and warming are known to decrease skeleton development and cause bleaching in adult corals (Foster et al., 2015), their combined effects on juvenile corals are less understood. Juvenile corals have been shown to be extremely sensitive to environmental stress than adults (Fabricius, 2005;Byrne and Przeslawski, 2013). The persistence of coral reefs is partially dependent on the successful dispersal and recruitment of early life phases of coral (Harrison, 2011). Based on the Intergovernmental Panel on Climate Change (IPCC, 2007) report, the SST is expected to increase by 2-3 • C, and surface ocean pH is expected to decline by ∼0.3 by 2100; this will result in a significant threat to juveniles, especially in symbiosis (Hughes et al., 2003;Cziesielski et al., 2019).
The relationship between the stony coral host and Symbiodiniaceae, including systematic identification and symbiont transmission, has been extensively studied. Seven genera (previously nine deeply divergent clades) within Symbiodiniaceae have been formally described based on their genetics and ecology: Symbiodinium (Clade A), Breviolum (Clade B), Cladocopium (Clade C), Durusdinium (Clade D), Effrenium (Clade E), Fugacium (Clade F), and Gerakladium (Clade G). Each genus includes several formally described molecular types (Lajeunesse et al., 2018). Moreover, the majority of reef-building corals (∼80%) establish symbiosis relationships with algal symbionts via horizontal transmission , acquiring symbionts at larval or juvenile stages that are the most suitable for the surrounding settlement environment (Schwarz et al., 1999). Despite the significance of the early life stages to overall coral reproduction, recruitment, and recovery, as well as the sensitivity of these life stages to environmental change, recent studies have primarily focused on the adaptation and response of adult corals to climate change (Silverstein et al., 2015;Davies et al., 2018). Specifically, the role of increased temperature and partial pressure of carbon dioxide (pCO 2 ) on symbiosis establishment in the juvenile stage has not been fully addressed.
It has been found that juvenile and larval corals may possess the ability to select and alter symbiont types and abundance in response to climate change. For example, juvenile Acropora tenuis and Acropora millepora preferentially select Durusdinium, and the relative abundance of Durusdinium spp. increased is high in A. tenuis juveniles under elevated temperature conditions (Abrego et al., 2012;Yorifuji et al., 2017). Similarly, Schnitzler et al. (2012) found that high temperature inhibits the ability of aposymbiotic Fungia scutaria larvae from establishing a symbiosis with Cladocopium spp. Overall, these studies suggest that corals harboring Durusdinium spp. may exhibit high thermal tolerance. Moreover, the acquisition of symbionts by Acropora digitifera juveniles decreased with an increase in seawater acidity (Suwa et al., 2010). Abrego et al. (2008) suggested that A. tenuis juveniles associated with Cladocopium spp. exhibit higher thermal tolerance than those associated with Durusdinium spp., indicating that the relationship between coral endosymbionts and the thermal tolerance of juvenile corals is not fully understood.
In this study, to develop a better understanding of the establishment of symbiosis and driving mechanisms during the onset of juvenile coral-algae symbiosis under climate change, we investigated the effect of elevated temperature (∼28.0 • C and ∼30.5 • C) and acidification (pH ∼8.12 and ∼7.85) for 33 days on the establishment of symbiosis in juvenile A. intermedia. We also measured the physiological parameters of survival, budding, symbiosis establishment, and green fluorescence in juvenile A. intermedia. Gene expression patterns of multiple genes in both coral host and symbionts are representatives of key cellular processes (e.g., energy metabolism and protein synthesis) occurring under different experimental conditions. Our results show the responses of juvenile A. intermedia to possible future thermal and pCO 2 stressors and provide valuable information on the onset of symbiosis and its underlying mechanisms. Furthermore, investigating the effects of elevated temperature and acidification on the onset of symbiosis in juvenile corals will provide new insights into the biological and ecological consequences of global climate change on scleractinian corals.

Coral Collection and Larval Culture
Five colonies of A. intermedia were collected on April 18, 2017, from the Luhuitou fringing reef (18 • 12 N, 109 • 28 E) in Sanya, Hainan Island, China, with the permission of the Administration of the Sanya Coral Reef National Nature Reserve.
The coral samples were immediately transferred to the Tropical Marine Biological Research Station in Hainan and placed in outdoor aquaria with flow-through seawater supply. Gametes were collected from different coral colonies between 21:00 and 21:30 h. Five days post-fertilization (20:00 h, April 23), ∼50-60 larvae were transferred to each Petri dish containing autoclaved chips of Porolithon onkodes, which is recognized as the settlement substrate for Acropora coral larvae (Heyward and Negri, 1999). We allowed the larvae to settle for 24 h and subsequently removed overlapping settlements and non-target recruits, leaving ∼20 recruits in each Petri dish. Additionally, three A. intermedia nubbins from different colonies were preserved in 90% alcohol for downstream symbiont analysis.
Specifically, three to four Petri dishes containing the settled recruits were placed in each flow-through tank (3 L h −1 ) containing 20 L of sand-filtered seawater pumped at a depth of ∼5 m from the nearby reef area. The seawater was mixed and pumped at a rate of 600 L h −1 with a submerged pump (Weipro, China), and the seawater pH was manipulated following a previously described protocol . The seawater pH was manipulated by pH controllers (accuracy: ±0.1 pH, Weipro, China) equipped with solenoid valves to regulate the injection of pure CO 2 and maintain the target pH. In contrast, the control treatment received a supply of outside air. The seawater temperature was maintained using a heater (Weipro, China) (accuracy: ±1 • C). The temperature was monitored at 30-min intervals using a HOBO Pendant, temperature data logger (Onset, Bourne, MA, United States). Lastly, the tanks were maintained under a cycle of 12 h darkness and 12 h light (light intensity ∼200 µmol photons m −2 s −1 ), and the light intensity was ∼200 µmol photons m −2 s −1 . During the 33-day experiment, the recruits were checked daily, and the tanks and tiles were cleaned every 4 days to remove turf algae.
At 5-day intervals, 50 ml seawater samples were collected from each tank and preserved in saturated HgCl 2 for the determination of total alkalinity (TA) with an automatic Total Alkalinity Titrator (AS-ALK2, Apollo, United States). Temperature, salinity, and pH were monitored and adjusted twice daily with a pH meter (SevenGo, Mettler Toledo, Switzerland). The pH probe was calibrated every other day utilizing a two-point calibration technique with 7.00 and 10.00 pH buffers. The carbonate parameters were calculated using CO2SYS (Lewis et al., 1998). The seawater conditions for each treatment are shown in Table 1.
At the end of the experiment, recruits in each Petri dish were checked under a dissecting microscope (SZX7, Olympus, Japan) and scored as alive or dead based on the presence of typical translucent tissue (Abrego et al., 2008). Survivorship of recruits was calculated for each replicate dish as the number of living recruits divided by the initial number of recruits. In addition, the symbiotic status of recruits was defined as follows: (1) aposymbiotic or low abundance of symbionts when the polyp tissue was translucent and lacked pigmentation, with the white skeleton clearly visible under the tissue, and (2) symbiotic when symbiont cells were clearly present and when the oral disk and tentacles of juvenile polyps were brown and pigmented (Abrego et al., 2012). It is important to note that although this visual assessment of symbiont uptake success could not preclude the possibility of symbiont infection, it did represent a good proxy for the establishment of symbiosis, where symbiont cells had not only entered the coral cells but also reproduced successfully in hospite; hereafter referred to as symbiosis establishment (Abrego et al., 2012;Nitschke et al., 2016;Ali et al., 2019). The symbiosis establishment rates were calculated as the proportion of pigmented juveniles relative to the total number of living recruits on the last day, and the budding state of each recruit was determined by recording the occurrence of new buds on the periphery of primary polyps. In addition, using a dissecting microscope, photographs of the polyps were analyzed to determine the budding rates (Humanes et al., 2016).
Furthermore, we found that the tissue of many recruits in the high-temperature treatments exhibited bleaching-induced fluorescence, and therefore, we also determined the percentage of this type of recruit using a dissecting microscope (Roth and Deheyn, 2013). At each census, the number of living recruits was recorded for each treatment. The remaining coral juveniles were preserved in liquid nitrogen and stored at −80 • C for transcriptome analysis.
Statistical analyses were conducted in IBM SPSS. The effects of temperature and pCO 2 on survival, budding, the onset of symbiosis, and green fluorescence were analyzed using a generalized linear model.

Transcriptome Sequencing
Total RNA of each coral juvenile sample was extracted using a TRIzol R Reagent RNA Isolation Kit (Invitrogen, Grand Island, NY, United States) following the manufacturer's instructions. After a quality check, 12 RNA samples were stored at −80 • C and transferred to Novogene Biological Information Technology Co., Ltd. (Beijing, China), for paired-end sequencing. All 12 libraries were sequenced on the Illumina HiSeq X Ten platform. Clean data were obtained by removing sequences containing adapters, reads containing poly-N, and low-quality reads from the raw data. In addition, we calculated the base quality scores (Q20 and Q30) and determined the GC (guanine-cytosine)-content (Supplementary Table 1).
program (Grabherr et al., 2011) with the default setting. A total of 103,285 unigenes were obtained. To identify the unigenes that were derived from either the coral host or the symbiont, we used BLASTN (coverage > 95%, identity > 90%, E < 1E −6 ) to map the database containing cnidarian and Symbiodiniaceae genomes and transcriptomes (Supplementary Table 2). Unigenes were mapped against database containing cnidarian genome and transcriptome using BLASTN, and the obtained sequences were called coral-all unigenes. The coral-all unigenes that could not be mapped to the Symbiodiniaceae database using BLASTN were considered coral unigenes and subsequently named as coral-specific unigenes. In addition, the unigenes were mapped against the database containing Symbiodiniaceae genome and transcriptome using BLASTN (coverage > 95%, identity > 90%, E < 1E −6 ), and the obtained sequences were called symbio-all unigenes. The symbio-all unigenes that could not be mapped to the cnidarian database using BLASTN were considered Symbiodiniaceae unigenes and subsequently named symbio-specific unigenes. The coral-specific and symbio-specific unigenes were used for downstream analysis.

Differential Expression Analysis
The reads per kilobase per million mapped reads (RPKM) value was used to assess unigene expression levels using RESM software (Mortazavi et al., 2008). The RPKM calculation formula is as follows: RPKM = (1,000,000 × C)/(N × L/1,000). Given RPKM (A) to be the expression of unigene A, C is the number of reads that are uniquely mapped to unigene A, N is the total number of reads that are uniquely mapped to all unigenes, and L is the length of unigene A. The RPKM method can eliminate the influence of different gene lengths and amount of sequencing data in the calculation of unigene expression. Additionally, unigenes were annotated by sequence similarity search using public databases: the Swiss-Prot, NCBI Nonredundant protein (Nr), and KOG databases using BLASTX, BLASTN, and BLASTP, respectively, with E value < 10 −5 . In addition, three comparisons were made to examine the changes in coral and symbiont gene expression between the treatments ATHC, HTAC, and HTHC and the control (ATAC): ATHC treatment vs. control; HTAC treatment vs. control, and HTHC treatment vs. control. Subsequently, we utilized the R package, DESeq2, to identify differentially expressed genes (DEGs) (Mohamed et al., 2019) with thresholds of log 2 (foldchange) > 2 and a false discovery rate (FDR) < 0.01.
Using a model based on the negative binomial distribution, DESeq2 provided statistical routines for determining differential expression in the digital gene expression data.

Identification of Symbiodiniaceae
To identify the Symbiodiniaceae genera present in juvenile A. intermedia, we conducted a local BLAST search of the symbiont transcriptome against a database containing chloroplast 23S rDNA and nuclear 28S rDNA sequences (Pochon et al., 2006). In contrast, to identify the Symbiodiniaceae genera present in adult A. intermedia, we performed highthroughput sequencing of the internal transcribed spacer 2 (ITS2) region of ribosomal RNA gene amplicons, following the protocols described by Zhou et al. (2017). Total DNA was extracted as described previously (Zhou et al., 2012). The purity of the DNA was measured using a NanoDrop spectrophotometer (Thermo Fisher Scientific, United States). The DNA samples were stored at −20 • C for further analysis. For library preparation, the DNA samples were amplified using polymerase chain reaction (PCR) primers employed by Lajeunesse and Trench (2000). The ITS2 amplicon was approximately 320 base pairs long. The primer sequences were 5 -AGATCGGAAGAGCGTCGTGTAGGGAAAGAGTGT N 6 GAA TTGCAGAACTCCGTG-3 (ITSintfor2) and 5 -AGATCGGA AGAGCACACGTCTGAACTCCAGTCACN 6 GGGATCCATAT GCTTAAGTTCAGCGGGT-3 (ITS2-reverse), targeting the ITS2 region of the ribosomal RNA gene for Symbiodiniaceae. Primers included Illumina library adapters (underlined) and a bar code (shown as N). Then, PCR was performed with 12.5 µl of PCR reagent (Bio-Rad, United States), 0.1 µM primer, 50 ng of DNA, and DNase-free water to make a total volume of 25 µl.
The following procedure was implemented: initial denaturation for 3 min at 94 • C, followed by 34 cycles at 98 • C for 10 s, 51 • C for 30 s, 68 • C for 30 s, and a final extension step of 5 min at 68 • C. The PCR products were validated using an Agilent 2100 Bioanalyzer (Agilent Technologies, Palo Alto, CA, United States) and quantified with a NanoDrop spectrophotometer. All qualified amplification products were mixed in equal amounts followed by sequencing on an Illumina MiSeq instrument (Illumina, San Diego, CA, United States) according to the manufacturer's instructions using a 300 × 2 paired-end configuration. The raw data were submitted to the National Center for Biotechnology Information (NCBI) Sequence Read Archive (PRJNA633301). Processing and analysis of the ITS2 sequences were conducted following our previously described protocol (Gong et al., 2018).

Response of Coral Juveniles to Elevated Temperature and Partial Pressure of Carbon Dioxide
Higher temperatures arrested symbiosis establishment (F 1,28 = 248.11, P < 0.001) and reduced budding in A. intermedia juveniles (F 1,28 = 8.868, P = 0.006), while the effects of reduction in pH were not significant (F 2,28 = 2.533, P = 0.123). The mean symbiosis establishment rates for ambient temperature treatments were ∼70% but 0 for high-temperature treatments ( Figure 1B). In addition, 13.48% and 4.09% of recruits successfully budded under the ATAC and ATHC treatments, respectively. In contrast, no recruits budded successfully under higher temperature treatments ( Figure 1C). Moreover, at the end of the 33-day experimental period, an exponential relationship (y = 0.0004e 6.43x , R = 0.72, P < 0.0001) existed between symbiosis establishment and budding rates (Figure 2A). All the data points in Figure 2A are from individual Petri dishes of the four treatments. The aposymbiotic polyp, symbiosis establishment, and budding are shown in Figures 2Ba-c, respectively.
The survival rates of coral recruits ranged from ∼89 to 96% across all treatments (Figure 1A), and the statistical results (Supplementary Table 3) show that the survival rates did not differ significantly between treatments. However, higher temperature significantly elevated green fluorescence in tissues (F 1,28 = 217.715, P < 0.001). The percentage of tissue expressing green fluorescence under ambient temperature treatments was ∼4%, but under high-temperature treatments, it was ∼50% (Figures 1D, 2B, d). At high temperatures, the number of corals exhibiting vibrantly green fluorescence increased 10-fold.

Assembly of the Reference Transcriptome
The obtained coral holobiont transcriptome was 88.19 Mb with an N50 size of 1,522 bp. The maximum, minimum, and mean lengths of the unigenes were 21,321,201, and 853 bp, respectively. The GC-content was 41.60%. The N50 of the coral transcriptome was 1,085 bp, with a mean unigene length of 697 bp, and a GCcontent of 39.90%. Lastly, the N50 of the symbiont transcriptome was 472 bp, with a mean unigene length of 434 bp and a GC-content of 57.11% (Supplementary Table 4). The raw data  have been submitted to the NCBI under accession numbers SRR10305223-SRR10305234.

Differential Gene Expression
Our analyses revealed 154 DEGs within the coral host. Specifically, under the ATHC, HTAC, and HTHC treatments, 3, 77, and 74 genes, respectively, were differentially expressed in the coral host. For the symbiont, 351 genes were differentially expressed (0, 133, and 218 under the ATHC, HTAC, and HTHC treatments, respectively) (Supplementary Figure 1). Although the number of DEGs increased under HTAC and HTHC treatments, the number of downregulated DEGs far exceeded that of upregulated ones.

Differentially Expressed Gene Functions
Unigenes were categorized based on functional similarity, while the search processes were based on sequence similarity. Unigenes were annotated by BLASTX against the Swiss-Prot, BLASTN against the Nr, and BLASTP against the KOG databases with E value < 10 −5 . The function of unigenes was grouped according to the public database search and literature (Desalvo et al., 2008). As a result, the functional categories of stress, transport protein, lipid or fatty acid metabolism, cell proliferation, fluorescent proteins, signaling, transcription factors, cell immune response, and others were related to the coral host (Figure 3 and Supplementary Figure 2). Using the same method, the DEG functions of the symbionts were also classified into 16 functional groups: stress, photosynthesis, Calvin cycle, glycolysis, gluconeogenesis, tricarboxylic acid (TCA) cycle, lipid or fatty acid metabolism, respiratory chain, protein metabolism, amino acid and vitamin metabolism, DNA replication, transcription, proton ion protein transportation, transport protein, flagellum-related factors, and others (Figure 4 and Supplementary Figure 2).

Identification of Symbiodiniaceae Communities
In juvenile A. intermedia, there were significant differences in the symbiont communities between treatments, suggesting that coral-algal symbiosis is highly flexible. Notably, these differences were attributed to changes in the relative abundance of existing Durusdinium spp. under high temperatures (Figure 5). The Symbiodiniaceae associated with the ATAC and ATHC treatments were mainly Symbiodinium spp. Of three replicates of the HTAC treatment, two were occupied by Durusdinium spp., whereas the other was not detected. At the subclade level, D1, D1a and D1.1, and D1.2 dominated in the HTAC treatment. Of three replicates of the HTHC treatment, two were occupied by both Durusdinium spp. and Symbiodinium spp., while the other was not detected. At the subclade level, the abundance of each of D1, D1a, A2-1, A2-2, and A3 was higher. In contrast, C3 and D1 were the most abundant Symbiodiniaceae species in adult A. intermedia (Figure 5).

DISCUSSION
Understanding the effects of global climate change on coral requires information not only concerning its the direct impact on physiological indicators but also on gene expression levels.
In this study, physiological responses and a large number of differentially expressed transcripts were observed in juvenile corals in the experimental conditions. Firstly, all aposymbiotic juveniles in the four treatments were of the same status before culturing. Secondly, at the end of the 33-day experiment, the juveniles had already shown different responses (e.g., symbiont community). As a result, all the different responses are likely attributable to the experimental conditions.
Ocean warming might be more deleterious for symbiosis establishment than acidification. Elevated temperature may not FIGURE 3 | Heatmap of annotated unigenes in the coral host in the ATAC (28 • C and pCO 2 ∼600 µatm), ATHC (28 • C and pCO 2 ∼1,200 µatm), HTAC (30.5 • C and pCO 2 ∼600 µatm), and HTHC (30.5 • C and pCO 2 ∼1,200 µatm) treatments. Note: numbers 1-7 represent fluorescent proteins, transport protein, cell proliferation, immunity-related factors, signal transduction, lipid or fatty acid metabolism, and stress category, respectively. The heatmap was drawn based on the log-transformed RPKM value of unigenes. The color bar in the heatmap represents the variation of the value of RPKM (from white to blue means the expression value of unigene from 0 to 0.4). High expression (blue color) represents the high value of RPKM, which means that the expression level of unigene among different groups or different unigenes is high. Unigenes were functionally categorized based on database searches and literature.
favor symbiosis onset, and this could affect asexual budding reproduction. Temperature stress reduces ribosomal protein gene expression in the symbionts and upregulates fluorescent protein gene expression in the coral host. This study yielded novel evidence of symbiont reproduction and coral host stress response under heat stress conditions. The results may provide a more comprehensive understanding of the influences of elevated temperature and pCO 2 on symbiosis establishment in juvenile corals.

Response of Symbiodiniaceae to Elevated Temperature
Although elevated temperature impeded the onset of symbiosis, there are symbionts in the juvenile corals based on the RNAseq data under ATAC and ATHC treatments, suggesting that the method of using a dissecting microscope to detect the symbiont infection may not be very accurate. Under ambient temperature conditions, the main components of the Symbiodiniaceae found in the recruits were Symbiodinium spp., but in response to thermal stress, Durusdinium spp. became predominant, indicating that juveniles are likely to become locally adapted through shifts in symbiont community composition. Juveniles possibly adapt to a high temperature by forming novel host-Symbiodiniaceae combinations and taking on symbionts that can tolerate warm conditions. Mounting evidence suggests that scleractinian corals may be able to choose algal symbiont types in response to increasing seawater temperatures (Abrego et al., 2009;Yorifuji et al., 2017) and ocean acidification . For adult corals, the dominant symbionts were Cladocopium and Durusdinium spp. Distinct symbiont communities were detected in juvenile and adult A. intermedia, suggesting that resilience to thermal stress is associated with a highly flexible symbiotic association.
Most of the downregulated genes identified in the symbionts were related to the encoding for ribosomal proteins (Supplementary Figure 3). Ribosomal proteins of Montastraea faveolata are downregulated in response to heat shock (Voolstra et al., 2009;Polato et al., 2010). In contrast, Kenkel et al. (2017) suggested that the upregulation of ribosomal proteins in symbionts resulted in increased cell growth and, thus, elevated net photosynthesis. Ribosomal function is associated with protein synthesis and, thus, cell growth through regulation of the cell cycle and cell size (Jorgensen and Tyers, 2004). Moreover, under high-temperature conditions, DNA replication-related genes were downregulated, probably further contributing to reduced symbiont cell growth and reproduction. Thus, the reduction in symbiont cell growth and development may have inhibited the successful establishment of symbiosis under high temperatures.
The transcriptome data revealed that the carbohydrate metabolism of the symbionts decreased because of high temperature. Some genes related to enzymes of photosynthesis, glycolysis, and TCA cycle were downregulated under hightemperature treatments (Figure 4). Regarding photosynthesis, FIGURE 4 | Heatmap of annotated unigenes related to photosynthesis (A) and metabolic processes (B) in the symbionts under the ATAC (28 • C and pCO 2 ∼600 µatm), ATHC (28 • C and pCO 2 ∼1,200 µatm), HTAC (30.5 • C and pCO 2 ∼600 µatm), and HTHC (30.5 • C and pCO 2 ∼1,200 µatm) treatments. Number 1-7 represents Calvin cycle, glycolysis, gluconeogenesis, TCA cycle, lipid or fatty acid metabolism, respiration chain and ATP synthesis category, respectively. The heatmap was drawn based on the log-transformed RPKM value of unigenes. The color bar in the heatmap represents the variation of the value of RPKM (from white to blue means the expression value of unigene from 0 to 1). High expression (blue color) represents the high value of RPKM, which means that the expression level of unigene among different groups or different unigenes is high. Unigenes were functionally categorized based on database searches and literature.
when symbionts are exposed to elevated temperatures, some critical proteins (e.g., light-harvesting protein and fucoxanthin-chlorophyll a-c binding protein F) are damaged, and this may result in the dysfunction of photosystem II (Warner et al., 1999). Still, heat stress may damage the thylakoid membranes, causing them to less produce ATP and nicotinamide adenine dinucleotide phosphate (NADPH) (Weis, 2008). Therefore, this result suggests that prolonged exposure to heat stress may cause symbiont photoinhibition, resulting in reduced carbon translocation from the symbiont to the coral host (Hillyer et al., 2017b). Moreover, pyruvate kinase (limited enzyme in glycolysis) is downregulated under high temperatures, less producing pyruvate ATP (Fernie et al., 2004). Pyruvate is then fed into the TCA cycle to generate less ATP. Regarding the reduction of glycolysis and photoinhibition, similar effects on symbionts are produced by marine pollutants such as 1,3,5-trinitro-1,3,5 triazine (Gust et al., 2014). A gene (unigene0050751) that encoded an essential enzyme [succinate dehydrogenase (ubiquinone) flavoprotein subunit, mitochondrial] within the TCA cycle was significantly downregulated [log 2 (foldchange) = −20.09] in the symbiont transcriptome. Therefore, there was a probable reduction in ATP production given that symbiont photoinhibition has been previously reported to be the result of a reduction in symbiont carbohydrate metabolism (Lesser, 1997).
Lastly, regarding the symbiont DEGs, it has been confirmed that the differentially expressed transcripts of C. goreaui and D. trenchii in juveniles were different (Yuyama et al., 2018). However, the determination of experimental conditions or symbiont species that account for these DEGs requires further investigation, since it cannot be addressed from the data obtained in the study -as in a previous study (Davies et al., 2018).

Survival, Symbiosis Establishment, and Asexual Budding
Our finding that the survivorship of juvenile A. intermedia remained ∼96% after 33-day exposure to high temperature is in agreement with results of Foster et al. (2015), who found no significant effect of elevated temperature on the post-settlement survival rate of juvenile A. spicifera. In the present study, the most plausible explanation for the lack of an elevated temperature effect on survivorship is unsuccessful symbiont infection under thermal stress. Yakovleva et al. (2009) suggested that symbionts can cause an increase in oxidative stress in coral larvae under heat stress, resulting in increased reactive oxygen species (ROS) and tissue damage, leading to reduced larval survival. Therefore, juvenile corals under high temperatures would be more resistant to temperature stress due to lower ROS production, thus favoring their survival. However, budding and symbiosis onset were significantly affected by high temperature. Budding is likely to be an energetically expensive process for newly settled coral and potentially strongly dependent on the successful establishment of symbiosis (Graham et al., 2013). Indeed, budding and symbiosis establishment rates formed a linear relationship (y = 0.0004e 6.43x , R = 0.72, P < 0.0001) in A. intermedia, possibly indicating that symbiosis establishment is positively correlated with budding in the coral host and that budding potentially benefits from symbiosis establishment.

Response of Fluorescent Protein Genes in Coral Host to Elevated Temperature
In this study, there was an increase in green fluorescent tissue in juvenile coral, and two genes (unigene103324 and 0000505) encoding fluorescent proteins (e.g., TPA: fluorescent protein 6, partial, and cyan fluorescent protein 2) were indeed upregulated 4.5-7.7-fold under elevated temperature conditions. Similarly, Yuyama et al. (2012) revealed that juvenile A. tenuis dominated by Durusdinium spp. highly expressed green fluorescent proteins (GFPs) under heat stress. Fluorescent proteins in corals exhibit ROS-scavenging ability and are generally regarded as one of the host-driven protective mechanisms that could contribute to the regulation of the bleaching response (Palmer et al., 2009;Bollati et al., 2020). The phenomenon of green fluorescence is highly expressed in response to high temperature in adult corals owing to a family of GFP-like proteins and is called colorful bleaching (Bollati et al., 2020). Thus, the overexpression of fluorescent proteins in coral that hosted Durusdinium spp. may enhance resistance and persistence under heat stress.

Response of Transport Protein Genes in Coral Host to Elevated Temperature
Several sets of genes encoding transport protein were downregulated in high temperatures, such as Tret1, SLC7A1, and SLC7A2. Tret1 codes for a "PREDICTED: solute carrier family 2, facilitated glucose transporter member 8-like" protein, while SLC7A1 and SLC7A2 code for "PREDICTED: cationic amino acid transporter 2-like isoform X1" and "PREDICTED: cationic amino acid transporter 2-like isoform, " respectively. These Red and green color scale represents upregulation and downregulation, respectively. AA, amino acid; Chl, chloroplast; Glc, glucose; mt, mitochondria; PA, pyruvate; TCA, tricarboxylic acid cycle. Italic letters around unigenes are annotation genes (i.e., RPL7A), and Arial letters are annotation proteins (i.e., cyan fluorescent protein). Unigenes in the schematic are selected by key genes, differential expression (up/down), and the value of RPKM (widely varied). The color variations in green (down expression) and orange (up expression) bar mean that the values of RPKM range from 0 to 10, 11 to 100, 101 to 1,000, and >1,000. Each gene box is equally divided into four pieces (A, C, H, and H+C), sequentially representing ATAC, ATHC, HTAC, and HTHC treatment, respectively. two amino acid transporters translocate L-arginine produced by symbiont algae from the symbiont to the coral host (Ferrier, 1994;Cérec et al., 2007), although glucose is the major photosynthetic product translocated from symbionts to coral host (Hillyer et al., 2017a). The onset of symbiosis was not successfully established under heat stress, suggesting a low production of glucose and amino acids. Our transcriptomic analyses did show that glycogen biosynthesis was downregulated in symbionts under thermal stress. Therefore, the overall pattern of low expression of Tret1, SLC7A1, and SLC7A2 may correspond to the reduced availability of glucose and amino acids during thermal stress. This low expression of transporter protein genes may simply reflect the dynamic interaction and coordination between coral host and symbionts under stress.

Metabolic Integration of Coral and Symbiodiniaceae
Our results provide molecular evidence that the number of symbionts in coral host is low under heat stress and, thus, support a model of metabolic integration, such as the one presented in Figure 6. In brief, the fact that symbionts cannot successfully establish symbiosis with coral host under high temperature may result from downregulation of several genes of ribosomal and DNA replication proteins in the symbionts, leading to low algal density. Meanwhile, heat stress negatively affects some vital enzymes or proteins in metabolic processes such as TCA, Calvin cycle, glycolysis, and photosynthesis. Due to the low algal population, the amount of fixed organic carbon obtained by the coral host is low. As a result, the genes related to glucose and amino acid transporters were downregulated. We hypothesize that these changes in symbiont metabolism would, in turn, result in fewer inorganic nutrients for the symbionts. In summary, the physiological and transcriptomic data demonstrated that ocean warming might be more deleterious to juvenile coral than ocean acidification.

DATA AVAILABILITY STATEMENT
The datasets generated for this study can be found in the all raw data from this study are available in the NCBI under the accession number SRR11802619-SRR11802621 and PRJNA633301.

AUTHOR CONTRIBUTIONS
YS, LJ, and HH conceived and designed the study. YS, LJ, SG, and XY conducted the experiments and performed the laboratory and transcriptomic analyses. YS analyzed the data and drafted the manuscript. LJ, SG, XL, YZ, MG, XY, GZ, JL, PQ, and HH contributed to the laboratory analyses and interpretation of the results. All authors contributed to the article and approved the submitted version.

ACKNOWLEDGMENTS
We thank the staff of the Tropical Marine Biological Research Station in Hainan for providing technical assistance and facilities for conducting the study.
Supplementary Figure 3 | Heatmap of annotated unigenes related to protein synthesis in the symbionts in the ATAC (28 • C and pCO 2 : ∼600 µatm), ATHC (28 • C and pCO 2 : ∼1200 µatm), HTAC (30.5 • C and pCO 2 : ∼600 µatm) and HTHC (30.5 • C and pCO 2 : ∼1200 µatm) treatments. The heatmap was drawn based on log-transformed RPKM value of unigenes. The color bar in heatmap represents the variation of the value of RPKM (from white to blue means the expression value of unigene from 0 to 0.4). High expression (blue color) represents the high value of RPKM, which means the expression level of unigene among different groups or different unigenes is high. Unigenes were functionally categorized based on database searches and literature.
Supplementary Table 2 | List of genomic and transcriptomic resources downloaded and used to filter the Acropora intermedia holobiont transcriptome into coral and symbiont compartments.
Supplementary Table 3 | Statistical results of generalized linear models and post-hoc comparisons (P-values) examining the effects of temperature (T) and pCO 2 on the survival, budding, symbiosis establishment, and green fluorescence of Acropora intermedia.