A genomic view of trophic and metabolic diversity in clade-specific Lamellodysidea sponge microbiomes

doi:10.21203/rs.2.17204/v1

Download PDF

Research

A genomic view of trophic and metabolic diversity in clade-specific Lamellodysidea sponge microbiomes

https://doi.org/10.21203/rs.2.17204/v1

This work is licensed under a CC BY 4.0 License

Journal Publication

published 23 Jun, 2020

Read the published version in Microbiome →

You are reading this older preprint version

Read the latest preprint version →

Background: Marine sponges and their microbiomes contribute significantly to carbon and nutrient cycling in global reefs, processing and remineralizing dissolved and particulate organic matter. Lamellodysidea herbacea sponges obtain additional energy from abundant photosynthetic Hormoscilla cyanobacterial symbionts, which also produce polybrominated diphenyl ethers (PBDEs) chemically similar to anthropogenic pollutants of environmental concern. Potential contributions of non-Hormoscilla bacteria to Lamellodysidea microbiome metabolism and the synthesis and degradation of additional secondary metabolites are currently unknown.

Results: This study has determined relative abundance, taxonomic novelty, metabolic capacities, and secondary metabolite potential in 21 previously uncharacterized, uncultured Lamellodysidea-associated microbial populations by reconstructing near-complete metagenome-assembled genomes (MAGs) to complement 16S rRNA gene amplicon studies. Microbial community structures aligned with sponge host subgroup phylogeny in 16 samples from four host clades collected from multiple sites in Guam over a three year period, including representatives of Alphaproteobacteria, Gammaproteobacteria, Oligoflexia, and Bacteroidetes as well as Cyanobacteria (Hormoscilla). Unexpectedly, microbiomes from one host clade also included Cyanobacteria from the chemically talented genus Prochloron, a common tunicate symbiont. Two novel Alphaprotobacteria MAGs encoded pathways diagnostic for methylotrophic metabolism as well as Type III secretion systems, and have been provisionally assigned to a new order, designated Candidatus Methylospongiales. MAGs from other taxonomic groups encoded light-driven energy production pathways using not only chlorophyll, but also bacteriochlorophyll and proteorhodopsin. Diverse heterotrophic capabilities favoring aerobic versus anaerobic conditions included pathways for degrading chitin, eukaryotic extracellular matrix polymers, phosphonates, dimethylsulfoniopropionate, trimethylamine, and benzoate. Genetic evidence identified an aerobic catabolic pathway for halogenated aromatics that may enable endogenous PBDEs to be used as a carbon and energy source.

Conclusions: The reconstruction of high quality MAGs from all microbial taxa comprising greater than 0.1% of the sponge microbiome enabled species-specific assignment of unique metabolic features that could not have been predicted from taxonomic data alone. This information will promote more representative models of marine invertebrate microbiome contributions to host bioenergetics, the identification of potential new sponge parasites and pathogens based on conserved metabolic and physiological markers, and a better understanding of biosynthetic and degradative pathways for secondary metabolites and halogenated compounds in sponge-associated microbiota.

General Microbiology

Lamellodysidea

sponge microbiome

cyanosponge

Hormoscilla

PBDE

Methylospongia

Prochloron

The prodigious seawater pumping capabilities of filter-feeding marine sponges contribute significantly to carbon cycling in global reef habitats, through the breakdown and remineralization of both dissolved and particulate forms of organic matter [1–3]. Nutrient recycling is facilitated by host-associated microbial communities comprised of diverse sponge-specific taxa [4], including Gamma- and Alphaproteobacteria, Acidobacteria, Actinobacteria, Chloroflexi, Cyanobacteria, Nitrospirae, Tectomicrobia, candidate phylum Poribacteria, and the archaeal Thaumarchaeota (reviewed in [5, 6]).

Some sponge taxa contain abundant cyanobacterial symbionts, whose photosynthetic activities are essential for host health and growth [7, 8]. Studies measuring potential contributions of photosynthetic symbionts to host energy metabolism have historically used Lamellodysidea herbacea as a prototypical model [8–10]. These studies have suggested that phototrophic cyanobacteria can contribute as much as 80% of the total sponge carbon budget, although significant variation has been observed among individual sponges [11].

The cyanobacterial symbionts of L. herbacea are dominated by a single taxonomic clade known as Hormosilla spongeliae (historically described under the genus names Oscillatoria and Phormidium), found in filamentous bundles called trichomes that can comprise up to 50% of sponge weight [8, 12]. Although Hormoscilla sponge symbionts have never been successfully cultured in the laboratory, two genomes of this bacterial clade have recently been assembled from trichome-enriched metagenomic samples [13]. 16S rRNA gene amplification studies have detected sponge-specific Alpha- and Gammaproteobacteria sequences associated with Lamellodysidea [14], but did not investigate how these non-Hormoscilla symbionts might contribute to metabolic activities and ecological interactions in the sponge microbiome.

Hormoscilla spongesymbionts have attracted particular attention for their ability to produce high concentrations of polybrominated diphenyl ethers (PBDEs), which are chemically similar to anthropogenic pollutants of environmental concern [15–19]. Laboratory studies have suggested that naturally occurring PBDEs might act as antimicrobial agents [20, 21] or as predator feeding deterrents [22], but the biologically relevant activities of these compounds in native marine habitats remain controversial.

PBDE product diversity among Lamellodysidea sponges is strongly correlated with the taxonomy of both sponge hosts and their Hormoscilla symbionts, suggesting potential co-evolutionary processes linking host genotype with microbiome and metabolome diversity [16]. The discovery that not all sponge-associated Hormoscilla strains produce PBDEs [15, 16] offers an opportunity to compare microbial community structure and metabolic capabilities in PBDE-producing versus non-producing sponge hosts. The current study addresses these questions through metagenomic reconstruction of 21 new sponge-associated microbial Metagenome Assembled Genomes (MAGs) from Hormoscilla-dominated Lamellodysidea sponge samples, enabling detailed exploration of their potential contributions to the overall microbial community.

Sponge samples and metagenomic assemblies

L. herbacea sponge samples were collected over a three-year period from four different coastal locations in Guam, including both summer and winter seasons (Supplementary Table 1, Additional file 1). Based on a phylogenetic tree of internal transcribed spacer (ITS) regions (Figure 1), these samples included representatives of all four previously described L. herbacea clades [16]. Clade Ia and Ib samples, whose microbial communities have not previously been described, were obtained at multiple different times and locations. Samples from clades II and III, previously analyzed by 16S rRNA gene amplification [14] are represented by tissue samples taken from different individual sponges collected on the same dates at a single site.

Assemblies for each metagenomic sample were initially performed using all available reads for that sample (Supplementary Table 2, Additional file 1). Combined properties of read coverage, nucleotide composition, and predicted protein matches to GenBank nr revealed a relatively small number of discrete scaffold clusters, mapping primarily to database sequences from Cyanobacteria, Bacteroidetes, Alphaproteobacteria, Gammaproteobacteria, and Oligoflexia (Figure 2). Scaffolds derived from host sponge DNA were present at similar coverage depths to bacterial clusters but were distinguishable by nucleotide composition and predicted matches to reference database sequences.

Metagenomic scaffold bins were used to recruit matching raw read pairs that were subsequently used in targeted sub-assemblies to produce 23 consensus population genomes (Table 1, Supplementary Table 3, Additional file 1). All of the assembled genomes except one (clade Ib: GM102ARS1) met or exceeded MIMAG standards for high or medium quality drafts [23]. CheckM analysis [24] estimated that 16 of the 23 MAGs were more than 90% complete, and 19 contained near full-length 16S rRNA gene sequences. The exceptional quality of these MAGs enabled detailed comparisons with both genomically sequenced relatives and environmental 16S rRNA gene surveys.

Taxonomic classification of assembled genomes

Taxonomic classifications for assembled MAGs were established using concatenated multi-locus phylogenetic trees of conserved proteins and matrices of 16S rRNA gene and average amino acid identity (AAI) scores (Figure 3, Supplementary Figures 1–5, Additional file 1). AAI scores were especially valuable in providing taxonomic classification levels for genomes with missing or incomplete 16S rRNA genes and those with few sequenced relatives, according to quantitative threshold ranges previously established for large numbers of database examples [25, 26].

Cyanobacteria genomes from the genus Hormoscilla were present in all samples, as expected, but an unanticipated genome from genus Prochloron (SP5CPC1; Supplementary Figure 1, Additional file 1) was also assembled in sample SP5 (clade Ia). Prochloron are well known tunicate symbionts [27], but have also been observed microbiologically on the surfaces of both marine sponges and sea cucumbers [28, 29]. No sponge-derived matches to Prochloron 16S rRNA genes or the SP5CPC1 genome were identified at 97% nucleotide identity in GenBank nr (Supplementary Table 4, Additional file 1) or the Sponge Microbiome Project Database [30, 31] (Supplementary Table 5, Additional file 1).

Genomes SP12BCY1, SP5BCY1, and GM202BCY1 were assigned to the Bacteroidetes family Cytophagaceae (Supplementary Figure 2, Additional file 1). Based on AAI scores greater than 98%, these three genomes represent members of the same species. Their closest cultured relative was Ekhidna lutea, a free-living aerobic heterotroph isolated from seawater [32]. 16S rRNA gene identities of 95% and AAI scores of 67% suggest membership in the same genus as E. lutea.

Four Alphaproteobacteria genomes (SP12ARB1, SP5ARB1, GM7ARB1, GM7ARB2) were classified in family Rhodobacteraceae (Supplementary Figure 3, Additional file 1). Their most closely related cultured isolate was Nioella sediminis, a free-living aerobic marine bacterium from a sister genus of the Roseobacteria clade [33, 34]. Alphaproteobacteria genomes SP5ARS3, GM102ARS1, and GM202ARS1 fell within family Rhodospirillaceae. Terrestrial plant-associated Azospirillum brasilense [35] was their closest sequenced relative. Other members of genus Azospirillum have been identified in environmental samples from marine habitats, but not yet cultured or sequenced. An AAI score of 87% suggests that SP5ARS3 should be classified in the genus Azospirillum, but GM102ARS1 and GM202ARS1, with AAI scores of 49–51%, represent a new genus within family Rhodospirillaceae. Alphaproteobacterial genomes SP12AHP2 and SP5AHP1 were classified in the Hyphomonadaceae family, as a sister genus to closest isolate relative Hellea balneolensis, an aerobic, heterotrophic bacterium isolated from surface seawater [36]. SP12AHP2 and SP5AHP1 have identical 16S rRNA genes and AAI scores of 99%, qualifying them as members of the same species in a new genus of family Hyphomonadaceae (Supplementary Figure 2B, Additional file 1).

GM7GCV1 was the only genome sequenced from Gammaproteobacteria order Cellovibrionales (Supplementary Figure 4, Additional file 1). Its closest sequenced relative, Halieaceae LZ–16–2, is an uncharacterized bacterium obtained from a mixed laboratory culture with the saxitoxin-producing dinoflagellate Alexandrium tamarense (NZ_RFLW00000000.1). AAI scores of 50–65% with other members of family Halieaceae suggest GM7GCV1 might be classified as a new genus within this family, but 16S rRNA gene identities of 88–89% imply a more distant relationship, potentially in a new family.

The remaining community genomes were too distant from reference database examples to allow precise taxonomic assignments. GM7GCR1 had no close matches among sequenced isolates, being most closely related to unclassified environmental Gammaproteobacteria with 16S rRNA gene nucleotide sequence identities of 90% and AAI scores below 45% (Supplementary Figure 3 and Supplementary Tables 4–5, Additional file 1). These results suggest it should be classified as either a new family within Chromatiales or a new order within Gammaproteobacteria.

Two Alphaproteobacteria genomes (GM202ARS2, GM7ARS4) fell outside any previously established orders, and could only be assigned at the class level, although their 55% AAI value when compared to each other implies they might be members of a single family (Supplementary Figure 2, Additional file 1). Genome SP5OBV1 was most closely related to Bdellovibrio bacteriovorus from the order Bdellovibrionales, but low 16S rRNA gene identity (82%) and AAI (40%) scores, suggest it might represent a new, previously unreported order within the recently described Proteobacteria class Oligoflexia (Supplementary Figure 5, Additional file 1) [37]. Verification of these taxonomic assignments will require the discovery of additional sequences for closely related taxa.

16S rRNA genes from Rhodobacteriaceae and Rhodospirillaceae genomes were conserved in samples obtained over three years in the current study, and also matched previously reported sequences from a 2005 amplification study of Guam Lamellodysidea sponges [14] at 98–100% identity (Supplementary Table 4, Additional file 1). No 16S rRNA gene matches to the Cytophagaceae, Cellvibrionales, Hyphomonadaceae, and Bdellovibrionales genomes from this study were detected in the 2005 Lamellodysidea microbiome study, but this may be due to the absence of samples from host clades Ia and Ib and the small number of clones analyzed (< 40 per sample) in this earlier work [14]. The closest GenBank 16S rRNA gene matches to MAGs in the current study were associated with Ircinia, Tethya. Stylissa, and Axinella sponges, at 90–96% identity levels (Supplementary Table 4, Additional file 1). Additional low abundance 16S rRNA gene matches were present at 97% identity in numerous other sponge genera from the Sponge Microbiome Project Database (Supplementary Table 5, Additional file 1). No species-level matches were detected in this database for 16S rRNA gene sequences from GM7GVC1 or SP5CPC1, but more distant matches were found at 95% identity, suggesting the presence of bacteria from the same genera but not the same species [38].

Relative abundance comparisons between samples

Microbial community structures in all samples were strongly correlated with host clade taxonomy (Figure 4). The most taxonomically diverse communities were observed in host clade Ia, which included two bacterial taxa that were absent from all other clades: alphaproteobacterial family Hyphomonadaecae and cyanobacterial genus Prochloron. Host clade II samples were notable for the consistent replacement of Alphaproteobacteria family Rhodobacteraceae with Rhodospirillaceae. Sponge metagenome 16S rRNA gene sequences from assembled MAGs accounted for 44–88% of all bacterial amplicon sequence variants (ASVs) in sponge samples, but only 0.1% in nearby seawater (Supplementary Figure 6, Additional file 1).

Relative abundances of the four MAGs with incomplete or missing 16S rRNA sequences (GM102ARS1, GM102CHS1, SP5ARS3, and GM7ARS4) were assessed by recruitment of unassembled metagenomic reads to assembled MAGs (Supplementary Figure 7, Additional file 1). Metagenomic read recruitment is known to underestimate relative abundance of genomes with smaller sizes, but 16S rRNA gene amplification over-reports taxa with multiple gene copies, precluding exact numerical agreement. Despite these differences, results from both procedures were consistent in showing the host clade-specificity of taxonomic compositions, including greater diversity and the unique presence of Prochoron and Hyphomonadaceae in clade 1a, as well as the replacement of Rhodobacteraceae with Rhodospirillaceae in clade II samples. These patterns were consistent over all available sample collection time points and locations tested.

Predicted bacterial lifestyles

Predicted protein annotations for the 23 microbial genomes in this study were used to explore both broad, community-wide patterns and detailed metabolic pathways specific to individual genomes (Figure 5). Shared features among all taxonomic groups included the presence of aerobic respiration, glycolysis, and TCA cycle enzymes. None of the assembled genomes contained genes encoding flagellar biosynthesis, although genes for gliding motility were found in all Cytophagaceae, and twitching motility in Cyanobacteria, GM7GCV1, and SP5OBV1 genomes (Supplementary Table 6, Additional file 1). Several genome groups contained expanded gene families encoding adhesive molecules with the potential to resist shear forces from high seawater flow rates, for example genes encoding cadherin and ankyrin domains, as well as type IV pilus structures. The absence of nitrogenase complex genes combined with the near-universal presence of ammonia transporters suggests ready availability of fixed nitrogen, consistent with detection of ammonia excretion in many sponge species [2].

Metabolic phenotype analysis of the assembled genomes identified pathways associated with phototrophic, methylotrophic, heterotrophic, and parasitic or pathogenic lifestyles. Phototrophic pathways included not only chlorophyll-based photosynthesis with carbon fixation via the Calvin-Benson Cycle in Hormoscilla and Prochloron cyanobacterial genomes, but also bacteriochlorophyll-mediated anoxygenic photosynthesis in alphaproteobacterial Rhodobacteriaceae, coupled with carbon fixation via the reductive TCA cycle. A third mode of light-driven energy production was identified in Bacteroidetes family Cytophagaceae genomes, supporting ATP synthesis through proteorhodopsin-generated proton motive force. Rhodobacteriaceae, Hyphomonadaceae, and Cytophagaceae genomes all contained carotenoid biosynthesis pathways, offering potential protection against free radicals generated during photosynthesis and/or exposure to ultraviolet radiation.

The GM7ARS4, GM202ARS2 and SP5GCR1 genomes all included a complete methylene-tetrahydromethanopterin dehydrogenase (dH4MPT)-dependent oxidation pathway, diagnostic for methylotrophic C–1 metabolism [39]. Each of these genomes also encoded a complete Type III secretion system (T3SS), often associated with eukaryotic pathogenicity. No pathways for dH4MPT-dependent oxidation or T3SS biosynthesis were present in other genomes recovered from this study. Based on these data, the taxonomic relationship between GM7ARS4 and GM202ARS2, and the absence of any other previously described bacteria from the same order, we have provisionally assigned these two genomes to a new order named Candidatus Methylospongiales.

T3SS operons in both SP5GCR1 and the two Methylospongiales genomes contained 15–18 proteins annotated as T3SS components, including base, inner rod, and needle proteins, pore-forming translocation proteins, chaperonins, ATPases and regulatory proteins. Although it was not possible to determine the nature of substrates being transported, T3SS operons also included matches to distinctive virulence-associated protein families (e.g. YscX), associated with human pathogens [40]. Placement of sponge microbiome YscR homologs in a reference tree of conserved, habitat-classified examples (Supplementary Figure 8, Additional file 1) [41], shows their closest relatives derive from extracellular bacteria associated with animal and insect hosts, but not plants, protists, or fungi. These results suggest possible roles for GM7ARS4, GM202ARS2, and SP5GCR1 in sponge host parasitism and/or pathogenicity. Genome sizes ranged from 1.6–1.8 Mbp, consistent with small sizes often observed in obligate symbionts and pathogenic bacteria. [42–44]

Evidence identifying other community members as potential obligate symbionts was constrained by assembled genome incompleteness and limited availability of well-characterized free-living relatives. However, supporting evidence was provided in some cases by consistent presence or absence of clade-specific diagnostic features in multiple closely-related MAGs from different samples. As an example, Hormoscilla genomes GUM202 and GUM007 have previously been shown to lack complete pathways for biotin synthesis [13]. This same deficiency was confirmed in Hormoscilla genomes from samples SP5, SP12, and GUM102, along with the presence of multiple transporters for importing this essential cofactor.

Genomic streamlining has previously been proposed as a shared feature of pelagic marine Rhodbacteriacae of the Roseobacter clade [45]. The genomes of GM7ARB1 and GM7ARB2, estimated to be 94–97% complete, are approximately 2.55 Mbp in size. This is smaller than the previously reported minimum genome size of 3.3 Mbp for Roseobacters, as well as closest free-living relative Nioella nitratereducens (4.0 Mbp). The four Rhodbacteriacae in the current study have 8–10% of their genomes devoted to transporter functions (Supplementary Table 6, Additional file 1), but lack genes for flagellar biosynthesis, and RuBisCO-mediated carbon fixation typically found in other family members. However, many other clade-specific characteristics have been preserved, including pathways for the synthesis of capsular polysaccharides and the degradation of phosphonates, urate and dimethylsulfoniopropionate (DSMP).

Marine Bacteroidetes from the Cytophaga-Flavobacteria group are noted for their role in degrading organic matter during phytoplankton blooms, but group members containing proteorhodopsin have consistently smaller genomes than close taxonomic relatives lacking this gene function [46]. Consistent with these observations, Cytophagaceae genomes SP12BCY2, SP5BCY1 and GM202BCY1 (~3 Mbp) are similar in size to free-living family members containing proteorhodopsin, but smaller than closest sequenced relative Ekhidna lutea (4.2 Mbp), which does not. SP12BCY2, SP5BCY1 and GM202BCY1 are enriched in signal peptide-containing peptidases (9–12 per genome) and glycosidases (3–4 per genome), suggesting the retention of conserved heterotrophic capabilities for degrading extracellular proteins and polysaccharides (Supplementary Table 6, Additional file 1).

All three of these genomes have retained Bacteroidetes-specific gliding motility and Por (Type IX) secretion system functions. The Type IX secretion system is both an essential component of gliding motility in non-pathogenic species [47, 48] and a virulence factor in human and fish pathogens [49–51]. A large number of novel cadherin domain-containing proteins (Supplementary Table 6, Additional file 1), with closest GenBank matches at 29–44% amino acid identity, may facilitate adhesion to sponge hosts and/or particulate matter. [52]

Potential symbiotic adaptations in Hyphomonadaceae genomes SP12HP1 and SP5AHP2 were suggested by the absence of both flagellar motility genes and the stalk formation pathway found in many other Hyphomonadaceae genomes [53]. However, genome sizes (2.4–2.5 Mbp) were only about 25% smaller than closest free-living relative Hellea balneolensis (3.2 Mbp). Lineage characteristic genes for heparinases, chondroitinases and chitinases, and signal-peptide containing glycosidases were retained in SP12HP1 and SP5AHP2, potentially facilitating the degradation of host extracellular matrix as well as organic marine particulates. The presence of typical Hyphomonadaceae pathways for the transport and degradation of aromatic compounds such as benzoic acid suggests retention of diverse heterotrophic capabilities.

Rhodospirillales genomes SP5ARS3, GUM102ARS1, and GM202ARS1 have much smaller genomes than their closest sequenced relative, the terrestrial soil bacterium Azospirillum brasilense (<3 versus 7.2 Mbp), but this difference could be due to expansion in Azospirillum rather than reduction in sponge-associated genomes. SP5ARS3, GUM102ARS1, and GM202ARS1 are unique among the genomes of this study in encoding pathways for carbon fixation via the Wood-Ljungdahl (Reductive Acetyl-CoA) pathway, the glyoxylate cycle, and trimethylamine degradation. The presence of nitrile hydratases and a large repertoire of transporters (9–12% of predicted coding sequences; Supplementary Table 6, Additional file 1) with predicted specificities for amino acids, oligopeptides, taurine, spermidine/putrescein, lipoproteins, ribose/xylose, and glycerol suggest metabolic versatility encompassing a wide variety of substrates, especially those containing amine groups. Genes encoding capsular biosynthesis may provide protection against host phagocytosis, viral attack, and/or hydrophobic toxins.

Marine relatives of genome GM7GCV1 include both free-living and host-associated species from Gammaproteobacteria family Cellovibronaceae, with AAI scores of 50% and 16S rRNA gene identities of 87–88%. All cultured examples have genome sizes of 4 Mbp or larger, including the dinoflagellate-associated Halieaceae LZ–16–2, GM7GCV1’s closest relative. The much smaller genome size of GM7GCV (1.8 Mbp) is typical for more distantly related, marine particle-associated Cellovibronaceae from the Tara Oceans project (e.g. TMED119 at AAI 45%, Supplementary Figure 3, Additional file 1), whose genome sizes range from 1.1 - 2.9 Mbp [54]. Factors responsible for smaller genome sizes in these uncultured, uncharacterized bacterial relatives are unknown.

The unavailability of sequenced genomes for relatives closer than order level makes it difficult to determine whether the SP5OBV1 genome lacks clade-specific features. However, it is unlikely that this species follows the same bacterivorous lifestyle as its distant cousins, not only because its genome size, estimated at 91% complete, is less than half as large as closest relative Bdellovibrio bacteriovorus (1.7 versus 3.8 Mbp), but also due to the absence of genes encoding flagellar motility and proteoglycan synthesis described as essential for this activity in other Bdellovibrionales species (reviewed in [55]). Fatty acid auxotrophy in SP5OBV1 is suggested by the absence of acetyl-CoA carboxylase and all other fatty acid biosynthesis pathway enzymes except FabG, coupled with the presence of complete pathways for lipoic acid metabolism and fatty acid degradation. Some mammalian pathogens are known to suppress endogenous biosynthesis while incorporating exogenous host fatty acids as a triclosan resistance mechanism, but previously reported genomic pathway deficiencies accompanying these adaptations have been much less extensive than those observed in SP5OVB1 [56].

Viral sequences and phage defense

Electron microscopic studies have identified phage-like particles associated with Lamellodysidea sponge samples [57], but these have not yet been characterized by sequence analysis. In this study, viral scaffold candidates detected by VirSorter [58] in assembled holobiont metagenomes were dominated by GenBank matches to double stranded DNA tailed bacteriophage (Caudovirales) of the Siphoviridae, Myoviridae, and Podoviridae lineages (Supplementary Figure 9, Additional file 1). Matches to Microviridae, archaeal Bicaudaviridae, and eukaryotic Baculoviridae and Herpes viruses were also detected at low levels. Genetic heterogeneity of integrated phage regions can make them difficult to capture in metagenomic assemblies, potentially underestimating occurrence. However, abundant CRISPR sequences, transposons, and restriction enzymes roughly proportional to genome size in assembled MAGs suggest a past history of recurring viral challenges (Supplementary Table 7, Additional file 1).

Integrated prophage genomes sometimes carry passenger genes of bacterial origin that can modify the phenotype of the host, resulting in improved fitness of infected cells, a process known as lysogenic conversion [59]. Lamellodysidea-associated phage candidate scaffolds included bacterial genes encoding DNA methyltransferases, glucanases, peptidases, and vitamin B–12 biosynthesis. One candidate Myoviridae scaffold, from sample SP5, encoded two key enzymes (PhnI and PhnJ) of the carbon-phosphonate lyase pathway [60]. This pathway has previously been demonstrated to be enriched under phosphate limiting conditions in metagenomic marine microbes [61], and could potentially assist in liberating phosphate from recalcitrant organic particles.

Secondary metabolite pathways

Secondary metabolite gene cluster candidates identified by antiSMASH [62] were most abundant in cyanobacterial taxa, with the largest number in previously described Hormoscilla genomes GM7CHS1 and GM202CHS1 (Supplementary Table 8, Additional file 1; [13]). Hormoscilla genomes SP12CHS1, SP5CHS1, and GM102CHS1 had fewer predicted clusters, but antiSMASH cluster detection sensitivity may have been reduced by scaffold fragmentation in these MAGs, which were assembled from Illumina reads without PacBio read supplementation.

The next most abundant source of biosynthetic gene clusters was Prochloron genome SP5CPC1, predicted to include two non-ribosomal peptide synthetases (NRPS), four terpene synthases, three ribosomally synthesized post-translationally modified peptides (RiPPs), and two flavin-dependent aromatic halogenases (Supplementary Table 8, Additional file 1). Although SP5CPC1 has a smaller genome size with fewer total clusters than previously sequenced Prochloron relatives [63], SP5CPC1 is also less complete, with short scaffolds and fragmented, incomplete pathway sequences that cannot easily be linked to specific molecular products. The aromatic halogenases were unrelated to those found in PBDE-producing strains of Hormoscilla spongeliae ([16], Supplementary Table 10B, Additional file 1), and putative RiPP clusters could only be classified as distantly related to non-cyanobactin bacteriocins and lassopeptides, based on HMM (Hidden Markov Model) pattern matches [62]. None of the putative RiPP clusters contained genes from the well-characterized patellamide pathway of Prochloron didimei [64].

Predicted biosynthetic clusters in non-cyanobacterial genomes were limited to terpenes and bacteriocin-related RiPPs, except for the SP12BCY1, SP5BCY1, and GUM202BCY1 genomes, which all contained an identical type III polyketide synthase (PKS). The closest database matches to this protein were from Cytophagaceae genera such as Ekhidna, Marinoscillum, and Pontibacter at 53–64% amino acid identity. Biosynthesis of plant-like flavonoids like those recently characterized in Flavibacteria cheonhonense [65] seems unlikely, because key pathway enzyme phenylalanine ammonia lyase is missing [66, 67]. However, type III PKS genes from the current study also matched CepAB genes from taxonomically distant sponge symbiont Entotheonella gemina (class Tectomicrobia) at 48% amino acid identity, suggesting phenolic lipids as a potential biosynthetic product [68].

Previously reported gene clusters encoding polybrominated compounds in assembled Hormoscilla genomes [13, 16] do not account for the full range of halogenated products described in Lamellodysidea sponges [69]. Additional diversity may be contributed by members of the broader microbial community, for example Hyphomonadaceae SP12AHP1 and SP5AHP2, which encode a highly expanded family of novel aromatic flavin-dependent halogenases (Supplementary Figure 10AB, Additional file 1). Several of these halogenases occur in sets of 2–4 tandem repeats, suggesting family expansion by gene duplication. Based on similarity to Pfam database model PF04820 [70], these proteins are annotated as flavin-dependent tryptophan halogenases. However, their closest match among experimentally characterized enzymes is the brvH gene product from Brevundimonas sp. BAL3, at 47% amino acid identity (Supplementary Figure 10B, Additional file 1). BrvH uses free indole rather than tryptophan as a substrate, preferentially incorporating bromine over chlorine into the C3 position [71],

It is possible that some of the additional diversity in halogenated compounds previously observed in Lamellodysidea sponges, including brominated phenols and catechols [72, 73], arises from degradative, rather than biosynthetic pathways. The SP12AHP1 and SP5AHP2 genomes each encode a complete 3,5-dichlorocatechol degradation pathway, characteristic of bacteria that use halogenated benzoates as sole carbon and energy sources [74]. Non-halogenated benzoate degradation pathways are common in other Hyphomonadaceae, but the SP12AHP1 and SP5AHP2 genomes are unique in containing chlorocatechol 1,2-dioxygenase, a key enzyme that cleaves chlorocatechol rings to 2,4-dichloro-cis,cis-muconate, followed by spontaneous dehalogenation during further processing by muconate cycloisomerase, carboxymethylenebutenolidase, and maleylacetate reductase [75]. The closest chlorocatechol 1,2-dioxygenase match in GenBank nr (64% amino acid identity) was found in Altererythrobacter marensis, a free-living coastal marine bacterium from the Alphaproteobacteria order Sphingomonadales [76]. Dichlorocatechol degradation pathways in the SP12AHP1 and SP5AHP2 genomes are clustered together in the same conserved gene order as A. marensis (Supplementary Figure 11, Additional file 1).

Plasmids containing the 3,5-dichlorocatechol degradation pathway are frequently exchanged between environmental bacteria in chloroaromatic contaminated environments (reviewed in [77]). Although the pathway is not located on a plasmid or within a genomic island in SP12AHP1, SP5APH2 or A. marensis, patchy phylogenetic distribution and order-level relationships between closest protein sequence relatives suggests historical dissemination by horizontal gene transfer (Supplementary Figure 12, Additional file 1). It is not known whether amino acid sequence similarity to experimentally characterized chlorocatechol dioxygenases might also capture activity towards brominated substrates, including PBDEs, but this seems like a reasonable hypothesis for future testing.

Sponge-associated microbial communities have historically been classified into two broad categories, designated HMA (high microbial abundance, high diversity) and LMA (low microbial abundance, low diversity) [78]. High sponge seawater pumping rates have been proposed to promote aerobic, LMA communities, with lower rates favoring more anaerobic HMA communities [2]. These hypotheses are supported by results in the current study, where multiple samples from Lamellodysidea herbacea sponges, noted for high seawater pumping rates [3], contain a predominance of bacteria favoring aerobic rather than anaerobic metabolic pathways, with low overall community diversity.

The Lamellodysidea herbacea sponge microbiome has previously been characterized almost exclusively in terms of functional activities associated with its most dominant symbiont, Hormoscilla spongeliae. The results of this study reveal previously unknown metabolic activities, interactions, and potential community contributions of non-Hormoscilla species. Low microbial community complexity facilitated effective taxonomic binning and metagenomic assembly, producing long scaffolds enabling operon-context assessment of functional gene activities. The assembly of multiple population genomes from closely related taxa in different, independent samples helped establish core functional characteristics for novel sponge-associated groups previously known only from 16S rRNA gene fragments, and encountered only rarely in other sponge species.

The GM7ARS4 and GM202ARS2 genomes, classified in a new Alphaproteobacteria order (Candidatus Methylospongiales), and SP5GCR1, representing a new Gammaproteobacteria family, are, to the best of our knowledge, the first reported candidate pathogens of Lamellodysidea herbacea. The unlikely combination of methylotrophy with type III secretion systems consistently shared in distantly related taxa suggests potential exploitation of a common niche opportunity. Although no close relatives of these organisms have yet been reported in other sponges, the three near-complete genomes obtained in this study should enable more effective future searches based on predicted metabolic and physiological features.

The genomes sequenced in this study show strong evidence of adaptation to their sponge-host environment. Abundant light availability in the reef habitat fuels energy generation through three different processes, using phototrophic pigments with widely varying absorption maxima: chlorophyll a (660 nm), bacteriochlorphyll a (800 nm), and proteorhodopsin (~520 nm). UV damage protection is provided by carotenoid pigments. The possibility that photoactive pigments with differing optimal absorption spectra might provide advantages within micro-niches associated with different layers of sponge tissue would be an interesting topic for future studies.

High seawater flow rates induced by host pumping action create strong shear forces that could potentially dislodge microbiome inhabitants from sponge tissues. The expanded families of adhesive molecules and type IV pili found in many Lamellodysidea-associated bacteria would be expected to counteract these disruptive forces. Oxygen depleted by microbial respiration should be rapidly replenished from seawater flow and cyanobacterial photosynthesis, consistent with an environment that appears to favor organisms using aerobic processes for heterotrophic degradation of dissolved and particulate organics, and potentially discriminating against symbionts requiring more anaerobic conditions. Nitrogenase enzymes requiring the exclusion of oxygen were not observed. They may however be unnecessary in this environment due to a combination of host ammonia excretion and the release of fixed nitrogen from the degradation of organic compounds [2].

Previous bioenergetic analyses of microbial contributions to Lamellodysidea sponge nutrition have been limited to relatively simplistic models based on chlorophyll-mediated photosynthesis by a single microorganism. Functional pathway predictions and abundance data describing the new phototrophic and heterotrophic symbionts identified in this study should allow the expansion of these models to include contributions from multiple metabolic sources, encompassing a more accurate sampling of microbiome community diversity.

Microbial community structures in all Lamellodysidea samples were strongly correlated with host clade taxonomy, whether assessed by metagenomic read recruitment or 16S rRNA gene analysis. The persistent co-occurrence of highly conserved Rhodobacteriaceae and Rhodospirillaceae species with Hormoscilla over an eleven year time period in Guam Lamellodysidea sponges (from 2005 in [14] to 2016 in the current work) suggests the likelihood of relatively stable, long-term host-microbe relationships. Previous analyses of Guam Lamellodysidea sponge microbiomes [14] failed to identify many of the taxa described in the current study, but examined 100-fold fewer 16S rRNA gene sequences, and did not include any representatives of host clades Ia or Ib. Further studies will be required to understand why host clade Ia contains some taxa, especially Hyphomonadaceae and Prochloron, that are absent in other clades, as well as why photoheterotrophic Rhodobacteriaceae were relatively abundant in clades Ia, Ib, and III, but nearly absent from clade II, where they were replaced by non-photosynthetic Rhodospirillaceae.

The biological functions and effects of locally concentrated PBDEs produced by Hormoscilla bacteria are unknown. Genomes SP12AHP and SP5AHP2 from clade Ia PBDE-producing sponges encode a potential mechanism for aerobic degradation of these compounds that may enable their use as nutritional energy sources. This proposed pathway includes a modified dioxygenase ring cleavage mechanism related to those used by other Hyphomonaceae for breaking down non-halogenated aromatic hydrocarbons.

A large number of biologically active secondary metabolite compounds have been isolated from Lamellodysidea herbacea sponges (reviewed in [79]), but distinguishing between those synthesized by the eukaryotic host versus microbial symbionts has proved challenging, especially in the absence of an assembled host genome. The unexpected discovery that the microbiomes of Lamellodysida herbacea clade Ia sponges contain representatives of the chemically talented genus Prochloron, in addition to Hormoscilla, provides new opportunities to identify microbial sources for both previously characterized and undiscovered natural product molecules.

The reconstruction of high quality MAGs representing all microbial taxa at relative abundances greater than 0.1% has enabled species-specific assignment of unique metabolic features that could not have been predicted from taxonomic data alone. This information will promote more representative models of marine invertebrate microbiome contributions to host bioenergetics, identification of new sponge pathogens based on conserved metabolic and physiological marker genes, and a better understanding of biosynthetic and degradative pathways for secondary metabolites and halogenated compounds in sponge-associated microbiota.

Sample collection and processing

Sponge samples used in this study, shown in Supplementary Table 1, Additional file 1, were collected between 2014-2016 from four different coastal sites in Guam. Frozen tissue samples were processed and prepared for DNA sequencing using Illumina HiSeq paired-end technology as described previously [13, 16]. Additional PacBio sequencing reads, obtained from cyanobacterial trichome enrichments of samples GUM7 and GUM202 [13], were used in the assembly of Hormoscilla spongeliae genomes from these samples. Sample host taxonomy was verified by amplification of ITS primer regions from bulk DNA, according to the protocol described in [80].

16S rRNA gene amplification was performed using V4-5 region primers 515F‐4Y (5′‐GTGYCAGCMGCCGCGGTAA) and 926R (5′‐CCGYCAATTYMTTTRAGTTT), as previously described [81]. Paired end reads were generating using Illumina MiSeq v3 (2 x 300 bp) sequencing. DADA2 v 1.12 [82] was used within QIIME2 [83] to trim, denoise, merge, remove chimeras, and cluster the assembled reads, with the following program parameters: --p-trunc-len-f 250 --p-trunc-len-r 180 --p-trim-left-f 20 --p-trim-left-r 20 --p-max-ee-f 3 --p-max-ee-r 5. Taxonomic classifications were assigned based on the SILVA v132 reference database [84], supplemented with 19 near-complete 16S rRNA genes obtained from targeted metagenomic assemblies in this study.

Metagenomic assembly

Illumina HiSeq 2500 reads (2 x 100 or 2 x150 PE) from metagenomic DNA libraries were quality filtered and trimmed using Trimmomatic version 0.339 [85]. Preliminary guide assemblies were created using IDBA-UD version 1.1.117 [86] with default parameters, followed by mapping of input reads to scaffolds using the end-to-end option of Bowtie2, version 2.218 [87]. Coverage depth and percentage of singleton versus paired reads were calculated from output of the idxstats module from samtools version 0.1.191 [88]. Preliminary scaffold bins were assigned based on percent GC, nucleotide composition, coverage depth, and taxonomic classification by DarkHorse version 2.0 [89, 90], as previously described [91]. Read subsets were supplemented as necessary to include both R1 and R2 members of each constituent read pair, then re-assembled.

Second round assemblies were performed using both IDBA-UD, and, alternatively Celera Assembler version 8.3 [92], configured with merSize =17, utgGenomeSize = 5Mb and utgErrorRate = 0.01. Second round assembly scaffolds were quality-screened by re-applying the same methods used for initial taxonomic binning (percent GC, coverage depth, and DarkHorse taxonomic classification), with an added requirement that each scaffold should contain more than 90% paired (versus singleton) reads. Final assembly completeness and contamination were assessed using CheckM version 1.07 with the default set of bacterial marker genes [24].

Hormoscilla genome assemblies have been previously reported for samples GUM202 and GUM007 using a combination of Illumina and PacBio reads from trichome enriched preparations [13]. In the current study, PacBio reads from Sample GUM007 were re-assembled without Illumina read supplementation using CANU version 1.8 with default parameters [93], which greatly improved the GUM007 Hormoscilla assembly (100% versus 94.5% completeness, 1.8% versus 2.2% contamination, and 6 versus 64 contigs). This approach could not, however, be applied to the GUM202 Hormoscilla genome, because the smaller number of available PacBio reads (31,034 versus 180,838) failed to provide sufficient coverage.

Phylogenetic placement

16S rRNA gene sequences for closely related species were obtained by blastn searches against SILVA database release 132 [94], and aligned with sequences extracted from assembled genomes using the SILVA Incremental Aligner (SINA) version 1.2.11 [84]. 16S rRNA gene trees were constructed from these alignments using FastTree version 2.1.8 [95] and visualized with FigTree version 1.4.3 [96]. Multi-locus gene trees were constructed using PhyloPhlAn version 0.99 [97]. Average amino acid identities were calculated using the online AAI-Matrix Genome-based distance matrix calculator [26].

Environmental database searches

16S rRNA genes identified in metagenome-assembled genomes (MAGs) from this study were compared to reference database sets by blastn searches, with a minimum cutoff alignment length of 360 nucleotides and maximum e-value of 1e-7. Quantitative abundances of blastn matches to OTUs from samples in the Sponge Microbiome Database were obtained by querying a custom MySQL database, constructed by reformatting Supplementary files S3, S4, and S5 from [30] into the schema shown in Supplementary Figure 13, Additional file 1.

Functional Annotation

All MAGs generated in this study were annotated at IMG-MER [98]. CRISPR repeat regions were identified using the MinCED (Mining CRISPRs in Environmental Datasets) program, version 0.3.0 [99]. Potentially over-represented protein functional families were identified using Hidden Markov Models (HMMs) from the Pfam-A version 32 [70] and TIGRFAM release 15 [100] databases. In cases where models with overlapping functional activities matched the same target protein, only the HMM with the highest bitscore was included in quantitative tallies, so that no protein was counted more than once. Candidate viral sequences were identified using VirSorter version 1.0.5 [58], and taxonomically classified using DarkHorse version 2.0 with a custom reference database consisting of all viral sequences in Genbank nr as of January 2, 2019. Biosynthetic clusters encoding potential secondary metabolites were identified using AntiSMASH version 5.0 [62].

PBDE polybrominated diphenyl ether; ITS internal transcribed spacer; Mbp, megabase pair; nt, nucleotide; MAG, metagenome associated genome; TCA tricarboxylic acid; NRPS non-ribosomal peptide synthase; HMM hidden markov model; PKS, polyketide synthase.

Ethics approval and consent to participate

Not applicable

Consent for publication

Not applicable

Availability of data and material

Sequences datasets generated and/or analyzed during the current study (sequencing reads, genome assemblies, and amplified ITS and 16S rRNA sequences) are available in the NCBI Genbank repository under BioProject ID PRJNA320446, and the IMG-MER database [101] under genome accession numbers 2643221529, 2747842407, 2758568580, 2758568581, 2767802581, 2808606420-34, 2814122959, 2830827416, and 3300008641.

Competing interests

The authors declare no conflicts of interest.

Funding

This work was supported by grants from NSF (OCE-1837116) and NIEHS (R01-ES030316 to BSM and EEA, and R00-ES026620 to V.A.).

Authors' contributions

S.P., V.A., B.S.M. and E.E.A. designed the study. S.P., J.M.B., A.O. and M.A.S. performed experiments. J.B. provided sponge samples, analytical tools and reagents. S.P. and E.E.A wrote the manuscript with input from all authors.

Acknowledgements

We thank our University of California, San Diego colleagues Dr. Alvaro Muñoz Plominsky for critical reading and suggestions for improvement of the manuscript and Katherine Bauman for assistance in sponge collections and DNA isolation.

1.Weisz, J. B., N. Lindquist, and C. S. Martens, Do associated microbial abundances impact marine demosponge pumping rates and tissue densities? Oecologia, 2008. 155(2): p. 367–76.

2.Maldonado, M., M. Ribes, and F. C. van Duyl, Nutrient Fluxes through Sponges: Biology, Budgets, and Ecological Implications. Advances in Sponge Science: Physiology, Chemical and Microbial Diversity, Biotechnology, 2012. 62: p. 113–182.

3.Biggerstaff, A., et al., Metabolic responses of a phototrophic sponge to sedimentation supports transitions to sponge-dominated reefs. Sci Rep, 2017. 7(1): p. 2725.

4.Taylor, M. W., et al., Soaking it up: the complex lives of marine sponges and their microbial associates. ISME J, 2007. 1(3): p. 187–90.

5.Pita, L., et al., The sponge holobiont in a changing ocean: from microbes to ecosystems. Microbiome, 2018 Mar 9. 6(1): p. Article number: 46.

6.Webster, N. S. and T. Thomas, The Sponge Hologenome. MBio, 2016. 7(2): p. e00135–16.

7.Bercerro MA, P. V., Effects of depth and light on secondary metabolites and cyanobacterial symbionts of the sponge Dysidea granulosa. Mar Ecol Prog Ser, 2004. 280: p. 115–128.

8.Thacker, R. W., Impacts of shading on sponge-Cyanobacteria symbioses: A comparison between host-specific and generalist associations. Integrative and Comparative Biology, 2005. 45(2): p. 369–376.

9.Hinde, R., F. Pironet, and M. A. Borowitzka, Isolation of Oscillatoria-Spongeliae, the Filamentous Cyanobacterial Symbiont of the Marine Sponge Dysidea-Herbacea. Marine Biology, 1994. 119(1): p. 99–104.

10.Usher, K. M., The ecology and phylogeny of cyanobacterial symbionts in sponges. Marine Ecology-an Evolutionary Perspective, 2008. 29(2): p. 178–192.

11.Freeman, C. J., et al., Quality or quantity: is nutrient transfer driven more by symbiont identity and productivity than by symbiont abundance? ISME J, 2013. 7(6): p. 1116–25.

12.Berthold, R. J., M. A. Borowitzka, and M. A. Mackay, The Ultrastructure of Oscillatoria-Spongeliae, the Blue-Green-Algal Endosymbiont of the Sponge Dysidea-Herbacea. Phycologia, 1982. 21(3): p. 327–335.

13.Schorn, M. A., et al., Comparative Genomics of Cyanobacterial Symbionts Reveals Distinct, Specialized Metabolism in Tropical Dysideidae Sponges. MBio, 2019. 10(3): p. e00821–19.

14.Ridley, C. P., D. John Faulkner, and M. G. Haygood, Investigation of Oscillatoria spongeliae-dominated bacterial communities in four dictyoceratid sponges. Appl Environ Microbiol, 2005. 71(11): p. 7366–75.

15.Ridley, C. P., et al., Speciation and biosynthetic variation in four dictyoceratid sponges and their cyanobacterial symbiont, Oscillatoria spongeliae. Chem Biol, 2005. 12(3): p. 397–406.

16.Agarwal, V., et al., Metagenomic discovery of polybrominated diphenyl ether biosynthesis by marine sponges. Nat Chem Biol, 2017. 13(5): p. 537–543.

17.Unson, M. D., N. D. Holland, and D. J. Faulkner, A Brominated Secondary Metabolite Synthesized by the Cyanobacterial Symbiont of a Marine Sponge and Accumulation of the Crystalline Metabolite in the Sponge Tissue. Marine Biology, 1994. 119(1): p. 1–11.

18.Hanif, N., et al., Polybrominated diphenyl ethers from the Indonesian sponge Lamellodysidea herbacea. J Nat Prod, 2007. 70(3): p. 432–5.

19.Liu, H., et al., Polybrominated Diphenyl Ethers: Structure Determination and Trends in Antibacterial Activity. J Nat Prod, 2016. 79(7): p. 1872–6.

20.Sharma, G. M. and B. Vig, Studies on.Antimicrobial Substances of Sponges. VI. Structures of 2 Antibacterial Substances Isolated from Marine Sponge Dysidea-Herbacea. Tetrahedron Letters, 1972(17): p. 1715–1718.

21.Sun, S., et al., Polybrominated diphenyl ethers with potent and broad spectrum antimicrobial activity from the marine sponge Dysidea. Bioorg Med Chem Lett, 2015. 25(10): p. 2181–3.

22.Becerro, M. A., et al., Biogeography of sponge chemical ecology: comparisons of tropical and temperate defenses. Oecologia, 2003. 135(1): p. 91–101.

23.Bowers, R. M., et al., Minimum information about a single amplified genome (MISAG) and a metagenome-assembled genome (MIMAG) of bacteria and archaea. Nat Biotechnol, 2017. 35(8): p. 725–731.

24.Parks, D. H., et al., CheckM: assessing the quality of microbial genomes recovered from isolates, single cells, and metagenomes. Genome Res, 2015. 25(7): p. 1043–55.

25.Rodriguez-R, L. M., and Konstantinos T. Konstantinidis, Bypassing Cultivation To Identify Bacterial Species. Microbe, 2014. 9(3): p. 111–8.

26.Rodriguez-R, L. and K. Konstantinidis, The enveomics collection: a toolbox for specialized analyses of microbial genomes and metagenomes. PeerJ Preprints, 2016. 4: p. e1900v1.

27.Lewin, R. A. and L. Cheng, Prochloron: a microbial enigma. 1989, New York: Chapman and Hall. xiii, 129 p., 26 p. of plates.

28.Cheng, L. and R. A. Lewin, Prochloron on Synaptula. Bulletin of Marine Science, 1984. 35(1): p. 95–98.

29.Parry, D. L., Prochloron on the Sponge Aplysilla Sp. Bulletin of Marine Science, 1986. 38(2): p. 388–390.

30.Thomas, T., et al., Diversity, structure and convergent evolution of the global sponge microbiome. Nat Commun, 2016. 7: p. 11870.

31.Moitinho-Silva, L., et al., The sponge microbiome project. Gigascience, 2017. 6(10): p. 1–7.

32.Alain, K., et al., Ekhidna lutea gen. nov., sp. nov., a member of the phylum Bacteroidetes isolated from the South East Pacific Ocean. Int J Syst Evol Microbiol, 2010. 60(Pt 12): p. 2972–8.

33.Liu, Y., et al., Nioella sediminis sp. nov., isolated from surface sediment and emended description of the genus Nioella. Int J Syst Evol Microbiol, 2017. 67(5): p. 1271–1274.

34.Rajasabapathy, R., et al., Nioella nitratireducens gen. nov., sp. nov., a novel member of the family Rhodobacteraceae isolated from Azorean Island. Antonie Van Leeuwenhoek, 2015. 107(2): p. 589–95.

35.Wisniewski-Dye, F., et al., Azospirillum genomes reveal transition of bacteria from aquatic to terrestrial environments. PLoS Genet, 2011. 7(12): p. e1002430.

36.Alain, K., et al., Hellea balneolensis gen. nov., sp. nov., a prosthecate alphaproteobacterium from the Mediterranean Sea. Int J Syst Evol Microbiol, 2008. 58(Pt 11): p. 2511–9.

37.Hahn, M. W., et al., Silvanigrella aquatica gen. nov., sp. nov., isolated from a freshwater lake, description of Silvanigrellaceae fam. nov. and Silvanigrellales ord. nov., reclassification of the order Bdellovibrionales in the class Oligoflexia, reclassification of the families Bacteriovoracaceae and Halobacteriovoraceae in the new order Bacteriovoracales ord. nov., and reclassification of the family Pseudobacteriovoracaceae in the order Oligoflexales. Int J Syst Evol Microbiol, 2017. 67(8): p. 2555–2568.

38.Yarza, P., et al., Uniting the classification of cultured and uncultured bacteria and archaea using 16S rRNA gene sequences. Nat Rev Microbiol, 2014. 12(9): p. 635–45.

39.Martinez-Gomez, N. C., S. Nguyen, and M. E. Lidstrom, Elucidation of the role of the methylene-tetrahydromethanopterin dehydrogenase MtdA in the tetrahydromethanopterin-dependent oxidation pathway in Methylobacterium extorquens AM1. J Bacteriol, 2013. 195(10): p. 2359–67.

40.Gurung, J. M., et al., Heterologous Complementation Studies With the YscX and YscY Protein Families Reveals a Specificity for Yersinia pseudotuberculosis Type III Secretion. Front Cell Infect Microbiol, 2018. 8: p. 80.

41.Nazir, R., et al., The Ecological Role of Type Three Secretion Systems in the Interaction of Bacteria with Fungi in Soil and Related Habitats Is Diverse and Context-Dependent. Front Microbiol, 2017. 8: p. 38.

42.Klasson, L. and S. G. Andersson, Evolution of minimal-gene-sets in host-dependent bacteria. Trends Microbiol, 2004. 12(1): p. 37–43.

43.Moran, N. A., J. P. McCutcheon, and A. Nakabachi, Genomics and evolution of heritable bacterial symbionts. Annu Rev Genet, 2008. 42: p. 165–90.

44.Moya, A., et al., Learning how to live together: genomic insights into prokaryote-animal symbioses. Nat Rev Genet, 2008. 9(3): p. 218–29.

45.Voget, S., et al., Adaptation of an abundant Roseobacter RCA organism to pelagic systems revealed by genomic and transcriptomic analyses. ISME J, 2015. 9(2): p. 371–84.

46.Gonzalez, J. M., et al., Genomics of the proteorhodopsin-containing marine flavobacterium Dokdonia sp. strain MED134. Appl Environ Microbiol, 2011. 77(24): p. 8676–86.

47.McBride, M. J. and D. Nakane, Flavobacterium gliding motility and the type IX secretion system. Curr Opin Microbiol, 2015. 28: p. 72–7.

48.Johnston, J. J., A. Shrivastava, and M. J. McBride, Untangling Flavobacterium johnsoniae Gliding Motility and Protein Secretion. J Bacteriol, 2018. 200(2).

49.Li, N., et al., The Type IX Secretion System Is Required for Virulence of the Fish Pathogen Flavobacterium columnare. Appl Environ Microbiol, 2017. 83(23).

50.Veith, P. D., et al., Type IX secretion: the generation of bacterial cell surface coatings involved in virulence, gliding motility and the degradation of complex biopolymers. Mol Microbiol, 2017. 106(1): p. 35–53.

51.Lasica, A.M., et al., The Type IX Secretion System (T9SS): Highlights and Recent Insights into Its Structure and Function. Front Cell Infect Microbiol, 2017. 7: p. 215.

52.Fraiberg, M., et al., Discovery and Characterization of Cadherin Domains in Saccharophagus degradans 2–40. Journal of Bacteriology, 2010. 192(4): p. 1066–1074.

53.Abraham, W. and M. Rohde, The Family Hyphomonadaceae, in The Prokaryotes—Alphaproteobacteria and Betaproteobacteria, Rosenberg E, et al., Editors. 2014, Springer: Berlin.

54.Tully, B. J., E. D. Graham, and J. F. Heidelberg, The reconstruction of 2,631 draft metagenome-assembled genomes from the global oceans. Sci Data, 2018. 5: p. 170203.

55.Sockett, R. E., Predatory lifestyle of Bdellovibrio bacteriovorus. Annu Rev Microbiol, 2009. 63: p. 523–39.

56.Yao, J. and C. O. Rock, How bacterial pathogens eat host lipids: implications for the development of fatty acid synthesis therapeutics. J Biol Chem, 2015. 290(10): p. 5940–6.

57.Pascelli, C., et al., Morphological characterization of virus-like particles in coral reef sponges. PeerJ, 2018. 6: p. e5625.

58.Roux, S., et al., VirSorter: mining viral signal from microbial genomic data. PeerJ, 2015. 3: p. e985.

59.Brussow, H., C. Canchaya, and W. D. Hardt, Phages and the evolution of bacterial pathogens: from genomic rearrangements to lysogenic conversion. Microbiol Mol Biol Rev, 2004. 68(3): p. 560–602, table of contents.

60.Seweryn, P., et al., Structural insights into the bacterial carbon-phosphorus lyase machinery. Nature, 2015. 525(7567): p. 68–72.

61.Martinez, A., et al., Metatranscriptomic and functional metagenomic analysis of methylphosphonate utilization by marine bacteria. Front Microbiol, 2013. 4: p. 340.

62.Blin, K., et al., antiSMASH 5.0: updates to the secondary metabolite genome mining pipeline. Nucleic Acids Res, 2019. 47(W1): p. W81-W87.

63.Donia, M. S., et al., Complex microbiome underlying secondary and primary metabolism in the tunicate-Prochloron symbiosis. Proc Natl Acad Sci U S A, 2011. 108(51): p. E1423–32.

64.Schmidt, E. W., et al., Patellamide A and C biosynthesis by a microcin-like pathway in Prochloron didemni, the cyanobacterial symbiont of Lissoclinum patella. Proc Natl Acad Sci U S A, 2005. 102(20): p. 7315–20.

65.Enisoglu-Atalay, V., et al., Chemical and molecular characterization of metabolites from Flavobacterium sp. PLoS One, 2018. 13(10): p. e0205817.

66.Moore, B. S., et al., Plant-like biosynthetic pathways in bacteria: from benzoic acid to chalcone. J Nat Prod, 2002. 65(12): p. 1956–62.

67.Xiang, L. and B. S. Moore, Biochemical characterization of a prokaryotic phenylalanine ammonia lyase. J Bacteriol, 2005. 187(12): p. 4286–9.

68.Reiter, S., et al., Characterization of an orphan type III polyketide synthase conserved in uncultivated ‘Entotheonella’ sponge symbionts Chembiochem, 2019 Aug 20 (Epub ahead of print).

69.Agarwal, V., et al., Complexity of Naturally Produced Polybrominated Diphenyl Ethers Revealed via Mass Spectrometry. Environmental Science & Technology, 2015. 49(3): p. 1339–1346.

70.Finn, R. D., et al., The Pfam protein families database: towards a more sustainable future. Nucleic Acids Res, 2016. 44(D1): p. D279–85.

71.Neubauer, P. R., et al., A flavin-dependent halogenase from metagenomic analysis prefers bromination over chlorination. PLoS One, 2018. 13(5): p. e0196797.

72.Agarwal, V. and B. S. Moore, Enzymatic synthesis of polybrominated dioxins from the marine environment. ACS Chem Biol, 2014. 9(9): p. 1980–4.

73.Fu, X. O., et al., Enzyme-Inhibitors - New and Known Polybrominated Phenols and Diphenyl Ethers from 4 Indo-Pacific Dysidea Sponges. Journal of Natural Products, 1995. 58(9): p. 1384–1391.

74.Miyazaki, R., et al., Comparative genome analysis of Pseudomonas knackmussii B13, the first bacterium known to degrade chloroaromatic compounds. Environ Microbiol, 2015. 17(1): p. 91–104.

75.Schwien, U., et al., Degradation of Chlorosubstituted Aromatic-Compounds by Pseudomonas Sp Strain-B13 - Fate of 3,5-Dichlorocatechol. Archives of Microbiology, 1988. 150(1): p. 78–84.

76.Seo, S. H. and S. D. Lee, Altererythrobacter marensis sp. nov., isolated from seawater. Int J Syst Evol Microbiol, 2010. 60(Pt 2): p. 307–11.

77.Phale, P. S., B. A. Shah, and H. Malhotra, Variability in Assembly of Degradation Operons for Naphthalene and its derivative, Carbaryl, Suggests Mobilization through Horizontal Gene Transfer. Genes (Basel), 2019. 10(8).

78.Hentschel, U., et al., Microbial diversity of marine sponges. Prog Mol Subcell Biol, 2003. 37: p. 59–88.

79.Huang, I. S. and P. V. Zimba, Cyanobacterial bioactive metabolites-A review of their chemistry and biology. Harmful Algae, 2019. 86: p. 139–209.

80.Thacker RW, S. S., Host specificity of the symbiotic cyanobacterium Oscillatoria spongeliae in marine sponges, Dysidea spp. Marine Biology, 2003. 142: p. 643–648.

81.Parada, A. E., D. M. Needham, and J. A. Fuhrman, Every base matters: assessing small subunit rRNA primers for marine microbiomes with mock communities, time series and global field samples. Environ Microbiol, 2016. 18(5): p. 1403–14.

82.Callahan, B. J., et al., DADA2: High-resolution sample inference from Illumina amplicon data. Nat Methods, 2016. 13(7): p. 581–3.

83.Bolyen, E., et al., Reproducible, interactive, scalable and extensible microbiome data science using QIIME 2. Nat Biotechnol, 2019. 37(8): p. 852–857.

84.Pruesse, E., J. Peplies, and F. O. Glockner, SINA: accurate high-throughput multiple sequence alignment of ribosomal RNA genes. Bioinformatics, 2012. 28(14): p. 1823–9.

85.Bolger, A.M., M. Lohse, and B. Usadel, Trimmomatic: a flexible trimmer for Illumina sequence data. Bioinformatics, 2014. 30(15): p. 2114–20.

86.Peng, Y., et al., IDBA-UD: a de novo assembler for single-cell and metagenomic sequencing data with highly uneven depth. Bioinformatics, 2012. 28(11): p. 1420–8.

87.Langmead, B. and S. L. Salzberg, Fast gapped-read alignment with Bowtie 2. Nat Methods, 2012. 9(4): p. 357–9.

88.Li, H., et al., The Sequence Alignment/Map format and SAMtools. Bioinformatics, 2009. 25(16): p. 2078–9.

89.Podell, S. DarkHorse2 GitHub Repository. 2017; Available from: https://github.com/spodell/Darkhorse2.

90.Podell, S. and T. Gaasterland, DarkHorse: a method for genome-wide prediction of horizontal gene transfer. Genome Biol, 2007. 8(2): p. R16.

91.Podell, S., et al., Assembly-driven community genomics of a hypersaline microbial ecosystem. PLoS One, 2013. 8(4): p. e61692.

92.Myers, E. W., et al., A whole-genome assembly of Drosophila. Science, 2000. 287(5461): p. 2196–204.

93.Koren, S., et al., Canu: scalable and accurate long-read assembly via adaptive k-mer weighting and repeat separation. Genome Res, 2017. 27(5): p. 722–736.

94.Quast, C., et al., The SILVA ribosomal RNA gene database project: improved data processing and web-based tools. Nucleic Acids Res, 2013. 41(Database issue): p. D590–6.

95.Price, M. N., P. S. Dehal, and A. P. Arkin, FastTree 2—approximately maximum-likelihood trees for large alignments. PLoS One, 2010. 5(3): p. e9490.

96.Rambaut, A. FigTree, version 1.4.3. In: Molecular evolytion phylogenetics and epidemiology. 2016; Available from: http://tree.bio.ed.ac.uk/software/figtree/ Accessed October, 2018.

97.Segata, N., et al., PhyloPhlAn is a new method for improved phylogenetic and taxonomic placement of microbes. Nat Commun, 2013. 4: p. 2304.

98.Markowitz, V. M., et al., IMG/M: the integrated metagenome data management and comparative analysis system. Nucleic Acids Res, 2012. 40(Database issue): p. D123–9.

99.Bland, C., et al., CRISPR recognition tool (CRT): a tool for automatic detection of clustered regularly interspaced palindromic repeats. BMC Bioinformatics, 2007. 8: p. 209.

100.Haft, D. H., et al., TIGRFAMs and Genome Properties in 2013. Nucleic Acids Res, 2013. 41(Database issue): p. D387–95.

101.Chen, I. A., et al., IMG/M: integrated genome and metagenome comparative data analysis system. Nucleic Acids Res, 2017. 45(D1): p. D507-D516.

102.Dewoody, R. S., P.M. Merritt, and M. M. Marketon, Regulation of the Yersinia type III secretion system: traffic control. Front Cell Infect Microbiol, 2013. 3: p. 4.

103.Agarwal, V., et al., Biosynthesis of polybrominated aromatic organic compounds by marine bacteria. Nat Chem Biol, 2014. 10(8): p. 640–7.

Due to technical limitations, table 1 is only available as a download in the supplemental files section

Download PDF

Journal Publication

published 23 Jun, 2020

Read the published version in Microbiome →

Editorial decision: Minor revision
10 Feb, 2020
Review #1 received at journal
15 Dec, 2019
Review #2 received at journal
06 Dec, 2019
Reviewer #2 agreed at journal
05 Dec, 2019
Reviewer #1 agreed at journal
28 Nov, 2019
Reviewers invited by journal
27 Nov, 2019
Editor assigned by journal
11 Nov, 2019
Submission checks completed at journal
10 Nov, 2019
Editor invited by journal
10 Nov, 2019
First submitted to journal
06 Nov, 2019

You are reading this older preprint version

Read the latest preprint version →

A genomic view of trophic and metabolic diversity in clade-specific Lamellodysidea sponge microbiomes

Status:

Journal Publication

Version 1

Abstract

Figures

background

results

Sponge samples and metagenomic assemblies

Taxonomic classification of assembled genomes

Relative abundance comparisons between samples

Predicted bacterial lifestyles

Viral sequences and phage defense

Secondary metabolite pathways

discussion

conclusions

methods

abbreviations

declarations

references

Table 1

Supplementary Files

Status:

Journal Publication

Version 1