16S rRNA Gene Amplicon Profiling of Anaerobic Bulking-Associated Prokaryotic Microbiota in a Mesophilic Expanded Granular Sludge Bed Reactor for Beverage Wastewater Treatment

Information regarding prokaryotic microbiota associated with anaerobic bulking is limited. Here, we provide 16S rRNA gene-based prokaryotic diversity profiles for anaerobic bulking and healthy granular sludge in a mesophilic expanded granular sludge bed (EGSB) reactor. These data were tabulated at the phylum level based on high-quality reads.

U pflow anaerobic sludge blanket (UASB) and expanded granular sludge bed (EGSB) reactors are core anaerobic wastewater treatment reactors for medium-and high-strength organic wastewater (1). Healthy granular sludge formation inside these reactors is essential. However, changes in granular sludge sedimentation due to the overgrowth of specific filamentous microorganisms contribute to anaerobic bulking, leading to the outflow of granular sludge from UASB and EGSB reactors (2)(3)(4). Anaerobic bulking, identified during the startup of a full-scale mesophilic (30 to 32°C) EGSB reactor (established in Shizuoka, Japan, in February 2014) for beverage wastewater treatment, was caused by certain filamentous microorganisms. 16S rRNA gene amplicon profiles of the microbiota were generated to determine biotechnological applications for reactor recovery and anaerobic bulking prevention.
Bulking sludge and two healthy granular sludges (500 ml each) were collected from a sampling port located 2 m from the EGSB reactor bottom (reactor volume, 177 m 3 ) ( Table 1). The composition of wastewater supplied to the EGSB reactor, including chemical oxygen demand (COD), was measured as reported previously (3,5,6). Most of the COD content in the wastewater comprised carbohydrates, proteins, and some organic acids, 2,880, 457, and 922 mg COD/liter, respectively. DNA was extracted from the sludges as described previously (5). PCR amplification of the V4 region of the 16S rRNA gene was performed using Blend Taq polymerase (Toyobo, Osaka, Japan) and a 515F/806R primer set (7). PCR product sequences were analyzed using the MiSeq platform and MiSeq reagent kit v2 (2 ϫ 300 bp; Illumina, San Diego, CA, USA) at the Bioengineering Lab. Co., Ltd. (Kanagawa, Japan) with six to seven technical replicates. The adaptor, index, and primer regions of each raw sequence read were trimmed using the FASTX-Toolkit v0.0.13 (8). Read sequences of Յ40 bp and with ambiguous bases and low-quality sequences (quality score, ՅQ20) were filtered out, together with their paired-end reads, using Sickle v1.33 (9). High-quality paired-end reads were merged using PEAR v0.9.10 with default settings (10). Merged sequences of Յ245 and Ն260 bp were discarded using SeqKit v0.8.0 (11). Operational taxonomic units (OTUs) of the sludge microbiota were classified using QIIME v1.9.1 (12) and the SILVA database (release 132) with 97% identity (13). The representative OTUs (Ͼ0.1%) at the phylum level within the domains Bacteria and Archaea are summarized in Table 1.
Data availability. The 16S rRNA gene amplicon data set was deposited in the DDBJ Sequence Read Archive (SRA) under accession number DRP005108 and SRA run accession numbers DRR180014 to DRR180032.

ACKNOWLEDGMENTS
This work was financially supported by JSPS KAKENHI grant 17H03333 (Grants-in-Aid for Scientific Research [B]) and JGC-S Scholarship Foundation Research grant program no. 1637.
We thank Masako Hamada for technical support at Toyohashi University of Technology.