Similar but different: characterization of dddD gene mediated dimethylsulfoniopropionatemetabolism among Coral-Associated Endozoicomonas

doi:10.21203/rs.3.rs-3251001/v1

Download PDF

Research Article

Similar but different: characterization of dddD gene mediated dimethylsulfoniopropionatemetabolism among Coral-Associated Endozoicomonas

https://doi.org/10.21203/rs.3.rs-3251001/v1

This work is licensed under a CC BY 4.0 License

Version 1

posted

You are reading this latest preprint version

Endozoicomonas, the predominant bacterium associated with corals, is believed to play crucial roles in maintaining coral health. Its role in dimethylsulfoniopropionate (DMSP) degradation has been a subject of discussion for over a decade. A previous study found Endozoicomonas degraded DMSP through dddD pathway. This process releases dimethyl sulfide (DMS), which is vital for corals coping with thermal stress. However, little is known about the related gene regulation and metabolic abilities of DMSP metabolism in Endozoicomonadaceae. In this study, we observed a distinct DMSP metabolic trend in two phylogenetically close dddD -harboring Endozoicomonas species confirming genetically by the comparative transcriptomic profiling and visualizing the change of DMSP stable isotopes in bacterial cell using NanoSIMS. Furthermore, we found that DMSP cleavage enzymes are ubiquitous in coral Endozoicomonas with a preference for possessing DddD lyase. We speculate harboring DMSP degrading genes enables Endozoicomonas to successfully colonize various coral species across the globe.

Teaser: Coral-associated Endozoicomonas are widespread DMSP degraders monopolized by the dddD gene, employing various degradation mechanisms.

dddD gene

DMSP cleavage enzyme

Endozoicomonas

coral-associated bacteria

Bacteria from the genus Endozoicomonas (Gammaproteobacteria, Oceanospirillales, Endozoicomonadaceae) are predominantly and commonly found in a wide array of animal hosts living in coral reefs (1–6) and are regarded in most reports as health-related microbes, especially in corals (7–9). Putative metabolic interactions between Endozoicomonas and coral hosts have been proposed based on genomic information (2, 10–12). For example, genes involved in amino acid metabolism (2, 11, 13–15) and the enrichment of Type III and Type IV secretion system effectors (11, 12, 16) suggest that the bacteria might play a role in supplying amino acids to the coral holobiont (17) and use effectors to communicate with other associated organisms, especially their coral hosts. Additionally, some eukaryotic-like proteins (ELPs) are found in Endozoicomonas (11, 16), implying the possibility of co-evolution with their eukaryotic hosts.

In addition to these genomic factors suggesting symbiosis, evidence from more than a decade ago indicates that Endozoicomonas could be a putative dimethylsulfoniopropionate (DMSP) degrader (18, 19). Only recently has this function been experimentally confirmed in three E. acroporae strains with the discovery of a complete dddD gene-based operon, combined with biochemical analyses (11). DMSP is a ubiquitous sulfur-containing compound in marine environments that, once broken down, produces dimethyl sulfide (DMS). DMS acts as a sulfur-containing aerosol particle indirectly affecting Earth's climate and linking the marine and terrestrial sulfur cycles by participating in cloud formation (20). Beyond this significant role in the geochemical cycle, DMSP also serves multiple other functions. In marine algae, DMSP functions as an osmolyte (21), cryoprotectant (22), scavenger of reactive oxygen species (ROS)(23), and antioxidant (24). The increase in DMSP concentration in stressed coral (25, 26) makes DMSP a stress indicator (27) and links to its function as an antioxidant (28).

Bacteria are key players in DMSP degradation, which involves two metabolic pathways — demethylation and cleavage. The demethylation pathway, triggered by enzymes DmdA, -B, -C, and -D, is dominant and widely distributed across bacterial taxa, providing both sulfur and carbon sources. In contrast, the cleavage pathway, in which bacteria retain carbon but release sulfur into the ambient environment, generates DMS. There are eight DMSP lysis enzymes that can be grouped into three categories based on the byproducts they produce (please read the detailed review for more information regarding bacterial DMSP catabolism (29)). The biggest category of DMSP lysis enzymes, which generates cytotoxic acrylate (30) contains 6 DMSP lyases, including 5 from the Cupin superfamily (DddL, Q, W, K, Y) and 1 (DddP) from the M24 metallopeptidase family. The second category, which produces acryloyl-CoA, consists of a novel ATP-dependent lyase (DddX) from the Acyl-CoA synthetase superfamily (31). The third category, which mediates the additional acyl coenzyme (CoA-transferase) prior to the lyase activity and generates 3-hydroxypropionate (3HP) (32, 33) is called the CoA- DMSP lysis pathway and contains 1 lyase (DddD) from the Type III acyl CoA transferase family. The DMSP catabolic operon, facilitated by the dddD gene, was first identified in Marinomonas MWYL (33), and consists of five genes, including detS, adh, mms,gcv, and dddD. First, DMSP is imported through detS, which belongs to the betaine-carnitine-choline-transporter (BCCT) family. It is then followed by the enzymes adh and mms belonging to families of Fe-containing alcohol dehydrogenases and methylmalonatesemialdehyde dehydrogenases respectively. DMS will not undergo cleavage without its transcriptional regulator, gcv, which belongs to the LysR family.

Among the eight lyase enzymes, the dddD gene has a higher abundance than dddP (34) and is thought to be the main DMS producer via DMSP cleavage in coastal seawater, particularly when lysed by a member of the order Oceanospirillales (35). However, until now most of the studies related to DMSP have focused on the pelagic SAR11 clades of Alphaproteobacteria and phytoplankton-related Rhodobacterales, which are the main degraders of DMSP via demethylation in the ocean (35). Little is known about the DMSP-degrading bacteria in other DMSP-producing hotspots, such as coral reefs. Being one of the predominant coral-associated bacteria from the order Oceanospirillales, E. acroporae harbors one copy of the functional DMSP lyase-DddD (11), suggesting that Endozoicomonas might play a role in DMSP degradation of coral holobionts. However, here is not enough evidence to represent the cleavage ability, features, and gene regulation of DMSP metabolism in Endozoicomonas.

In this study, we successfully isolated a novel dddD gene harboring DMSP degrader, known as Candidatus Endozoicomonas ruthgatesiae strain 8E (hereafter “8E”). Genomes of dddD gene harboring Endozoicomonas were screened for DddD protein phylogenetic characteristic.

Through direct measurement of DMS by GC-MS, we observed a distinct metabolic profile within two Endozoicomonas species and visually confirmed this for 8E on a single-cell level using NanoSIMs Gene expression results indicated that 8E utilizes DMSP as a carbon source, as evidenced by the upregulation of main energy and carbohydrate related gene, whereas no similar upregulation was observed in E. acorporae. Based on these results, we propose that the distinct gene regulation and metabolic profiles observed in the utilization of DMSP by the two bacterial species reflect different functional roles, including differences in adaptation and evolution within DMSP metabolism. Our findings illuminate the mysterious role of Endozoicomonas in DMSP metabolism. Furthermore, based on the remarkable DMSP cleavage capacity in 8E and the metabolic differences between two Endozoicomonas species, we propose that contribution and functionality of Endozoicomonas on the coral sulfur cycle is more significant and diverse than previously thought.

Phenotypic and taxonomic traits of a novel coral-associated bacterium Candidatus Endozoicomonas ruthgatesiae 8E

8E was isolated from the tissue-mucus mixture of Acropora sp. coral from Kenting, Taiwan. Bacteria are rod-shaped (Fig. 1a), but their size and morphology change on different media (Fig. 1a, b, c). Peritrichous flagella-like and fimbria-like structures were only observed in minimum medium cultures with 0.1 mM DMSP (Fig. 1b) but not in Mmbv4. Otherwise, a few small granule-like structures in the middle of the cells have been observed (Fig. 1d). A similar structure has been found in another Endozoicomonas species (15). (Other detailed morphological and specific biochemical features are described in Supplementary Information and Table S1).

The phylogenetic analysis based on the 16S rRNA gene sequences (Fig. 1e) and the 92 bacterial core gene sets (Fig. 1f) of Endozoicomonadaceae both indicated that the closest species to 8E was E. numazuensis HC50^T (97.40% identity). An average nucleotide identity of 80.56% (ANI) (Fig. 1g) and an average amino acid identity of 81.35% (AAI) (Fig. 1h) also supported the result that E. numazuensis HC50^T is phylogenetically close to 8E.

The objective threshold for defining a species requires a 16S rRNA similarity of 97% (36)ANI (Average Nucleotide Identity) greater than 95%, and an Amino Acid Identity (AAI) greater than 62% (37, 38). 8E shared more than 80% of AAI with E. numazuensis HC50^T and E arenosclerae CABS5572^T, but only around 62% of AAI with other available Endozoicomonas genomes. This wide range of 62–100% suggests interspecies variation and puts 8E at the edge of a new genus classification (37, 38). These results, as well as the phenotypic characteristics, 16S rRNA gene phylogenetic tree and ANI calculation, lead us to propose that 8E is a novel species, Candidatus Endozoicomonas ruthgatesiae, in the genus Endozoicomonas.

Genome assembly and features of 8E

We recovered a single contig with a length of 7,377,917 bp using a hybrid sequencing strategy of ONT and Illumina sequencing platforms. This contig was larger than most of the other genomes detected in the genus Endozoicomonas (Table S2). Based on the analysis by CheckM, only 2 of 507 markers (c_Gammaproteobacteria [UID4444]) were missing, 502 genes were detected once, and 3 genes appeared twice, which resulted in 99.14% genome completeness, 0.54% contamination, and 0.00% strain heterogeneity. Overall, the quality of the 8E genome is ranked as “finished” according to the Genomic Standards Consortium (GSC) (39).

To clarify genome characteristics, comparative genomic analyses were conducted to detect unique features of 8E. Seven Endozoicomonas type strain genomes were included in the genomic comparison. Ten regions in the genome of 8E shared low similarity with other Endozoicomonas species (Fig. 2). In these 10 regions, we identified different clusters of genes which included WD40 domain protein, ankyrin related protein, midasin AAA ATPase, chromosome segregation protein/ATPase, glycotransferase family, and tandem repeat. Notably, these copies of genes usually gather adjacently in the 8E, especially the Eukaryotic-like proteins (ELPs) which we define as proteins found in the bacterial genome that may have originally in eukaryotes. We detected 160 WD40 domain proteins (WD40 superfamily) in the 8E genome (Table S2). Most of these were disseminated into 4 regions (Fig. 2) located at 0.197–0.331 Mbps (47 copies), 3.49–3.64 Mbps (34 copies), 4.62–4.78 Mbps (38 copies), and 5.50–5.79 Mbps (23 copies) (Fig. 2; Fig. S1a). A total of 66 gene copies of ankyrin-related proteins were found and annotated as ankyrin repeat (3 copies) or ankyrin-like protein (63 copies). Unlike the WD40 proteins, most of the ankyrin-related proteins were clustered in a single region around 0.436–0.744 Mbps (Fig. 2; Fig. S1b). Both types of protein mentioned above shared similar arrangements with E. acroporae Acr-14^T (ring 7) which had 23 WD40 proteins and 97 ankyrin repeat proteins (Table S2). Midasin AAA ATPase is a member of the MDN1 superfamily, which is a group of ELPs that are rarely observed in bacterial genomes. Midasin AAA ATPase was found to have 23 copies in the region from 0.966 to 1.23 Mbps (Fig. 2; Fig. S1c), suggesting that it may be functionally involved in ribosome maturation. It is not surprising that Endozoicomonas contained ELPs, but through gene annotation analysis using conserved domain search and EffectiveDB we found that the ELPs occupy up to 8.73% of the total 8E genome (Fig. S3) (E. acroporae 5.51%; E. montiporae 3.76%). The unexpectedly high percentage of ELPs in the genome is also found in another coral isolate E. marisrubri sp. 6C (16). However, the biological, mechanical, and original eukaryotic sources of these high proportions of ELPs in Endozoicomonas are still unclear (detailed descriptions of the functional annotation and genomic features of 8E are in “Supplementary Information”). No signs of genome reduction were found in the genome of 8E but genome expansion. We found multiple copies of chromosome segregation ATPases (33 copies), chromosome segregation proteins (16 copies), and glycosyltransferase (6 copies) which all copies are located neighboring. except that, there are a high density of various sized DNA Tandem repeats (1,126 TRs) ranging from 30 to 2,000 bp scattering between 4.94–5.07Mbp (74 TRs) and 5.84–6.31Mbp (361 TRs). 8E is a second big-genome Endozoicomonas (> 7MB), and with the advancements and widespread use of sequencing, we anticipate that more and more big genome Endozoicomonas will be sequenced. It is worth thoroughly investigating the genome evolution mechanism in those animal -associated bacteria with large genomes.

dddD gene- mediated DMSP lyases operon

Since 8E was isolated from the DMSP-rich coral genus Acropora, we paid attention to functions related to DMSP metabolism, as DMSP is an important molecule in the coral holobiont relative to scavenging free radicals in stressful conditions (40). Notably, the gene dddD encodes CoA-transferase lyase, which is a key participant in in the DMSP cleavage pathways. dddD was found in the 8E genome along with a complete DMSP lyase operon. We downloaded and annotated all Endozoicomonas genomic assembly on NCBI for the DddD protein sequence similarity comparison (Table S5). We built a DddD protein phylogenetic tree for Proteobacteria to locate the position of 8E in other DddD-harboring Proteobacteria (Fig. 3a). All 7 Endozoicomonas DddD sequences were clustered with other Gammaproteobacteria (Fig. 3a). Among them, 7 Endozoicomonas species separated into two different clades. Clade A contained 4 Endozoicomonas species, including E. acroporae, and clade B contained 3 Endozoicomonas species, including 8E (Fig. 3a). The DddD protein of 8E had a lower similarity with clade B Endozoicomonas (> 67.9% identity) than with Marinobacter zhenjiangenesis and Pseudomonas sp. SJZ079 (> 75.5% identity) (Table S3), suggesting that the two clades of Endozoicomonas strains might have different horizontal transfer sources. On the other hand, the gene arrangement and gene content of the dddD-mediated DMSP cleavage operon in the two clades are different (Fig. 3b). In clade A, E. acroporae and Endozoicomonas sp. G2_1 have the same gene arrangement, and both strains contain adhA (long-chain-alcohol dehydrogenase A). Clade B Endozoicomonas contains adhB (long-chain-alcohol dehydrogenase B), and their gene arrangement rather than close to phylogenetic close Endozoicomonas in Clade A but more similar to that of Pseudomonas guineae (Fig. 3b). Even through most of Endozoicomonas strains we included in the analysis are unculturable, we speculated that species placed in two different clades might possibly display similar DMSP metabolic features as E. acroporae or 8E.

8E is a potent DMS producer

To detect and compare the DMSP metabolic ability in 8E and other bacteria, we measured the production of DMS and consumption of DMSP in 8E, E. acroporae Acr-14^T, and a common species of DMSP degrader in the Rosebacter clade, Ruegeria atlantica DSM5823^T. The second time points of two Endozoicomonas represent the earliest detectable time of DMS after the addition of DMSP (Fig. 4a). Based on our observations using GC-MS, the overall trend of DMSP metabolism in the three bacteria can be categorized into two types, rapid emission, and gradual emission. At the first sampling point (6 h), the highest DMS peaks were detected from the culture of 8E with a consumption of 34.8% (from 0.1mM to 0.0655 \(\pm 0.00309 \text{m}\text{M}\)) of the DMSP in the medium (Fig. 4a). After the bacteria entered the stationary phase (Fig. 4b), the DMS_(g) concentration gradually decreased (from 5.84 to 5.53 nmol/ml).

Meanwhile, a different DMS production profile was observed in E. acroporae and R. atlantica. These species initially consumed more than half (0.0729 mM and 0.0876 mM) of the DMSP in the medium but emitted 5 and 9 times lower DMS_(g) concentration (1.07 and 0.056 nmol/ml respectively) than 8E (Fig. 4a). Additionally, the DMS concentration in these cultures continued to accumulate over time, reaching 56.1 nmol/ml at 24 h for E. acroporae and 11.15 nmol/ml at 24 h for R. atlantica.

To further confirm the DMSP cleavage of the rapid emission type in 8E, we used two different stable isotopes of DMSP (³⁴S-DMSP and 3¹³C -DMSP) (Fig. S4), and visualized the bacterial element ratios using NanoSIMS (Fig. 4c). ³⁴S was labelled to monitor the DMS emission, and 3¹³C -DMSP was used to monitor DMSP in the cell. At least 37 bacterial cells (regions of interest, ROIs) were analyzed in each isotopic treatment. The strong DMS emission caused no obvious ³⁴S enrichment in 8E. The highest ratio of ³⁴S/³²S was 4.50 (nature \(\delta\) ³⁴S/³²S \(\approx\)4.25) (Fig. 4d) after 6 h. In contrast, we observed enrichment of ¹³C at 6 h and 24 h, with 67 ROIs showing \(\delta\) ¹³C/¹²C higher than the nature \(\delta\) ¹³C/¹²C (nature value of \(\delta\) ¹³C/¹²C \(\approx\) 1.07) and 4 ROIs had two times higher than the natural \(\delta\) ¹³C/¹²C value (Fig. 4d). However, for both ³⁴S and ¹³C, the enrichment of isotope signals was not equally distributed in every cell. Instead, a few cells appeared brighter and had a stronger signal than the surrounding cells (Fig. 4c). Otherwise, the cell morphology was diverse. Cells with higher \(\delta\) ¹³C/¹²C values were mostly short and rod-shaped (Fig. 4c).

Transcriptional responses of Endozoicomonas species to DMSP

Based on the variation of DMSP utilization in the two Endozoicomonas species (Fig. 4a), we hypothesized that these two species might also display different gene regulations, especially in DMSP metabolism. We performed RNA-seq and collected transcripts at the time points when the DMS signal was first detected (8E at 6 h; E. acroporae at 8 h, both with and without addition of DMSP) (Fig. 4a). High-quality (Phred score > 35) sequences were also collected (Table S4). More than 10 million reads were obtained for all treatments, and more than 90% of the total reads were mapped to the two bacteria genomes. In total, around 90% of mapped reads were assigned to unique features (Table S4).

To understand the associated and subsequent responses of DMSP metabolism, we created a map of a putative metabolism scenario derived from all known intermediates of the dddD-centric process and initiated enzymatic reactions (Fig. 5).

The map starts from the dddD-mediated DMSP metabolic pathway and expands with each product. DMSP is a source of both sulfur and carbon in the demethylation pathway and is used only as a carbon source in cleavage enzyme-containing bacteria, including Endozoicomonas. The acetate produced from the dddD-mediated DMSP metabolic operon suggests a direct role of DMSP in bacterial central carbon metabolism. However, two bacteria displayed distinct transcript regulation after DMSP consumption, especially in energy and carbohydrate metabolism (Fig. 5, Fig. S3). Most of the genes related to the TCA cycle, glycolysis and gluconeogenesis were upregulated in 8E but not in E. acroporae, suggesting that DMSP may suppress the energy metabolism in E. acroporae. The upregulated transcript in amino acid metabolism was also differently in the two bacteria, valine and isoleucine biosynthesis pathway in 8E and proline (Fig. S3a), polyamine biosynthesis in E. acroporae which the average normalized count of the genes in. the pathway are 10 times more than the control group (Fig. S3b).

The metabolic pathway not directly linked with the DMSP operon also differs between the two bacteria. Pantothenate and aerobactin biosynthesis were significantly upregulated in E. acroporae. The transporter-related genes were also affected by DMSP addition; the reads of genes related to branched-chain and general L-amino acid ABC transporters were 10 times higher in 8E bacteria treated with DMSP than in 8E cultures that received no DMSP treatment (Fig. S3a). Further, in E. acroporae we found that four transporters (molybdate, methyl-galactoside, Fe (III) hydroxamate, and zine) were upregulated (Fig. S3b). Notably, secretion system structural genes and genes governing cell processes were expressed distinctly between the two bacteria. The increased of normalized read in type I secretion system, sec-SRP, quorum sensing, and flagellar wsd only observed in 8E (Fig. S3a). On the other hand, in E. acroporae, the type III secretion system and twin-arginine targeting (Tat) were found upregulated (Fig. S3b).

Endozoicomonas is a core bacterial group in the microbiota of coral and has been functionally associated with coral health in many studies over the past decade. However, their ecological functions are mostly unknown. In the present study, we discovered a novel, ubiquitous Endozoicomonas species, 8E, which possesses potent DMSP degradation activity via the DddD protein-mediated cleavage pathway, compared to E. acroporae (11) and the common oceanic DMSP utilizer, Ruegeria atlantica (41). DMSP cleavage activity was previously thought to be rare in Endozoicomonas, but the identification and proof of function in this study revealed that DMSP degradation is a common ecological function of Endozoicomonas species in coral reefs. Using comparative transcriptomic analysis, we uncorved variation in the regulation of dddD gene mediated DMSP metabolism in Endozoicomonas. This included pathways directly related and indirectly related, as well as other bacteria. These findings suggest that the underlying utilization mechanism of DMSP in bacteria is more complex, diverse, and crucial than previously thought.

Divergence in metabolic activities of DMSP degradation in two Endozoicomonas species

Within 12 hours of incubation, we clearly observed different profiles of DMSP degradation and uptake between 8E and E. acroporae. 8E displayed rapid and efficient DMSP cleavage, produced the highest DMS concentration (79.92 nmol/ml) within 6 hours of incubation, and consumed less than 35% of the supplied DMSP (Fig. 4a). Contrastingly, E. acroporae consumed more than 70% of the supplied DMSP in the first 8 hours and produced only a small amount of DMS (14.55 nmol/ml) (Fig. 4b). The underlying molecular mechanism causing the metabolic differences is still unclear, but three possible causes based on these results include: (i) Molecular variations at the gene and enzyme levels (i.e., variations in the dddD gene cluster and differentiations in enzymatic kinetics and activity); (ii) Variations in the physiological purpose of DMSP metabolism; (iii) More than one DMSP lyase gene participated in the DMSP degradation.

The DddD protein in the two Endozoicomonas species shares relatively low sequence identities with M. zhenjiangenesis and Pseudomonas sp. SJZ079 (Table S3), indicating the presence of molecular variations at the gene and enzyme levels. Furthermore, the phylogenetic relationship of the DddD protein and operon arrangement shows incongruence with bacterial taxonomy and phylogeny (Fig. 3a), indicating that the dddD gene clusters in both bacteria were horizontally acquired from different sources. Therefore, the different profiles of DMSP degradation may be caused by primary variations in the metabolic features of uniquely-sourced dddD genes or operons. Furthermore, the enzymatic kinetics and activity of the DddD protein between the two bacteria would be variable. We propose that the DddD enzyme might play a key role in the different metabolic profiles of the bacteria. The efficiency of DMSP cleavage activity was high in 8E (confirmed by the NanoSIMS results) which produced a DMS concentration 5 times higher than that of E. acroporae with only half the DMSP consumption (Fig. 4). This difference in observed efficiency suggests that the DddD lyase in 8E might have a high V_max value (indicating higher DMSP cleavage activity) and a low K_m value (indicating higher affinity to DMSP) than that in E. acroporae. Further investigation is required to confirm enzymatic characters. Cloning the purified DddD protein into a clear system such as E. coli and observing heterogeneous expression may be a direction for future study.

Apart from variations in the operons or genes, we also propose that variations in the physiological purpose of DMSP metabolism could lead to varying profiles of DMS production and DMSP uptake in the two bacteria. The upregulation of genes involved in the TCA cycle and amino acid metabolism in 8E following the addition of DMSP likely indicates that DMSP is used as a substitute carbon source for energy generation in poor nutrient conditions. On the other hand, the energy-related genes were commonly downregulated when DMSP was added in E. acroporae, suggesting that DMSP is used for non-energy production metabolisms. Similarly, we found that the proliferation efficiency of 8E was two times higher than E. acroporae with the addition of DMSP in the specific minimum medium, suggesting that DMSP may be used as an alternative energy source for 8E to propagate. The rate of proliferation of E. acroporae did not increase regardless of whether DMSP was added to the medium, indicating that the bacteria was under stressful conditions. Hence, E. acroporae likely uses DMSP to cope with stressful conditions. The bacterium might store DMSP and gradually cleave it to produce DMS. This strategy is common for DMSP users. Microorganisms accumulate DMSP for its beneficial antioxidant (24), osmolytic and cryoprotectant (42) properties, which might help soothe stressful physiology. Overall, though DMSP might mainly act as an antioxidant or other beneficial anti-stress molecule in E. acroporae, it could act as an energy source in 8E. Many studies have discussed how DMSP plays different physiological roles in different organisms (24, 29, 43). However, details comparative studies on the subject are rare, especially concerning only dddD-mediated DMSP utilization in closely related organisms.

After findings indicating that DMSP is not only the precursor of atmospheric trace gas DMS, but also plays multiple roles in different marine organisms (21, 42), increasing numbers of DMSP lyase have been identified. To date, nine DMSP lyase (including algal Alma1) have been found and categorized into three different lysis pathways based on their byproducts. There is a possibility of a novel degradation enzyme participating in Endozoicomonas DMSP degradation, as more than half of the genes (~ 2500 genes) are hypothetical. A future dddD gene knock-out experiment could be applied to clarify the hypothesis despite the lack of a successful gene knockout experiment conducted on Endozoicomonas to date.

Differential gene expression in DMSP pathways and relevant metabolisms in Endozoicomonas

According to the biochemical assay and the distinct gene regulation in two Endozoicomonas species that is associated with energy- and carbohydrate-related metabolic pathways, we believe the two species use DMSP for different biological purposes. This difference is also found when examining the pathway not directly linked with DMSP metabolism (layer 3, Fig. 5, S3). Genes for flagella, chemotaxis and quorum-sensing related pathways were upregulated in 8E. We observed a similar trend of gene regulation when the coral mucus-associated bacteria Roseobacter clades were exposed to DMSP (44–48). A flagella-like structure was observed on 8E during incubation in 0.1 mM DMSP in minimal medium (Fig. 1b), suggesting that flagella related physiological responses including chemotaxis and quorum-sensing are coupled with DMSP metabolism. However, a similar coupling response was not observed in E. acroporae, indicating that the gene upregulation concurrence of flagella, chemotaxis and quorum-sensing only appears in some DMSP bacterial degraders.

Additionally, the upregulated genes involved in aerobactin (a bacterial iron chelating agent, also called siderophore) biosynthesis and the Fe (III) hydroxamate transporter may increase resistance to oxidative stress in E. acroporae but not in 8E. For example, siderophores might increase resistance to oxidative stress in E. coli (49) and in the human pathogen Yersinia pseudotuberculosis (50). We speculate that aerobactin, a citrate-hydroxamate siderophore, could be another way for E. acroporae to mitigate oxidative stress.

dddD -mediated DMSP metabolism in coral reefs is more ecologically crucial than previously thought

The dominant DMSP lyase in coastal seawater is DddD protein, which is mainly found in Oceanospirillales (class Gammaproteobacteria). However, the ecological roles and influences of Oceanospirillales in DMSP-abundant coral reefs are still unclear. Endozoicomonas, being one predominant group of Oceanospirillales in coral, has long been thought to participate in coral sulfur cycling via DMSP metabolism (19, 51). With the framework of additional Endozoicomonas genomes sequenced in recent years, we inspected the gene contents of our sample and identified 10 Endozoicomonas species that contain DMSP cleavage genes (Table S6), indicating that DMSP degradation is a common eco-physiological function in Endozoicomonas and suggesting that this bacterial group plays a more crucial role in the sulfur cycle of coral holobionts or coral reefs than previously thought.

The DMS concentration produced by both Endozoicomonas species detected in our experiment was higher than that of R. atlantica. In out experiment, the cleavage ability of Endozoicomonas were outcompete to the widely documented DMSP degraders Roseobacter. In contrast, Ruegeria can tolerate a high DMSP concentration (5 mM) incubation and displayed higher DMS production at this high DMSP concentration than at a low DMSP concentration (0.1 mM) feeding incubation (52). In most cases, DMSP concentration in coral is less than 0.1 mM (40, 53, 54), and the highest DMSP concentration of 0.177 mM was detected in the coral species Acropora millepora under thermal stress (55). Our observations of the low DMS production rate of R. atlantica matched results from previous studies showing that Ruegeria species can shift the DMSP metabolic pathway from demethylation to cleavage along with an increase in DMSP concentration (from 0.1 mM to 5 mM DMSP) (52). Hence, we suggest that the Ruegeria in the stressed corals mostly metabolize DMSP via a non-DMS generated demethylation pathway. Unlike Ruegeria, Endozoicomonas cannot tolerate high DMSP concentrations, but can cope with 0.1 mM DMSP (Fig. S6), suggesting that a species of Endozoicomonas, rather than Ruegeria, could be the main producer of DMS in stressed corals.

We noticed that most Endozoicomonas genomes (9/10) accommodate the dddD gene (Table S6) rather than other DMSP cleavage genes. Compared to other DMSP cleavage pathways, dddD is unique and advantageous in the DMSP cleavage pathway as it produces 3-hydroxiopropionate instead of cytotoxic compound acrylate (30). DMSP cleavage does not accumulate harmful compounds in the cell, and also releases the potent antioxidant DMS, indicating that the selection of the dddD gene in Endozoicomonas may be associated with the establishment of a bacteria-coral mutualism. We performed a simple distribution analysis of 8E across the Indo-Pacific Ocean (Fig. S7) and found that 8E are widely distributed and have a similar distribution with E. acroporae (11). We suspect that the ability to cleave DMSP, enabled by the dddD gene, might be one of the reasons for the widespread appearance of Endozoicomonas, especially in the Acropora coral.

De novo Endozoicomonas isolation, genome sequencing and assembly.

The significant DMSP-producing coral genus, Acropora sp., (56) in Kenting, Taiwan, was selected to isolate a new species of Endozoicomonas. The methods used to isolate and culture this species have been previously described (6). A comprehensive description of the morphological, physiological, and biochemical characteristics of strain 8E is provided in the supplementary Information. The same method for bacterial genomic DNA isolation as described in (11) was employed in this study. High-quality genomic DNA was subjected to two different sequencing techniques. Illumina MiSeq paired-end sequencing was conducted at Yourgene Biosciences, Taiwan, while Nanopore sequencing was performed at the NGS High Throughput Genomics Core, Biodiversity Research Center, Academia Sinica, Taiwan. 8E genome assembly was obtained by combining high-accuracy Illumina short reads and Nanopore sequencing long reads. After conducting a quality check using nanoQC (57), the nanopore reads were assembled using Flye v2.9 (58). The resulting draft contigs were polished using Illumina reads through the utilization of the software, Pilon (59). After four rounds of polishing, we achieved 0 corrected SNPs, 0 ambiguous bases, 0 small insertions, and 0 deletions. This generated a high-quality contig with a length of 7,377,917 bp. CheckM V1.1.3 (60) was utilized to assess completeness, contamination, and heterogeneity. The method used to annotate the 8E genome is described in detail in “Supplementary Information.”

Phylogenetically comparing the dddD gene within Endozoicomonas to that of other bacteria

All Endozoicomonas dddD genes and operons used in this study were annotated using Prokka and a web-based CD-Search (61). To compare the position of the DddD protein within Endozoicomonas with that of other Proteobacteria, 67 DddD protein sequences were downloaded from the NCBI (keyword: dddD CoA-transferase/lyase) for phylogenetic tree construction. The collected DddD protein sequences were initially aligned using MAFFT v7. Subsequently, a maximum-likelihood tree was constructed using IQ-TREE (62) with 1000 replicates. Meanwhile, for the comparison of gene arrangement within the operon, we selected eight bacterial strains that contained the DddD protein from two different clades which were close to E. acorporae Acr-14^T and 8E, respectively. These bacterial strains included three from the class Gammaproteobacteria (Pseudomonas guineae, Marinobacter shejiangensis, and Marinomonas pollencensis), and 4 species from the same genus, Endozoicomonas (Endozoicomonas sp. ONNA1, Endozoicomonas sp. YOMI1, Endozoicomonas sp. G2-1 and E. acroporae Acr-14^T). The genomes used in this study were downloaded from the NCBI database and are listed in Table S5.

DMSP consumption and DMS quantified assay.

To ensure an equal initial bacterial count at the start of the experiment across different treatments, the initial bacterial population was standardized (OD600 = 0.6). Bacterial cultures were initially enriched in MMB medium supplemented with 0.1% glucose for 24 hours. After that, residual medium was removed via centrifuge, and the bacteria were then resuspended twice using fresh medium. After medium replacement, the bacteria were cultured in sterile 300 ml gas-tight vials, with rubber stoppers on the tops for sample collection (the bacterial concentration was measured and marked as BA; before acclimation). Before the addition of DMSP, the bacteria were acclimated in the fresh medium for 2 hours (h). We began collecting gas from the headspace of the bottles after DMSP was added (first measurement marked as 0 h).

Two different mediums, with and without subtract combination, were applied. As Endozoicomonas cannot grow in minimal medium with only DMSP, 0.2% casamino acid was added to all minimal medium treatments. The only difference between the control group (labeled as: mini) and the DMSP treatment group (label as: mini + DMSP) is the presence of 0.1 mM DMSP. The second medium used in this study was a rich medium, MMbv4. Unlike the minimal medium, no additional casamino acid was used in any of the MMbv4 medium treatment groups.

DMS (g) measurement used in the assay followed the protocol described in (11). Collection time points for the two Endozoicomonas species were chosen based on the initial emission time of DMS determined through our preliminary test (E. acroporae 8 h, 8E 6 h). We sampled R. atlantica DSS-3^T at both the 6 and 8 h time points for DMSP cleavage trend comparison.

DMSP concentration in the medium was determined via alkaline lysis. At each time point, 1.2 ml of the bacterial culture was drawn out from the gas-tight vials using a sterile syringe. The cultures were then centrifuged at 4500 rpm for 5 min, after which the resulting supernatant was collected and mixed with 100 ml 10M NaOH in a 20ml gas-tight vial. Afterwards, the vials were incubated at 200 rpm, 25 \(℃\) for 30 min to ensure the full lysis of DMSP to DMS. The target peak areas from both bacterial cleavage DMS in the headspace and alkaline lysis DMS from DMSP in the medium were quantified using a standard DMS/DMSP prepared calibration curve (detailed information can be found in Supplementary Information). Bacterial growth was monitored using OD600 (Analytik Jena, ScanDrop 250, Germany) and flow cytometry.

Stable-isotopic DMSP synthesis

To track the products derived from DMSP, a three-carbon 13 labeled [1,2,3-¹³C₃] stable-isotopic DMSP was essential (Fig. S4b). However, no such product was commercially available. We synthesized the stable isotopes by following a synthesis protocol of [1,2,3-¹³C₃] DMSP described in (63, 64). In brief, DMS was converted from DMSP hydrochloride (Supelco, USA) by alkaline hydrolysis. The resulting DMS was converted to [1,2,3-¹³C₃] DMSP hydrochloride via a Michael addition to [¹³C₃] acrylic acid (99 atom%, Sigma-Aldrich) under acidic conditions and washed with methylene chloride to remove any excess reactants. The structure and accurate mass of [1,2,3-¹³C₃] DMSP hydrochloride was confirmed by NMR and HRMS (Fig. S5). [³⁴S] DMSP (Fig. S4a) was provided by Professor William B. Whitman from the Department of Microbiology, University of Georgia.

Stable-isotopic labeled DMSP assay and NanoSIMS analysis

In the stable-isotopic labeled DMSP assay, the sampling time points, and bacterial incubation conditions were the same as those described for the DMSP consumption and DMS quantification assay. The two different stable isotopes labeled DMSP (Fig. S4) generated above were used as the DMSP source for further NanoSIMS visualization and analysis.

We filtered 10⁶ bacterial cellsthrough an Au-Pd (80:20) coated polycarbonate filter (Millipore, pore size 0.2 µm). Filters were washed by a series of ethanol solutions from 0%, 25%, 50%, 75–100% (v/v) for dehydration. After that, the isotopic compositions of bacterial cells on these filters were analyzed using a NanoSIMS 50 L (Cameca-Ametek, Gennevilliers, France) housed at Academia Sinica, Taiwan. Secondary ions of ¹²C₂⁻, ¹²C¹³C⁻, ¹²C¹⁴N⁻, ³²S⁻ and ³⁴S⁻ were collected simultaneously by multiple electron multipliers. Random areas on the polycarbonate filter were presputtered by 174–202 pA Cs⁺ ion beam for 1 minute and analyzed by 12.5–13.4 pA Cs⁺ ion beam for ~ 45–60 min to obtain 50 × 50 µm square image data at 512 × 512 pixels resolution. 30 mm and 150 mm of entrance slit, and aperture slit respectively were used to reach a mass resolving power of 5000 for 32S⁻. The images and isotope ratios of the cells were processed by L’Image software (developed by Larry Nittler, School of Earth and Space Exploration, Arizona State University). More than 37 cells were analyzed for each treatment. More than 2 × 10⁵ counts of ¹²C₂ and more than 10⁴ counts of ³²S⁻ were collected for each bacterial cell.

Acknowledgments

We thank the two microscopic facilities at Academia Sinica: the IPMB Cell Biology Core EM Division for the technical support in electronic microscopic imaging; and the Institute of Astronomy and Astrophysics NanoSIMS lab for stable-isotopic imaging analysis. We thank The NGS Core Facility at the Biodiversity Research Center, Academia Sinica for technical support of the high-throughput sequencing. Finally, we thank the research group of Dr. Huang from National Tsing Hua University and Prof. William B. Whitman from the Department of Microbiology, University of Georgia, for providing us with stable isotopes.

Funding:

Investigator Award [AS-IA-109-L05] Academia Sinica, Taipei, Taiwan

National Science and Technology council [MOST-110-2611-M-001-004]

Author contributions:

Conceptualization: Y.J.C , S.L.T

Performed research: Y.J.C, Y.F.C, S.P.P, C.Y.L, P.W.C, N.W

Data analysis: Y.J.C, S.S.Y.H, P.Y.L

Contributed data or analysis tools: Y.J.C, T.C.H, W.N.J, Y.L, D.C.L, Y.W.H

Supervision: S.L.T

Writing—original draft: Y.J.C , S.L.T

Writing—review & editing: Y.J.C, S.L.T

Competing interests: None of the authors have conflicts of interests to disclose.

Data and materials availability:

The bacterial genome assembly and the sequencing data of RNA-seq are available on NCBI under BioProject PRJNA946181. Genome assembly of 8E is available with accession number CP120717. All of the genome sequencing used in this study is listed in Supplementary Table S5. Previously published data used for analyses in this study are listed in “Supplementary Information”

S. Li, T. Young, S. Archer, K. Lee, S. Sharma, A. C. Alfaro, Mapping the Green-Lipped Mussel (Perna canaliculus) Microbiome: A Multi-Tissue Analysis of Bacterial and Fungal Diversity. Current microbiology 79, 1–8 (2022).
M. J. Neave, C. T. Michell, A. Apprill, C. R. Voolstra, Endozoicomonas genomes reveal functional adaptation and plasticity in bacterial strains symbiotically associated with diverse marine hosts. Scientific reports 7, 40579 (2017).
M. Nishijima, K. Adachi, A. Katsuta, Y. Shizuri, K. Yamasato, Endozoicomonas numazuensis sp. nov., a gammaproteobacterium isolated from marine sponges, and emended description of the genus Endozoicomonas Kurahashi and Yokota 2007. International journal of systematic and evolutionary microbiology 63, 709–714 (2013).
R. E. Pike, B. Haltli, R. G. Kerr, Description of Endozoicomonas euniceicola sp nov and Endozoicomonas gorgoniicola sp nov., bacteria isolated from the octocorals Eunicea fusca and Plexaura sp., and an emended description of the genus Endozoicomonas. International Journal of Systematic and Evolutionary Microbiology 63, 4294–4302 (2013).
L. Schreiber, K. U. Kjeldsen, M. Obst, P. Funch, A. Schramm, Description of Endozoicomonas ascidiicola sp. nov., isolated from Scandinavian ascidians. Systematic and applied microbiology 39, 313–318 (2016).
S. Y. Sheu, K. R. Lin, M. Y. Hsu, D. S. Sheu, S. L. Tang, W. M. Chen, Endozoicomonas acroporae sp nov., isolated from Acropora coral. International Journal of Systematic and Evolutionary Microbiology 67, 3791–3797 (2017).
J. L. Meyer, V. J. Paul, M. Teplitski, Community shifts in the surface microbiomes of the coral Porites astreoides with unusual lesions. PLoS One 9, e100316 (2014).
K. M. Morrow, D. G. Bourne, C. Humphrey, E. S. Botté, P. Laffy, J. Zaneveld, S. Uthicke, K. E. Fabricius, N. S. Webster, Natural volcanic CO2 seeps reveal future trajectories for host–microbial associations in corals and sponges. The ISME journal 9, 894–908 (2015).
Z. Qin, K. Yu, S. Chen, B. Chen, Q. Yao, X. Yu, N. Pan, X. Wei, Significant Changes in Bacterial Communities Associated with Pocillopora Corals Ingestion by Crown-of-Thorns Starfish: An Important Factor Affecting the Coral’s Health. Microorganisms 10, 207 (2022).
M. J. Neave, C. T. Michell, A. Apprill, C. R. Voolstra, Whole-genome sequences of three symbiotic Endozoicomonas strains. Genome Announc. 2, e00802-00814 (2014).
K. Tandon, C. Y. Lu, P. W. Chiang, N. Wada, S. H. Yang, Y. F. Chan, P. Y. Chen, H. Y. Chang, Y. J. Chiou, M. S. Chou, W. M. Chen, S. L. Tang, Comparative genomics: Dominant coral-bacterium Endozoicomonas acroporae metabolizes dimethylsulfoniopropionate (DMSP). Isme Journal, (2020).
J.-Y. Ding, J.-H. Shiu, W.-M. Chen, Y.-R. Chiang, S.-L. Tang, Genomic insight into the host–endosymbiont relationship of endozoicomonas montiporae CL-33T with its Coral Host. Frontiers in microbiology 7, 251 (2016).
K. Tandon, Y.-J. Chiou, S.-P. Yu, H. J. Hsieh, C.-Y. Lu, M.-T. Hsu, P.-W. Chiang, H.-J. Chen, N. Wada, S.-L. Tang, Microbiome restructuring: dominant coral bacterium Endozoicomonas species display differential adaptive capabilities to environmental changes. bioRxiv, (2021).
K. Tandon, Y.-J. Chiou, S.-P. Yu, H. J. Hsieh, C.-Y. Lu, M.-T. Hsu, P.-W. Chiang, H.-J. Chen, N. Wada, S.-L. Tang, Microbiome restructuring: dominant coral bacterium Endozoicomonas species display differential adaptive capabilities to environmental changes. bioRxiv, (2022).
K. Tandon, Y.-J. Chiou, S.-P. Yu, H. J. Hsieh, C.-Y. Lu, M.-T. Hsu, P.-W. Chiang, H.-J. Chen, N. Wada, S.-L. Tang, Microbiome Restructuring: Dominant Coral Bacterium Endozoicomonas Species Respond Differentially to Environmental Changes. Msystems, e00359-00322 (2022).
C. Pogoreutz, C. A. Oakley, N. Rädecker, A. Cárdenas, G. Perna, N. Xiang, L. Peng, S. K. Davy, D. K. Ngugi, C. R. Voolstra, Coral holobiont cues prime Endozoicomonas for a symbiotic lifestyle. The ISME Journal, 1–13 (2022).
L. M. Fitzgerald, A. M. SZMANT, Biosynthesis of ‘essential’amino acids by scleractinian corals. Biochemical Journal 322, 213–221 (1997).
J.-B. Raina, D. Tapiolas, B. L. Willis, D. G. Bourne, Coral-associated bacteria and their role in the biogeochemical cycling of sulfur. Applied and environmental microbiology 75, 3492–3501 (2009).
E. Ransome, S. J. Rowley, S. Thomas, K. Tait, C. B. Munn, Disturbance to conserved bacterial communities in the cold-water gorgonian coral Eunicella verrucosa. FEMS Microbiology Ecology 90, 404–416 (2014).
M. O. Andreae, Ocean-atmosphere interactions in the global biogeochemical sulfur cycle. Marine Chemistry 30, 1–29 (1990).
G. Kirst, Osmotic adjustment in phytoplankton and macroalgae: the use of dimethylsulfoniopropionate (DMSP). Biological and environmental chemistry of DMSP and related sulfonium compounds, 121–129 (1996).
G. Kirst, C. Thiel, H. Wolff, J. Nothnagel, M. Wanzek, R. Ulmke, Dimethylsulfoniopropionate (DMSP) in icealgae and its possible biological role. Marine Chemistry 35, 381–388 (1991).
A. M. Theseira, D. A. Nielsen, K. Petrou, Uptake of dimethylsulphoniopropionate (DMSP) reduces free reactive oxygen species (ROS) during late exponential growth in the diatom Thalassiosira weissflogii grown under three salinities. Marine Biology 167, 127 (2020).
W. Sunda, D. Kieber, R. Kiene, S. Huntsman, An antioxidant function for DMSP and DMS in marine algae. Nature 418, 317–320 (2002).
C. Aguilar, J.-B. Raina, C. A. Motti, S. Fôret, D. C. Hayward, B. Lapeyre, D. G. Bourne, D. J. Miller, Transcriptomic analysis of the response of Acropora millepora to hypo-osmotic stress provides insights into DMSP biosynthesis by corals. BMC genomics 18, 1–14 (2017).
S. G. Gardner, M. R. Nitschke, J. O'Brien, C. A. Motti, J. R. Seymour, P. J. Ralph, K. Petrou, J.-B. Raina, Increased DMSP availability during thermal stress influences DMSP-degrading bacteria in coral mucus. Frontiers in Marine Science, (2022).
S. G. Gardner, D. A. Nielsen, O. Laczka, R. Shimmon, V. H. Beltran, P. J. Ralph, K. Petrou, Dimethylsulfoniopropionate, superoxide dismutase and glutathione as stress response indicators in three corals under short-term hyposalinity stress. Proceedings of the Royal Society B: Biological Sciences 283, 20152418 (2016).
E. S. Deschaseaux, G. B. Jones, M. A. Deseo, K. M. Shepherd, R. Kiene, H. Swan, P. L. Harrison, B. D. Eyre, Effects of environmental factors on dimethylated sulfur compounds and their potential role in the antioxidant system of the coral holobiont. Limnology and Oceanography 59, 758–768 (2014).
C. R. Reisch, M. A. Moran, W. B. Whitman, Bacterial catabolism of dimethylsulfoniopropionate (DMSP). Frontiers in microbiology 2, 172 (2011).
E. Yoshii, Cytotoxic effects of acrylates and methacrylates: relationships of monomer structures and cytotoxicity. Journal of Biomedical Materials Research: An Official Journal of The Society for Biomaterials and The Japanese Society for Biomaterials 37, 517–524 (1997).
C.-Y. Li, X.-J. Wang, X.-L. Chen, Q. Sheng, S. Zhang, P. Wang, M. Quareshy, B. Rihtman, X. Shao, C. Gao, A novel ATP dependent dimethylsulfoniopropionate lyase in bacteria that releases dimethyl sulfide and acryloyl-CoA. Elife 10, e64045 (2021).
U. Alcolombri, P. Laurino, P. Lara-Astiaso, A. Vardi, D. S. Tawfik, DddD is a CoA-transferase/lyase producing dimethyl sulfide in the marine environment. Biochemistry 53, 5473–5475 (2014).
J. D. Todd, R. Rogers, Y. G. Li, M. Wexler, P. L. Bond, L. Sun, A. R. Curson, G. Malin, M. Steinke, A. W. Johnston, Structural and regulatory genes required to make the gas dimethyl sulfide in bacteria. Science 315, 666–669 (2007).
Y. Cui, S.-K. Wong, R. Kaneko, A. Mouri, Y. Tada, I. Nagao, S.-J. Chun, H.-G. Lee, C.-Y. Ahn, H.-M. Oh, Distribution of dimethylsulfoniopropionate degradation genes reflects strong water current dependencies in the Sanriku coastal region in Japan: from mesocosm to field study. Frontiers in microbiology 11, 1372 (2020).
J. Liu, C.-X. Xue, J. Wang, A. T. Crombie, O. Carrión, A. W. Johnston, J. C. Murrell, J. Liu, Y. Zheng, X.-H. Zhang, Oceanospirillales containing the DMSP lyase DddD are key utilisers of carbon from DMSP in coastal seawater. Microbiome 10, 1–21 (2022).
E. Stackebrandt, B. M. GOEBEL, Taxonomic note: a place for DNA-DNA reassociation and 16S rRNA sequence analysis in the present species definition in bacteriology. International journal of systematic and evolutionary microbiology 44, 846–849 (1994).
M. R. Olm, A. Crits-Christoph, S. Diamond, A. Lavy, P. B. Matheus Carnevali, J. F. Banfield, Consistent metagenome-derived metrics verify and delineate bacterial species boundaries. Msystems 5, e00731-00719 (2020).
M. Kim, H.-S. Oh, S.-C. Park, J. Chun, Towards a taxonomic coherence between average nucleotide identity and 16S rRNA gene sequence similarity for species demarcation of prokaryotes. International journal of systematic and evolutionary microbiology 64, 346–351 (2014).
R. M. Bowers, N. C. Kyrpides, R. Stepanauskas, M. Harmon-Smith, D. Doud, T. Reddy, F. Schulz, J. Jarett, A. R. Rivers, E. A. Eloe-Fadrosh, Minimum information about a single amplified genome (MISAG) and a metagenome-assembled genome (MIMAG) of bacteria and archaea. Nature biotechnology 35, 725–731 (2017).
J.-B. Raina, D. M. Tapiolas, S. Forêt, A. Lutz, D. Abrego, J. Ceh, F. O. Seneca, P. L. Clode, D. G. Bourne, B. L. Willis, DMSP biosynthesis by an animal and its role in coral thermal stress response. Nature 502, 677–680 (2013).
J. M. Gonzalez, J. S. Covert, W. B. Whitman, J. R. Henriksen, F. Mayer, B. Scharf, R. Schmitt, A. Buchan, J. A. Fuhrman, R. P. Kiene, Silicibacter pomeroyi sp. nov. and Roseovarius nubinhibens sp. nov., dimethylsulfoniopropionate-demethylating bacteria from marine environments. International journal of systematic and evolutionary microbiology 53, 1261–1269 (2003).
U. Karsten, K. Kück, C. Vogt, G. Kirst, "Dimethylsulfoniopropionate production in phototrophic organisms and its physiological functions as a cryoprotectant" in Biological and environmental chemistry of DMSP and related sulfonium compounds (Springer, 1996), pp. 143–153.
U. Karsten, G. Kirst, C. Wiencke, Dimethylsulphoniopropionate (DMSP) accumulation in green macioalgae from polar to temperate regions: interactive effects of light versus salinity and light versus temperature. Polar biology 12, 603–607 (1992).
T. R. Miller, K. Hnilicka, A. Dziedzic, P. Desplats, R. Belas, Chemotaxis of Silicibacter sp. strain TM1040 toward dinoflagellate products. Applied and environmental microbiology 70, 4692–4701 (2004).
T. R. Miller, R. Belas, Motility is involved in Silicibacter sp. TM1040 interaction with dinoflagellates. Environmental microbiology 8, 1648–1659 (2006).
M. Landa, A. S. Burns, S. J. Roth, M. A. Moran, Bacterial transcriptome remodeling during sequential co-culture with a marine dinoflagellate and diatom. The ISME journal 11, 2677–2690 (2017).
J. Li, W. Kuang, L. Long, S. Zhang, Production of quorum-sensing signals by bacteria in the coral mucus layer. Coral Reefs 36, 1235–1241 (2017).
W. M. Johnson, M. C. Kido Soule, E. B. Kujawinski, Evidence for quorum sensing and differential metabolite production by a marine bacterium in response to DMSP. The ISME journal 10, 2304–2316 (2016).
C. Adler, N. S. Corbalan, D. R. Peralta, M. F. Pomares, R. E. de Cristóbal, P. A. Vincent, The alternative role of enterobactin as an oxidative stress protector allows Escherichia coli colony development. PloS one 9, e84734 (2014).
C. Li, D. Pan, M. Li, Y. Wang, L. Song, D. Yu, Y. Zuo, K. Wang, Y. Liu, Z. Wei, Aerobactin-mediated iron acquisition enhances biofilm formation, oxidative stress resistance, and virulence of Yersinia pseudotuberculosis. Frontiers in Microbiology 12, (2021).
J.-B. Raina, D. Tapiolas, B. L. Willis, D. G. Bourne, Coral-associated bacteria and their role in the biogeochemical cycling of sulfur. Appl. Environ. Microbiol. 75, 3492–3501 (2009).
J. S. Wirth, T. Wang, Q. Huang, R. H. White, W. B. Whitman, Dimethylsulfoniopropionate sulfur and methyl carbon assimilation in Ruegeria species. MBio 11, e00329-00320 (2020).
A. D. Broadbent, G. B. Jones, R. J. Jones, DMSP in corals and benthic algae from the Great Barrier Reef. Estuarine, Coastal and Shelf Science 55, 547–555 (2002).
D. M. Tapiolas, J.-B. Raina, A. Lutz, B. L. Willis, C. A. Motti, Direct measurement of dimethylsulfoniopropionate (DMSP) in reef-building corals using quantitative nuclear magnetic resonance (qNMR) spectroscopy. Journal of Experimental Marine Biology and Ecology 443, 85–89 (2013).
S. G. Gardner, M. R. Nitschke, J. O'Brien, C. A. Motti, J. R. Seymour, P. J. Ralph, K. Petrou, J.-B. Raina, Increased DMSP availability during thermal stress influences DMSP-degrading bacteria in coral mucus. Frontiers in Marine Science, 1684 (2022).
K. L. Van Alstyne, P. Schupp, M. Slattery, The distribution of dimethylsulfoniopropionate in tropical Pacific coral reef invertebrates. Coral Reefs 25, 321–327 (2006).
W. De Coster, S. D’Hert, D. T. Schultz, M. Cruts, C. Van Broeckhoven, NanoPack: visualizing and processing long-read sequencing data. Bioinformatics 34, 2666–2669 (2018).
M. Kolmogorov, D. M. Bickhart, B. Behsaz, A. Gurevich, M. Rayko, S. B. Shin, K. Kuhn, J. Yuan, E. Polevikov, T. P. Smith, metaFlye: scalable long-read metagenome assembly using repeat graphs. Nature Methods 17, 1103–1110 (2020).
B. J. Walker, T. Abeel, T. Shea, M. Priest, A. Abouelliel, S. Sakthikumar, C. A. Cuomo, Q. Zeng, J. Wortman, S. K. Young, Pilon: an integrated tool for comprehensive microbial variant detection and genome assembly improvement. PloS one 9, e112963 (2014).
D. H. Parks, M. Imelfort, C. T. Skennerton, P. Hugenholtz, G. W. Tyson, CheckM: assessing the quality of microbial genomes recovered from isolates, single cells, and metagenomes. Genome research 25, 1043–1055 (2015).
A. Marchler-Bauer, S. Lu, J. B. Anderson, F. Chitsaz, M. K. Derbyshire, C. DeWeese-Scott, J. H. Fong, L. Y. Geer, R. C. Geer, N. R. Gonzales, CDD: a Conserved Domain Database for the functional annotation of proteins. Nucleic acids research 39, D225-D229 (2010).
B. Q. Minh, H. A. Schmidt, O. Chernomor, D. Schrempf, M. D. Woodhams, A. von Haeseler, R. Lanfear, IQ-TREE 2: New Models and Efficient Methods for Phylogenetic Inference in the Genomic Era. Molecular Biology and Evolution 37, 1530–1534 (2020).
S. T. Chambers, C. M. Kunin, D. Miller, A. Hamada, Dimethylthetin can substitute for glycine betaine as an osmoprotectant molecule for Escherichia coli. Journal of Bacteriology 169, 4845–4847 (1987).
J. S. Wirth, W. B. Whitman, An efficient method for synthesizing dimethylsulfonio-34S‐propionate hydrochloride from 34S8. Journal of Labelled Compounds and Radiopharmaceuticals 62, 52–58 (2019).

Supplementary Information not available with this version.

Download PDF

Version 1

posted

You are reading this latest preprint version

Similar but different: characterization of dddD gene mediated dimethylsulfoniopropionatemetabolism among Coral-Associated Endozoicomonas

Status:

Version 1

Abstract

Figures

Introduction

Results

Genome assembly and features of 8E

8E is a potent DMS producer

Discussion

Materials and Methods

Stable-isotopic DMSP synthesis

Stable-isotopic labeled DMSP assay and NanoSIMS analysis

Declarations

References

Supplementary Informations

Status:

Version 1