Community and Proteomic Analysis of Anaerobic Consortia Converting Tetramethylammonium to Methane.
Tetramethylammonium hydroxide ([[(C[H.sub.3]).sub.4]N][OH]) is a developing fluid used in photolithography processes in semiconductor, thin-film transistor liquid crystal display, and light-emitting diode manufacturing industries. The waste stream generated usually contains high concentrations of tetramethylammonium (QMA) and has been treated by the anaerobic process in serially combined with the activated sludge process  or with the processes for the nitrogen removal [2, 3] to meet the discharge standards of carbon and nitrogen. Other than the advantages such as recovery of methane energy, low-nutrient requirements and less sludge production, the anaerobic process is a key to release ammonium, which efficacy greatly influences the performance of downstream aerobic/anoxic reactors. Because of its ability to manage high concentrations of methylated compounds and ammonia, the anaerobic process as a prior treatment unit can accommodate high organic loadings to minimize the land use for the treatment plant .
Under the obligately anaerobic conditions, QMA with four methyl moieties serves as an excellent substrate for methane production through methanogenesis pathways. In the growth of methylotrophic methanogens with QMA, it has been proposed that the electrons required for reduction of the methyl groups to methane are obtained from the oxidation of an additional methyl group to C[O.sub.2] (1). During direct methylotrophic methanogenesis, the methyl groups of the compounds are transferred to the thiol group of coenzyme M (HS-CoM) by substrate-specific methyltransferases and then produce methane [5, 6]. Several marine methylotrophic methanogens under genera Methanococcoides and Methanolobus have been obtained in pure cultures to prove the direct methanogenesis from QMA [6-8]. In contrast to this mode of growth, the quaternary amines can be converted to methane in the cocultures of methanogens and sulfate-reducing bacteria . Because the reaction from QMA degradation to [H.sub.2] is energetically unfavorable ([DELTA]G[degrees]' of (2) is +319.9 kJ/reaction at the standard conditions of 1M, 25[degrees]C, 1 atm, and pH 7), it requires the hydrogen-scavenging methanogenesis (3) or homoacetogenesis (4) as a coupling reaction to proceed. In the anaerobic reactor where hydrogen partial pressure is usually sufficiently low for syntrophic growth ((2) + (3)/(4)), the corresponding bacterial and archaeal populations can interact closely in a complex manner in the sludge consortia [10, 11]. Rather, complex microbial populations are expected to participate in the indirect methanogenesis route. However, the microbial community structure to the indirect route of QMA degradation has not been reported, and to date, relevant knowledge remains relatively limited.
[(C[H.sub.3]).sub.4][N.sup.+] + 3[H.sub.2]O [right arrow] 3C[H.sub.4] + HC[O.sub.3.sup.-] + N[H.sub.4.sup.+] + [H.sup.+] ([DELTA]G[degrees]' = -86.2 kJ/reaction) (1)
[(C[H.sub.3]).sub.4][N.sup.+] + 12[H.sub.2]O [right arrow] 4HC[O.sub.3.sup.-] + 12[H.sub.2] + N[H.sub.4.sup.+] + 4[H.sup.+] ([DELTA]G[degrees]' = +319.9 kJ/reaction) (2)
HC[O.sub.3.sup.-] + 4[H.sub.2] + [H.sup.+] [right arrow] C[H.sub.4] + 3[H.sub.2O] ([DELTA]G[degrees]' = -135.6 kJ/reaction) (3)
2HC[O.sub.3.sup.-] + 4[H.sub.2] + [H.sup.+] [right arrow] C[H.sub.3]CO[O.sup.-] + 4[H.sub.2]O ([DELTA]G[degrees]' = -104.6 kJ/reaction) (4)
Free energy change ([DELTA]G[degrees]') was calculated at the standard conditions (1 M, 1 atm, 25[degrees]C, and pH 7).
Various types of anaerobic reactors, such as upflow anaerobic blanket (UASB)  and complete-mixing suspended sludge (CMSS)  reactors, have been used for the treatment of QMA-containing wastewater. The two types of reactors differ in mixing conditions (complete mixing versus plug flow) and growth models (dispersed growth versus attached growth). Because reactor performance is highly associated with microbial activity, clarifying the relevant population structure and their nexus inside anaerobic reactors can provide key information for increasing the efficiency of organic matter decomposition and methane recovery, as well as improving reactor performance. In the present study, the sludge samples from QMA-degrading CMSS and UASB reactors are analyzed using multiple PCR-based techniques and the shotgun proteomic approach , with attempt to gain insight into QMA-degrading communities in the reactor environment and the metabolic pathways used by predominant populations.
2. Materials and Methods
2.1. Sludge Samples. Anaerobic sludge from a laboratory-scale CMSS system (sample CMJP) and a full-scale UASB reactor (samples UASB1a-UASB5a) was studied. The laboratory-scale CMSS system consisting of a 10 L gas-tight continuously stirred tank reactor, and two clarifier tanks were operated at ambient temperature (23[degrees]C). The sludge from an anaerobic digester in Kitakyushu, Japan was first enriched with QMA (60 g/L) and then approximately 4000 mg COD/L of enriched sludge was introduced into the reactor. As detailed parameters in a previous study , QMA (343 mg/L) was fed as the sole substrate continuously with a stepwise decrement in hydraulic retention time (HRT) to accommodate an increasing volumetric loading of 0.6 g C/L/d. Upon sampling at day 147 of reactor operation, the system displayed excellent performance in total organic carbon (TOC) removal efficiency (>90%). The full-scale UASB reactor in Tainan, Taiwan had a working volume of approximately 1000 [m.sup.3] and has been operated for treating QMA-rich wastewater with HRT of 1.3-2 days without temperature control (18~26[degrees]C) for more than three years . The influent mainly contained QMA with the concentrations of 455~1100mg/L, which was equivalent to ~85% of total organic carbon and 98% of total organic nitrogen in wastewater. Sludge samples, UASB1aUASB5a, from the full-scale UASB reactor were taken from the sampling ports located at heights of 30 cm, 100 cm, 200 cm, 300 cm, and 400 cm (Table 1). Upon sampling, the UASB reactor achieved a TOC removal efficiency of 85-90% with slightly high pH (~7.5) in the effluent and high methane content (>90%) in the biogas stream. From each sample, 20-40 mL of sludge was preserved at -80[degrees]C in a freezer until use. The concentrations of QMA were analyzed by an ion chromatograph DX-120 (Thermo Fisher, CA, USA).
2.2. Protein Extraction and In-Sol Digestion. Proteins were recovered from the sludge following the protocol of freeze-thaw and acetone precipitation described elsewhere . The protein extracts (60 [micro]L) were first diluted with 50 mM ammonium bicarbonate and then reduced with 10 mM dithiothreitol at 25[degrees]C for 60 min, followed by cysteine-blocking with 40 mM iodoacetamide at 25[degrees]C for 30 min. Samples were digested with sequencing grade modified porcine trypsin (Promega, Madison, WI, USA) at 25[degrees]C for 16 hours. The peptides were then desalted and dried by vacuum centrifugation and store at -80[degrees]C until mass spectrometry analysis.
2.3. Online Two-Dimensional LC-MS/MS Analysis. The analysis of proteolytic peptides was performed on a nanoflow-HPLC system (Thermo Finnigan, Surveyor MS Pump Plus, Thermo Scientific) with a 2D linear ion trap mass spectrometer (LTQ-Orbitrap ELITE; Thermo Fisher, San Jose, CA, USA). Briefly, the desalted peptide mixtures were reconstituted in HPLC buffer A (30% acetonitrile/0.1% formic acid) and loaded onto a homemade column (Luna SCX 5 [micro]m, 0.5 x 100 mm). The peptides were then fractionated to 22 fractions by eluting with 0 to 100% HPLC buffer B (0.5 M ammonium chloride/30% acetonitrile/0.1% formic acid) using online 2D-HPLC (Dionex Ultimate 3000, Thermo Fisher). Each SCX fraction was diluted in-line prior to trap onto a reverse-phase column (Zorbax 300SB-C18, 0.3 x5mm; Agilent Technologies). The peptides were then separated on a HydroRP column (2.5 [micro]m, 75 [micro]m I.D. x 24 cm with a 15 fm tip) using a multistep gradient of HPLC buffer C (99.9% acetonitrile/0.1% formic acid) for 65 minutes with a flow rate of 0.25 f L/min.
The full-scan MS was performed in the Orbitrap over a range of m/z 350 to 2000 and a resolution of 60, 000 at m/z 400. The 20 data-dependent MS/MS scan events were followed by one MS scan for the 20 most abundant precursor ions in the preview MS scan. The m/z values selected for MS/MS were dynamically excluded for 40 seconds with a relative mass window of 15 ppm. The electrospray voltage was set to 2.0 kV, and the temperature of the capillary was set to 200[degrees]C. MS and MS/MS automatic gain control were set to 1000 ms (full scan) and 200 ms (MS/MS) or 2 x 106 ions (full scan) and 3 x 103 ions (MS/MS) for maximum accumulated time or ions, respectively.
2.4. Protein Identification and Annotation. The mass spectrometry data analysis was carried out using Proteome Discoverer software (version 1.4, Thermo Fisher Scientific). The MS/MS spectra were searched against the UniProt database using the Mascot search engine (Matrix Science, London, UK; version 2.5). To detail the expressed function of predominant populations, specific genome sequences were searched separately. For peptide identification, 10ppm mass tolerance was permitted for intact peptide masses, and 0.5 Da for CID fragment ions with allowance for one missed cleavage made from the trypsin digestion: oxidized methionine and acetyl (protein N-terminal) as variable modifications and carbamidomethyl (cysteine) as the fixed modifications. Peptide-spectrum matches were then filtered based on high confidence and Mascot search engine rank 1 of peptide identification to ensure an overall false discovery rate below 0.01. To ensure the protein identification, only proteins with at least two peptides identified in the LC-MS/MS analysis were included for analysis. The annotation of the identified proteins was performed with the updated Clusters of Orthologous Group (COG) (https://www.ncbi.nlm.nih.gov/COG/) and database of Kyoto Encyclopedia of Genes and Genomes (KEGG) (http://www.kegg.jp/kegg/).
2.5. DNA Recovery and Quantitative PCR. The samples were washed thrice with phosphate buffer (pH 7.4) and subjected to DNA recovery with a MoBio PowerSoil DNA isolation kit (Carlsbad, CA). The DNA quality was verified spectrophotometrically, and the concentrations of double-stranded DNA were determined with a PicoGreen quantitation reagent (Molecular Probes, Oregon, USA) in a TBS-380 Mini-Fluorometer (Turner BioSystem, CA, USA). An SYBR Green quantitative PCR (Q-PCR) was performed using a CFX96 real-time PCR detection system (BioRad, USA), with bacteria-specific primers (8F, 5'-AGA GTT TGA TCC TGG CTC AG-3'; 518R, 5'-GWA TTA CCG CGG CKG CTG-3') and archaea-specific primers
CGG GGG AGC AC-3'; 1059R, 5'-GCC ATG CAC CWC CTC T-3') [14, 15], respectively, to determine the quantities of bacterial and archaeal 16S rRNA genes in the sludge. The Q-PCR experiment was performed in triplicate, with each reaction solution (20 [micro]L) containing 200 nM of each primer, 10 [micro]L of 2X SYBR Green Supermix (BioRad), and 5-30 ng of the DNA template. The thermal program was set to 95[degrees]C for 3 min, followed by 45 cycles at 95[degrees]C for 5 s, 56[degrees]C for 30 s (for bacteria), or 54[degrees]C for 30 s (for archaea) and 72[degrees]C for 15 s; at this step of each cycle, the fluorescent products were monitored at a wavelength of 530 nm. The accuracy of the Q-PCR assay was confirmed through melting curve analysis and agarose gel electrophoresis. Calibration curves were obtained using 10-fold serial dilutions of known concentrations of cloned DNA samples that were prepared in the laboratory.
2.6. PCR Amplification, Clone Library Construction, and Phylogenetic Analysis of Archaeal 16S rRNA Genes. The archaeal communities obtained from various heights (UASB1a-UASB5a) in a methanogenic bioreactor were investigated by denaturing gradient gel electrophoresis(DGGE-) cloning approaches through amplification of the 16S rRNA by using the primers Arc344f-mod and Arch958r-mod . The amplicons were cloned to Escherichia coli DH5[alpha] cells using a pGEM-T easy vector (Promega, Madison, USA). The clones with the correct DNA inserts were screened through SP6/T7 PCR. Sequence types among the positive clones were detected by DGGE on a DCode System (BioRad). The DGGE profiles of the individual clones were compared with each other and with the DGGE profile of the entire community ; the profiles were obtained through seminested PCR amplification with the primers Arc344f-mod-GC and 524F-10-ext-rv (5'-TTA CCG CGG CTG RCA-3') . The sequences of the representative clones, referring to the phylotype, were obtained through the Sanger sequencing method and deposited in the GenBank database under the accession numbers KU569977-KU569987. The obtained sequences were compared with sequences in the GenBank database by using the BLAST search tool (http://blast.ncbi.nlm.nih.gov/Blast.cgi). A phylogenetic tree was constructed with the neighborjoining method, and bootstrap resampling analysis was performed with 1000 replicates in the MEGA6 program .
2.7. 16S rRNA Gene High-Throughput Sequencing. The compositions of archaeal and bacterial populations were analyzed using a high-throughput sequencing method. The 16S rRNA genes of archaea were amplified using a barcoded fusion archaea-specific primer set (Arc344f-mod/Arch958r-mod) . The PCR reaction mixtures comprised 10 ng of genomic DNA, 200 nM (each) of each primer, 0.2 mM dNTP (each), 1x FastStart Buffer, and 1.25 U of FastStart HiFi Polymerase (Roche, Indiana, USA) in a volume of 24 [micro]L. Reactions were conducted in a 9700 thermal cycler (Applied Biosystems, Foster, CA, USA) by denaturing DNA for 3 min at 94[degrees]C, followed by 35 cycles of 15 s at 94[degrees]C, 45 s at 52[degrees]C, and 1 min at 72[degrees]C, and a final extension at 72[degrees]C for 8 min. The PCR amplicons were purified using Agencourt AMPure XP Reagent (Beckman Coulter Inc., Beverly, Massachusetts, USA) and subsequently quantified using an Agilent Bioanalyzer (Agilent Technologies Inc., Deutschland, Waldbronn, Germany). Equimolar amounts of the PCR amplicons were mixed in a single tube, and pyrosequencing was performed using a 454 GS Junior, according to the manufacturer's instruction (Roche Applied Science, Branford, CT, USA).
For bacteria, the hypervariable region of the 16S rRNA gene was amplified using a barcoded fusion primer set comprising Pro341F and Pro805R . The amplification solution (20 [micro]L) contained 10 [micro]L of 2X Phusion HF master mix (New England BioLabs, Ipswich, MA, USA), 0.5 [micro]M of each primer (with customized barcodes present on both primers for multiplex sequencing), and 50-150ng of DNA template. The PCR conditions were as follows: initial denaturation at 95[degrees]C for 2 min, 30 cycles of annealing starting at 65[degrees]C (ending at 55[degrees]C) for 15 s, and extension at 68[degrees]C for 30 s. The amplicons in the triplicate samples were pooled and purified using an AMPureXP PCR Purification Kit (Agencourt, Brea, CA, USA) and quantified using a Qubit dsDNA HS Assay Kit on a Qubit 2.0 Fluorometer (Invitrogen, Carlsbad, CA, USA). To sequence the library preparation, Illumina adapters were attached to the amplicons by using an Illumina TruSeq DNA sample preparation kit, v2. Purified libraries were used for cluster generation and sequencing on the Illumina Miseq platform.
2.8. Analysis of High-Throughput Sequencing Data. Sequence reads were sorted on the basis of their respective barcodes into individual libraries. The sequences of primers, barcodes, and adaptors were then trimmed, and reads shorter than 150 bp and those containing ambiguous nucleotides were removed. The reverse sequences were complemented on the RDP Pipeline Initial Process. The qualified sequences were then aligned based on the RDP infernal and assigned phylogenetically using the RDP classifier at a 95% confidence level . The cluster files were subsequently employed to generate rarefaction curves, which showed an accumulating trend among OTUs defined at the 97% similarity level relative to the total number of sequence reads. To normalize the uneven sequencing effects, the OTU table was 3X randomly rarefied subset of 12, 000 and 20, 000 sequences per sample for archaeal and bacterial libraries, respectively. The Pearson correlation between specific methanogen genera was obtained using Microsoft Excel 2010. The sequencing data were submitted to the EMBL-EBI European Nucleotide Archive under the study accession number PRJEB13976 and PRJEB14151.
3.1. Microbial Diversity. The microbial communities of QMA-degrading sludge sample (CMJP) from a CMSS reactor and samples (UASB1a-UASB5a) from various heights in a full-scale UASB reactor were studied using tag-encoded amplicon sequencing of archaeal and bacterial 16S rRNA gene sequences. To assess the internal (within-sample) complexity of individual communities, various microbial diversity indices were compared. As shown in Table 1, the number of archaeal OTUs of UASB sludge displayed a cardinal distribution with the low values of 282-447 at the bottom and top of the reactor, and the highest value of up to 710 at 200 cm. According to an analysis of the Chao 1 index, species richness between the sludge samples collected at reactor heights of 100 cm and 400 cm was 1341 and 500, respectively, suggesting a 2.7-fold difference in the degree of complexity of archaeal diversity. The values of OTUs (218-299) and Chao 1 (270-384) for bacterial populations appeared to be lower than archaeal populations, even though a higher number of sequence reads was analyzed. The Shannon-Weaver index values of the samples taken at 100-300 cm were similar and were higher than the corresponding values of samples taken at heights of 30 and 400 cm. The values of Shannon-Weaver index for bacterial populations were significantly higher than archaeal populations (t-test, p = 0.01). The results showed that the archaeal populations displayed higher species richness but lower species diversity, as compared to the bacterial populations in the UASB reactor. Still, the species richness and diversity of archaeal populations and bacterial diversity in most samples of UASB reactor were higher than in the CMJP sample, which might be attributed to the influence of substrate complexity in real wastewater.
3.2. Microbial Community Structures. In the archaeal library of CMJP, more than 99.9% of the sequence reads were assigned to methanogenic archaea under Euryarchaeota and distributed with sequence abundance of 97.7, 2.2, and 0.04% in the order of Methanosarcinales, Methanobacteriales, and Methanomicrobiales, respectively. As shown in Figure 1(a), the members of the genera Methanomethylovorans (85.6%) and Methanolobus (8.3%) of the order Methanosarcinales account for 93.9% of the total reads, representing the most predominant types of methanogenic archaea in the CMSS reactor. A relatively low fraction of the sequences corresponded to acetotrophic and hydrogenotrophic methanogens, including the aceticlastic Methanosaeta (0.59%) and hydrogenotrophic methanogens such as Methanobacterium (1.93%), Methanolinea (0.01%), and Methanospirillum (0.01%). Accordingly, the results suggested that in the CMSS reactor, methane formation following the methylotrophic pathway was more critical than that following the acetotrophic and hydrogenotrophic pathways. Notably, approximately 2.7% of the total reads were assigned to the ANME-3 group, an archaeal population related phylogenetically to the order Methanosarcinales and thought to be capable of oxidizing methane to C[O.sub.2] anaerobically . Also, a low abundance of sequence reads (0.09%) was assigned to Thermoplasmata, whose members were recently recognized as methylotrophic methanogens .
The bacterial community structure in CMJP was relatively complex, as 25 bacterial phyla (reads > 0.01%) were detected. As shown in Figure 1(b), the phyla Proteobacteria, Bacteroidetes, Firmicutes, and Ignavibacteriae accounted for 31.0, 29.5, 17.1, and 16.6% of the total sequence reads, representing the top four predominant bacterial groups in the QMA-degrading consortium. Within the phyla, Gammaproteobacteria (accounting for 66.2% of Proteobacteria-related sequences) and the unclassified Bacteroidetes and Bacteroidia (accounting for 55.7 and 31.3% of Bacteroidetes-related sequences, resp.) were highly predominant. Figure 1(c) shows a detailed analysis of the sequences at the genus level. The dominant bacterial genera with a sequence abundance of >5% included Ignavibacterium (16.6%), Saccharofermentans (13.3%), Thiovirga (11.2%), Halomonas (5.4%), and the unclassified Bacteroidales group (23.8%). Members of the genus Thiovirga are chemolithoautotrophic for sulfur oxidation , whereas the other three known genera are usually organotrophic. Their presence has been associated with the anaerobic fermentation of organic substances [24-27].
As shown in Figure 2, the QMA profile and archaeal populations distributed distinctively in a full-scale UASB reactor. The decreasing trend of QMA concentrations with the reactor height revealed that the QMA was mostly degraded in the sludge bed (~1.0m) (Figure2(a)). In addition to the constantly low abundance of Methanococcales, the detected populations related to Methanosarcinales, Methanobacteriales, and Thermoplasmatales displayed a specific spatial distribution. To further elucidate the distribution, the genera data were analyzed in detail. As shown in Figure 2(b), at least three genera of methylotrophic methanogens were detected. Among them, the relative sequence abundance of Methanomethylovorans and Methanolobus in the samples remained low (<3.4%) in comparison with Methanosarcina. The relative sequence abundance of Methanosarcina varied markedly with the reactor height, being higher in sequence abundance at the bottom (74.0%) and top (89.5%) than at the middle (3.8-37.2%) of the USAB reactor. Notably, this spatial distribution of Methanosarcina appeared to exhibit an inverse relationship with aceticlastic Methanosaeta (0.08-6.6%) and hydrogenotrophic Methanobacterium (3.3-76.5%), with high Pearson's coefficients (r) of -0.961 and -0.968, respectively, but the distribution of Methanosarcina was associated positively with methylotrophic Methanomethylovorans (r = 0.992) (Figure S1). Except for Methanobacterium, all other detected hydrogenotrophic methanogens, including Methanolinea, Methanoregula, Methanoculleus, Methanospirillum, Methanobrevibacter, Methanosphaera, and Methanococcus, were present in relatively low abundance (approximately 0.77% in total). These results suggested that the archaeal populations exhibited a distinct spatial distribution in the UASB reactor, even though in which the upflow of inflow wastewater and gaseous products provided vertical mixing to the sludge.
The bacterial libraries of UASB sludge comprised totally 30 bacterial phyla (reads > 0.01%). As shown in Figure S2(a), the phyla Proteobacteria (12.4-58.6%), Bacteroidetes (15.8-27.1%), Firmicutes (15.2-19.4%), Synergistetes (2.6-8.4%), and Caldiserica (1.2-9.2%) represented the top five predominant bacterial groups in the sludge of UASB reactor. Figure S2(b) showed a detailed analysis of the abundant (>1%) sequence types at the genus level. Several sulfate-reducing bacteria such as Sulfurospirillum, and Desulfovibrio, homoacetogenic Acetobacterium, Geobacter, and Anaerovorax in the samples UASB1a (sampling height, 30 cm) and UASB5a (sampling height, 400 cm) shared similar sequence abundance higher than those in the sample UASB3a (sampling height, 200 cm). As shown in Figure S3, 11 bacterial phyla were shared by the UASB and CMSS sludge, while 14 and 9 bacterial phyla were specific to UASB and CMSS sludge, respectively. The difference of bacterial populations between the reactors might be partially attributed to the effects of influent substrates and lysed biomass.
3.3. Spatial Distribution of Archaeal Communities as Revealed by 16S rRNA Gene DGGE-Cloning. To validate the observed distinctive spatial distribution of archaeal populations, we employed the DGGE-cloning approach with the five samples from the UASB reactor. In total, 87 of 500 clones (100 clones with each sample) with 16S rRNA gene inserts with a unique banding profile were sequenced in the DGGE analysis and classified into 11 phylotypes within the phylum Euryarchaeota. As shown in Figure 3, most of these phylotypes were closely affiliated with four methanogen genera, namely, Methanosarcina (46% of total clones), Methanobacterium (27% of total clones), Methanosaeta (21% of total clones), and Methanomethylovorans (6% of total clones). The results from the DGGE-cloning analysis validated that all the three types of methanogen were predominant in the UASB reactor. The detected methanogens had a similar spatial distribution along the reactor height, as revealed by the pyrosequencing results. Figure S4 showed that Methanosarcina appeared to be more abundant (70-80% for bottom/top versus, 10-42% for middle) in the sludge samples taken at 0.3 m (bottom) and 4 m (top) of the UASB reactor, whereas the hydrogenotrophic Methanobacterium (6-7% for bottom/top versus 39-43% for middle) and the aceticlastic Methanosaeta (0-5% for bottom/top versus 15-48% for middle) populations were more abundant in the sludge samples taken at 1-3 m. As compared with other samples of UASB, the methylotrophic Methanomethylovorans was abundant at the top of the UASB reactor because 19% of total clones were detected in the UASB5a sample. Because of low number of clones screened, the DGGE-cloning analysis, however, did not detect some methanogens (e.g., Methanolobus) of low abundance, as suggested by pyrosequencing analysis. No sequence related to nonmethanogenic archaeon, Thermogymnomonas, and unclassified Thermoplasmatales was detected using DGGE-cloning analysis, but the relevant sequences could be retrieved using pyrosequencing. These dissimilarities are attributable to the biases associated with the PCR and cloning procedures, as well as sequencing depth [16, 28].
3.4. Quantitative Analysis of Archaeal and Bacterial Populations. As revealed by Q-PCR, Figure 4(a) showed the distribution of archaeal and bacterial 16S rRNA gene copies of the sludge samples in this study. For samples UASB1a, UASB5a, and CMJP, the results showed that the 16S rRNA gene copies of archaeal populations (3.43~4.54 x [10.sup.6] copies/ng-DNA for UASB1a, and UASB5a; 1.23 x [10.sup.6] copies/ng-DNA for CMJP) were present in numerically higher proportions than in bacterial populations (4.31~5.18 x [10.sup.5] copies/ng-DNA; 1.56 x [10.sup.5] copies/ng-DNA). Because the copy numbers per archaeal and bacterial genome are, respectively, 2.5 and 3.6 on averages , the corresponding cell numbers can be estimated accordingly, which leads to an archaeal to bacterial cell ratio in the range of 9.6-15.3 (Figure 4(b)). For the sludge samples from the middle of UASB reactor, the archaeal and bacterial populations shared a similar abundance with a cell ratio estimated to be 1.1-1.2, which was lower than that in the bottom and top of the UASB reactor.
3.5. Proteomic Analysis of QMA-Degrading Consortium. Since the QMA as the sole substrate was applied to the reactor, the CMJP sample was very suitable for studying microbial activities associated with the QMA degradation using the proteomic analysis. Online two-dimensional proteomic analysis identified at least 748 protein entries with the prokaryotic protein database from UniProtKB/TrEMBL. The numbers of archaea-related proteins (706 proteins, 94.4% of total identified proteins) were much higher than the bacterial counterparts (42 proteins, 5.6% of total identified proteins). Among the archaea-related proteins, 667 proteins were related to 15 genera of cultivable archaeal methanogens. As shown in Figure S5(a), M. hollandica was the most dominant population that produced the highest number (523 proteins), followed by Methanosarcina (60 proteins), Methanolobus (29 proteins), uncultured archaeon (14 proteins), and Methanococcoides (12 proteins) in order. These top 5 abundant archaeal populations accounted for protein abundance up to 93.9% of total archaeal proteins detected. The detailed analysis of the archaeal proteins showed that a large number of proteins associated with methanogenesis of methylamines and methanol were produced by the methylotrophic methanogens, including Methanosarcina, Methanolobus, Methanohalophilus, and Methanosalsum (Table S1), and M. hollandica (Table S2). As shown in Figure S5(b), the number of identified proteins assigned to every bacterial species detected was relatively low (<3 proteins per species). The overall results of proteomic analysis were in well consistent with the community analysis data (Figure 1(a)) and suggested that the methylotrophic methanogens (especially M. hollandica) exhibited vigorous activities in CMJP sludge.
3.6. Protein Expression of Methanomethylovorans hollandica. The protein expression of M. hollandica was closely elaborated by matching the mass spectrometry data with the protein sequences coded by the M. hollandica genome. Approximately, 33.3% (841 genes) of coding genes predicted in the M. hollandica genome were converted to proteins. The annotation of the identified proteins with updated COG obtained 592 COGs in 21 categories and showed that the proteins with functions assigned to the categories J (translation and ribosomal structure and biogenesis), C (Energy conservation), E (amino acid transport and metabolism), and H (coenzyme transport and metabolism) were numerically abundant (Figure S6). It was observed that the proteins involved in various methanogenesis pathways were highly expressed (Figure 5, Table S2). The expression of abundant proteins mttBC/mtbA for trimethylamine, mtbBC/ mtbA for dimethylamine, mtmBC/mtbA for monomethylamine, and mtaABC for methanol revealed the activities of M. hollandica in the transfer of methyl group from methanol and methylamines to form methyl-S-CoM, which was further reduced to methane by the methyl-CoM reductase (mcrABDG). The proteomic results showed all proteins (fwdABCDFG, ftr, mch, mtd, mer, and mtrABDEFGH) needed in the pathway of [H.sub.2]/C[O.sub.2] methanogenesis, while the acetyl-CoA synthetase (ACSS) and acetyl-CoA decarbonylase/synthase complex (cdhABCDE) that participated in the pathway of acetotrophic methanogenesis were identified, too (Figure 5). In addition, the abundant proteins matched to ROS scavengers like superoxide dismutase (L0KZ58) and catalase (L0L0M1) of Methanomethylovorans, as that archaeon may be under oxidative stress in the reactor loaded with QMA. The detection of abundant NAD(P)Hnitrite reductase (L0KW74) and hydroxylamine reductase (L0KZ19 and L0KU78) and glutamate dehydrogenase (L0KYE3) and glutamine synthetase (L0KT19) (Table S2) suggested the activities of M. hollandica in the regulation of nitrogen metabolism.
The overall results of this study revealed a central role of the Methanosarcinales in the conversion of QMA to methane in anaerobic reactors. To date, only marine methanogens have been obtained in pure cultures with QMA as the sole energy source. Methanococcoides methylutens was the first marine methanogen reported to grow with QMA [6-8]. This genus of marine methanogen, however, was not detected using PCR-based molecular methods but was detectable using a shotgun proteomic approach in the present study. The low abundance was likely due to its deteriorated growth in low-salt environments (<0.3-0.4 M sodium). Recently, Methanolobus vulcani obtained from brackish river sediment was also shown to grow with QMA in a wider range of salt concentrations (0.05-0.94 M sodium) . In this study, the Methanolobus spp. could also inhabit the lowsalt conditions, since the reactor was loaded with low sodium medium (approximately 10.4 mM).
The Methanomethylovorans populations are methylotrophic and frequently thrive in terrestrial (freshwater) ecosystems over a broad temperature range [30-35]. The proteomics data of this study suggested the methylotrophic activities of M. hollandica associated with QMA degradation (Figure 5). This result was in well accordance with our previous study  and suggested a demethylation of the QMA degradation by M. hollandica. However, the proteomic data did not uncover the protein entries for the demethylation of QMA. To search possible candidates, the amino acid sequences of proteins (MtqABC) that had been detected for the demethylation of QMA in Methanococcoides sp.  were matched to the sequence of M. hollandica genome. Three methyl-Co(III) methylamine-specific corrinoid protein: CoM methyltransferase (Metho_0037, Metho_0007, Metho_0355) of M. hollandica shared high-sequence similarity (69-74%) with the MtqA of Methanococcoides sp. but only one (mtbA, accession: L0KUG6/Metho_0007, Table S2) was detected in this study. This observation suggested that the M. hollandica used the same methyltransferase (mtbA) to mediate the formation of methyl-SCoM for tetra and tri-, di-, and mono-methylamines (Figure 5). However, likely because of high protein-substrate specificity, the M. hollandica genome in the database did not contain any homolog to the MtqBC. To identify the corresponding genes/proteins responsible for QMA demethylation, the genome sequences of the exact degrader strains are needed accordingly.
As predicted with the KEGG pathway database, the M. hollandica genome possessed the gene sets involved in the methanogenesis pathways from acetate. Interestingly, our proteomic data revealed that the M. hollandica produced abundant acetyl-CoA synthetase (ACSS) and acetyl-CoA decarbonylase/synthase complex (cdhCDE) for the acetotrophic pathway. Because the M. hollandica is not aceticlastic , one possible explanation was that a fraction of methylS-CoM molecules might be converted to acetate through a reverse acetotrophic pathway. The acetate production was known to be produced only from CO and formate by Methanosarcina acetivorans . Until recently, the homoacetogenesis (namely, acetate production from [H.sub.2]/C[O.sub.2]) was recognized in the new archaeal phylum Bathyarchaeota, in which some members were suggested to be capable of methylotrophic methanogenesis . Alternatively, the observation might account for the autotrophic carbon assimilation (Figure 5), because the CO dehydrogenase (L0KWW2), which was the key enzyme in the reductive acetyl-CoA pathway (Wood-Ljungdahl pathway) was also expressed by M. hollandica. The trait of acetate production/assimilation would be beneficial to M. hollandica to conserve more energy in the metabolism of QMA for thriving in energylimited ecosystems.
Unlike Methanomethylovorans, high species diversity of Methanosarcina was observed in both reactors. The sludge samples taken from the UASB reactor harbored at least five types of Methanosarcina 16S rRNA gene sequences. In the phylogenetic tree (Figure 3), the closest relative species, M. mazei, M. siciliae, and M. barkeri, have been reported to use methylamines for methane formation in pure cultures [38-41]. The proteomic analysis also detected the three and other three Methanosarcina species to function methylotrophically in the CMSS reactor. Our results showed that Methanomethylovorans and Methanosarcina outperformed other methylotrophic methanogens in dominating the CMSS system and the sludge bed of the UASB reactor, respectively. This suggests that the occurrence and dominance of methylotrophic methanogen populations could be influenced by the type of reactor and how it was operated. The methylotrophic Methanosarcina and Methanomethylovorans spp. compete with each other (depending on substrate utilization kinetics) and occupy the niches of high- and low-QMA concentrations, respectively . Our findings were consistent with this argument. The CMSS system with its complete mixing conferred a consistently low QMA concentration (2-5mg C/L), facilitating the dominance of the Methanomethylovorans population. In contrast to the uniformity with the CMSS reactor, the QMA profiles (48~162mg C/L) in the UASB reactor, of which the flow regime was closer to a plug-flow type, were at higher concentrations at the inlet of the reactor than at higher levels, resulting in the formation of a favorable environment for Methanosarcina in the sludge bed. However, the concentration effects cannot account for the constantly low abundance of Methanomethylovorans and Methanolobus in the UASB reactor or the dominance of Methanosarcina in the zone near the outlet (height, 400 cm). In a previous study, quantitative analysis of samples taken from a UASB reactor on different dates yielded a similar distribution of methylotrophic methanogens but with higher abundance of Methanomethylovorans than Methanosarcina . This suggested dynamic competition between Methanomethylovorans and Methanosarcina. The dominance of Methanosarcina in the zone (400 cm) close to reactor outlet might be somehow attributed to its high ability to resist oxidative stress in oxic environments .
Notably, the Methanobacterium spp. also dominated the UASB reactor, particularly in the center space (100200 cm) (approximately 76.5% of total archaeal 16S rRNA sequence), suggesting that in addition to methylotrophic methanogenesis by Methanosarcina, a major fraction of methane in that zone could be formed through a hydrogenotrophic methanogenesis pathway. Interestingly, the percentage was extraordinarily higher than that of the hydrogenotrophic methanogens detected in the CMSS reactor (2.0%) and in the bottom and top of UASB reactor, as well as other anaerobic reactors degrading industrial wastewaters and sludge biomass (<45.6%) [35, 44]. This unusual was also related to the contrast of archaeal to bacterial cell ratios analyzed by using Q-PCR. However, the absolute abundance of Methanobacterium still remained low in the center of the UASB reactor (Figure 4). Reasonably, the [H.sub.2] source to the growth of hydrogenotrophic Methanobacterium spp. could be attributed to the bacterial degradation of decayed sludge biomass and fermentable substrates in wastewater (~15% of TOC) .
Besides, the finding might be indicative of the methyl oxidation of methylamines (2) to [H.sub.2] and C[O.sub.2], which were converted to methane (3). Finke and coworkers (2007) proposed that the methylotrophic Methanosarcina can produce [H.sub.2] through shifting the metabolism to oxidize more methyl groups of methylated substrates to C[O.sub.2]. The production and diffusive loss of [H.sub.2] from the Methanosarcina cell could be achieved by growing with the hydrogen-scavenging methanogens to facilitate low [H.sub.2] partial pressure conditions . Consequently, the observed distribution of methanogens suggested a mutual interaction between methylotrophic Methanosarcina and hydrogen-scavenging Methanobacterium for the [H.sub.2] transfer reaction to occur. Such interaction was not recognized in the CMSS reactor (dispersed cell growth) but was specific to the UASB reactor (attached cell growth) in this study and in the previous studies with marine sediment , methanol/acetate-degrading cocultures , and methanethiol-degrading granular sludge . Considering absolute abundance of Methanobacterium and assumed conditions inside the UASB reactor (QMA= 1 mM, bicarbonate = 0.1 M, ammonium = 0.1 M, pH = 7, [H.sub.2] partial pressure = [10.sup.-5] atm, and methane = 0.75 atm), the QMA degradation (2) became thermodynamically favorable ([DELTA]G' = -33.9 kJ/ reaction). Since the degree of diverting the methyl oxidation into [H.sub.2] production was dependent on [H.sub.2] partial pressure , it is speculated that the syntrophic degradation of methylated compounds by Methanosarcina spp. could be facilitated at low [H.sub.2] partial pressure conditions. The Methanosarcina spp. may be superior to M. hollandica to produce extracellular [H.sub.2] during metabolizing methylated compounds, because the Methanosarcina genomes possessed the gene encoding membrane-bound hydrogenase [30, 47], but M. hollandica lacked the kind of genes in the genome. Methanosarcina spp. have the potential to use the hydrogenotrophic pathway for methane production. However, whether simultaneous methylotrophic and hydrogenotrophic activities could occur in the same organism remains unclear, and further study is required to confirm this.
Aceticlastic Methanosaeta with low-medium sequence abundance (4.7-6.6%, Figure 2(b)) inhabited the zone (1-2 m height in UASB reactor) where methane was largely produced from [H.sub.2]/C[O.sub.2]. The growth of Methanosaeta in the CMSS and UASB reactors might not have been inhibited by ammonia toxicity because the ammonia concentrations (approximately 15-54 mM) resulting from the degradation of QMA were markedly lower than the reported inhibitory levels (>75 mM) [48, 49]. In the UASB reactor, the acetate sources could be partially attributable to fermentation of the decay of sludge biomass and substrates in sanitary wastewater, as well as the potential acetogenesis of methylotrophic methanogens. It was also likely the contribution from the bacterial homoacetogenesis. This was further supported by the results of the high-throughput analysis of bacterial 16S rRNA genes, which showed detectable abundance of homoacetogens such as Acetobacterium (1.0-4.4%) and Clostridium (0.05-1.48%) (Figure S2(b)). These observations may in part account for higher relative abundance of Methanosaeta in the UASB reactor (~6.61%) than the CMSS reactor (0.59%). Since the activity of homoacetogenic populations in the conversion of methylated compounds required close association with hydrogenotrophic methanogens , the nexus of actual microbial interactions and metabolisms in the UASB reactor was far complicated and should be further studied in the future.
In the present study, the microbial community structures of methanogenic sludge samples obtained from QMA-degrading CMSS and UASB reactors were contrasted. The overall results showed high degree of diversity of freshwater methylotrophic methanogens to sustain the anaerobic degradation of QMA to methane. In particular, the methylotrophic Methanomethylovorans sp. highly dominated in the CMSS reactor, whereas diverse methylotrophic Methanosarcina spp. possibly in the association with hydrogenotrophic Methanobacterium distributed abundantly in the UASB reactor, which suggested the effects of the reactor configuration and operation on the QMA-degrading communities. This finding advances the understanding of methanogenic degradation of quaternary amines in the engineering environments and may facilitate improving reactor technology used in the anaerobic treatment of relevant wastewater.
Please Note: Illustration(s) are not available due to copyright restrictions.
Conflicts of Interest
The authors declare that they have no conflicts of interest.
This work was supported by SAS-MOST Joint Research Cooperation (Bilateral project Slovakia-Taiwan) under Grant nos. MOST 104-2923-E-006-001-MY3 and SASMOST JRP 2014/3. The authors thank Ms. Chu G.-L. and Mr. Kuo H.-C. for experimental assistance.
Table S1: annotation of proteins related to the conversion of methylamines by the methanogens in the CMJP sample analyzed in this study. Table S2: annotation of proteins related to the conversion of methylamines, methanogenesis, and nitrogen metabolism in Methanomethylovorans hollandica detected in the CMJP sample in this study. Figure S1: (a) correlation of Methanosarcina versus Methanosaeta and Methanomethylovorans. (b) Correlation of Methanosarcina versus Methanobacterium. The lowercase "r" denotes the Pearson coefficient. Figure S2: high-throughput sequencing analysis of bacterial 16S rRNA sequences of sludge samples taken from different heights (30 cm, 200 cm, and 400 cm) of a full-scale UASB reactor treating tetramethylammoniumrich wastewater. Figure S3: Venn diagram analysis of bacterial phyla detected with CMJP and UASB samples. Figure S4: quantitative distribution of archaeal populations in the full-scale UASB reactor for treating tetramethylammoniumrich wastewater as revealed by a DGGE-cloning approach. Figure S5: numerical distribution and accumulated relative abundance of the identified proteins from (a) archaeal and (b) bacterial members in the sludge sampled from a tetramethylammonium-degrading CMSS reactor. Figure S6: distribution of the identified proteins assigned to the functional categories based on cluster of orthologous group (COG) classification. (Supplementary Materials)
 T. H. Hu, L. M. Whang, P. W. Liu et al., "Biological treatment of TMAH (tetra-methyl ammonium hydroxide) in a full-scale TFT-LCD wastewater treatment plant," Bioresource Technology, vol. 113, pp. 303-310, 2012.
 H. L. Lin, S. K. Chen, Y. W. Huang, W. C. Chen, W. C. Chien, and S. S. Cheng, "A combined upflow anaerobic sludge bed, aerobic, and anoxic fixed-bed reactor system for removing tetramethylammonium hydroxide and nitrogen from light-emitting diode wastewater," Environmental Technology, vol. 37, no. 11, pp. 1325-1333, 2016.
 H.-L. Lin, B.-K. Chen, H.-P. Hsia et al., "Use of two-stage biological process in treating thin film transistor liquid crystal display wastewater of tetramethylammonium hydroxide, " Sustainable Environment Research, vol. 21, no. 3, pp. 155-160, 2011.
 K.-F. Chang, S.-Y. Yang, H.-S. You, and J. R. Pan, "Anaerobic treatment of tetramethylammonium hydroxide (TMAH) containing wastewater," IEEE Transactions on Semiconductor Manufacturing, vol. 21, no. 3, pp. 486-491, 2008.
 J. A. Krzycki, "Function of genetically encoded pyrrolysine in corrinoid-dependent methylamine methyltransferases, " Current Opinion in Chemical Biology, vol. 8, no. 5, pp. 484-491, 2004.
 T. Ticak, D. Hariraju, M. B. Arcelay, B. A. Arivett, S. E. Fiester, and D. J. Ferguson Jr., "Isolation and characterization of a tetramethylammonium-degrading Methanococcoides strain and a novel glycine betaine-utilizing Methanolobus strain, " Archives of Microbiology, vol. 197, no. 2, pp. 197-209, 2015.
 S. Asakawa, K. Sauer, W. Liesack, and R. K. Thauer, "Tetramethylammonium: coenzyme M methyltransferase system from Methanococcoides sp," Archives of Microbiology, vol. 170, no. 4, pp. 220-226, 1998.
 K. Tanaka, "Anaerobic degradation of tetramethylammonium by a newly isolated marine methanogen," Journal of Fermentation and Bioengineering, vol. 78, no. 5, pp. 386-388, 1994.
 G. M. King, "Utilization of hydrogen, acetate, and "noncompetitive"; substrates by methanogenic bacteria in marine sediments," Geomicrobiology Journal, vol. 3, no. 4, pp. 275-306, 1984.
 J. Weijma and A. J. Stams, "Methanol conversion in high-rate anaerobic reactors," Water Science and Technology, vol. 44, no. 8, pp. 7-14, 2001.
 N. Finke, T. M. Hoehler, and B. B. Jorgensen, "Hydrogen 'leakage' during methanogenesis from methanol and methylamine: implications for anaerobic carbon degradation pathways in aquatic sediments," Environmental Microbiology, vol. 9, no. 4, pp. 1060-1071, 2007.
 B. Liu, K. Yoshinaga, J.-H. Wu et al., "Kinetic analysis of biological degradation for tetramethylammonium hydroxide (TMAH) in the anaerobic activated sludge system at ambient temperature," Biochemical Engineering Journal, vol. 114, pp. 42-49, 2016.
 H.-J. Huang, W.-Y. Chen, and J.-H. Wu, "Total protein extraction for metaproteomics analysis of methane producing biofilm: the effects of detergents," International Journal of Molecular Sciences, vol. 15, no. 6, pp. 10169-10184, 2014.
 M. Haakensen, C. Dobson, H. Deneer, and B. Ziola, "Real-time PCR detection of bacteria belonging to the Firmicutes phylum, " International Journal of Food Microbiology, vol. 125, no. 3, pp. 236-241, 2008.
 Y. Ji, H. Kim, H. Park et al., "Modulation of the murine microbiome with a concomitant anti-obesity effect by Lactobacillus rhamnosus GG and Lactobacillus sakei NR28," Beneficial Microbes, vol. 3, no. 1, pp. 13-22, 2012.
 A. C. Pires, D. F. Cleary, A. Almeida et al., "Denaturing gradient gel electrophoresis and barcoded pyrosequencing reveal unprecedented archaeal diversity in mangrove sediment and rhizosphere samples," Applied and Environmental Microbiology, vol. 78, no. 16, pp. 5520-5528, 2012.
 A. Puskarova, M. Buckova, B. Habalova, L. Krakova, A. Makova, and D. Pangallo, "Microbial communities affecting albumen photography heritage: a methodological survey, " Scientific Reports, vol. 6, no. 1, article 20810, 2016.
 K. Tamura, G. Stecher, D. Peterson, A. Filipski, and S. Kumar, "MEGA6: molecular evolutionary genetics analysis version 6.0," Molecular Biology and Evolution, vol. 30, no. 12, pp. 2725-2729, 2013.
 S. Takahashi, J. Tomita, K. Nishioka, T. Hisada, and M. Nishijima, "Development of a prokaryotic universal primer for simultaneous analysis of bacteria and archaea using next-generation sequencing," PLoS One, vol. 9, no. 8, article e105592, 2014.
 Q. Wang, G. M. Garrity, J. M. Tiedje, and J. R. Cole, "Naive Bayesian classifier for rapid assignment of rRNA sequences into the new bacterial taxonomy," Applied and Environmental Microbiology, vol. 73, no. 16, pp. 5261-5267, 2007.
 S. Scheller, M. Goenrich, R. Boecher, R. K. Thauer, and B. Jaun, "The key nickel enzyme of methanogenesis catalyses the anaerobic oxidation of methane," Nature, vol. 465, no. 7298, pp. 606-608, 2010.
 M. Poulsen, C. Schwab, B. B. Jensen et al., "Methylotrophic methanogenic Thermoplasmata implicated in reduced methane emissions from bovine rumen," Nature Communications, vol. 4, p. 1428, 2013.
 T. Ito, K. Sugita, I. Yumoto, Y. Nodasaka, and S. Okabe, "Thiovirga sulfuroxydans gen. nov., sp. nov., a chemolithoautotrophic sulfur-oxidizing bacterium isolated from a microaerobic waste-water biofilm," International Journal of Systematic and Evolutionary Microbiology, vol. 55, no. 3, pp. 1059-1064, 2005.
 S. Chen, L. Niu, and Y. Zhang, "Saccharofermentans acetigenes gen. nov., sp. nov., an anaerobic bacterium isolated from sludge treating brewery wastewater," International Journal of Systematic and Evolutionary Microbiology, vol. 60, no. 12, pp. 2735-2738, 2010.
 J. A. Mata, J. Martinez-Canovas, E. Quesada, and V. Bejar, "A detailed phenotypic characterisation of the type strains of Halomonas species," Systematic and Applied Microbiology, vol. 25, no. 3, pp. 360-375, 2002.
 T. Iino, K. Mori, Y. Uchino, T. Nakagawa, S. Harayama, and K. Suzuki, "Ignavibacterium album gen. nov., sp. nov., a moderately thermophilic anaerobic bacterium isolated from microbial mats at a terrestrial hot spring and proposal of Ignavibacteria classis nov., for a novel lineage at the periphery of green sulfur bacteria," International Journal of Systematic and Evolutionary Microbiology, vol. 60, no. 6, pp. 1376-1382, 2010.
 Z. Liu, N. U. Frigaard, K. Vogl et al., "Complete genome of Ignavibacterium album, a metabolically versatile, flagellated, facultative anaerobe from the phylum Chlorobi," Frontiers in Microbiology, vol. 3, p. 185, 2012.
 L. J. Forney, X. Zhou, and C. J. Brown, "Molecular microbial ecology: land of the one-eyed king," Current Opinion in Microbiology, vol. 7, no. 3, pp. 210-220, 2004.
 S. F. Stoddard, B. J. Smith, R. Hein, B. R. Roller, and T. M. Schmidt, "rrnDB: improved tools for interpreting rRNA gene abundance in bacteria and archaea and a new foundation for future development," Nucleic Acids Research, vol. 43, no. D1, pp. D593-D598, 2015.
 Y. Liu and W. B. Whitman, "Metabolic, phylogenetic, and ecological diversity of the methanogenic archaea," Annals of the New York Academy of Sciences, vol. 1125, no. 1, pp. 171-189, 2008.
 B. P. Lomans, R. Maas, R. Luderer et al., "Isolation and characterization of Methanomethylovorans hollandica gen. nov., sp. nov., isolated from freshwater sediment, a methylotrophic methanogen able to grow on dimethyl sulfide and methanethiol, " Applied and Environmental Microbiology, vol. 65, no. 8, pp. 3641-3650, 1999.
 K. Roest, M. Altinbas, P. L. Paulo et al., "Enrichment and detection of microorganisms involved in direct and indirect methanogenesis from methanol in an anaerobic thermophilic bioreactor," Microbial Ecology, vol. 50, no. 3, pp. 440-446, 2005.
 B. Jiang, S. N. Parshina, W. van Doesburg, B. P. Lomans, and A. J. Stams, "Methanomethylovorans thermophila sp. nov., a thermophilic, methylotrophic methanogen from an anaerobic reactor fed with methanol," International Journal of Systematic and Evolutionary Microbiology, vol. 55, no. 6, pp. 2465-2470, 2005.
 F. A. de Bok, R. C. van Leerdam, B. P. Lomans et al., "Degradation of methanethiol by methylotrophic methanogenic archaea in a lab-scale upflow anaerobic sludge blanket reactor, " Applied and Environmental Microbiology, vol. 72, no. 12, pp. 7540-7547, 2006.
 J. H. Wu, H. P. Chuang, M. H. Hsu, and W. Y. Chen, "Use of a hierarchical oligonucleotide primer extension approach for multiplexed relative abundance analysis of methanogens in anaerobic digestion systems," Applied and Environmental Microbiology, vol. 79, no. 24, pp. 7598-7609, 2013.
 M. Rother and W. W. Metcalf, "Anaerobic growth of Methanosarcina acetivorans C2A on carbon monoxide: an unusual way of life for a methanogenic archaeon," Proceedings of the National Academy of Sciences of the United States of America, vol. 101, no. 48, pp. 16929-16934, 2004.
 Y. He, M. Li, V. Perumal et al., "Genomic and enzymatic evidence for acetogenesis among multiple lineages of the archaeal phylum Bathyarchaeota widespread in marine sediments, " Nature Microbiology, vol. 1, no. 6, article 16035, 2016.
 H. Hippe, D. Caspari, K. Fiebig, and G. Gottschalk, "Utilization of trimethylamine and other N-methyl compounds for growth and methane formation by Methanosarcina barkeri, " Proceedings of the National Academy of Sciences of the United States of America, vol. 76, no. 1, pp. 494-498, 1979.
 K. R. Sowers, S. F. Baron, and J. G. Ferry, "Methanosarcina acetivorans sp. nov., an acetotrophic methane-producing bacterium isolated from marine sediments," Applied and Environmental Microbiology, vol. 47, no. 5, pp. 971-978, 1984.
 S. Ni, C. R. Woese, H. C. Aldrich, and D. R. Boone, "Transfer of Methanolobus siciliae to the genus Methanosarcina, naming it Methanosarcina siciliae, and emendation of the genus Methanosarcina," International Journal of Systematic Bacteriology, vol. 44, no. 2, pp. 357-359, 1994.
 K. Veit, C. Ehlers, and R. A. Schmitz, "Effects of nitrogen and carbon sources on transcription of soluble methyltransferases in Methanosarcina mazei strain Go1," Journal of Bacteriology, vol. 187, no. 17, pp. 6147-6154, 2005.
 L. M. Whang, T. H. Hu, P. W. Liu et al., "Molecular analysis of methanogens involved in methanogenic degradation of tetra-methylammonium hydroxide in full-scale bioreactors, " Applied Microbiology and Biotechnology, vol. 99, no. 3, pp. 1485-1497, 2015.
 R. Angel, D. Matthies, and R. Conrad, "Activation of methanogenesis in arid biological soil crusts despite the presence of oxygen," PLoS One, vol. 6, no. 5, article e20453, 2011.
 T. Narihiro, T. Terada, A. Ohashi et al., "Quantitative detection of culturable methanogenic archaea abundance in anaerobic treatment systems using the sequence-specific rRNA cleavage method," The ISME Journal, vol. 3, no. 5, pp. 522-535, 2009.
 G. M. King, M. J. Klug, and D. R. Lovley, "Metabolism of acetate, methanol, and methylated amines in intertidal sediments of lowes cove, maine," Applied and Environmental Microbiology, vol. 45, no. 6, pp. 1848-1853, 1983.
 T. J. Phelps, R. Conrad, and J. G. Zeikus, "Sulfate-dependent interspecies [H.sub.2] transfer between Methanosarcina barkeri and Desulfovibrio vulgaris during coculture metabolism of acetate or methanol," Applied and Environmental Microbiology, vol. 50, no. 3, pp. 589-594, 1985.
 R. K. Thauer, A. K. Kaster, H. Seedorf, W. Buckel, and R. Hedderich, "Methanogenic archaea: ecologically relevant differences in energy conservation," Nature Reviews. Microbiology, vol. 6, no. 8, pp. 579-591, 2008.
 M. Westerholm, L. Leven, and A. Schnurer, "Bioaugmentation of syntrophic acetate-oxidizing culture in biogas reactors exposed to increasing levels of ammonia," Applied and Environmental Microbiology, vol. 78, no. 21, pp. 7619-7625, 2012.
 J. J. Werner, M. L. Garcia, S. D. Perkins et al., "Microbial community dynamics and stability during an ammonia-induced shift to syntrophic acetate oxidation," Applied and Environmental Microbiology, vol. 80, no. 11, pp. 3375-3383, 2014.
Wei-Yu Chen, (1) Lucia Krakova, (2) Jer-Horng Wu, (1) Domenico Pangallo, (2) Lenka Jeszeova, (2) Bing Liu, (3) and Hidenari Yasui (3)
(1) Department of Environmental Engineering, National Cheng Kung University, No. 1, University Road, East District, Tainan City 701, Taiwan
(2) Institute of Molecular Biology, Slovak Academy of Sciences, Dubravska Cesta 21, 84551 Bratislava, Slovakia
(3) Faculty of Environmental Engineering, The University of Kitakyushu, 1-1 Hibikino, Wakamatsu, Kitakyushu 808-0135, Japan
Correspondence should be addressed to Jer-Horng Wu; email@example.com
Received 23 June 2017; Revised 25 September 2017; Accepted 31 October 2017; Published 17 December 2017
Academic Editor: William B. Whitman
Caption: Figure 1: High-throughput sequencing analysis of archaeal and bacterial 16S rRNA sequences of sludge samples taken from a tetramethylammonium-degrading CMSS reactor. (a) Phylogenetic distribution of the detected archaeal populations. "Others" included the low abundant populations, namely, Methanolinea (0.01%), Methanospirillum (0.01%), Methanobacteriaceae (unclassified) (0.04%), and Thermoplasmata (0.09%). (b) Phylogenetic distribution of the detected bacterial populations at the phylum and class levels. "Others" included 10 bacterial phylum taxa with each abundance <0.05%. (c) Distribution of bacterial genera with relative sequence abundance >1%.
Caption: Figure 2: Tetramethylammonium (QMA) profile (a) and quantitative distribution of archaeal populations in the full-scale UASB reactor for treating tetramethylammonium-rich wastewater as revealed by the pyrosequencing 16S rRNA gene amplicons (sequence abundance >0.8%) (b). "Other [H.sub.2] methanogens" included the low abundant populations, namely, Methanolinea (<0.07%), Methanoculleus (<0.01%), Methanospirillum (<0.03%), Methanobrevibacter (<0.03%), Methanosphaera (<0.01%), and Methanococcus (<0.01%).
Caption: Figure 3: Phylogenetic tree and distribution of archaeal 16S rRNA genes retrieved from the sludge taken at various heights (0.3-4 m) of the full-scale UASB reactor for treating tetramethylammonium-rich wastewater. The 16S rRNA gene sequences were obtained through the DGGE-cloning approach, and the representative clones were shown in bold in the tree. The values in the square brackets indicated the numbers of clones that were similar to the representative phylotypes (>99% sequence similarity). The 16S rRNA gene sequence of Methanopyrus kandleri (NC003551) was used as an outgroup to root the tree. Bootstrap values of 50-75% and >75% were indicated with black circles and squares, respectively, at the nodes of the tree. The scale bar represented the estimated number of nucleotide changes per sequence position.
Caption: Figure 4: (a) 16S rRNA gene copies (on average) of archaeal and bacterial populations and (b) the estimated ratios of archaeal to bacterial populations in methanogenic sludge from tetramethylammonium-degrading CMSS (sample CMJP) and UASB reactors (samples UASB1a-UASB5a). Coefficients of variance (CV) for archaeal and bacterial Q-PCR experiments (triplicate analysis per sample) are 1.4~35.1% and 2.3~12.8%, respectively.
Caption: Figure 5: Proteins expressed in the pathways of methylotrophic, hydrogenotrophic, and acetotrophic methanogenesis of Methanomethylovorans hollandica in the conversion of tetramethylammonium (QMA) to methane. TMA: trimethylamine; DMA: dimethylamine; MMA: monomethylamine; MeOH: methanol; MFR: methanofuran; [H.sub.4]MPT: tetrahydromethanopterin; HS-CoM coenzyme M; HS-CoB: coenzyme B; CoM-S-S-CoB: heterodisulfide of HS-CoM and HS-CoB; [F.sub.420][H.sub.2]: reduced coenzyme [F.sub.420]; [Fd.sub.red]: reduced ferredoxin; [Fd.sub.ox]: oxidized ferredoxin; [H.sub.2]ase: hydrogenase. The proteins detected in the study were highlighted in red and referred to in Table S2.
Table 1: Reads and diversity indices obtained from the high-throughput analysis of archaeal and bacterial 16S rRNA gene amplicons for sludge samples taken from tetramethylammonium-degrading CMSS and UASB reactors. Sample ID Reactor Target Height (cm) Reads (1) UASB1a UASB Archaea 30 12,000 UASB2a UASB Archaea 100 12,000 UASB3a UASB Archaea 200 12,000 UASB4a UASB Archaea 300 12,000 UASB5a UASB Archaea 400 12,000 UASB1a UASB Bacteria 30 20,000 UASB3a UASB Bacteria 200 20,000 UASB5a UASB Bacteria 400 20,000 CMJP CMSS Archaea NA 12,000 CMJP CMSS Bacteria NA 20,000 Sample ID Richness (2) Chao1 (3) Shannon (4) UASB1a 447 839 3.30 UASB2a 703 1341 3.94 UASB3a 710 1174 4.13 UASB4a 641 1196 4.09 UASB5a 282 500 2.04 UASB1a 218 270 5.02 UASB3a 249 266 5.88 UASB5a 299 384 4.37 CMJP 250 367 2.75 CMJP 232 375 2.54 (1) Number of qualified 16S rRNA sequence reads. (2) Number of detected OTUs (97% similarity). (3) Chao1 richness estimate; a higher number represents a higher degree of diversity. (4) Shannon-Weiner index; a higher number represents a higher degree of diversity. NA: not applicable.
|Printer friendly Cite/link Email Feedback|
|Title Annotation:||Research Article|
|Author:||Chen, Wei-Yu; Krakova, Lucia; Wu, Jer-Horng; Pangallo, Domenico; Jeszeova, Lenka; Liu, Bing; Yasui,|
|Date:||Jan 1, 2018|
|Previous Article:||Enhancement of Organic Matter Removal in an Integrated Biofilm-Membrane Bioreactor Treating High-Salinity Wastewater.|
|Next Article:||Diversity and Niche of Archaea in Bioremediation.|