Survivorship and growth
A. cf. hyacinthus survivorship after 12 months was not statistically different amongst coral groups (wild vs nursery grown corals, 100.0 ± 0.0 vs 77.8 ± 11.1 %, respectively; n = 6 for each coral group) (Wilcoxon rank-sum test, p = 0.2500, Fig. 2A). However, higher relative areal growth rate (ΔGA; % growth in cm2 y-1) over 12 months was evident for the nursery corals (308.9 ± 52.8) compared to the wild corals (12.4 ± 11.5) (t(5.0) = 4.2, p = 0.0083, Fig. 2B). This outcome was also observed in the relative increase in area covered by coral tissue (%), (112.4 ± 11.5 vs 309.0 ± 54.0 %, for wild vs nursery corals, respectively; n = 6 for each coral group) (Fig. 2C; t(5.0) = 4.7, p = 0.0055). Similarly, relative linear growth rate (ΔGL; % growth in cm y-1) was higher for nursery compared to wild corals (84.4 ± 11.2 vs 4.5 ± 5.0; t(5.0) = 4.3, p = 0.0076, Fig. S3A). Amongst the six colonies examined, highest areal growth rates in wild and nursery corals were consistently observed for colony #1 (50.7 and 432.0 cm2 y-1, respectively; see example in Fig. S3B-C). Overall, areal and linear growth rates observed for nursery corals were ~20-25 times higher compared against wild corals.
Physical appearance
Multiple traits were examined to describe bio-physical colouration (red [R], green [G], blue [B], hue, saturation, brightness, and bleaching state; Fig. S4). Initial RGB values for wild corals (WT0) were 103.2 ± 10.0, 89.0 ± 8.5, and 68.0 ± 9.8 for RGB, respectively (n = 6). Higher RGB and brightness values correspond to a paler coloration (i.e. closer to white). R and G values were slightly increased (138.5 ± 6.0 and 97.8 ± 3.4), and B values decreased (42.2 ± 3.9) in wild corals (n = 6) after 12 months (WT12). In contrast, the nursery corals (NT12) became darker after 12 months (RGB of 98.5 ± 5.1, 67.3 ± 4.6, and 27.8 ± 2.6, respectively; n = 6) compared to the other two coral groups. Differences in the RGB, and brightness values among WT12 and NT12 coral groups were evident (see RM-ANOVA in Fig. S4A-D) indicating different colouration between these two coral groups. We also observed similar saturation and hue levels among WT12 and NT12 (see RM-ANOVA in Fig. S4E-F). All coral groups were in a non-bleached stage (average score of 6.0 ± 0.0 for all three coral groups) based on the coral colour reference card42 (See Fig. S4G). In summary, at 12 months post-fragmentation, corals propagated within the nursery (NT12) were darker than those retained in the wild (WT12) compared to the start (WT0).
Photobiology
Whilst nursery-propagation resulted in darker colonies after 12 months, photobiological traits for NT12 vs WT12 were generally similar; of note, however, photobiological characteristics varied as a function of time whereby both NT12 and WT12 were largely different compared to WT0 (for 8 out of 12 photobiological traits; EK, Fq´/Fm´MAX, [1-C]MAX, [1-Q]MAX, symbiont cell density, and total pigment, Chl a, and Chl c2 density). Specifically, the light saturation coefficient (EK, μmol photons m–2 s–1) was higher for WT12 and NT12 (292.8 ± 64.4 and 328.2 ± 49.7, respectively) compared to WT0 (126.4 ± 16.9) (Fig. 3A; one-way RM ANOVA, F(1.4,7.2) = 6.2, p = 0.0337; Tukey’s test WT0 vs WT12, and WT0 vs NT12; p = 0.0447, and 0.0310, respectively). Maximum photochemical efficiency of PSII (Fq´/Fm´MAX, dimensionless) was lower for WT12 and NT12 (0.6412 ± 0.0130 and 0.6368 ± 0.0162) compared to WT0 (0.7043 ± 0.0118) (Fig. 3B; one-way RM ANOVA, F(1.3,6.5) = 12.2, p = 0.0091; Tukey’s test WT0 vs WT12, and WT0 vs NT12, p = 0.0134, and 0.0438, respectively). When examining the fluorescence quenching parameters, whereby photochemical quenching capacity was reduced (i.e. [1-C]MAX values were higher) for WT12 and NT12 (0.2330 ± 0.0422 and 0.3037 ± 0.0342) compared to WT0 (0.1214 ± 0.0198) (one-way RM ANOVA, F(2.0,9.8) = 16.8, p = 0.0007; Tukey’s test WT0 vs WT12 (p = 0.0350) and WT0 vs NT12 (p = 0.0070); Fig. 3C). Similarly, non-photochemical quenching capacity was also reduced ([1-Q]MAX values were higher) for WT12 and NT12 (0.8539 ± 0.0345 and 0.8123 ± 0.0387) vs WT0 (0.6879 ± 0.0255) (Fig. 3D; one-way RM ANOVA, F(1.7,8.6) = 7.1, p = 0.0169; Tukey’s test WT0 vs WT12, p = 0.0488).
Symbiont cell densities (cells x 106 cm-2) were higher for WT12 and NT12 (0.5600 ± 0.1037 and 0.7317 ± 0.1224, n = 6, respectively) compared to the initial wild colonies (WT0) (0.2525 ± 0.0368, n = 4) (Fig. 3E; one-way ANOVA, F(2,13) = 2.5, p = 0.0331; Tukey’s test WT0 vs NT12; p = 0.0264). No differences were found among coral groups for total pigment per cell, although lowest values were observed for NT12 (8.7 ± 1.1), followed by WT12 (10.5 ± 2.2) and highest values for WT0 (13.4 ± 2.2) (Fig. 3F-H; one-way ANOVA, F(2,13) = 1.345, p = 0.2945), consistent with higher light-acclimated symbiont photobiology (above). However, pigment per coral surface area demonstrated that NT12 (but not WT12) A. cf. hyacinthus was more pigmented than WT0 (Fig. 3I-K), consistent with the darker coloration observed for the NT12 corals (Fig. S4). Specifically, total pigment, and Chl a and Chl c2 densities (per coral surface area) were higher for NT12 (5.9 ± 0.7, 4.2 ± 0.6, and 1.7 ± 0.2, respectively) vs WT0 (3.2 ± 0.3, 2.2 ± 0.2, and 1.0 ± 0.1, respectively; Tukey’s test, p = 0.0117, 0.0275, and 0.0026, respectively). Intriguingly, the chlorophyll ratio (a:c2) remained similar over time and environments (WT0 (2.3 ± 0.2), WT12 (2.8 ±0.2) and NT12 (2.6 ± 0.3); one-way ANOVA, F(2,13) = 0.3, p = 0.3543, Fig. 3L). In summary, after 12 months of growth, colonies on both the nursery and reef contained more symbiont cells (and more pigments) (typical of lower-light adapted corals) but were counterintuitively acclimated to higher light intensities.
Symbiodiniaceae ITS2 identity
All samples showed an association with Cladocopium; however, the major ITS2 type profiles were highly variable across all groups (WT0, WT12, and NT12; see Fig. S5), with only one major ITS2 type profile (C50a/C3k/C50c/C3-C3b-C50f) present in all three groups. All four initial samples had unique major ITS2 type profiles that included C21/C3, C40/3, C3-C21-C3k-C3at-C3b-C3av-C3dp, and C50a/C3k/C50c/C3-C3b-C50f. The wild colonies after 12 months also all had distinct major ITS2 type profiles. The highest association of WT12 (3 out of 6 samples) was with Cladocopium of the C3k radiation, albeit three distinct C3k ITS2-type profiles observed (representative of different genotypes). Of all groups, the nursery corals had the most consistent major ITS2 type profiles, with 5 of the 6 samples associating with Cladocopium of the C3k radiation (4 unique type profiles). In summary, after 12 months of growth, both the wild and nursery colonies had distinct major ITS2 profiles relative to the initial wild colony major ITS2 type profiles.
Metabolism
Surface area-normalised gross photosynthesis rate (PG, µmol O2 cm-2 h-1) of NT12 (932.5 ± 73.8) was similar to WT0 (932.4 ± 108.8) and ~40% higher than WT12 (690.0 ± 57.5) but overall not statistically different across groups (one-way RM ANOVA, F(1.6,7.9) = 3.4, p = 0.0941; Fig. 4A). Similarly, no differences were found in net photosynthesis (PN, µmol O2 cm-2 h-1) among coral groups (one-way RM ANOVA, F(1.1,5.6) = 1.9, p = 0.2277; Fig. 4B), where PN of NT12 (278.5 ± 51.1) was slightly lower than WT0 and WT12 (442.2 ± 93.8 and 332.4 ± 33.1, respectively). In contrast to PG and PN, respiration rates (R, µmol O2 cm-2 h-1) were higher for NT12 (654.0 ± 35.4) than both WT0 and WT12 (490.2 ± 24.5 and 357.6 ± 50.6, respectively) (one-way RM ANOVA, F(1.5,7.4) = 19.3, p = 0.0017; Tukey’s test, p = 0.0107 for WT0 vs NT12 and, p = 0.0092 for WT12 vs NT12; Fig. 4C). Consequently, the PG:R ratio was lower for NT12 (~1.4) than WT0 (~1.9) or WT12 (~2.1) but not statistically distinguishable between coral groups (one-way RM ANOVA, F(1.5,7.6) = 4.2, p = 0.0674; Fig. 4D). Light-dependent calcification rate (GL, µmol CaCO3 cm-2 h-1) for NT12 (0.3704 ± 0.0811; n =4) was lower than WT0 (0.6726 ± 0.1347) and higher than WT12 (0.1112 ± 0.0411; n = 4) but again not statistically different (Fig. 4E; one-way ANOVA, F(2,11) = 6.6, p = 0.0130; but note Tukey’s test detected differences only for WT0 vs WT12, p = 0.0108). Significant differences were observed for corresponding dark calcification (GD, µmol CaCO3 cm-2 h-1) across groups. In descending order, WT12, NT12, and WT0 were -0.3876 ± 0.1178, 0.1010 ± 0.1043, and 0.2250 ± 0.0654, respectively; Fig. 4F; one-way RM ANOVA, F(1.2,5.9) = 11.9, p = 0.0124), particularly Tukey’s test detected differences between WT12 vs WT0 and NT12, p = 0.0393 and 0.0022, respectively. The GL:GD ratio remained constant for all coral groups (Fig. 4G; Kruskal-Wallis, H = 6.3, p = 0.0327), and Dunn’s test revealed differences between WT0 and WT12 (p = 0.0373).
Energy reserves and tissue biomass
No differences were found for any of the 28 energetic characteristics when nursery-propagating corals were compared to their wild counterparts (NT12 vs WT12) (Fig. 5, Fig. S7 and Fig. S8). However, a time effect was observed for carbohydrate, lipid, and protein content (mg:mg dry weight [DW]) and total energy reserves (J:g DW), whereby wild corals at T12 exhibited more energetic reserves to total tissue biomass compared to WT0 (Tukey’s test WT0 vs WT12; p = 0.0747, 0.0421, 0.0110, and 0.0053, respectively) (Fig. 5A-D). Specifically, carbohydrate concentrations for WT0, WT12, and NT12 were 0.0101 ± 0.0022, 0.0252 ± 0.0023, 0.0293 ± 0.0058, respectively; lipid content was 0.0984 ± 0.0107, 0.1356 ± 0.0098, and 0.1339 ± 0.0071; and protein content 0.0254 ± 0.0034, 0.0407 ± 0.0035, and 0.0415 ± 0.0019, respectively. Energy reserves determined from the carbohydrate, lipid, and protein contents were lower for WT0 (4671.0 ± 468.5) compared to WT12 and NT12 (6770.0 ± 414.9, and 6794 ± 222.3, respectively). Note for comparability, Fig. S6A-D presents the same data normalised to ash-free dry weight and surface area, but where no differences were detected amongst coral groups. Surface area-normalised total dry weight tissue biomass (mg DW cm-2; Fig. 5E) — but not the ash-free total dry weight biomass (mg DW cm-2; Fig. 5F) — was greater for WT0 (21.4 ± 2.3) than WT12 and NT12 (12.1 ± 1.0, and 14.5 ± 1.0). Thus, despite the greater growth of A. cf. hyacinthus for nursery compared to reef colonies over 12 months, all colonies exhibited less tissue biomass (but more energy reserves) at 12 months compared to initial samples. A strong time effect (but not location effect) was also observed when profiling lipids in more detail via 19 fatty acid molecules (Fig. S7). In general, differences in MUFA and PUFA (mono- and poly-unsaturated fatty acids) (p = 0.0324 and 0.0099, respectively) but not in SFA (saturated fatty acids; but note p = 0.0501) were found for both nursery and reef corals at T12 compared to T0 (Fig. S8).
Elemental content
After 12 months propagation, carbon-to-nitrogen ratios (C:N) for WT0, WT12, and NT12 were similar (7.11 ± 0.18, 7.02 ± 0.24, and 6.98 ± 0.54, respectively) (Fig. 6A; ANOVA, F(2,12)= 0.0238, p = 0.9765). In contrast, carbon-to-phosphorus (C:P) and nitrogen-to-phosphorus (N:P) ratios for WT0 (1332.9 ± 223.5, and 189.9 ± 35.7; Fig. 6B) were 100% greater than WT12 (679.6 ± 72.0, and 94.8 ± 7.6; n = 5) and NT12 (604.7 ± 95.9, and 84.4 ± 8.55), respectively (Fig. 6C). This change over time appeared driven by lower P content in WT0 (0.0208 ± 0.0028 mmol P g sample-1; Fig. S9D). Furthermore, C:N:P was highest for WT0 (1333:190:1) compared to W12 and NT12 (680:95:1 and 605:84:1, respectively), and all were significantly higher than the Redfield ratio of 106:16:143.
In general, no differences were detected over time (WT0 vs WT12) or at 12 months post-propagation (WT12 vs NT12) for C, N, P, and 17 major and trace elements (Fig. S9), except for Cu content of NT12 (0.8836 ± 0.1097) which was similar to WT0 (0.9543 ± 0.0624) and higher than WT12 (0.5567 ± 0.0434); one-way ANOVA, F(2,13) = 6.8, p = 0.0096; Tukey’s test WT0 vs WT12, and WT12 vs NT12; p = 0.0159, and 0.0260, respectively; Fig. S9L.
Skeletal properties
Skeletal traits appeared largely unchanged after 12-months of propagation (WT12 vs NT12) (bulk volume, bulk density, biomineral density, pore volume, apparent (internal) porosity, hardness, and colony mass per area; Fig. 7A-G). The only exception was hardness (HD) (Fig. 7F; one-way RM ANOVA, F(1.4,7.2) = 8.2, p = 0.0182, Tukey’s test WT0 vs WT12, p = 0.0488) where values were higher for WT0 (38.8 ± 1.5) than both WT12 or NT12 after 12 months (31.5 ± 1.1, and 35.3 ± 0.9, respectively). Thus, overall, propagation in nurseries compared to the reef did not result in a difference in skeletal properties despite the greater growth rates observed for nursery colonies (Fig. 2B).
Multi-trait assessment and “trait redundancy”
We finally performed a series of principal components analysis (PCA) to reduce the trait space and identify the main drivers of variance from amongst the 90 traits examined (Table S2). Specifically, a PCA was first conducted for each of the following seven biological properties: physical appearance, photobiology, metabolism, energy reserves, FAME, elementome, and skeletal traits; to identify the two main contributing traits to the variation (Table S3A) beyond survivorship and growth, and in doing so where “trait redundancy” (i.e. traits that have similar contribution to the separation of coral groups) existed within biological properties. This is observed via overlapped trait-based vectors and correlated variables in Fig. S10 and Table S3A, respectively. A total of 14 (of 90) traits (brightness, blue, symbiont cell density, maximum photochemical efficiency of PSII [Fq´/Fm´MAX], dark and light calcifications [GD and GL, respectively], total energy reserves, ash-free dry weight biomass, methyl arachidate, methyl myristoleate, potassium [K], strontium [Sr], pore volume, and biomineral density) were identified as the strongest drivers of variation from the individual PCAs (Table S4), and were therefore subsequently selected for the final PCA (the so-called multitraits, Fig. 8). We identified three clusters corresponding to the three groups (WT0, WT12, and NT12), despite some overlap; and where WT0 was somewhat separated from WT12 and NT12 (Fig. 8). Overall, 51.6% of the coral trait variance between coral groups was accounted by the first and second principal components collectively (PC1 and PC2, respectively; Table S4B). PC1 accounted for 34.5% of the total trait variance of the coral holobiont, with the following top 4 loadings: methyl arachidate (15.41%), total energy reserves (14.49%), blue (11.47%), and K (10.47%) (see Table S4A for contributions by the rest of the traits) contributing the largest loadings to this vector (Fig. 8). A further 17.1% of the total trait variance was accounted for by PC2, with brightness (25.77%), dark calcification (15.80%), ash-free dry weight biomass (11.51%), and symbiont cell density (8.23%) followed by other traits (see Table S4A) contributing the largest loadings to this vector (Fig. 8). Differences between coral groups were dependent on both PC1 and PC2, with PC resolving separation between all three coral groups. Here, WT0 separated from NT12 by PC1 (p = 0.0237; Table S5A), whilst all three coral groups were separated from each other by PC2 (WT0 - WT12, WT0 - NT12, and WT12 – NT12; p = 0.025, 0.027, and 0.0067, respectively). However, the one-way PERMANOVA between coral groups for the main 14 traits did not return differences between coral groups (p = 0.3706, Table S6).