Diuron tolerance and potential degradation by pelagic microbiomes in the Great Barrier Reef lagoon

Diuron is a herbicide commonly used in agricultural areas where excess application causes it to leach into rivers, reach sensitive marine environments like the Great Barrier Reef (GBR) lagoon and pose risks to marine life. To investigate the impact of diuron on whole prokaryotic communities that underpin the marine food web and are integral to coral reef health, GBR lagoon water was incubated with diuron at environmentally-relevant concentration (8 µg/L), and sequenced at specific time points over the following year. 16S rRNA gene amplicon profiling revealed no significant short- or long-term effect of diuron on microbiome structure. The relative abundance of prokaryotic phototrophs was not significantly altered by diuron, which suggests that they were largely tolerant at this concentration. Assembly of a metagenome derived from waters sampled at a similar location in the GBR lagoon did not reveal the presence of mutations in the cyanobacterial photosystem that could explain diuron tolerance. However, resident phages displayed several variants of this gene and could potentially play a role in tolerance acquisition. Slow biodegradation of diuron was reported in the incubation flasks, but no correlation with the relative abundance of heterotrophs was evident. Analysis of metagenomic reads supports the hypothesis that previously uncharacterized hydrolases carried by low-abundance species may mediate herbicide degradation in the GBR lagoon. Overall, this study offers evidence that pelagic phototrophs of the GBR lagoon may be more tolerant of diuron than other tropical organisms, and that heterotrophs in the microbial seed bank may have the potential to degrade diuron and alleviate local anthropogenic stresses to inshore GBR ecosystems.


INTRODUCTION
Coral reefs are very rich and diverse ecosystems, though due to both local and global anthropogenic disturbances, they are in a state of gradual decline (Pandolfi et al., 2003).
A recent study found that microbial community variations in the GBR lagoon are primarily driven by riverine effluents (Angly et al., 2016) and a year-long seawater simulation experiment presented evidence that microorganisms play a role in the biodegradation of a wide range of PSII herbicides in this region (Mercurio et al., 2015). Although the community-wide effects of diuron on marine microorganisms are not characterized, this research suggests that pelagic microbiomes may protect coral reef and seagrass ecosystems by degrading this herbicide. In the present study, we hypothesized that diuron significantly affects the composition of these microbiomes by specifically: (i) inhibiting phototrophs in the short-term; and (ii) sustaining the long-term growth of selected heterotrophs that have the potential to metabolize it, leading to the herbicide's disappearance. To investigate these hypotheses, we collected samples from the simulation study of Mercurio et al. (2015) and an inshore GBR location, and characterized their microbiomes (Archaea and Bacteria) using high-throughput 16S rRNA amplicon profiling and metagenomics. Mercurio et al. (2015) collected surface seawater (24 L) in sterile containers in the tropical dry season (15 May 2012) from Cape Ferguson, QLD, Australia (latitude −19.2673297, longitude 147.0591537) (Fig. S1), a site that is ∼17 km downstream from the Haughton River mouth and where diuron is consistently reported in the wet season (Lewis et al., 2009;Kennedy et al., 2012a;Kennedy et al., 2012b). The investigators passed seawater through 20 µm impact filters, dispensed it in 500 mL glass flasks and incubated it during 365 days on a shaking platform at 25 • C either in the dark or in the light (12:12 light day cycle with 40 µmol photons m −2 s −1 ), and with or without amendment of diuron (at the ecologically relevant concentration of 8 µg/L) (Lewis et al., 2009;Kennedy et al., 2012a;Kennedy et al., 2012b) (Figs. S2A and S2B). The investigators performed each experimental treatment in triplicate (12 flasks in total) and monitored diuron concentration for each flask over the life of the incubation. For more details of the experimental set up and diuron degradation results, see Mercurio et al. (2015).

Cape Ferguson diuron incubations
In the present study, we collected subsamples (705 µL) from each flask of the Cape Ferguson incubation experiment at day 0, 2, 7, 28, 120 and 365 for 16S rRNA gene analysis. Each subsample was collected with a pipette after thoroughly shaking the flask, 5 µL were stained with 5 µL of DAPI (1 µg/mL) in the dark for 15 min, and observed with a Nikon Ci-L epifluorescence microscope (Fig. S2C). This confirmed the presence of DNA-containing cells, a prerequisite for sequencing.
Total DNA was extracted by first centrifuging each sample at 13,000 × g for 30 min. Each resulting pellet was then resuspended in 20 µL microLYSIS-Plus DNA release buffer (Microzone, West Sussex, UK) for 30 min at room temperature and incubated using a thermal cycler following the manufacturer's tough cell lysis protocol (65 • C for 15 min; 96 • C for 2 min; 65 • C for 4 min; 96 • C for 1 min; 65 • C for 1 min; 96 • C for 30 s). Control of the absence of microLYSIS-Plus reagent contamination by foreign DNA was carried out by adding a blank sample, containing only the microLYSIS-Plus buffer (without template DNA).
Microbial amplicons were generated by PCR-amplifying the V6-V8 variable regions of the 16S rRNA gene in the total DNA using a universal primer set targeting Archaea and Bacteria (iTAG926F and iTAG1392wR primers) according to the protocol in Dove et al. (2013). These amplicons were paired-end sequenced on an Illumina MiSeq instrument at the Australian Centre for Ecogenomics (ACE).

Statistical analysis
The significance of changes in taxon relative abundance between sampling points were evaluated using LEfSe (Segata et al., 2011). The effects of incubation regimen on microbial community structure was assessed by Hellinger-transforming the microbial profiles and analyzing them using the R language (R Foundation for Statistical Computing, Vienna, Austria, 0000), specifically using the capscale(), adonis() and rda() functions of the vegan packages (Dixon, 2003), for PCoA, PERMANOVA and RDA analysis respectively.

Dunk Island metagenome preparation
An additional seawater sample was collected in the tropical dry season (13 October 2009), north of Dunk Island, QLD, Australia (latitude −17.9242918, longitude 146.1429637) (Fig. S1). This site is ∼15 km downstream from the Tully River mouth, and exposed to diuron and other PSII herbicides in comparable concentrations as the Cape Ferguson site (Lewis et al., 2009;Kennedy et al., 2012a;Kennedy et al., 2012b). A 20 L volume was taken from a depth of 5 m and pre-filtered through a 2.7 µm Whatman GF/D filter and a 1.6 µm Whatman GF/A filter to remove particles and most eukaryotic microorganisms. The filtrate was then passed through a 0.22 µm Millipore Express Plus filter to capture the bacterial and archaeal fraction. The filters were folded in half, cells inward, added to a tube containing 20 mL of lysis buffer (40 mM Na 2 EDTA, 50 mM Tris pH 8.3 and 0.73 M sucrose, sterilized), stored shipboard at −20 • C and transferred to −80 • C on land.
DNA was extracted from the filter using a modified method from Suzuki et al. (2004). In brief, the filter was thawed on ice, added 6 mL of lysis buffer with 5 mg/mL lysozyme and the tube was incubated for 30 min at 37 • C, while rotating at 10 rpm. Proteinase K (1.1 mg/mL final concentration) and 10% sodium dodecyl sulfate (1.1% final concentration) were added and the sample was incubated at 55 • C for 2 h, with rotation. The lysate was split in half and DNA was extracted from each using two rounds of phenol:chloroform:isoamyl alcohol (25:24:1, pH 8.0), then one round of chloroform:isoamyl alcohol (24:1). Aqueous phases were pooled and frozen overnight at −20 • C. The aqueous phase was then cleaned by passage of 15 mL at a time through Amicon Ultra-15 100 kDa spin unit (EMD Millipore, Billerica, MA, USA). The filter was washed once with 8 mL of Tris EDTA buffer (TE, 10 mM, pH 8.0) and recovered with 50 µL of TE (1 mM, pH 8.0). DNA was then further cleaned by precipitation with 70% ethanol, the pellet was washed once with 70% ethanol, air dried, and resuspended in 100 µL TE, for a total yield of 37 µg DNA. The resulting DNA was sequenced on an Illumina (Solexa) Genome Analyzer II instrument at the University of Arizona, producing 25.4 million pairs of 101 bp long reads.
Under the assumption that each distinct identified PuhAB protein is encoded by a different species (the three known Puh proteins are encoded on three distinct genomes), we approximated this species' relative abundance as: where R is the number of reads matching puhAB, M is the number of screened metagenomic reads (22,927,633), L is the average read length L (93.7 bp), P is the average length of the puhAB genes (1,376 bp), G is the average genome length in marine microbiomes (2.58 Mbp) (Angly et al., 2009), and S is the number of species in inshore GBR water column (643 OTUs in the diuron incubation experiment).

Contig-centric metagenomic screening for photosystem genes
For this analysis, the Dunk Island metagenomic raw read pairs were cleaned with TRIMMOMATIC by removing Illumina adapters, deleting reads with uncalled bases, truncating their 5 end to a final length of 80 bp, and removing smaller reads. The data were assembled using IDBA-UD (Peng et al., 2012), and the resulting scaffolds translated into their six possible reading frames. The hmmsearch tool of HMMER3 (Eddy, 2011) was employed to look for photosystem B proteins in these translated scaffolds using the TIGR001151 PsbA hidden Markov Model profile of TIGRFAMs (Haft, Selengut & White, 2003). A maximum E-value of 1e-50 was used to retrieve significant matches and their alignment was visualized in Jalview. The taxonomic affiliation of the scaffolds matching PsbA was determined by best BLASTN similarity against the NCBI nr database (minimum identity of 70% over a minimum alignment length of 1,200 bp, i.e., the length of PsbA + 40 amino acids). Nesoni (https://github.com/Victorian-Bioinformatics-Consortium/nesoni) and SHRiMP (Rumble et al., 2009) were used to map the Illumina reads against the metagenomic scaffolds and call single nucleotide polymorphisms (SNPs). Mercurio et al. (2015) collected seawater during the tropical dry season at Cape Ferguson, an inshore region of the GBR (Figs. S1 and S2), to conduct a year-long diuron incubation experiment. In the present study, the 16S rRNA amplicon sequencing of 72 samples taken at set time points from the incubation flasks generated a total of 4.96 million read pairs (NCBI accession PRJNA276057). Processing through the Hitman bioinformatic pipeline resulted in 3.83 million high-quality sequences (77.3% of the initial amount). Rarefaction at a depth of 20,000 counts per sample provided a sequencing depth-independent view of the diversity of the samples (Fig. S3, Table S1), collectively containing 4,743 distinct OTUs (97% identity level).

Microbial dynamics in Cape Ferguson diuron incubation
The taxonomic affiliation conducted by best global alignment against the Greengenes database and subsequent gene-copy number correction (Angly et al., 2014) permitted estimation of changes in the relative abundance of prokaryotic taxa over time (Fig. S4). When averaging the replicates (Fig. 1, Fig. S5), the most abundant taxa at the start of the incubation (day 0) included the orders Rickettsiales (19%) and Synechococcales (14% average relative abundance), from the Cyanobacteria and Proteobacteria phyla, respectively. The microbiomes were marked by a succession of various taxa over time, as seen in previous work (Fierer et al., 2010). For example, Sphingomonadales increased significantly from an initial average of 0.38% relative abundance (day 0) to dominate the communities with 30% at day 2 (LEfSe; α < 0.05). Rhodobacterales-affiliated sequences increased significantly, reaching a maximum relative abundance of 36% on average a week after the start of the incubation (day 7) (LEfSe; α < 0.05), and subsequently significantly decreased until day 120 (LEfSe; α < 0.05). At the end of the incubation experiment (day 365), Oceanospirillales were very abundant in the control flasks exposed to light, while Thiotrichales dominated the samples incubated in the dark (both control and diuron-treated).
Three predominant OTUs (OTU 12, 13 and 20) characteristic of the flasks kept in the dark could not be assigned to a taxonomic group. Further identification efforts using the RDP Classifier and Silva's SINA suggest that they all belong to the Proteobacteria phylum, more precisely to the Salinisphaera, Coxiella and GR-WP33-30 taxa (Table S2). The genus Salinisphaera includes a recently sequenced species that is adapted to environments with fluctuating conditions (Antunes et al., 2011), while the genus Coxiella contains a single species that is highly resistant to environmental stresses such as temperature, osmotic pressure and ultraviolet radiation (Voth & Heinzen, 2007), and representatives of the order GR-WP33-30 were detected in uranium mines (Selenska-Pobell & Radeva, 2004). The robustness of these taxa may be responsible for their success in the dark and likely oligotrophic conditions of the incubation flasks.

Effect of diuron on microbial profiles
The diuron measurements made by Mercurio et al. (2015) in the incubation flasks ranged from an initial 8.77 µg/L (dark conditions, replicate R4) down to 3.78 µg/L (light conditions, replicate R3) after one year of incubation. Here, we included the diuron concentration of each individual flask as an input for a constrained ordination (Fig. 2), which demonstrated a significant influence of incubation time and light exposure, but not of diuron concentrations on the microbial profiles (PERMANOVA, p < 0.05). Dissection of the differences between diuron-treated and control flasks for each individual sampling day, PCoA (Fig. S6) confirmed that diuron did not affect microbiome composition significantly (PERMANOVA, p < 0.05).

Resistance of photosynthetic bacteria to diuron
Some phototrophic bacteria are inhibited by diuron, while others are insensitive. For example, photosystems I and II exist in Cyanobacteria and vascular plants, and physiological experiments have demonstrated binding of diuron on the cyanobacterial photosystem II, leading to photosynthesis inhibition (Allen et al., 1983;Gadkari, 1988;Brusslan & Haselkorn, 1989;Deng, Gao & Sun, 2012). Conversely, diuron does not bind to the photosynthetic reaction center of purple bacteria and they may remain unaffected (Sinning, 1992). We thus hypothesized that the majority of phototrophic prokaryotes in the incubation experiment would be affected by diuron toxicity, resulting in their rapid decline.

Figure 2 OTU-level RDA of the microbiomes (Hellinger-based) in the Cape Ferguson diuron incubations.
OTUs are indicated by a red cross and the Greengenes taxonomic affiliation of the most discriminating is shown. Circle size is proportional to incubation time (sampling day). Environmental factors are green arrows depicting light amount in the light and dark treatments, incubation time and measured diuron concentration. Asterisks denote environmental factors that are statistically significant (PER-MANOVA; p < 0.05). Samples from day 2 and 7 were omitted from this analysis because diuron concentration was not measured on these days.
Cyanobacteria and purple bacteria such as Rhodobacteraceae were prevalent in the incubation flasks, but their relative abundance did not decline between day 0 and 28 (Fig. 2,  Figs. S6A-S6D), despite the presence of more than 8.45 µg/L diuron on average during this period, a concentration that markedly inhibits the photosynthesis of diatoms and green algae (Magnusson, Heimann & Negri, 2008;Magnusson et al., 2012). This supports previous reports that Cyanobacteria are less sensitive to PSII herbicides than eukaryotic phototrophs (Lürling & Roessink, 2006). The relative insensitivity of Cyanobacteria in our dataset could be explained by pollution-induced community tolerance (PICT) following chronic exposure to herbicides, which was previously reported for biofilms in a French river (Tlili et al., 2008;Tlili et al., 2011) and for periphyton (a mixture of detritus, algae and microorganisms growing on submerged surfaces) in the GBR lagoon (Magnusson et al., 2012) and in a Swedish fjord (Molander & Blanck, 1992). The mechanism underpinning this tolerance for diuron is not yet elucidated, but may be related to the evolution and enrichment of high-turnover variants of the PsbA protein upon which diuron and other PSII herbicides such as irgarol 1051 bind (Eriksson et al., 2009;Deng, Gao & Sun, 2012).
Metagenomic analysis was undertaken to explore the presence of PsbA variants and the potential for diuron resistance in the GBR lagoon. Since the samples collected during the Cape Ferguson incubation experiment contained too little biomass for comprehensive metagenomic sequencing, we prepared a metagenome from a sample collected during the tropical dry season at Dunk Island, another inshore GBR location (NCBI accession SRR1819825). Weather, river effluent and diuron exposure data indicate that the Dunk Island and Cape Ferguson samples were both representative of the GBR lagoon during the dry season, when the effects of riverine floodwaters are minimal (Text S1, Tables S3-S6), and therefore comparable (Angly et al., 2016). The Dunk Island metagenome was assembled into ∼74,000 scaffolds (771 bp average length, 879 bp N50) from which all putative PsbA protein sequences were identified (Fig. 3). The introduction of Val219 and Ser264 mutations in PsbA confers PSII herbicide resistance (Bettini et al., 1987;Mengistu et al., 2000) and mutations in the PsbA PEST domain (rich in amino-acids P, E, S and T) were previously correlated with resistance in the environment (Eriksson et al., 2009). But none of these mutations were detected in the Dunk Island scaffolds of Cyanobacteria (Synechococcus and Prochlorococcus). Further, a total of ∼24,900 reads mapped onto the 12.9 kb long Prochlorococcus scaffold, but no SNPs could be identified within the psbA gene.
PEST sequence type 11 (RETTENESANAGYK), representing a PEST type hypothesized to confer irgarol tolerance to the Swedish fjord microbiomes (Eriksson et al., 2009), was detected in the metagenomic scaffolds (Fig. 3). However, BLASTN analyses suggest that this sequence was part of a eukaryotic genome (100% query coverage and 100% identity to Micromonas, a member of the Prasinophyceae class) (Worden et al., 2009). The Prasinophyceae are composed of unicellular photosynthetic green algae and Micromonas (<2 µm) could have passed through the 1.6 µm wide pores of the filters used during metagenome preparation. Similarly, the Swedish fjord sequencing read containing PEST type 11 (accession AM933747) best matched a eukaryotic genome from another subdivision of the Prasinophyceae (100% query coverage and 99% identity to the Pycnococcaceae family). This suggests that PEST type 11 is a general feature of the Prasinophyceae genome and, at least in the present study, not an adaptive mutation of Bacteria and Archaea to protect against PSII herbicides.
BLASTN investigation of the metagenomic scaffolds from Dunk Island encoding PsbA revealed that they were not only prokaryotic and eukaryotic. Most of them (14 out of 24) were of viral origin (Fig. 3), and 9 out of the 11 PEST types identified in the present study matched some carried by cyanophages. This large PEST type diversity and the propensity of phages to transfer genes to and from their hosts (lateral gene transfer) raises the possibility that cyanophages steer the stability of PsbA in their hosts (Lindell et al., 2004;Zeidner et al., 2005;Sharon et al., 2007). Future research should consider how this may influence the tolerance of phototrophs to PSII herbicides.

Potential for degradation by heterotrophic bacteria
Another hypothesis formulated in this study was that specific heterotrophic populations would carry genes for the degradation of diuron and take advantage of this resource, leading to their increase in relative abundance over time. During the one-year Cape Ferguson incubation experiment, Mercurio et al. (2015) reported 15-31% diuron degradation. They attributed this slow degradation in part to prokaryotic breakdown but, in the present work, we found no significant association between diuron-treated incubation flasks and heterotrophic abundance between day 28 and 365 (Fig. 2,. While this evidence goes against our hypothesis of rapid heterotrophic degradation, a similar marine incubation study also detected a lack of diuron degradation over a shorter 42 d timeframe (Thomas, McHugh & Waldock, 2002).
To explore the reasons for the lack of rapid heterotrophic degradation of diuron in the incubation experiment and the potential for diuron degradation in the GBR lagoon at large, we looked for phenylurea hydrolase genes, known to degrade diuron (Turnbull et al., 2001;Khurana et al., 2009), in the reads of the Dunk Island metagenome. Two reads had BLASTX similarities to a protein database covering the metal-dependent amidohydrolase superfamily. These reads were more similar to PuhB than to other PSII hydrolases and proteins from the same superfamily, with a high 51% amino acid identity over 54 amino acids (Fig. S7) indicative of the presence of potential phenylurea hydrolase homologs (Rost, 1999). Despite the short length of these metagenomic reads and the potential sequencing errors they contain, the proteins identified here may represent novel phenylurea or other PSII hydrolases, whose existence has previously been suggested (Pesce et al., 2012). Further research will be needed to characterize the sequence, structure and function of this protein and thus confirm this hypothesis. The five top scoring similarities of one of these metagenomic reads (TBLASTX, ≥98% query cover, E value ≤ 2e-22) suggest that a bacterium from the Bacteroidetes phylum (Flavobacteriia or Cytophagia order) encodes this PuhB-like protein, while the five top scoring similarities for the other read (TBLASTX, ≥97% identity, E value ≤ 5e-12) did not agree on a precise taxonomic origin. Further, calculations (see 'Materials & Methods' section) indicate that this putative hydrolase could be present in low-abundance species, in the tail of the microbial rank-abundance curve (∼0.01% relative abundance). Overall, these findings suggest that the marine microbial seed bank, ''a reservoir of dormant individuals that can potentially be resuscitated'' (Lennon & Jones, 2011), may have a potential for herbicide degradation.
Microbial enrichment studies that reported rapid biodegradation of diuron were conducted with a rich substrate or supplemented with alternative sources of carbon and nitrogen, sometimes under the form of soil or sediments (Widehem et al., 2002;Sørensen, Albers & Aamand, 2008). From this evidence, we conclude that resources may be a limiting factor for marine heterotrophs to express their diuron-degrading potential in often oligotrophic marine waters (Schaffelke et al., 2012), as is the case in marine incubations performed without supplementation (Thomas, McHugh & Waldock, 2002;Mercurio et al., 2015). GBR microorganisms are thought to metabolize nutrients from land runoff at inshore sites (Alongi & McKinnon, 2005) and, given that these sites receive high diuron and nutrient input during the wet season (Packett et al., 2009) perhaps along with diuron-degrading species, we predict that heterotrophic diuron degradation may be enhanced episodically in the GBR lagoon.

CONCLUSIONS
This study used amplicon and metagenomic sequencing to evaluate the effects of a PSII herbicide on the composition of entire prokaryotic communities, rather than selected species. It provides a baseline for future research on the impacts of herbicides on the marine ecosystem by suggesting that the effects of the PSII herbicide diuron on communities of GBR near-shore pelagic prokaryotes are limited. Metagenomic evidence suggests that prokaryotic heterotrophs in the marine water column may encode potential new herbicide hydrolase genes, though their expression may be limited by scarce environmental resources in the dry season. The apparent tolerance of marine pelagic phototrophs to diuron may have been due to the acquisition of a resistance mechanism following regular exposure to this herbicide. While no PEST sequence mutations in Cyanobacteria could explain this resistance in the present study, resident phages carried various PEST sequence types and could act as a reservoir. In summary, many components of coral reef ecosystems are stressed by herbicides from land runoff, but in contrast, the pelagic microbiome that underpins the marine food web and is integral to reef functioning, may represents an important buffer that mitigates the impacts of local anthropogenic and natural stresses on coral reefs.
• Virginia Rich and Philip Mercurio conceived and designed the experiments, performed the experiments, reviewed drafts of the paper.
• Hemerson Tonin performed the experiments, reviewed drafts of the paper.
• David G. Bourne, Andrew P. Negri and Gene W. Tyson conceived and designed the experiments, reviewed drafts of the paper.

DNA Deposition
The following information was supplied regarding the deposition of DNA sequences: PRJNA276057 SRR1819825

Data Availability
The following information was supplied regarding data availability: All raw data (DNA sequences and associated metadata) were submitted to GenBank.

Supplemental Information
Supplemental information for this article can be found online at http://dx.doi.org/10.7717/ peerj.1758#supplemental-information.