Widespread occurrence of chitinase-encoding genes suggests the Endozoicomonadaceae family as a key player in chitin processing in the marine benthos

Introduction

Chitin, the unbranched polymer of β−1,4-N-acetylglucosamine, is the second-most abundant natural polysaccharide on earth after cellulose, with global production rates estimated at 1012–1014 tonnes per year [1, 2]. It is a major constituent of fungal cell walls, marine diatoms, crustaceans, and zooplankton, and difficult to decompose in nature [1, 3]. The water insolubility of chitin contributes to the formation of ‘marine snow’ in marine ecosystems [4]. Nevertheless, chitin particles rarely accumulate on the seafloor as they are degraded by bacteria which use chitin as carbon and nitrogen sources [5]. Filter- and suspension-feeding animals, such as sponges and corals, likely process numerous chitin-rich particles present in the water column and may thus host a microbiome well adapted to chitinous food processing [6].

The gammaproteobacterial family Endozoicomonadaceae, first described in 2018, has a worldwide distribution and lives associated with marine eukaryotic organisms, including corals, sponges, mollusks, bryozoans, ascidians, and fishes [7]. This family currently comprises three validly published genera: Endozoicomonas, Parendozoicomonas and Kistimonas [7]. Based on metagenome-assembled-genome (MAGs) data, we recently proposed the novel genus Candidatus Gorgonimonas along with two candidate species, Gorgonimonas eunicellae and Gorgonimonas leptogorgiae, found to be specifically associated with temperate octocorals (Octocorallia) [8]. Endozoicomonadaceae symbionts can represent over 90% of a coral bacterial community [9] and are presumed indicators of coral health [9, 10]. Owing to their abundance in corals and ubiquity in marine ecosystems, the roles played by Endozoicomonadaceae symbionts of marine animals are of increasing research interest yet their precise functions across multiple symbioses are still poorly understood. A major limitation is the small number of cultured representatives available [11] and the apparent difficulty to manipulate these symbionts in the laboratory. Genome analyses suggested that Endozoicomonadaceae symbionts participate in amino acid and B vitamin supply [8, 12,13,14], micronutrient acquisition [8] and sulfur cycling [15]. Recent metagenomics and genomics data also indicate a role in carbohydrate degradation [14], including chitin [8]. However, studies performed so far were mostly based on the inspection of one or few Endozoicomonadaceae species, while family-wide, comparative surveys with multiple genera are lacking. Hence, this study inspected 42 Endozoicomonadaceae genomes (i.e., all genomes publicly available until May 2022), including cultured and uncultured lineages and all formally described genera, to determine the distribution and diversity of endo-chitinase-encoding genes and whether chitin hydrolysis is indeed a widespread capacity in this family. The full methodology employed in this study is documented in Additional file 1. Details on the genomes examined here and their annotations are available in Additional file 2.

Results and discussion

Genes coding for endo-chitinases (EC 3.2.1.14), the enzymes that cleave chitin polymers into oligomers [16], were found on 32 of 42 Endozoicomonadaceae genomes (Fig. 1 and Table S1), including representatives of all formally described genera and all Candidatus Gorgonimonas MAGs. Several genomes harbored more than one endo-chitinase encoding gene, resulting in 57 such genes detected across 32 genomes with a coding potential for chitin hydrolysis (RASTtk annotation, see also Table S2 for details). Of the 57 genes, 40 genes were found to code for glycoside hydrolase family GH18 endo-chitinases, the most common type of chitinase in bacteria [1]. In addition, we found four GH19-type endo-chitinases, one dual GH18-GH19 endo-chitinase and two GH6-type enzymes possibly indicating cellulase activity, while no GH domain was detected on the remaining 10 sequences. Twenty-nine Endozoicomonadaceae genomes also had genes coding for exo-chitinases/β-N-acetylglucosaminidases (EC 3.2.1.52, RASTtk annotation), which can cleave chito-oligomers to produce the mono-sugar N-acetylglucosamine [16]. Moreover, 26 genomes harbored polysaccharide deacetylases, indicating that Endozoicomonadaceae may also convert chitin into chitosan [1], and 13 had chitin-binding protein-encoding genes, revealing indeed a versatile chitin degradation machinery in this marine bacterial family.

Fig. 1: Chitin-degradation related genes present in the Endozoicomonadaceae family.
figure 1

The Maximum Likelihood phylogenomic tree of the Endozoicomonadaceae family on the left was constructed with the Insert Genomes into Species Tree—v.2.2.0 app of DOE Systems Biology Knowledgebase (KBase). This tree presents 42 genomes, including 19 metagenome-assembled genomes (MAGs, in bold), two single amplified genomes (SAGs, in italics), and 21 genomes from cultured representatives. Gray circles on the branches indicate bootstrap support of ≥75%. The tree is drawn to scale and was edited in iTOL. The columns next to the tree indicate the number of chitin-degradation related proteins/protein domains present on each genome. Chitin-degradation related genes were annotated using RASTtk v.2.0 (Rapid Annotation using Subsystem Technology toolkit; black circles), COG (Cluster of Orthologous Groups of proteins, empty triangles) and Pfam (protein families, empty circles) databases, respectively. Ca.—Candidatus. More information on the genomes used here, including their genome assembly accession numbers, is available in Table S1 of Additional File 2. Genome-wide COG and Pfam annotations of the 42 Endozoicomonadaceae genomes are available in Tables S3 and S4. Further details on the chitin degradation-related entries shown here can be found in Table S2 (RASTtk) and Table S5 (COGs, Pfams).

Full size image

Phylogenomic inference of the Endozoicomonadaceae family (Fig. 1) showed that strains without endo-chitinase encoding genes were usually not of the same species as those that carried endo-chitinases, indicating interspecific trait variation. Habitat origin, on the other hand, did not seem to be a determining factor, since both, genomes with and without endo-chitinases, derived from a variety of host animals and even non-host environments (Table S1). Of the 10 genomes that lacked endo-chitinase genes, five were MAGs and five were genomes from cultured representatives (Fig. 1). Likewise, genomes from cultured representatives with >99% genome completeness (e.g., Endozoicomonas sp. G2 1, Endozoicomonadaceae bacterium GTF-13, Kistimonas asteriae KMD 001) were found not to carry an endo-chitinase whereas multiple MAGs with 60–70% genome completeness did possess endo-chitinases (Table S1), indicating that neither cultivation status nor assembly completeness were determining factors for endo-chitinase presence across the surveyed genomes.

Nevertheless, the number of chitin-degradation related entries varied somewhat depending on the annotation tool/database used. For example, Clusters of Orthologous Groups of proteins (COG) and Protein families (Pfam) annotations identified 41 and 58 endo-chitinase domains (compared to 57 identified by RASTtk) on 26 out of 42 Endozoicomonadaceae genomes (Fig. 1, Tables S3–S5). The number of predicted exo-chitinases was lower with COG (n = 9) and Pfam (n = 16) annotations than with RASTtk (n = 29), while the number of chitin/carbohydrate binding domains (25 and 28) and polysaccharide deacetylases (73 and 53) was higher with COG and Pfam than with RASTtk (13 and 26). Quite congruent were the endo-chitinase predictions for the 11 Candidatus Gorgonimonas MAGs with all annotation schemes predicting one endo-chitinase per MAG (Fig. 1 and Table S1).

A phylogenetic analysis (Fig. 2 and Table S6) of 37 full-length, translated amino acid sequences showed that endo-chitinases of the 11 Candidatus Gorgonimonas MAGs formed a separate, well-supported cluster, unique of octocorals, and congruent with the distinct phylogenomic grouping of this candidate genus within the family (Fig. 1). Resembling species-level phylogeny, Candidatus Gorgonimonas endo-chitinases formed two subclusters (Fig. 2), indicating that each Candidatus Gorgonimonas species has its specific endo-chitinase. These endo-chitinases differed in amino acid (524 versus 577 aa) sequence lengths, predicted molecular weight (54.89 versus 61.29 kDa) and isoelectric point (7.08 versus 4.77), and by the presence of two Domains of Unknown Function (DUFA5011) in Candidatus Gorgonimonas leptogorgiae, which were absent in Candidatus Gorgonimonas eunicellae. Following a Blastp search, the closest endo-chitinases present in public databases were GH18 and DUF5011 domain-containing proteins from Photobacterium spp. (Vibrionaceae, Gammaproteobacteria). However, amino acid sequence similarity was low (~50%, see Table S6 for details), underlining the uniqueness of Candidatus Gorgonimonas endo-chitinases within the Endozoicomonadaceae family and Gammaproteobacteria in general.

Fig. 2: Phylogenetic analysis of endo-chitinases (EC 3.2.1.14) in the Endozoicomonadaceae family.
figure 2

The tree was produced using Mega11 with 37 endo-chitinase amino acid sequences, translated from full-length endo-chitinase gene sequences obtained from Endozoicomonadaceae genomes. Endo-chitinases from a MAG are highlighted in bold, the remaining endo-chitinase genes derived from cultured representatives. Sequences were aligned using MAFFT (Multiple Alignment using Fast Fourier Transform) v. 7.48. The evolutionary history was inferred using the Maximum Likelihood method based on a Whelan And Goldman (WAG) model. The tree with the highest log likelihood (−10396.47) is shown. A discrete Gamma distribution was used to model evolutionary rate differences among sites (5 categories (+G, parameter = 50.1306)). The rate variation model allowed for some sites to be evolutionarily invariable ([+I], 0.44% sites). The tree is drawn to scale, with branch lengths measured in the number of substitutions per site. There were a total of 340 amino acid positions in the final dataset. The tree was constructed with 1000 bootstrap repetitions, and values on the branches indicate bootstrap support of ≥60%. Blue stars on cluster nodes highlight genus-specific (Candidatus Gorgonimonas and Endozoicomonas, respectively) endo-chitinase clusters. Protein domains present on the endo-chitinase sequences are indicated by colored dots as identified in the legend below the tree. Black stars indicate the presence of a signal peptide and white stars the prediction of an active site. The theoretical average molecular weight (in kDa) and isoelectric point (pI) of each endo-chitinase (after removal of the corresponding signal peptide domain) is displayed on the right. The colored background shading indicates the origin of the Endozoicomonadaceae genomes as detailed in the legend below the tree. More information about the endo-chitinase sequences shown here is available in Table S6 of Additional file 2.

Full size image

Eighty-four percent of the Endozoicomonadaceae endo-chitinases depicted in Fig. 2 contained a signal peptide domain, responsible for transporting the enzyme outside of the cell [17] where endo-chitinases typically degrade the large, water-insoluble chitin fibers [1]. Moreover, 86% of these Endozoicomonadaceae endo-chitinases had an active site where the chitinous substrate is bound and the catalytic reaction initiated [18] (Fig. 2), suggesting that the endo-chitinases are functional. Several Endozoicomonadaceae genomes harbored two or three distinct endo-chitinases that differed significantly in their protein domains, molecular weight (from 20 kDa to 118 kDa) and isoelectric points (from 4.15 to 9.04) (Fig. 2). Endozoicomonas montiporae, Endozoicomonas sp. SM1973 and Paraendozoicomonas haliclonae, for example, harbored both GH18 and GH19-type chitinases. Possibly, these differences in protein sequence domains, molecular weight, and isoelectric points translate into distinct substrate affinities and/or enzymatic properties, such as different temperature and/or pH optima, or heavy metal tolerance, as recently observed in chitinases from Vibrio species [19]. Such chitinolytic plasticity could provide an advantage for Endozoicomonadaceae bacteria to cope with environmental changes in their natural habitat and with unique physical-chemical conditions and substrate quality and availability across different host-microbe symbioses.

We demonstrate that chitinase-encoding genes are widespread among Endozoicomonadaceae family members and frequently accompanied by auxiliary genes related to chitin breakdown, suggesting that these bacteria are well-adapted to process chitin-containing matrices in diverse marine micro-niches. Endozoicomonadaceae symbionts are thus likely important contributors to chitin turnover in filter- and suspension-feeding marine animals, capable of providing energetic, carbon and nitrogen-rich mono- and oligo-sugars to the holobiont. Given their ubiquity in the marine benthos, Endozoicomonadaceae bacteria could play key roles in carbon and nitrogen cycling on reefs and in marine ecosystems at large.

Studies of the kinetics and biochemical properties of endo-chitinases from the Endozoicomonadaceae family, including heterologously expressed endo-chitinases from uncultured organisms, such as Candidatus Gorgonimonas symbionts of octocorals, are needed to deepen our understanding of their physiology and function. This knowledge may open future opportunities for the circular blue bioeconomy sector since chitinases and chitin-derived products find various biotechnological and pharmaceutical applications [2, 3]. Moreover, future in-vivo studies using stable isotope probing will be helpful to determine the fate of chitin and to reveal chitin-associated cross-feeding mechanisms and food webs in coral and sponge holobionts. Such insights will contribute to a more holistic understanding of symbiont ecology and carbon and nitrogen fluxes in the marine realm.

Related Articles

Protein signatures predict coral resilience and survival to thermal bleaching events

Coral bleaching events from thermal stress are increasing globally in duration, frequency, and intensity. While bleaching can cause mortality, some corals survive, reacquire symbionts, and recover. We experimentally bleached Montipora capitata to examine molecular and physiological differences between corals that recover (resilient) and those that die (susceptible). Corals were collected and monitored for eight months post-bleaching to identify genets with long-term resilience. Using an integrated systems-biology approach that included quantitative proteomics, 16S rRNA sequencing to characterize the coral microbiome, total coral lipids, symbiont community composition and density, we explored molecular-level mechanisms of tolerance in corals pre- and post-bleaching. Prior to thermal stress, resilient corals have a more diverse microbiome and abundant proteins essential for carbon acquisition, symbiont retention, and pathogen resistance. Protein signatures of susceptible corals showed early symbiont rejection and utilized urea for carbon and nitrogen. Our results reveal molecular factors for surviving bleaching events and identify diagnostic protein biomarkers for reef management and restoration.

Microbiome-based interventions to modulate gut ecology and the immune system

The gut microbiome lies at the intersection between the environment and the host, with the ability to modify host responses to disease-relevant exposures and stimuli. This is evident in how enteric microbes interact with the immune system, e.g., supporting immune maturation in early life, affecting drug efficacy via modulation of immune responses, or influencing development of immune cell populations and their mediators. Many factors modulate gut ecosystem dynamics during daily life and we are just beginning to realise the therapeutic and prophylactic potential of microbiome-based interventions. These approaches vary in application, goal, and mechanisms of action. Some modify the entire community, such as nutritional approaches or faecal microbiota transplantation, while others, such as phage therapy, probiotics, and prebiotics, target specific taxa or strains. In this review, we assessed the experimental evidence for microbiome-based interventions, with a particular focus on their clinical relevance, ecological effects, and modulation of the immune system.

Type 2 immunity in allergic diseases

Significant advancements have been made in understanding the cellular and molecular mechanisms of type 2 immunity in allergic diseases such as asthma, allergic rhinitis, chronic rhinosinusitis, eosinophilic esophagitis (EoE), food and drug allergies, and atopic dermatitis (AD). Type 2 immunity has evolved to protect against parasitic diseases and toxins, plays a role in the expulsion of parasites and larvae from inner tissues to the lumen and outside the body, maintains microbe-rich skin and mucosal epithelial barriers and counterbalances the type 1 immune response and its destructive effects. During the development of a type 2 immune response, an innate immune response initiates starting from epithelial cells and innate lymphoid cells (ILCs), including dendritic cells and macrophages, and translates to adaptive T and B-cell immunity, particularly IgE antibody production. Eosinophils, mast cells and basophils have effects on effector functions. Cytokines from ILC2s and CD4+ helper type 2 (Th2) cells, CD8 + T cells, and NK-T cells, along with myeloid cells, including IL-4, IL-5, IL-9, and IL-13, initiate and sustain allergic inflammation via T cell cells, eosinophils, and ILC2s; promote IgE class switching; and open the epithelial barrier. Epithelial cell activation, alarmin release and barrier dysfunction are key in the development of not only allergic diseases but also many other systemic diseases. Recent biologics targeting the pathways and effector functions of IL4/IL13, IL-5, and IgE have shown promising results for almost all ages, although some patients with severe allergic diseases do not respond to these therapies, highlighting the unmet need for a more detailed and personalized approach.

A Consensus Statement on establishing causality, therapeutic applications and the use of preclinical models in microbiome research

The gut microbiome comprises trillions of microorganisms and profoundly influences human health by modulating metabolism, immune responses and neuronal functions. Disruption in gut microbiome composition is implicated in various inflammatory conditions, metabolic disorders and neurodegenerative diseases. However, determining the underlying mechanisms and establishing cause and effect is extremely difficult. Preclinical models offer crucial insights into the role of the gut microbiome in diseases and help identify potential therapeutic interventions. The Human Microbiome Action Consortium initiated a Delphi survey to assess the utility of preclinical models, including animal and cell-based models, in elucidating the causal role of the gut microbiome in these diseases. The Delphi survey aimed to address the complexity of selecting appropriate preclinical models to investigate disease causality and to study host–microbiome interactions effectively. We adopted a structured approach encompassing a literature review, expert workshops and the Delphi questionnaire to gather insights from a diverse range of stakeholders. Experts were requested to evaluate the strengths, limitations, and suitability of these models in addressing the causal relationship between the gut microbiome and disease pathogenesis. The resulting consensus statements and recommendations provide valuable insights for selecting preclinical models in future studies of gut microbiome-related diseases.

Single cell dynamics and nitrogen transformations in the chain forming diatom Chaetoceros affinis

Colony formation in phytoplankton is often considered a disadvantage during nutrient limitation in aquatic systems. Using stable isotopic tracers combined with secondary ion mass spectrometry (SIMS), we unravel cell-specific activities of a chain-forming diatom and interactions with attached bacteria. The uptake of 13C-bicarbonate and15N-nitrate or 15N-ammonium was studied in Chaetoceros affinis during the stationary growth phase. Low cell-to-cell variance of 13C-bicarbonate and 15N-nitrate assimilation within diatom chains prevailed during the early stationary phase. Up to 5% of freshly assimilated 13C and 15N was detected in attached bacteria within 12 h and supported bacterial C- and N-growth rates up to 0.026 h−1. During the mid-stationary phase, diatom chain-length decreased and 13C and 15N-nitrate assimilation was significantly higher in solitary cells as compared to that in chain cells. During the late stationary phase, nitrate assimilation ceased and ammonium assimilation balanced C fixation. At this stage, we observed highly active cells neighboring inactive cells within the same chain. In N-limited regimes, bacterial remineralization of N and the short diffusion distance between neighbors in chains may support surviving cells. This combination of “microbial gardening” and nutrient transfer within diatom chains represents a strategy which challenges current paradigms of nutrient fluxes in plankton communities.

Responses

Your email address will not be published. Required fields are marked *