Metagenome-based diversity analyses suggest a strong locality signal for bacterial communities associated with oyster aquaculture farms in Ofunato Bay

Abstract

Ofunato Bay, in Japan, is the home of buoy-and-rope-type oyster aquaculture activities. Since the oysters filter suspended materials and excrete organic matters into the seawater, bacterial communities residing in its vicinity may show dynamic changes depending on the oyster culture activities. We employed a shotgun metagenomic technique to study bacterial communities near oyster aquaculture facilities at the center of the bay (KSt. 2) and compared the results with those of two other localities far from the station, one to the northeast (innermost bay, KSt. 1) and the other to the southwest (bay entrance, KSt. 3). Seawater samples were collected every month from January to December 2015 from the surface (1 m) and deeper (8 or 10 m) layers of the three locations, and the sequentially filtered fraction on 0.2-μm membranes was sequenced on an Illumina MiSeq system. The acquired reads were uploaded to MG-RAST for KEGG functional abundance analysis, while taxonomic analyses at the phylum and genus levels were performed using MEGAN after parsing the BLAST output. Discrimination analyses were then performed using the ROC-AUC value of the cross validation, targeting the depth (shallow or deep), locality [(KSt. 1 + KSt. 2) vs. KSt 3; (KSt. 1 + KSt. 3) vs. KSt. 2 or the (KSt. 2 + KSt. 3) vs. KSt. 1] and seasonality (12 months). The matrix discrimination analysis on the adjacent 2 continuous seasons by ROC-AUC, which was based on the datasets that originated from different depths, localities and months, showed the strongest discrimination signal on the taxonomy matrix at the phylum level for the datasets from July to August compared with those from September to June, while the KEGG matrix showed the strongest signal for the datasets from March to June compared with those from July to February. Then, the locality combination was subjected to the same ROC-AUC discrimination analysis, resulting in significant differences between KSt. 2 and KSt. 1 + KSt. 3 on the KEGG matrix. These results suggest that aquaculture activities markedly affect bacterial functions.