- Research
- Open access
- Published:
DNA from non-viable bacteria biases diversity estimates in the corals Acropora loripes and Pocillopora acuta
Environmental Microbiome volume 18, Article number: 86 (2023)
Abstract
Background
Nucleic acid-based analytical methods have greatly expanded our understanding of global prokaryotic diversity, yet standard metabarcoding methods provide no information on the most fundamental physiological state of bacteria, viability. Scleractinian corals harbour a complex microbiome in which bacterial symbionts play critical roles in maintaining health and functioning of the holobiont. However, the coral holobiont contains both dead and living bacteria. The former can be the result of corals feeding on bacteria, rapid swings from hyper- to hypoxic conditions in the coral tissue, the presence of antimicrobial compounds in coral mucus, and an abundance of lytic bacteriophages.
Results
By combining propidium monoazide (PMA) treatment with high-throughput sequencing on six coral species (Acropora loripes, A. millepora, A. kenti, Platygyra daedalea, Pocillopora acuta, and Porites lutea) we were able to obtain information on bacterial communities with little noise from non-viable microbial DNA. Metabarcoding of the 16S rRNA gene showed significantly higher community evenness (85%) and species diversity (31%) in untreated compared with PMA-treated tissue for A. loripes only. While PMA-treated coral did not differ significantly from untreated samples in terms of observed number of ASVs, > 30% of ASVs were identified in untreated samples only, suggesting that they originated from cell-free/non-viable DNA. Further, the bacterial community structure was significantly different between PMA-treated and untreated samples for A. loripes and P. acuta indicating that DNA from non-viable microbes can bias community composition data in coral species with low bacterial diversity.
Conclusions
Our study is highly relevant to microbiome studies on coral and other host organisms as it delivers a solution to excluding non-viable DNA in a complex community. These results provide novel insights into the dynamic nature of host-associated microbiomes and underline the importance of applying versatile tools in the analysis of metabarcoding or next-generation sequencing data sets.
Background
Nucleic acid-based analytical methods, in particular 16S ribosomal RNA (rRNA) gene metabarcoding, have greatly expanded our understanding of global prokaryotic diversity. However, standard metabarcoding methods provide no information on the metabolic activity or physiological state of bacteria. Even the most fundamental physiological state, viability, cannot be assessed with conventional DNA-targeted methods such as PCR [1, 2].
One system with a complex microbiome is the coral holobiont [3], which includes the coral animal along with its diverse set of microbial symbionts [4, 5]. Of particular interest is the coral-associated bacterial community, which has been suggested to aid in the defence against pathogens [6, 7], nutrient cycling [8,9,10,11,12,13], the biosynthesis of essential amino acids [14, 15], and protection against climate change [16]. However, corals feeding on bacteria [17], rapid swings from hyper- to hypoxic conditions in the coral tissue [18], the presence of antimicrobial compounds in coral mucus [6], and an abundance of lytic bacteriophages [19,20,21,22,23] can contribute to a considerable prevalence of non-viable microbes and extracellular DNA (exDNA). exDNA and non-viable cells can persist in environmental samples for hours to months [24, 25]. This can bias or overwhelm analyses seeking to understand the diversity of coral-associated viable bacteria in the holobiont. While there are 304 studies using 16S rRNA gene metabarcoding to evaluate bacterial communities in coral samples (from Web of Science Core Collection search across all fields on 20 Nov. 2023 for: coral AND (metagenom* OR metabarcoding OR amplicon) AND bacteria), none have assessed bacterial community viability.
One method used for community viability characterization is high-throughput sequencing combined with a propidium monoazide (PMA) treatment to deplete signals from exDNA [26, 27]. PMA is not permeable to intact cell membranes; thus, the dye only interacts with DNA in membrane-compromised cells (i.e., dead or damaged cells) or exDNA. Once the PMA is inside the cells it intercalates into the DNA and, after photoactivation, is crosslinked to the DNA, preventing PCR amplification [28]. PMA has been used in combination with PCR and 16S rRNA gene metabarcoding to selectively detect viable microbes without signals from damaged cells and exDNA in soil [26], sewage sludge [29], and guts of mammals [30,31,32,33] or fish [34, 35]. In these systems, DNA from non-viable microbes represented 9–40% of the amplified community and inflated community diversity by 25–55%. Overall, exDNA or DNA from non-viable microbes has been shown to obscure experimental treatment effects, spatiotemporal patterns, and relationships between microbial taxa and environmental conditions.
Here, we characterize the viable bacterial communities from six Great Barrier Reef (GBR)-sourced coral species, Acropora loripes, A. millepora, A. kenti, Platygyra daedalea, Pocillopora acuta, and Porites lutea.
Methods
Coral sampling
All corals used in this study were part of the coral stock at the Australian Institute of Marine Science's (AIMS) National Sea Simulator (SeaSim) in Townsville, Queensland, where they were housed in outdoor mesocosms with natural light conditions and a flow-through system with a full seawater exchange every two hours. Details about collection dates and locations are in Additional file 1: Table S1. The daily temperature profile followed the average temperature profile recorded at Davies Reef weather station. Corals were maintained on a daily regimen of 0.5 Artemia nauplii ml−1 and 2000 cells ml−1 of a mixed-species microalgae solution. In June 2022, ~ 2 cm fragments were cut from five different positions in each of five colonies (distinct genotypes) of A. loripes, P. daedalea, P. acuta, and P. lutea without removing the colonies from the aquarium (n = 25 each). Four genotypes of A. kenti and A. millepora were available in the SeaSim at the time of sampling (Additional file 1: Table S1), thus 20 instead of 25 samples were collected for each of these two species. GBR acroporids have recently undergone taxonomic revision; phylogenomic reconstructions consistently places ‘A. tenuis’ sampled from the GBR with other specimens from eastern Australia, identified as A. kenti [36]. Following this formal taxonomic revision, we have adopted this nomenclature. To minimise cross-contamination, nitrile gloves were discarded between each sampling, and all collection equipment (bone cutters and forceps) was sequentially sterilized in 10% sodium hypochlorite, reverse osmosis water, 80% ethanol, with a final wash in 0.22 µm filtered seawater (fSW) as described in Damjanovic et al. [37].
Fragments were gently rinsed in fSW to remove loosely associated bacteria. Using an airbrush, tissue was removed from the fragment into a volume of 10 ml fSW, which was transferred to a 15 ml polypropylene tube. During the tissue airbrushing, 12 blanks were collected to account for contamination that might be introduced at this stage. Coral homogenate and blank samples were centrifuged (5250 × g for 10 min at room temperature (RT)), supernatant removed, and pelleted cells resuspended in 1000 µl 0.22 µm filter sterilized 3× concentrated phosphate buffered saline. Aliquots of 400 µl were transferred to two separate 1.5 ml tubes; the “A” tube to receive PMA treatment and the “B” tube as untreated. The untreated “B” tube was centrifuged (5000×g for 10 min at RT), supernatant removed, and pelleted cells stored at − 20 °C until shipment.
PMA treatment
Extracellular DNA (exDNA) and DNA from membrane-compromised bacterial cells was removed by treatment of coral homogenates and tissue airbrushed blanks with the viability dye PMAxx™ (Biotium, Fremont, CA, USA). PMAxx is a DNA-intercalating agent that forms photo-induced crosslinks making the bound DNA inaccessible for downstream molecular applications. PMAxx was added to the “A” fraction at a final concentration of 25 µM by adding 1 µl of 10 mM PMAxx stock to the 400 µl aliquot. This was followed by 10 min incubation in a rotary mixer inside a dark room at RT. Photoactivation was performed by using the PMA-Lite™ LED Photolysis Device (Biotium) with the exposure time set to 30 min. During this exposure, samples were mixed every 3 min to ensure adequate photoactivation. After photoactivation, the PMA-treated “A” tube was centrifuged (5000 × g for 10 min at RT), supernatant removed, and pelleted cells stored at − 20 °C until shipment.
The efficacy of PMAxx in cnidarian tissue was confirmed prior to coral collection at AIMS using a model for corals, the sea anemone Exaiptasia diaphana [38]. Triplicate heat-killed anemone samples (250 µl homogenized tissue heated to 90 °C for 60 min and plated to confirm no bacterial growth) were spiked with either 3.7 × 106 heat-killed or viable Endozoicomonas sp. cells (strain ALC013, provided by Cecilie Goetze). Additional triplicate heat-killed anemone samples were included without addition of bacteria. Subsequently, these samples (n = 9) were split into two aliquots of which one was treated with PMAxx as described above and one left untreated. DNA was extracted from these samples (as below) and bacterial load quantified using droplet digital PCR (ddPCR) targeting the bacterial 16S rRNA genes.
DNA extraction and library prep
PMA-treated and untreated samples (n = 302) were shipped to The University of Melbourne for DNA extractions, which were completed as described by Dungan et al. [39]. Blank DNA extractions (n = 9) were included to identify reagent contamination.
Extracted DNA was amplified using bacterial primers targeting the V5–V6 regions of the 16S rRNA gene: 784F [5′ GTGACCTATGAACTCAGGAGTCAGGATTAGATACCCTGGTA 3′], 1061R [5′ CTGAGACTTGCACATCGCAGC-CRRCACGAGCTGACGAC 3′] [40] with overhang adapters (underlined). Negative template PCR controls (n = 7) were included to test for potential contamination. Triplicate PCRs were carried out in 15 µl reactions containing 1 × UCP Multiplex PCR Master Mix (Qiagen, Venlo, The Netherlands), 0.3 μmol l−1 each of the forward and reverse primers, and 1 µl of template DNA (1:40 dilution) or nuclease-free water (negative template control). The PCR cycling conditions consisted of an initial heating step at 95 °C for 3 min; 18 cycles of 95 °C for 15 s, 55 °C for 30 s and 72 °C for 30 s; and a final extension step of 72 °C for 7 min. Triplicate PCR products were then pooled; successful DNA extraction was confirmed by agarose gel electrophoresis.
A volume of 20 µl of each PCR product pool was purified by size-selection using Nucleomag NGS Clean-up and Size Select beads (Scientifix, Clayton, VIC, Australia). The purified DNA was resuspended in 40 µl of nuclease-free water. Indexing PCRs were created by combining 10 μl of purified DNA with 10 μl 2 × MyTaq HS Mix polymerase (Bioline, Narellan, NSW, Australia) and 1 μl (5 μM) of forward and reverse indexing primers. The PCR cycling conditions consisted of an initial heating step at 95 °C for 3 min; 24 cycles of 95 °C for 15 s, 60 °C for 30 s and 72 °C for 30 s; and a final extension step of 72 °C for 7 min. For a subset of randomly chosen samples, product size was confirmed by agarose gel electrophoresis. Barcoded PCR amplicons were pooled by species; pools were purified with a bead clean up. Each pool was checked for quality and quantity (2200 TapeStation; Agilent Technologies, Mulgrave, VIC, Australia) to guide pool normalization at equimolar concentrations, then sequenced on a single Illumina MiSeq run using v3 (2 × 300 bp) reagents at the Walter and Eliza Hall Institute, Melbourne, Australia.
Metabarcoding
Raw demultiplexed 16S rRNA gene sequences were imported into QIIME2 v2022.11 [41] where primers and overhang adapters were removed [using cutadapt v2.6; 42]. The data were then filtered, denoised using the pseudo-pooling method, and chimera checked [using DADA2 v1.22.0; 43] to generate amplicon sequence variants (ASVs). Taxonomy for each ASV was assigned against a SILVA database (version 138) trained with a naïve Bayes classifier against the same V5-V6 region targeted for sequencing [44]. A phylogenetic tree was produced in QIIME2 by aligning ASVs using the PyNAST method [45] with mid-point rooting.
Droplet digital PCR
ddPCR was used to track changes in the absolute abundance of bacteria and Endozoicomonas spp. between PMA-treated and untreated samples for A. loripes only. ddPCRs for each sample (n = 50) were prepared in an initial volume of 72 µl comprising 37.5 µl EvaGreen Supermix (QX200; Bio-Rad, South Granville, NSW, Australia), 31.5 µl sterile water, and 3 µl diluted DNA template. The solution was mixed and then split into three 24 µl aliquots, one for host cell quantification (A. loripes β-actin gene), one for all bacteria, and a third to target Endozoicomonas spp. The β-actin gene occurs as a single copy in A. loripes and is used here as an internal control gene. Primers to target the A. loripes β-actin gene were designed using Primer3 [v4.1.0; 46] using the β-actin gene sequence from the existing draft genome [47]. One microlitre each of the appropriate 5 µM forward and reverse primers (Table 1) were added to each reaction aliquot, giving final primer concentrations of 200 nM and volumes of 25 µl. From each 25 µl volume, 20 µl was loaded into a DG8 cartridge (1864008; Bio-Rad), followed by 70 µl of droplet generation oil for EvaGreen (1864005; Bio-Rad), and droplets were generated in a droplet- generator (QX200; Bio-Rad). A volume of 40 µl of generated droplets per reaction was then transferred to a 96-well plate and foil-sealed (1814040; Bio-Rad) with a thermal plate sealer (PX1; Bio-Rad). No-template controls and duplicate samples were included in each plate. Thermal cycler settings were: one cycle at 95.0 °C for 6 min; 40 cycles at 95 °C for 15 s + 56 °C for 15 s + 72 °C for 15 s; one cycle at 72 °C for 4 min; 12 °C hold following Hartman et al. [48].
Droplets were read on a Bio-Rad QX200 droplet reader, and fluorescence data were analysed in QuantaSoft v1.7.4.0917. To quantify the change in Endozoicomonas spp. abundance in coral cells, the 16S rRNA gene of Endozoicomonas was measured using Endozoicomonas-specific primers and normalized using the coral β-actin gene. To quantify the variation in total bacterial abundance in coral, the 16S rRNA gene in bacteria was quantified using bacterial primers and normalized with the β-actin gene. To determine the dynamics in relative abundance of Endozoicomonas to the total bacterial community, the Endozoicomonas 16S rRNA gene was measured using Endozoicomonas-specific primers and normalized by the 16S rRNA gene of total bacteria.
Statistical analysis
All data were analyzed in R [v4.2.3; 49]. For metabarcoding data, ASV, taxonomy, metadata and phylogenetic tree files were imported into R and combined into a phyloseq object [50]. Contaminant ASVs were identified and removed sequentially from the dataset according to their abundance in the PCR negative controls, extraction blanks, and tissue processing blanks relative to the samples using the prevalence method in the R package decontam with p = 0.1 [51].
α-diversity metrics (observed ASVs, Shannon’s and inverse Simpson’s indices) were calculated from a rarefied dataset by species using the ‘estimate_richness’ function in the R package phyloseq [50]. These and the ddPCR data were then analyzed using linear mixed effects models, with PMA treatment and coral genotype as fixed effects, and aquarium as a random effect, using the R package nlme [52]. Where appropriate, post hoc comparisons were performed using Tukey’s honestly significant difference test in the R package emmeans [53]. Box plots were created with ggplot2 [54], plotrix [55], and gridExtra [56] by merging data by PMA treatment. Summary data was produced using the function ‘summarySE’ from the package Rmisc [57].
β-diversity was evaluated using a weighted Unifrac distance matrix. The difference in community compositions among groups (PMA treatment and coral genotype) was calculated using ‘adonis2’ (a modified version of a permutational multivariate analysis of variance (PERMANOVA)) in the vegan package in R [58] with ‘Tank’ assigned as a random effect and 9999 permutations. Where the assumption of homogeneity of dispersion was not met, a permutation test for multivariate dispersion to account for unequal dispersion between groups was included. Where appropriate, Holm corrected pairwise comparisons were computed using the package pairwiseAdonis [59].This method tests the hypothesis that the bacterial communities from different groups (e.g., PMA treatment and genotype) are significantly different from each other in terms of their composition. PERMANOVA uses a distance matrix as input and is often more robust to non-normality and heteroscedasticity than traditional ANOVA. The results of the PERMANOVA were cross validated with and visualized using principal coordinates analysis (PCoA) to gain a more robust understanding of the underlying patterns. Where there were significant differences in community composition between PMA treatments for a given species, ASVs that were significantly associated with treatment groups were identified using the indicspecies function ‘multipatt’ [60] with specificity (At) and fidelity (Bt) set to 0.6. Barplots visualizing overall microbiota composition were made using ggplot2 [54] by agglomerating taxa at the genus or ASV level based on relative abundances. Venn diagrams were created using packages Eulerr [61] and microbiome [62] to identify ASVs in PMA-treated or untreated samples for all corals and for each species.
For all tests p-values < 0.05 were considered significant.
Results
Method validation
Prior to collecting coral samples, we examined the efficacy of PMAxx treatment in removing DNA from non-viable bacteria in a complex cnidarian matrix. We used heat-killed anemone tissue as our model system to mimic the complexity and opacity of homogenized coral tissue, while removing the potential influence of a living bacterial community. Endozoicomonas sp. was used to supplement the tissue and provide signal as it is one of the most ubiquitous bacterial symbionts of corals. Spiking of heat-killed anemone tissue samples with heat-killed or live Endozoicomonas showed that PMAxx treatment effectively eliminated DNA from non-viable bacteria; there was a substantial reduction in copies of the 16S rRNA gene (per µl) in ddPCR (using general bacteria probes) when PMA was applied to anemone tissue spiked with heat-killed bacteria (mean ± SE, 68.2 ± 26.3) compared to the untreated fraction (2358.7 ± 225.8), equivalent to a signal reduction of 97.1%. This copy number was largely recovered when anemone tissue was spiked with live Endozoicomonas in parallel with PMA treatment (Additional file 2: Fig. S1), though signal after bacterial addition with PMA treatment was only half that of the untreated fraction suggesting some bacteria in the added culture may have not been viable.
Metabarcoding data
Sequencing produced 10.7 M reads across the six coral species (n = 290), sampling blanks (n = 12), extraction blanks (n = 9), and PCR negative control samples (n = 7). After merging, denoising and chimera filtering, 6.8 M reads remained. ASVs with fewer than 10 reads were removed from the dataset as they likely represent sequencing errors. Decontam identified eight putative contaminant ASVs from PCR amplification (0.02% of total reads), eight putative contaminant ASVs from DNA extraction (0.07% of total reads), and 11 putative contaminant ASVs from tissue blasting blanks (0.50% of total reads) (Additional file 1: Table S2). Nine of 290 samples (four A. loripes, one A. millepora, four P. daedalea) were removed from analysis as they had < 500 reads. After all filtering steps, there were 8,499 ASVs across the remaining coral samples (n = 281).
Coral species were subset into separate phyloseq objects for analysis, and independently rarefied (Table 2) based on rarefaction curves for species richness (Additional file 2: Fig. S2). PMA treatment had a significant effect on α-diversity indices for the microbiota present in A. loripes only. The evenness (inverse Simpson’s index) in untreated (mean ± SE) (9.20 ± 1.12) compared with PMA-treated (4.97 ± 0.55) samples showed an 85% inflation as result of non-viable DNA. Likewise, Shannon’s Index of diversity (2.49 ± 0.10) in untreated samples was overestimated by 31% compared to the PMA treatment (1.90 ± 0.12) (Fig. 1; Table 2).
For A. kenti and P. daedalea, there was a significant effect of host genotype on all α-diversity parameters (Table 2; Additional file 2: Fig. S3), with one genotype largely driving the differences (A for A. kenti and B for P. daedalea). Only Shannon’s index of diversity was significantly different by genotype for bacterial communities in A. millepora (Table 2). While PMA-treated corals did not differ significantly from untreated samples in terms of mean ASV richness for any species, of the 8,499 ASVs in the full coral dataset, > 30% were identified in untreated samples only, alluding that they originated from exDNA or non-viable DNA (Fig. 2A).
PERMANOVA analysis revealed a significant effect of PMA treatment on bacterial community structure for A. loripes and P. acuta samples, regardless of genotype (Table 2), which is illustrated by the distribution of samples in the PCoA of weighted unifrac distances and community composition (Fig. 3). The communities of both PMA-treated and untreated A. loripes were dominated by Endozoicomonas spp.; however, the relative abundance of Endozoicomonas shifted from 29.5% in the untreated to 78.2% in the PMA-treated samples (Fig. 3A). Untreated A. loripes had 24-fold greater relative abundance of an unknown Rickettsiales (26.4% untreated, 1.1% PMA-treated) and 5 × more Candidatus Fritschea, a member of the Simkaniaceae family (3.7% untreated, 0.7% PMA-treated). PMA-treated samples had sevenfold more Vibrio (0.4% untreated, 3.0% PMA-treated) and saw the emergence of Acinetobacter (0.002% untreated, 4.0% PMA-treated). For P. acuta (Fig. 3C), PMA treatment resulted in a loss of Reinekea (15.7% untreated, 0% PMA-treated) and Nonlabens (8.2% untreated, 2.0% PMA-treated), with increases in Salinipshaera (1.8% untreated, 8.8% PMA-treated), an unknown Oceanospirillales (2.3% untreated, 5.4% PMA-treated) and the appearance of XBB1006 (0% untreated, 6.1% PMA-treated). The bacterial β-diversity of A. millepora, A. kenti, P. daedalea, and P. lutea was not impacted by PMA treatment, though there were differences in the presence of individual ASVs by PMA treatment (Fig. 2) and all but A. millepora had a significant impact of host genotype on bacterial community composition (Table 2, Additional file 2: Fig. S4. Additional file 1: Table S3).
The indicator species analysis provided two ASVs characteristic of the microbiomes of the PMA-treated or untreated samples for A. loripes and one ASV for P. acuta (Table 3). ASV001 was identified as an indicator for the PMA-treatment in A. loripes; BLASTn output for this sequence provided Endozocomonas coralli and E. atrinae as top hits at 99.22% identity. There was a 2.3-fold increase in ASV001 from the untreated to the PMA-treated condition. This ASV was also found, albeit in low relative abundance, in A. millepora (Additional file 2: Fig. S5A-B). ASV002 is an unknown member of the class Alphaproteobacteria, potentially belonging to the order Rickettsiales or Rhodospirialles. This ASV was an indicator of the untreated condition in A. loripes where it decreased 23-fold when samples were PMA-treated (Table 3). Only one ASV was identified as an indicator taxon for P. acuta, ASV003. This ASV was a 100% match to Reinekea thalattae and, while it made up 15.1% of the abundance of all bacteria in the untreated condition, it virtually disappeared with PMA-treatment (0.002% relative abundance).
A. loripes had 72 Endozoicomonas ASVs between all the samples (Additional file 2: Fig. S5A); ASV001 was the only variant identified as being indicative of PMA treatment. Other coral species were also dominated by Endozoicomonas including A. millepora (70 ASVs), A. kenti (33 ASVs), and P. daedalea (33 ASVs) (Additional file 2: Fig. S5). Four of the P. daedalea genotypes were dominated by one ASV (ASV023; Additional file 2: Fig. S5D), an unknown Rhodanobacteraceae who’s top BLASTn hit (98.07% identity) was Fulvimonas yonginensis.
ddPCR
Given the significant effect of PMA-treatment on α-and β-diversity for bacterial communities associated with A. loripes, particularly Endozoicomonas spp., these samples were further processed to quantify bacterial load. For each sample, three data points were collected using ddPCR: 1) bacterial 16S rRNA gene copies per coral cell; 2) Endozoicomonas spp. 16S rRNA gene copies per coral cell; and 3) the relative abundance of Endozoicomonas. Contrary to expectations, PMA-treated samples had significantly more bacterial gene copies (F(1,36) = 38.74, p < 0.0001) and Endozoicomonas (F(1,37) = 15.74, p = 0.0003) per coral cell, but there was no difference in the ratio of Endozoicomonas to bacterial signal between the two groups (Fig. 4). There was no effect of A. loripes genotype on any of the ddPCR data.
Discussion
Viability testing shows significant differences in bacteria associated with low diversity corals
Metabarcoding is often non- or semi-quantitative. DNA extraction and amplification bias, 16S rRNA gene copy number variation, bioinformatics classification accuracy, and (non)viability all influence the degree to which metabarcoding truly quantifies a community’s taxonomic profile. The use of 16S rRNA/16S rRNA gene ratio to identify viable bacteria in the coral system is gradually appearing [16, 65]. However, this approach shows inconsistent results in other systems [66]. Though not interpreted as such, Sun et al. [67] indirectly looked at the variation in coral-associated bacterial communities based on viability when comparing metatranscriptome and metagenome data, finding significant differences in the microbial community composition at DNA and RNA levels. This study is the first to apply PMA in combination with metabarcoding to coral tissue samples and shows that DNA from non-viable microbes significantly inflates community evenness (85%) and species diversity (31%) in A. loripes and bacterial community composition for A. loripes and P. acuta. This exaggeration in diversity with the presence of non-viable bacteria has previously been reported in a range of systems from rabbits to rainwater (Additional file 1: Table S4), but for the first time here for corals.
PMA significantly altered α- and β- bacterial diversity metrics in A. loripes and β-diversity in P. acuta. Average bacterial richness in adult corals can range from 20 to > 500 ASVs [68, 69]. A. loripes and P. acuta are on the lower end of this spectrum with (mean ± SE) 24 ± 2 and 20 ± 2 ASVs per sample, respectively. The variation in response to PMA treatment by coral species could be explained by differential and specific antibacterial capabilities of the surface mucus layer [70], whereby some species may be more proficient at degrading environmentally sourced bacteria. Alternatively, coral species may have unequal expression of immune pathway genes, such as those involved in the Toll-like and nucleotide-binding oligomerization domain-like receptor signalling pathways, resulting in a range of responses to bacterial exposure. Low diversity corals may have a more active role in selecting their bacterial symbionts and have immune mechanisms in place to influence community composition. Further, since we did not quantify the number of bacteria per coral cell for all species, it is possible that low diversity corals had fewer bacteria all together and that eliminating exDNA and DNA from dead bacteria had a greater influence on the abundance of other viable members.
These results correspond with previous studies that show bacterial communities post-PMA treatment had significantly reduced α-[29, 31, 33, 34] or β-[31, 34, 71] diversity indices, suggesting that standard 16S rRNA gene metabarcoding methods miscalculate coral-associated bacterial diversity for low diversity systems. For A. loripes the PCoA (Fig. 3B) shows greater variability in community composition for untreated compared to PMA-treated corals. This suggests that DNA from non-viable microbes may obscure effects of treatments between bacteria and environmental conditions for this species by exaggerating diversity.
The differences between untreated and PMA-treated bacterial communities were driven by a few taxa
For A. loripes, the reduction in diversity and evenness of bacterial communities was largely reflected in the abundance of Endozoicomonas spp. As non-viable bacteria and exDNA signal were removed, Endozoicomonas increased from 29.5% to 78.2% relative abundance in the untreated to PMA-treated condition. This suggests that in hospite Endozoicomonas are alive with few dead cells. Correspondingly, indicator species analysis identified ASV001, an Endozoicomonas sp., indicative of PMA treatment. In addition to ASV001, there were 71 other Endozoicomonas ASVs in A. loripes (Additional file 2: Fig. S5A). The dominance of Endozoicomonas spp. in the viable bacterial community of A. loripes adds a valuable component to other work showing Endozoicomonas are a key symbiotic partner in this species [37, 72], other acroporids [73,74,75,76,77,78], and a global distribution of corals [11, 79,80,81,82].
Bacterial viability studies have found that previously identified “key” bacterial taxa were not active in the community sampled [34]. Here were saw a depletion of Rickettsiales in the viable bacterial community of A. loripes, a member of which was identified as an indicator taxon for the untreated samples (ASV002). Rickettsiales are consistently found in Acropora spp. [83,84,85,86,87] and corals globally [88]. This taxon is thought to be associated with white band disease [89] with genomic signatures of pathogenicity [88], though Casas et al. [83] found Rickettsiales to be widespread in both healthy and diseased A. cervicornis and A. palmata in the Caribbean. The absence of Rickettsiales in the viable bacterial community of A. loripes in contrast to their abundance in the untreated condition implies that they were predominately dead cells or exDNA. Their abundance could be explained by changes in environmental conditions that allow Rickettsiales to proliferate, coral immune responses, acquisition of Rickettsiales exDNA from the external environment, or feeding. Coral-associated Rickettsiales thrive under high nutrient conditions [90]; Rickettsiales could have bloomed under previous conditions and crashed in low nutrient load tanks. Alternatively, it is possible that acroporids have developed strategies to prevent infection by Rickettsiales leading to an abundance of dead cells. Members of Rickettsiales have been identified as infecting coral hosts through horizontal transmission [91], suggesting that corals may acquire this taxon from the surrounding seawater, possibly though heterotrophic feeding. Corals are known to heterotrophically feed on available bacteria [92,93,94] and Rickettsiales are found in tropical marine waters surrounding coral reefs [95].
ASV003 (Reinekea thalattae) was identified as an indicator taxon for untreated P. acuta, where it made up 15.1% of the microbiome and was depleted with PMA-treatment (0.002% relative abundance). Reinekea spp. have been previously reported in the coral-associated microbiome [96], potentially associated with disease [97], but there is otherwise little information about this taxon and its relationship to the coral holobiont.
Not all coral species harbour a distinct viable microbial community, but there are differences
Four of the six coral species tested did not have significant differences by PMA treatment for bacterial α-or β-diversity. However, there were clear differences in the ASVs found in PMA-treated versus untreated samples, where 24.8% (A. kenti)–36.1% (A. millepora) of the unique reads were found only in untreated samples. Since this value considers all ASVs in all samples (Fig. 2A) or all ASVs for each coral species (Fig. 2B-G), as opposed to the average number of ASV’s per sample (α-diversity richness measure), this suggests that, particularly for the rare/low abundant taxa, the ASVs are different between PMA-treated and untreated fractions. One hypothesis is that specific more highly abundant bacterial sequences are preferentially amplified in the untreated samples. A reduction in the abundance of these sequences following PMA treatment would then allow primers to bind to low abundant templates that were previously unbound. This is in line with observations that PMA treatment increases evenness and diversity without a change in the number of ASVs, a phenomenon driven by increased detection of rare taxa [30, 31, 98]. We also found that that coral microbiota varied by genotype for A. kenti, P. daedalea, and P. lutea, which has previously been reported [78, 87, 99].
Coral tissue-associated bacteria can form dense clusters termed cell-associated microbial aggregates (CAMAs). CAMAs have previously been visualized in all coral genera used in this study—Acropora, Platygyra, Pocillopora, and Porites [79, 80, 100,101,102,103,104,105,106,107,108,109]. One consistent trend is that Endozoicomas, a widespread coral symbiont [110], forms CAMAs. No studies have yet visualized CAMAs in A. kenti, but the consistency and abundance of Endozoicomonas in this study and others suggest that CAMAs can occur in this species. CAMA formation in P. daedalea has not been described since Work TM and Aeby GS [104] but the high abundance (> 75%) of ASV023 in P. daedalea (Additional file 2: Fig. S5D) is consistent with the presence of CAMAs. ASV023 is an unknown Rhodanobacteraceae who’s top BLASTn hit (98.07% identity) was Fulvimonas yonginensis. This family has been found in the mucus of A. millepora [111], but otherwise little literature exists on its relationship with coral.
Eukaryotic DNA influences ddPCR results and can skew normalized bacterial results
ddPCR, qPCR, or fluorescence activated cell sorting (FACS), are elegant approaches for quantitative microbiome profiling, which can complement metabarcoding findings and can be used to track the amount of microbial DNA derived from non-viable bacteria. However, it is important to realize that technical sources of variability may introduce substantial additional biases depending on the quantification method used [33]. As with qPCR, ddPCR uses Taq polymerase in a standard PCR to amplify the target DNA. The ddPCR technology, however, partitions the PCR into thousands of droplets (individual reaction vessels) prior to amplification and acquires the data at the reaction end point. This enables more precise and reproducible data and direct quantification without the need for standard curves [112]. ddPCR results with A. loripes show that the number of bacteria per coral cell and Endozoicomonas per coral cell were significantly greater in the PMA-treated samples, with no difference in relative abundance of Endozoicomonas. This finding is surprising because the metabarcoding data showed an increase in Endozoicomonas relative abundance following PMA treatment. Given that PMA treatment blocks non-viable bacterial DNA from being amplified, the opposite result would be expected. Because PMA also targets eukaryotic DNA, it is probable that there was a substantial number of damaged host cells during the tissue removal via air brushing that were labelled by PMA. This would in turn shift the proportion of bacteria relative to coral cells to be greater in the PMA treatment and result in no difference when looking at the ratio of Endozoicomonas to all bacteria. Possible ways to mitigate this issue are to use normalization approaches that are not dependent on DNA, such as surface area of the collected fragment or by the number of coral polyps. While it was not the motivation of this study, one advantage of the host DNA depletion by PMA is the application of shotgun metagenomics. Eukaryotic DNA can dominate coral samples, obscuring changes in microbial populations because few DNA sequence reads are from the microbial component [113, 114]. PMA treatment after osmotic lysis of human saliva samples has been a successful approach to remove host-derived metagenomic sequences [115] and could be an option to deplete coral DNA from coral tissue as well. In addition to PCR, successful PMA treatment has been reported in combination with loop-mediated isothermal amplification [116], multiple displacement amplification [117], flow cytometry [33, 118], and metagenomics [115, 118, 119]. There are other promising applications such as PMA in combination with fluorescence in situ hybridization (PMA-FISH) [120] that can be explored in future work.
Despite the many achievements of PMA-based methods in microbial studies (Additional file 1: Table S4), as this is the first time PMA has been used in a coral study, there are limitations that future studies can address. These limitations include creating optimal PMA application conditions to decrease the likelihood of false positives (dead cells not being bound by PMA) or false negatives (living cells being labelled with PMA). Factors such as light penetration [121], which would need to be adjusted based on Symbiodiniaceae density, PMA concentration, and PMA responsiveness to gram positive versus gram negative bacteria [122], would potentially need to be tailored for each coral species tested. Further, PMA is limited by its requirement of a damaged membrane for penetration. Not all dead bacteria have a damaged membrane, and in the case of this study, sample processing can damage membranes of living eukaryotic cells.
Conclusions
Coral microbiomes are dynamic ecosystems, composed of tens to hundreds of unique microbial taxa. Metabarcoding of 16S rRNA genes has been widely adopted to understand this diversity more generally, making the profiling of existing microbiomes in different host species a common analysis in the coral field [123]. The function of these microbiomes depends on which members of the community are alive or dead. While only the bacterial communities of A. loripes and P. acuta were significantly changed by PMA, any removal of non-viable signal has the potential to increase knowledge of the potentially metabolically active microbes. These results provide novel insights into the dynamic nature of host-associated microbiomes and underline the importance of applying versatile tools in the analysis of metabarcoding or next-generation sequencing data sets.
Availability of data and materials
The datasets generated and analysed during the current study are available under NCBI BioProject ID PRJNA971764. QIIME2 and R code can be found at https://github.com/adungan31/bacteria_viability.
References
Emerson JB, Adams RI, Roman CMB, Brooks B, Coil DA, Dahlhausen K, Ganz HH, Hartmann EM, Hsu T, Justice NB, et al. Schrodinger’s microbes: Tools for distinguishing the living from the dead in microbial ecosystems. Microbiome. 2017;5:86.
Cangelosi GA, Meschke JS. Dead or alive: molecular assessment of microbial viability. Appl Environ Microbiol. 2014;80:5884–91.
Rohwer F, Seguritan V, Azam F, Knowlton N. Diversity and distribution of coral-associated bacteria. Mar Ecol Prog Ser. 2002;243:10.
Ricci F, Rossetto Marcelino V, Blackall LL, Kühl M, Medina M, Verbruggen H. Beneath the surface: community assembly and functions of the coral skeleton microbiome. Microbiome. 2019;7:159.
Blackall LL, Wilson B, van Oppen MJ. Coral-the world’s most diverse symbiotic ecosystem. Mol Ecol. 2015;24:5330–47.
Ritchie KB. Regulation of microbial populations by coral surface mucus and mucus-associated bacteria. Mar Ecol Prog Ser. 2006;322:1–14.
Morrow KM, Muller E, Lesser MP. How Does the Coral Microbiome Cause, Respond to, or Modulate the Bleaching Process? In: van Oppen MJH, Lough JM, editors. Coral Bleaching: Patterns, Processes, Causes and Consequences. Cham: Springer International Publishing; 2018. p. 153–88.
Ceh J, Kilburn MR, Cliff JB, Raina J-B, van Keulen M, Bourne DG. Nutrient cycling in early coral life stages: Pocillopora damicornis larvae provide their algal symbiont (Symbiodinium) with nitrogen acquired from bacterial associates. Ecol Evol. 2013;3:2393–400.
Lema KA, Clode PL, Kilburn MR, Thornton R, Willis BL, Bourne DG. Imaging the uptake of nitrogen-fixing bacteria into larvae of the coral Acropora millepora. ISME J. 2016;10:1804–8.
Kimes NE, Van Nostrand JD, Weil E, Zhou J, Morris PJ. Microbial functional structure of Montastraea faveolata, an important Caribbean reef-building coral, differs between healthy and yellow-band diseased colonies. Environ Microbiol. 2010;12:541–56.
Bourne DG, Morrow KM, Webster NS. Insights into the coral microbiome: underpinning the health and resilience of reef ecosystems. Annu Rev Microbiol. 2016;70:317–40.
Raina JB, Tapiolas DM, Willis BL, Bourne DG. Coral-associated bacteria and their role in the biogeochemical cycling of sulfur. Appl Environ Microbiol. 2009;75:3492–501.
Tandon K, Lu CY, Chiang PW, Wada N, Yang SH, Chan YF, Chen PY, Chang HY, Chiou YJ, Chou MS, et al. Comparative genomics: dominant coral-bacterium Endozoicomonas acroporae metabolizes dimethylsulfoniopropionate (DMSP). ISME J. 2020;14:1290–303.
Matthews JL, Raina J-B, Kahlke T, Seymour JR, van Oppen MJH, Suggett DJ. Symbiodiniaceae-bacteria interactions: rethinking metabolite exchange in reef-building corals as multi-partner metabolic networks. Environ Microbiol. 2020;22:1675–87.
Robbins SJ, Singleton CM, Chan CX, Messer LF, Geers AU, Ying H, Baker A, Bell SC, Morrow KM, Ragan MA, et al. A genomic view of the reef-building coral Porites lutea and its microbial symbionts. Nat Microbiol. 2019;4:2090–100.
Ziegler M, Seneca FO, Yum LK, Palumbi SR, Voolstra CR. Bacterial community dynamics are linked to patterns of coral heat tolerance. Nat Commun. 2017;8:14213.
Meunier V, Bonnet S, Pernice M, Benavides M, Lorrain A, Grosso O, Lambert C, Houlbreque F. Bleaching forces coral’s heterotrophy on diazotrophs and Synechococcus. ISME J. 2019;13:2882–6.
Linsmayer LB, Deheyn DD, Tomanek L, Tresguerres M. Dynamic regulation of coral energy metabolism throughout the diel cycle. Sci Rep. 2020;10:19881.
Sweet M, Bythell J. The role of viruses in coral health and disease. J Invertebr Pathol. 2017;147:136–44.
Wood-Charlson EM, Weynberg KD, Suttle CA, Roux S, van Oppen MJ. Metagenomic characterization of viral communities in corals: mining biological signal from methodological noise. Environ Microbiol. 2015;17:3440–9.
van Oppen MJH, Leong J, Gates RD. Coral-virus interactions: A double-edged sword? Symbiosis. 2009;47:1–8.
Messyasz A, Rosales SM, Mueller RS, Sawyer T, Correa AMS, Thurber AR, Vega Thurber R. Coral bleaching phenotypes associated with differential abundances of nucleocytoplasmic large DNA viruses. Front Mar Sci. 2020;7:555474.
Weynberg KD, Laffy PW, Wood-Charlson EM, Turaev D, Rattei T, Webster NS, van Oppen MJH. Coral-associated viral communities show high levels of diversity and host auxiliary functions. PeerJ. 2017;5:e4054.
Agnelli A, Ascher J, Corti G, Ceccherini MT, Pietramellara G, Nannipieri P. Purification and isotopic signatures (δ 13C, δ 15N, Δ14C) of soil extracellular DNA. Biol Fertil Soils. 2007;44:353–61.
Gutierrez-Cacciabue D, Cid AG, Rajal VB. How long can culturable bacteria and total DNA persist in environmental waters? The role of sunlight and solid particles. Sci Total Environ. 2016;539:494–502.
Carini P, Marsden PJ, Leff JW, Morgan EE, Strickland MS, Fierer N. Relic DNA is abundant in soil and obscures estimates of soil microbial diversity. Nat Microbiol. 2016;2:16242.
Kumar SS, Ghosh AR. Assessment of bacterial viability: a comprehensive review on recent advances and challenges. Microbiology (Reading). 2019;165:593–610.
Nocker A, Mazza A, Masson L, Camper AK, Brousseau R. Selective detection of live bacteria combining propidium monoazide sample treatment with microarray technology. J Microbiol Methods. 2009;76:253–61.
Ni J, Hatori S, Wang Y, Li YY, Kubota K. Uncovering viable microbiome in anaerobic sludge digesters by propidium monoazide (PMA)-PCR. Microb Ecol. 2020;79:925–32.
Stinson LF, Keelan JA, Payne MS. Characterization of the bacterial microbiome in first-pass meconium using propidium monoazide (PMA) to exclude nonviable bacterial DNA. Lett Appl Microbiol. 2019;68:378–85.
Fu X, Zeng B, Wang P, Wang L, Wen B, Li Y, Liu H, Bai S, Jia G. Microbiome of total versus live bacteria in the gut of rex rabbits. Front Microbiol. 2018;9:733.
Loublier C, Taminiau B, Heinen J, Lecoq L, Amory H, Daube G, Cesarini C. Evaluation of bacterial composition and viability of equine feces after processing for transplantation. Microorganisms. 2023;11:231.
Galazzo G, van Best N, Benedikter BJ, Janssen K, Bervoets L, Driessen C, Oomen M, Lucchesi M, van Eijck PH, Becker HEF, et al. How to count our microbes? The effect of different quantitative microbiome profiling approaches. Front Cell Infect Microbiol. 2020;10:403.
Legrand T, Wos-Oxley ML, Wynne JW, Weyrich LS, Oxley APA. Dead or alive: microbial viability treatment reveals both active and inactive bacterial constituents in the fish gut microbiota. J Appl Microbiol. 2021;131:2528–38.
Dvergedal H, Sandve SR, Angell IL, Klemetsdal G, Rudi K. Association of gut microbiota with metabolism in juvenile Atlantic salmon. Microbiome. 2020;8:160.
Bridge T, Cowman P, Bonito V, Head C, Hung JY-h, Rongo T, Baird A, Harii S, Sinniger F, Halafihi T, Quattrini A. A tenuis relationship: traditional taxonomy obscures systematics and biogeography of the ‘Acropora tenuis’ (Scleractinia: Acroporidae) species complex. Zool J Linn Soc. 2023.
Damjanovic K, Blackall LL, Peplow LM, van Oppen MJH. Assessment of bacterial community composition within and among Acropora loripes colonies in the wild and in captivity. Coral Reefs. 2020;39:1245–55.
Dungan AM, Hartman LM, Tortorelli G, Belderok R, Lamb AM, Pisan L, McFadden GI, Blackall LL, van Oppen MJH. Exaiptasia diaphana from the great barrier reef: a valuable resource for coral symbiosis research. Symbiosis. 2020;80:195–206.
Dungan AM, van Oppen MJH, Blackall LL. Short-term exposure to sterile seawater reduces bacterial community diversity in the sea anemone, Exaiptasia diaphana. Front Mar Sci. 2021;7:599314.
Andersson AF, Lindberg M, Jakobsson H, Backhed F, Nyren P, Engstrand L. Comparative analysis of human gut microbiota by barcoded pyrosequencing. PLoS ONE. 2008;3:e2836.
Bolyen E, Rideout JR, Dillon MR, Bokulich NA, Abnet CC, Al-Ghalith GA, Alexander H, Alm EJ, Arumugam M, Asnicar F, et al. Reproducible, interactive, scalable and extensible microbiome data science using QIIME 2. Nat Biotechnol. 2019;37:852–7.
Martin M. Cutadapt removes adapter sequences from high-throughput sequencing reads. EMBnet J. 2011;2011(17):3.
Callahan BJ, McMurdie PJ, Rosen MJ, Han AW, Johnson AJ, Holmes SP. DADA2: High-resolution sample inference from Illumina amplicon data. Nat Methods. 2016;13:581–3.
Bokulich NA, Kaehler BD, Rideout JR, Dillon M, Bolyen E, Knight R, Huttley GA, Gregory Caporaso J. Optimizing taxonomic classification of marker-gene amplicon sequences with QIIME 2’s q2-feature-classifier plugin. Microbiome. 2018;6:90.
Caporaso JG, Bittinger K, Bushman FD, DeSantis TZ, Andersen GL, Knight R. PyNAST: a flexible tool for aligning sequences to a template alignment. Bioinformatics. 2010;26:266–7.
Kõressaar T, Lepamets M, Kaplinski L, Raime K, Andreson R, Remm M. Primer3_masker: integrating masking of template sequence with primer design software. Bioinformatics. 2018;34:1937–8.
Salazar OR, Arun NP, Cui G, Bay LK, van Oppen MJH, Webster NS, Aranda M. The coral Acropora loripes genome reveals an alternative pathway for cysteine biosynthesis in animals. Sci Adv. 2022;8:eabq0304.
Hartman LM, Blackall LL, van Oppen MJH. Antibiotics reduce bacterial load in Exaiptasia diaphana, but biofilms hinder its development as a gnotobiotic coral model. Access Microbiol 2022; 4.
R Core Team. R: a language and environment for statistical computing. Vienna: R Foundation for Statistical Computing, Austria; 2019.
McMurdie PJ, Holmes S. phyloseq: an R package for reproducible interactive analysis and graphics of microbiome census data. PLoS ONE. 2013;8:e61217.
Davis NM, Proctor DM, Holmes SP, Relman DA, Callahan BJ. Simple statistical identification and removal of contaminant sequences in marker-gene and metagenomics data. Microbiome. 2018;6:1–14.
Pinheiro J, Bates D, DebRoy S, Sarkar D, Team RC. Linear and nonlinear mixed effects models. R Package Version. 2007;3:1–89.
Searle SR, Speed FM, Milliken GA. Population marginal means in the linear model: an alternative to least squares means. Am Stat. 1980;34:216–21.
Wickham H. ggplot2: elegant graphics for data analysis version 3.3.1. New York: Springer; 2020.
Lemon J, Bolker B, Oom S, Klein E, Rowlingson B, Wickham H, Tyagi A, Eterradossi O, Grothendieck G, Toews M. Package ‘plotrix.’ Vienna: R Development Core Team; 2015.
Auguie B, Antonov A, Auguie MBJMFfGG. Package ‘gridExtra’. 2017.
Hope RM, Hope MRM, Collate'CI R. Package ‘Rmisc’. 2013, 101:2.
Oksanen J. Vegan: community ecology package. R package version 1.8–5. http://www.cran.r-project.org 2007.
Arbizu PM. JShgcpp: pairwiseAdonis: Pairwise multilevel comparison using adonis R package version 0.3. 2019.
Cáceres MD, Legendre P. Associations between species and groups of sites: indices and statistical inference. Ecology. 2009;90:3566–74.
Larsson J. Eulerr: area-proportional Euler diagrams with ellipses. 2018.
Lahti L, Shetty S. Introduction to the microbiome R package. 2018.
Sogin ML, Morrison HG, Huber JA, Welch DM, Huse SM, Neal PR, Arrieta JM, Herndl GJ. Microbial diversity in the deep sea and the underexplored “rare biosphere.” Proc Natl Acad Sci. 2006;103:12115–20.
Shiu JH, Yu SP, Fong CL, Ding JY, Tan CJ, Fan TY, Lu CY, Tang SL. Shifting in the Dominant Bacterial Group Endozoicomonas Is Independent of the Dissociation With Coral Symbiont Algae. Front Microbiol. 2020;11:1791.
Zou Y, Chen Y, Wang L, Zhang S, Li J. Differential responses of bacterial communities in coral tissue and mucus to bleaching. Coral Reefs. 2022;41:951–60.
Mei R, Narihiro T, Nobu MK, Kuroda K, Liu WT. Evaluating digestion efficiency in full-scale anaerobic digesters by identifying active microbial populations through the lens of microbial activity. Sci Rep. 2016;6:34090.
Sun F, Yang H, Wang G, Shi Q. Combination analysis of metatranscriptome and metagenome reveal the composition and functional response of coral symbionts to bleaching during an El Niño event. Front Microbiol. 2020;11:448.
Marchioro GM, Glasl B, Engelen AH, Serrao EA, Bourne DG, Webster NS, Frade PR. Microbiome dynamics in the tissue and mucus of acroporid corals differ in relation to host and environmental parameters. PeerJ. 2020;8:e9644.
Damjanovic K, van Oppen MJH, Menéndez P, Blackall LL. Experimental inoculation of coral recruits with marine bacteria indicates scope for microbiome manipulation in Acropora tenuis and Platygyra daedalea. Front Microbiol. 2019;10:1702.
Rivera-Ortega J, Thomé PE. Contrasting antibacterial capabilities of the surface mucus layer from three symbiotic cnidarians. Front Mar Sci. 2018;5:392.
Reyneke B, Waso M, Ndlovu T, Clements T, Havenga B, Khan S, Khan W. EMA- Versus PMA-amplicon-based sequencing to elucidate the viable bacterial community in rainwater. Water Air Soil Pollut. 2022;233:103.
Chan WY, Peplow LM, Menendez P, Hoffmann AA, van Oppen MJH. The roles of age, parentage and environment on bacterial and algal endosymbiont communities in Acropora corals. Mol Ecol. 2019;28:3830–43.
Tandon K, Chiou YJ, Yu SP, Hsieh HJ, Lu CY, Hsu MT, Chiang PW, Chen HJ, Wada N, Tang SL. Microbiome restructuring: dominant coral bacterium Endozoicomonas species respond differentially to environmental changes. mSystems. 2022;7:e0035922.
Ochsenkuhn MA, Schmitt-Kopplin P, Harir M, Amin SA. Coral metabolite gradients affect microbial community structures and act as a disease cue. Commun Biol. 2018;1:184.
Gignoux-Wolfsohn SA, Aronson FM, Vollmer SV. Complex interactions between potentially pathogenic, opportunistic, and resident bacteria emerge during infection on a reef-building coral. FEMS Microbiol Ecol. 2017;93:fix080.
Ide K, Nakano Y, Ito M, Nishikawa Y, Fujimura H, Takeyama H. Coral co-culture alters bacterial community structure by altering rare bacterial species. Res Sq 2022.
Quigley KM, Alvarez Roa C, Torda G, Bourne DG, Willis BL. Co-dynamics of Symbiodiniaceae and bacterial populations during the first year of symbiosis with Acropora tenuis juveniles. Microbiologyopen. 2020;9:e959.
Glasl B, Smith CE, Bourne DG, Webster NS. Disentangling the effect of host-genotype and environment on the microbiome of the coral Acropora tenuis. PeerJ. 2019;7:e6377.
Wada N, Hsu MT, Tandon K, Hsiao SS, Chen HJ, Chen YH, Chiang PW, Yu SP, Lu CY, Chiou YJ, et al. High-resolution spatial and genomic characterization of coral-associated microbial aggregates in the coral Stylophora pistillata. Sci Adv. 2022;8:eabo2431.
Neave MJ, Rachmawati R, Xun L, Michell CT, Bourne DG, Apprill A, Voolstra CR. Differential specificity between closely related corals and abundant Endozoicomonas endosymbionts across global scales. ISME J. 2017;11:186–200.
Ding J-Y, Shiu J-H, Chen W-M, Chiang Y-R, Tang S-L. Genomic insight into the host-endosymbiont relationship of Endozoicomonas montiporae CL-33T with its coral host. Front Microbiol. 2016;7:251.
Hochart C, Paoli L, Ruscheweyh HJ, Salazar G, Boissin E, Romac S, Poulain J, Bourdin G, Iwankow G, Moulin C, et al. Ecology of endozoicomonadaceae in three coral genera across the Pacific Ocean. Nat Commun. 2023;14:3037.
Casas V, Kline DI, Wegley L, Yu Y, Breitbart M, Rohwer F. Widespread association of a Rickettsiales-like bacterium with reef-building corals. Environ Microbiol. 2004;6:1137–48.
Libro S, Vollmer SV. Genetic signature of resistance to white band disease in the caribbean staghorn coral Acropora cervicornis. PLoS ONE. 2016;11:e0146636.
Ide K, Nishikawa Y, Kogawa M, Iwamoto E, Samuel AZ, Nakano Y, Takeyama H. High-Quality Draft Genome Sequence of a Rickettsiales Bacterium Found in Acropora tenuis Coral from Okinawa, Japan. Microbiol Resour Announc 2020, 9.
Godoy-Vitorino F, Ruiz-Diaz CP, Rivera-Seda A, Ramirez-Lugo JS, Toledo-Hernandez C. The microbial biosphere of the coral Acropora cervicornis in Northeastern Puerto Rico. PeerJ. 2017;5:e3717.
Miller N, Maneval P, Manfrino C, Frazer TK, Meyer JL. Spatial distribution of microbial communities among colonies and genotypes in nursery-reared Acropora cervicornis. PeerJ. 2020;8:e9635.
Klinges JG, Rosales SM, McMinds R, Shaver EC, Shantz AA, Peters EC, Eitel M, Worheide G, Sharp KH, Burkepile DE, et al. Phylogenetic, genomic, and biogeographic characterization of a novel and ubiquitous marine invertebrate-associated Rickettsiales parasite, Candidatus Aquarickettsia rohweri, gen. nov., sp. nov. ISME J. 2019;13:2938–53.
Gignoux-Wolfsohn SA, Vollmer SV. Identification of candidate coral pathogens on white band disease-infected staghorn coral. PLoS ONE. 2015;10:e0134416.
Shaver EC, Shantz AA, McMinds R, Burkepile DE, Vega Thurber RL, Silliman BR. Effects of predation and nutrient enrichment on the success and microbiome of a foundational coral. Ecology. 2017;98:830–9.
Baker LJ, Reich HG, Kitchen SA, Grace Klinges J, Koch HR, Baums IB, Muller EM, Thurber RV. The coral symbiont Candidatus Aquarickettsia is variably abundant in threatened Caribbean acroporids and transmitted horizontally. ISME J. 2022;16:400–11.
Meunier V, Bonnet S, Pernice M, Benavides M, Lorrain A, Grosso O, Lambert C, Houlbrèque F. Bleaching forces coral’s heterotrophy on diazotrophs and Synechococcus. ISME J. 2019;13:2882–6.
Tremblay P, Naumann MS, Sikorski S, Grover R, Ferrier-Pagès C. Experimental assessment of organic carbon fluxes in the scleractinian coral Stylophora pistillata during a thermal and photo stress event. Mar Ecol Prog Ser. 2012;453:63–77.
Ducklow HW. The biomass, production and fate of bacteria in coral reefs. In: Dubinsky Z, editor. Ecosystems of the world: coral reefs, vol. 25. Amsterdam: Elsevier; 1990. p. 265–89.
Angly FE, Heath C, Morgan TC, Tonin H, Rich V, Schaffelke B, Bourne DG, Tyson GW. Marine microbial communities of the Great Barrier Reef lagoon are influenced by riverine floodwaters and seasonal weather events. PeerJ. 2016;4:e1511.
Pratte ZA, Richardson LL. Possible links between white plague-like disease, scleractinian corals, and a cryptochirid gall crab. Dis Aquat Organ. 2016;122:153–61.
Sweet M, Bythell J. White syndrome in Acropora muricata: nonspecific bacterial infection and ciliate histophagy. Mol Ecol. 2015;24:1150–9.
Rogers GB, Cuthbertson L, Hoffman LR, Wing PA, Pope C, Hooftman DA, Lilley AK, Oliver A, Carroll MP, Bruce KD, van der Gast CJ. Reducing bias in bacterial community analysis of lower respiratory infections. ISME J. 2013;7:697–706.
Dube CE, Ziegler M, Merciere A, Boissin E, Planes S, Bourmaud CA, Voolstra CR. Naturally occurring fire coral clones demonstrate a genetic and environmental basis of microbiome composition. Nat Commun. 2021;12:6402.
Wada N, Ishimochi M, Matsui T, Pollock FJ, Tang SL, Ainsworth TD, Willis BL, Mano N, Bourne DG. Characterization of coral-associated microbial aggregates (CAMAs) within tissues of the coral Acropora hyacinthus. Sci Rep. 2019;9:14662.
Ainsworth TD, Fine M, Blackall LL, Hoegh-Guldberg O. Fluorescence in situ hybridization and spectral imaging of coral-associated bacterial communities. Appl Environ Microbiol. 2006;72:3016–20.
Ainsworth TD, Hoegh-Guldberg O. Bacterial communities closely associated with coral tissues vary under experimental and natural reef conditions and thermal stress. Aquat Biol. 2009;4:289–96.
Bayer T, Neave MJ, Alsheikh-Hussain A, Aranda M, Yum LK, Mincer T, Hughen K, Apprill A, Voolstra CR. The microbiome of the Red Sea coral Stylophora pistillata is dominated by tissue-associated Endozoicomonas bacteria. Appl Environ Microbiol. 2013;79:4759–62.
Work TM, Aeby GS. Microbial aggregates within tissues infect a diversity of corals throughout the Indo-Pacific. Mar Ecol Prog Ser. 2014;500:1–9.
Work TM, Aeby GS, Lasne G, Tribollet A. Gross and microscopic pathology of hard and soft corals in New Caledonia. J Invertebr Pathol. 2014;120:50–8.
Aeby GS, Shore A, Jensen T, Ziegler M, Work T, Voolstra CR. A comparative baseline of coral disease in three regions along the Saudi Arabian coast of the central Red Sea. PLoS ONE. 2021;16:e0246854.
Rodriguez-Villalobos JC, Rocha-Olivares A, Work TM, Calderon-Aguilera LE, Caceres-Martinez JA. Gross and microscopic pathology of lesions in Pocillopora spp from the subtropical eastern Pacific. J Invertebr Pathol. 2014;120:9–17.
Work TM, Aeby GS, Hughen KA. Gross and microscopic lesions in corals from micronesia. Vet Pathol. 2016;53:153–62.
Maire J, Tandon K, Collingro A, van de Meene A, Damjanovic K, Gotze CR, Stephenson S, Philip GK, Horn M, Cantin NE, et al. Colocalization and potential interactions of Endozoicomonas and chlamydiae in microbial aggregates of the coral Pocillopora acuta. Sci Adv. 2023;9:eadg0773.
Neave MJ, Apprill A, Ferrier-Pages C, Voolstra CR. Diversity and function of prevalent symbiotic marine bacteria in the genus Endozoicomonas. Appl Microbiol Biotechnol. 2016;100:8315–24.
Gardner SG, Nitschke MR, O’Brien J, Motti CA, Seymour JR, Ralph PJ, Petrou K, Raina J-B. Increased DMSP availability during thermal stress influences DMSP-degrading bacteria in coral mucus. Front Mar Sci. 2022;9:912862.
Kim TG, Jeong S-Y, Cho K-S. Comparison of droplet digital PCR and quantitative real-time PCR for examining population dynamics of bacteria in soil. Appl Microbiol Biotechnol. 2014;98:6105–13.
Rosales SM, Huebner LK, Clark AS, McMinds R, Ruzicka RR, Muller EM. Bacterial metabolic potential and micro-eukaryotes enriched in stony coral tissue loss disease lesions. Front Mar Sci. 2022;8:776859.
Dungan AM, Tandon K, Jameson V, Gotze CR, Blackall LL, van Oppen MJH. A targeted approach to enrich host-associated bacteria for metagenomic sequencing. FEMS Microb 2023 (in press, accepted 15/11/23).
Marotz CA, Sanders JG, Zuniga C, Zaramela LS, Knight R, Zengler K. Improving saliva shotgun metagenomics by chemical host DNA depletion. Microbiome. 2018;6:42.
Chen S, Wang F, Beaulieu JC, Stein RE, Ge B. Rapid detection of viable salmonellae in produce by coupling propidium monoazide with loop-mediated isothermal amplification. Appl Environ Microbiol. 2011;77:4008–16.
Dean FB, Nelson JR, Giesler TL, Lasken RS. Rapid amplification of plasmid and phage DNA using Phi 29 DNA polymerase and multiply-primed rolling circle amplification. Genome Res. 2001;11:1095–9.
Marotz C, Morton JT, Navarro P, Coker J, Belda-Ferre P, Knight R, Zengler K. Quantifying live microbial load in human saliva samples over time reveals stable composition and dynamic load. mSystems. 2021;6:10–1128.
Weinmaier T, Probst AJ, La Duc MT, Ciobanu D, Cheng JF, Ivanova N, Rattei T, Vaishampayan P. A viability-linked metagenomic analysis of cleanroom environments: eukarya, prokaryotes, and viruses. Microbiome. 2015;3:62.
Mohapatra BR, La Duc MT. Rapid detection of viable Bacillus pumilus SAFR-032 encapsulated spores using novel propidium monoazide-linked fluorescence in situ hybridization. J Microbiol Methods. 2012;90:15–9.
Papanicolas LE, Wang Y, Choo JM, Gordon DL, Wesselingh SL, Rogers GB. Optimisation of a propidium monoazide based method to determine the viability of microbes in faecal slurries for transplantation. J Microbiol Methods. 2019;156:40–5.
Wang Y, Yan Y, Thompson KN, Bae S, Accorsi EK, Zhang Y, Shen J, Vlamakis H, Hartmann EM, Huttenhower C. Whole microbial community viability is not quantitatively reflected by propidium monoazide sequencing approach. Microbiome. 2021;9:17.
Silva DP, Epstein HE, Vega Thurber RL. Best practices for generating and analyzing 16S rRNA amplicon data to track coral microbiome dynamics. Front Microbiol. 2022;13:1007877.
Acknowledgements
The corals used in this study were originally collected from the sea estates of the Wuthathi, Bindal, Wulgurukaba, and Manbarra People. We acknowledge their contributions to this research not only by way of the coral samples, but also through Indigenous Knowledge and custodianship of the land and sea country on which we work. We are grateful to Lonidas Koukoumaftsis, Dr. Jose Montalvo-Proano, Dr. Thomas (Ed) Roberts, Rachel Neil, and Hugo Scharfenstein for coral material.
Funding
This research was funded by a University of Melbourne Early Career Researcher Grant (to AMD) and the Australian Research Council (FL180100036 to MJHvO; DP210100630 to MJHvO and LLB). Funding bodies had no influence in the design of the study, the collection, analysis, and interpretation of data, or in writing the manuscript.
Author information
Authors and Affiliations
Contributions
AMD conceived the project; AMD, LLB, and MJHvO designed the study; AMD, LG, and AW collected and processed samples through PMA treatment and metabarcoding library prep; ECF completed ddPCR; AMD analysed and interpreted the data; AMD wrote the first draft; all authors substantively revised the manuscript and approve of the submitted version.
Corresponding author
Ethics declarations
Ethics approval and consent to participate
Not applicable.
Consent for publication
Not applicable.
Competing interests
The authors declare that they have no competing interests.
Additional information
Publisher's Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Supplementary Information
Additional file 1.
Supplementary Tables S1–S4. Supplementary Table S1: Collection and sampling details for all coral colonies used in this study. Supplementary Table S2: Putative contaminants in the metabarcoding dataset identified with decontam. Supplementary Table S3: β-diversity pairwise comparisons by genotype for A. kenti, P. daedalea, and P. lutea. Supplementary Table S4: Review of studies examining non-viable versus viable bacteria with PMA in mixed communities.
Additional file 2.
Supplementary Figures S1–S5. Supplementary Figure S1: Box plots showing 16S rRNA gene copies μl−1 of method validation samples from ddPCR. Supplementary Figure S2: Rarefaction curves for each of six coral species sampled. Supplementary Figure S3: Alpha diversity indices by coral genotype for A. kenti and P. daedalea. Supplementary Figure S4: Relative abundance of bacterial genera for A. millepora, A. kenti, P. daedalea, and P. lutea by genotype. Supplementary Figure S5: Relative abundance of Endozoicomonas and unknown Rhodanobacteraceae ASVs by coral genotype.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.
About this article
Cite this article
Dungan, A.M., Geissler, L., Williams, A.S. et al. DNA from non-viable bacteria biases diversity estimates in the corals Acropora loripes and Pocillopora acuta. Environmental Microbiome 18, 86 (2023). https://doi.org/10.1186/s40793-023-00541-6
Received:
Accepted:
Published:
DOI: https://doi.org/10.1186/s40793-023-00541-6