Prokaryotic communities of the French Polynesian sponge Dactylospongia metachromia display a site-specific and stable diversity during an aquaculture trial

Prokaryotic communities of the French Polynesian sponge Dactylospongia metachromia display a site-specific and stable diversity during an aquaculture trial | Research Square window.SnipcartSettings = { analytics: { enabled: false } }; (function() { var accessVector = localStorage.getItem('access_vector') || ''; window.dataLayer = window.dataLayer || []; if (accessVector) { window.dataLayer.push({ user: { profile: { profileInfo: { snid: accessVector } } } }); } })(); (function(w,d,s,l,i){w[l]=w[l]||[];w[l].push({'gtm.start':new Date().getTime(),event:'gtm.js'});var f=d.getElementsByTagName(s)[0],j=d.createElement(s),dl=l!='dataLayer'?'&l='+l:'';j.async=true;j.src='https://www.googletagmanager.com/gtm.js?id='+i+dl;f.parentNode.insertBefore(j,f);})(window,document,'script','dataLayer','GTM-K279D39R'); Browse Preprints In Review Journals COVID-19 Preprints AJE Video Bytes Research Tools Research Promotion AJE Professional Editing AJE Rubriq About Preprint Platform In Review Editorial Policies Our Team Advisory Board Help Center Sign In Submit a Preprint Cite Share Download PDF Research Article Prokaryotic communities of the French Polynesian sponge Dactylospongia metachromia display a site-specific and stable diversity during an aquaculture trial Mathilde Maslin, Benoît Paix, Niels Windt, Rohani Ambo-Rappe, and 4 more This is a preprint; it has not been peer reviewed by a journal. https://doi.org/ 10.21203/rs.3.rs-4002744/v1 This work is licensed under a CC BY 4.0 License Status: Published Journal Publication published 11 Apr, 2024 Read the published version in Antonie van Leeuwenhoek → Version 1 posted 7 You are reading this latest preprint version Abstract Dynamics of microbiomes through time are fundamental regarding survival and resilience of their hosts when facing environmental alterations. As for marine species with commercial applications, such as marine sponges, assessing the temporal change of prokaryotic communities allows us to better consider the adaptation of sponges to aquaculture designs. The present study aims to investigate the factors shaping the microbiome of the sponge Dactylospongia metachromia (Order Dictyoceratida, family Thorectidae), in a context of aquaculture development in French Polynesia (Rangiroa, Tuamotu archipelago). A temporal approach targeting explants collected during farming trials revealed a relative high stability of the prokaryotic diversity, meanwhile a complementary biogeographical study confirmed a spatial specificity amongst samples at different longitudinal scales. Results from this additional spatial analysis confirmed that differences in prokaryotic communities might first be explained by environmental changes (mainly temperature and salinity), while no significant effect of the host phylogeny was observed. The core community of D. metachromia is thus characterized by a high spatiotemporal constancy, which is a good prospect for the sustainable exploitation of this species towards drug development. Indeed, a microbiome stability across locations and throughout the farming process, as evidenced by our results, should go against a negative influence of sponge translocation during in situ aquaculture. holobiont marine sponges microbiome farming trials biogeography French Polynesia Figures Figure 1 Figure 2 Figure 3 Figure 4 Introduction Sponges are ubiquitous filter-feeders involved in key ecological processes such as nutrient cycling, energy transfer towards higher trophic levels, substratum stability, and habitat for a wide range of associated marine organisms (Bell 2008 ; Pawlik and McMurray 2020 ). They have been farmed by humans since Greek antiquity, alternatively used as hygienic tools, household devices, medical instruments, artistic material, and biopolymers (Ehrlich and Worch 2007 ; Voultsiadou 2007 ). Over the past decades, sponges have been increasingly studied as prolific sources of nature-based pharmaceuticals (Sipkema et al. 2005a ; Lowe et al. 2016 ; Maslin et al. 2021 ). For instance, sponge farming programs are an important practice to grow enough biomass to extract their natural products of interest while conserving natural stocks (Duckworth 2009 ; Bierwirth et al. 2022 ). To this end, in situ and ex situ trials have been conducted in a wide range of geographical zones, for example with Neopetrosia sp., Stylissa massa and Callyspongia biru in the Indo-Pacific region (de Voogd 2007 ; Schiefenhövel and Kunzmann 2012 ), Axinella damicornis , Ircinia variabilis , Chondrosia reniformis , and Crambe crambe in the Mediterranean Sea (van Treeck et al. 2003 ; Ternon et al. 2017 ), and Discodermia dissoluta in the Caribbean (Ruiz et al. 2013 ). Increased specialized metabolite yields were for instance observed in Negombata magnifica from the Red Sea (Hadas et al. 2005 ), Dysidea avara from the Mediterranean (De Caralt et al. 2010 ), as well as some New Zealand species (Duckworth and Battershill 2003a , b ), leading to important economical perspectives (Sipkema et al. 2005a ; Osinga et al. 2010 ). Sponges host a large diversity of microorganisms, which are often described as host-specific (Thomas et al. 2016 ). Together, they form a functional entity named the “sponge holobiont”, whereby both host and microbiome display tight interactions contributing to their overall functioning (Pita et al. 2018 ; Slaby et al. 2019 ). The sponge holobiont concept has allowed to better understand the role of sponge prokaryotic communities in the production of various biosynthetic pathways (Wilson et al. 2014 ; Harvey et al. 2015 ). As an example, sponge-associated bacteria can contribute to the production of bioactive compounds (Piel et al. 2004 ; Indraningrat et al. 2016 ), with chemical properties acting against biofoulers (Aguila-Ramírez et al. 2014 ), predators (e.g. fishes; Thoms et al. 2007 ; Slaby et al. 2019 ) or competitors (e.g. corals competing for space; Thoms et al. 2007 ). Sponge microbiomes may also contribute to the nutrient cycling and energy transfer roles of their host in their environment (de Goeij et al. 2017 ; Pita et al. 2018 ). In a global context of environmental changes such as ocean acidification, immune functions associated with the microbiome can influence the survivorship and the resilience of the sponge holobiont to future climate scenarios (Posadas et al. 2022 ). Considering these key roles, the dynamics of the microbiome diversity appears as an essential factor to be monitored during transplantation and farming programs (Mohamed et al. 2008a , b ; Thomas et al. 2016 ). Such investigations provide a better holistic understanding of the parameters linked to the sponge health and thus, to the farming performance itself under a context of environmental pressures (Pita et al. 2018 ). To date, these dynamics have been monitored only under few aquaculture trials. Alpha-diversity increases of bacterial communities from Mycale laximissa and Ircinia strobilina (both collected from the Gulf of Mexico) were observed under closed aquaculture systems (Mohamed et al. 2008a , b ). Another example is related to the transplantation of Gelliodes obtusa to fish farms in Philippines and revealed a stable bacterial community while facing eutrophic pressures (e.g. high concentrations of ammonium and phosphate; Baquiran and Conaco 2018 ). The Indo-Pacific sponge Dactylospongia metachromia (Demospongiae: Dictyoceratida: Thorectidae) was part of a farming program for the extraction of several bioactive compounds in French Polynesia (Maslin 2022 ). This sponge species is known to produce two secondary metabolites of interest named ilimaquinone (Capon and MacLeod 1987 ) and 5-epi-ilimaquinone (Carte et al. 1985 ). Promising pharmaceutical properties were evidenced for these two compounds, including anticancer, antimicrobial, anti-HIV, antiparasitic, anti-diabetes and anti-inflammatory properties (Radeke et al. 1999 ; Lu et al. 2007 ; Kochanowska et al. 2008 ; Leary et al. 2009 ; Yasuhara-Bell and Lu 2010 ; Pereira et al. 2013 ; Daletos et al. 2014 ; Hitora et al. 2021 ). This species has a wide distribution range from the Red Sea to French Polynesia (Hooper et al. 2002 ; Lim et al. 2016 ; Cleary et al. 2020 ), suggesting an important adaptation potential to different environments. However, the microbial diversity associated with D. metachormia is still underexplored to date. Using a denaturing gradient gel electrophoresis fingerprinting approach, Jeong and Park ( 2013 ) described the bacterial community from Micronesian D. metachromia samples, while a Next Generation Sequencing method was undertaken by Cleary et al. ( 2020 ) for two D. metachromia samples from the Red Sea. In both studies, the most represented bacterial phyla were Actinobacteria, Proteobacteria, Chloroflexi, Acidobacteria, and Gemmatimmonadetes which are commonly abundant in so-called high microbial abundance (HMA) sponges. These sponges harbor prokaryotic communities with higher diversities compared to Low Microbial Abundance (LMA) sponges (Moitinho-Silva et al. 2017 ). As the temporal dynamics of D. metachromia microbiome have never been investigated before, the present study aimed to monitor its diversity during the farming experiment conducted in French Polynesia (Maslin 2022 ). Subsequently, the objective was also to identify the major environmental factors (i.e. temperature, salinity, chlorophyll a and several nutrients) shaping the potentially-observed microbial changes. Considering the good survival and growth rates of the sponge host observed during trials in sites where it naturally occurs, we hypothesized a temporally stable microbiome. Additionally, a multi-scale biogeographical approach was conducted to have a broader understanding of the relative importance of these abiotic parameters, as well as the role of the host phylogeny. The core community associated with D. metachromia was also analyzed in order to look for essential sponge-related symbionts that could either be shaped by such variables or be broadly stable over time and space. Materials and Methods Sampling and biological materials from the farming trials The farming trials were conducted at Rangiroa, an atoll from the Tuamotu archipelago. The reef site named “Avatoru” (Fig. 1 ) was chosen, in consultation with the maritime authorities, to perform the aquaculture trials since D. metachromia was found naturally growing in that area. The sponge exploitation remained under the regulations of a permanent operating license from the Direction des Ressources Marines (DRM, Tahiti) during the whole study. The aquaculture frames were directly anchored to the reef as supports of sponge explants coming from nearby parental specimens. These structures consisted of three horizontal, table-shaped iron frames submerged between 16 and 18m depth where the abundance of wild D. metachromia populations was found the highest ( Figure S1 ). Four different sponge explants (roughly 4*4*4cm) were excised in December 2019 from the same original sponge. An additional piece was directly sampled as an uncultivated control, named “Donor”. The four explants were attached to the same frame using a nylon rope directly threaded through their spongin skeleton ( Figure S2 ). Their sampling was then carried out using SCUBA at four distinct times, spaced by 3 months. During every field survey, survival and growth rates of the new sponge generations were calculated and one explant originally coming from the same donor sponge was removed to be subsampled in triplicates. To do so, pieces of D. metachromia specimens were manually excised with a knife to prevent squeezing the animals and transferred into a plastic bag filled with the surrounding seawater. Only small fragments of about 2*2cm were removed from each selected individual, leaving sufficient pinacoderm (i.e. outer layer of sponge cells) intact to allow regeneration after damage. Samples were named according to their farming duration as follows: “9 months” for those collected in September 2020, “12 months” for December 2020, “15 months” for March 2021 and “18 months” for May 2021 ( Table S1 ). Subsamples were cut using sterilized tweezers and scalpels for the barcoding analysis of the sponge and the metabarcoding study of the prokaryotic community. They were subsequently preserved in sterilized plastic vials filled with 96% ethanol and conserved at -20°C until DNA extraction. Sampling and biological materials for the biogeographical complementary study Different geographical scales were considered for the biogeographical study. A first approach was to conduct a global comparison, with samples collected from four different regions: the Red Sea (Saudi Arabia), the Indian Ocean (Maldives), the Makassar Strait (Indonesia), and the South Pacific Ocean (French Polynesia) at different sampling times from 2014 to 2019 (Fig. 1 A, Table S2 ). Then, a regional scale was considered by targeting three different French Polynesian archipelagos. The sampling of French Polynesian sponges was organized during an oceanographic campaign led by the French “Institut de Recherche pour le Développement” (IRD) from 19th October to 23th November 2018 onboard the Alis oceanographic vessel. Subsamples were taken from donor sponge populations located at 20m depth on the barrier reef of six different islands, encompassing more than 4000km 2 and covering three neighboring archipelagos: Rangiroa, Makemo, Raroia, Tematangi (Tuamotu), Mangareva (Gambier) and Tetiaroa (Society) (Fig. 1 B, Table S2 and S3 ). To analyze differences at a higher geographical resolution, considering differences in terms of prevalent currents and oceanic exposition, four specific sites were chosen within the atoll of Rangiroa (Fig. 1 C). Two locations named “Tiputa” and “Avatoru” were located on the external reef on the Western sides of both Tiputa and Avatoru passes. The two remaining areas named “Lagon 1” and “Lagon 2” were located within the lagoon, both facing the Avatoru pass in order to benefit from the semi-oceanic conditions allowed by the daily water exchange. The sampling and subsampling were then processed as previously described for the study of farming trials. Environmental parameters Different devices were used to record temperature and salinity conditions throughout the farming experiment in Rangiroa. First, a temperature data logger (iBee 22L Thermo Button, Plug and Track, Proges Plus, France) was directly attached to one of the frames supporting the explants, to record seasonal temperature variations every hour with an accuracy of 0,1°C. The ThermoTrack PC V8 software was used for device calibration and data processing. Then, for salinity measurements in Practical Salinity Unit (PSU), the AquaSonde-5000 submersible multi-parameter probe (Aquaread, SDEC, France) was deployed using SCUBA at the beginning of each field survey. It was previously calibrated in the laboratory according to the manufacturer's recommendations and using the SondeLink-PC software. Lastly, for the remaining French Polynesian islands investigated in this study except Makemo due to a lack of working probes, temperature and salinity were automatically registered during each stopover using factory sensors directly operated onboard the Alis vessel. For nutrient analyses, a subsample of 20mL of seawater was fixed with 20µL of HgCl 2 in polyethylene scintillation vials (Dominique Dutscher SAS). The vials were vigorously shaken, and then kept at 4°C and away from light exposure. Nutrient content analyses of dissolved phosphate (PO 4 3− ), nitrogen oxides (NO x ), and orthosilicic acid (Si[OH] 4 ) were performed by the Laboratoire des Moyens Analytiques (LAMA, Nouméa, New Caledonia) using a SEAL AA3 HR AutoAnalyzer. For chlorophyll a measurement, a subsample of 500 mL of seawater was filtered on a Whatman® GF/F glass-fiber filter (0,7 nominal pore size; 47 mm diameter) to prevent the loss of picoplankton (Aminot and Rey 2001 ). Filters were then placed in aluminum foils and stored at -20°C. Extraction of chlorophyll a from phytoplanktonic cells was performed using 90% acetone after mechanical grinding of the filters into glass tubes. The 6-mL extracts were then stored overnight in darkness at 4°C and centrifuged the next day a few hours prior to analysis. Concentrations were measured by a Trilogy® Laboratory Fluorometer (Turner Designs) using the chlorophyll a Non-Acidification module. Calibration of the fluorometer was made using a Prim’Light SECOMAM® spectrophotometer. Barcoding of sponge samples For phylogenetic analyses, DNA was extracted from all samples (except those from Lagon 1 and Lagon 2) using the Qiagen DNeasy Blood & Tissue kit and following manufacturer protocols for spin column extractions. The sponge 28S rRNA gene was amplified using the primers 28S-C2-fwd (5′-GAAAAGAACTTTGRARAGAGAGT-3′) and 28S-D2-rev (5′-TCCGTGTTTCAAGACGGG-3′) (Chombard et al. 1998 ) in a 25 µl reaction containing: 7.9µl distilled water (sterile milliQ); 5µl 5x Thermo Fisher Phire Green; 1.0µl 25mM MgCl 2 ; 1.0µl 10mg.ml − 1 Life BSA; 5.0µl 5x Qiagen Q-solution; 1.3µl 10 pMol.µl − 1 forward primer; 1.3µl 10 pMol.µl − 1 reverse primer; 0.5µl 2.5mM dNTP; 0.5µl Thermo Fisher Phire II HS Taq and 1.5µl 0.05-0.5ng.µl − 1 DNA template, with the following amplification parameters: initial denature of 98℃ for 30 sec, followed by 40 cycles of 98℃ for 10 sec, 52.5℃ for 10 sec and 72℃ for 15 sec and a final extension of 72℃ for 5 min. Sponge mitochondrial cytochrome oxidase subunit 1 (COI mtDNA) gene was amplified using the primers dgLCO1490 (5′-GGTCAACAAATCATAAAGAYATYGG-3’) and dgHCO2198 (5′- TAAACTTCAGGGTGACCAAARAAYCA-3′) (Meyer et al. 2005 ) in a 25 µl reaction containing: 12.9µl distilled water (sterile milliQ); 5µl 5x Thermo Fisher Phire Green; 1.0µl 25mM MgCl 2 ; 1.0µl 10mg.ml − 1 Life BSA; 1.3µl 10 pMol.µl − 1 forward primer; 1.3µl 10 pMol.µl − 1 reverse primer; 0.5µl 2.5mM dNTP; 0.5µl Thermo Fisher Phire II HS Taq and 1.5µl 0.05-0.5ng.µl − 1 DNA template, with the following amplification parameters: initial denature of 98℃ for 30 sec, followed by 40 cycles of 98℃ for 10 sec, 55℃ for 10 sec and 72℃ for 15 sec and a final extension of 72℃ for 5 min. PCR product yield was checked using Invitrogen E-gels 96 2% agarose. PCR products were sent to Baseclear B.V. in the Netherlands for forward and reverse sanger sequencing using the same primers than those of the PCR amplification. The resulting sequences were analyzed and edited in Geneious Prime 19.2.3. Reads were assembled using a De Novo Assembly. Contigs were stripped of primer sequences and analyzed for reverse complements and quality. Sequences were aligned in Geneious using MAFFT v7.450 (Katoh et al. 2002 ). A phylogenetic tree was constructed using MrBayes 3.2.6 (Huelsenbeck et al. 2001 ). A GTR substitution model with ‘invgamma’ rate variation was preferred, setting the chain length to 1100000 with a subsampling frequency of 200 and a burn-in length of 100000. A first Mantel test was conducted to test the correlation between the phylogenetic distance matrix (based on the 28S rRNA gene) and the geographical distances matrix. The latter was obtained using the GPS coordinates of each sampling site and the earth.dist() function from the “fossil” R package (Vavrek 2011 ). DNA extractions, library preparation and high throughput sequencing of 16S rRNA gene amplicons. For the microbiome analysis, DNA was extracted using the FastDNA™ SPIN Kit for Soil (MP Biomedicals, Inc.) following the manufacturer instructions. Sponge tissues were cut into small pieces of approximately 3*1*0.5mm using sterilized tweezers and scalpel blades. An extraction blank, in which no tissue was added to the tube, was also included. The library preparation was conducted through a two-step PCR protocol for all samples, together with the extraction blank and two negative controls (mQ water instead of template DNA). For the first PCR, the V4-V5 regions of the 16S rRNA gene were targeted and amplified (Parada et al. 2016 ). We used 515F-Y/926R primers and the KAPA HiFi HotStart Ready Mix PCR Kit (Roche Molecular Systems, Inc.). Reactions were performed in a T100 Thermal Cycler (Bio-Rad, Hercules, CA, United States). The following thermal cycling scheme was set up: initial denaturation at 95°C for 3 min, 30 cycles of denaturation at 98°C for 20 s, annealing at 50°C for 30 s, followed by extension at 72°C for 30 s. The final extension was carried out at 72°C for 5 min. PCR products from the samples were checked using an E-Gel™ (agarose gels at 2%) and the absence of amplification was validated for the negative controls and the extraction blank. PCR products were then cleaned with NucleoMag NGS-Beads (bead volume at 0.9 times the total volume of the sample, Macherey Nagel, Düren, Germany) using the VP 407AM-N 96 Pin Magnetic Bead Extractor stamp (V&P Scientific, San Diego, CA, United States). Through a second PCR, 3 µL of the cleaned PCR products were then amplified and labeled using the MiSeq Nextera XTDNA library preparation kit (Illumina, San Diego, CA, United States), with the same thermal cycling scheme limited to 8 cycles. PCR products were then analyzed with the Fragment Analyzer Agilent 5300 using the DNF-910-33 dsDNA Reagent Kit (35–1,500 bp) protocol (Agilent Technologies, Santa Clara, CA, United States) to confirm successful labeling of the DNA fragments. Negative controls and extraction blanks remained negatives. The pooling at the equimolar concentration was performed with QIAgility (Qiagen, Hilden, Germany). The final pool was then cleaned with NucleoMag NGSBeads, eluted in Milli-Q and the DNA concentration was verified using Tapestation 4150 (Kit HSD 5000, Agilent Technologies, Santa Clara, CA, United States). The sequencing was performed on an Illumina MiSeq V3 PE300 platform at BaseClear B.V. (Leiden, the Netherlands). 16S rRNA gene metabarcoding data processing The raw reads were first treated by BaseClear B.V. for demultiplexing (using bcl2fastq version 2.20, Illumina) and filtering based on two quality controls (using Illumina Chastity filtering, and a PhiX control signal filtering). The following reads were then processed with the DADA2 workflow, allowing an inference to Amplicon Sequence Variant (ASV) (Callahan et al. 2016a , 2017 ). The “dada2” R package was used following the workflow described in Callahan et al. ( 2016b ) and guidelines described in the online tutorial ( https://benjjneb.github.io/dada2/tutorial.html ). The parameters used for filtering and trimming reads were the following ones: truncation length of 290 and 240 base pairs for forward and reverse reads, respectively, maxN = 0, maxEE = 2 and truncQ = 2. The error learning step was based on 153 188 440 and 126 776 640 bases for forward and reverse reads, respectively. After the construction of the ASV table, chimeric sequences were filtered and taxonomic assignment was performed using the Silva v138 reference database (Quast et al. 2013 ). The ASV and taxonomy tables produced by the pipeline were then combined into a phyloseq object, together with the sample metadata table using the “phyloseq” R package (McMurdie and Holmes 2013 ). The dataset was then filtered by removing all sequences affiliated to Eukaryota, chloroplast and mitochondria. Data was then decontaminated with two negative controls and the extraction blank used as control samples, using the “decontam” R package (Davis et al. 2018 ). The dataset was then divided into two subsets, one with the samples for the farming experiment and the second one with those for the biogeographical study. The sampling time (Donor, 9 months, 12 months, 15 months, and 18 months) and sampling sites were used as factors for the statistical analyses ( Table S1 and S2 ). Metabarcoding and statistical analyses After the filtering of reads affiliated to Eukaryota, chloroplast and mitochondria (representing together an average of 0.27% of all reads per sample, SD = ± 0.34%), rarefaction curves performed with the “vegan” R package (Oksanen et al. 2019 ), reached a plateau for all samples ( Figure S3 ). Since the default DADA2 pipeline is performing a removal of singletons and thus considering the potential loss of very rare taxa, these rarefactions curves and the corresponding α -diversity metrics should be interpreted with caution. Despite the removal of these singletons, we remain confident about the relevant coverage of the richness displayed in the whole study. The α -diversity measures were estimated with Chao1 (estimated richness), Pielou (evenness) and Shannon (both richness and evenness) indices using the “ade4” and “phyloseq” R packages (Dray and Dufour 2007 ; McMurdie and Holmes 2013 ) and the rarefied datasets (rarefaction performed to the minimum library size, i.e. 31559 and 34103 reads, for the farming experiment and the biogeographical datasets respectively). Depending on the normal distribution assessed using a Shapiro test, the significant effect of the sampling times and sites on the diversity metrics was tested using ANOVA followed by HSD Tukey’s tests as parametric tests, as well as Kruskal-Wallis followed by pairwise Wilcoxon tests as non-parametric ones. Following recommendations for compositional approaches from McMurdie and Holmes ( 2014 ) and Gloor et al. ( 2017 ), all other analyses were conducted without rarefaction, using the phyloseq R package and the datasets normalized to the total number of sequences per sample. For both datasets, the β -diversity of all samples was analyzed with a non-metric multidimensional scaling (NMDS) using the Bray-Curtis dissimilarity. Differences of β -diversity between groups were statistically checked with one-way PERMANOVA tests followed by pairwise Adonis tests. Distance-based Redundancy Analysis (db-RDA) using the Bray-Curtis dissimilarity was conducted to investigate the effect of environmental parameters on the β -diversity of the prokaryotic community of D. metachromia , using the “vegan” R package. A first analysis was performed with all samples from the farming dataset and tested with all available environmental parameters (temperature, salinity, PO 4 3− , NO x , Si[OH] 4 , and chlorophyll a ). Since some parameters were not described for all sampling sites from the biogeographical study, the analysis was conducted on a restricted dataset with a selection of samples for which the temperature, salinity, PO 4 3− , NO x and Si(OH) 4 were measured ( Table S2 and S3 ). This restricted biogeographical dataset includes only samples from French Polynesian sites (except Makemo). The ordiR2step () function was used to verify that all these environmental factors can be selected as significant constraints to be used as explanatory variables. The overall significance of the two models was tested using the anova.cca () function with 999 permutations. Following db-RDA, variance partitioning analyses (Borcard et al. 1992 ; Cleary et al. 2022 ) were conducted using the varpart () function from the “vegan” package for both the farming and the restricted biogeographical datasets. These analyses were performed to assess the relative percentage of the variance of the β -diversity explained by the selected environmental parameters, using the formula (~ temperature, ~ salinity, ~ PO 4 3− + NO x + Si[OH] 4 ) for both datasets. Additionally, the variance of the β -diversity explained by environmental parameters was also assessed through a variance partitioning analysis combining the geographical distances. The geographical distances matrix was transformed using the pcnm function and the analysis was conducted using the following formula: (~ temperature + salinity + PO 4 3− + NO x + Si[OH] 4 , ~ pcnm [geographical.dist] $ vectors). For the biogeographical study, Mantel tests (Spearman correlation, 9999 permutations) were conducted within all samples to evaluate the correlation between the β -diversity distance matrix and the geographical distance matrix. Excepted for samples not barcoded (Lagon 1 and Lagon 2, Fig. 1 , Table S2 ), an additional Mantel test was also conducted with the biogeographical dataset to assess the correlation between the β -diversity distance matrix and the phylogenetic distance matrix. The "core community" was obtained using the “microbiome” R package (Lahti et al. 2017 ) with the core () function: a prevalence threshold of 90% in all samples (including both geographical and farming datasets) and an abundance threshold of 0.01% were kept. Percentages of core ASVs numbers and sequences were estimated by comparison with the total community and statistically tested using Shapiro, ANOVA and HSD Tukey’s tests. Results Environmental conditions, and sponge monitoring during farming experiment The water temperature increased from 27.3 to 28.3°C for farmed explants of 9 and 12 months respectively, corresponding to the beginning of the wet and warm tropical season (i.e. September to December) ( Figure S4A ). The peak of 29.3°C was reached for “15 months” samples, followed three months later by a slight decrease of -0.5°C. Indeed, the “18 months” sampling corresponds to the start of the dry and cold season in the tropics (i.e. March-April). The salinity varied between 35.04 and 32 PSU, between December 2019 and March 2021 (i.e. “Donor” and “15 months” respectively, Figure S4A ). Yet, the average salinity dropped in May 2021 (i.e. “18 months”) to 29.9 PSU, corresponding to the return of an El Niño–Southern Oscillation climate pattern (World Meteorological Association 2022 ). Regarding the nutrients, no clear trend was observed in the phosphate plot, as the averages remained very similar throughout the farming period. Nitrogen oxides and orthosilicic acid seem to vary in opposite ways ( Figure S4A ), (NO x ) concentrations being the lowest for farmed samples of “12 months” and “15 months” whereas peaks in (Si[OH] 4 ) averages are noticed for the same months (i.e. December 2020 and March 2021, respectively). Without including the Red Sea, the Maldives, Indonesia and the atoll of Makemo (within French Polynesia) since no environmental parameters were measured for these sites ( Table S2 and S3 ), spatial analyses of all environmental parameters revealed significant differences between French Polynesian islands (Kruskal-Wallis: p < 0.02; Table S4 ). Lower sea temperatures for the southern part of the Tuamotu and the Gambier archipelagos (i.e. Tematangi and Mangareva) were observed compared to warmer northern islands featuring on average 4 to 5°C temperatures ( Figure S4B ). The atoll of Rangiroa was characterized by significantly lower salinity levels (averages ranging between 31 and 33 PSU depending on the location) in contrast to the ones registered for Raroia and Tematangi, usually above 36 PSU (Pairwise Wilcoxon test, Figure S4B ). Phosphate and silicate concentrations appear to be relatively similar for all sampled sites, except for the atoll of Raroia, characterized by significantly higher concentrations compared to all sites except Mangareva and Tetiaroa (Pairwise Wilcoxon test, Figure S4B ). Finally, nitrogen oxide values fall along a north-west/south-east gradient from 0.4µg.L − 1 in Tetiaroa to 0.25µg.L − 1 in Raroia, 0.1µg.L − 1 in Tematangi and barely 0.05µg.L − 1 in Mangareva. High and stable survival rates were observed in Avatoru throughout the farming experiment hereby addressed, ranging between 69–71% ( Figure S5A ). For volumes, a higher peak average reaching 99 ± 59 cm 3 after 12 months of farming can be observed, but no significant changes were determined ( Figure S5B , Kruskal-Wallis, p > 0.05). Our results display an annual growth rate of 100%, as the standardized initial size of each explant was set close to 50 cm 3 . Sponge phylogeny Phylogenetic trees based on 28S rRNA gene and COI mtDNA sequences have been constructed using Dactylospongia metachromia samples collected from the farming trials and additional material from the sponge collections of Naturalis Biodiversity Center. The phylogeny of 28S rRNA genes shows some intraspecific variation with several clades forming within D. metachromia . Despite a significant correlation between the distance matrix based on 28S rRNA gene sequences and the geographical distances, for the whole dataset (Mantel test: r = 0.1785; p = 0.0241), no clear geographical pattern was noticed within the French Polynesian sites ( Figure S6 ). More precisely, these samples formed two different clades, with (i) Makemo, Mangareva, and Raroia samples more related to Indonesia and Red Sea samples, and (ii) samples from Rangiroa, Tematangi, and Tetiaroa. Moreover, COI mtDNA shows no intraspecific variation between all samples, with a tree forming a single clade. Prokaryotic diversity of farmed samples For the farming experiment, no significant differences were observed for Shannon, Chao1, and Pielou indices (Fig. 2 A, 2 B and 2 C, respectively) between samples collected at different time points (“donor samples”, and samples farmed during “9-”, “12-”, “15-” and “18-months”, ANOVA: p > 0.05 for the three indices, Table S5 ). As for the NMDS analysis of β -diversity, distinct clusters of samples appear only for “15 months'' and “18 months” samples, while no clear patterns were observed for “Donor”, “9 months” and “12 months” samples (Fig. 2 D). Despite a significant difference observed with the PERMANOVA test ( p = 0.003), no significant pairwise differences between the sample times were noticed (pairwise Adonis: p > 0.05 for all comparisons). The ASVs scores of the NMDS analysis revealed major ASVs (with relative abundances > 2% on average) affiliated with several phyla including Crenarchaeota (e.g. ASV1), Chloroflexi (e.g. ASV5 and AVV10), the PAUC34f clade (e.g. ASV9), Entheonellaeota (e.g. ASV4) and Acidobacteriota (e.g. ASV7) (Fig. 2 E). No clear clusters of ASVs associated with a specific phylum were observed. However, one minor ASV located on the left side of the NMDS and affiliated to the Acidobacteriota appears to be specifically associated with the cluster of “18 months” samples. The db-RDA analysis shows three distinct clusters (Fig. 2 F): (i) one gathering “Donor” and “15 months” samples, positively associated to the temperature and to a lesser extent chlorophyll a , (ii) one gathering “9 months” and “12 months”, positively and negatively associated to the salinity and temperature, respectively, and (iii) a last one gathering “18 months” samples, positively associated to PO 4 3− and NO x concentrations and negatively to the salinity (ANOVA.CCA: p = 0.002, Table S6B ). The ASVs scores of the db-RDA analysis did not reveal ASV clusters specifically associated with environmental parameters (Fig. 2 G). However, ASV5 affiliated to the Chloroflexi which happens to be amongst the dominant ones, was found positively correlated to the temperature. The variance partitioning analysis with environmental parameters revealed that 3% of the β -diversity variance is explained by temperature, 10% by salinity, and 6% by nutrients (PO 4 3− , NO x and Si[OH] 4 together); while the percentage explained by the combination of salinity and nutrients reach 6%. Temperature, salinity and nutrients combined explain 5% of the variance ( Figure S7A ). At the family level, the prokaryotic community composition did not reveal any major changes between sample times (for the farming experiment, Figure S7B ). Most of the families observed were affiliated with three phyla: Acidobacteriota, Chloroflexi, and Proteobacteria. The main families observed in the whole dataset belong to an unidentified family from the Dehalococcoidia class (Chloroflexi phylum), the Nitrosopumilaceae family (Crenarchaeota phylum), the Caldilineaceae family (Chloroflexi phylum), an unaffiliated family from the candidate phylum PAUC34f and two unaffiliated families from the Acidobacteriae and Vicinamibacteria classes (Acidobacteria phylum). The above are respectively representing 11.5%, 6.7%, 6.0%, 6.0%, 5.5% and 4.9% of the community (SD: ±4.5, ± 5.6, ± 2.8, ± 2.3, ± 2.3 and ± 1.8, respectively). Prokaryotic diversity of the biogeographical samples When focusing on the biogeographical study ( Figure S8A ), significant differences between sampling sites were observed for the three α -diversity indices (ANOVA: p < 0.001, p = 0.001 for Shannon and Pielou, respectively; Kruskal-Wallis: p = 0.007 for Chao1; Table S5 ). For the Shannon and Chao1 indices, the highest values were observed for Tiputa samples ( Figure S8A ). For the Pielou index, a more stable tendency was discerned, despite lower values observed for Raroia compared to Rangiroa (Avatoru) and Temantagi (HSD Tukey’s test, Figure S8A ). The NMDS plot of the samples from the biogeographical study (Fig. 3 A) showed a distinct clustering according to sampling sites (PERMANOVA: p < 0.001, Table S6A ), with main differences observed in the first axis (NMDS1) between French Polynesian and other sites (Red Sea, Indonesia and Maldives). Moreover, differences were observed within the French Polynesian sites on the second axis (NMDS2), with a clear clustering observed between each island (Raroia, Temantangi, Mangareva, Tetiaroa and Rangiroa). When focusing on a smaller geographical scale, with sampling sites from the Rangiroa atoll (Tiputa, Lagon1, Lagon2 and Avatoru), no significant differences of β -diversity were observed, except between Lagon1 and Tiputa (pairwise Adonis test). When focusing on the NMDS plot displaying the ASV scores of the biogeographical study (Fig. 3 B), major ASVs were found to be notably affiliated to the Crenarchaeota (e.g. ASV1), Acidobacteriota, Chloroflexi and PAUC34f clades. While no specific patterns of ASV clusters could be associated with a specific group of samples, ASV41 affiliated with the Crenarchaeota was found to be closely related to the Red Sea samples. Considering environmental parameters as explanatory variables (available only for French Polynesian samples, Fig. 3 C and 3 D), the db-RDA was conducted using the Bray-Curtis dissimilarity. The analysis suggested that temperature and salinity were the main parameters explaining β -diversity differences, together with PO 4 3− concentration to a lesser extent (ANOVA.CCA: p < 0.001, Table S4B ). Temperature and salinity appears positively and negatively associated with Rangiroa samples, respectively, while the opposite situation appears for Tetiaroa, Tematangi and Mangareva samples (Fig. 3 C). Additionally, PO 4 3− concentration was found positively associated with the β -diversity of Raroia samples. No specific patterns were observed within the ASV distribution on the db-RDA plot. However, two major ASVs can be noticed: (i) ASV1 affiliated to the Crenarchaeota and positively associated with the salinity, (ii) ASV5 affiliated to the Chloroflexi, mainly occurring within the Rangiroa samples and positively associated with the temperature. A first variance partitioning analysis ( Figure S8B ) revealed that temperature, salinity, and nutrients [combining PO 4 3− , NO x and Si(OH) 4 ] explain 2%, 6% and 5% of the β -diversity variance, respectively. More importantly, temperature and salinity combined explain 18% of the variance, while these percentages reach 3% and 8% when nutrients are combined either with the temperature or the salinity, respectively. For the biogeographical study, a significant correlation was observed between the biogeographical and β -diversity (Bray-Curtis) dissimilarities (Mantel test: r = 0.8722; p < 0.001). Furthermore, a variance partition analysis was conducted with environmental parameters (temperature, salinity and nutrients combined) and the geographical distances as explanatory factors of the β -diversity variance ( Figure S8C ). Results revealed that 9% of the variance is explained by the environmental parameters, 10% by the geographical distances, and 35% is explained by both datasets together. The residual variance (variance left unexplained) reached 46% for this analysis. Additionally, when comparing distance matrices from the sponge phylogeny and β -diversity analyses, no significant correlations were observed between the two datasets (Mantel test: r = 0.1213; p = 0.0693). Finally, the prokaryotic community composition of the biogeographical samples at the family level did not reveal any major changes across space ( Figure S8D ), with similar percentages of relative abundances observed for the temporal study. Core community of all D. metachromia samples The core community of D. metachromia was defined by keeping only ASVs occurring at least in 90% of all samples (samples from both biogeographical and farming studies together). A total of 18 ASVs were determined as part of the core community ( Table S7 ), representing 0.3% of the total number of all ASVs in both datasets combined. The relative percentage of numbers of core ASVs per sample ( Figure S9 ), ranging between 4.1% and 8.8% (average 6.6%, SD ± 1.1), was not significantly different among sampling times for the farming experiment (ANOVA: p = 0.03, but Tukey’s test: p > 0.05) or among sampling sites for the biogeographical approach (ANOVA: p > 0.05). The relative percentage of sequences of core ASVs (Fig. 4 ), ranging between 5.9% and 14.1% (average 10.6%, SD: ±2.1), did not show significant differences between sampling times (for the farming experiment, ANOVA: p > 0.05) or sampling sites (for the biogeographical study, ANOVA: p > 0.05). These sequences of the 18 core ASVs were affiliated to 13 genera (belonging mostly to unidentified families). The most abundant core ASVs were affiliated to the phyla Acidobacteriota, Chloroflexi, Dadabacteria, Myxococcota, Proteobacteria and the candidate phylum PAUC34f, representing in average 2.3%, 1.7%, 1.8%, 1.4%, 1.4%, 1.5%, of the total community, respectively (SD: ±0.7, ± 0.8, ± 0.7, ± 0.7, ± 0.6, ± 1.0, respectively). Some of these phyla were represented by different classes such as Acidobacteriae, Thermoanaerobaculia, and Vicinamibacteria (Acidobacteria), or Alpha- and Gamma-proteobacteria (Proteobacteria). Discussion As a model of interest for the production of pharmacological promising natural products (such as ilimaquinone and 5-epi-ilimaquinone), Dactylospongia metachromia has been the subject of an aquaculture program in French Polynesia (Maslin 2022 ). The production of sponge natural products and, more generally, the whole sponge physiology, is known to be intimately linked to its microbiome (Indraningrat et al. 2016 ; Bibi et al. 2020 ). When considering the launching of a sponge aquaculture industry, it thus appears essential to study the effects of transplantation and environmental changes on the whole holobiont of the targeted species (Thomas et al. 2016 ; Pita et al. 2018 ). A temporal study revealing an overall stability of the holobiont during farming experiment, despite a moderate effect of environmental changes The analysis of the prokaryotic diversity of D. metachromia was conducted on a temporal dimension relative to the duration of the farming experiment targeting sponge explants on a barrier reef area in Rangiroa, named Avatoru. No significant differences of α - diversity between the different months of farming were obtained, with a higher intra-individual variability suggesting an overall stability of the prokaryotic community. Similarly, no clear differences of β -diversity within the overall community were provided by the NMDS analysis. However, several environmental parameters including temperature, salinity and nutrients appear to have a slight impact within the community, notably impacting a small fraction of it [i.e. two ASVs affiliated to Crenaracheaota (ASV1) and Chloroflexi (ASV5)]. The prokaryotic community associated with marine sponges is commonly described as geographically and temporally stable. This is explained by the importance of symbiotic taxa displaying specific relationships with their sponge hosts compared to other environmental samples such as seawater, sediments or other holobionts (Hentschel et al. 2002 ; Taylor et al. 2004 ; Hill et al. 2006 ; Schmitt et al. 2012b ; Chen et al. 2020 ; Hardoim et al. 2021 ). The HMA-LMA specificity is also linked to distinct prokaryotic composition (Campana et al. 2021b ; Lesser et al. 2022 ), while no clear evidence of stability changes was observed over seasons and years (Erwin et al. 2015 ). Studies have shown that notable differences in microbiota composition at multiple taxonomic levels exist between healthy and diseased sponges or under environmental stress (Webster et al. 2008 ). Sponges belonging to the order Dictyoceratida have been reported as particularly sensitive to pathogen epidemics (Smith 1941 ; Wulff 2001 ). Still, identifying the temporally stable microbial community found in the present work, no potentially harmful bacteria taxa were highlighted in the analysis of the community composition. Zaneveld et al. ( 2017 ) introduced the Anna Karenina principle for sponge-holobionts, stating that an increased dispersion of the β -diversity through time could be linked to multiple stress scenarios led by the emergence of diseases. In this study, the opposite situation seems to occur, with a less dispersed β -diversity for the last two sampling times (15 and 18 months). While the transplantation of wild sponge explants farmed in situ to closed or semi-closed systems has demonstrated important changes in the microbiota (Mohamed et al. 2008a , b ), our study indicates that wild transfers can only induce moderate effects on the sponge holobiont. Such stability throughout natural contexts was also observed for other sponges, namely Tethya bergquistae (Cárdenas et al. 2014 ), Aplysina cavernicola (Thoms et al. 2003 ), Aplysina archeri (Hunting et al. 2013 ) and Hymeniacidon heliophila (Weigel and Erwin 2017 ). Additionally, a similar temporal stability was also recently described within Hymeniacidon perlevis and Suberites massa over seasonal changes (Lamb and Watts 2023 ). Interestingly enough, despite highly heterogeneous ratios between ilimaquinone and 5-epi-ilimaquinone within the explants, no significant variation was either obtained in their concentrations over time or specimens during farming experiments (Maslin 2022 ). Together with the stable α - and β -diversity observed over time, these results corroborate the absence of a negative influence of the aquaculture system during the farming period. This finding allows us to consider a chemical exploitation in line with the positive enlargement of the sponge explants in order to seek an optimal sampling time. Nevertheless, being of bacterial origin or not, the average concentrations of secondary metabolites found in the extracts of D. metachromia are still insufficient to match the needs of the pharmaceutical market towards new drug synthesis (Sipkema et al. 2005b ; Malve 2016 ). A complementary biogeographical approach disclosing site-specific communities and confirming the importance of temperature and salinity Values of α -diversity metrics were found to slightly differ between sampling sites, especially for the Red Sea samples characterized by a lower Shannon index compared to Rangiroa and Tematangi atolls. Such differences for this particular region might enlighten host-specific relationships, since for other sponges such as Xestospongia spp. the α -diversity did not reveal lower values for samples from Saudi Arabia compared to the Australian coast or the Indo-Pacific region (Thomas et al. 2016 ; Swierts et al. 2018 ). However, comparison of prokaryotic richness between different studies should be interpreted with caution since many methodological biases can be involved, from the molecular preparation to the data processing. The β -diversity of the prokaryotic community associated with D. metachromia samples was found to be highly site-specific, whether at the global scale (when considering samples from French Polynesia together with Red Sea, Indonesia and Maldives) or within a particular region (when comparing the different French Polynesian islands). However, this specificity was found to be locally limited when assessing different sampling sites from Rangiroa. Samples from the outer part of the reef of this atoll (Tiputa and Avatoru sites) were exposed to higher currents compared to the ones located within the lagoon and only relying on daily water exchanges through the passes (sites “Lagon” 1 and “Lagon 2”). Such differences in terms of oceanographic conditions do not appear to be involved in the variations of the prokaryotic communities. A similar result was observed for the sponge Cliona delitrix , characterized by a prokaryotic β -diversity only correlated to large geographical distances (> 1000 km), whereas no significant differences were observed more locally (< 300 km) (Easson et al. 2020 ). Many studies have provided insights on correlation between sponge phylogeny and microbiota. In some cases, dominant symbiont selection is believed to be related to sponge host identity (Easson and Thacker 2014 ; O’Brien et al. 2020 ). Several holobiont features like stress flexibility or immunity can be decisive in species adaptation, especially when considering future ocean scenarios with perturbed conditions (Posadas et al. 2022 ). In this sense, tropical patterns in microbial communities amongst marine sponge orders have been revealed, most probably driven not only by vertical transmission but also by abiotic factors (Taylor et al. 2007 ; Schmitt et al. 2012b ). Our study revealed significant correlations between geographical and β -diversity distances, suggesting an influence of both environmental parameters and/or phylogenetic differences between the sponge samples. However, no clear geographical pattern was observed within the phylogenetic trees and no significant correlations were observed between the host phylogeny and the β -diversity distances. This result suggests that the prokaryotic community could rather be shaped by the environment than by the host specificity. Similar results were observed for Xestospongia spp. collected from the Indo-Pacific region, the Mozambique Channel, and the Red Sea, with a prokaryotic community more influenced by the biogeography than by the host-specificity (Swierts et al. 2018 ). Other studies suggest a more important effect of the host-specificity (Meyer and Kuever 2008 ; Busch et al. 2022 ), especially when smaller geographical scales are considered (Cleary et al. 2018 ), or when the phylogenetic comparisons are covering a large sponge diversity with higher taxonomic differences (i.e. family or order levels; Yang et al. 2019 ). As studied for C. delitrix , more insight on the potential effect of host-specificity should be considered at the level of the population genetics (Easson et al. 2020 ). In our study, such level may justify why 10% of the β -diversity variance is explained exclusively by the geographical distance independently from the environmental factors. However, a large proportion of the variance associated with the environmental parameters, whether combined with the geographical distances (35%) or not (9%), confirm their predominant role driving the community structure. When considering the relative effect of the different environmental parameters measured, the biogeographical study confirmed that temperature and salinity constitute two key abiotic factors shaping the β -diversity of the prokaryotic community. The effect of these parameters appears to be more important on the biogeographical study compared to the temporal one, probably in line with the greater variability of these factors through space rather than time. The temperature and salinity were found to have opposite effects on the community structure, explaining their high proportion of the shared variance. The negative relationship betwen temperature and salinity has been considered to play an important role within deep-sea sponge-holobionts (Busch et al. 2022 ), as reported for Geodia barretti (Strand et al. 2017 ; Steffen et al. 2022 ). In these studies, such parameters are consequently often interpreted carefully due to their collinearity with other depth-related factors, such as hydrostatic pressure or oxygen level. However, all samples in our study were collected from shallow waters with an average of 20m (SD: ±5m), suggesting a limited role of the depth in the combined effect of temperature and salinity. Experimental studies have proven that an increase in temperature is one of the main factors impacting sponge holobionts health (Webster et al. 2008 ; Pantile and Webster 2011 ; Pita et al. 2018 ; Slaby et al. 2019 ). Nevertheless, it has also been observed that bacterial communities within Ircinia spp. can be highly stable under large fluctuations of temperature occurring in situ (Erwin et al. 2012 ) or in controlled conditions (Pita et al. 2013 ). As for salinity changes, only few studies have looked for specific effects on sponges, especially at the holobiont scale. Although decreasing salinity gradients were linked to a decline in growth of Cymbastela concentrica , the abundance of its photosymbionts was not significantly affected (Roberts et al. 2006 ). Additionally, a stability of the host health and prokaryotic communities was demonstrated for diverse HMA and LMA species exposed to salinity disturbances over time (Glasl et al. 2018 ). Stability and role of the core community The core concept within holobionts has been defined as the spatio-temporal stable part of the microbiome community, allowing to identify potential key symbionts displaying strong associations with their hosts (Astudillo-García et al. 2017 ; Risely 2020 ; Neu et al. 2021a ). In particular, identifying the dynamics and stability of these core taxa is relevant to better understand the health and resilience of holobionts facing stress or damage, such as in a context of climate change or wild harvest for pharmaceutical valorization (Hutchins et al. 2019 ; Risely 2020 ; Neu et al. 2021a ). However, the parameters to identify the core microbiome can vary a lot depending on the studies, featuring a broad set of context-based parameters, and the lack of a consensus definition of the core-community, often leading to different interpretations (Astudillo-García et al. 2017 ; Astudillo-García et al. 2020 ). In our study, we defined the spatio-temporal core community of D. metachromia by identifying core ASVs with a prevalence of 90% in all samples. This parameter allowed a conservative approach to keep only persistent ASVs across space and time, but also to compare our results with previous studies using this approach (Thomas et al. 2016 ; Astudillo‐García et al. 2020). The threshold returned 18 core ASVs, representing a relatively small set of the overall community in terms of average percentages of sequences (10.6%) and ASV numbers (6.6%). With a threshold of 85% of occurrence, the core microbiome described within 5 sponge species was found to vary from 5 to 20 core OTUs (Thomas et al. 2016 ). Furthermore, a similar threshold (90% of occurrence, with an abundance threshold of 0.1%) returned 4 core OTUs for 20 host sponge species (Astudillo-García et al. 2020 ). Together with these results, our work suggests that a limited number of taxa are core members of sponge microbiomes. However, the number of core ASVs or OTUs depends on sampling size, especially since larger spatio-temporal scales tend to decrease the richness of the core microbiome (Turnbaugh et al. 2007 ). The core bacterial community of D. metachromia was found stable, with no significant changes of relative percentages of sequences and ASV numbers over time and space. Transplantation, thus, seems to have had a very limited impact on the structure of the potential key symbiotic community of the targeted sponge species. Also, the high growth rates of farmed explants could be linked to the overall stability of the core taxa, both indicating a steady fitness of the holobiont. These findings are in agreement with previous studies analyzing sponge-associated core microbiomes, either in tropical regions (Campana et al. 2021b ; Leal et al. 2022 ), temperate waters (Lamb and Watts 2023 ) or polar ecosystems (Happel et al. 2022 ), suggesting that changing environmental conditions might have a relatively low impact in the shaping of the core sponge microbiome. The 18 core ASVs identified belong to 13 genera mostly affiliated to Acidobacteriota, Chloroflexi, PAUC34f and Proteobacteria. Such results, returning rather small globally-distributed taxa, are similar to previous works focusing on core communities defined from multiple species coming from very distinct locations (Schmitt et al. 2012a , b ). The core taxa of D. metachromia within the phylum Acidobacteriota, were assigned to Vicinamibacterale, Thermoanaerobaculales and the uncultured group PAUC26f. Marine Thermoanaerobaculales are not strictly sponge-related, as they have also been extracted from molluscs (Neu et al. 2021b ; Banker et al. 2022 ). Likewise, a core ASV was affiliated to the class bacteriap25 (Myxococcales; Proteobacteria), which is found in soils and terrestrial plants such as some non-indigenous Asteraceae that have been introduced in French Polynesia (Landwehr et al. 2016 ; Monteiro 2020 ). One of the main core ASVs belongs to the candidate phylum PAUC34f, which is described as part of the most important clades directly associated with the cycling of dissolved organic matter (DOM) within the Caribbean sponge Plakortis angulospiculatus (Campana et al. 2021a ). Within the Nitrospinota phylum, a core ASV was affiliated to the candidate clade P9X2b3D02, which was recently also observed in Aplysina caissara restricted to southern Brazil (Hardoim et al. 2021 ). This clade did not occur in the marine sediment sampled nearby or in other coastal sponge species from the Western Atlantic. Along with some Chloroflexi and Dadabacteria taxa, Nitrospinota bacteria are known to participate in the sulfur cycling and are capable of synthesizing various B-vitamins or essential amino acids (Bayer et al. 2018 ; Engelberts et al. 2020 ). Additionally, two classes (S085 and SAR202 clades) from the phylum Chloroflexi were represented within the core ASVs. Clade SAR202 has been observed in the water column of the aphotic zone and is believed to take part in the oxidation of sulfite particles (Mehrshad et al. 2018 ). Recently, this clade was also found to dominate the microbiome of sponges living in extreme conditions, with large variation of pH values, temperature levels and oxygen concentrations in a semi-enclosed reef lagoon (Maggioni et al. 2023b ). Such dominance of SAR202 is suggested to play an important role for providing energy to the host through their DOM assimilation ability, thus helping with its adaptation in a changing environment (Maggioni et al. 2023a , b). In respect to its potential role in DOM turnover and nutrient dynamics, especially in low-light environments, SAR202 might be providing a competitive advantage to the sponge for DOM assimilation and cycling (Bayer et al. 2018 ; Morganti et al. 2022 ). Finally, the most abundant core ASVs belong to the order Dadabacteriales (Dadabacteria). Dadabacteria have been observed from the open ocean, from living systems related to the presence of organic carbon (including marine sponges) or from hydrothermal structures. This order notably includes obligate aerobic and photo-heterotrophic taxa; their ecological distribution throughout the water column could actually be driven by phototrophy (Graham and Tully 2021 ). In sponges, they have only been recently discovered (Schuster et al. 2016 ; Gavriilidou et al. 2023 ), and their ecological involvements remain partially unaddressed. Yet, they appear to degrade microbial dissolved organic matter like phospholipids and complex carbohydrates, but also to produce important compounds for the host-interaction such as vitamins and cofactors (Hug et al. 2016 ; Gavriilidou et al. 2023 ). Conclusion Environmental changes often play a determinant role on the prokaryotic community of a sponge holobiont (Morrow et al. 2016 ; Batista et al. 2018 ; Busch et al. 2022 ). Consequently, it appears essential to better investigate the factors that could shape sponge holobiont dynamics, especially when a commercial outlet for pharmaceutics is tested during farming trials. An overall stable α - and β -diversity of the prokaryotic community was observed for the sponge Dactylospongia metachromia whose explants were farmed for 18 months on the outer reef of Rangiroa (Tuamotu, French Polynesia). These results are in line with the good viability, high survival rates and constant biomass gain of the new sponge generations observed in situ. However, some parameters such as temperature and salinity appear to partially impact a small fraction of the community. The moderated importance of these two abiotic factors was confirmed through a biogeographical study considering a broad regional scale (i.e. Red Sea and Indo-Pacific regions). Such a spatial investigation also revealed that both environmental variables are predominant over the effects of host-phylogeny differences. Hence, salinity and temperature need to be considered for their potential influence on the production of bioactive metabolites by the sponge-microorganism complex (Brinkmann et al. 2017 ). The core microbiota of D. metachromia is composed of bacterial phyla widely distributed among sponges (except the poorly characterized Dadabacteria phylum) and these taxa are known for their DOM assimilation abilities, sulfur cycling and vitamin production. Such findings highlight their global importance for the holobiont functioning. Such taxa do not seem to be strongly affected either by geographical and environmental characteristics or by the duration of the farming trials. Indeed, core communities of sponge holobionts commonly appear to be highly similar alongside a broad spatio-temporal gradient. However, thermal changes have previously been linked to bacterial shifts in sponge-associated microbiome, global issues such as a general rise in sea temperatures are likely to affect such bacterial communities through the disappearance of certain symbionts for the benefit of others (Lemoine et al. 2007 ; Webster et al. 2008 ). In the context of climate change, it becomes ever more important nowadays to keep extending our knowledge of the potential effects of ocean variables on a large diversity of sponge holobionts such as the one targeted by this work, especially considering the stability and function of their core community. In addition, our findings are of interest in a growing context of marine-based chemical substances valuation, notably leading the path to a potential exploitation of this promising sponge species in other geographical regions. Declarations Acknowledgments We want to acknowledge the availability and assistance provided by the whole crew of the Alis vessel, including all scientific members who helped with sample collection and analysis during the “TUAM’2018” oceanographic mission. We are also grateful to the DRM headquarters and field operators, either coming from Tahiti or based in Rangiroa, for their help during the set-up and monitoring of the aquaculture trials. Special thanks to all UMR 241 EIO lab technicians who assisted in the organization of the field missions and the preparation of general equipment. Lastly, we want to show our great appreciation to IRD researcher Dr. Sylvain Petek, who contributed to the general organization of the REDAME project and performed some key steps of the chemical analyses. Funding MM was granted a PhD fellowship at the University of French Polynesia (French "Ministère de l'Enseignement Supérieur, de la Recherche et de l’Innovation"). The oceanographic campaign was majoritarily funded by the IRD and their LEMAR lab team. Additional funds were allocated by the GENAVIR research vessel company and the UPF. Field campaigns in Rangiroa and environmental analyzes within French Polynesia were funded through the REDAME project, financially led by the "Direction des Ressources Marines" (grant number 08858/VP/DRMM) and the Government of French Polynesia (grant number 3299/MTF). The lab work was funded by the NWO-VIDI with project number 16.161.301. Data accessibility COI and 28S rRNA gene sequences were deposited and are publicly available in GenBank under the accession numbers OM980557-OM980598 ( Table S3 ). 16S rRNA gene sequences were deposited and are publicly available in the NCBI Sequences Read Archive (SRA) under the BioProject ID PRJNA917856, accession number. The R script used for all the 16S rRNA gene metabarcoding analysis can be found at https://github.com/BenoitPAIX/Dactylospongia_metachromia_microbiome Author contributions CD, NGM and RH contributed to the initial design and set-up of the REDAME project which financially supported this study. CD and MM collected and prepared the environmental samples during the oceanic campaign in French Polynesian archipelagos. MM performed the installation of farming structures and the monitoring of the survival, growth and environmental data of the sponges in Rangiroa. As for the microbiological analyses, MM and NV conceived the study design and protocols. MM provided the samples from French Polynesia. NV and RA provided the samples from the Red Sea, the Maldives and Indonesia. NW and MM extracted the DNA. BP performed the library preparation for the 16S rRNA gene sequencing. NW and MM prepared the samples for the sponge barcoding. BP analyzed the 16S rRNA gene metabarcoding data. NW and MM analyzed the barcoding data. BP created the graphics. BP and MM prepared the first draft of the manuscript. All authors read, contributed to the revision, and approved the submitted version. References Aguila-Ramírez RN, Hernández-Guerrero CJ, González-Acosta B, et al (2014) Antifouling activity of symbiotic bacteria from sponge Aplysina gerardogreeni . International Biodeterioration & Biodegradation 90:64–70. https://doi.org/10.1016/j.ibiod.2014.02.003 Aminot A, Rey F (2001) Chlorophyll a: Determination by spectroscopic methods. International Council for the Exploration of the Sea (ICES) Astudillo-García C, Bell JJ, Montoya JM, et al (2020) Assessing the strength and sensitivity of the core microbiota approach on a highly diverse sponge reef. Environ Microbiol 22:3985–3999. https://doi.org/10.1111/1462-2920.15185 Astudillo-García C, Bell JJ, Webster NS, et al (2017) Evaluating the core microbiota in complex communities: A systematic investigation. Environ Microbiol 19:1450–1462. https://doi.org/10.1111/1462-2920.13647 Banker RMW, Lipovac J, Stachowicz JJ, Gold DA (2022) Sodium molybdate does not inhibit sulfate-reducing bacteria but increases shell growth in the Pacific oyster Magallana gigas . PLoS One 17:e0262939. https://doi.org/10.1371/journal.pone.0262939 Baquiran JIP, Conaco C (2018) Sponge-microbe partnerships are stable under eutrophication pressure from mariculture. Marine Pollution Bulletin 136:125–134. https://doi.org/10.1016/j.marpolbul.2018.09.011 Batista D, Costa R, Carvalho AP, et al (2018) Environmental conditions affect activity and associated microorganisms of marine sponges. Marine Environmental Research 142:59–68. https://doi.org/10.1016/j.marenvres.2018.09.020 Bayer K, Jahn MT, Slaby BM, et al (2018) Marine sponges as Chloroflexi hot spots: Genomic insights and high-resolution visualization of an abundant and diverse symbiotic clade. mSystems 3:e00150-18. https://doi.org/10.1128/mSystems.00150-18 Bell JJ (2008) The functional roles of marine sponges. Estuarine, Coastal and Shelf Science 79:341–353. https://doi.org/10.1016/j.ecss.2008.05.002 Bibi F, Yasir M, Al-Sofyani A, et al (2020) Antimicrobial activity of bacteria from marine sponge Suberea mollis and bioactive metabolites of Vibrio sp. EA348. Saudi J Biol Sci 27:1139–1147. https://doi.org/10.1016/j.sjbs.2020.02.002 Bierwirth J, Mantas TP, Villechanoux J, Cerrano C (2022) Restoration of marine sponges—what can we learn from over a century of experimental cultivation? Water 14:1055. https://doi.org/10.3390/w14071055 Borcard D, Legendre P, Drapeau P (1992) Partialling out the Spatial Component of Ecological Variation. Ecology 73:1045–1055. https://doi.org/10.2307/1940179 Brinkmann CM, Marker A, Kurtböke Dİ (2017) An Overview on Marine Sponge-Symbiotic Bacteria as Unexhausted Sources for Natural Product Discovery. Diversity 9:40. https://doi.org/10.3390/d9040040 Busch K, Slaby BM, Bach W, et al (2022) Biodiversity, environmental drivers, and sustainability of the global deep-sea sponge microbiome. Nat Commun 13:5160. https://doi.org/10.1038/s41467-022-32684-4 Callahan BJ, McMurdie PJ, Holmes SP (2017) Exact sequence variants should replace operational taxonomic units in marker-gene data analysis. ISME J 11:2639–2643. https://doi.org/10.1038/ismej.2017.119 Callahan BJ, McMurdie PJ, Rosen MJ, et al (2016a) DADA2: High-resolution sample inference from Illumina amplicon data. Nat Methods 13:581–583. https://doi.org/10.1038/nmeth.3869 Callahan BJ, Sankaran K, Fukuyama JA, et al (2016b) Bioconductor Workflow for Microbiome Data Analysis: from raw reads to community analyses Campana S, Busch K, Hentschel U, et al (2021a) DNA-stable isotope probing (DNA-SIP) identifies marine sponge-associated bacteria actively utilizing dissolved organic matter (DOM). Environ Microbiol 23:4489–4504. https://doi.org/10.1111/1462-2920.15642 Campana S, Demey C, Busch K, et al (2021b) Marine sponges maintain stable bacterial communities between reef sites with different coral to algae cover ratios. FEMS Microbiol Ecol 97:fiab115. https://doi.org/10.1093/femsec/fiab115 Capon RJ, MacLeod JK (1987) Revision of the absolute stereochemistry of ilimaquinone. The Journal of Organic Chemistry 52:5059–5060 Cárdenas CA, Bell JJ, Davy SK, et al (2014) Influence of environmental variation on symbiotic bacterial communities of two temperate sponges. FEMS Microbiology Ecology 88:516–527. https://doi.org/10.1111/1574-6941.12317 Carte B, Rose CB, Faulkner DJ (1985) 5-Epi-Ilimaquinone, a metabolite of the sponge Fenestraspongia sp. The Journal of Organic Chemistry 50:2785–2787 Chen ML, Becraft ED, Pachiadaki M, et al (2020) Hiding in plain sight: The globally distributed bacterial candidate Phylum PAUC34f. Front Microbiol 11:376. https://doi.org/10.3389/fmicb.2020.00376 Chombard C, Boury-Esnault N, Tillier S (1998) Reassessment of homology of morphological characters in tetractinellid sponges based on molecular data. Systematic Biology 47:351–366. https://doi.org/10.1080/106351598260761 Cleary DFR, Polónia ARM, de Voogd NJ (2018) Prokaryote composition and predicted metagenomic content of two Cinachyrella Morphospecies and water from West Papuan Marine Lakes. FEMS Microbiology Ecology 94:. https://doi.org/10.1093/femsec/fix175 Cleary DFR, Polónia ARM, Reijnen BT, et al (2020) Prokaryote communities inhabiting endemic and newly discovered sponges and octocorals from the Red Sea. Microb Ecol 80:103–119. https://doi.org/10.1007/s00248-019-01465-w Cleary DFR, Polónia ARM, Swierts T, et al (2022) Spatial and environmental variables structure sponge symbiont communities. Molecular Ecology 31:4932–4948. https://doi.org/10.1111/mec.16631 Daletos G, de Voogd NJ, Müller WEG, et al (2014) Cytotoxic and protein kinase inhibiting nakijiquinones and nakijiquinols from the sponge Dactylospongia metachromia . J Nat Prod 77:218–226. https://doi.org/10.1021/np400633m Davis NM, Proctor DM, Holmes SP, et al (2018) Simple statistical identification and removal of contaminant sequences in marker-gene and metagenomics data. Microbiome 6:226. https://doi.org/10.1186/s40168-018-0605-2 De Caralt S, Sánchez-Fontenla J, Uriz MJ, Wijffels RH (2010) In situ aquaculture methods for Dysidea avara (Demospongiae, Porifera) in the northwestern mediterranean. Marine Drugs 8:1731–1742. https://doi.org/10.3390/md8061731 de Goeij JM, Lesser MP, Pawlik JR (2017) Nutrient fluxes and ecological functions of coral reef sponges in a changing ocean. In: Carballo JL, Bell JJ (eds) Climate Change, Ocean Acidification and Sponges: Impacts Across Multiple Levels of Organization. Springer International Publishing, Cham, pp 373–410 de Voogd NJ (2007) The mariculture potential of the Indonesian reef-dwelling sponge Callyspongia ( Euplacella ) biru : Growth, survival and bioactive compounds. Aquaculture 262:54–64. https://doi.org/10.1016/j.aquaculture.2006.09.028 Dray S, Dufour A-B (2007) The ade4 package: implementing the duality diagram for ecologists. J Stat Softw 22:1–20. https://doi.org/10.18637/jss.v022.i04 Duckworth A (2009) Farming sponges to supply bioactive metabolites and bath sponges: a review. Mar Biotechnol 11:669–679. https://doi.org/10.1007/s10126-009-9213-2 Duckworth A, Battershill C (2003a) Sponge aquaculture for the production of biologically active metabolites: the influence of farming protocols and environment. Aquaculture 221:311–329. https://doi.org/10.1016/S0044-8486(03)00070-X Duckworth AR, Battershill CN (2003b) Developing farming structures for production of biologically active sponge metabolites. Aquaculture 217:139–156. https://doi.org/10.1016/S0044-8486(02)00038-8 Easson CG, Chaves-Fonnegra A, Thacker RW, Lopez JV (2020) Host population genetics and biogeography structure the microbiome of the sponge Cliona delitrix . Ecology and Evolution 10:2007. https://doi.org/10.1002/ece3.6033 Easson CG, Thacker RW (2014) Phylogenetic signal in the community structure of host-specific microbiomes of tropical marine sponges. Front Microbiol 5:532. https://doi.org/10.3389/fmicb.2014.00532 Ehrlich H, Worch H (2007) Sponges as natural composites: from biomimetic potential to development of new biomaterials. In: In: Custódio MR, Lôbo-Hajdu G, Hajdu E, Muricy G (eds) Porifera Research. Biodiversity, Innovation and Sustainability. Livros de Museu Nacional 28, Rio de Janeiro, pp 303–312 Engelberts JP, Robbins SJ, de Goeij JM, et al (2020) Characterization of a sponge microbiome using an integrative genome-centric approach. ISME J 14:1100–1110. https://doi.org/10.1038/s41396-020-0591-9 Erwin PM, Coma R, López-Sendino P, et al (2015) Stable symbionts across the HMA-LMA dichotomy: low seasonal and interannual variation in sponge-associated bacteria from taxonomically diverse hosts. FEMS Microbiol Ecol 91:fiv115. https://doi.org/10.1093/femsec/fiv115 Erwin PM, Pita L, López-Legentil S, Turon X (2012) Stability of sponge-associated bacteria over large seasonal shifts in temperature and irradiance. Applied and Environmental Microbiology 78:7358–7368. https://doi.org/10.1128/AEM.02035-12 Gavriilidou A, Avcı B, Galani A, et al (2023) Candidatus Nemesobacterales is a sponge-specific clade of the candidate phylum Desulfobacterota adapted to a symbiotic lifestyle. ISME J 17:1808–1818. https://doi.org/10.1038/s41396-023-01484-z Glasl B, Smith CE, Bourne DG, Webster NS (2018) Exploring the diversity-stability paradigm using sponge microbial communities. Sci Rep 8:8425. https://doi.org/10.1038/s41598-018-26641-9 Gloor GB, Macklaim JM, Pawlowsky-Glahn V, Egozcue JJ (2017) Microbiome datasets are compositional: and this is not optional. Front Microbiol 8:2224. https://doi.org/10.3389/fmicb.2017.02224 Graham ED, Tully BJ (2021) Marine Dadabacteria exhibit genome streamlining and phototrophy-driven niche partitioning. ISME J 15:1248–1256. https://doi.org/10.1038/s41396-020-00834-5 Hadas E, Shpigel M, Ilan M (2005) Sea ranching of the marine sponge Negombata magnifica (Demospongiae, Latrunculiidae) as a first step for latrunculin B mass production. Aquaculture 244:159–169. https://doi.org/10.1016/j.aquaculture.2004.11.052 Happel L, Rondon R, Font A, et al (2022) Stability of the microbiome of the sponge Mycale ( Oxymycale ) acerata in the Western Antarctic Peninsula. Front Microbiol 13:827863. https://doi.org/10.3389/fmicb.2022.827863 Hardoim CCP, Ramaglia ACM, Lôbo-Hajdu G, Custódio MR (2021) Community composition and functional prediction of prokaryotes associated with sympatric sponge species of southwestern Atlantic coast. Sci Rep 11:9576. https://doi.org/10.1038/s41598-021-88288-3 Harvey AL, Edrada-Ebel R, Quinn RJ (2015) The re-emergence of natural products for drug discovery in the genomics era. Nat Rev Drug Discov 14:111–129. https://doi.org/10.1038/nrd4510 Hentschel U, Hopke J, Horn M, et al (2002) Molecular evidence for a uniform microbial community in sponges from different oceans. Appl Environ Microbiol 68:4431–4440. https://doi.org/10.1128/AEM.68.9.4431-4440.2002 Hill M, Hill A, Lopez N, Harriott O (2006) Sponge-specific bacterial symbionts in the Caribbean sponge, Chondrilla nucula (Demospongiae, Chondrosida). Marine Biology 148:1221–1230. https://doi.org/10.1007/s00227-005-0164-5 Hitora Y, Sejiyama A, Honda K, et al (2021) Fluorescent image-based high-content screening of extracts of natural resources for cell cycle inhibitors and identification of a new sesquiterpene quinone from the sponge, Dactylospongia metachromia . Bioorg Med Chem 31:115968. https://doi.org/10.1016/j.bmc.2020.115968 Hooper JNA, Van Soest RWM, Willenz P (eds) (2002) Systema Porifera. Springer US, Boston, MA Huelsenbeck JP, Ronquist F, Nielsen R, Bollback JP (2001) Bayesian Inference of Phylogeny and Its Impact on Evolutionary Biology. Science 294:2310–2314. https://doi.org/10.1126/science.1065889 Hug LA, Thomas BC, Sharon I, et al (2016) Critical biogeochemical functions in the subsurface are associated with bacteria from new phyla and little studied lineages. Environ Microbiol 18:159–173. https://doi.org/10.1111/1462-2920.12930 Hunting ER, Franken O, Knopperts F, et al (2013) Substrate as a driver of sponge distributions in mangrove ecosystems. Marine Ecology Progress Series 486:133–141. https://doi.org/10.3354/meps10376 Hutchins DA, Jansson JK, Remais JV, et al (2019) Climate change microbiology — problems and perspectives. Nat Rev Microbiol 17:391–396. https://doi.org/10.1038/s41579-019-0178-5 Indraningrat AAG, Smidt H, Sipkema D (2016) Bioprospecting sponge-associated microbes for antimicrobial compounds. Mar Drugs 14:87. https://doi.org/10.3390/md14050087 Jeong IH, Park JS (2013) Bacterial diversity of the South Pacific sponge, Dactylospongia metachromia based on DGGE fingerprinting. The Korean Journal of Microbiology Katoh K, Misawa K, Kuma K, Miyata T (2002) MAFFT: a novel method for rapid multiple sequence alignment based on fast Fourier transform. Nucleic Acids Res 30:3059–3066. https://doi.org/10.1093/nar/gkf436 Kochanowska AJ, Rao KV, Childress S, et al (2008) Secondary metabolites from three Florida sponges with antidepressant activity. J Nat Prod 71:186–189. https://doi.org/10.1021/np070371u Lahti L, Shetty S, Blake T, Salojarvi J (2017) Tools for microbiome analysis in R. Version 1:28 Lamb CE, Watts JEM (2023) Microbiome species diversity and seasonal stability of two temperate marine sponges Hymeniacidon perlevis and Suberites massa . Environmental Microbiome 18:52. https://doi.org/10.1186/s40793-023-00508-7 Landwehr W, Wolf C, Wink J (2016) Actinobacteria and Myxobacteria-two of the most important bacterial resources for novel antibiotics. Curr Top Microbiol Immunol 398:273–302. https://doi.org/10.1007/82_2016_503 Leal CV, Avelino-Alves D, Salazar V, et al (2022) Sponges present a core prokaryotic community stable across Tropical Western Atlantic. Science of The Total Environment 835:155145. https://doi.org/10.1016/j.scitotenv.2022.155145 Leary D, Vierros M, Hamon G, et al (2009) Marine genetic resources: A review of scientific and commercial interest. Marine Policy 33:183–194. https://doi.org/10.1016/j.marpol.2008.05.010 Lemoine N, Buell N, Hill A, Hill M (2007) Assessing the utility of sponge microbial symbiont communities as models to study global climate change: a case study with Halichondria bowerbanki . In: In: Custódio MR, Lôbo-Hajdu G, Hajdu E, Muricy G (eds) Porifera Research. Biodiversity, Innovation and Sustainability. Livros de Museu Nacional 28, Rio de Janeiro, pp 419–425 Lesser MP, Sabrina Pankey M, Slattery M, et al (2022) Microbiome diversity and metabolic capacity determines the trophic ecology of the holobiont in Caribbean sponges. ISME COMMUN 2:1–12. https://doi.org/10.1038/s43705-022-00196-3 Lim SC, Putchakarn S, Thai M-Q, et al (2016) Inventory of sponge fauna from the Singapore Strait to Taiwan Strait along the western coastline of the South China Sea Lowe B, Venkatesan J, Ehrlich H, Kim S-K (2016) Global Constraints, Prospects, and Perspectives of Marine Sponge Research. In: Pallela R, Ehrlich H (eds) Marine Sponges: Chemicobiological and Biomedical Applications. Springer India, New Delhi, pp 25–35 Lu P-H, Chueh S-C, Kung F-L, et al (2007) Ilimaquinone, a marine sponge metabolite, displays anticancer activity via GADD153-mediated pathway. European Journal of Pharmacology 556:45–54. https://doi.org/10.1016/j.ejphar.2006.10.061 Maggioni F, Bell JJ, Pujo-Pay M, et al (2023a) Sponge organic matter recycling: Reduced detritus production under extreme environmental conditions. Marine Pollution Bulletin 190:114869. https://doi.org/10.1016/j.marpolbul.2023.114869 Maggioni F, Stenger P-L, Jourand P, Majorel C (2023b) The phylum Chloroflexi and their SAR202 clade dominate the microbiome of two marine sponges living in extreme environmental conditions. Marine Ecology 44:e12757. https://doi.org/10.1111/maec.12757 Malve H (2016) Exploring the ocean for new drug developments: Marine pharmacology. J Pharm Bioallied Sci 8:83–91. https://doi.org/10.4103/0975-7406.171700 Maslin M (2022) Étude de la ressource en éponge Dactylospongia metachromia pour une production durable. Theses, Université de la Polynésie française Maslin M, Gaertner-Mazouni N, Debitus C, et al (2021) Marine sponge aquaculture towards drug development: An ongoing history of technical, ecological, chemical considerations and challenges. Aquaculture Reports 21:100813. https://doi.org/10.1016/j.aqrep.2021.100813 McMurdie PJ, Holmes S (2013) phyloseq: an R package for reproducible interactive analysis and graphics of microbiome census data. PLoS One 8:e61217. https://doi.org/10.1371/journal.pone.0061217 McMurdie PJ, Holmes S (2014) Waste not, want not: why rarefying microbiome data is inadmissible. PLoS Comput Biol 10:e1003531. https://doi.org/10.1371/journal.pcbi.1003531 Mehrshad M, Rodriguez-Valera F, Amoozegar MA, et al (2018) The enigmatic SAR202 cluster up close: shedding light on a globally distributed dark ocean lineage involved in sulfur cycling. ISME J 12:655–668. https://doi.org/10.1038/s41396-017-0009-5 Meyer B, Kuever J (2008) Phylogenetic Diversity and Spatial Distribution of the Microbial Community Associated with the Caribbean Deep-water Sponge Polymastia cf. corticata by 16S rRNA, aprA, and amoA Gene Analysis. Microb Ecol 56:306–321. https://doi.org/10.1007/s00248-007-9348-5 Meyer CP, Geller JB, Paulay G (2005) Fine scale endemism on coral reefs: archipelagic differentiation in turbinid gastropods. Evolution 59:113–125 Mohamed NM, Enticknap JJ, Lohr JE, et al (2008a) Changes in bacterial communities of the marine sponge Mycale laxissima on transfer into aquaculture. Applied and Environmental Microbiology 74:1209–1222. https://doi.org/10.1128/AEM.02047-07 Mohamed NM, Rao V, Hamann MT, et al (2008b) Monitoring bacterial diversity of the marine sponge Ircinia strobilina upon transfer into aquaculture. Applied and Environmental Microbiology 74:4133–4143. https://doi.org/10.1128/AEM.00454-08 Moitinho-Silva L, Steinert G, Nielsen S, et al (2017) Predicting the HMA-LMA status in marine sponges by machine learning. Front Microbiol 8:752. https://doi.org/10.3389/fmicb.2017.00752 Monteiro LCP (2020) Rhizospheric microbial communities of plants coexisting under different edaphic conditions Morganti TM, Slaby BM, de Kluijver A, et al (2022) Giant sponge grounds of Central Arctic seamounts are associated with extinct seep life. Nat Commun 13:638. https://doi.org/10.1038/s41467-022-28129-7 Morrow KM, Fiore CL, Lesser MP (2016) Environmental drivers of microbial community shifts in the giant barrel sponge, Xestospongia muta , over a shallow to mesophotic depth gradient. Environmental Microbiology 18:2025–2038. https://doi.org/10.1111/1462-2920.13226 Neu AT, Allen EE, Roy K (2021a) Defining and quantifying the core microbiome: Challenges and prospects. Proceedings of the National Academy of Sciences 118:e2104429118. https://doi.org/10.1073/pnas.2104429118 Neu AT, Torchin ME, Allen EE, Roy K (2021b) Microbiome divergence of marine gastropod species separated by the Isthmus of Panama. 2021.07.08.451645 O’Brien PA, Tan S, Yang C, et al (2020) Diverse coral reef invertebrates exhibit patterns of phylosymbiosis. ISME J 14:2211–2222. https://doi.org/10.1038/s41396-020-0671-x Oksanen J, Blanchet FG, Friendly M, et al (2019) vegan: Community ecology package Osinga R, Sidri M, Cerig E, et al (2010) Sponge aquaculture trials in the East-Mediterranean Sea: new approaches to earlier ideas. The Open Marine Biology Journal 4: Pantile R, Webster N (2011) Strict thermal threshold identified by quantitative PCR in the sponge Rhopaloeides odorabile . Marine Ecology Progress Series 431:97–105. https://doi.org/10.3354/meps09128 Parada AE, Needham DM, Fuhrman JA (2016) Every base matters: assessing small subunit rRNA primers for marine microbiomes with mock communities, time series and global field samples. Environ Microbiol 18:1403–1414. https://doi.org/10.1111/1462-2920.13023 Pawlik JR, McMurray SE (2020) The emerging ecological and biogeochemical importance of sponges on coral reefs. Annual Review of Marine Science 12:315–337. https://doi.org/10.1146/annurev-marine-010419-010807 Pereira DM, Valentão P, Andrade PB (2013) Chap. 7 - Lessons from the Sea: Distribution, SAR, and Molecular Mechanisms of Anti-inflammatory Drugs from Marine Organisms. In: Atta-ur-Rahman (ed) Studies in Natural Products Chemistry. Elsevier, pp 205–228 Piel J, Hui D, Wen G, et al (2004) Antitumor polyketide biosynthesis by an uncultivated bacterial symbiont of the marine sponge Theonella swinhoei . Proceedings of the National Academy of Sciences 101:16222–16227. https://doi.org/10.1073/pnas.0405976101 Pita L, Erwin PM, Turon X, López-Legentil S (2013) Till Death Do Us Part: Stable Sponge-Bacteria Associations under Thermal and Food Shortage Stresses. PLOS ONE 8:e80307. https://doi.org/10.1371/journal.pone.0080307 Pita L, Rix L, Slaby BM, et al (2018) The sponge holobiont in a changing ocean: from microbes to ecosystems. Microbiome 6:46. https://doi.org/10.1186/s40168-018-0428-1 Posadas N, Baquiran JIP, Nada MAL, et al (2022) Microbiome diversity and host immune functions influence survivorship of sponge holobionts under future ocean conditions. ISME J 16:58–67. https://doi.org/10.1038/s41396-021-01050-5 Quast C, Pruesse E, Yilmaz P, et al (2013) The SILVA ribosomal RNA gene database project: improved data processing and web-based tools. Nucleic Acids Res 41:D590–D596. https://doi.org/10.1093/nar/gks1219 Radeke HS, Digits CA, Casaubon RL, Snapper ML (1999) Interactions of (–)-ilimaquinone with methylation enzymes: implications for vesicular-mediated secretion. Chemistry & Biology 6:639–647. https://doi.org/10.1016/S1074-5521(99)80115-X Risely A (2020) Applying the core microbiome to understand host–microbe systems. Journal of Animal Ecology 89:1549–1558. https://doi.org/10.1111/1365-2656.13229 Roberts DE, Davis AR, Cummins SP (2006) Experimental manipulation of shade, silt, nutrients and salinity on the temperate reef sponge Cymbastela concentrica . Marine Ecology Progress Series 307:143–154. https://doi.org/10.3354/meps307143 Ruiz C, Valderrama K, Zea S, Castellanos L (2013) Mariculture and natural production of the antitumoural (+)-discodermolide by the Caribbean marine sponge Discodermia dissoluta . Mar Biotechnol 15:571–583. https://doi.org/10.1007/s10126-013-9510-7 Schiefenhövel K, Kunzmann A (2012) Sponge farming trials: survival, attachment, and growth of two Indo-Pacific sponges, Neopetrosia sp. and Stylissa massa . Journal of Marine Sciences 2012:e417360. https://doi.org/10.1155/2012/417360 Schmitt S, Hentschel U, Taylor MW (2012a) Deep sequencing reveals diversity and community structure of complex microbiota in five Mediterranean sponges. Hydrobiologia 687:341–351. https://doi.org/10.1007/s10750-011-0799-9 Schmitt S, Tsai P, Bell J, et al (2012b) Assessing the complex sponge microbiota: core, variable and species-specific bacterial communities in marine sponges. ISME J 6:564–576. https://doi.org/10.1038/ismej.2011.116 Schuster A, Knapp I, Pomponi S, et al (2016) Validation of RAD-Seq for sponge phylogenetics Sipkema D, Franssen MCR, Osinga R, et al (2005a) Marine Sponges as Pharmacy. Mar Biotechnol 7:142–162. https://doi.org/10.1007/s10126-004-0405-5 Sipkema D, Osinga R, Schatton W, et al (2005b) Large-scale production of pharmaceuticals by marine sponges: Sea, cell, or synthesis? Biotechnology and Bioengineering 90:201–222. https://doi.org/10.1002/bit.20404 Slaby BM, Franke A, Rix L, et al (2019) Marine sponge holobionts in health and disease. In: Li Z (ed) Symbiotic Microbiomes of Coral Reefs Sponges and Corals. Springer Netherlands, Dordrecht, pp 81–104 Smith FGW (1941) Sponge disease in British Honduras, and its transmission by water currents. Ecology 22:415–421. https://doi.org/10.2307/1930719 Steffen K, Indraningrat AAG, Erngren I, et al (2022) Oceanographic setting influences the prokaryotic community and metabolome in deep-sea sponges. Sci Rep 12:3356. https://doi.org/10.1038/s41598-022-07292-3 Strand R, Whalan S, Webster NS, et al (2017) The response of a boreal deep-sea sponge holobiont to acute thermal stress. Sci Rep 7:1660. https://doi.org/10.1038/s41598-017-01091-x Swierts T, Cleary DFR, de Voogd NJ (2018) Prokaryotic communities of Indo-Pacific giant barrel sponges are more strongly influenced by geography than host phylogeny. FEMS Microbiology Ecology 94:. https://doi.org/10.1093/femsec/fiy194 Taylor MW, Hill RT, Piel J, et al (2007) Soaking it up: the complex lives of marine sponges and their microbial associates. ISME J 1:187–190. https://doi.org/10.1038/ismej.2007.32 Taylor MW, Schupp PJ, Dahllöf I, et al (2004) Host specificity in marine sponge-associated bacteria, and potential implications for marine microbial diversity. Environ Microbiol 6:121–130. https://doi.org/10.1046/j.1462-2920.2003.00545.x Ternon E, Perino E, Manconi R, et al (2017) How environmental factors affect the production of Guanidine Alkaloids by the Mediterranean sponge Crambe crambe . Marine Drugs 15:181. https://doi.org/10.3390/md15060181 Thomas T, Moitinho-Silva L, Lurgi M, et al (2016) Diversity, structure and convergent evolution of the global sponge microbiome. Nat Commun 7:11870. https://doi.org/10.1038/ncomms11870 Thoms C, Horn M, Wagner M, et al (2003) Monitoring microbial diversity and natural product profiles of the sponge Aplysina cavernicola following transplantation. Marine Biology 142:685–692. https://doi.org/10.1007/s00227-002-1000-9 Thoms C, Schupp PJ, Custódio MR, et al (2007) Chemical defense strategies in sponges: a review. Porifera research: biodiversity, Innovation and sustainability 28:627–637 Turnbaugh PJ, Ley RE, Hamady M, et al (2007) The Human Microbiome Project. Nature 449:804–810. https://doi.org/10.1038/nature06244 van Treeck P, Eisinger M, Müller J, et al (2003) Mariculture trials with Mediterranean sponge species: The exploitation of an old natural resource with sustainable and novel methods. Aquaculture 218:439–455. https://doi.org/10.1016/S0044-8486(03)00010-3 Vavrek M (2011) Fossil: Palaeoecological and Palaeogeographical Analysis Tools. Palaeontologia Electronica 14:16 Voultsiadou E (2007) Sponges: an historical survey of their knowledge in Greek antiquity. Journal of the Marine Biological Association of the United Kingdom 87:1757–1763. https://doi.org/10.1017/S0025315407057773 Webster NS, Cobb RE, Negri AP (2008) Temperature thresholds for bacterial symbiosis with a sponge. ISME J 2:830–842. https://doi.org/10.1038/ismej.2008.42 Weigel BL, Erwin PM (2017) Effects of reciprocal transplantation on the microbiome and putative nitrogen cycling functions of the intertidal sponge, Hymeniacidon heliophila . Sci Rep 7:43247. https://doi.org/10.1038/srep43247 Wilson MC, Mori T, Rückert C, et al (2014) An environmental bacterial taxon with a large and distinct metabolic repertoire. Nature 506:58–62. https://doi.org/10.1038/nature12959 World Meteorological Association (2022) State of the Global Climate 2021 Wulff J (2001) Assessing and monitoring coral reef sponges: Why and how? Bulletin of Marine Science 69:831–846 Yang Q, Franco CMM, Lin H-W, Zhang W (2019) Untapped sponge microbiomes: structure specificity at host order and family levels. FEMS Microbiology Ecology 95:fiz136. https://doi.org/10.1093/femsec/fiz136 Yasuhara-Bell J, Lu Y (2010) Marine compounds and their antiviral activities. Antiviral Research 86:231–240. https://doi.org/10.1016/j.antiviral.2010.03.009 Zaneveld JR, McMinds R, Vega Thurber R (2017) Stress and stability: applying the Anna Karenina principle to animal microbiomes. Nat Microbiol 2:1–8. https://doi.org/10.1038/nmicrobiol.2017.121 Additional Declarations No competing interests reported. Supplementary Files SIDactylospongiamicrobiomeAntonievanLeeuwenhoek.docx Cite Share Download PDF Status: Published Journal Publication published 11 Apr, 2024 Read the published version in Antonie van Leeuwenhoek → Version 1 posted Editorial decision: Revision requested 25 Mar, 2024 Reviews received at journal 24 Mar, 2024 Reviewers agreed at journal 04 Mar, 2024 Reviewers invited by journal 04 Mar, 2024 Submission checks completed at journal 02 Mar, 2024 Editor assigned by journal 02 Mar, 2024 First submitted to journal 01 Mar, 2024 You are reading this latest preprint version Research Square lets you share your work early, gain feedback from the community, and start making changes to your manuscript prior to peer review in a journal. As a division of Research Square Company, we’re committed to making research communication faster, fairer, and more useful. We do this by developing innovative software and high quality services for the global research community. Our growing team is made up of researchers and industry professionals working together to solve the most critical problems facing scientific publishing. Also discoverable on Platform About Our Team In Review Editorial Policies Advisory Board Help Center Resources Author Services Accessibility API Access RSS feed Manage Cookie Preferences © Research Square 2026 | ISSN 2693-5015 (online) Privacy Policy Terms of Service Do Not Sell My Personal Information {"props":{"pageProps":{"initialData":{"identity":"rs-4002744","acceptedTermsAndConditions":true,"allowDirectSubmit":false,"archivedVersions":[],"articleType":"Research Article","associatedPublications":[],"authors":[{"id":276378091,"identity":"3a1327d5-6dc1-441a-990a-12f32f1d7493","order_by":0,"name":"Mathilde Maslin","email":"","orcid":"","institution":"Univ. Polynesie Française","correspondingAuthor":false,"prefix":"","firstName":"Mathilde","middleName":"","lastName":"Maslin","suffix":""},{"id":276378092,"identity":"eeb38bae-613f-4ab2-9993-3ada76d66640","order_by":1,"name":"Benoît Paix","email":"data:image/png;base64,iVBORw0KGgoAAAANSUhEUgAAAZAAAAAyAQMAAABI0h/eAAAABlBMVEX///8AAABVwtN+AAAACXBIWXMAAA7EAAAOxAGVKw4bAAAA50lEQVRIiWNgGAWjYBACAxCRACKYGRgfACkePuK0JIC1MIM4PGxEaYFYw8AmASYJaTFn7z344OEPhjz+dh6zyq85djJsDMwPH93Ao8Wy51yyAdBhxRKH2dJuy25LBjqMzdg4B5/DbuSYSQC1JDYcZj52W3IbM1ALD5s0Xi3330C0zD/M2FYsua2eCC03eCBaNgBtYfy47TBhLZY9OcYGCWkSxYaH2ZKlGbcd52FjJuAXc/Yzhg9/2NjkyZ0/Y/jx57Zqe3725oeP8WmBAokEEMnMAyYJKwcDsBbGH0SqHgWjYBSMgpEFAJsuQQZzR/ntAAAAAElFTkSuQmCC","orcid":"","institution":"Naturalis Biodiversity Center","correspondingAuthor":true,"prefix":"","firstName":"Benoît","middleName":"","lastName":"Paix","suffix":""},{"id":276378093,"identity":"b03b0fb8-3df0-47f3-8233-c42d4b9746f2","order_by":2,"name":"Niels Windt","email":"","orcid":"","institution":"Naturalis Biodiversity Center","correspondingAuthor":false,"prefix":"","firstName":"Niels","middleName":"","lastName":"Windt","suffix":""},{"id":276378094,"identity":"304d8d58-6569-49b1-9948-43b575f416d9","order_by":3,"name":"Rohani Ambo-Rappe","email":"","orcid":"","institution":"Hasanuddin University","correspondingAuthor":false,"prefix":"","firstName":"Rohani","middleName":"","lastName":"Ambo-Rappe","suffix":""},{"id":276378095,"identity":"2a8a02a3-1e05-44cf-93eb-0e507c1f655f","order_by":4,"name":"Cécile Debitus","email":"","orcid":"","institution":"IRD, Univ Brest, CNRS, Ifremer, LEMAR","correspondingAuthor":false,"prefix":"","firstName":"Cécile","middleName":"","lastName":"Debitus","suffix":""},{"id":276378096,"identity":"c99ef94d-b833-400d-a069-b48576f71e8c","order_by":5,"name":"Nabila Gaertner-Mazouni","email":"","orcid":"","institution":"Univ. Polynesie Française","correspondingAuthor":false,"prefix":"","firstName":"Nabila","middleName":"","lastName":"Gaertner-Mazouni","suffix":""},{"id":276378097,"identity":"c955982a-9afc-4f67-a495-a54d83c84ad8","order_by":6,"name":"Raimana Ho","email":"","orcid":"","institution":"Univ. Polynesie Française","correspondingAuthor":false,"prefix":"","firstName":"Raimana","middleName":"","lastName":"Ho","suffix":""},{"id":276378098,"identity":"cd4b95dc-fe66-4ac7-a247-7b74f268a332","order_by":7,"name":"Nicole J. Voogd","email":"","orcid":"","institution":"Naturalis Biodiversity Center","correspondingAuthor":false,"prefix":"","firstName":"Nicole","middleName":"J.","lastName":"Voogd","suffix":""}],"badges":[],"createdAt":"2024-03-01 09:25:09","currentVersionCode":1,"declarations":"","doi":"10.21203/rs.3.rs-4002744/v1","doiUrl":"https://doi.org/10.21203/rs.3.rs-4002744/v1","draftVersion":[],"editorialEvents":[{"content":"https://doi.org/10.1007/s10482-024-01962-0","type":"published","date":"2024-04-11T15:01:50+00:00"}],"editorialNote":"","failedWorkflow":false,"files":[{"id":52304567,"identity":"c505e5af-829a-4ab0-8bf8-eb1cb229eb29","added_by":"auto","created_at":"2024-03-08 19:05:16","extension":"jpg","order_by":1,"title":"Figure 1","display":"","copyAsset":false,"role":"figure","size":650223,"visible":true,"origin":"","legend":"\u003cp\u003eLocations of the sampling sites of the study. \u003cstrong\u003eA\u003c/strong\u003e: map of all sampling sites (Mercator projection). \u003cstrong\u003eB\u003c/strong\u003e: map of the six French Polynesian islands of the study (Tetiaroa, Rangiroa, Raroia, Makemo Tematangi and Mangareva). \u003cstrong\u003eC\u003c/strong\u003e: map of sampling sites located in Rangiroa atoll (Avatoru, Tiputa, Lagon 1 and Lagon 2). Farming trials were conducted on the Western reef side of Rangiroa, roughly 1km away from the Avatoru channel.\u003c/p\u003e","description":"","filename":"Figure1.jpg","url":"https://assets-eu.researchsquare.com/files/rs-4002744/v1/927633d4d2597c15ac1ddfc8.jpg"},{"id":52304563,"identity":"b059efbd-25bf-4a21-818e-20e857512072","added_by":"auto","created_at":"2024-03-08 19:05:16","extension":"jpg","order_by":2,"title":"Figure 2","display":"","copyAsset":false,"role":"figure","size":1212809,"visible":true,"origin":"","legend":"\u003cp\u003eDiversity of the prokaryotic communities associated with \u003cem\u003eDactylospongia metachromia\u003c/em\u003esamples collected during the temporal study. \u003cstrong\u003eA\u003c/strong\u003e, \u003cstrong\u003eB\u003c/strong\u003e and \u003cstrong\u003eC\u003c/strong\u003e: boxplots representing the \u003cem\u003eα\u003c/em\u003e-diversity metrics (\u003cstrong\u003eA\u003c/strong\u003e for Shannon, \u003cstrong\u003eB\u003c/strong\u003e for Chao1 and \u003cstrong\u003eC\u003c/strong\u003e for Pielou indices). \u003cstrong\u003eD\u003c/strong\u003e and \u003cstrong\u003eE\u003c/strong\u003e: NMDS plots representing the \u003cem\u003eβ\u003c/em\u003e-diversity (Bray-Curtis distances). \u003cstrong\u003eF\u003c/strong\u003e and \u003cstrong\u003eG\u003c/strong\u003e: db-RDA plots representing the environmental parameters as explanatory variables of the \u003cem\u003eβ\u003c/em\u003e-diversity (Bray-Curtis distances). \u003cstrong\u003eD\u003c/strong\u003e and \u003cstrong\u003eF\u003c/strong\u003e: plots of the sample scores of the NMDS and db-RDA, respectively. \u003cstrong\u003eE\u003c/strong\u003e and \u003cstrong\u003eG\u003c/strong\u003e: plots of the ASVs scores of the NMDS and db-RDA, respectively.\u003c/p\u003e","description":"","filename":"Figure2.jpg","url":"https://assets-eu.researchsquare.com/files/rs-4002744/v1/530b97e67ea3a96d1017284f.jpg"},{"id":52304798,"identity":"569e8158-a071-485b-8fcc-7680a7c5d077","added_by":"auto","created_at":"2024-03-08 19:13:16","extension":"jpg","order_by":3,"title":"Figure 3","display":"","copyAsset":false,"role":"figure","size":703876,"visible":true,"origin":"","legend":"\u003cp\u003e\u003cem\u003eβ\u003c/em\u003e-diversity of the prokaryotic communities associated with \u003cem\u003eDactylospongia metachromia\u003c/em\u003e samples from the biogeographical study. \u003cstrong\u003eA\u003c/strong\u003e and \u003cstrong\u003eB\u003c/strong\u003e: NMDS plot. \u003cstrong\u003eC\u003c/strong\u003e and \u003cstrong\u003eD\u003c/strong\u003e: db-RDA plots representing the environmental parameters as explanatory variables of the \u003cem\u003eβ\u003c/em\u003e-diversity. \u003cstrong\u003eA\u003c/strong\u003e and \u003cstrong\u003eC\u003c/strong\u003e: plots of the sample scores of the NMDS and db-RDA, respectively. \u003cstrong\u003eB\u003c/strong\u003e and \u003cstrong\u003eD\u003c/strong\u003e: plots of the ASVs scores of the NMDS and db-RDA, respectively\u003c/p\u003e","description":"","filename":"Figure3.jpg","url":"https://assets-eu.researchsquare.com/files/rs-4002744/v1/3c88906dc44cc43bc5d5ebc0.jpg"},{"id":52304566,"identity":"1ea8b505-de75-4787-a467-2f305d811f6b","added_by":"auto","created_at":"2024-03-08 19:05:16","extension":"jpg","order_by":4,"title":"Figure 4","display":"","copyAsset":false,"role":"figure","size":700122,"visible":true,"origin":"","legend":"\u003cp\u003eBarplots of the relative percentages of the core prokaryotic community composition (at the family level) of \u003cem\u003eDactylospongia metachromia\u003c/em\u003e. Core ASVs are present at least in 90% of the samples. \u003cstrong\u003eA\u003c/strong\u003e: samples from the farming trials (Donor: samples collected in December 2019; 9-, 12-, 15- and 18-months samples: cultured samples respectively collected in September 2020, December 2020, March 2021 and May 2021). \u003cstrong\u003eB:\u003c/strong\u003esamples from the biogeographical study (abbreviations used: Mal: Maldives, Mak: Makemo).\u003c/p\u003e","description":"","filename":"Figure4.jpg","url":"https://assets-eu.researchsquare.com/files/rs-4002744/v1/e2cd9727a47654d169495ddf.jpg"},{"id":54713169,"identity":"3fa4945c-ce69-4189-b86a-09a93f7a8647","added_by":"auto","created_at":"2024-04-15 15:14:51","extension":"pdf","order_by":0,"title":"","display":"","copyAsset":false,"role":"manuscript-pdf","size":1336128,"visible":true,"origin":"","legend":"","description":"","filename":"manuscript.pdf","url":"https://assets-eu.researchsquare.com/files/rs-4002744/v1/8c1d9c1a-1ee4-43c0-ae0d-b81feb932f2b.pdf"},{"id":52304564,"identity":"e7ee4642-76c6-4eff-ba2c-daea4b609699","added_by":"auto","created_at":"2024-03-08 19:05:16","extension":"docx","order_by":6,"title":"","display":"","copyAsset":false,"role":"supplement","size":7024349,"visible":true,"origin":"","legend":"","description":"","filename":"SIDactylospongiamicrobiomeAntonievanLeeuwenhoek.docx","url":"https://assets-eu.researchsquare.com/files/rs-4002744/v1/09ceb043579e1b46940645cd.docx"}],"financialInterests":"No competing interests reported.","formattedTitle":"Prokaryotic communities of the French Polynesian sponge Dactylospongia metachromia display a site-specific and stable diversity during an aquaculture trial","fulltext":[{"header":"Introduction","content":"\u003cp\u003eSponges are ubiquitous filter-feeders involved in key ecological processes such as nutrient cycling, energy transfer towards higher trophic levels, substratum stability, and habitat for a wide range of associated marine organisms (Bell \u003cspan citationid=\"CR9\" class=\"CitationRef\"\u003e2008\u003c/span\u003e; Pawlik and McMurray \u003cspan citationid=\"CR96\" class=\"CitationRef\"\u003e2020\u003c/span\u003e). They have been farmed by humans since Greek antiquity, alternatively used as hygienic tools, household devices, medical instruments, artistic material, and biopolymers (Ehrlich and Worch \u003cspan citationid=\"CR39\" class=\"CitationRef\"\u003e2007\u003c/span\u003e; Voultsiadou \u003cspan citationid=\"CR127\" class=\"CitationRef\"\u003e2007\u003c/span\u003e). Over the past decades, sponges have been increasingly studied as prolific sources of nature-based pharmaceuticals (Sipkema et al. \u003cspan citationid=\"CR111\" class=\"CitationRef\"\u003e2005a\u003c/span\u003e; Lowe et al. \u003cspan citationid=\"CR71\" class=\"CitationRef\"\u003e2016\u003c/span\u003e; Maslin et al. \u003cspan citationid=\"CR77\" class=\"CitationRef\"\u003e2021\u003c/span\u003e). For instance, sponge farming programs are an important practice to grow enough biomass to extract their natural products of interest while conserving natural stocks (Duckworth \u003cspan citationid=\"CR34\" class=\"CitationRef\"\u003e2009\u003c/span\u003e; Bierwirth et al. \u003cspan citationid=\"CR11\" class=\"CitationRef\"\u003e2022\u003c/span\u003e). To this end, \u003cem\u003ein situ\u003c/em\u003e and \u003cem\u003eex situ\u003c/em\u003e trials have been conducted in a wide range of geographical zones, for example with \u003cem\u003eNeopetrosia\u003c/em\u003e sp., \u003cem\u003eStylissa massa\u003c/em\u003e and \u003cem\u003eCallyspongia biru\u003c/em\u003e in the Indo-Pacific region (de Voogd \u003cspan citationid=\"CR32\" class=\"CitationRef\"\u003e2007\u003c/span\u003e; Schiefenh\u0026ouml;vel and Kunzmann \u003cspan citationid=\"CR107\" class=\"CitationRef\"\u003e2012\u003c/span\u003e), \u003cem\u003eAxinella damicornis\u003c/em\u003e, \u003cem\u003eIrcinia variabilis\u003c/em\u003e, \u003cem\u003eChondrosia reniformis\u003c/em\u003e, and \u003cem\u003eCrambe crambe\u003c/em\u003e in the Mediterranean Sea (van Treeck et al. \u003cspan citationid=\"CR125\" class=\"CitationRef\"\u003e2003\u003c/span\u003e; Ternon et al. \u003cspan citationid=\"CR120\" class=\"CitationRef\"\u003e2017\u003c/span\u003e), and \u003cem\u003eDiscodermia dissoluta\u003c/em\u003e in the Caribbean (Ruiz et al. \u003cspan citationid=\"CR106\" class=\"CitationRef\"\u003e2013\u003c/span\u003e). Increased specialized metabolite yields were for instance observed in \u003cem\u003eNegombata magnifica\u003c/em\u003e from the Red Sea (Hadas et al. \u003cspan citationid=\"CR47\" class=\"CitationRef\"\u003e2005\u003c/span\u003e), \u003cem\u003eDysidea avara\u003c/em\u003e from the Mediterranean (De Caralt et al. \u003cspan citationid=\"CR30\" class=\"CitationRef\"\u003e2010\u003c/span\u003e), as well as some New Zealand species (Duckworth and Battershill \u003cspan citationid=\"CR35\" class=\"CitationRef\"\u003e2003a\u003c/span\u003e, \u003cspan citationid=\"CR36\" class=\"CitationRef\"\u003eb\u003c/span\u003e), leading to important economical perspectives (Sipkema et al. \u003cspan citationid=\"CR111\" class=\"CitationRef\"\u003e2005a\u003c/span\u003e; Osinga et al. \u003cspan citationid=\"CR93\" class=\"CitationRef\"\u003e2010\u003c/span\u003e).\u003c/p\u003e \u003cp\u003eSponges host a large diversity of microorganisms, which are often described as host-specific (Thomas et al. \u003cspan citationid=\"CR121\" class=\"CitationRef\"\u003e2016\u003c/span\u003e). Together, they form a functional entity named the \u0026ldquo;sponge holobiont\u0026rdquo;, whereby both host and microbiome display tight interactions contributing to their overall functioning (Pita et al. \u003cspan citationid=\"CR100\" class=\"CitationRef\"\u003e2018\u003c/span\u003e; Slaby et al. \u003cspan citationid=\"CR113\" class=\"CitationRef\"\u003e2019\u003c/span\u003e). The sponge holobiont concept has allowed to better understand the role of sponge prokaryotic communities in the production of various biosynthetic pathways (Wilson et al. \u003cspan citationid=\"CR130\" class=\"CitationRef\"\u003e2014\u003c/span\u003e; Harvey et al. \u003cspan citationid=\"CR50\" class=\"CitationRef\"\u003e2015\u003c/span\u003e). As an example, sponge-associated bacteria can contribute to the production of bioactive compounds (Piel et al. \u003cspan citationid=\"CR98\" class=\"CitationRef\"\u003e2004\u003c/span\u003e; Indraningrat et al. \u003cspan citationid=\"CR59\" class=\"CitationRef\"\u003e2016\u003c/span\u003e), with chemical properties acting against biofoulers (Aguila-Ram\u0026iacute;rez et al. \u003cspan citationid=\"CR1\" class=\"CitationRef\"\u003e2014\u003c/span\u003e), predators (e.g. fishes; Thoms et al. \u003cspan citationid=\"CR123\" class=\"CitationRef\"\u003e2007\u003c/span\u003e; Slaby et al. \u003cspan citationid=\"CR113\" class=\"CitationRef\"\u003e2019\u003c/span\u003e) or competitors (e.g. corals competing for space; Thoms et al. \u003cspan citationid=\"CR123\" class=\"CitationRef\"\u003e2007\u003c/span\u003e). Sponge microbiomes may also contribute to the nutrient cycling and energy transfer roles of their host in their environment (de Goeij et al. \u003cspan citationid=\"CR31\" class=\"CitationRef\"\u003e2017\u003c/span\u003e; Pita et al. \u003cspan citationid=\"CR100\" class=\"CitationRef\"\u003e2018\u003c/span\u003e). In a global context of environmental changes such as ocean acidification, immune functions associated with the microbiome can influence the survivorship and the resilience of the sponge holobiont to future climate scenarios (Posadas et al. \u003cspan citationid=\"CR101\" class=\"CitationRef\"\u003e2022\u003c/span\u003e). Considering these key roles, the dynamics of the microbiome diversity appears as an essential factor to be monitored during transplantation and farming programs (Mohamed et al. \u003cspan citationid=\"CR83\" class=\"CitationRef\"\u003e2008a\u003c/span\u003e, \u003cspan citationid=\"CR84\" class=\"CitationRef\"\u003eb\u003c/span\u003e; Thomas et al. \u003cspan citationid=\"CR121\" class=\"CitationRef\"\u003e2016\u003c/span\u003e). Such investigations provide a better holistic understanding of the parameters linked to the sponge health and thus, to the farming performance itself under a context of environmental pressures (Pita et al. \u003cspan citationid=\"CR100\" class=\"CitationRef\"\u003e2018\u003c/span\u003e). To date, these dynamics have been monitored only under few aquaculture trials. Alpha-diversity increases of bacterial communities from \u003cem\u003eMycale laximissa\u003c/em\u003e and \u003cem\u003eIrcinia strobilina\u003c/em\u003e (both collected from the Gulf of Mexico) were observed under closed aquaculture systems (Mohamed et al. \u003cspan citationid=\"CR83\" class=\"CitationRef\"\u003e2008a\u003c/span\u003e, \u003cspan citationid=\"CR84\" class=\"CitationRef\"\u003eb\u003c/span\u003e). Another example is related to the transplantation of \u003cem\u003eGelliodes obtusa\u003c/em\u003e to fish farms in Philippines and revealed a stable bacterial community while facing eutrophic pressures (e.g. high concentrations of ammonium and phosphate; Baquiran and Conaco \u003cspan citationid=\"CR6\" class=\"CitationRef\"\u003e2018\u003c/span\u003e).\u003c/p\u003e \u003cp\u003eThe Indo-Pacific sponge \u003cem\u003eDactylospongia metachromia\u003c/em\u003e (Demospongiae: Dictyoceratida: Thorectidae) was part of a farming program for the extraction of several bioactive compounds in French Polynesia (Maslin \u003cspan citationid=\"CR76\" class=\"CitationRef\"\u003e2022\u003c/span\u003e). This sponge species is known to produce two secondary metabolites of interest named ilimaquinone (Capon and MacLeod \u003cspan citationid=\"CR20\" class=\"CitationRef\"\u003e1987\u003c/span\u003e) and 5-epi-ilimaquinone (Carte et al. \u003cspan citationid=\"CR22\" class=\"CitationRef\"\u003e1985\u003c/span\u003e). Promising pharmaceutical properties were evidenced for these two compounds, including anticancer, antimicrobial, anti-HIV, antiparasitic, anti-diabetes and anti-inflammatory properties (Radeke et al. \u003cspan citationid=\"CR103\" class=\"CitationRef\"\u003e1999\u003c/span\u003e; Lu et al. \u003cspan citationid=\"CR72\" class=\"CitationRef\"\u003e2007\u003c/span\u003e; Kochanowska et al. \u003cspan citationid=\"CR62\" class=\"CitationRef\"\u003e2008\u003c/span\u003e; Leary et al. \u003cspan citationid=\"CR67\" class=\"CitationRef\"\u003e2009\u003c/span\u003e; Yasuhara-Bell and Lu \u003cspan citationid=\"CR134\" class=\"CitationRef\"\u003e2010\u003c/span\u003e; Pereira et al. \u003cspan citationid=\"CR97\" class=\"CitationRef\"\u003e2013\u003c/span\u003e; Daletos et al. \u003cspan citationid=\"CR28\" class=\"CitationRef\"\u003e2014\u003c/span\u003e; Hitora et al. \u003cspan citationid=\"CR53\" class=\"CitationRef\"\u003e2021\u003c/span\u003e). This species has a wide distribution range from the Red Sea to French Polynesia (Hooper et al. \u003cspan citationid=\"CR54\" class=\"CitationRef\"\u003e2002\u003c/span\u003e; Lim et al. \u003cspan citationid=\"CR70\" class=\"CitationRef\"\u003e2016\u003c/span\u003e; Cleary et al. \u003cspan citationid=\"CR26\" class=\"CitationRef\"\u003e2020\u003c/span\u003e), suggesting an important adaptation potential to different environments. However, the microbial diversity associated with \u003cem\u003eD. metachormia\u003c/em\u003e is still underexplored to date. Using a denaturing gradient gel electrophoresis fingerprinting approach, Jeong and Park (\u003cspan citationid=\"CR60\" class=\"CitationRef\"\u003e2013\u003c/span\u003e) described the bacterial community from Micronesian \u003cem\u003eD. metachromia\u003c/em\u003e samples, while a Next Generation Sequencing method was undertaken by Cleary et al. (\u003cspan citationid=\"CR26\" class=\"CitationRef\"\u003e2020\u003c/span\u003e) for two \u003cem\u003eD. metachromia\u003c/em\u003e samples from the Red Sea. In both studies, the most represented bacterial phyla were Actinobacteria, Proteobacteria, Chloroflexi, Acidobacteria, and Gemmatimmonadetes which are commonly abundant in so-called high microbial abundance (HMA) sponges. These sponges harbor prokaryotic communities with higher diversities compared to Low Microbial Abundance (LMA) sponges (Moitinho-Silva et al. \u003cspan citationid=\"CR85\" class=\"CitationRef\"\u003e2017\u003c/span\u003e).\u003c/p\u003e \u003cp\u003eAs the temporal dynamics of \u003cem\u003eD. metachromia\u003c/em\u003e microbiome have never been investigated before, the present study aimed to monitor its diversity during the farming experiment conducted in French Polynesia (Maslin \u003cspan citationid=\"CR76\" class=\"CitationRef\"\u003e2022\u003c/span\u003e). Subsequently, the objective was also to identify the major environmental factors (i.e. temperature, salinity, chlorophyll \u003cem\u003ea\u003c/em\u003e and several nutrients) shaping the potentially-observed microbial changes. Considering the good survival and growth rates of the sponge host observed during trials in sites where it naturally occurs, we hypothesized a temporally stable microbiome. Additionally, a multi-scale biogeographical approach was conducted to have a broader understanding of the relative importance of these abiotic parameters, as well as the role of the host phylogeny. The core community associated with \u003cem\u003eD. metachromia\u003c/em\u003e was also analyzed in order to look for essential sponge-related symbionts that could either be shaped by such variables or be broadly stable over time and space.\u003c/p\u003e"},{"header":"Materials and Methods","content":"\u003cdiv id=\"Sec3\" class=\"Section2\"\u003e \u003ch2\u003eSampling and biological materials from the farming trials\u003c/h2\u003e \u003cp\u003eThe farming trials were conducted at Rangiroa, an atoll from the Tuamotu archipelago. The reef site named \u0026ldquo;Avatoru\u0026rdquo; (Fig.\u0026nbsp;\u003cspan refid=\"Fig1\" class=\"InternalRef\"\u003e1\u003c/span\u003e) was chosen, in consultation with the maritime authorities, to perform the aquaculture trials since \u003cem\u003eD. metachromia\u003c/em\u003e was found naturally growing in that area. The sponge exploitation remained under the regulations of a permanent operating license from the Direction des Ressources Marines (DRM, Tahiti) during the whole study.\u003c/p\u003e \u003cp\u003eThe aquaculture frames were directly anchored to the reef as supports of sponge explants coming from nearby parental specimens. These structures consisted of three horizontal, table-shaped iron frames submerged between 16 and 18m depth where the abundance of wild \u003cem\u003eD. metachromia\u003c/em\u003e populations was found the highest (\u003cb\u003eFigure \u003cspan refid=\"MOESM1\" class=\"InternalRef\"\u003eS1\u003c/span\u003e\u003c/b\u003e). Four different sponge explants (roughly 4*4*4cm) were excised in December 2019 from the same original sponge. An additional piece was directly sampled as an uncultivated control, named \u0026ldquo;Donor\u0026rdquo;. The four explants were attached to the same frame using a nylon rope directly threaded through their spongin skeleton (\u003cb\u003eFigure S2\u003c/b\u003e). Their sampling was then carried out using SCUBA at four distinct times, spaced by 3 months. During every field survey, survival and growth rates of the new sponge generations were calculated and one explant originally coming from the same donor sponge was removed to be subsampled in triplicates. To do so, pieces of \u003cem\u003eD. metachromia\u003c/em\u003e specimens were manually excised with a knife to prevent squeezing the animals and transferred into a plastic bag filled with the surrounding seawater. Only small fragments of about 2*2cm were removed from each selected individual, leaving sufficient pinacoderm (i.e. outer layer of sponge cells) intact to allow regeneration after damage.\u003c/p\u003e \u003cp\u003eSamples were named according to their farming duration as follows: \u0026ldquo;9 months\u0026rdquo; for those collected in September 2020, \u0026ldquo;12 months\u0026rdquo; for December 2020, \u0026ldquo;15 months\u0026rdquo; for March 2021 and \u0026ldquo;18 months\u0026rdquo; for May 2021 (\u003cb\u003eTable \u003cspan refid=\"MOESM1\" class=\"InternalRef\"\u003eS1\u003c/span\u003e\u003c/b\u003e). Subsamples were cut using sterilized tweezers and scalpels for the barcoding analysis of the sponge and the metabarcoding study of the prokaryotic community. They were subsequently preserved in sterilized plastic vials filled with 96% ethanol and conserved at -20\u0026deg;C until DNA extraction.\u003c/p\u003e \u003cp\u003e \u003c/p\u003e \u003c/div\u003e \u003cdiv id=\"Sec4\" class=\"Section2\"\u003e \u003ch2\u003eSampling and biological materials for the biogeographical complementary study\u003c/h2\u003e \u003cp\u003eDifferent geographical scales were considered for the biogeographical study. A first approach was to conduct a global comparison, with samples collected from four different regions: the Red Sea (Saudi Arabia), the Indian Ocean (Maldives), the Makassar Strait (Indonesia), and the South Pacific Ocean (French Polynesia) at different sampling times from 2014 to 2019 (Fig.\u0026nbsp;\u003cspan refid=\"Fig1\" class=\"InternalRef\"\u003e1\u003c/span\u003eA, \u003cb\u003eTable S2\u003c/b\u003e). Then, a regional scale was considered by targeting three different French Polynesian archipelagos.\u003c/p\u003e \u003cp\u003eThe sampling of French Polynesian sponges was organized during an oceanographic campaign led by the French \u0026ldquo;Institut de Recherche pour le D\u0026eacute;veloppement\u0026rdquo; (IRD) from 19th October to 23th November 2018 onboard the \u003cem\u003eAlis\u003c/em\u003e oceanographic vessel. Subsamples were taken from donor sponge populations located at 20m depth on the barrier reef of six different islands, encompassing more than 4000km\u003csup\u003e2\u003c/sup\u003e and covering three neighboring archipelagos: Rangiroa, Makemo, Raroia, Tematangi (Tuamotu), Mangareva (Gambier) and Tetiaroa (Society) (Fig.\u0026nbsp;\u003cspan refid=\"Fig1\" class=\"InternalRef\"\u003e1\u003c/span\u003eB, \u003cb\u003eTable S2 and S3\u003c/b\u003e).\u003c/p\u003e \u003cp\u003eTo analyze differences at a higher geographical resolution, considering differences in terms of prevalent currents and oceanic exposition, four specific sites were chosen within the atoll of Rangiroa (Fig.\u0026nbsp;\u003cspan refid=\"Fig1\" class=\"InternalRef\"\u003e1\u003c/span\u003eC). Two locations named \u0026ldquo;Tiputa\u0026rdquo; and \u0026ldquo;Avatoru\u0026rdquo; were located on the external reef on the Western sides of both Tiputa and Avatoru passes. The two remaining areas named \u0026ldquo;Lagon 1\u0026rdquo; and \u0026ldquo;Lagon 2\u0026rdquo; were located within the lagoon, both facing the Avatoru pass in order to benefit from the semi-oceanic conditions allowed by the daily water exchange.\u003c/p\u003e \u003cp\u003eThe sampling and subsampling were then processed as previously described for the study of farming trials.\u003c/p\u003e \u003c/div\u003e \u003cdiv id=\"Sec5\" class=\"Section2\"\u003e \u003ch2\u003eEnvironmental parameters\u003c/h2\u003e \u003cp\u003eDifferent devices were used to record temperature and salinity conditions throughout the farming experiment in Rangiroa. First, a temperature data logger (iBee 22L Thermo Button, Plug and Track, Proges Plus, France) was directly attached to one of the frames supporting the explants, to record seasonal temperature variations every hour with an accuracy of 0,1\u0026deg;C. The ThermoTrack PC V8 software was used for device calibration and data processing. Then, for salinity measurements in Practical Salinity Unit (PSU), the AquaSonde-5000 submersible multi-parameter probe (Aquaread, SDEC, France) was deployed using SCUBA at the beginning of each field survey. It was previously calibrated in the laboratory according to the manufacturer's recommendations and using the SondeLink-PC software. Lastly, for the remaining French Polynesian islands investigated in this study except Makemo due to a lack of working probes, temperature and salinity were automatically registered during each stopover using factory sensors directly operated onboard the \u003cem\u003eAlis\u003c/em\u003e vessel.\u003c/p\u003e \u003cp\u003eFor nutrient analyses, a subsample of 20mL of seawater was fixed with 20\u0026micro;L of HgCl\u003csub\u003e2\u003c/sub\u003e in polyethylene scintillation vials (Dominique Dutscher SAS). The vials were vigorously shaken, and then kept at 4\u0026deg;C and away from light exposure. Nutrient content analyses of dissolved phosphate (PO\u003csub\u003e4\u003c/sub\u003e\u003csup\u003e3\u0026minus;\u003c/sup\u003e), nitrogen oxides (NO\u003csub\u003ex\u003c/sub\u003e), and orthosilicic acid (Si[OH]\u003csub\u003e4\u003c/sub\u003e) were performed by the Laboratoire des Moyens Analytiques (LAMA, Noum\u0026eacute;a, New Caledonia) using a SEAL AA3 HR AutoAnalyzer.\u003c/p\u003e \u003cp\u003eFor chlorophyll \u003cem\u003ea\u003c/em\u003e measurement, a subsample of 500 mL of seawater was filtered on a Whatman\u0026reg; GF/F glass-fiber filter (0,7 nominal pore size; 47 mm diameter) to prevent the loss of picoplankton (Aminot and Rey \u003cspan citationid=\"CR2\" class=\"CitationRef\"\u003e2001\u003c/span\u003e). Filters were then placed in aluminum foils and stored at -20\u0026deg;C. Extraction of chlorophyll \u003cem\u003ea\u003c/em\u003e from phytoplanktonic cells was performed using 90% acetone after mechanical grinding of the filters into glass tubes. The 6-mL extracts were then stored overnight in darkness at 4\u0026deg;C and centrifuged the next day a few hours prior to analysis. Concentrations were measured by a Trilogy\u0026reg; Laboratory Fluorometer (Turner Designs) using the chlorophyll \u003cem\u003ea\u003c/em\u003e Non-Acidification module. Calibration of the fluorometer was made using a Prim\u0026rsquo;Light SECOMAM\u0026reg; spectrophotometer.\u003c/p\u003e \u003c/div\u003e \u003cdiv id=\"Sec6\" class=\"Section2\"\u003e \u003ch2\u003eBarcoding of sponge samples\u003c/h2\u003e \u003cp\u003eFor phylogenetic analyses, DNA was extracted from all samples (except those from Lagon 1 and Lagon 2) using the Qiagen DNeasy Blood \u0026amp; Tissue kit and following manufacturer protocols for spin column extractions. The sponge 28S rRNA gene was amplified using the primers 28S-C2-fwd (5\u0026prime;-GAAAAGAACTTTGRARAGAGAGT-3\u0026prime;) and 28S-D2-rev (5\u0026prime;-TCCGTGTTTCAAGACGGG-3\u0026prime;) (Chombard et al. \u003cspan citationid=\"CR24\" class=\"CitationRef\"\u003e1998\u003c/span\u003e) in a 25 \u0026micro;l reaction containing: 7.9\u0026micro;l distilled water (sterile milliQ); 5\u0026micro;l 5x Thermo Fisher Phire Green; 1.0\u0026micro;l 25mM MgCl\u003csub\u003e2\u003c/sub\u003e; 1.0\u0026micro;l 10mg.ml\u003csup\u003e\u0026minus;\u0026thinsp;1\u003c/sup\u003e Life BSA; 5.0\u0026micro;l 5x Qiagen Q-solution; 1.3\u0026micro;l 10 pMol.\u0026micro;l\u003csup\u003e\u0026minus;\u0026thinsp;1\u003c/sup\u003e forward primer; 1.3\u0026micro;l 10 pMol.\u0026micro;l\u003csup\u003e\u0026minus;\u0026thinsp;1\u003c/sup\u003e reverse primer; 0.5\u0026micro;l 2.5mM dNTP; 0.5\u0026micro;l Thermo Fisher Phire II HS Taq and 1.5\u0026micro;l 0.05-0.5ng.\u0026micro;l\u003csup\u003e\u0026minus;\u0026thinsp;1\u003c/sup\u003e DNA template, with the following amplification parameters: initial denature of 98℃ for 30 sec, followed by 40 cycles of 98℃ for 10 sec, 52.5℃ for 10 sec and 72℃ for 15 sec and a final extension of 72℃ for 5 min. Sponge mitochondrial cytochrome oxidase subunit 1 (COI mtDNA) gene was amplified using the primers dgLCO1490 (5\u0026prime;-GGTCAACAAATCATAAAGAYATYGG-3\u0026rsquo;) and dgHCO2198 (5\u0026prime;- TAAACTTCAGGGTGACCAAARAAYCA-3\u0026prime;) (Meyer et al. \u003cspan citationid=\"CR82\" class=\"CitationRef\"\u003e2005\u003c/span\u003e) in a 25 \u0026micro;l reaction containing: 12.9\u0026micro;l distilled water (sterile milliQ); 5\u0026micro;l 5x Thermo Fisher Phire Green; 1.0\u0026micro;l 25mM MgCl\u003csub\u003e2\u003c/sub\u003e; 1.0\u0026micro;l 10mg.ml\u003csup\u003e\u0026minus;\u0026thinsp;1\u003c/sup\u003e Life BSA; 1.3\u0026micro;l 10 pMol.\u0026micro;l\u003csup\u003e\u0026minus;\u0026thinsp;1\u003c/sup\u003e forward primer; 1.3\u0026micro;l 10 pMol.\u0026micro;l\u003csup\u003e\u0026minus;\u0026thinsp;1\u003c/sup\u003e reverse primer; 0.5\u0026micro;l 2.5mM dNTP; 0.5\u0026micro;l Thermo Fisher Phire II HS Taq and 1.5\u0026micro;l 0.05-0.5ng.\u0026micro;l\u003csup\u003e\u0026minus;\u0026thinsp;1\u003c/sup\u003e DNA template, with the following amplification parameters: initial denature of 98℃ for 30 sec, followed by 40 cycles of 98℃ for 10 sec, 55℃ for 10 sec and 72℃ for 15 sec and a final extension of 72℃ for 5 min. PCR product yield was checked using Invitrogen E-gels 96 2% agarose.\u003c/p\u003e \u003cp\u003ePCR products were sent to Baseclear B.V. in the Netherlands for forward and reverse sanger sequencing using the same primers than those of the PCR amplification. The resulting sequences were analyzed and edited in Geneious Prime 19.2.3. Reads were assembled using a De Novo Assembly. Contigs were stripped of primer sequences and analyzed for reverse complements and quality. Sequences were aligned in Geneious using MAFFT v7.450 (Katoh et al. \u003cspan citationid=\"CR61\" class=\"CitationRef\"\u003e2002\u003c/span\u003e). A phylogenetic tree was constructed using MrBayes 3.2.6 (Huelsenbeck et al. \u003cspan citationid=\"CR55\" class=\"CitationRef\"\u003e2001\u003c/span\u003e). A GTR substitution model with \u0026lsquo;invgamma\u0026rsquo; rate variation was preferred, setting the chain length to 1100000 with a subsampling frequency of 200 and a burn-in length of 100000.\u003c/p\u003e \u003cp\u003eA first Mantel test was conducted to test the correlation between the phylogenetic distance matrix (based on the 28S rRNA gene) and the geographical distances matrix. The latter was obtained using the GPS coordinates of each sampling site and the \u003cem\u003eearth.dist()\u003c/em\u003e function from the \u0026ldquo;fossil\u0026rdquo; R package (Vavrek \u003cspan citationid=\"CR126\" class=\"CitationRef\"\u003e2011\u003c/span\u003e).\u003c/p\u003e \u003cp\u003e \u003cb\u003eDNA extractions, library preparation and high throughput sequencing of 16S rRNA gene amplicons.\u003c/b\u003e \u003c/p\u003e \u003cp\u003eFor the microbiome analysis, DNA was extracted using the FastDNA\u0026trade; SPIN Kit for Soil (MP Biomedicals, Inc.) following the manufacturer instructions. Sponge tissues were cut into small pieces of approximately 3*1*0.5mm using sterilized tweezers and scalpel blades. An extraction blank, in which no tissue was added to the tube, was also included.\u003c/p\u003e \u003cp\u003eThe library preparation was conducted through a two-step PCR protocol for all samples, together with the extraction blank and two negative controls (mQ water instead of template DNA). For the first PCR, the V4-V5 regions of the 16S rRNA gene were targeted and amplified (Parada et al. \u003cspan citationid=\"CR95\" class=\"CitationRef\"\u003e2016\u003c/span\u003e). We used 515F-Y/926R primers and the KAPA HiFi HotStart Ready Mix PCR Kit (Roche Molecular Systems, Inc.). Reactions were performed in a T100 Thermal Cycler (Bio-Rad, Hercules, CA, United States). The following thermal cycling scheme was set up: initial denaturation at 95\u0026deg;C for 3 min, 30 cycles of denaturation at 98\u0026deg;C for 20 s, annealing at 50\u0026deg;C for 30 s, followed by extension at 72\u0026deg;C for 30 s. The final extension was carried out at 72\u0026deg;C for 5 min.\u003c/p\u003e \u003cp\u003ePCR products from the samples were checked using an E-Gel\u0026trade; (agarose gels at 2%) and the absence of amplification was validated for the negative controls and the extraction blank. PCR products were then cleaned with NucleoMag NGS-Beads (bead volume at 0.9 times the total volume of the sample, Macherey Nagel, D\u0026uuml;ren, Germany) using the VP 407AM-N 96 Pin Magnetic Bead Extractor stamp (V\u0026amp;P Scientific, San Diego, CA, United States). Through a second PCR, 3 \u0026micro;L of the cleaned PCR products were then amplified and labeled using the MiSeq Nextera XTDNA library preparation kit (Illumina, San Diego, CA, United States), with the same thermal cycling scheme limited to 8 cycles. PCR products were then analyzed with the Fragment Analyzer Agilent 5300 using the DNF-910-33 dsDNA Reagent Kit (35\u0026ndash;1,500 bp) protocol (Agilent Technologies, Santa Clara, CA, United States) to confirm successful labeling of the DNA fragments. Negative controls and extraction blanks remained negatives.\u003c/p\u003e \u003cp\u003eThe pooling at the equimolar concentration was performed with QIAgility (Qiagen, Hilden, Germany). The final pool was then cleaned with NucleoMag NGSBeads, eluted in Milli-Q and the DNA concentration was verified using Tapestation 4150 (Kit HSD 5000, Agilent Technologies, Santa Clara, CA, United States). The sequencing was performed on an Illumina MiSeq V3 PE300 platform at BaseClear B.V. (Leiden, the Netherlands).\u003c/p\u003e \u003c/div\u003e \u003cdiv id=\"Sec7\" class=\"Section2\"\u003e \u003ch2\u003e16S rRNA gene metabarcoding data processing\u003c/h2\u003e \u003cp\u003eThe raw reads were first treated by BaseClear B.V. for demultiplexing (using bcl2fastq version 2.20, Illumina) and filtering based on two quality controls (using Illumina Chastity filtering, and a PhiX control signal filtering). The following reads were then processed with the DADA2 workflow, allowing an inference to Amplicon Sequence Variant (ASV) (Callahan et al. \u003cspan citationid=\"CR16\" class=\"CitationRef\"\u003e2016a\u003c/span\u003e, \u003cspan citationid=\"CR15\" class=\"CitationRef\"\u003e2017\u003c/span\u003e). The \u0026ldquo;dada2\u0026rdquo; R package was used following the workflow described in Callahan et al. (\u003cspan citationid=\"CR17\" class=\"CitationRef\"\u003e2016b\u003c/span\u003e) and guidelines described in the online tutorial (\u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://benjjneb.github.io/dada2/tutorial.html\u003c/span\u003e\u003cspan address=\"https://benjjneb.github.io/dada2/tutorial.html\" targettype=\"URL\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e). The parameters used for filtering and trimming reads were the following ones: truncation length of 290 and 240 base pairs for forward and reverse reads, respectively, maxN\u0026thinsp;=\u0026thinsp;0, maxEE\u0026thinsp;=\u0026thinsp;2 and truncQ\u0026thinsp;=\u0026thinsp;2. The error learning step was based on 153 188 440 and 126 776 640 bases for forward and reverse reads, respectively. After the construction of the ASV table, chimeric sequences were filtered and taxonomic assignment was performed using the Silva v138 reference database (Quast et al. \u003cspan citationid=\"CR102\" class=\"CitationRef\"\u003e2013\u003c/span\u003e).\u003c/p\u003e \u003cp\u003eThe ASV and taxonomy tables produced by the pipeline were then combined into a phyloseq object, together with the sample metadata table using the \u0026ldquo;phyloseq\u0026rdquo; R package (McMurdie and Holmes \u003cspan citationid=\"CR78\" class=\"CitationRef\"\u003e2013\u003c/span\u003e). The dataset was then filtered by removing all sequences affiliated to Eukaryota, chloroplast and mitochondria. Data was then decontaminated with two negative controls and the extraction blank used as control samples, using the \u0026ldquo;decontam\u0026rdquo; R package (Davis et al. \u003cspan citationid=\"CR29\" class=\"CitationRef\"\u003e2018\u003c/span\u003e). The dataset was then divided into two subsets, one with the samples for the farming experiment and the second one with those for the biogeographical study. The sampling time (Donor, 9 months, 12 months, 15 months, and 18 months) and sampling sites were used as factors for the statistical analyses (\u003cb\u003eTable \u003cspan refid=\"MOESM1\" class=\"InternalRef\"\u003eS1\u003c/span\u003e and S2\u003c/b\u003e).\u003c/p\u003e \u003c/div\u003e \u003cdiv id=\"Sec8\" class=\"Section2\"\u003e \u003ch2\u003eMetabarcoding and statistical analyses\u003c/h2\u003e \u003cp\u003eAfter the filtering of reads affiliated to Eukaryota, chloroplast and mitochondria (representing together an average of 0.27% of all reads per sample, SD\u0026thinsp;=\u0026thinsp;\u0026plusmn;\u0026thinsp;0.34%), rarefaction curves performed with the \u0026ldquo;vegan\u0026rdquo; R package (Oksanen et al. \u003cspan citationid=\"CR92\" class=\"CitationRef\"\u003e2019\u003c/span\u003e), reached a plateau for all samples (\u003cb\u003eFigure S3\u003c/b\u003e). Since the default DADA2 pipeline is performing a removal of singletons and thus considering the potential loss of very rare taxa, these rarefactions curves and the corresponding \u003cem\u003eα\u003c/em\u003e-diversity metrics should be interpreted with caution. Despite the removal of these singletons, we remain confident about the relevant coverage of the richness displayed in the whole study.\u003c/p\u003e \u003cp\u003eThe \u003cem\u003eα\u003c/em\u003e-diversity measures were estimated with Chao1 (estimated richness), Pielou (evenness) and Shannon (both richness and evenness) indices using the \u0026ldquo;ade4\u0026rdquo; and \u0026ldquo;phyloseq\u0026rdquo; R packages (Dray and Dufour \u003cspan citationid=\"CR33\" class=\"CitationRef\"\u003e2007\u003c/span\u003e; McMurdie and Holmes \u003cspan citationid=\"CR78\" class=\"CitationRef\"\u003e2013\u003c/span\u003e) and the rarefied datasets (rarefaction performed to the minimum library size, i.e. 31559 and 34103 reads, for the farming experiment and the biogeographical datasets respectively). Depending on the normal distribution assessed using a Shapiro test, the significant effect of the sampling times and sites on the diversity metrics was tested using ANOVA followed by HSD Tukey\u0026rsquo;s tests as parametric tests, as well as Kruskal-Wallis followed by pairwise Wilcoxon tests as non-parametric ones. Following recommendations for compositional approaches from McMurdie and Holmes (\u003cspan citationid=\"CR79\" class=\"CitationRef\"\u003e2014\u003c/span\u003e) and Gloor et al. (\u003cspan citationid=\"CR45\" class=\"CitationRef\"\u003e2017\u003c/span\u003e), all other analyses were conducted without rarefaction, using the phyloseq R package and the datasets normalized to the total number of sequences per sample. For both datasets, the \u003cem\u003eβ\u003c/em\u003e-diversity of all samples was analyzed with a non-metric multidimensional scaling (NMDS) using the Bray-Curtis dissimilarity. Differences of \u003cem\u003eβ\u003c/em\u003e-diversity between groups were statistically checked with one-way PERMANOVA tests followed by pairwise Adonis tests.\u003c/p\u003e \u003cp\u003eDistance-based Redundancy Analysis (db-RDA) using the Bray-Curtis dissimilarity was conducted to investigate the effect of environmental parameters on the \u003cem\u003eβ\u003c/em\u003e-diversity of the prokaryotic community of \u003cem\u003eD. metachromia\u003c/em\u003e, using the \u0026ldquo;vegan\u0026rdquo; R package. A first analysis was performed with all samples from the farming dataset and tested with all available environmental parameters (temperature, salinity, PO\u003csub\u003e4\u003c/sub\u003e\u003csup\u003e3\u0026minus;\u003c/sup\u003e, NO\u003csub\u003ex\u003c/sub\u003e, Si[OH]\u003csub\u003e4\u003c/sub\u003e, and chlorophyll \u003cem\u003ea\u003c/em\u003e). Since some parameters were not described for all sampling sites from the biogeographical study, the analysis was conducted on a restricted dataset with a selection of samples for which the temperature, salinity, PO\u003csub\u003e4\u003c/sub\u003e\u003csup\u003e3\u0026minus;\u003c/sup\u003e, NO\u003csub\u003ex\u003c/sub\u003e and Si(OH)\u003csub\u003e4\u003c/sub\u003e were measured (\u003cb\u003eTable S2\u003c/b\u003e and \u003cb\u003eS3\u003c/b\u003e). This restricted biogeographical dataset includes only samples from French Polynesian sites (except Makemo). The \u003cem\u003eordiR2step\u003c/em\u003e() function was used to verify that all these environmental factors can be selected as significant constraints to be used as explanatory variables. The overall significance of the two models was tested using the \u003cem\u003eanova.cca\u003c/em\u003e() function with 999 permutations.\u003c/p\u003e \u003cp\u003eFollowing db-RDA, variance partitioning analyses (Borcard et al. \u003cspan citationid=\"CR12\" class=\"CitationRef\"\u003e1992\u003c/span\u003e; Cleary et al. \u003cspan citationid=\"CR27\" class=\"CitationRef\"\u003e2022\u003c/span\u003e) were conducted using the \u003cem\u003evarpart\u003c/em\u003e() function from the \u0026ldquo;vegan\u0026rdquo; package for both the farming and the restricted biogeographical datasets. These analyses were performed to assess the relative percentage of the variance of the \u003cem\u003eβ\u003c/em\u003e-diversity explained by the selected environmental parameters, using the formula (~\u0026thinsp;temperature, ~ salinity, ~ PO\u003csub\u003e4\u003c/sub\u003e \u003csup\u003e3\u0026minus;\u003c/sup\u003e + NO\u003csub\u003ex\u003c/sub\u003e + Si[OH]\u003csub\u003e4\u003c/sub\u003e) for both datasets. Additionally, the variance of the \u003cem\u003eβ\u003c/em\u003e-diversity explained by environmental parameters was also assessed through a variance partitioning analysis combining the geographical distances. The geographical distances matrix was transformed using the \u003cem\u003epcnm\u003c/em\u003e function and the analysis was conducted using the following formula: (~\u0026thinsp;temperature\u0026thinsp;+\u0026thinsp;salinity\u0026thinsp;+\u0026thinsp;PO\u003csub\u003e4\u003c/sub\u003e \u003csup\u003e3\u0026minus;\u003c/sup\u003e + NO\u003csub\u003ex\u003c/sub\u003e + Si[OH]\u003csub\u003e4\u003c/sub\u003e, ~ \u003cem\u003epcnm\u003c/em\u003e[geographical.dist]\u003cspan\u003e$\u003c/span\u003evectors).\u003c/p\u003e \u003cp\u003eFor the biogeographical study, Mantel tests (Spearman correlation, 9999 permutations) were conducted within all samples to evaluate the correlation between the \u003cem\u003eβ\u003c/em\u003e-diversity distance matrix and the geographical distance matrix. Excepted for samples not barcoded (Lagon 1 and Lagon 2, Fig.\u0026nbsp;\u003cspan refid=\"Fig1\" class=\"InternalRef\"\u003e1\u003c/span\u003e, \u003cb\u003eTable S2\u003c/b\u003e), an additional Mantel test was also conducted with the biogeographical dataset to assess the correlation between the \u003cem\u003eβ\u003c/em\u003e-diversity distance matrix and the phylogenetic distance matrix.\u003c/p\u003e \u003cp\u003eThe \"core community\" was obtained using the \u0026ldquo;microbiome\u0026rdquo; R package (Lahti et al. \u003cspan citationid=\"CR63\" class=\"CitationRef\"\u003e2017\u003c/span\u003e) with the \u003cem\u003ecore\u003c/em\u003e() function: a prevalence threshold of 90% in all samples (including both geographical and farming datasets) and an abundance threshold of 0.01% were kept. Percentages of core ASVs numbers and sequences were estimated by comparison with the total community and statistically tested using Shapiro, ANOVA and HSD Tukey\u0026rsquo;s tests.\u003c/p\u003e \u003c/div\u003e"},{"header":"Results","content":"\u003cdiv id=\"Sec10\" class=\"Section2\"\u003e \u003ch2\u003eEnvironmental conditions, and sponge monitoring during farming experiment\u003c/h2\u003e \u003cp\u003eThe water temperature increased from 27.3 to 28.3\u0026deg;C for farmed explants of 9 and 12 months respectively, corresponding to the beginning of the wet and warm tropical season (i.e. September to December) (\u003cb\u003eFigure S4A\u003c/b\u003e). The peak of 29.3\u0026deg;C was reached for \u0026ldquo;15 months\u0026rdquo; samples, followed three months later by a slight decrease of -0.5\u0026deg;C. Indeed, the \u0026ldquo;18 months\u0026rdquo; sampling corresponds to the start of the dry and cold season in the tropics (i.e. March-April). The salinity varied between 35.04 and 32 PSU, between December 2019 and March 2021 (i.e. \u0026ldquo;Donor\u0026rdquo; and \u0026ldquo;15 months\u0026rdquo; respectively, \u003cb\u003eFigure S4A\u003c/b\u003e). Yet, the average salinity dropped in May 2021 (i.e. \u0026ldquo;18 months\u0026rdquo;) to 29.9 PSU, corresponding to the return of an El Ni\u0026ntilde;o\u0026ndash;Southern Oscillation climate pattern (World Meteorological Association \u003cspan citationid=\"CR131\" class=\"CitationRef\"\u003e2022\u003c/span\u003e).\u003c/p\u003e \u003cp\u003eRegarding the nutrients, no clear trend was observed in the phosphate plot, as the averages remained very similar throughout the farming period. Nitrogen oxides and orthosilicic acid seem to vary in opposite ways (\u003cb\u003eFigure S4A\u003c/b\u003e), (NO\u003csub\u003ex\u003c/sub\u003e) concentrations being the lowest for farmed samples of \u0026ldquo;12 months\u0026rdquo; and \u0026ldquo;15 months\u0026rdquo; whereas peaks in (Si[OH]\u003csub\u003e4\u003c/sub\u003e) averages are noticed for the same months (i.e. December 2020 and March 2021, respectively).\u003c/p\u003e \u003cp\u003eWithout including the Red Sea, the Maldives, Indonesia and the atoll of Makemo (within French Polynesia) since no environmental parameters were measured for these sites (\u003cb\u003eTable S2 and S3\u003c/b\u003e), spatial analyses of all environmental parameters revealed significant differences between French Polynesian islands (Kruskal-Wallis: \u003cem\u003ep\u003c/em\u003e\u0026thinsp;\u0026lt;\u0026thinsp;0.02; \u003cb\u003eTable S4\u003c/b\u003e). Lower sea temperatures for the southern part of the Tuamotu and the Gambier archipelagos (i.e. Tematangi and Mangareva) were observed compared to warmer northern islands featuring on average 4 to 5\u0026deg;C temperatures (\u003cb\u003eFigure S4B\u003c/b\u003e). The atoll of Rangiroa was characterized by significantly lower salinity levels (averages ranging between 31 and 33 PSU depending on the location) in contrast to the ones registered for Raroia and Tematangi, usually above 36 PSU (Pairwise Wilcoxon test, \u003cb\u003eFigure S4B\u003c/b\u003e). Phosphate and silicate concentrations appear to be relatively similar for all sampled sites, except for the atoll of Raroia, characterized by significantly higher concentrations compared to all sites except Mangareva and Tetiaroa (Pairwise Wilcoxon test, \u003cb\u003eFigure S4B\u003c/b\u003e). Finally, nitrogen oxide values fall along a north-west/south-east gradient from 0.4\u0026micro;g.L\u003csup\u003e\u0026minus;\u0026thinsp;1\u003c/sup\u003e in Tetiaroa to 0.25\u0026micro;g.L\u003csup\u003e\u0026minus;\u0026thinsp;1\u003c/sup\u003e in Raroia, 0.1\u0026micro;g.L\u003csup\u003e\u0026minus;\u0026thinsp;1\u003c/sup\u003e in Tematangi and barely 0.05\u0026micro;g.L\u003csup\u003e\u0026minus;\u0026thinsp;1\u003c/sup\u003e in Mangareva.\u003c/p\u003e \u003cp\u003eHigh and stable survival rates were observed in Avatoru throughout the farming experiment hereby addressed, ranging between 69\u0026ndash;71% (\u003cb\u003eFigure S5A\u003c/b\u003e). For volumes, a higher peak average reaching 99\u0026thinsp;\u0026plusmn;\u0026thinsp;59 cm\u003csup\u003e3\u003c/sup\u003e after 12 months of farming can be observed, but no significant changes were determined (\u003cb\u003eFigure S5B\u003c/b\u003e, Kruskal-Wallis, \u003cem\u003ep\u003c/em\u003e\u0026thinsp;\u0026gt;\u0026thinsp;0.05). Our results display an annual growth rate of 100%, as the standardized initial size of each explant was set close to 50 cm\u003csup\u003e3\u003c/sup\u003e.\u003c/p\u003e \u003c/div\u003e \u003cdiv id=\"Sec11\" class=\"Section2\"\u003e \u003ch2\u003eSponge phylogeny\u003c/h2\u003e \u003cp\u003ePhylogenetic trees based on 28S rRNA gene and COI mtDNA sequences have been constructed using \u003cem\u003eDactylospongia metachromia\u003c/em\u003e samples collected from the farming trials and additional material from the sponge collections of Naturalis Biodiversity Center. The phylogeny of 28S rRNA genes shows some intraspecific variation with several clades forming within \u003cem\u003eD. metachromia\u003c/em\u003e. Despite a significant correlation between the distance matrix based on 28S rRNA gene sequences and the geographical distances, for the whole dataset (Mantel test: r\u0026thinsp;=\u0026thinsp;0.1785; \u003cem\u003ep\u003c/em\u003e\u0026thinsp;=\u0026thinsp;0.0241), no clear geographical pattern was noticed within the French Polynesian sites (\u003cb\u003eFigure S6\u003c/b\u003e). More precisely, these samples formed two different clades, with (i) Makemo, Mangareva, and Raroia samples more related to Indonesia and Red Sea samples, and (ii) samples from Rangiroa, Tematangi, and Tetiaroa. Moreover, COI mtDNA shows no intraspecific variation between all samples, with a tree forming a single clade.\u003c/p\u003e \u003c/div\u003e \u003cdiv id=\"Sec12\" class=\"Section2\"\u003e \u003ch2\u003eProkaryotic diversity of farmed samples\u003c/h2\u003e \u003cp\u003eFor the farming experiment, no significant differences were observed for Shannon, Chao1, and Pielou indices (Fig.\u0026nbsp;\u003cspan refid=\"Fig2\" class=\"InternalRef\"\u003e2\u003c/span\u003eA, \u003cspan refid=\"Fig2\" class=\"InternalRef\"\u003e2\u003c/span\u003eB and \u003cspan refid=\"Fig2\" class=\"InternalRef\"\u003e2\u003c/span\u003eC, respectively) between samples collected at different time points (\u0026ldquo;donor samples\u0026rdquo;, and samples farmed during \u0026ldquo;9-\u0026rdquo;, \u0026ldquo;12-\u0026rdquo;, \u0026ldquo;15-\u0026rdquo; and \u0026ldquo;18-months\u0026rdquo;, ANOVA: \u003cem\u003ep\u003c/em\u003e\u0026thinsp;\u0026gt;\u0026thinsp;0.05 for the three indices, \u003cb\u003eTable S5\u003c/b\u003e).\u003c/p\u003e \u003cp\u003eAs for the NMDS analysis of \u003cem\u003eβ\u003c/em\u003e-diversity, distinct clusters of samples appear only for \u0026ldquo;15 months'' and \u0026ldquo;18 months\u0026rdquo; samples, while no clear patterns were observed for \u0026ldquo;Donor\u0026rdquo;, \u0026ldquo;9 months\u0026rdquo; and \u0026ldquo;12 months\u0026rdquo; samples (Fig.\u0026nbsp;\u003cspan refid=\"Fig2\" class=\"InternalRef\"\u003e2\u003c/span\u003eD). Despite a significant difference observed with the PERMANOVA test (\u003cem\u003ep\u003c/em\u003e\u0026thinsp;=\u0026thinsp;0.003), no significant pairwise differences between the sample times were noticed (pairwise Adonis: \u003cem\u003ep\u003c/em\u003e\u0026thinsp;\u0026gt;\u0026thinsp;0.05 for all comparisons). The ASVs scores of the NMDS analysis revealed major ASVs (with relative abundances\u0026thinsp;\u0026gt;\u0026thinsp;2% on average) affiliated with several phyla including Crenarchaeota (e.g. ASV1), Chloroflexi (e.g. ASV5 and AVV10), the PAUC34f clade (e.g. ASV9), Entheonellaeota (e.g. ASV4) and Acidobacteriota (e.g. ASV7) (Fig.\u0026nbsp;\u003cspan refid=\"Fig2\" class=\"InternalRef\"\u003e2\u003c/span\u003eE). No clear clusters of ASVs associated with a specific phylum were observed. However, one minor ASV located on the left side of the NMDS and affiliated to the Acidobacteriota appears to be specifically associated with the cluster of \u0026ldquo;18 months\u0026rdquo; samples.\u003c/p\u003e \u003cp\u003e \u003c/p\u003e \u003cp\u003eThe db-RDA analysis shows three distinct clusters (Fig.\u0026nbsp;\u003cspan refid=\"Fig2\" class=\"InternalRef\"\u003e2\u003c/span\u003eF): (i) one gathering \u0026ldquo;Donor\u0026rdquo; and \u0026ldquo;15 months\u0026rdquo; samples, positively associated to the temperature and to a lesser extent chlorophyll \u003cem\u003ea\u003c/em\u003e, (ii) one gathering \u0026ldquo;9 months\u0026rdquo; and \u0026ldquo;12 months\u0026rdquo;, positively and negatively associated to the salinity and temperature, respectively, and (iii) a last one gathering \u0026ldquo;18 months\u0026rdquo; samples, positively associated to PO\u003csub\u003e4\u003c/sub\u003e\u003csup\u003e3\u0026minus;\u003c/sup\u003e and NO\u003csub\u003ex\u003c/sub\u003e concentrations and negatively to the salinity (ANOVA.CCA: \u003cem\u003ep\u003c/em\u003e\u0026thinsp;=\u0026thinsp;0.002, \u003cb\u003eTable S6B\u003c/b\u003e). The ASVs scores of the db-RDA analysis did not reveal ASV clusters specifically associated with environmental parameters (Fig.\u0026nbsp;\u003cspan refid=\"Fig2\" class=\"InternalRef\"\u003e2\u003c/span\u003eG). However, ASV5 affiliated to the Chloroflexi which happens to be amongst the dominant ones, was found positively correlated to the temperature.\u003c/p\u003e \u003cp\u003eThe variance partitioning analysis with environmental parameters revealed that 3% of the \u003cem\u003eβ\u003c/em\u003e-diversity variance is explained by temperature, 10% by salinity, and 6% by nutrients (PO\u003csub\u003e4\u003c/sub\u003e\u003csup\u003e3\u0026minus;\u003c/sup\u003e, NO\u003csub\u003ex\u003c/sub\u003e and Si[OH]\u003csub\u003e4\u003c/sub\u003e together); while the percentage explained by the combination of salinity and nutrients reach 6%. Temperature, salinity and nutrients combined explain 5% of the variance (\u003cb\u003eFigure S7A\u003c/b\u003e).\u003c/p\u003e \u003cp\u003eAt the family level, the prokaryotic community composition did not reveal any major changes between sample times (for the farming experiment, \u003cb\u003eFigure S7B\u003c/b\u003e). Most of the families observed were affiliated with three phyla: Acidobacteriota, Chloroflexi, and Proteobacteria. The main families observed in the whole dataset belong to an unidentified family from the Dehalococcoidia class (Chloroflexi phylum), the Nitrosopumilaceae family (Crenarchaeota phylum), the Caldilineaceae family (Chloroflexi phylum), an unaffiliated family from the candidate phylum PAUC34f and two unaffiliated families from the Acidobacteriae and Vicinamibacteria classes (Acidobacteria phylum). The above are respectively representing 11.5%, 6.7%, 6.0%, 6.0%, 5.5% and 4.9% of the community (SD: \u0026plusmn;4.5, \u0026plusmn;\u0026thinsp;5.6, \u0026plusmn;\u0026thinsp;2.8, \u0026plusmn;\u0026thinsp;2.3, \u0026plusmn;\u0026thinsp;2.3 and \u0026plusmn;\u0026thinsp;1.8, respectively).\u003c/p\u003e \u003c/div\u003e \u003cdiv id=\"Sec13\" class=\"Section2\"\u003e \u003ch2\u003eProkaryotic diversity of the biogeographical samples\u003c/h2\u003e \u003cp\u003eWhen focusing on the biogeographical study (\u003cb\u003eFigure S8A\u003c/b\u003e), significant differences between sampling sites were observed for the three \u003cem\u003eα\u003c/em\u003e-diversity indices (ANOVA: \u003cem\u003ep\u003c/em\u003e\u0026thinsp;\u0026lt;\u0026thinsp;0.001, \u003cem\u003ep\u003c/em\u003e\u0026thinsp;=\u0026thinsp;0.001 for Shannon and Pielou, respectively; Kruskal-Wallis: \u003cem\u003ep\u003c/em\u003e\u0026thinsp;=\u0026thinsp;0.007 for Chao1; \u003cb\u003eTable S5\u003c/b\u003e). For the Shannon and Chao1 indices, the highest values were observed for Tiputa samples (\u003cb\u003eFigure S8A\u003c/b\u003e). For the Pielou index, a more stable tendency was discerned, despite lower values observed for Raroia compared to Rangiroa (Avatoru) and Temantagi (HSD Tukey\u0026rsquo;s test, \u003cb\u003eFigure S8A\u003c/b\u003e).\u003c/p\u003e \u003cp\u003eThe NMDS plot of the samples from the biogeographical study (Fig.\u0026nbsp;\u003cspan refid=\"Fig3\" class=\"InternalRef\"\u003e3\u003c/span\u003eA) showed a distinct clustering according to sampling sites (PERMANOVA: \u003cem\u003ep\u003c/em\u003e\u0026thinsp;\u0026lt;\u0026thinsp;0.001, \u003cb\u003eTable S6A\u003c/b\u003e), with main differences observed in the first axis (NMDS1) between French Polynesian and other sites (Red Sea, Indonesia and Maldives). Moreover, differences were observed within the French Polynesian sites on the second axis (NMDS2), with a clear clustering observed between each island (Raroia, Temantangi, Mangareva, Tetiaroa and Rangiroa). When focusing on a smaller geographical scale, with sampling sites from the Rangiroa atoll (Tiputa, Lagon1, Lagon2 and Avatoru), no significant differences of \u003cem\u003eβ\u003c/em\u003e-diversity were observed, except between Lagon1 and Tiputa (pairwise Adonis test). When focusing on the NMDS plot displaying the ASV scores of the biogeographical study (Fig.\u0026nbsp;\u003cspan refid=\"Fig3\" class=\"InternalRef\"\u003e3\u003c/span\u003eB), major ASVs were found to be notably affiliated to the Crenarchaeota (e.g. ASV1), Acidobacteriota, Chloroflexi and PAUC34f clades. While no specific patterns of ASV clusters could be associated with a specific group of samples, ASV41 affiliated with the Crenarchaeota was found to be closely related to the Red Sea samples.\u003c/p\u003e \u003cp\u003e \u003c/p\u003e \u003cp\u003eConsidering environmental parameters as explanatory variables (available only for French Polynesian samples, Fig.\u0026nbsp;\u003cspan refid=\"Fig3\" class=\"InternalRef\"\u003e3\u003c/span\u003eC and \u003cspan refid=\"Fig3\" class=\"InternalRef\"\u003e3\u003c/span\u003eD), the db-RDA was conducted using the Bray-Curtis dissimilarity. The analysis suggested that temperature and salinity were the main parameters explaining \u003cem\u003eβ\u003c/em\u003e-diversity differences, together with PO\u003csub\u003e4\u003c/sub\u003e\u003csup\u003e3\u0026minus;\u003c/sup\u003e concentration to a lesser extent (ANOVA.CCA: \u003cem\u003ep\u003c/em\u003e\u0026thinsp;\u0026lt;\u0026thinsp;0.001, \u003cb\u003eTable S4B\u003c/b\u003e). Temperature and salinity appears positively and negatively associated with Rangiroa samples, respectively, while the opposite situation appears for Tetiaroa, Tematangi and Mangareva samples (Fig.\u0026nbsp;\u003cspan refid=\"Fig3\" class=\"InternalRef\"\u003e3\u003c/span\u003eC). Additionally, PO\u003csub\u003e4\u003c/sub\u003e\u003csup\u003e3\u0026minus;\u003c/sup\u003e concentration was found positively associated with the \u003cem\u003eβ\u003c/em\u003e-diversity of Raroia samples. No specific patterns were observed within the ASV distribution on the db-RDA plot. However, two major ASVs can be noticed: (i) ASV1 affiliated to the Crenarchaeota and positively associated with the salinity, (ii) ASV5 affiliated to the Chloroflexi, mainly occurring within the Rangiroa samples and positively associated with the temperature.\u003c/p\u003e \u003cp\u003eA first variance partitioning analysis (\u003cb\u003eFigure S8B\u003c/b\u003e) revealed that temperature, salinity, and nutrients [combining PO\u003csub\u003e4\u003c/sub\u003e\u003csup\u003e3\u0026minus;\u003c/sup\u003e, NO\u003csub\u003ex\u003c/sub\u003e and Si(OH)\u003csub\u003e4\u003c/sub\u003e] explain 2%, 6% and 5% of the \u003cem\u003eβ\u003c/em\u003e-diversity variance, respectively. More importantly, temperature and salinity combined explain 18% of the variance, while these percentages reach 3% and 8% when nutrients are combined either with the temperature or the salinity, respectively.\u003c/p\u003e \u003cp\u003eFor the biogeographical study, a significant correlation was observed between the biogeographical and \u003cem\u003eβ\u003c/em\u003e-diversity (Bray-Curtis) dissimilarities (Mantel test: r\u0026thinsp;=\u0026thinsp;0.8722; \u003cem\u003ep\u003c/em\u003e\u0026thinsp;\u0026lt;\u0026thinsp;0.001). Furthermore, a variance partition analysis was conducted with environmental parameters (temperature, salinity and nutrients combined) and the geographical distances as explanatory factors of the \u003cem\u003eβ\u003c/em\u003e-diversity variance (\u003cb\u003eFigure S8C\u003c/b\u003e). Results revealed that 9% of the variance is explained by the environmental parameters, 10% by the geographical distances, and 35% is explained by both datasets together. The residual variance (variance left unexplained) reached 46% for this analysis. Additionally, when comparing distance matrices from the sponge phylogeny and \u003cem\u003eβ\u003c/em\u003e-diversity analyses, no significant correlations were observed between the two datasets (Mantel test: r\u0026thinsp;=\u0026thinsp;0.1213; \u003cem\u003ep\u003c/em\u003e\u0026thinsp;=\u0026thinsp;0.0693).\u003c/p\u003e \u003cp\u003eFinally, the prokaryotic community composition of the biogeographical samples at the family level did not reveal any major changes across space (\u003cb\u003eFigure S8D\u003c/b\u003e), with similar percentages of relative abundances observed for the temporal study.\u003c/p\u003e \u003cp\u003e \u003cb\u003eCore community of all\u003c/b\u003e \u003cb\u003eD. metachromia\u003c/b\u003e \u003cb\u003esamples\u003c/b\u003e\u003c/p\u003e \u003cp\u003eThe core community of \u003cem\u003eD. metachromia\u003c/em\u003e was defined by keeping only ASVs occurring at least in 90% of all samples (samples from both biogeographical and farming studies together).\u003c/p\u003e \u003cp\u003e \u003c/p\u003e \u003cp\u003eA total of 18 ASVs were determined as part of the core community (\u003cb\u003eTable S7\u003c/b\u003e), representing 0.3% of the total number of all ASVs in both datasets combined. The relative percentage of numbers of core ASVs per sample (\u003cb\u003eFigure S9\u003c/b\u003e), ranging between 4.1% and 8.8% (average 6.6%, SD\u0026thinsp;\u0026plusmn;\u0026thinsp;1.1), was not significantly different among sampling times for the farming experiment (ANOVA: \u003cem\u003ep\u003c/em\u003e\u0026thinsp;=\u0026thinsp;0.03, but Tukey\u0026rsquo;s test: \u003cem\u003ep\u003c/em\u003e\u0026thinsp;\u0026gt;\u0026thinsp;0.05) or among sampling sites for the biogeographical approach (ANOVA: \u003cem\u003ep\u003c/em\u003e\u0026thinsp;\u0026gt;\u0026thinsp;0.05). The relative percentage of sequences of core ASVs (Fig.\u0026nbsp;\u003cspan refid=\"Fig4\" class=\"InternalRef\"\u003e4\u003c/span\u003e), ranging between 5.9% and 14.1% (average 10.6%, SD: \u0026plusmn;2.1), did not show significant differences between sampling times (for the farming experiment, ANOVA: \u003cem\u003ep\u003c/em\u003e\u0026thinsp;\u0026gt;\u0026thinsp;0.05) or sampling sites (for the biogeographical study, ANOVA: \u003cem\u003ep\u003c/em\u003e\u0026thinsp;\u0026gt;\u0026thinsp;0.05). These sequences of the 18 core ASVs were affiliated to 13 genera (belonging mostly to unidentified families). The most abundant core ASVs were affiliated to the phyla Acidobacteriota, Chloroflexi, Dadabacteria, Myxococcota, Proteobacteria and the candidate phylum PAUC34f, representing in average 2.3%, 1.7%, 1.8%, 1.4%, 1.4%, 1.5%, of the total community, respectively (SD: \u0026plusmn;0.7, \u0026plusmn;\u0026thinsp;0.8, \u0026plusmn;\u0026thinsp;0.7, \u0026plusmn;\u0026thinsp;0.7, \u0026plusmn;\u0026thinsp;0.6, \u0026plusmn;\u0026thinsp;1.0, respectively). Some of these phyla were represented by different classes such as Acidobacteriae, Thermoanaerobaculia, and Vicinamibacteria (Acidobacteria), or Alpha- and Gamma-proteobacteria (Proteobacteria).\u003c/p\u003e \u003c/div\u003e"},{"header":"Discussion","content":"\u003cp\u003eAs a model of interest for the production of pharmacological promising natural products (such as ilimaquinone and 5-epi-ilimaquinone), \u003cem\u003eDactylospongia metachromia\u003c/em\u003e has been the subject of an aquaculture program in French Polynesia (Maslin \u003cspan citationid=\"CR76\" class=\"CitationRef\"\u003e2022\u003c/span\u003e). The production of sponge natural products and, more generally, the whole sponge physiology, is known to be intimately linked to its microbiome (Indraningrat et al. \u003cspan citationid=\"CR59\" class=\"CitationRef\"\u003e2016\u003c/span\u003e; Bibi et al. \u003cspan citationid=\"CR10\" class=\"CitationRef\"\u003e2020\u003c/span\u003e). When considering the launching of a sponge aquaculture industry, it thus appears essential to study the effects of transplantation and environmental changes on the whole holobiont of the targeted species (Thomas et al. \u003cspan citationid=\"CR121\" class=\"CitationRef\"\u003e2016\u003c/span\u003e; Pita et al. \u003cspan citationid=\"CR100\" class=\"CitationRef\"\u003e2018\u003c/span\u003e).\u003c/p\u003e \u003cp\u003e \u003cb\u003eA temporal study revealing an overall stability of the holobiont during farming experiment, despite a moderate effect of environmental changes\u003c/b\u003e \u003c/p\u003e \u003cp\u003eThe analysis of the prokaryotic diversity of \u003cem\u003eD. metachromia\u003c/em\u003e was conducted on a temporal dimension relative to the duration of the farming experiment targeting sponge explants on a barrier reef area in Rangiroa, named Avatoru. No significant differences of \u003cem\u003eα\u003c/em\u003e- diversity between the different months of farming were obtained, with a higher intra-individual variability suggesting an overall stability of the prokaryotic community. Similarly, no clear differences of \u003cem\u003eβ\u003c/em\u003e-diversity within the overall community were provided by the NMDS analysis. However, several environmental parameters including temperature, salinity and nutrients appear to have a slight impact within the community, notably impacting a small fraction of it [i.e. two ASVs affiliated to Crenaracheaota (ASV1) and Chloroflexi (ASV5)].\u003c/p\u003e \u003cp\u003eThe prokaryotic community associated with marine sponges is commonly described as geographically and temporally stable. This is explained by the importance of symbiotic taxa displaying specific relationships with their sponge hosts compared to other environmental samples such as seawater, sediments or other holobionts (Hentschel et al. \u003cspan citationid=\"CR51\" class=\"CitationRef\"\u003e2002\u003c/span\u003e; Taylor et al. \u003cspan citationid=\"CR119\" class=\"CitationRef\"\u003e2004\u003c/span\u003e; Hill et al. \u003cspan citationid=\"CR52\" class=\"CitationRef\"\u003e2006\u003c/span\u003e; Schmitt et al. \u003cspan citationid=\"CR109\" class=\"CitationRef\"\u003e2012b\u003c/span\u003e; Chen et al. \u003cspan citationid=\"CR23\" class=\"CitationRef\"\u003e2020\u003c/span\u003e; Hardoim et al. \u003cspan citationid=\"CR49\" class=\"CitationRef\"\u003e2021\u003c/span\u003e). The HMA-LMA specificity is also linked to distinct prokaryotic composition (Campana et al. \u003cspan citationid=\"CR19\" class=\"CitationRef\"\u003e2021b\u003c/span\u003e; Lesser et al. \u003cspan citationid=\"CR69\" class=\"CitationRef\"\u003e2022\u003c/span\u003e), while no clear evidence of stability changes was observed over seasons and years (Erwin et al. \u003cspan citationid=\"CR41\" class=\"CitationRef\"\u003e2015\u003c/span\u003e).\u003c/p\u003e \u003cp\u003eStudies have shown that notable differences in microbiota composition at multiple taxonomic levels exist between healthy and diseased sponges or under environmental stress (Webster et al. \u003cspan citationid=\"CR128\" class=\"CitationRef\"\u003e2008\u003c/span\u003e). Sponges belonging to the order Dictyoceratida have been reported as particularly sensitive to pathogen epidemics (Smith \u003cspan citationid=\"CR114\" class=\"CitationRef\"\u003e1941\u003c/span\u003e; Wulff \u003cspan citationid=\"CR132\" class=\"CitationRef\"\u003e2001\u003c/span\u003e). Still, identifying the temporally stable microbial community found in the present work, no potentially harmful bacteria taxa were highlighted in the analysis of the community composition. Zaneveld et al. (\u003cspan citationid=\"CR135\" class=\"CitationRef\"\u003e2017\u003c/span\u003e) introduced the Anna Karenina principle for sponge-holobionts, stating that an increased dispersion of the \u003cem\u003eβ\u003c/em\u003e-diversity through time could be linked to multiple stress scenarios led by the emergence of diseases. In this study, the opposite situation seems to occur, with a less dispersed \u003cem\u003eβ\u003c/em\u003e-diversity for the last two sampling times (15 and 18 months). While the transplantation of wild sponge explants farmed \u003cem\u003ein situ\u003c/em\u003e to closed or semi-closed systems has demonstrated important changes in the microbiota (Mohamed et al. \u003cspan citationid=\"CR83\" class=\"CitationRef\"\u003e2008a\u003c/span\u003e, \u003cspan citationid=\"CR84\" class=\"CitationRef\"\u003eb\u003c/span\u003e), our study indicates that wild transfers can only induce moderate effects on the sponge holobiont. Such stability throughout natural contexts was also observed for other sponges, namely \u003cem\u003eTethya bergquistae\u003c/em\u003e (C\u0026aacute;rdenas et al. \u003cspan citationid=\"CR21\" class=\"CitationRef\"\u003e2014\u003c/span\u003e), \u003cem\u003eAplysina cavernicola\u003c/em\u003e (Thoms et al. \u003cspan citationid=\"CR122\" class=\"CitationRef\"\u003e2003\u003c/span\u003e), \u003cem\u003eAplysina archeri\u003c/em\u003e (Hunting et al. \u003cspan citationid=\"CR57\" class=\"CitationRef\"\u003e2013\u003c/span\u003e) and \u003cem\u003eHymeniacidon heliophila\u003c/em\u003e (Weigel and Erwin \u003cspan citationid=\"CR129\" class=\"CitationRef\"\u003e2017\u003c/span\u003e). Additionally, a similar temporal stability was also recently described within \u003cem\u003eHymeniacidon perlevis\u003c/em\u003e and \u003cem\u003eSuberites massa\u003c/em\u003e over seasonal changes (Lamb and Watts \u003cspan citationid=\"CR64\" class=\"CitationRef\"\u003e2023\u003c/span\u003e).\u003c/p\u003e \u003cp\u003eInterestingly enough, despite highly heterogeneous ratios between ilimaquinone and 5-epi-ilimaquinone within the explants, no significant variation was either obtained in their concentrations over time or specimens during farming experiments (Maslin \u003cspan citationid=\"CR76\" class=\"CitationRef\"\u003e2022\u003c/span\u003e). Together with the stable \u003cem\u003eα\u003c/em\u003e- and \u003cem\u003eβ\u003c/em\u003e-diversity observed over time, these results corroborate the absence of a negative influence of the aquaculture system during the farming period. This finding allows us to consider a chemical exploitation in line with the positive enlargement of the sponge explants in order to seek an optimal sampling time. Nevertheless, being of bacterial origin or not, the average concentrations of secondary metabolites found in the extracts of \u003cem\u003eD. metachromia\u003c/em\u003e are still insufficient to match the needs of the pharmaceutical market towards new drug synthesis (Sipkema et al. \u003cspan citationid=\"CR112\" class=\"CitationRef\"\u003e2005b\u003c/span\u003e; Malve \u003cspan citationid=\"CR75\" class=\"CitationRef\"\u003e2016\u003c/span\u003e).\u003c/p\u003e \u003cdiv id=\"Sec15\" class=\"Section2\"\u003e \u003ch2\u003eA complementary biogeographical approach disclosing site-specific communities and confirming the importance of temperature and salinity\u003c/h2\u003e \u003cp\u003eValues of \u003cem\u003eα\u003c/em\u003e-diversity metrics were found to slightly differ between sampling sites, especially for the Red Sea samples characterized by a lower Shannon index compared to Rangiroa and Tematangi atolls. Such differences for this particular region might enlighten host-specific relationships, since for other sponges such as \u003cem\u003eXestospongia\u003c/em\u003e spp. the \u003cem\u003eα\u003c/em\u003e-diversity did not reveal lower values for samples from Saudi Arabia compared to the Australian coast or the Indo-Pacific region (Thomas et al. \u003cspan citationid=\"CR121\" class=\"CitationRef\"\u003e2016\u003c/span\u003e; Swierts et al. \u003cspan citationid=\"CR117\" class=\"CitationRef\"\u003e2018\u003c/span\u003e). However, comparison of prokaryotic richness between different studies should be interpreted with caution since many methodological biases can be involved, from the molecular preparation to the data processing.\u003c/p\u003e \u003cp\u003eThe \u003cem\u003eβ\u003c/em\u003e-diversity of the prokaryotic community associated with \u003cem\u003eD. metachromia\u003c/em\u003e samples was found to be highly site-specific, whether at the global scale (when considering samples from French Polynesia together with Red Sea, Indonesia and Maldives) or within a particular region (when comparing the different French Polynesian islands). However, this specificity was found to be locally limited when assessing different sampling sites from Rangiroa. Samples from the outer part of the reef of this atoll (Tiputa and Avatoru sites) were exposed to higher currents compared to the ones located within the lagoon and only relying on daily water exchanges through the passes (sites \u0026ldquo;Lagon\u0026rdquo; 1 and \u0026ldquo;Lagon 2\u0026rdquo;). Such differences in terms of oceanographic conditions do not appear to be involved in the variations of the prokaryotic communities. A similar result was observed for the sponge \u003cem\u003eCliona delitrix\u003c/em\u003e, characterized by a prokaryotic \u003cem\u003eβ\u003c/em\u003e-diversity only correlated to large geographical distances (\u0026gt;\u0026thinsp;1000 km), whereas no significant differences were observed more locally (\u0026lt;\u0026thinsp;300 km) (Easson et al. \u003cspan citationid=\"CR37\" class=\"CitationRef\"\u003e2020\u003c/span\u003e).\u003c/p\u003e \u003cp\u003eMany studies have provided insights on correlation between sponge phylogeny and microbiota. In some cases, dominant symbiont selection is believed to be related to sponge host identity (Easson and Thacker \u003cspan citationid=\"CR38\" class=\"CitationRef\"\u003e2014\u003c/span\u003e; O\u0026rsquo;Brien et al. \u003cspan citationid=\"CR91\" class=\"CitationRef\"\u003e2020\u003c/span\u003e). Several holobiont features like stress flexibility or immunity can be decisive in species adaptation, especially when considering future ocean scenarios with perturbed conditions (Posadas et al. \u003cspan citationid=\"CR101\" class=\"CitationRef\"\u003e2022\u003c/span\u003e). In this sense, tropical patterns in microbial communities amongst marine sponge orders have been revealed, most probably driven not only by vertical transmission but also by abiotic factors (Taylor et al. \u003cspan citationid=\"CR118\" class=\"CitationRef\"\u003e2007\u003c/span\u003e; Schmitt et al. \u003cspan citationid=\"CR109\" class=\"CitationRef\"\u003e2012b\u003c/span\u003e).\u003c/p\u003e \u003cp\u003eOur study revealed significant correlations between geographical and \u003cem\u003eβ\u003c/em\u003e-diversity distances, suggesting an influence of both environmental parameters and/or phylogenetic differences between the sponge samples. However, no clear geographical pattern was observed within the phylogenetic trees and no significant correlations were observed between the host phylogeny and the \u003cem\u003eβ\u003c/em\u003e-diversity distances. This result suggests that the prokaryotic community could rather be shaped by the environment than by the host specificity. Similar results were observed for \u003cem\u003eXestospongia\u003c/em\u003e spp. collected from the Indo-Pacific region, the Mozambique Channel, and the Red Sea, with a prokaryotic community more influenced by the biogeography than by the host-specificity (Swierts et al. \u003cspan citationid=\"CR117\" class=\"CitationRef\"\u003e2018\u003c/span\u003e). Other studies suggest a more important effect of the host-specificity (Meyer and Kuever \u003cspan citationid=\"CR81\" class=\"CitationRef\"\u003e2008\u003c/span\u003e; Busch et al. \u003cspan citationid=\"CR14\" class=\"CitationRef\"\u003e2022\u003c/span\u003e), especially when smaller geographical scales are considered (Cleary et al. \u003cspan citationid=\"CR25\" class=\"CitationRef\"\u003e2018\u003c/span\u003e), or when the phylogenetic comparisons are covering a large sponge diversity with higher taxonomic differences (i.e. family or order levels; Yang et al. \u003cspan citationid=\"CR133\" class=\"CitationRef\"\u003e2019\u003c/span\u003e). As studied for \u003cem\u003eC. delitrix\u003c/em\u003e, more insight on the potential effect of host-specificity should be considered at the level of the population genetics (Easson et al. \u003cspan citationid=\"CR37\" class=\"CitationRef\"\u003e2020\u003c/span\u003e). In our study, such level may justify why 10% of the \u003cem\u003eβ\u003c/em\u003e-diversity variance is explained exclusively by the geographical distance independently from the environmental factors. However, a large proportion of the variance associated with the environmental parameters, whether combined with the geographical distances (35%) or not (9%), confirm their predominant role driving the community structure.\u003c/p\u003e \u003cp\u003eWhen considering the relative effect of the different environmental parameters measured, the biogeographical study confirmed that temperature and salinity constitute two key abiotic factors shaping the \u003cem\u003eβ\u003c/em\u003e-diversity of the prokaryotic community. The effect of these parameters appears to be more important on the biogeographical study compared to the temporal one, probably in line with the greater variability of these factors through space rather than time. The temperature and salinity were found to have opposite effects on the community structure, explaining their high proportion of the shared variance. The negative relationship betwen temperature and salinity has been considered to play an important role within deep-sea sponge-holobionts (Busch et al. \u003cspan citationid=\"CR14\" class=\"CitationRef\"\u003e2022\u003c/span\u003e), as reported for \u003cem\u003eGeodia barretti\u003c/em\u003e (Strand et al. \u003cspan citationid=\"CR116\" class=\"CitationRef\"\u003e2017\u003c/span\u003e; Steffen et al. \u003cspan citationid=\"CR115\" class=\"CitationRef\"\u003e2022\u003c/span\u003e). In these studies, such parameters are consequently often interpreted carefully due to their collinearity with other depth-related factors, such as hydrostatic pressure or oxygen level. However, all samples in our study were collected from shallow waters with an average of 20m (SD: \u0026plusmn;5m), suggesting a limited role of the depth in the combined effect of temperature and salinity. Experimental studies have proven that an increase in temperature is one of the main factors impacting sponge holobionts health (Webster et al. \u003cspan citationid=\"CR128\" class=\"CitationRef\"\u003e2008\u003c/span\u003e; Pantile and Webster \u003cspan citationid=\"CR94\" class=\"CitationRef\"\u003e2011\u003c/span\u003e; Pita et al. \u003cspan citationid=\"CR100\" class=\"CitationRef\"\u003e2018\u003c/span\u003e; Slaby et al. \u003cspan citationid=\"CR113\" class=\"CitationRef\"\u003e2019\u003c/span\u003e). Nevertheless, it has also been observed that bacterial communities within \u003cem\u003eIrcinia\u003c/em\u003e spp. can be highly stable under large fluctuations of temperature occurring \u003cem\u003ein situ\u003c/em\u003e (Erwin et al. \u003cspan citationid=\"CR42\" class=\"CitationRef\"\u003e2012\u003c/span\u003e) or in controlled conditions (Pita et al. \u003cspan citationid=\"CR99\" class=\"CitationRef\"\u003e2013\u003c/span\u003e).\u003c/p\u003e \u003cp\u003eAs for salinity changes, only few studies have looked for specific effects on sponges, especially at the holobiont scale. Although decreasing salinity gradients were linked to a decline in growth of \u003cem\u003eCymbastela concentrica\u003c/em\u003e, the abundance of its photosymbionts was not significantly affected (Roberts et al. \u003cspan citationid=\"CR105\" class=\"CitationRef\"\u003e2006\u003c/span\u003e). Additionally, a stability of the host health and prokaryotic communities was demonstrated for diverse HMA and LMA species exposed to salinity disturbances over time (Glasl et al. \u003cspan citationid=\"CR44\" class=\"CitationRef\"\u003e2018\u003c/span\u003e).\u003c/p\u003e \u003c/div\u003e \u003cdiv id=\"Sec16\" class=\"Section2\"\u003e \u003ch2\u003eStability and role of the core community\u003c/h2\u003e \u003cp\u003eThe core concept within holobionts has been defined as the spatio-temporal stable part of the microbiome community, allowing to identify potential key symbionts displaying strong associations with their hosts (Astudillo-Garc\u0026iacute;a et al. \u003cspan citationid=\"CR4\" class=\"CitationRef\"\u003e2017\u003c/span\u003e; Risely \u003cspan citationid=\"CR104\" class=\"CitationRef\"\u003e2020\u003c/span\u003e; Neu et al. \u003cspan citationid=\"CR89\" class=\"CitationRef\"\u003e2021a\u003c/span\u003e). In particular, identifying the dynamics and stability of these core taxa is relevant to better understand the health and resilience of holobionts facing stress or damage, such as in a context of climate change or wild harvest for pharmaceutical valorization (Hutchins et al. \u003cspan citationid=\"CR58\" class=\"CitationRef\"\u003e2019\u003c/span\u003e; Risely \u003cspan citationid=\"CR104\" class=\"CitationRef\"\u003e2020\u003c/span\u003e; Neu et al. \u003cspan citationid=\"CR89\" class=\"CitationRef\"\u003e2021a\u003c/span\u003e). However, the parameters to identify the core microbiome can vary a lot depending on the studies, featuring a broad set of context-based parameters, and the lack of a consensus definition of the core-community, often leading to different interpretations (Astudillo-Garc\u0026iacute;a et al. \u003cspan citationid=\"CR4\" class=\"CitationRef\"\u003e2017\u003c/span\u003e; Astudillo-Garc\u0026iacute;a et al. \u003cspan citationid=\"CR3\" class=\"CitationRef\"\u003e2020\u003c/span\u003e). In our study, we defined the spatio-temporal core community of \u003cem\u003eD. metachromia\u003c/em\u003e by identifying core ASVs with a prevalence of 90% in all samples. This parameter allowed a conservative approach to keep only persistent ASVs across space and time, but also to compare our results with previous studies using this approach (Thomas et al. \u003cspan citationid=\"CR121\" class=\"CitationRef\"\u003e2016\u003c/span\u003e; Astudillo‐Garc\u0026iacute;a et al. 2020).\u003c/p\u003e \u003cp\u003eThe threshold returned 18 core ASVs, representing a relatively small set of the overall community in terms of average percentages of sequences (10.6%) and ASV numbers (6.6%). With a threshold of 85% of occurrence, the core microbiome described within 5 sponge species was found to vary from 5 to 20 core OTUs (Thomas et al. \u003cspan citationid=\"CR121\" class=\"CitationRef\"\u003e2016\u003c/span\u003e). Furthermore, a similar threshold (90% of occurrence, with an abundance threshold of 0.1%) returned 4 core OTUs for 20 host sponge species (Astudillo-Garc\u0026iacute;a et al. \u003cspan citationid=\"CR3\" class=\"CitationRef\"\u003e2020\u003c/span\u003e). Together with these results, our work suggests that a limited number of taxa are core members of sponge microbiomes. However, the number of core ASVs or OTUs depends on sampling size, especially since larger spatio-temporal scales tend to decrease the richness of the core microbiome (Turnbaugh et al. \u003cspan citationid=\"CR124\" class=\"CitationRef\"\u003e2007\u003c/span\u003e).\u003c/p\u003e \u003cp\u003eThe core bacterial community of \u003cem\u003eD. metachromia\u003c/em\u003e was found stable, with no significant changes of relative percentages of sequences and ASV numbers over time and space. Transplantation, thus, seems to have had a very limited impact on the structure of the potential key symbiotic community of the targeted sponge species. Also, the high growth rates of farmed explants could be linked to the overall stability of the core taxa, both indicating a steady fitness of the holobiont. These findings are in agreement with previous studies analyzing sponge-associated core microbiomes, either in tropical regions (Campana et al. \u003cspan citationid=\"CR19\" class=\"CitationRef\"\u003e2021b\u003c/span\u003e; Leal et al. \u003cspan citationid=\"CR66\" class=\"CitationRef\"\u003e2022\u003c/span\u003e), temperate waters (Lamb and Watts \u003cspan citationid=\"CR64\" class=\"CitationRef\"\u003e2023\u003c/span\u003e) or polar ecosystems (Happel et al. \u003cspan citationid=\"CR48\" class=\"CitationRef\"\u003e2022\u003c/span\u003e), suggesting that changing environmental conditions might have a relatively low impact in the shaping of the core sponge microbiome. The 18 core ASVs identified belong to 13 genera mostly affiliated to Acidobacteriota, Chloroflexi, PAUC34f and Proteobacteria. Such results, returning rather small globally-distributed taxa, are similar to previous works focusing on core communities defined from multiple species coming from very distinct locations (Schmitt et al. \u003cspan citationid=\"CR108\" class=\"CitationRef\"\u003e2012a\u003c/span\u003e, \u003cspan citationid=\"CR109\" class=\"CitationRef\"\u003eb\u003c/span\u003e).\u003c/p\u003e \u003cp\u003eThe core taxa of \u003cem\u003eD. metachromia\u003c/em\u003e within the phylum Acidobacteriota, were assigned to Vicinamibacterale, Thermoanaerobaculales and the uncultured group PAUC26f. Marine Thermoanaerobaculales are not strictly sponge-related, as they have also been extracted from molluscs (Neu et al. \u003cspan citationid=\"CR90\" class=\"CitationRef\"\u003e2021b\u003c/span\u003e; Banker et al. \u003cspan citationid=\"CR5\" class=\"CitationRef\"\u003e2022\u003c/span\u003e). Likewise, a core ASV was affiliated to the class bacteriap25 (Myxococcales; Proteobacteria), which is found in soils and terrestrial plants such as some non-indigenous Asteraceae that have been introduced in French Polynesia (Landwehr et al. \u003cspan citationid=\"CR65\" class=\"CitationRef\"\u003e2016\u003c/span\u003e; Monteiro \u003cspan citationid=\"CR86\" class=\"CitationRef\"\u003e2020\u003c/span\u003e).\u003c/p\u003e \u003cp\u003eOne of the main core ASVs belongs to the candidate phylum PAUC34f, which is described as part of the most important clades directly associated with the cycling of dissolved organic matter (DOM) within the Caribbean sponge \u003cem\u003ePlakortis angulospiculatus\u003c/em\u003e (Campana et al. \u003cspan citationid=\"CR18\" class=\"CitationRef\"\u003e2021a\u003c/span\u003e). Within the Nitrospinota phylum, a core ASV was affiliated to the candidate clade P9X2b3D02, which was recently also observed in \u003cem\u003eAplysina caissara\u003c/em\u003e restricted to southern Brazil (Hardoim et al. \u003cspan citationid=\"CR49\" class=\"CitationRef\"\u003e2021\u003c/span\u003e). This clade did not occur in the marine sediment sampled nearby or in other coastal sponge species from the Western Atlantic. Along with some Chloroflexi and Dadabacteria taxa, Nitrospinota bacteria are known to participate in the sulfur cycling and are capable of synthesizing various B-vitamins or essential amino acids (Bayer et al. \u003cspan citationid=\"CR8\" class=\"CitationRef\"\u003e2018\u003c/span\u003e; Engelberts et al. \u003cspan citationid=\"CR40\" class=\"CitationRef\"\u003e2020\u003c/span\u003e). Additionally, two classes (S085 and SAR202 clades) from the phylum Chloroflexi were represented within the core ASVs. Clade SAR202 has been observed in the water column of the aphotic zone and is believed to take part in the oxidation of sulfite particles (Mehrshad et al. \u003cspan citationid=\"CR80\" class=\"CitationRef\"\u003e2018\u003c/span\u003e). Recently, this clade was also found to dominate the microbiome of sponges living in extreme conditions, with large variation of pH values, temperature levels and oxygen concentrations in a semi-enclosed reef lagoon (Maggioni et al. \u003cspan citationid=\"CR74\" class=\"CitationRef\"\u003e2023b\u003c/span\u003e). Such dominance of SAR202 is suggested to play an important role for providing energy to the host through their DOM assimilation ability, thus helping with its adaptation in a changing environment (Maggioni et al. \u003cspan citationid=\"CR73\" class=\"CitationRef\"\u003e2023a\u003c/span\u003e, b). In respect to its potential role in DOM turnover and nutrient dynamics, especially in low-light environments, SAR202 might be providing a competitive advantage to the sponge for DOM assimilation and cycling (Bayer et al. \u003cspan citationid=\"CR8\" class=\"CitationRef\"\u003e2018\u003c/span\u003e; Morganti et al. \u003cspan citationid=\"CR87\" class=\"CitationRef\"\u003e2022\u003c/span\u003e).\u003c/p\u003e \u003cp\u003eFinally, the most abundant core ASVs belong to the order Dadabacteriales (Dadabacteria). Dadabacteria have been observed from the open ocean, from living systems related to the presence of organic carbon (including marine sponges) or from hydrothermal structures. This order notably includes obligate aerobic and photo-heterotrophic taxa; their ecological distribution throughout the water column could actually be driven by phototrophy (Graham and Tully \u003cspan citationid=\"CR46\" class=\"CitationRef\"\u003e2021\u003c/span\u003e). In sponges, they have only been recently discovered (Schuster et al. \u003cspan citationid=\"CR110\" class=\"CitationRef\"\u003e2016\u003c/span\u003e; Gavriilidou et al. \u003cspan citationid=\"CR43\" class=\"CitationRef\"\u003e2023\u003c/span\u003e), and their ecological involvements remain partially unaddressed. Yet, they appear to degrade microbial dissolved organic matter like phospholipids and complex carbohydrates, but also to produce important compounds for the host-interaction such as vitamins and cofactors (Hug et al. \u003cspan citationid=\"CR56\" class=\"CitationRef\"\u003e2016\u003c/span\u003e; Gavriilidou et al. \u003cspan citationid=\"CR43\" class=\"CitationRef\"\u003e2023\u003c/span\u003e).\u003c/p\u003e \u003c/div\u003e"},{"header":"Conclusion","content":"\u003cp\u003eEnvironmental changes often play a determinant role on the prokaryotic community of a sponge holobiont (Morrow et al. \u003cspan citationid=\"CR88\" class=\"CitationRef\"\u003e2016\u003c/span\u003e; Batista et al. \u003cspan citationid=\"CR7\" class=\"CitationRef\"\u003e2018\u003c/span\u003e; Busch et al. \u003cspan citationid=\"CR14\" class=\"CitationRef\"\u003e2022\u003c/span\u003e). Consequently, it appears essential to better investigate the factors that could shape sponge holobiont dynamics, especially when a commercial outlet for pharmaceutics is tested during farming trials.\u003c/p\u003e \u003cp\u003eAn overall stable \u003cem\u003eα\u003c/em\u003e- and \u003cem\u003eβ\u003c/em\u003e-diversity of the prokaryotic community was observed for the sponge \u003cem\u003eDactylospongia metachromia\u003c/em\u003e whose explants were farmed for 18 months on the outer reef of Rangiroa (Tuamotu, French Polynesia). These results are in line with the good viability, high survival rates and constant biomass gain of the new sponge generations observed \u003cem\u003ein situ.\u003c/em\u003e\u003c/p\u003e \u003cp\u003eHowever, some parameters such as temperature and salinity appear to partially impact a small fraction of the community. The moderated importance of these two abiotic factors was confirmed through a biogeographical study considering a broad regional scale (i.e. Red Sea and Indo-Pacific regions). Such a spatial investigation also revealed that both environmental variables are predominant over the effects of host-phylogeny differences. Hence, salinity and temperature need to be considered for their potential influence on the production of bioactive metabolites by the sponge-microorganism complex (Brinkmann et al. \u003cspan citationid=\"CR13\" class=\"CitationRef\"\u003e2017\u003c/span\u003e).\u003c/p\u003e \u003cp\u003eThe core microbiota of \u003cem\u003eD. metachromia\u003c/em\u003e is composed of bacterial phyla widely distributed among sponges (except the poorly characterized Dadabacteria phylum) and these taxa are known for their DOM assimilation abilities, sulfur cycling and vitamin production. Such findings highlight their global importance for the holobiont functioning. Such taxa do not seem to be strongly affected either by geographical and environmental characteristics or by the duration of the farming trials. Indeed, core communities of sponge holobionts commonly appear to be highly similar alongside a broad spatio-temporal gradient.\u003c/p\u003e \u003cp\u003eHowever, thermal changes have previously been linked to bacterial shifts in sponge-associated microbiome, global issues such as a general rise in sea temperatures are likely to affect such bacterial communities through the disappearance of certain symbionts for the benefit of others (Lemoine et al. \u003cspan citationid=\"CR68\" class=\"CitationRef\"\u003e2007\u003c/span\u003e; Webster et al. \u003cspan citationid=\"CR128\" class=\"CitationRef\"\u003e2008\u003c/span\u003e). In the context of climate change, it becomes ever more important nowadays to keep extending our knowledge of the potential effects of ocean variables on a large diversity of sponge holobionts such as the one targeted by this work, especially considering the stability and function of their core community. In addition, our findings are of interest in a growing context of marine-based chemical substances valuation, notably leading the path to a potential exploitation of this promising sponge species in other geographical regions.\u003c/p\u003e"},{"header":"Declarations","content":"\u003cp\u003eAcknowledgments\u003c/p\u003e\n\u003cp\u003eWe want to acknowledge the availability and assistance provided by the whole crew of the \u003cem\u003eAlis\u0026nbsp;\u003c/em\u003evessel, including all scientific members who helped with sample collection and analysis during the \u0026ldquo;TUAM\u0026rsquo;2018\u0026rdquo; oceanographic mission. We are also grateful to the DRM headquarters and field operators, either coming from Tahiti or based in Rangiroa, for their help during the set-up and monitoring of the aquaculture trials. Special thanks to all UMR 241 EIO lab technicians who assisted in the organization of the field missions and the preparation of general equipment. Lastly, we want to show our great appreciation to IRD researcher Dr. Sylvain Petek, who contributed to the general organization of the REDAME project and performed some key steps of the chemical analyses.\u003c/p\u003e\n\u003cp\u003eFunding\u0026nbsp;\u003c/p\u003e\n\u003cp\u003eMM was granted a PhD fellowship at the University of French Polynesia (French \u0026quot;Minist\u0026egrave;re de l\u0026apos;Enseignement Sup\u0026eacute;rieur, de la Recherche et de l\u0026rsquo;Innovation\u0026quot;). The oceanographic campaign was majoritarily funded by the IRD and their LEMAR lab team. Additional funds were allocated by the GENAVIR research vessel company and the UPF. Field campaigns in Rangiroa and environmental analyzes within French Polynesia were funded through the REDAME project, financially led by the \u0026quot;Direction des Ressources Marines\u0026quot; (grant number 08858/VP/DRMM) and the Government of French Polynesia (grant number 3299/MTF). The lab work was funded by the NWO-VIDI with project number 16.161.301.\u003c/p\u003e\n\u003cp\u003eData accessibility\u003c/p\u003e\n\u003cp\u003eCOI and 28S rRNA gene sequences were deposited and are publicly available in GenBank under the accession numbers OM980557-OM980598 (\u003cstrong\u003eTable S3\u003c/strong\u003e).\u003c/p\u003e\n\u003cp\u003e16S rRNA gene sequences were deposited and are publicly available in the NCBI Sequences Read Archive (SRA) under the BioProject ID PRJNA917856, accession number.\u003c/p\u003e\n\u003cp\u003eThe R script used for all the 16S rRNA gene metabarcoding analysis can be found at https://github.com/BenoitPAIX/Dactylospongia_metachromia_microbiome\u003c/p\u003e\n\u003cp\u003eAuthor contributions\u003c/p\u003e\n\u003cp\u003eCD, NGM and RH contributed to the initial design and set-up of the REDAME project which financially supported this study. CD and MM collected and prepared the environmental samples during the oceanic campaign in French Polynesian archipelagos. MM performed the installation of farming structures and the monitoring of the survival, growth and environmental data of the sponges in Rangiroa. As for the microbiological analyses, MM and NV conceived the study design and protocols. MM provided the samples from French Polynesia. NV and RA provided the samples from the Red Sea, the Maldives and Indonesia. NW and MM extracted the DNA. BP performed the library preparation for the 16S rRNA gene sequencing. NW and MM prepared the samples for the sponge barcoding. BP analyzed the 16S rRNA gene metabarcoding data. NW and MM analyzed the barcoding data. BP created the graphics. BP and MM prepared the first draft of the manuscript. All authors read, contributed to the revision, and approved the submitted version. \u003c/p\u003e"},{"header":"References","content":"\u003col\u003e\u003cli\u003e\u003cspan\u003eAguila-Ram\u0026iacute;rez RN, Hern\u0026aacute;ndez-Guerrero CJ, Gonz\u0026aacute;lez-Acosta B, et al (2014) Antifouling activity of symbiotic bacteria from sponge \u003cem\u003eAplysina gerardogreeni\u003c/em\u003e. International Biodeterioration \u0026amp; Biodegradation 90:64\u0026ndash;70. \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1016/j.ibiod.2014.02.003\u003c/span\u003e\u003cspan address=\"10.1016/j.ibiod.2014.02.003\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eAminot A, Rey F (2001) Chlorophyll a: Determination by spectroscopic methods. International Council for the Exploration of the Sea (ICES)\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eAstudillo-Garc\u0026iacute;a C, Bell JJ, Montoya JM, et al (2020) Assessing the strength and sensitivity of the core microbiota approach on a highly diverse sponge reef. Environ Microbiol 22:3985\u0026ndash;3999. \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1111/1462-2920.15185\u003c/span\u003e\u003cspan address=\"10.1111/1462-2920.15185\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eAstudillo-Garc\u0026iacute;a C, Bell JJ, Webster NS, et al (2017) Evaluating the core microbiota in complex communities: A systematic investigation. Environ Microbiol 19:1450\u0026ndash;1462. \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1111/1462-2920.13647\u003c/span\u003e\u003cspan address=\"10.1111/1462-2920.13647\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eBanker RMW, Lipovac J, Stachowicz JJ, Gold DA (2022) Sodium molybdate does not inhibit sulfate-reducing bacteria but increases shell growth in the Pacific oyster \u003cem\u003eMagallana gigas\u003c/em\u003e. PLoS One 17:e0262939. \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1371/journal.pone.0262939\u003c/span\u003e\u003cspan address=\"10.1371/journal.pone.0262939\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eBaquiran JIP, Conaco C (2018) Sponge-microbe partnerships are stable under eutrophication pressure from mariculture. Marine Pollution Bulletin 136:125\u0026ndash;134. \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1016/j.marpolbul.2018.09.011\u003c/span\u003e\u003cspan address=\"10.1016/j.marpolbul.2018.09.011\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eBatista D, Costa R, Carvalho AP, et al (2018) Environmental conditions affect activity and associated microorganisms of marine sponges. Marine Environmental Research 142:59\u0026ndash;68. \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1016/j.marenvres.2018.09.020\u003c/span\u003e\u003cspan address=\"10.1016/j.marenvres.2018.09.020\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eBayer K, Jahn MT, Slaby BM, et al (2018) Marine sponges as Chloroflexi hot spots: Genomic insights and high-resolution visualization of an abundant and diverse symbiotic clade. mSystems 3:e00150-18. \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1128/mSystems.00150-18\u003c/span\u003e\u003cspan address=\"10.1128/mSystems.00150-18\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eBell JJ (2008) The functional roles of marine sponges. Estuarine, Coastal and Shelf Science 79:341\u0026ndash;353. \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1016/j.ecss.2008.05.002\u003c/span\u003e\u003cspan address=\"10.1016/j.ecss.2008.05.002\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eBibi F, Yasir M, Al-Sofyani A, et al (2020) Antimicrobial activity of bacteria from marine sponge Suberea mollis and bioactive metabolites of \u003cem\u003eVibrio\u003c/em\u003e sp. EA348. Saudi J Biol Sci 27:1139\u0026ndash;1147. \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1016/j.sjbs.2020.02.002\u003c/span\u003e\u003cspan address=\"10.1016/j.sjbs.2020.02.002\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eBierwirth J, Mantas TP, Villechanoux J, Cerrano C (2022) Restoration of marine sponges\u0026mdash;what can we learn from over a century of experimental cultivation? Water 14:1055. \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.3390/w14071055\u003c/span\u003e\u003cspan address=\"10.3390/w14071055\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eBorcard D, Legendre P, Drapeau P (1992) Partialling out the Spatial Component of Ecological Variation. Ecology 73:1045\u0026ndash;1055. \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.2307/1940179\u003c/span\u003e\u003cspan address=\"10.2307/1940179\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eBrinkmann CM, Marker A, Kurtb\u0026ouml;ke Dİ (2017) An Overview on Marine Sponge-Symbiotic Bacteria as Unexhausted Sources for Natural Product Discovery. Diversity 9:40. \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.3390/d9040040\u003c/span\u003e\u003cspan address=\"10.3390/d9040040\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eBusch K, Slaby BM, Bach W, et al (2022) Biodiversity, environmental drivers, and sustainability of the global deep-sea sponge microbiome. Nat Commun 13:5160. \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1038/s41467-022-32684-4\u003c/span\u003e\u003cspan address=\"10.1038/s41467-022-32684-4\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eCallahan BJ, McMurdie PJ, Holmes SP (2017) Exact sequence variants should replace operational taxonomic units in marker-gene data analysis. ISME J 11:2639\u0026ndash;2643. \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1038/ismej.2017.119\u003c/span\u003e\u003cspan address=\"10.1038/ismej.2017.119\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eCallahan BJ, McMurdie PJ, Rosen MJ, et al (2016a) DADA2: High-resolution sample inference from Illumina amplicon data. Nat Methods 13:581\u0026ndash;583. \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1038/nmeth.3869\u003c/span\u003e\u003cspan address=\"10.1038/nmeth.3869\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eCallahan BJ, Sankaran K, Fukuyama JA, et al (2016b) Bioconductor Workflow for Microbiome Data Analysis: from raw reads to community analyses\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eCampana S, Busch K, Hentschel U, et al (2021a) DNA-stable isotope probing (DNA-SIP) identifies marine sponge-associated bacteria actively utilizing dissolved organic matter (DOM). Environ Microbiol 23:4489\u0026ndash;4504. \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1111/1462-2920.15642\u003c/span\u003e\u003cspan address=\"10.1111/1462-2920.15642\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eCampana S, Demey C, Busch K, et al (2021b) Marine sponges maintain stable bacterial communities between reef sites with different coral to algae cover ratios. FEMS Microbiol Ecol 97:fiab115. \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1093/femsec/fiab115\u003c/span\u003e\u003cspan address=\"10.1093/femsec/fiab115\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eCapon RJ, MacLeod JK (1987) Revision of the absolute stereochemistry of ilimaquinone. The Journal of Organic Chemistry 52:5059\u0026ndash;5060\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eC\u0026aacute;rdenas CA, Bell JJ, Davy SK, et al (2014) Influence of environmental variation on symbiotic bacterial communities of two temperate sponges. FEMS Microbiology Ecology 88:516\u0026ndash;527. \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1111/1574-6941.12317\u003c/span\u003e\u003cspan address=\"10.1111/1574-6941.12317\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eCarte B, Rose CB, Faulkner DJ (1985) 5-Epi-Ilimaquinone, a metabolite of the sponge \u003cem\u003eFenestraspongia\u003c/em\u003e sp. The Journal of Organic Chemistry 50:2785\u0026ndash;2787\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eChen ML, Becraft ED, Pachiadaki M, et al (2020) Hiding in plain sight: The globally distributed bacterial candidate Phylum PAUC34f. Front Microbiol 11:376. \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.3389/fmicb.2020.00376\u003c/span\u003e\u003cspan address=\"10.3389/fmicb.2020.00376\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eChombard C, Boury-Esnault N, Tillier S (1998) Reassessment of homology of morphological characters in tetractinellid sponges based on molecular data. Systematic Biology 47:351\u0026ndash;366. \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1080/106351598260761\u003c/span\u003e\u003cspan address=\"10.1080/106351598260761\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eCleary DFR, Pol\u0026oacute;nia ARM, de Voogd NJ (2018) Prokaryote composition and predicted metagenomic content of two \u003cem\u003eCinachyrella\u003c/em\u003e Morphospecies and water from West Papuan Marine Lakes. FEMS Microbiology Ecology 94:. \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1093/femsec/fix175\u003c/span\u003e\u003cspan address=\"10.1093/femsec/fix175\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eCleary DFR, Pol\u0026oacute;nia ARM, Reijnen BT, et al (2020) Prokaryote communities inhabiting endemic and newly discovered sponges and octocorals from the Red Sea. Microb Ecol 80:103\u0026ndash;119. \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1007/s00248-019-01465-w\u003c/span\u003e\u003cspan address=\"10.1007/s00248-019-01465-w\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eCleary DFR, Pol\u0026oacute;nia ARM, Swierts T, et al (2022) Spatial and environmental variables structure sponge symbiont communities. Molecular Ecology 31:4932\u0026ndash;4948. \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1111/mec.16631\u003c/span\u003e\u003cspan address=\"10.1111/mec.16631\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eDaletos G, de Voogd NJ, M\u0026uuml;ller WEG, et al (2014) Cytotoxic and protein kinase inhibiting nakijiquinones and nakijiquinols from the sponge \u003cem\u003eDactylospongia metachromia\u003c/em\u003e. J Nat Prod 77:218\u0026ndash;226. \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1021/np400633m\u003c/span\u003e\u003cspan address=\"10.1021/np400633m\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eDavis NM, Proctor DM, Holmes SP, et al (2018) Simple statistical identification and removal of contaminant sequences in marker-gene and metagenomics data. Microbiome 6:226. \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1186/s40168-018-0605-2\u003c/span\u003e\u003cspan address=\"10.1186/s40168-018-0605-2\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eDe Caralt S, S\u0026aacute;nchez-Fontenla J, Uriz MJ, Wijffels RH (2010) \u003cem\u003eIn situ\u003c/em\u003e aquaculture methods for \u003cem\u003eDysidea avara\u003c/em\u003e (Demospongiae, Porifera) in the northwestern mediterranean. Marine Drugs 8:1731\u0026ndash;1742. \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.3390/md8061731\u003c/span\u003e\u003cspan address=\"10.3390/md8061731\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003ede Goeij JM, Lesser MP, Pawlik JR (2017) Nutrient fluxes and ecological functions of coral reef sponges in a changing ocean. In: Carballo JL, Bell JJ (eds) Climate Change, Ocean Acidification and Sponges: Impacts Across Multiple Levels of Organization. Springer International Publishing, Cham, pp 373\u0026ndash;410\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003ede Voogd NJ (2007) The mariculture potential of the Indonesian reef-dwelling sponge \u003cem\u003eCallyspongia\u003c/em\u003e (\u003cem\u003eEuplacella\u003c/em\u003e) \u003cem\u003ebiru\u003c/em\u003e: Growth, survival and bioactive compounds. Aquaculture 262:54\u0026ndash;64. \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1016/j.aquaculture.2006.09.028\u003c/span\u003e\u003cspan address=\"10.1016/j.aquaculture.2006.09.028\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eDray S, Dufour A-B (2007) The ade4 package: implementing the duality diagram for ecologists. J Stat Softw 22:1\u0026ndash;20. \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.18637/jss.v022.i04\u003c/span\u003e\u003cspan address=\"10.18637/jss.v022.i04\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eDuckworth A (2009) Farming sponges to supply bioactive metabolites and bath sponges: a review. Mar Biotechnol 11:669\u0026ndash;679. \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1007/s10126-009-9213-2\u003c/span\u003e\u003cspan address=\"10.1007/s10126-009-9213-2\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eDuckworth A, Battershill C (2003a) Sponge aquaculture for the production of biologically active metabolites: the influence of farming protocols and environment. Aquaculture 221:311\u0026ndash;329. \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1016/S0044-8486(03)00070-X\u003c/span\u003e\u003cspan address=\"10.1016/S0044-8486(03)00070-X\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eDuckworth AR, Battershill CN (2003b) Developing farming structures for production of biologically active sponge metabolites. Aquaculture 217:139\u0026ndash;156. \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1016/S0044-8486(02)00038-8\u003c/span\u003e\u003cspan address=\"10.1016/S0044-8486(02)00038-8\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eEasson CG, Chaves-Fonnegra A, Thacker RW, Lopez JV (2020) Host population genetics and biogeography structure the microbiome of the sponge \u003cem\u003eCliona delitrix\u003c/em\u003e. Ecology and Evolution 10:2007. \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1002/ece3.6033\u003c/span\u003e\u003cspan address=\"10.1002/ece3.6033\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eEasson CG, Thacker RW (2014) Phylogenetic signal in the community structure of host-specific microbiomes of tropical marine sponges. Front Microbiol 5:532. \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.3389/fmicb.2014.00532\u003c/span\u003e\u003cspan address=\"10.3389/fmicb.2014.00532\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eEhrlich H, Worch H (2007) Sponges as natural composites: from biomimetic potential to development of new biomaterials. In: In: Cust\u0026oacute;dio MR, L\u0026ocirc;bo-Hajdu G, Hajdu E, Muricy G (eds) Porifera Research. Biodiversity, Innovation and Sustainability. Livros de Museu Nacional 28, Rio de Janeiro, pp 303\u0026ndash;312\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eEngelberts JP, Robbins SJ, de Goeij JM, et al (2020) Characterization of a sponge microbiome using an integrative genome-centric approach. ISME J 14:1100\u0026ndash;1110. \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1038/s41396-020-0591-9\u003c/span\u003e\u003cspan address=\"10.1038/s41396-020-0591-9\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eErwin PM, Coma R, L\u0026oacute;pez-Sendino P, et al (2015) Stable symbionts across the HMA-LMA dichotomy: low seasonal and interannual variation in sponge-associated bacteria from taxonomically diverse hosts. FEMS Microbiol Ecol 91:fiv115. \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1093/femsec/fiv115\u003c/span\u003e\u003cspan address=\"10.1093/femsec/fiv115\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eErwin PM, Pita L, L\u0026oacute;pez-Legentil S, Turon X (2012) Stability of sponge-associated bacteria over large seasonal shifts in temperature and irradiance. Applied and Environmental Microbiology 78:7358\u0026ndash;7368. \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1128/AEM.02035-12\u003c/span\u003e\u003cspan address=\"10.1128/AEM.02035-12\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eGavriilidou A, Avcı B, Galani A, et al (2023) Candidatus Nemesobacterales is a sponge-specific clade of the candidate phylum Desulfobacterota adapted to a symbiotic lifestyle. ISME J 17:1808\u0026ndash;1818. \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1038/s41396-023-01484-z\u003c/span\u003e\u003cspan address=\"10.1038/s41396-023-01484-z\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eGlasl B, Smith CE, Bourne DG, Webster NS (2018) Exploring the diversity-stability paradigm using sponge microbial communities. Sci Rep 8:8425. \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1038/s41598-018-26641-9\u003c/span\u003e\u003cspan address=\"10.1038/s41598-018-26641-9\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eGloor GB, Macklaim JM, Pawlowsky-Glahn V, Egozcue JJ (2017) Microbiome datasets are compositional: and this is not optional. Front Microbiol 8:2224. \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.3389/fmicb.2017.02224\u003c/span\u003e\u003cspan address=\"10.3389/fmicb.2017.02224\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eGraham ED, Tully BJ (2021) Marine Dadabacteria exhibit genome streamlining and phototrophy-driven niche partitioning. ISME J 15:1248\u0026ndash;1256. \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1038/s41396-020-00834-5\u003c/span\u003e\u003cspan address=\"10.1038/s41396-020-00834-5\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eHadas E, Shpigel M, Ilan M (2005) Sea ranching of the marine sponge \u003cem\u003eNegombata magnifica\u003c/em\u003e (Demospongiae, Latrunculiidae) as a first step for latrunculin B mass production. Aquaculture 244:159\u0026ndash;169. \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1016/j.aquaculture.2004.11.052\u003c/span\u003e\u003cspan address=\"10.1016/j.aquaculture.2004.11.052\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eHappel L, Rondon R, Font A, et al (2022) Stability of the microbiome of the sponge \u003cem\u003eMycale\u003c/em\u003e (\u003cem\u003eOxymycale\u003c/em\u003e) \u003cem\u003eacerata\u003c/em\u003e in the Western Antarctic Peninsula. Front Microbiol 13:827863. \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.3389/fmicb.2022.827863\u003c/span\u003e\u003cspan address=\"10.3389/fmicb.2022.827863\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eHardoim CCP, Ramaglia ACM, L\u0026ocirc;bo-Hajdu G, Cust\u0026oacute;dio MR (2021) Community composition and functional prediction of prokaryotes associated with sympatric sponge species of southwestern Atlantic coast. Sci Rep 11:9576. \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1038/s41598-021-88288-3\u003c/span\u003e\u003cspan address=\"10.1038/s41598-021-88288-3\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eHarvey AL, Edrada-Ebel R, Quinn RJ (2015) The re-emergence of natural products for drug discovery in the genomics era. Nat Rev Drug Discov 14:111\u0026ndash;129. \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1038/nrd4510\u003c/span\u003e\u003cspan address=\"10.1038/nrd4510\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eHentschel U, Hopke J, Horn M, et al (2002) Molecular evidence for a uniform microbial community in sponges from different oceans. Appl Environ Microbiol 68:4431\u0026ndash;4440. \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1128/AEM.68.9.4431-4440.2002\u003c/span\u003e\u003cspan address=\"10.1128/AEM.68.9.4431-4440.2002\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eHill M, Hill A, Lopez N, Harriott O (2006) Sponge-specific bacterial symbionts in the Caribbean sponge, \u003cem\u003eChondrilla nucula\u003c/em\u003e (Demospongiae, Chondrosida). Marine Biology 148:1221\u0026ndash;1230. \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1007/s00227-005-0164-5\u003c/span\u003e\u003cspan address=\"10.1007/s00227-005-0164-5\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eHitora Y, Sejiyama A, Honda K, et al (2021) Fluorescent image-based high-content screening of extracts of natural resources for cell cycle inhibitors and identification of a new sesquiterpene quinone from the sponge, \u003cem\u003eDactylospongia metachromia\u003c/em\u003e. Bioorg Med Chem 31:115968. \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1016/j.bmc.2020.115968\u003c/span\u003e\u003cspan address=\"10.1016/j.bmc.2020.115968\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eHooper JNA, Van Soest RWM, Willenz P (eds) (2002) Systema Porifera. Springer US, Boston, MA\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eHuelsenbeck JP, Ronquist F, Nielsen R, Bollback JP (2001) Bayesian Inference of Phylogeny and Its Impact on Evolutionary Biology. Science 294:2310\u0026ndash;2314. \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1126/science.1065889\u003c/span\u003e\u003cspan address=\"10.1126/science.1065889\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eHug LA, Thomas BC, Sharon I, et al (2016) Critical biogeochemical functions in the subsurface are associated with bacteria from new phyla and little studied lineages. Environ Microbiol 18:159\u0026ndash;173. \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1111/1462-2920.12930\u003c/span\u003e\u003cspan address=\"10.1111/1462-2920.12930\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eHunting ER, Franken O, Knopperts F, et al (2013) Substrate as a driver of sponge distributions in mangrove ecosystems. Marine Ecology Progress Series 486:133\u0026ndash;141. \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.3354/meps10376\u003c/span\u003e\u003cspan address=\"10.3354/meps10376\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eHutchins DA, Jansson JK, Remais JV, et al (2019) Climate change microbiology \u0026mdash; problems and perspectives. Nat Rev Microbiol 17:391\u0026ndash;396. \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1038/s41579-019-0178-5\u003c/span\u003e\u003cspan address=\"10.1038/s41579-019-0178-5\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eIndraningrat AAG, Smidt H, Sipkema D (2016) Bioprospecting sponge-associated microbes for antimicrobial compounds. Mar Drugs 14:87. \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.3390/md14050087\u003c/span\u003e\u003cspan address=\"10.3390/md14050087\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eJeong IH, Park JS (2013) Bacterial diversity of the South Pacific sponge, \u003cem\u003eDactylospongia metachromia\u003c/em\u003e based on DGGE fingerprinting. The Korean Journal of Microbiology\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eKatoh K, Misawa K, Kuma K, Miyata T (2002) MAFFT: a novel method for rapid multiple sequence alignment based on fast Fourier transform. Nucleic Acids Res 30:3059\u0026ndash;3066. \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1093/nar/gkf436\u003c/span\u003e\u003cspan address=\"10.1093/nar/gkf436\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eKochanowska AJ, Rao KV, Childress S, et al (2008) Secondary metabolites from three Florida sponges with antidepressant activity. J Nat Prod 71:186\u0026ndash;189. \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1021/np070371u\u003c/span\u003e\u003cspan address=\"10.1021/np070371u\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eLahti L, Shetty S, Blake T, Salojarvi J (2017) Tools for microbiome analysis in R. Version 1:28\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eLamb CE, Watts JEM (2023) Microbiome species diversity and seasonal stability of two temperate marine sponges \u003cem\u003eHymeniacidon perlevis\u003c/em\u003e and \u003cem\u003eSuberites massa\u003c/em\u003e. Environmental Microbiome 18:52. \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1186/s40793-023-00508-7\u003c/span\u003e\u003cspan address=\"10.1186/s40793-023-00508-7\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eLandwehr W, Wolf C, Wink J (2016) Actinobacteria and Myxobacteria-two of the most important bacterial resources for novel antibiotics. Curr Top Microbiol Immunol 398:273\u0026ndash;302. \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1007/82_2016_503\u003c/span\u003e\u003cspan address=\"10.1007/82_2016_503\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eLeal CV, Avelino-Alves D, Salazar V, et al (2022) Sponges present a core prokaryotic community stable across Tropical Western Atlantic. Science of The Total Environment 835:155145. \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1016/j.scitotenv.2022.155145\u003c/span\u003e\u003cspan address=\"10.1016/j.scitotenv.2022.155145\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eLeary D, Vierros M, Hamon G, et al (2009) Marine genetic resources: A review of scientific and commercial interest. Marine Policy 33:183\u0026ndash;194. \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1016/j.marpol.2008.05.010\u003c/span\u003e\u003cspan address=\"10.1016/j.marpol.2008.05.010\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eLemoine N, Buell N, Hill A, Hill M (2007) Assessing the utility of sponge microbial symbiont communities as models to study global climate change: a case study with \u003cem\u003eHalichondria bowerbanki\u003c/em\u003e. In: In: Cust\u0026oacute;dio MR, L\u0026ocirc;bo-Hajdu G, Hajdu E, Muricy G (eds) Porifera Research. Biodiversity, Innovation and Sustainability. Livros de Museu Nacional 28, Rio de Janeiro, pp 419\u0026ndash;425\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eLesser MP, Sabrina Pankey M, Slattery M, et al (2022) Microbiome diversity and metabolic capacity determines the trophic ecology of the holobiont in Caribbean sponges. ISME COMMUN 2:1\u0026ndash;12. \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1038/s43705-022-00196-3\u003c/span\u003e\u003cspan address=\"10.1038/s43705-022-00196-3\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eLim SC, Putchakarn S, Thai M-Q, et al (2016) Inventory of sponge fauna from the Singapore Strait to Taiwan Strait along the western coastline of the South China Sea\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eLowe B, Venkatesan J, Ehrlich H, Kim S-K (2016) Global Constraints, Prospects, and Perspectives of Marine Sponge Research. In: Pallela R, Ehrlich H (eds) Marine Sponges: Chemicobiological and Biomedical Applications. Springer India, New Delhi, pp 25\u0026ndash;35\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eLu P-H, Chueh S-C, Kung F-L, et al (2007) Ilimaquinone, a marine sponge metabolite, displays anticancer activity via GADD153-mediated pathway. European Journal of Pharmacology 556:45\u0026ndash;54. \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1016/j.ejphar.2006.10.061\u003c/span\u003e\u003cspan address=\"10.1016/j.ejphar.2006.10.061\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eMaggioni F, Bell JJ, Pujo-Pay M, et al (2023a) Sponge organic matter recycling: Reduced detritus production under extreme environmental conditions. Marine Pollution Bulletin 190:114869. \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1016/j.marpolbul.2023.114869\u003c/span\u003e\u003cspan address=\"10.1016/j.marpolbul.2023.114869\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eMaggioni F, Stenger P-L, Jourand P, Majorel C (2023b) The phylum Chloroflexi and their SAR202 clade dominate the microbiome of two marine sponges living in extreme environmental conditions. Marine Ecology 44:e12757. \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1111/maec.12757\u003c/span\u003e\u003cspan address=\"10.1111/maec.12757\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eMalve H (2016) Exploring the ocean for new drug developments: Marine pharmacology. J Pharm Bioallied Sci 8:83\u0026ndash;91. \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.4103/0975-7406.171700\u003c/span\u003e\u003cspan address=\"10.4103/0975-7406.171700\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eMaslin M (2022) \u0026Eacute;tude de la ressource en \u0026eacute;ponge \u003cem\u003eDactylospongia metachromia\u003c/em\u003e pour une production durable. Theses, Universit\u0026eacute; de la Polyn\u0026eacute;sie fran\u0026ccedil;aise\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eMaslin M, Gaertner-Mazouni N, Debitus C, et al (2021) Marine sponge aquaculture towards drug development: An ongoing history of technical, ecological, chemical considerations and challenges. Aquaculture Reports 21:100813. \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1016/j.aqrep.2021.100813\u003c/span\u003e\u003cspan address=\"10.1016/j.aqrep.2021.100813\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eMcMurdie PJ, Holmes S (2013) phyloseq: an R package for reproducible interactive analysis and graphics of microbiome census data. PLoS One 8:e61217. \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1371/journal.pone.0061217\u003c/span\u003e\u003cspan address=\"10.1371/journal.pone.0061217\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eMcMurdie PJ, Holmes S (2014) Waste not, want not: why rarefying microbiome data is inadmissible. PLoS Comput Biol 10:e1003531. \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1371/journal.pcbi.1003531\u003c/span\u003e\u003cspan address=\"10.1371/journal.pcbi.1003531\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eMehrshad M, Rodriguez-Valera F, Amoozegar MA, et al (2018) The enigmatic SAR202 cluster up close: shedding light on a globally distributed dark ocean lineage involved in sulfur cycling. ISME J 12:655\u0026ndash;668. \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1038/s41396-017-0009-5\u003c/span\u003e\u003cspan address=\"10.1038/s41396-017-0009-5\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eMeyer B, Kuever J (2008) Phylogenetic Diversity and Spatial Distribution of the Microbial Community Associated with the Caribbean Deep-water Sponge \u003cem\u003ePolymastia\u003c/em\u003e cf. \u003cem\u003ecorticata\u003c/em\u003e by 16S rRNA, aprA, and amoA Gene Analysis. Microb Ecol 56:306\u0026ndash;321. \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1007/s00248-007-9348-5\u003c/span\u003e\u003cspan address=\"10.1007/s00248-007-9348-5\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eMeyer CP, Geller JB, Paulay G (2005) Fine scale endemism on coral reefs: archipelagic differentiation in turbinid gastropods. Evolution 59:113\u0026ndash;125\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eMohamed NM, Enticknap JJ, Lohr JE, et al (2008a) Changes in bacterial communities of the marine sponge \u003cem\u003eMycale laxissima\u003c/em\u003e on transfer into aquaculture. Applied and Environmental Microbiology 74:1209\u0026ndash;1222. \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1128/AEM.02047-07\u003c/span\u003e\u003cspan address=\"10.1128/AEM.02047-07\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eMohamed NM, Rao V, Hamann MT, et al (2008b) Monitoring bacterial diversity of the marine sponge \u003cem\u003eIrcinia strobilina\u003c/em\u003e upon transfer into aquaculture. Applied and Environmental Microbiology 74:4133\u0026ndash;4143. \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1128/AEM.00454-08\u003c/span\u003e\u003cspan address=\"10.1128/AEM.00454-08\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eMoitinho-Silva L, Steinert G, Nielsen S, et al (2017) Predicting the HMA-LMA status in marine sponges by machine learning. Front Microbiol 8:752. \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.3389/fmicb.2017.00752\u003c/span\u003e\u003cspan address=\"10.3389/fmicb.2017.00752\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eMonteiro LCP (2020) Rhizospheric microbial communities of plants coexisting under different edaphic conditions\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eMorganti TM, Slaby BM, de Kluijver A, et al (2022) Giant sponge grounds of Central Arctic seamounts are associated with extinct seep life. Nat Commun 13:638. \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1038/s41467-022-28129-7\u003c/span\u003e\u003cspan address=\"10.1038/s41467-022-28129-7\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eMorrow KM, Fiore CL, Lesser MP (2016) Environmental drivers of microbial community shifts in the giant barrel sponge, \u003cem\u003eXestospongia muta\u003c/em\u003e, over a shallow to mesophotic depth gradient. Environmental Microbiology 18:2025\u0026ndash;2038. \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1111/1462-2920.13226\u003c/span\u003e\u003cspan address=\"10.1111/1462-2920.13226\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eNeu AT, Allen EE, Roy K (2021a) Defining and quantifying the core microbiome: Challenges and prospects. Proceedings of the National Academy of Sciences 118:e2104429118. \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1073/pnas.2104429118\u003c/span\u003e\u003cspan address=\"10.1073/pnas.2104429118\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eNeu AT, Torchin ME, Allen EE, Roy K (2021b) Microbiome divergence of marine gastropod species separated by the Isthmus of Panama. 2021.07.08.451645\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eO\u0026rsquo;Brien PA, Tan S, Yang C, et al (2020) Diverse coral reef invertebrates exhibit patterns of phylosymbiosis. ISME J 14:2211\u0026ndash;2222. \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1038/s41396-020-0671-x\u003c/span\u003e\u003cspan address=\"10.1038/s41396-020-0671-x\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eOksanen J, Blanchet FG, Friendly M, et al (2019) vegan: Community ecology package\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eOsinga R, Sidri M, Cerig E, et al (2010) Sponge aquaculture trials in the East-Mediterranean Sea: new approaches to earlier ideas. The Open Marine Biology Journal 4:\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003ePantile R, Webster N (2011) Strict thermal threshold identified by quantitative PCR in the sponge \u003cem\u003eRhopaloeides odorabile\u003c/em\u003e. Marine Ecology Progress Series 431:97\u0026ndash;105. \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.3354/meps09128\u003c/span\u003e\u003cspan address=\"10.3354/meps09128\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eParada AE, Needham DM, Fuhrman JA (2016) Every base matters: assessing small subunit rRNA primers for marine microbiomes with mock communities, time series and global field samples. Environ Microbiol 18:1403\u0026ndash;1414. \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1111/1462-2920.13023\u003c/span\u003e\u003cspan address=\"10.1111/1462-2920.13023\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003ePawlik JR, McMurray SE (2020) The emerging ecological and biogeochemical importance of sponges on coral reefs. Annual Review of Marine Science 12:315\u0026ndash;337. \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1146/annurev-marine-010419-010807\u003c/span\u003e\u003cspan address=\"10.1146/annurev-marine-010419-010807\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003ePereira DM, Valent\u0026atilde;o P, Andrade PB (2013) Chap. 7 - Lessons from the Sea: Distribution, SAR, and Molecular Mechanisms of Anti-inflammatory Drugs from Marine Organisms. In: Atta-ur-Rahman (ed) Studies in Natural Products Chemistry. Elsevier, pp 205\u0026ndash;228\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003ePiel J, Hui D, Wen G, et al (2004) Antitumor polyketide biosynthesis by an uncultivated bacterial symbiont of the marine sponge \u003cem\u003eTheonella swinhoei\u003c/em\u003e. Proceedings of the National Academy of Sciences 101:16222\u0026ndash;16227. \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1073/pnas.0405976101\u003c/span\u003e\u003cspan address=\"10.1073/pnas.0405976101\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003ePita L, Erwin PM, Turon X, L\u0026oacute;pez-Legentil S (2013) Till Death Do Us Part: Stable Sponge-Bacteria Associations under Thermal and Food Shortage Stresses. PLOS ONE 8:e80307. \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1371/journal.pone.0080307\u003c/span\u003e\u003cspan address=\"10.1371/journal.pone.0080307\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003ePita L, Rix L, Slaby BM, et al (2018) The sponge holobiont in a changing ocean: from microbes to ecosystems. Microbiome 6:46. \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1186/s40168-018-0428-1\u003c/span\u003e\u003cspan address=\"10.1186/s40168-018-0428-1\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003ePosadas N, Baquiran JIP, Nada MAL, et al (2022) Microbiome diversity and host immune functions influence survivorship of sponge holobionts under future ocean conditions. ISME J 16:58\u0026ndash;67. \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1038/s41396-021-01050-5\u003c/span\u003e\u003cspan address=\"10.1038/s41396-021-01050-5\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eQuast C, Pruesse E, Yilmaz P, et al (2013) The SILVA ribosomal RNA gene database project: improved data processing and web-based tools. Nucleic Acids Res 41:D590\u0026ndash;D596. \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1093/nar/gks1219\u003c/span\u003e\u003cspan address=\"10.1093/nar/gks1219\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eRadeke HS, Digits CA, Casaubon RL, Snapper ML (1999) Interactions of (\u0026ndash;)-ilimaquinone with methylation enzymes: implications for vesicular-mediated secretion. Chemistry \u0026amp; Biology 6:639\u0026ndash;647. \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1016/S1074-5521(99)80115-X\u003c/span\u003e\u003cspan address=\"10.1016/S1074-5521(99)80115-X\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eRisely A (2020) Applying the core microbiome to understand host\u0026ndash;microbe systems. Journal of Animal Ecology 89:1549\u0026ndash;1558. \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1111/1365-2656.13229\u003c/span\u003e\u003cspan address=\"10.1111/1365-2656.13229\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eRoberts DE, Davis AR, Cummins SP (2006) Experimental manipulation of shade, silt, nutrients and salinity on the temperate reef sponge \u003cem\u003eCymbastela concentrica\u003c/em\u003e. Marine Ecology Progress Series 307:143\u0026ndash;154. \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.3354/meps307143\u003c/span\u003e\u003cspan address=\"10.3354/meps307143\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eRuiz C, Valderrama K, Zea S, Castellanos L (2013) Mariculture and natural production of the antitumoural (+)-discodermolide by the Caribbean marine sponge \u003cem\u003eDiscodermia dissoluta\u003c/em\u003e. Mar Biotechnol 15:571\u0026ndash;583. \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1007/s10126-013-9510-7\u003c/span\u003e\u003cspan address=\"10.1007/s10126-013-9510-7\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eSchiefenh\u0026ouml;vel K, Kunzmann A (2012) Sponge farming trials: survival, attachment, and growth of two Indo-Pacific sponges, \u003cem\u003eNeopetrosia\u003c/em\u003e sp. and \u003cem\u003eStylissa massa\u003c/em\u003e. Journal of Marine Sciences 2012:e417360. \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1155/2012/417360\u003c/span\u003e\u003cspan address=\"10.1155/2012/417360\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eSchmitt S, Hentschel U, Taylor MW (2012a) Deep sequencing reveals diversity and community structure of complex microbiota in five Mediterranean sponges. Hydrobiologia 687:341\u0026ndash;351. \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1007/s10750-011-0799-9\u003c/span\u003e\u003cspan address=\"10.1007/s10750-011-0799-9\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eSchmitt S, Tsai P, Bell J, et al (2012b) Assessing the complex sponge microbiota: core, variable and species-specific bacterial communities in marine sponges. ISME J 6:564\u0026ndash;576. \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1038/ismej.2011.116\u003c/span\u003e\u003cspan address=\"10.1038/ismej.2011.116\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eSchuster A, Knapp I, Pomponi S, et al (2016) Validation of RAD-Seq for sponge phylogenetics\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eSipkema D, Franssen MCR, Osinga R, et al (2005a) Marine Sponges as Pharmacy. Mar Biotechnol 7:142\u0026ndash;162. \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1007/s10126-004-0405-5\u003c/span\u003e\u003cspan address=\"10.1007/s10126-004-0405-5\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eSipkema D, Osinga R, Schatton W, et al (2005b) Large-scale production of pharmaceuticals by marine sponges: Sea, cell, or synthesis? Biotechnology and Bioengineering 90:201\u0026ndash;222. \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1002/bit.20404\u003c/span\u003e\u003cspan address=\"10.1002/bit.20404\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eSlaby BM, Franke A, Rix L, et al (2019) Marine sponge holobionts in health and disease. In: Li Z (ed) Symbiotic Microbiomes of Coral Reefs Sponges and Corals. Springer Netherlands, Dordrecht, pp 81\u0026ndash;104\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eSmith FGW (1941) Sponge disease in British Honduras, and its transmission by water currents. Ecology 22:415\u0026ndash;421. \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.2307/1930719\u003c/span\u003e\u003cspan address=\"10.2307/1930719\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eSteffen K, Indraningrat AAG, Erngren I, et al (2022) Oceanographic setting influences the prokaryotic community and metabolome in deep-sea sponges. Sci Rep 12:3356. \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1038/s41598-022-07292-3\u003c/span\u003e\u003cspan address=\"10.1038/s41598-022-07292-3\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eStrand R, Whalan S, Webster NS, et al (2017) The response of a boreal deep-sea sponge holobiont to acute thermal stress. Sci Rep 7:1660. \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1038/s41598-017-01091-x\u003c/span\u003e\u003cspan address=\"10.1038/s41598-017-01091-x\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eSwierts T, Cleary DFR, de Voogd NJ (2018) Prokaryotic communities of Indo-Pacific giant barrel sponges are more strongly influenced by geography than host phylogeny. FEMS Microbiology Ecology 94:. \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1093/femsec/fiy194\u003c/span\u003e\u003cspan address=\"10.1093/femsec/fiy194\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eTaylor MW, Hill RT, Piel J, et al (2007) Soaking it up: the complex lives of marine sponges and their microbial associates. ISME J 1:187\u0026ndash;190. \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1038/ismej.2007.32\u003c/span\u003e\u003cspan address=\"10.1038/ismej.2007.32\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eTaylor MW, Schupp PJ, Dahll\u0026ouml;f I, et al (2004) Host specificity in marine sponge-associated bacteria, and potential implications for marine microbial diversity. Environ Microbiol 6:121\u0026ndash;130. \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1046/j.1462-2920.2003.00545.x\u003c/span\u003e\u003cspan address=\"10.1046/j.1462-2920.2003.00545.x\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eTernon E, Perino E, Manconi R, et al (2017) How environmental factors affect the production of Guanidine Alkaloids by the Mediterranean sponge \u003cem\u003eCrambe crambe\u003c/em\u003e. Marine Drugs 15:181. \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.3390/md15060181\u003c/span\u003e\u003cspan address=\"10.3390/md15060181\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eThomas T, Moitinho-Silva L, Lurgi M, et al (2016) Diversity, structure and convergent evolution of the global sponge microbiome. Nat Commun 7:11870. \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1038/ncomms11870\u003c/span\u003e\u003cspan address=\"10.1038/ncomms11870\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eThoms C, Horn M, Wagner M, et al (2003) Monitoring microbial diversity and natural product profiles of the sponge Aplysina cavernicola following transplantation. Marine Biology 142:685\u0026ndash;692. \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1007/s00227-002-1000-9\u003c/span\u003e\u003cspan address=\"10.1007/s00227-002-1000-9\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eThoms C, Schupp PJ, Cust\u0026oacute;dio MR, et al (2007) Chemical defense strategies in sponges: a review. Porifera research: biodiversity, Innovation and sustainability 28:627\u0026ndash;637\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eTurnbaugh PJ, Ley RE, Hamady M, et al (2007) The Human Microbiome Project. Nature 449:804\u0026ndash;810. \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1038/nature06244\u003c/span\u003e\u003cspan address=\"10.1038/nature06244\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003evan Treeck P, Eisinger M, M\u0026uuml;ller J, et al (2003) Mariculture trials with Mediterranean sponge species: The exploitation of an old natural resource with sustainable and novel methods. Aquaculture 218:439\u0026ndash;455. \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1016/S0044-8486(03)00010-3\u003c/span\u003e\u003cspan address=\"10.1016/S0044-8486(03)00010-3\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eVavrek M (2011) Fossil: Palaeoecological and Palaeogeographical Analysis Tools. Palaeontologia Electronica 14:16\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eVoultsiadou E (2007) Sponges: an historical survey of their knowledge in Greek antiquity. Journal of the Marine Biological Association of the United Kingdom 87:1757\u0026ndash;1763. \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1017/S0025315407057773\u003c/span\u003e\u003cspan address=\"10.1017/S0025315407057773\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eWebster NS, Cobb RE, Negri AP (2008) Temperature thresholds for bacterial symbiosis with a sponge. ISME J 2:830\u0026ndash;842. \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1038/ismej.2008.42\u003c/span\u003e\u003cspan address=\"10.1038/ismej.2008.42\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eWeigel BL, Erwin PM (2017) Effects of reciprocal transplantation on the microbiome and putative nitrogen cycling functions of the intertidal sponge, \u003cem\u003eHymeniacidon heliophila\u003c/em\u003e. Sci Rep 7:43247. \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1038/srep43247\u003c/span\u003e\u003cspan address=\"10.1038/srep43247\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eWilson MC, Mori T, R\u0026uuml;ckert C, et al (2014) An environmental bacterial taxon with a large and distinct metabolic repertoire. Nature 506:58\u0026ndash;62. \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1038/nature12959\u003c/span\u003e\u003cspan address=\"10.1038/nature12959\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eWorld Meteorological Association (2022) State of the Global Climate 2021\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eWulff J (2001) Assessing and monitoring coral reef sponges: Why and how? Bulletin of Marine Science 69:831\u0026ndash;846\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eYang Q, Franco CMM, Lin H-W, Zhang W (2019) Untapped sponge microbiomes: structure specificity at host order and family levels. FEMS Microbiology Ecology 95:fiz136. \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1093/femsec/fiz136\u003c/span\u003e\u003cspan address=\"10.1093/femsec/fiz136\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eYasuhara-Bell J, Lu Y (2010) Marine compounds and their antiviral activities. Antiviral Research 86:231\u0026ndash;240. \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1016/j.antiviral.2010.03.009\u003c/span\u003e\u003cspan address=\"10.1016/j.antiviral.2010.03.009\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eZaneveld JR, McMinds R, Vega Thurber R (2017) Stress and stability: applying the Anna Karenina principle to animal microbiomes. Nat Microbiol 2:1\u0026ndash;8. \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1038/nmicrobiol.2017.121\u003c/span\u003e\u003cspan address=\"10.1038/nmicrobiol.2017.121\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e\u003c/ol\u003e"}],"fulltextSource":"","fullText":"","funders":[],"hasAdminPriorityOnWorkflow":false,"hasManuscriptDocX":true,"hasOptedInToPreprint":true,"hasPassedJournalQc":"","hasAnyPriority":false,"hideJournal":false,"highlight":"","institution":"","isAcceptedByJournal":true,"isAuthorSuppliedPdf":false,"isDeskRejected":"","isHiddenFromSearch":false,"isInQc":false,"isInWorkflow":false,"isPdf":false,"isPdfUpToDate":true,"isWithdrawnOrRetracted":false,"journal":{"display":true,"email":"[email protected]","identity":"antonie-van-leeuwenhoek","isNatureJournal":false,"hasQc":true,"allowDirectSubmit":false,"externalIdentity":"anto","sideBox":"Learn more about [Antonie van Leeuwenhoek](https://www.springer.com/journal/10482)","snPcode":"10482","submissionUrl":"https://submission.nature.com/new-submission/10482/3","title":"Antonie van Leeuwenhoek","twitterHandle":"","acdcEnabled":true,"dfaEnabled":true,"editorialSystem":"em","reportingPortfolio":"Springer Hybrid","inReviewEnabled":true,"inReviewRevisionsEnabled":false},"keywords":"holobiont, marine sponges, microbiome, farming trials, biogeography, French Polynesia","lastPublishedDoi":"10.21203/rs.3.rs-4002744/v1","lastPublishedDoiUrl":"https://doi.org/10.21203/rs.3.rs-4002744/v1","license":{"name":"CC BY 4.0","url":"https://creativecommons.org/licenses/by/4.0/"},"manuscriptAbstract":"\u003cp\u003eDynamics of microbiomes through time are fundamental regarding survival and resilience of their hosts when facing environmental alterations. As for marine species with commercial applications, such as marine sponges, assessing the temporal change of prokaryotic communities allows us to better consider the adaptation of sponges to aquaculture designs. The present study aims to investigate the factors shaping the microbiome of the sponge \u003cem\u003eDactylospongia metachromia\u003c/em\u003e (Order Dictyoceratida, family Thorectidae), in a context of aquaculture development in French Polynesia (Rangiroa, Tuamotu archipelago). A temporal approach targeting explants collected during farming trials revealed a relative high stability of the prokaryotic diversity, meanwhile a complementary biogeographical study confirmed a spatial specificity amongst samples at different longitudinal scales. Results from this additional spatial analysis confirmed that differences in prokaryotic communities might first be explained by environmental changes (mainly temperature and salinity), while no significant effect of the host phylogeny was observed. The core community of \u003cem\u003eD. metachromia\u003c/em\u003e is thus characterized by a high spatiotemporal constancy, which is a good prospect for the sustainable exploitation of this species towards drug development. Indeed, a microbiome stability across locations and throughout the farming process, as evidenced by our results, should go against a negative influence of sponge translocation during \u003cem\u003ein situ\u003c/em\u003e aquaculture.\u003c/p\u003e","manuscriptTitle":"Prokaryotic communities of the French Polynesian sponge Dactylospongia metachromia display a site-specific and stable diversity during an aquaculture trial","msid":"","msnumber":"","nonDraftVersions":[{"code":1,"date":"2024-03-08 19:05:11","doi":"10.21203/rs.3.rs-4002744/v1","editorialEvents":[{"type":"communityComments","content":0},{"type":"decision","content":"Revision requested","date":"2024-03-25T15:08:30+00:00","index":"","fulltext":""},{"type":"editorInvitedReview","content":"","date":"2024-03-24T07:37:50+00:00","index":"hide","fulltext":""},{"type":"reviewerAgreed","content":"697a0a5a-2006-4d32-98d4-d34bd28558d9","date":"2024-03-04T15:01:59+00:00","index":"hide","fulltext":""},{"type":"reviewersInvited","content":"","date":"2024-03-04T07:46:36+00:00","index":"","fulltext":""},{"type":"checksComplete","content":"","date":"2024-03-02T07:18:37+00:00","index":"","fulltext":""},{"type":"editorAssigned","content":"","date":"2024-03-02T07:18:37+00:00","index":"","fulltext":""},{"type":"submitted","content":"Antonie van Leeuwenhoek","date":"2024-03-01T09:20:19+00:00","index":"","fulltext":""}],"status":"published","journal":{"display":true,"email":"[email protected]","identity":"antonie-van-leeuwenhoek","isNatureJournal":false,"hasQc":true,"allowDirectSubmit":false,"externalIdentity":"anto","sideBox":"Learn more about [Antonie van Leeuwenhoek](https://www.springer.com/journal/10482)","snPcode":"10482","submissionUrl":"https://submission.nature.com/new-submission/10482/3","title":"Antonie van Leeuwenhoek","twitterHandle":"","acdcEnabled":true,"dfaEnabled":true,"editorialSystem":"em","reportingPortfolio":"Springer Hybrid","inReviewEnabled":true,"inReviewRevisionsEnabled":false}}],"origin":"","ownerIdentity":"82878eba-8d47-4c18-ba2f-fc881d17cf50","owner":[],"postedDate":"March 8th, 2024","published":true,"recentEditorialEvents":[],"rejectedJournal":[],"revision":"","amendment":"","status":"published-in-journal","subjectAreas":[],"tags":[],"updatedAt":"2024-04-15T15:12:23+00:00","versionOfRecord":{"articleIdentity":"rs-4002744","link":"https://doi.org/10.1007/s10482-024-01962-0","journal":{"identity":"antonie-van-leeuwenhoek","isVorOnly":false,"title":"Antonie van Leeuwenhoek"},"publishedOn":"2024-04-11 15:01:50","publishedOnDateReadable":"April 11th, 2024"},"versionCreatedAt":"2024-03-08 19:05:11","video":"","vorDoi":"10.1007/s10482-024-01962-0","vorDoiUrl":"https://doi.org/10.1007/s10482-024-01962-0","workflowStages":[]},"version":"v1","identity":"rs-4002744","journalConfig":"researchsquare"},"__N_SSP":true},"page":"/article/[identity]/[[...version]]","query":{"redirect":"/article/rs-4002744","identity":"rs-4002744","version":["v1"]},"buildId":"qtupq5eGEP_6zYnWcrvyt","isFallback":false,"isExperimentalCompile":false,"dynamicIds":[84888],"gssp":true,"scriptLoader":[]}