{'fraction_non_alphanumeric': 0.0389287173581962, 'fraction_numerical': 0.02364325247947699, 'mean_word_length': 5.183875043068795, 'pattern_counts': {'":': 0, '<': 0, '<?xml version=': 0, '>': 0, 'https://': 2, 'lorem ipsum': 0, 'www.': 2, 'xml': 0}, 'pii_count': 0, 'substrings_counts': 0, 'word_list_counts': {'cursed_substrings.json': 1, 'profanity_word_list.json': 0, 'sexual_word_list.json': 0, 'zh_pornsignals.json': 0}}
Microbiome Restructuring: Dominant Coral Bacterium Endozoicomonas Species Respond Differentially to Environmental Changes

ABSTRACT Bacteria in the coral microbiome play a crucial role in determining coral health and fitness, and the coral host often restructures its microbiome composition in response to external factors. An important but often neglected factor determining this microbiome restructuring is the ability of microbiome members to respond to changes in the environment. To address this issue, we examined how the microbiome structure of Acropora muricata corals changed over 9 months following a reciprocal transplant experiment. Using a combination of metabarcoding, genomics, and comparative genomics approaches, we found that coral colonies separated by a small distance harbored different dominant Endozoicomonas-related phylotypes belonging to two different species, including a novel species, “Candidatus Endozoicomonas penghunesis” 4G, whose chromosome-level (complete) genome was also sequenced in this study. Furthermore, the two dominant Endozoicomonas species had different potentials to scavenge reactive oxygen species, suggesting potential differences in responding to the environment. Differential capabilities of dominant members of the microbiome to respond to environmental change can (i) provide distinct advantages or disadvantages to coral hosts when subjected to changing environmental conditions and (ii) have positive or negative implications for future reefs. IMPORTANCE The coral microbiome has been known to play a crucial role in host health. In recent years, we have known that the coral microbiome changes in response to external stressors and that coral hosts structure their microbiome in a host-specific manner. However, an important internal factor, the ability of microbiome members to respond to change, has been often neglected. In this study, we combine metabarcoding, culturing, and genomics to delineate the differential ability of two dominant Endozoicomonas species, including a novel “Ca. Endozoicomonas penghunesis” 4G, to respond to change in the environment following a reciprocal transplant experiment.

Inner Bay (here Inner Bay) of Penghu has suffered substantial marine biodiversity losses, including significant damage to marine aquaculture, wild fisheries, and coral bleaching due to extreme weather events in the winter (39). Furthermore, based on regional news and government reports, domestic sewage dumping, the presence of a shipping port, and aquaculture practices have increased the concentrations of nitrogen and ammonia in the calmer waters of the Inner Bay compared to the Outer Bay region, threatening corals. The contrasting local environmental conditions between the Inner and Outer Bays make them excellent sites to study the response of locally acclimated dominant coral microbial communities, e.g., Endozoicomonas, as members of this genus show lower relative abundance in degraded or anthropogenically impacted reefs compared to the one that is healthy (28). Further, contrasting conditions in the sites also help trace coral microbiome restructuring at a fine scale and test the differential capability of dominant coral-associated bacteria to respond to environmental change.
We examined how the microbial community restructures after reciprocal transplant and whether the Endozoicomonas phylotypes have differential responses to change in local conditions using colonies of the coral Acropora muricata (genus Acropora) in the Penghu Archipelago, Taiwan. These coral species in the Penghu Archipelago have been reported to harbor Endozoicomonas as their dominant bacteria (24). We conducted a longitudinal (9-month) in situ reciprocal transplant experiment with repeated sampling, where coral colonies from the semiclosed Inner Bay were transplanted into the open ocean region of the Outer Bay and vice versa. Furthermore, we aimed to isolate, culture, and characterize dominant Endozoicomonas members to provide genomic insights into how bacteria adapt to these environments.

RESULTS
Sampling and sequencing overview. We collected a total of 110 coral and 12 seawater samples from the experiment, of which 10 coral fragments were removed (all from the Inner Bay) during the experiment (marked with an X in Fig. 2A and B), as they appeared to be dead. These dead samples were only used to help contrast with the microbial community compositions of healthy corals and were later removed before downstream analyses, including a and b diversity analyses. At the end of the experiment, we had 100 coral and 12 seawater samples. A total of 2,015,935 high-quality reads (an average of 16,524 reads per sample) were obtained after removing chimeras and poor-quality reads from the 110 coral and 12 seawater samples. These reads were denoised into 2,064 zero-radius operational taxonomic units (zOTUs). Healthy corals (n = 100) had 1,815,002 reads (range, 6,092 to 67,598) and 1,966 zOTUs, and seawater (n = 12) had 149,638 reads (range, 6,897 to 19,777) and 1,521 zOTUs.
Coral and seawater microbiomes differ in bacterial diversity and compositions. The coral and seawater samples were significantly different in bacterial diversity and evenness, measured through zOTU richness (see Fig. S1A in the supplemental material), Shannon (Fig. S1B) and Chao1 (Fig. S1C) diversities, and inverse Simpson evenness (Fig. S1D). The seawater samples had more than twice the zOTU richness and Shannon and Chao1 diversities of the coral samples. August samples showed an unusual alpha diversity pattern, particularly in seawater samples; this could be because the samples were collected during heavy rainfall (63 days). We observed an increase in richness and Chao1 between O!I and OC samples, but no apparent differences were observed between I!O and IC samples. In terms of corals at different locations, there was no significant difference in calculated alpha diversity measures between control and transplant samples from the Inner and Outer Bays (Fig. S1A to D).
Proteobacteria (specifically class Alphaproteobacteria) (average relative abundance, 51.71%) and Bacteroidetes (22.76%) (class Flavobacteriia) were the dominant phyla in seawater samples across all time points, followed by Cyanobacteria (17.75%), which was particularly abundant in Inner Bay samples. There was heavy rainfall during the week in August that samples were collected, and we noticed a higher abundance of Marinimicrobia (class Marinimicrobia SAR406 clade), Planctomycetes, and Verrucomicrobia in those samples, which might explain the different bacterial diversity and evenness results obtained that month. On the contrary, coral samples were dominated by Proteobacteria (specifically class Gammaproteobacteria) (73.08%), Chlamydiae (13.70%), and Tenericutes (8.99%) (class Mollicutes). Dead coral samples had a bacterial community composition similar to that of seawater samples ( Fig. 2A and Fig. S2). Changes in the coral microbial community throughout the reciprocal transplant. Proteobacteria was the dominant phylum across all sample groups (control, IC and OC; transplant, I!O and O!I) from the two locations throughout the experiment. We observed shifts in the microbial community of control samples (IC and OC) from April to August and December. The relative abundance of Chlamydiae (family Simkaniaceae), the second abundant phylum in the Outer Bay (OC), with all zOTUs (10 in count) belonging to "Unclassified Simakaniaceae," decreased from April to August before increasing slightly in December. Tenericutes (specifically Mollicutes), the second dominant phylum in the Inner Bay control samples (IC), increased in abundance over time, peaking in December. One zOTU annotated as "Unclassified Entomoplasmatales" had the highest abundance among different members of Tenericutes, including zOTUs belonging to Acholeplasma, Mycoplasma, "Candidatus Bacilloplasma," and "Candidatus Hepatoplasma." In July, we observed a sudden spike in Verrucomicrobia abundance in two samples from IC, but soon after in August the community composition became similar to that in June. We also observed patterns of community dynamics in Chlamydiae and Tenericutes in cross-transplant samples (I!O and O!I) over the sampling period. Chlamydiae became the most dominant group in O!I (May) samples, but its abundance decreased sharply thereafter, whereas in I!O samples, Chlamydiae and Tenericutes both remained stable, with Chlamydiae being dominant in May and June and Tenericutes being dominant in July and August ( Fig. 2A and Fig. S2).
At the genus taxonomic rank, Endozoicomonas species were the most dominant. Forty-eight zOTUs (out of 2,064) were taxonomically classified as Endozoicomonas. These 48 zOTUs accounted for an average of ;54% relative abundance in coral fragments and 0% in sea water samples. Of these 48 zOTUs, 13 contributed ;90% of the total Endozoicomonas abundance. Interestingly, IC and OC samples harbored different dominant Endozoicomonas zOTUs, with zOTU1 and zOTU2 being dominant in OC and zOTU7 and zOTU9 in IC. Across the sampling time, we also observed shifts in the dominant Endozoicomonas phylotypes. zOTU2 was dominant from April to June, whereas zOTU1 became dominant in July-December OC samples. In IC samples, zOTU7 was dominant from April to May, but after that its relative abundance declined (Fig. 2B). It is also worth noting that a significant decline in the Endozoicomonas abundance was observed in IC samples in July, August, and December (Fig. 2B), suggestive of location dependence.
In cross-transplant samples, the Endozoicomonas phylotypes from OC remained resistant to change when transplanted in the Inner Bay (O!I), with zOTU2 being dominant across all sampling times. For I!O transplanted samples, however, instead of Endozoicomonas phylotypes from IC, we observed that zOTU1, another dominant Endozoicomonas phylotype in OC samples, was dominant (Fig. 2B), suggesting that the phylotypes had different robustness under different environmental scale disturbances.
Location-dependent robustness in the coral microbiome. The dispersion of homogeneity analysis identified that the bacterial community in corals with the Inner Bay as the final location (IC and O!I) were significantly different from each other (analysis of variance [ANOVA], F = 9.23, P , 0.001), whereas samples whose final destination was the Outer Bay (OC and I!O) had no significant difference (ANOVA, F = 1.98, P . 0.05), indicating that the microbiome had location specificity. Therefore, samples whose final destinations were the Inner Bay and Outer Bay were analyzed independently to test for differences in community composition between the control and transplant groups. Ordination analysis using nonmultidimensional scaling (nMDS), followed by permutational multivariate analysis of variance (PERMANOVA), identified the significant influence of coral sample, sampling month, and their combined effect (interaction term) (Fig. 3A). Ellipses with a 95% confidence interval suggested that samples for which the Outer Bay was their final location (OC and I!O) were more similar to each other than to samples for which the Inner Bay was their final location (IC and O!I). These findings support locational variability and differential robustness in the coral microbiome ( Fig. 3A and B). Transplanted samples (O!I) clustered tightly compared to IC samples, indicating less variability after transplantation in the transplant samples.
However, highly overlapping ellipses were observed for OC and I!O samples, suggesting a highly similar microbial community in the control (OC) and transplanted samples from the Inner Bay (I!O) ( Fig. 3A and B).
Differentially abundant microbiome dominated by Endozoicomonas-related phylotypes. Linear discriminant analysis effect size (LEfSe) analysis identified differentially abundant zOTUs of different taxa across all the sampling groups: nine zOTUs were differentially abundant in OC samples, 13 in I!O samples, and 16 in IC and O!I samples ( Fig. 3C and D). Interestingly, all the differentially abundant zOTUs in the OC samples belonged to Endozoicomonas, but zOTUs belonging to diverse taxa, including the BD1-7 clade (Gammaproteobacteria), Entomoplasmatales (phylum Tenericutes; class Mollicutes), and Alteromonadaceae (class Gammaproteobacteria), were differentially abundant in I!O samples (Fig. 3C). Similarly, out of the 16 zOTUs that were differentially abundant in O!I samples, 13 were Endozoicomonas; IC samples also had zOTUs belonging to diverse taxa that were differentially abundant, including Surface 1_ge (class Alphaproteobacteria), Synechococcus (class Cyanobacteria), and others (Fig. 3D).
Description of "Ca. Endozoicomonas penghunesis" 4G. "Ca. Endozoicomonas penghunesis" 4G is a Gram-negative, facultatively anaerobic, and slightly motile bacterium that forms beige-colored colonies (size, 2.14 by 0.66 mm) and is slightly susceptible to the antibiotics streptomycin and ampicillin. No catalase enzymatic activity was reported for this bacterium, but the bacterial culture was trypsin and oxidase positive (Table S2). This new bacterial species was identified to tolerate the specified temperature (15 to 35°C) and salinity (5 to 30 PSU) ranges (Table S3). The pH range for growth was pH 6.0 to 10.0, with optimal growth observed at a slightly alkaline pH (pH 8.0). Scanning electron microscopy (SEM) and transmission electron microscopy (TEM) analyses identified rod-shaped cells (Fig. S4C) surrounded by a possible mucus lining (Fig. S4A) and with structures that appeared to be granules or vacuoles in the cell (Fig. S4B). A 16S rRNA gene sequence blast search identified the closest cultured relative to be Endozoicomonas montiporae CL-33 (GenBank accession no. CP013251) with 96.17% identity; based on a species identity cutoff of 97%, this suggests that the bacterium is a novel Endozoicomonas species. Further, average nucleotide identity (ANI) analysis on genome sequences confirmed our assertion that "Ca. Endozoicomonas penghunesis" 4G is a novel species (Fig. S3E).
Genome assembly features of "Ca. Endozoicomonas penghunesis" 4G. The genome of "Ca. Endozoicomonas penghunesis" 4G was first assembled using Nanopore reads and later polished with quality filtered Illumina reads, resulting in a single contig of 6,004,453 bp and N 50 of 6,004,453. Genome completeness, contamination, and strain heterogeneity were estimated to be 97.52, 0.98, and 0%, respectively. Out of 573 single-copy marker genes (c_Gammaproteobacteria) from the checkM database (40), 493 genes were present only once, five single-copy markers were duplicated, and nine were missing. Based on the criteria of "minimum information for single amplified and metagenome-assembled genome of bacteria" (41), our genome can be considered "finished." The GC content of the genome was 49.1%, which is similar to that of other Endozoicomonas species. Genomic features of "Ca. Endozoicomonas penghunesis" 4G. A total of 5,019 genes and 4,913 coding DNA sequences (CDS) were predicted from the genome. We annotated seven copies of 16S, nine of 23S, eight of 5S rRNA genes, and 80 tRNAs. A sequence similarity analysis of 16S rRNA gene copies revealed that all copies were at least 98.76% identical, with four copies .99.28% identical and two 100% identical to each other (Fig. S3C). Copy-1 of the 16S rRNA gene was used as a representative sequence to classify the closest relative of dominant Endozoicomonas zOTUs identified in this study (Fig. 4A). There were no CRISPR elements, and only one prophage was identified in the genome. Out of the 5,019 genes predicted, more than 50% (2,721) were annotated to be hypothetical. Since Endozoicomonas species have been exclusively isolated from their marine eukaryotic hosts, including "Ca. Endozoicomonas penghunesis" 4G, we searched for eukaryote-like proteins (ELPs) and identified 43 WD40 domain proteins (WD40), four ankyrin repeat proteins (ARPs), and 12 tetratricopeptide repeat proteins (TRPs). Almost all the WD40 domain-containing proteins were arranged consecutively ( Fig. 4B) and flanked by transposons. Most (27 out of 41) of the WD40 domain proteins were annotated as TolB protein from the Tol-Pal system, which is important for maintaining cellular integrity, and others (14) were classified as hypothetical proteins. A wide array of secretory proteins was also annotated with 248 type III secretion system (T3SS), 50 type IV secretion system (T4SS), and 10 type VI secretion system (T6SS) effectors annotated from the proteome.
Metabolic repertoire of "Ca. Endozoicomonas penghunesis" 4G. RAST classified only 36% (1,756) of the total genes into subsystems. Subsystems (i) carbohydrates and (ii) cofactors, vitamins, prosthetic groups, and pigments had the highest number of annotated genes, 270 and 238, respectively (Fig. S5). In the stress response subsystem, 108 genes were annotated, most of which were related to oxidative stress response (46 genes), followed by heat shock response (18) and detoxification response (16). Interestingly, within the osmotic stress response, we identified genes for betaine transport via ATP-binding cassette transporter BetS (high-affinity choline uptake protein BetS), arranged in tandem with an L-proline glycine betaine ABC transport system permease (ProV and OusW) (Fig. 4B). Multiple copies of superoxide dismutase, alkyl hydroperoxide reductase, and methionine sulfoxide reductase genes were also identified in the genome.
"Ca. Endozoicomonas penghunesis" 4G had genes encoding essential amino acids and pathways including glycolysis and tricarboxylic acid cycle; genes for converting nitrate to nitrite (NapAB, KEGG Ontology no. K02567) and ammonia to L-glutamate were identified, but none were related to the conversion of nitrite to ammonia. Assimilatory and dissimilatory sulfate reduction and oxidation pathways were also completely absent. Furthermore, no genes related to the uptake of extracellular taurine or its metabolism to sulfite were identified. Interestingly, siroheme biosynthesis and sirohemedependent anaerobic sulfite reduction operons were present in "Ca. Endozoicomonas penghunesis" 4G (Fig. 4B). We also identified genes arranged in an operon-like manner for anaerobic glycerol degradation. Comparing E. acroporae and "Ca. Endozoicomonas penghunesis" physiological and genomic features. Endozoicomonas phylotypes belonging to E. acroporae and "Ca. Endozoicomonas penghunesis" 4G were dominant in colonies of coral Acropora muricata in the Outer and Inner Bay, respectively ( Fig. 2B and 4A). This selective dominance could be attributed to multiple factors, including bacterial physiological and genetic repertoire. Therefore, we compared the two species; "Ca. Endozoicomonas penghunesis" had a wider growth temperature range than E. acroporae and was slightly motile. E. acroporae, on the other hand, was nonmotile (Table S2) and had a wider salinity and growth pH range (Table S3). A wider growth temperature range of "Ca. Endozoicomonas penghunesis" and its dominance in the Inner Bay aligns with the high variation of temperature fluctuations (from summer to winter) observed in the calmer waters of the Inner Bay. Comparing the metabolic repertoire of the two species, we found that genes for dimethylsulfoniopropionate (DMSP) metabolism and dimethyl sulfoxide (DMSO) reduction were absent from "Ca. Endozoicomonas penghunesis," but E. acroporae had a complete operon for DMSP metabolism, as reported in our previous study (36), as well as genes for DMSO reduction. The lack of a potent oxidative stress response gene repertoire could be the reason for the loss of abundance of "Ca. Endozoicomonas penghunesis" in the control (IC) and transplant samples (I!O) in summer (June-August). Furthermore, the robustness of E. acroporae in the control (OC) and transplant samples (O!I) throughout the year could be due to the species' ability to remove oxidative stress (which increased in summer) more efficiently via DMSP metabolism and the presence of catalase activity. This was further confirmed by using the 2,2 diphenyl-1-picrylhydrazyl (DPPH) assay as a proxy for measuring the ROS scavenging ability, where E. acroporae was identified as a strong free radical scavenger compared to "Ca. Endozoicomonas penghunesis" in ambient (25°C) as well as higher temperature conditions (30°C) (Fig. S6).
Coral mortality in the Inner Bay. We observed coral mortality exclusively in the Inner Bay for both IC and O!I samples ( Fig. 2A and B, marked X). Grazing by Drupella cornus was also observed in the samples from the Inner Bay during the sampling in May and continued until December (Fig. S7).

DISCUSSION
This study aimed to test the differential capability of Endozoicomonas, one of the most dominant bacterial groups in the coral microbiome toward environmental change. We analyzed the microbiome dynamics of the common Indo-Pacific coral A. muricata over 9 months following a reciprocal transplant experiment at the finest resolution of zOTUs. We identified that different Endozoicomonas phylotypes in the A. muricata coral colonies belong to two dominant species, including a novel species, and have differential adaptation capabilities, with one species more resilient to environmental change than the other. Our results shed light on an often-neglected factor when determining variations in community composition: bacterial species/strains adapt differently when coral hosts are subjected to biotic and abiotic stressors. Furthermore, we also isolated, cultured, and sequenced a single chromosome-level genome of one of the dominant phylotypes belonging to a novel Endozoicomonas species, "Ca. Endozoicomonas penghunesis" 4G, to ascertain the ecological and functional role of this bacterium and add to the growing knowledge and genome data sets of this key microbe in coral reefs.
Acropora muricata microbiome was dominated by members of class Gammaproteobacteria (phylum Proteobacteria) (Silva v132 grouped Gammaproteobacteria into Betaproteobacteria), particularly by Endozoicomonas-related phylotypes ( Fig. 2A and B and Fig. S2). Members of the genus Endozoicomonas are often found to be the dominant group in the Response of Endozoicomonas to Environmental Change mSystems microbiome of several coral species, e.g., members of Acropora, Pocillopora, and Stylophora (3,25,31,42), and have been proposed to play a significant role in coral health and protection (2,43) and coral sulfur cycling (36,44). Another dominant bacterial group, Simkaniaceae (phylum Chlamydiae, class Chlamydiae), was described as an obligate intracellular bacterium, but its function has remained enigmatic (42,45). Like Endozoicomonas, Simkaniaceae-related phylotypes were also recently found to be abundant in healthy corals from the reefs in Florida, but their abundance decreased in corals suffering from stony coral tissue loss disease (46). Members of class Mollicutes, particularly zOTUs related to Entomoplasmatales and Mycoplasmatales, are suggested to be mutualistic or commensal bacteria in temperate and deep-sea gorgonians and cold-water Scleractinia corals, but their specific function remains unknown (47,48). Overall, the microbial community composition in the coral colonies of the control group remained stable throughout the experiment timeline, with only transient differences observed between them across sampling times ( Fig. 2A and B). Spatial and temporal fluctuations in the microbiome were observed throughout the experiment, and community structure ( Fig. 2A) and ordination analysis showed that community varies across temporal and spatial scales (Fig. 3A). Various degrees of overlap between the samples from the same location suggested that the microbes show different scales of variability and that microbial structure is a function of the local environment. Previous studies have also shown that the microbiome varies spatially due to differences in the sites' local environments (19,49). Ocean currents are believed to have a homogenizing effect on the microbial communities, and the same coral species separated by hundreds to thousands of kilometers have been found to have similar microbiome compositions (1, 50). Our results were in contrast to this observation, as a high site-to-site variation was observed at a relatively small scale. One potential reason for this variation is high abiotic and anthropogenic pollution in the Inner Bay compared to the Outer Bay; similar results were obtained in earlier studies from the Red Sea (28) and Singapore (51).
Utilizing the zOTU approach for metabarcoding data analysis and focusing on the dominant bacterial genus, we identified that colonies of A. muricata from the Inner and Outer Bays were not dominated by a single Endozoicomonas phylotype but had several differentially abundant phylotypes associated with them (Fig. 2B). These results are similar to the multiple phylotypes that were reported to be dominant in colonies of the same coral genus as identified earlier (42). However, the dominant Endozoicomonas phylotypes identified in our study were different in coral colonies in the two locations, with Inner Bay colonies harboring a novel species "Ca. Endozoicomonas penghunesis" and Outer Bay colonies harboring E. acroporae-related phylotypes (Fig. 4A). This result was intriguing, as corals of the Acropora genus are known to have a strong influence on their microbial community composition (19). Another observation that arises from these results is how the single-nucleotide variation approach utilized to obtain ASVs or zOTUs potentially can lead to increased diversity (richness) estimates, which rely on ASV or zOTU counts, especially with bacterial groups known to have more than one copy of nonidentical 16S rRNA genes in their genome. In our study this is true in the case of the genus Endozoicomonas, as members of this genus are known to harbor more than one copy of 16S rRNA and complete genomes of Endozoicomonas montiporae CL-33 T have seven copies (52), similar to "Ca. Endozoicomonas penghunesis," which are not all identical (Fig. S3); hence, we used a phylogenetic approach to assign taxonomy to these Endozoicomonas phylotypes.
Members of the coral holobiont potentially engage in complex interactions to maintain the health and fitness of the coral host, and external stressors could disturb these interactions by influencing the composition of the holobiont. To overcome the influence of the external stressors, a genomic adaptation of the holobiont members (in our case, bacteria) could play an important role in host survival, since a bacterium is likely to be not as well adapted to a new niche as resident strains, unless it has the genetic capability to mitigate the new stressors (53,54). During our reciprocal transplant experiment, we observed that E. acroporae-related zOTUs remained dominant in the control (OC) and the transplanted samples (O!I), with only a change in the dominant zOTU from zOTU1 to zOTU2, whereas those related to "Ca. Endozoicomonas penghunesis" 4G were only dominant in April and May months in the control samples of the Inner Bay (IC) (Fig. 2B). There are a few possible reasons for this observation. One is that E. acroporae may be more resilient and better adapted to diverse conditions like increased oxidative stress during summers, as evident from the DPPH assay compared to "Ca. Endozoicomonas penghunesis" 4G. However, further investigation is required to say this with greater confidence, as the DPPH assay only measures polyphenolic compounds and, hence, cannot be taken as representative of the overall bacterial antioxidant capacity. Another potential explanation is the restructuring of the microbiome by the host A. muricata. The latter seems more plausible with evidence from recent studies on Acropora.
In addition, variation in the abundance of the different Endozoicomonas phylotypes is potentially analogous to the abundance of genotypes of different endosymbiotic algae Symbiodinium, often found in coral colonies (55). Several coral species can perform "symbiont shuffling or switching" to select for the more thermotolerant genotype of endosymbiotic algae in response to thermal stress (56)(57)(58). We observed microbial shuffling for Endozoicomonas phylotypes to a certain degree in our transplant samples. zOTU1 (potentially E. acroporae (Fig. 4A)) became the dominant Endozoicomonas phylotype in I!O samples, although zOTU7 and zOTU9 (potentially "Ca. Endozoicomonas penghunesis" 4G [ Fig. 4A]) were the dominant phylotypes in IC samples (Fig. 2B). It is important to note here that no microbial shuffling was observed in the case of O!I versus OC samples, as phylogenetic analysis identified both zOTU1 and -2 as E. acroporae (Fig. 4A). Furthermore, the DPPH assay identified a higher potential for E. acroporae to scavenge ROS than "Ca. Endozoicomonas penghunesis" (Fig. S6), so shuffling microbial community might confer selective advantages to the coral host but further exploration is required.
Understanding the functional and ecological roles of coral-associated microbes in reefs has become critical to developing an intervention for coral reef protection, such as developing a coral probiotic (2,59). These interventions require in-depth information about members of the coral holobiont. In the current study, we isolated, cultured, and sequenced the complete genome of a dominant Endozoicomonas phylotype identified in the metabarcoding data analysis. Phylogenetic analysis identified that the dominant zOTUs (zOTU1, zOTU2, zOTU10, zOTU11, zOTU13, and zOTU18) from the Outer Bay were closer to a previously characterized species, E. acroporae (60), whose genome was sequenced earlier (36,61). On the other hand, the Inner Bay-dominant zOTUs (zOTU7, zOTU9, zOTU16, and zOTU17) were closest to the novel species "Ca. Endozoicomonas penghunesis" 4G isolated and characterized in this study. Genomic analysis of "Ca. Endozoicomonas penghunesis" revealed features similar to those of other Endozoicomonas species, i.e., the large genome size (;6.00 Mb), many coding genes (4,913), and complete pathways for essential amino acids, suggesting a free-living life stage. Furthermore, a comparative genomics study identified that Endozoicomonas species are capable of differential functional specificity, and different genotypes may play disparate metabolic roles in their hosts (37). This is true for sulfur metabolism, where E. acroporae is the only known Endozoicomonas species capable of metabolizing dimethylsulfoniopropionate (DMSP) to dimethylsulfide (DMS) (36). Genes for DMSP metabolism operon were not present in the genome of "Ca. Endozoicomonas penghunesis." However, transporters (three copies) for glycine-betaine, another osmolyte, were identified in the genome of "Ca. Endozoicomonas penghunesis." Other Endozoicomonas species have also been identified to have the ability to scavenge glycine-betaine through transporters (62), potentially alleviating oxidative stress. Identification of putative siroheme-dependent anaerobic sulfite reduction operon was interesting as this process facilitates growth under anaerobic conditions (B 12 -dependent anaerobic growth) by oxidizing 1,2-propanediol with tetrathionate as an electron acceptor (63). Physiological tests also showed that "Ca. Endozoicomonas penghunesis" is a facultative anaerobe; however, more functional evidence is required to confirm this outcome and the advantage (if any) it provides to the bacterium and coral host that maintains it.
The antioxidant-producing ability of marine invertebrate-associated bacteria had been wildly reported (64,65). Living in highly variable environments, marine bacteria have developed different strategies to protect themselves from reduced oxygen intermediates and increase environmental adaptability (66). Compared to "Ca. Endozoicomonas penghunesis" 4G, E. acroporae Acr-14 T showed stronger free radical scavenging ability at higher temperatures, no matter in cell lysate or supernatant (see Fig. S6 in the supplemental material). Based on our results, we suggest that the dominance of E. acroporae-related phylotypes after reciprocal is due to higher antioxidant abilities, especially in summer. However, the interactions between coral and bacteria are complicated, and it is hard to clarify the specific function of bacterial antioxidants and their relationship with niche selection. Therefore, more focused studies are required to discern the putative benefits of the coral microbiome concerning oxidative stress management for the host.
We observed coral mortality during our experiment exclusively in the Inner Bay; the corallivorous snail Drupella cornus, which exclusively feeds on living tissue, grazed there (Fig. S7). These gastropods occur throughout the shallow waters of the Indo-Pacific region (67). Outbreaks of this corallivorous marine gastropod have been recorded in different parts of the Gulf of Eilat, Israel (68), and the Great Barrier Reef, Australia (69). Coral feeding gastropods of Drupella species show a strong preference for preying on Acroporids (70) and are known to be efficient vectors for brown band disease in corals (69,71). Although no Drupella species outbreaks to date have been recorded in Taiwan's coral reefs and no visible signs of brown band disease were observed in our study, it is important to keep monitoring the corals in the Penghu Archipelago for signs of climate change and disease outbreaks in the near future.
Conclusions. A variety of factors, many of which are external, are known to influence the coral microbiome composition and its dynamics. However, an important internal factor, the adaptation capability of microbiome members, which governs the survival of a bacterium in a niche, has been overlooked. Using a combination of metabarcoding, genomic, and comparative genomic approaches, we showed that members of the dominant bacterial group Endozoicomonas are capable of sustaining and proliferating in a new niche following a reciprocal transplant experiment. Our ability to isolate and culture one of the dominant bacterial species, "Ca. Endozoicomonas penghunesis" 4G, builds on our knowledge of these important bacterial groups in the coral holobiont. Furthermore, we address critical aspects of using zOTUs/ASVs to estimate bacterial richness using metabarcoding data, which can result in often falsely inflated diversity estimates, especially in the case of microbes harboring more than one copy of nonidentical 16S rRNA gene, e.g., Endozoicomonas. In summary, we conclude that different members of the coral holobiont belonging to the same bacterial group can have differential adaptation capabilities, and this internal factor should also be considered when devising interventions to protect coral reefs, like developing a coral probiotic.

MATERIALS AND METHODS
Study design and experimental setup. Five colonies of Acropora muricata (40 by 40 cm) were collected at a depth of 3 m from the Outer Bay (O) (N23°33.097' E119°38.3359) and Inner Bay (I) (N23°3 1.8539 E119°33.629') along the reef adjacent to the coast of the Penghu Archipelago, Taiwan (Fig. 1A). Mother coral colonies were first collected in April. These acted as controls for the native coral microbiome in the study sites (Fig. 1A). Later, mother colonies were fragmented into two halves (approximately 20 by 20 cm each). Coral fragments from mother colonies were either cross swapped (I!O or O!I) or transplanted in their original location (IC or OC) (Fig. 1B). Coral fragments from each colony that remained in their original location acted as controls to measure any change in the microbiome due to the transplant procedure and change in the microbiome based on colony age and experimental time. Coral fragments were glued onto the reef with epoxy putty.
Sampling timeline and sample collection. Study sites were visited every month from April to August 2018 and then in December 2018 to check the status of transplanted fragments and collect samples. In total, we collected 122 samples, including seawater samples (1 liter) at each time point and location. Two-by 2-cm fragments were taken from each colony, rinsed with filtered seawater, and stored in 99% ethanol at 220°C before DNA extraction.
DNA extraction and 16S rRNA gene amplicon sequencing. Frozen coral fragments were sprayed (70 lb/in 2 ) with ;15 ml 1Â TE buffer (10 mM Tris-HCl, 1 mM EDTA, pH 8) and then placed into sterile Ziploc bags. Total DNA was extracted using a modified cetyltrimethylammonium bromide (CTAB) Response of Endozoicomonas to Environmental Change mSystems method (72). Coral tissue sample DNA was extracted with conventional chloroform-isoamyl alcohol (24:1) and phenol-chloroform-isoamyl alcohol (25:24:1) step and isopropanol precipitation method. The DNA pellet was rinsed with 70% ethanol and then dissolved in 50 ml double-distilled water and stored at 220°C. DNA concentration was determined using a NanoDrop 1000 spectrophotometer (Thermo Fisher Scientific, Waltham, MA) and Quant-iT dsDNA HS (high-sensitivity) assay kit. Seawater samples were processed similarly with the modified CTAB method (72). For DNA library construction, 968F (59-AAC GCG AAG AAC CTT AC-39) (73) and 1391R (59-ACG GGC GGT GWG TRC-39) (74) universal primers were used to amplify the bacterial V6-V8 hypervariable region of the 16S rRNA gene from the total DNA from samples using PCR. For PCRs, 50 to 150 ng of template DNA was used. PCR was performed in 50-ml reaction volumes, consisting of 1.5 U TaKaRa Ex Taq (TaKaRa Bio, Otsu, Japan), 1Â TaKaRa Ex Taq buffer, 0.2 mM deoxynucleotide triphosphate mixture (dNTPs), 0.2 mM forward and reverse primers, and template DNA. The PCR conditions consisted of an initial denaturing step at 95°C for 5 min, followed by 30 cycles at 94°C for 30 s, 52°C for 20 s, and 72°C for 45 s and a final extension at 72°C for 10 min. The amplified product was visually confirmed using 1.5% agarose gel with a 5-ml PCR product. Target bands (;420 bp) were cut and eluted using a QIAEX II gel extraction kit (Qiagen, Valencia, CA, USA).
Each bacterial V6-V8 amplicon was tagged with a unique barcode sequence by designing tag primers with 4-base overhangs at 59 ends. The sample-specific tagging reaction was performed with a 5cycle PCR, with a reaction program of initial denaturation at 94°C for 3 min, followed by denaturation at 94°C for 30 s, annealing at 52°C for 20 s, extension at 72°C for 45 s, and a final extension at 72°C for 10 min. The amplified product was purified using the QIAquick PCR purification kit (Qiagen, Valencia, CA). Purified products were pooled into four independent libraries and sequenced with Illumina MiSeq paired-end sequencing (2 Â 300 bp) at Yourgene Biosciences, Taiwan. No kit and PCR negative controls were sequenced in this study.
Sequence data processing and analysis. Paired-end raw reads obtained from Illumina sequencing were merged using USEARCH v11 (75) with the parameters minovlen = 16, maxdiffs = 30, and pctid = 80. Merged reads were sorted, quality filtered, and trimmed using Mothur v1.3.81 (76). Reads 400 to 470 bp long with an average quality of .25 were kept. Chimeric reads were inspected and eliminated with UCHIME (77) by USEARCH v11. Qualified reads were retained for subsequent analysis. High-quality reads were denoised using UNOISE3 (78), and zero-radius operational taxonomic units (zOTUs), which are equivalent to exact sequence variants, were obtained. The denoised sequences were aligned against the SILVA128 (79, 80) rRNA database for a taxonomic assignment up to the genus level using Mothur on a per-sample basis with a pseudobootstrap cutoff of 80%.
Statistical analyses. All statistical analyses and graphs were generated in R (R Core Team, 2020). Stacked bar plots were obtained by converting absolute abundance profiles into relative abundances. Abundance profiles were processed with the R packages phyloseq (81), vegan (82), ggplot2 (83), pheatmap (84), and microbiomeMarker (85) for downstream analyses and visualization. Alpha diversity analysis was conducted after rarifying the samples to an even depth of 5,704 reads using the estimate_richness function from phyloseq. Alpha diversity metrics were compared using analysis of variance (ANOVA) and Tukey's post hoc tests using vegan package P value correction for multiple testing. Multivariate analysis was performed after square root transforming the zOTU count data. Betadisper function was used to calculate the multivariate dispersion of samples (Bray-Curtis distance) between sample groups. Homogeneity of multivariate dispersion was tested with ANOVA. Nonmultidimensional scaling (nMDS) was performed to compare community compositions using the Bray-Curtis distance metric between sample groups. Permutational multivariate analysis of variance (PERMANOVA) with the "adonis" function (with 9,999 permutations) was used to statistically test for differences in community compositions between the back and cross transplant samples for each location as dispersion was significantly different between groups. Linear discriminant analysis effect size (LEfSe) implemented in the microbiomeMarker package in R was used to identify shifts in zOTUs between back-and cross-transplant samples for each location with a log Linear discriminant analysis (LDA) cutoff of 3 (Kruskal-Wallis test, P , 0.05). z score-transformed abundance profiles of marker zOTUs identified from LEfSe were visualized with a heatmap via pheatmap.
Environmental parameters. The water temperatures of the Outer (O) and Inner (I) Bays of the Penghu Archipelago, Taiwan, were obtained from May 2018 through December 2018 using temperature data loggers (HOBO Pendant, Onset Corp.) located ;3 m deep, close to target colonies, and recording temperatures every 30 min. Abiotic factors, including NH 3 , NO 3 , and PO 4 , were measured with a LaMotte 1910 SMART 3 colorimeter; pH was measured with a HORIBA LAQUA act water quality meter, and salinity was measured with an ATAGO master refractometer.
Bacteria isolation and culturing. "Ca. Endozoicomonas penghunesis" 4G was isolated from the coral Acropora muricata off the coast of the Inner Bay, Penghu Archipelago, Taiwan (GPS location, N23°3 1.8519 E119°33.6319). Coral tissue and mucus were sprayed with TE buffer (10 mM Tris-HCl, 1 mM EDTA, pH 8) and serially diluted to 10 24 . All dilutions were plated on modified marine broth version 4 (MMBv4 agar) (52) and incubated at 25°C. Each colony was screened first by the following primers: bacterial universal forward 27F (5'-AGA GTT TGA TCM TGG CTC AG-39) and Endozoicomonas-specific reverse En771R (59-TCA GTG TCA RRC CTG AGT GT-39) (86). Endozoicomonas species 16S rRNA gene V1-V4 region was PCR amplified by 35 cycles with a denaturing step at 94°C for the 30 s, followed by annealing at 54°C for 30 s and an extension step at 72°C for 45 s. PCR product was checked on a 1.5% agarose gel after electrophoresis. All samples with bands ;750 bp long were then subcultured in MMB medium. Fulllength 16S rRNA genes were amplified by universal bacterial primer 27F (59-AGA GTT TGA TCC TGG CTC AG-39) and 1541R (59-AAG GAG GTG ATC CAG CC-39). The full-length 16S rRNA PCR was amplified using 30 cycles with a denaturing step at 95°C for 30 s, annealing at 55°C for 30 s, and a final extension at 72°C for 90 s. Amplified products with target bands (;1,465 bp) were cut and later sequenced by Sanger sequencing (3730 DNA analyzer; Thermo) from Genomics, Taipei, Taiwan. Chromatograms obtained were manually checked and sequences were trimmed. The final length of the high-quality trimmed sequence was ;600 bp. Sequences with #98% identity to 16S rRNA genes of type strains from genus Endozoicomonas were deemed new candidates for novel Endozoicomonas species.
Physiological characterization. "Ca. Endozoicomonas penghunesis" 4G was cultivated on MMBv4 medium (52) (see Table S1 in the supplemental material) for the enrichment, and a broad range of physiological characterizations was performed. The optimum salinity was tested on MMB medium with NaCl concentrations adjusted as required (0.5% and 1.0 to ;4.0%, wt/vol, in increments of 1.0%).
The growth temperature range was tested at 4°C and 10 to 40°C (at 5°C intervals). The pH tolerance was determined using the following buffers: pH 4.0 to 7.0, HCl; pH 7.0 to 10.0, NaOH (at 1.0 pH unit intervals).
Three physiological tests (pH, salinity, and temperature) were measured based on the turbidity (at the optical density at 600 nm) of cultures grown at different pH values, NaCl concentrations, and temperatures, respectively. Commercial API 20NE kits (bioMérieux, France) were used to test the ability to metabolize different carbon substrates per the manufacturer's protocol. Additional carbon utilization was evaluated in the modified marine medium (see details in Table S2). Bacterial motility was tested in marine broth semisolid agar (0.5% agarose). The Gram stain kit (Fluka, England) was used to distinguish bacterial Gram reactions. Relation to oxygen was determined after incubating "Ca. Endozoicomonas penghunesis" 4G on MMB agar in the 2.5-liter Oxoid AnaroGen system (Thermo) and cultured at 25°C for 7 days. Oxidase and catalase activity was tested independently by adding 35% H 2 O 2 and 0.1% tetramethyl-phenylenediamine dihydrochloride, respectively. An antibiotic sensitivity test was performed after spreading bacteria on an MMB plate, with each disc containing different antibiotics (10 mg streptomycin and 10 mg ampicillin). The results were observed after 5 days of incubation at 25°C, and sensitivity was measured based on the distance from the discs to the edge of the clear zone. Bacteria were scored as sensitive if the diameter was greater than 2 mm, slightly sensitive if the diameter was 1 to 2 mm, and resistant otherwise.
Morphological characterization. The morphology of "Ca. Endozoicomonas penghunesis" 4G, including colony shape and color, was observed by a stereomicroscope (EZ4; Leica, Germany). General cell structure and cell inner structure were studied by transmission electron microscopy (TEM). The bacterial shape on a single colony was observed by scanning electron microscopy (SEM). TEM and SEM observations were made after bacteria were cultured in MMB for 1 day and MMB agar (1.5%) for 3 days, respectively. Colonies were incubated at 25°C.
For the "Ca. Endozoicomonas penghunesis" 4G thin section, bacteria were first centrifuged at 2,500 Â g for 5 min, and bacterial pellets were collected and fixed in 2.5% glutaraldehyde and 4% paraformaldehyde in a 0.1 M sodium phosphate buffer (pH 7.0) at room temperature for 1 h. After three 20-min buffer rinses, the samples were postfixed in 1% OsO 4 in the same buffer for 1 h at room temperature and then rinsed again as described above. Samples were dehydrated in an alcohol series, embedded in Spurr's resin (EMS), and sectioned with a Leica EM UC6 ultramicrotome (Leica, Germany). The ultrathin sections (70 to 90 nm) were stained with 5% uranyl acetate in 50% methanol and 0.4% lead citrate in 0.1 N sodium hydroxide. An FEI G2 Tecnai Spirit Twin TEM was used at 80 kV for viewing, and images were captured with a Gatan Orius charge-coupled device (CCD) camera.
The colony of "Ca. Endozoicomonas penghunesis" 4G was observed using cryo-SEM (FEI Quanta 200 SEM/Quorum Cryo System PP2000TR). The MMBv4 agar plate containing a single colony of "Ca. Endozoicomonas penghunesis" 4G was sectioned into 1 mm by 1 mm, loaded onto the medium-containing stub, and then frozen with liquid nitrogen slush. The frozen sample was transferred to the sample preparation chamber at 2160°C. After 5 min, the temperature was raised to 285°C, and the samples were etched for 20 min. After coating at 2130°C, the samples were transferred to the SEM chamber and observed at 2160°C and 20 kV.
The general cell morphology was studied by negative staining and observed under TEM. Bacteria were enriched in MMB for 1 day before adding a fixative solution (2.5% glutaraldehyde plus 4% paraformaldehyde-0.1 M PBS) at 37°C for 10 min. To reduce the background signal of TEM observation, MMB was replaced first by PBS and then by sterilized H 2 O twice, and the bacteria were mounted onto growdischarge carbon-Formvar grids. Bacteria were stained by 2% phosphotungstate for 1 s, and finally the sample was rinsed with sterilized H 2 O twice and viewed under an FEI G2 Tecnai Spirit Twin TEM at 80 KV. The images were then captured with a Gatan Orius CCD camera.
DPPH radical scavenging assay. The ability of free radical scavenging of two bacterial species was tested by a 2,2 diphenyl-1-picrylhydrazyl (DPPH) assay. The two Endozoicomonas species were incubated in MMbv4 medium with 0.1% glucose at 25°C and 30°C for 24 h, and the bacterial supernatant was harvested by centrifugation at 13,000 rpm for 5 min at 4°C. To collect cell lysate, cell pellets were resuspended by 0.8 ml 50 mM Tris-HCl lysis buffer (50 mM HEPES, 2% Triton X-100, 5% glycerol, and 40 mg 180-mm plastic beads) and loaded into a 2-mL O-ring-capped plastic microcentrifuge tube. The cell pellets with lysis buffer were vortexed at maximum intensity for 5 min and incubated on ice for 1 min. The step was repeated twice. After that, the cell lysates were centrifuged at 13,000 rpm for 10 min at 4°C to pellet down the cell debris and beads. The hydrogen atom donating ability of bacterial cultural supernatant and bacterial cell lysates was determined by decolorization of methanol solution of DPPH from violet color to yellow color in the presence of antioxidants. The 0.2 mM DPPH solution in 100% methanol was prepared, and 100 ml of the sample (supernatant or cell lysates) was mixed with 100 ml DPPH in 96-well plates in triplicate. Blanks were 100 ml sample mixed with 100 ml DPPH solvent (100% methanol), and control was 100 ml sample solvent (MMbv4 as the supernatant solvent and Tris-HCl lysis buffer as cell lysate sample) mixed with 100 ml DPPH solution. The reaction mixtures were incubated in the dark at 25°C for 30 min. The absorbance of the mixture was measured by a microplate reader (SpectraMax M2; Molecular Devices) at 517 nm. The DPPH scavenging activity was defined as DPPH free radical scavenging % ð Þ ¼ 1 2 ABS sample 2 ABS blank ABS control 2 ABS blank Â 100 (1) where ABS sample is the absorbance of the mixture of DPPH solution with samples, ABS blank is the absorbance of the mixture of samples and DPPH solvent, and ABS control is the absorbance of the mixture of sample solvent and DPPH solution.
Phylogenetic analysis of Endozoicomonas species zOTUs. To phylogenetically place the dominant Endozoicomonas zOTUs identified and the 16S rRNA gene of "Ca. Endozoicomonas penghunesis" 4G and to identify their closest neighbor within genus Endozoicomonas and its cultured isolates, representative 16S rRNA sequences from type strains (12 total) and one outgroup Halospina denitrificans HGD were downloaded from the NCBI taxonomy database (https://www.ncbi.nlm.nih.gov/taxonomy). Sequences were aligned using the RNA homology search tool cmalign (87) from the infernal package, and the CM models for domain bacteria were acquired from the rfam database (88). A maximum-likelihood phylogeny tree was built using the IQ-TREE web server (89) with 1,000 bootstraps and the best model selection enabled (best model, K2P1I1G4). The tree was finally visualized and edited in the laboratory-licensed version of iTOL v4 (90).
Long-and short-read paired-end sequencing and genome assembly. Long reads obtained from nanopore sequencing were first quality checked with nanoqc (91) and only high-quality paired-end reads were used for genome assembly using metaFlye (92) with default settings. Illumina reads (2Â 300) were first quality checked with FastQC (https://www.bioinformatics.babraham.ac.uk/projects/fastqc/), and then the adapters were removed and the reads trimmed with AdapterRemoval v2 (93). High quality (phred, .30) and trimmed paired-end reads were used to polish the crude nanopore assembly with four rounds of pilon (94) with default settings.
Genome annotation. The assembled genome was first checked for completeness, contamination, and heterogeneity using CheckM (40). The E. acroporae Acr-14 T genome (61) was assembled previously in our laboratory. Protein predictions in two Endozoicomonas species were performed with Prodigal in Prokka (95) with default settings to keep gene calls preserved for further functional categories analysis. A rapid annotation using subsystem technology (RAST) server (96) was used to obtain higher-order subsystem-level features. The "reconstruct pathway" approach in blastKOALA v2.2 (97) was used to obtain KEGG Ontology (KO) terms and in-depth annotation of the proteome. CRISPRcasFinder (98) was used to access the CRISPR-spacer. Eukaryote-like proteins (ELPs) were searched from a Batch Web-CD search against the CDD database (99), with minimum E value of 1e25 and maximum hit number set to 50. Circular genomic map of "Ca. Endozoicomonas penghunesis" 4G was visualized by CGView Server beta (100). Average nucleotide identity (ANI) analysis was performed to delineate the genome identity of "Ca. Endozoicomonas penghunesis" 4G against other Endozoicomonas species genomes which are publicly available from National Center for Biotechnology Information (NCBI) using the ANI calculator from the enveomics collection (http://enve-omics.ce.gatech.edu/ani/) (101).
Data availability. All sequencing data generated in the manuscript was submitted under the BioProject no. PRJNA758232, and the "Ca. Endozoicomonas penghunesis" genome was made available under GenBank accession no. SAMN21016876.

SUPPLEMENTAL MATERIAL
Supplemental material is available online only.