Science of the Total Environment 652 (2019) 709–717

Contents lists available at ScienceDirect

Science of the Total Environment

journal homepage: www.elsevier.com/locate/scitotenv

Characterization of microbial community and main functional groups of

prokaryotes in thermophilic anaerobic co-digestion of food waste and
paper waste
Lu Li, Yu Qin, Zhe Kong, Jing Wu, Kengo Kubota, Yu-You Li ⁎
Department of Civil and Environmental Engineering, Graduate School of Engineering, Tohoku University, 6-6-06 Aza-Aoba, Aramaki, Aoba Ward, Sendai, Miyagi 980-8579, Japan

H I G H L I G H T S G R A P H I C A L A B S T R A C T

• Effect of paper waste on microbial com-

munity in thermophilic digestion was
clarified.
• Functional groups of prokaryotes
changed greatly as the paper waste
ratio increased.
• Community of cellulose-degrading bac-
teria changed greatly due to declined
pH value.
• Hydrogenotrophic methanogenesis was
enhanced when the paper waste ratio
≥ 50%.
• The predominant microorganism re-
sponse for each metabolic step was
summarized.

a r t i c l e i n f o a b s t r a c t

Article history: The thermophilic anaerobic co-digestion of food waste and paper waste was successfully operated with a 0% to
Received 19 August 2018 70% fraction of paper waste. The variation of functional microbial community was investigated by 16S rRNA
Received in revised form 17 October 2018 gene analysis. The results indicated that the hydrolyzing bacterial community changed from carbohydrate/
Accepted 21 October 2018
protein-degrading bacteria to cellulose-degrading bacteria when the paper waste ratio was higher than 50%. Sig-
Available online 22 October 2018
nificant changes in the taxon responsible for cellulose degradation were found depending on the paper waste
Editor: Jay Gan fraction. Cellulose-degrading bacteria outcompeted lactic acid bacteria in the degradation of monosaccharide,
resulting in a decline in the proportion of lactic acid bacteria and the absence of an accumulation of lactic acid.
Keywords: At high paper waste ratios, because the cellulose-degrading bacteria, such as Defluviitoga tunisiensis, were
Thermophilic anaerobic co-digestion more likely to degrade monosaccharides directly to acetate and hydrogen rather than to propionate and butyrate,
Food waste the abundance of syntrophs was reduced. The variation of those bacteria with high H2-producing ability signifi-
Paper waste cantly influenced the proportion of hydrogenotrophic archaea. The change in the microbial community as the
Microbial community paper waste fraction increased was illustrated with regard to anaerobic degradation steps.
Cellulose degradation
© 2018 Elsevier B.V. All rights reserved.

1. Introduction

Thanks to its high loading tolerance, energy recovery in the form of

methane and its low generation of biomass, anaerobic digestion has
⁎ Corresponding author. been widely applied in the treatment of the organic fraction of munici-
E-mail address: [email protected] (Y.-Y. Li). pal solid waste (Banks et al., 2011). As a mono-substrate, the digestion

https://doi.org/10.1016/j.scitotenv.2018.10.292
0048-9697/© 2018 Elsevier B.V. All rights reserved.
710 L. Li et al. / Science of the Total Environment 652 (2019) 709–717

of paper waste (PW), one of the biggest components of municipal solid metabolic pathway for the co-digestion of FW and PW in the thermo-
waste, has been shown to be difficult (Jokela et al., 1997; Kamali et al., philic anaerobic digestion system was proposed.
2016). While PW contains a large amount of carbohydrate, it lacks nutri-
ents: an excessive dosage of PW in the anaerobic digestion system in-
creases the C/N ratio. Because of the lack of basic raw materials for the 2. Materials and methods
synthesis of amino acids and proteins (McCarty, 1964), this increase in
the PW ratio significantly inhibits the assimilation of microorganisms. 2.1. Experimental apparatus and preparation of feedstock
On the other hand, the carbohydrate in PW is in the form of cellulose
and hemicellulose: unlike starch or glucose, these are referred to as re- A thermophilic anaerobic digestion system was operated for the
fractory carbohydrates due to their high degree of polymerization co-digestion of FW and PW. As illustrated in Fig. 1, the system con-
(Leschine, 1995). Besides, lignin, another major component in PW, has tains a substrate tank in which the substrate was stored at 4 °C,
also been classified as a refractory organic matter for degradation. and a continuously stirred tank reactor with an effective volume of
These factors make it difficult for PW to be degraded anaerobically 3 L. The reactor was operated under the thermophilic condition of
and to be consumed by microorganisms. 55 °C, an organic loading rate of approximately 14.0 g COD L−1 d−1
Conversely, food waste (FW) is classified as a nitrogen-rich biomass and a hydraulic retention time of 30 days. Biogas was collected and
due to its high protein and amino acid content (Banks et al., 2011). It is measured by a wet tip gas meter (W-NK-0.5B, Shinagawa, Japan).
therefore reasonable to assume that a supply of essential nitrogen ele- Water bath circulations were introduced into the system to main-
ment from FW has potential to foster higher effectiveness in the anaer- tain the temperature for the substrate tank and the reactor. Seed
obic co-digestion of PW and FW, with a higher harvest than the sole sludge was collected from a local wastewater treatment plant
digestion of PW. As a matter of fact, when collecting municipal waste, which deals with municipal wastewater. This system was operated
PW and FW are always mixed; for example, a mixture waste of food re- for 485 days and the PW ratio in substrate was changed at the
siduals and tissue is the common garbage from public restaurants. Be- 127th, 202nd, 279th, 366th and 404th day.
cause the separation of these two organic wastes is difficult, the direct Initially, the reactor was dosed with FW as the sole substrate. When
co-digestion of these two wastes provides a convenient way to treat the monitoring parameters of the system became stable, the PW ratio in
both PW and FW. the mixture substrate was elevated to 20% and then to 40%, 50%, 60%,
Although most previous studies on anaerobic digestion have been 70% of the total solid (TS) in the following runs. The artificial FW was
performed under the mesophilic condition of 35 °C (Kong et al., 2018), prepared monthly and contained the mixture of fruit (30%), vegetables
attention is increasingly turning towards thermophilic digestion (36%), staple food (20%), meat and eggs, and all of the organic waste
(Micolucci et al., 2016; Riya et al., 2018). The reported advantages of was crushed into the slurry with a TS of approximately 100 g L−1. The
thermophilic digestion include the much higher permissible organic PW was prepared by mixing office paper, newsprint paper and toilet
loading rate and a higher methane production than that of mesophilic paper at a 1:1:1 ratio according to their air-dried weights. The office
digestion (Jang et al., 2016; Labatut et al., 2014). As the degradation of paper and newsprint paper had been cut into tatters approximately
cellulose is considered the most important point in the anaerobic diges- 2 mm × 10 mm in size by a paper shredder. The mixed PW was stirred
tion of PW (Valdez-Vazquez et al., 2005; Yuan et al., 2014), it is well- with water to make a bucket of PW slurry with the same TS of 100 g L−1.
understood that to improve the cellulose-degrading rate, the optimum In order to maintain the stable growth of microorganisms, FeCl2·H2O,
temperature must be determined. It was recently reported that the CoCl2·6H2O and NiCl2·6H2O were added to the substrate to maintain
highest hydrogen production activity of microorganisms using cellulose the trace elements of 10 mg L−1 Fe, 1 mg L−1 Co and 1 mg L−1 Ni
was at a thermophilic temperature of 55 °C rather than the mesophilic (Qiang et al., 2012).
temperature of 35 °C or a hyper-thermophilic temperature of 70 °C
(Gadow et al., 2016). The results of this comparative study also indi-
2.2. Chemical analysis and calculation
cated that the activation energy of cellulose reached its lowest point
when the cellulose-degrading consortium was cultured at a tempera-
The pH was measured by an HM-30R pH meter (TOA-DKK, Japan).
ture of 55 °C.
Alkalinity was measured by the titrimetric method and the ammonia ni-
Because of the principal role played by anaerobic microorganisms
trogen concentration was determined by 7100 Capillary Electrophoresis
in methane fermentation, an investigation of the composition and
(Agilent Technologies, USA). Chemical oxygen demand (COD) was ana-
variation of microorganisms can lead to a better understanding of
lyzed according to standard methods (Franson et al., 1994). Carbohy-
the mechanisms of co-digestion of FW and PW. While the archaeal
drate, protein and lipid concentrations in both the substrate tank and
communities in anaerobic digestion have been revealed, the specific
the reactor were determined with the same method described in a pre-
functions of bacteria in the co-digestion system of FW and PW re-
vious study (Qin et al., 2018). The concentrations of VFAs (acetate, pro-
main largely unknown. While the variations in the microbial com-
pionate, butyrate, isobutyrate, valerate, isovalerate and hexanoate) and
munity of thermophilic digestion with FW or PW feed as the mono-
lactic acid were measured by a 6890 Series GC system (Agilent, USA). El-
substrate have been investigated (Kamali et al., 2016; Li et al.,
ement components of FW and PW were analyzed by Elementar system
2011; Zamanzadeh et al., 2016; Zamorano-López et al., 2018), no fur-
and all the analytical reagents were purchased from Wako Co. Ltd.,
ther detailed studies have been published on the variations of bacte-
Japan. The volume of generated gas was recorded by a wet gas-meter
rial and archaeal communities as the PW ratio in the thermophilic
and then calculated to the standard condition of 273.15 K, 101325 Pa.
anaerobic digestion system increases for the co-digestion of high-
Biogas components were measured by a GC-8A gas chromatograph
strength FW and PW.
(Shimadzu, Japan) with a 0.4 mL injection volume.
In this study, a thermophilic reactor was set up for the co-digestion
Remove efficiency was calculated based on the following equation:
of FW and PW. One of the objectives of this work was to investigate
the shift in the structure of the microbial community in response to
the long-term performance of the reactor with the increase in the PW C sub −C rea
ratio. Another aim was to seek and identify the functional microorgan- Removal efficiency ð%Þ ¼  100% ð1Þ
C sub
isms which play important roles in cellulose degradation, protein degra-
dation, lactic acid production and together with the syntrophs,
homoacetogens and methanogens which work in the process of me- where Csub: the concentration of composition in the substrate tank; Crea:
thanogenic degradation. Based on these outcomes, a reasonable the concentration of composition in the reactor.
 L. Li et al. / Science of the Total Environment 652 (2019) 709–717 711

Fig. 1. Schematic diagram of lab-scale thermophilic continuously stirred tank reactor (CSTR) system.

Hydrolysis efficiency was calculated based on the following equa- to the manufacturers' instructions. Purified DNA was measured by
tion: Qubit 3.0 ® (Life technologies, USA) and diluted to 4 nmol L−1. Barcode
  added PCR products were sequenced by Illumina Miseq platform.
−1
Hydrolysis efficiency g COD L−1 d
ðC −SC sub Þ−ðC rea −SC rea Þ 2.4. Sequencing data analysis
¼ sub  100% ð2Þ
ðC sub −SC sub Þ
The sequence reads were processed with QIIME (version 1.8.0)
where Csub: the concentration of composition in the substrate tank; (Caporaso et al., 2010), which is widely used in bioinformatics. Raw
SCsub: the concentration of soluble part composition in the substrate data were filtered to remove poor-quality sequences, including se-
tank; Crea: the concentration of composition in the reactor; SCrea: the quences which were longer than 480 bp or shorter than 400 bp based
concentration of soluble part composition in the reactor. on our primers. Operation taxonomic units (OTUs) were generated
The hydrolysis rate was calculated based on the following equation: based on 97% similarity, and the clean sequences were aligned against
  the most recent vision of the Greengenes database, 13_8, (http://
−1
Hydrolysis rate g COD L−1 d qiime.org/home_static/dataFiles.html) before removing the Chimeras
ðTCODsub −SCODsub Þ−ðTCODrea −SCODrea Þ with ChimeraSlayer. Singleton OTUs, which only appeared once in all
¼  100% ð3Þ of the samples, were removed, and sequences were randomly selected
HRT
to unify the sequence number of each sample to 26,000 to remove the
where TCODsub: total COD of substrate; SCODsub: soluble COD of sub- bias of sequencing depth. The Chao1 richness index, the Shannon diver-
strate, TCODrea: total COD of reactor; SCODrea: soluble COD of reactor. sity index and the Simpson diversity index were calculated based on
When the hydrolysis rate of carbohydrate and protein was calculated, clustered sequence data.
the COD was conformed from the corresponding concentration with
the conversion factors of 1.07 and 1.51. 3. Results and discussions

2.3. Microbial population analysis 3.1. Variation of parameters and performance of reactor

15 mL of sludge samples were collected from the outlet of the reac- Based on the results of the elemental analysis, the molecular formu-
tor during the stable phase of each run, respectively on Day 97 (Run I), las of FW and PW in this study were calibrated as C20H33O12N and
181 (Run II), 275 (Run III), 350 (Run IV), 396 (Run V) and 453 (Run VI) C324H497O308N with a C/N ratio of 17.1 and 277.9, respectively. As the
for microbial analysis. Samples were centrifuged at 10,000 rpm for PW ratio increased, the nitrogen in the substrate decreased correspond-
10 min at 4 °C, and 0.15 g of the sediment was used for DNA extraction ingly, resulting in an increase in the C/N ratio from 17.1 to 45.0 as the
by ISOIL for Beads Beating kit (Nippon gene, Japan), then concentrations PW ratio increased from 0% to 70%. The reactor parameters during the
and quality were measured with NanoDrop 2000 (Nanodrop Inc., USA) long-term operation also showed significant variations with the in-
and diluted to 10 ng μL−1 for PCR. V3-V4 fragment of 16S rRNA gene creased PW content. As shown in Table 1, pH decreased from 8.33 to
were amplified with former primer 341F (5′-CCT AYG GGR BGC ASC 7.11, and the concentration of total ammonia nitrogen reduced greatly
AG-3′) and mixed reverse primer 806R & 806R-P (5′-GGA CTA CHV from 2.25 g N L−1 to 0.11 g N L−1, which resulted in a lowering of the
GGG THT CTA AT-3′ & 5′-GGA CTA CCA GGG TAT CTA AG-3′ in ratio of alkalinity from 8.31 g CaCO3 L−1 to 2.02 g CaCO3 L−1. A gradual increase
30:1) (Matsubayashi et al., 2017). The PCR condition was as follows: in the concentration of COD and carbohydrate in the reactor was noted,
30 cycles at 94 °C for 5 s, 50 °C for 5 s, 68 °C for 10 s, and a final extension while the concentration of protein decreased. The lipid concentration of
at 68 °C for 7 min with Low DNA Ex Taq® (TaKaRa, Japan). The purity all the samples was under 0.19 g L−1, the VFAs were under 991 mg L−1
and concentration of the amplified DNA was verified by DNA chip and the gradual decrease in VFAs indicated no accumulation. No lactic
with Bioanalyzer 7500 (Agilent Technologies, USA) and then purified acid was detected from the reactor. The total refractory carbohydrate
with Agencourt® AMPure® XP (Beckman Coulter, Inc., USA) according (including hemicellulose and cellulose) and lignin content accounted
712 L. Li et al. / Science of the Total Environment 652 (2019) 709–717

Table 1
Summary of the important operational parameters in thermophilic anaerobic digestion reactor.

Parameters PW ratio (%)

0 20 40 50 60 70

Influent C/N 17.1 20.4 25.8 29.9 35.8 45.0

pH 8.33 ± 0.10 8.29 ± 0.15 7.94 ± 0.20 7.49 ± 0.09 7.42 ± 0.08 7.11 ± 0.12
Ammonia-N (g L−1) 2.25 ± 0.37 1.86 ± 0.12 1.23 ± 0.30 0.52 ± 0.05 0.33 ± 0.10 0.11 ± 0.07
T-Alkalinity (g-CaCO3 L−1) 8.31 ± 0.25 7.94 ± 0.23 5.48 ± 0.59 3.37 ± 0.11 2.82 ± 0.37 2.02 ± 0.28
COD (g L−1) 28.26 ± 1.57 25.50 ± 2.21 27.97 ± 2.98 30.40 ± 1.87 31.10 ± 2.13 33.08 ± 3.82
Carbohydrate (g L−1) 2.42 ± 0.34 3.11 ± 0.51 5.08 ± 0.54 6.65 ± 0.35 7.17 ± 0.95 9.11 ± 1.57
Protein (g L−1) 11.47 ± 0.96 10.66 ± 0.85 12.08 ± 1.11 9.12 ± 0.81 9.05 ± 1.27 8.91 ± 1.21
Lipid (g L−1) 0.17 ± 0.03 0.19 ± 0.05 0.17 ± 0.02 0.16 ± 0.02 N.D.a N.D.
VFA (mg L−1) 991 ± 12 63 ± 3 128 ± 10 12 ± 3 40 ± 16 190 ± 36
GPR (L L−1 d−1) 2.46 ± 0.32 2.27 ± 0.22 2.30 ± 0.58 2.40 ± 0.22 2.25 ± 0.27 2.07 ± 0.23
CH4 content (%) 61.43 ± 2.75 58.38 ± 2.07 56.89 ± 2.26 53.89 ± 2.05 53.58 ± 1.37 51.51 ± 2.10
CO2 content (%) 38.01 ± 2.80 40.72 ± 1.00 43.81 ± 2.40 45.89 ± 0.80 46.42 ± 1.37 48.20 ± 2.54

COD: total chemical oxygen demand; VFA: volatile fatty acids; GPR: gas production rate.
a
N.D. means not detected.

for 23.38% of the TS in FW, while it was as high as 88.68% of the TS in the total of 25 phyla (account for 98.7% in total sequence) and 40 classes
PW (Table S1). A high gas production rate of 2.46 L L−1 d−1 was ob- (account for 98.7% in total sequence) were obtained and classified
served when the PW ratio was 0%, and this was maintained at above from all the samples: 24 phyla were bacteria while only 1 phylum was
2.07 L L−1 d−1 when the PW ratio was 70%. identified as archaea. The relative proportion of archaea ranged be-
As shown in Fig. 2, the removal efficiency of COD maintained at tween 3.6% and 7.3%, as shown in Table S3.
around 80% in each PW ratio. Carbohydrate removal efficiency was
high at 95.78% when the PW ratio was 0%, and while it slightly de- 3.2.1. Changes in functional bacterial groups
creased with increasing PW ratio, it did not fall below 88.67%. The hy- The community variation trend of bacteria differed significantly be-
drolysis efficiency of COD and carbohydrate slightly fluctuated. tween PW40 and PW50, as shown in Fig. 3. Based on the Student's t-test
However, the protein removal efficiency, around 50% at first, fell to of these two groups (Group 1 include PW0, PW20 and PW40, Group2
just 15% when the PW ratio was 70%. It should be noted that despite include PW50, PW60 and PW70), the p-value is 0.013, which is b0.05,
the slight decrease in the removal efficiencies of carbohydrate, the hy- indicating a statistically significant difference between the low PW
drolysis rate of carbohydrate rose from 1.96 g COD L−1 d−1 to ratio sample group (b50%) and the high PW ratio sample group
2.53 g COD L−1 d−1 along with the increase in the PW ratio. (≥50%). Phylum Firmicutes was overwhelmingly predominant in all
samples, accounting for over 70% in the low PW ratio samples, but rap-
3.2. Overall analysis of bacterial and archaea community idly reduced to about 40% in high PW ratio samples. A decreasing trend
in the relative abundances of the major orders in Firmicutes was noted,
High-throughput sequencing technology based on an Illumina including Thermoanaerobacterales, MBA08, Clostridiales, Lactobacillales
MiSeq platform was used to reveal the overall taxonomic composition and OPB54, the exception was SHA98, which increased greatly in high
of the microbial community. As shown in Table S2, the coverage of PW ratio samples. The second predominant phylum, Thermotogae,
each sample was over 99.3%, exhibiting good coverage in this study. which accounted for about 10% in low PW ratio samples increased to
The Chao1 index, Shannon diversity index and Simpson diversity about 20% in high PW ratio samples. These two phyla have been re-
index were calculated. Among them, the Shannon diversity index de- ported to be predominant in thermophilic anaerobic digestion (Jang
creased from 4.71 to 3.86 with the increase of the PW ratio from 0% to et al., 2016). The third predominant candidate phylum, ‘Atribacteria’,
60% and eventually increased to 4.38 when PW ratio was 70%. The also increased from 5.4% to 24.3% before decreasing to 11.8% in the
changes in the Simpson diversity index exhibited the same trend. A last run. The phylum Synergistetes changed from 3.6% to 8.5%, without

Fig. 2. (a) Removal efficiencies of COD, carbohydrate and protein; (b) hydrolysis efficiencies of COD, carbohydrate and protein; (c) hydrolysis rates of COD, carbohydrate and protein.
 L. Li et al. / Science of the Total Environment 652 (2019) 709–717 713

Fig. 3. Relative proportions of predominant bacterial phyla and predominant orders in phylum Firmicutes.

showing a regular trend. The fifth most abundant phylum in the sam- ethanol (Sekiguchi et al., 2006). The total relative proportion of
ples, Bacteroidetes, increased from 0.3% to 8.7% across the runs. Other syntrophs slightly decreased from 2.7% to 0.9%.
phyla in this study were Acidobacteria, Actinobacteria, Planctomycetes,
Armatimonadetes, Spirochaetes, Caldiserica, Cyanobacteria, Chlorobi,
Fusobacteria, Nitrospirae and some other candidatus phyla. Because all 3.2.2. Changes of methanogenic archaea
of these phyla had an abundance of b0.05% in all the sequences, they The archaeal community structure was much simpler than the bac-
were included in the “Others” category. teria community, since all the sequences were identified as only 6 OTUs
In order to illustrate the ecological function of the community in and all of them belonged to phylum Euryarchaeota as given in Table S5.
more detail, the OTUs which accounted for N0.5% in each sample (se- Among them, denovo3242 was the most predominant OTU, accounting
quence number N 130) were selected and aligned with the NCBI data- for N79.6% of the total archaea, and had 100% similarity with the
base and the Greengenes database. The sum of all of the revealed hydrogenotrophic archaea Methanothermobacter thermautotrophicus
genera accounted for N90% in each sample. These bacterial genera (Wasserfallen et al., 2000). Beside this OTU, there were 4 other OTUs,
were divided into several functional groups according to the role played Methanothermobacter tenebrarum, Methanoculleus thermophilus,
in the hydrolysis, acidogenesis and acetogenesis processes. The details Methanolinea mesophila and Methanomassiliicoccus luminyensis, which
are listed in Table 2. A total of 9 genera belonging to the phylum have been classified as hydrogenotrophic archaea (Dridi et al., 2012;
Firmicutes and Thermotogae were classified as cellulose-degrading bac- Nakamura et al., 2013; Sakai et al., 2012; Zhu et al., 2011). The relative
teria (CDB) based on a sequence identity of above 95% with the cultured proportion of hydrogenotrophic archaea decreased from 6.0% to 2.9%
strains. As shown in Fig. S1(a), the abundance of CDB was at its lowest at at low PW ratios and increased from 2.9% to 6.8% at high PW ratios, as
20.0% in PW20 and peaked at 30.5% in PW50. It has been reported that shown in Fig. S1(b). The predominance of hydrogenotrophic species
all of these genera have (hemi) cellulose-degrading ability (Ben Hania in the archaea community in this study is in good agreement with the
et al., 2012a; Park et al., 2014; Shiratori et al., 2009). Worthy of particu- results reported from previous studies focused on thermophilic meth-
lar attention is that the great increase in the population of denovo975 ane fermentation (Chachkhiani et al., 2004; Hori et al., 2006; Yu et al.,
and denovo6440 with the increasing PW ratio, from 0% to 19.4% and 2014). The second most predominant species was Methanosarcina
from 4.7% to 10.9% respectively (referencing to Table S4). Denovo975 thermophila, which can utilize both acetate and methyl compounds
had an identity of 90% with the species of Caloranaerobacter azorensis, (Aceti and Ferry, 1988). The relative abundance of Methanosarcina
and denovo6440 had an identity of 91% with the strain Candidatus thermophila was stable at each PW ratio and ranged from 0.4% to 0.6%.
OP9 SCG 110722-44 according to the BLAST results. Note that both of
the Caloranaerobacter azorensis species and the strain Candidatus OP9
SCG 110722-44 have been confirmed to have cellulose-degrading abil- 3.3. CDB and protein-degrading bacteria responding to the increased PW
ity (Dodsworth et al., 2013; Wery et al., 2001). Thus, we inferred that ratio
these two OTUs are also responsible for cellulose degradation, and clas-
sified them as putative CDB. Genera of Coprothermobacter in phylum The predominant hydrolyzing groups, which play the important
Firmicutes and Anaerobaculum in phylum Synergistetes were identified roles in the hydrolysis step, were the CDB and the protein-degrading
as protein- or peptides-degrading bacteria. With the decline in the pro- bacteria. The relative abundance of both shifted significantly. The in-
tein concentration in the substrate, the total abundance of protein- crease of PW not only elevated the amount of carbohydrates in the
degrading bacteria decreased steadily from 23.6% to 10.5%, as expected. substrate, but also raised the refractory proportion of carbohy-
Genera Lactobacillus and Leuconostoc were identified as lactic acid bac- drates, such as cellulose. A 25.6% increase in the relative abundance
teria and their total proportions decreased from 7.7% to 1.0% throughout of the putative CDB resulting in a 28.2% increase in the total CDB.
the whole process. As an essential group in anaerobic digestion, Conversely, because the sole contribution of protein was the FW in
syntrophs, which are symbiotic with hydrogenotrophic archaea, play a this study, when the dosage of FW decreased, the protein concentra-
role in oxidizing fatty-acid and generating acetate and hydrogen simul- tion in the substrate also dropped, which led to a 12.7% reduction in
taneously. These were identified as Pelotomaculum, Syntrophaceticus the proteolyzer. The increase in the CDB and decrease in the
and Tepidanaerobacter in this study. Among them, S. schinkii can only proteolyzers resulted in a shift of the metabolism type in the entire
use acetate as an electron donor to generate hydrogen (Westerholm anaerobic digestion system from carbohydrate/protein-degradation
et al., 2010) while T. syntrophicus can also use both lactic acid and to cellulose-degradation.
714 L. Li et al. / Science of the Total Environment 652 (2019) 709–717

3.4. Variation of CDB effecting on the lactic acid bacteria and syntrophs
Optimum pH

6.0–7.7

5.7–6.1
7.3–7.5
7.0–9.0

6.6–7.3
Unlike previous studies on the mono-digestion of cellulose or lignin

6.9

6.5

7.5

7.8
(Callaghan et al., 2002; Yen and Brune, 2007), stable operation was

7
achieved in the co-digestion of PW and FW in this study with a high car-
Ethanol

bohydrate removal efficiency above 90.0% even at a high organic load-

ing rate of 14 g L−1 d−1. The high removal efficiency of carbohydrate
+

+
+
+
+
+

+
was mainly attributed to the overwhelmingly predominance of CDB in
Acetate

the bacterial community. As shown in Fig. S1(a), the relative abundance

of CDB and putative CDB increased by 28.2% when the PW ratio in-
+
+
+
+
+
+
+

+
+

+
+

+
creased from 0% to 70%. The high occupancy of CDB ensured the hydro-
H2 + CO2

lysis of carbohydrate, enabling a stable carbohydrate removal efficiency

during the entire operation. Apart from the increase in the total abun-
dance of CDB, other significant changes in specific species when the
+
+
+
+
+
+
+

+
+

+
+
+
PW ratio was elevated from 40% to 50% greatly influenced the following
Lactic acid

acidogenesis and acetogenesis processes. As shown in Fig. 4, while

Products

Defluviitoga tunisiensis, Candidatus OP9 SCG 110722-44, Clostridium

+

+
+

+
+
putrefaciens, Ruminiclostridium sp., Herbinix hemicellulosilytica and
Defluviitalea sp. were predominant at low PW ratios, the predominant
Function of closest relatives

Monosaccharide-degrading

species shifted to Defluviitoga tunisiensis, Candidatus OP9 SCG 110722-

44, Caloranaerobacter sp. and Ruminiclostridium thermocellum at high
Lactic acid bacteria
Lactic acid bacteria
Lactic acid bacteria
Lactic acid bacteria
Lactic acid bacteria
Lactic acid bacteria

PW ratios. Those CDB which were predominant in high PW ratios

Protein-degrading
Protein-degrading

have been reported to grow at an optimum pH at a neutral condition

b7.0 (Ben Hania et al., 2012b; Johnson et al., 1982; Wery et al., 2001),
Syntrophs
Syntrophs
Syntrophs

while at low PW ratios, the optimum pH for the predominant CDB has
been reported to be relative alkalinity, and over 7.0 (Bouanane-
CDB
CDB
CDB
CDB
CDB
CDB
CDB
CDB
CDB

Darenfed et al., 2011). This is completely different from that in high

PW ratios. Therefore, the variation in the CDB species was considered
Identity (%)
Variation of relative proportion, phylogenetic classification and metabolism descriptions of functional bacteria at different paper waste ratios.

to correspond with the decrease in the pH value from 8.3 to 7.1 in this
study.
100

100

100

100
100
100
100
100

100
95

99
95
97
96
98
95

95

98

99
97
96

The variation of the CDB affected the abundance of lactic acid bacte-
ria. All the lactic acid bacteria in this study were reported to use mono-
Accession number

saccharide as the sole carbon source, while the CDB in this study were
all reported to have monosaccharide-degrading ability. The increase in
MH298782

MG825725

MG674712
KM036189

MF179627

MF093231
MF443391
NR125455

NR136763

NR108234
NR074653

NR074997

NR041320
NR116297
NR040966
LN881581
KT274710

CP016502

CP016502
KF766957

the CDB led to competition with lactic acid bacteria because the CDB
AB10748

consumed the majority of monosaccharide, leaving little monosaccha-

ride for the lactic acid bacteria. Also, with the increase in the PW ratio,
cellulose gradually became the main carbohydrate component rather
Coprothermobacter proteolyticus
Ruminiclostridium thermocellum

Tepidanaerobacter syntrophicus
Ruminiclostridium stercorarium

Pelotomaculum isophthalicicum

than the monosaccharide in the substrate. This greatly reduced the

Defluviitalea isophthalicicum
Ruminiclostridium caenicola
Herbinix hemicellulosilytica

availability of monosaccharide for lactic acid bacteria, resulting eventu-

Fervidobacterium riparium
Hydrogenispora ethanolica

Syntrophaceticus schinkii
Clostridium putrefaciens

CDB: cellulose-degrading bacteria. The sequence number of each sample was unified to 26,000.

ally in a sharp decline in lactic acid bacteria occupancy. Lactic acid is

Lactobacillus helveticus
Defluviitoga tunisiensis

Anaerobaculum mobile
Herbivorax saccincola

Lactobacillus hilgardii

Leuconostoc gelidum
Lactobacillus rossiae

usually reported as one of the organic acids produced in anaerobic di-

Lactobacillus brevis
Lactobacillus sakei
Closest relatives

gestion. The accumulation of lactic acid further acidizes the system

and inhibits methanogenesis. The lactic acid bacteria were also reported
to generate bacteriocins which inhibit the growth of other bacteria
(Noike et al., 2002). In this regard, the decrease in the lactic acid bacteria
inversely benefited the growth of CDB and other bacteria.
The CDB also reduced the proportion of the syntrophs. The CDB at
4622

1088

1909
204

296

832

193

135
70

10

10
77

71

86

13

the high PW ratios, such as Defluviitoga tunisiensis, were more likely to

0

4
0

1
0
1

degrade monosaccharide directly to acetate and hydrogen without pro-

5222

1108

1036
1768

ducing propionate and butyrate (Ben Hania et al., 2012a) than the CDB
196

203
60

12

33
36
48
4

2
6
0

3
0
6

0
0
0

at low PW ratios. The decrease in the VFA concentration also resulted in

the decline in the abundance of syntrophs which utilize VFAs as a sub-
5387

1065
1334

3964
2164
119

175

119
156
232
50

34

64

65

strate to generate acetate and hydrogen. The relative proportion of

8

0
0

1
0
0

syntrophs at high PW ratios decreased to less than half that at low PW

1989

1330

1824
1410

4222
1341
190

278

915

200
147

384

ratios.
40

87
85

75

68
0
0

4
0
0
2736

3855
1736
1055

3.5. Hydrogen production and consumption

PW ratio (%)

178
456

250

926
513

491

132

153
20

69

43

18

84

49
0

High H 2 -producing ability has been reported at a high C/N

3170
1071

5197

1162
599
421
212
103

850

489
258

136
296
276
32

31
26
15

67
16

ratio of 47:1 (Lin, 2004) as opposed to the ratio of 25:1 in nor-

0

mal anaerobic digestion (Parkin and Owen, 1986). This was in

denovo6151
denovo2999

denovo3430
denovo3241
denovo5581
denovo1702

denovo6439

denovo3720
denovo4182
denovo4667
denovo5881
denovo4772

denovo4262

denovo4664
denovo942
denovo473

denovo433

denovo471

denovo623

denovo474

accordance with a C/N ratio of 45:1 when PW was 70% in this

denovo6

research.
OTU ID
Table 2

As explained, the CDB at high PW ratios were capable of directly

degrading monosaccharide to acetate rather than propionate and
 L. Li et al. / Science of the Total Environment 652 (2019) 709–717 715

Fig. 4. Variation of the main identified bacterial and archaea species along with the increase in the paper waste ratio.

butyrate, thus enabling them to produce more hydrogen. This is based perspectives, both the H2-producing and H2-consuming processes
on the following equation: were enhanced at high PW ratios.

C 6 H12 O6 þ 2H 2 O→2CH3 COOH þ 4H2 þ 2CO2 ð4Þ 3.6. Variation of metabolic pathways

C 6 H12 O6 þ 2H 2 →C 2 H 5 COOH þ 2H 2 O ð5Þ In this study, great differences were found in the metabolic path-
ways for low PW ratios and high PW ratios due to the change of the cor-
C 6 H12 O6 →C 3 H7 COOH þ 2H2 þ 2CO2 ð6Þ responding functional bacteria and archaea. The details are provided in
Fig. 5, where the microorganisms which increased in abundance are
This was especially relevant to Defluviitoga tunisiensis, which has been marked in red and those which decreased in abundance are marked in
reported to have high H2-producing ability when fed with cellulose green. Note that the abundance of the microorganisms is roughly indi-
(Maus et al., 2015). With the increase in the total number of H2- cated by the size of the wording. As the condition changed from a low
producing bacteria, the hydrogen production capacity of the system was PW ratio to a high PW ratio, bacteria evolved either in the direction of
increased, and the hydrogenotrophic archaea increased correspondingly. degrading cellulose rather than carbohydrate and protein, or by
It was reported that with the increasing content of H2 in the anaero- degrading monosaccharides directly to acetate and hydrogen. The ar-
bic digestion system, the metabolic pathway from propionate to acetate chaea community changed along with the H2-producing bacteria. That
becomes inhibited, leading to VFA accumulation (Lin, 2004). In this re- is, due to the formation of a new colony structure, when the PW ratio
search, the increasing abundance of hydrogenotrophic archaea elimi- ≥ 50%, the thermophilic anaerobic digestion system obtained apprecia-
nates the adverse effects of excessively produced H2 from cellulose, ble stability and high cellulose tolerance. The microbial community
ensuring that the operation of the system remains stable. From these has been reported to maintain a functional redundancy which can

Fig. 5. Functional groups of microorganisms in the hydrolysis, acidogenesis, acetogenesis and methanogenesis processes. The microorganisms which increased in abundance are marked in
red and those which decreased in abundance are marked in green.
716 L. Li et al. / Science of the Total Environment 652 (2019) 709–717

continue to conduct normal functions even when the community struc- Chachkhiani, M., Dabert, P., Abzianidze, T., Partskhaladze, G., Tsiklauri, L., Dudauri, T.,
Godon, J.J., 2004. 16S rDNA characterisation of bacterial and archaeal communities
ture is forced to change due to environmental stress (Berga et al., 2012). during start-up of anaerobic thermophilic digestion of cattle manure. Bioresour.
The variations of the functional microorganisms in this system and the Technol. 93, 227–232. https://doi.org/10.1016/j.biortech.2003.11.005.
main metabolic pathways illustrate the stability and high cellulose tol- Dodsworth, J.A., Blainey, P.C., Murugapiran, S.K., Swingley, W.D., Ross, C.A., Tringe, S.G.,
Chain, P.S.G., Scholz, M.B., Lo, C.-C., Raymond, J., Quake, S.R., Hedlund, B.P., 2013.
erance the of thermophilic anaerobic co-digestion system. Single-cell and metagenomic analyses indicate a fermentative and saccharolytic life-
style for members of the OP9 lineage. Nat. Commun. 4. https://doi.org/10.1038/
ncomms2884.
4. Conclusions Dridi, B., Fardeau, M.-L., Ollivier, B., Raoult, D., Drancourt, M., 2012. Methanomassiliicoccus
luminyensis gen. nov., sp. nov., a methanogenic archaeon isolated from human faeces.
The change in the functional bacteria and archaea community Int. J. Syst. Evol. Microbiol. 62, 1902–1907. https://doi.org/10.1099/ijs.0.033712-0.
Franson, M.A.H., Association, A.P.H, Association, A.W.W, Federation, W.E, 1994. Standard
during thermophilic anaerobic co-digestion of food waste and paper Methods for the Examination of Water and Wastewater 18th Edition Supplement.
waste was investigated by 16S rRNA gene analysis. The hydrolysis bac- American Public Health Association, American Water Works Association and Water
teria shifted from carbohydrate/protein-degrading bacteria to cellulose- Environment Federation, Washington, D.C Published jointly by.
Gadow, S.I., Jiang, H., Li, Y.-Y., 2016. Characterization and potential of three temperature
degrading bacteria. The evolution of CDB was observed from low to high ranges for hydrogen fermentation of cellulose by means of activity test and 16s
paper waste ratios, enhancing the cellulose tolerance and helping to rRNA sequence analysis. Bioresour. Technol. 209, 80–89. https://doi.org/10.1016/j.
maintain the stability of thermophilic co-digestion of food waste and biortech.2016.02.098.
Hori, T., Haruta, S., Ueno, Y., Ishii, M., Igarashi, Y., 2006. Dynamic transition of a methan-
paper waste. CDB outcompeted lactic acid bacteria in the degradation ogenic population in response to the concentration of volatile fatty acids in a thermo-
of monosaccharides, reducing the proportion of lactic acid bacteria. philic anaerobic digester. Appl. Environ. Microbiol. 72, 1623–1630. https://doi.org/
CDB in high paper waste ratios directly hydrolyzed monosaccharide to 10.1128/AEM.72.2.1623-1630.2006.
Jang, H.M., Ha, J.H., Kim, M.-S., Kim, J.-O., Kim, Y.M., Park, J.M., 2016. Effect of increased
acetate and hydrogen, avoiding the accumulation of propionate and bu-
load of high-strength food wastewater in thermophilic and mesophilic anaerobic
tyrate. The increase in CDB, which obtained high H2-producing ability, co-digestion of waste activated sludge on bacterial community structure. Water
eventually enhanced H2-production, also resulting in the predominant Res. 99, 140–148. https://doi.org/10.1016/j.watres.2016.04.051.
abundance of hydrogenotrophic archaea at a high paper waste ratio. Johnson, E.A., Sakajoh, M., Halliwell, G., Madia, A., Demain, A.L., 1982. Saccharification of
complex cellulosic substrates by the cellulase system from Clostridium thermocellum.
Appl. Environ. Microbiol. 43, 1125–1132.
Jokela, J., Rintala, J., Oikari, A., Reinikainen, O., Mutka, K., Nyrönen, T., 1997. Aerobic
Acknowledgements composting and anaerobic digestion of pulp and paper mill sludges. Water Sci.
Technol. 36, 181–188. https://doi.org/10.1016/S0273-1223(97)00680-X (Sludge
This work was also supported by Japan Society for the Promotion of Rheology).
Kamali, M., Gameiro, T., Costa, M.E.V., Capela, I., 2016. Anaerobic digestion of pulp and
Science (JSPS) KAKENHI Grant No. 18J11397. We gratefully acknowl-
paper mill wastes – an overview of the developments and improvement opportuni-
edge the financial support from the Japan Science and Technology ties. Chem. Eng. J. 298, 162–182. https://doi.org/10.1016/j.cej.2016.03.119.
Agency (JST) in the Japanese-Chinese Research Cooperative Program Kong, Z., Li, L., Kurihara, R., Kubota, K., Li, Y.-Y., 2018. Anaerobic treatment of N,N-
on “Research and Development to Find Solutions to Environmental dimethylformamide-containing wastewater by co-culturing two sources of inocu-
lum. Water Res. 139, 228–239. https://doi.org/10.1016/j.watres.2018.03.078.
and Energy Issues in Urban Areas” (16769220A) and a Grant in-Aid of Labatut, R.A., Angenent, L.T., Scott, N.R., 2014. Conventional mesophilic vs. thermophilic
Tohoku University Division for Interdisciplinary Advanced Research anaerobic digestion: a trade-off between performance and stability? Water Res. 53,
and Education. 249–258. https://doi.org/10.1016/j.watres.2014.01.035.
Leschine, S.B., 1995. Cellulose degradation in anaerobic environments. Annu. Rev.
Microbiol. 49, 399–426. https://doi.org/10.1146/annurev.mi.49.100195.002151.
Appendix A. Supplementary data Li, Y., Park, S.Y., Zhu, J., 2011. Solid-state anaerobic digestion for methane production from
organic waste. Renew. Sust. Energ. Rev. 15, 821–826. https://doi.org/10.1016/j.
rser.2010.07.042.
Supplementary data to this article can be found online at https://doi. Lin, C., 2004. Carbon/nitrogen-ratio effect on fermentative hydrogen production by mixed
org/10.1016/j.scitotenv.2018.10.292. microflora. Int. J. Hydrog. Energy 29, 41–45. https://doi.org/10.1016/S0360-3199(03)
00083-1.
Matsubayashi, M., Shimada, Y., Li, Y.-Y., Harada, H., Kubota, K., 2017. Phylogenetic diver-
References sity and in situ detection of eukaryotes in anaerobic sludge digesters. PLoS ONE 12,
e0172888. https://doi.org/10.1371/journal.pone.0172888.
Aceti, D.J., Ferry, J.G., 1988. Purification and characterization of acetate kinase from Maus, I., Cibis, K.G., Wibberg, D., Winkler, A., Stolze, Y., König, H., Pühler, A., Schlüter, A.,
acetate-grown Methanosarcina thermophila. Evidence for regulation of synthesis. 2015. Complete genome sequence of the strain Defluviitoga tunisiensis L3, isolated
J. Biol. Chem. 263, 15444–15448. from a thermophilic, production-scale biogas plant. J. Biotechnol. 203, 17–18.
Banks, C.J., Chesshire, M., Heaven, S., Arnold, R., 2011. Anaerobic digestion of source- https://doi.org/10.1016/j.jbiotec.2015.03.006.
segregated domestic food waste: performance assessment by mass and energy balance. McCarty, P.L., 1964. Anaerobic waste treatment fundamentals, part III, toxic materials and
Bioresour. Technol. 102, 612–620. https://doi.org/10.1016/j.biortech.2010.08.005. their control. Public Works 95, 91–94.
Ben Hania, W., Godbane, R., Postec, A., Hamdi, M., Ollivier, B., Fardeau, M.-L., 2012a. Micolucci, F., Gottardo, M., Cavinato, C., Pavan, P., Bolzonella, D., 2016. Mesophilic and
Defluviitoga tunisiensis gen. nov., sp. nov., a thermophilic bacterium isolated from a thermophilic anaerobic digestion of the liquid fraction of pressed biowaste for high
mesothermic and anaerobic whey digester. Int. J. Syst. Evol. Microbiol. 62, energy yields recovery. Waste Manag. 48, 227–235. https://doi.org/10.1016/j.
1377–1382. https://doi.org/10.1099/ijs.0.033720-0. wasman.2015.09.031.
Ben Hania, W., Godbane, R., Postec, A., Hamdi, M., Ollivier, B., Fardeau, M.-L., 2012b. Nakamura, Kohei, Takahashi, A., Mori, C., Tamaki, H., Mochimaru, H., Nakamura, Kazunori,
Defluviitoga tunisiensis gen. nov., sp. nov., a thermophilic bacterium isolated from a Takamizawa, K., Kamagata, Y., 2013. Methanothermobacter tenebrarum sp. nov., a
mesothermic and anaerobic whey digester. Int. J. Syst. Evol. Microbiol. 62, hydrogenotrophic, thermophilic methanogen isolated from gas-associated formation
1377–1382. https://doi.org/10.1099/ijs.0.033720-0. water of a natural gas field. Int. J. Syst. Evol. Microbiol. 63, 715–722. https://doi.org/
Berga, M., Székely, A.J., Langenheder, S., 2012. Effects of disturbance intensity and fre- 10.1099/ijs.0.041681-0.
quency on bacterial community composition and function. PLoS ONE 7, e36959. Noike, T., Takabatake, H., Mizuno, O., Ohba, M., 2002. Inhibition of hydrogen fermentation
https://doi.org/10.1371/journal.pone.0036959. of organic wastes by lactic acid bacteria. Int. J. Hydrog. Energy 27, 1367–1371. https://
Bouanane-Darenfed, A., Fardeau, M.-L., Grégoire, P., Joseph, M., Kebbouche-Gana, S., doi.org/10.1016/S0360-3199(02)00120-9 (BIOHYDROGEN).
Benayad, T., Hacene, H., Cayol, J.-L., Ollivier, B., 2011. Caldicoprobacter algeriensis sp. Park, J.-H., Lee, S.-H., Yoon, J.-J., Kim, S.-H., Park, H.-D., 2014. Predominance of cluster I
nov. a new thermophilic anaerobic, Xylanolytic bacterium isolated from an Algerian Clostridium in hydrogen fermentation of galactose seeded with various heat-
hot spring. Curr. Microbiol. 62, 826–832. https://doi.org/10.1007/s00284-010-9789- treated anaerobic sludges. Bioresour. Technol. 157, 98–106. https://doi.org/10.1016/
9. j.biortech.2014.01.081.
Callaghan, F.J., Wase, D.A.J., Thayanithy, K., Forster, C.F., 2002. Continuous co-digestion of Parkin, Gene F., Owen, William F., 1986. Fundamentals of anaerobic digestion of wastewa-
cattle slurry with fruit and vegetable wastes and chicken manure. Biomass Bioenergy ter Sludges. J. Environ. Eng. 112, 867–920. https://doi.org/10.1061/(ASCE)0733-9372
22, 71–77. https://doi.org/10.1016/S0961-9534(01)00057-5. (1986)112:5(867).
Caporaso, J.G., Kuczynski, J., Stombaugh, J., Bittinger, K., Bushman, F.D., Costello, E.K., Qiang, H., Lang, D.-L., Li, Y.-Y., 2012. High-solid mesophilic methane fermentation of food
Fierer, N., Peña, A.G., Goodrich, J.K., Gordon, J.I., Huttley, G.A., Kelley, S.T., Knights, waste with an emphasis on iron, cobalt, and nickel requirements. Bioresour. Technol.
D., Koenig, J.E., Ley, R.E., Lozupone, C.A., McDonald, D., Muegge, B.D., Pirrung, M., 103, 21–27. https://doi.org/10.1016/j.biortech.2011.09.036.
Reeder, J., Sevinsky, J.R., Turnbaugh, P.J., Walters, W.A., Widmann, J., Yatsunenko, T., Qin, Y., Wu, J., Xiao, B., Hojo, T., Li, Y.-Y., 2018. Biogas recovery from two-phase anaerobic
Zaneveld, J., Knight, R., 2010. QIIME allows analysis of high-throughput community digestion of food waste and paper waste: optimization of paper waste addition. Sci.
sequencing data. Nat. Methods 7, 335–336. https://doi.org/10.1038/nmeth.f.303. Total Environ. 634, 1222–1230. https://doi.org/10.1016/j.scitotenv.2018.03.341.
 L. Li et al. / Science of the Total Environment 652 (2019) 709–717 717

Riya, S., Suzuki, K., Meng, L., Zhou, S., Terada, A., Hosomi, M., 2018. The influence of the thermophilic bacterium isolated from a deep-sea hydrothermal vent. Int. J. Syst.
total solid content on the stability of dry-thermophilic anaerobic digestion of rice Evol. Microbiol. 51, 1789–1796. https://doi.org/10.1099/00207713-51-5-1789.
straw and pig manure. Waste Manag. https://doi.org/10.1016/j.wasman.2018.02.033. Westerholm, M., Roos, S., Schnürer, A., 2010. Syntrophaceticus schinkii gen. nov., sp. nov.,
Sakai, S., Ehara, M., Tseng, I.-C., Yamaguchi, T., Brauer, S.L., Cadillo-Quiroz, H., Zinder, S.H., an anaerobic, syntrophic acetate-oxidizing bacterium isolated from a mesophilic an-
Imachi, H., 2012. Methanolinea mesophila sp. nov., a hydrogenotrophic methanogen aerobic filter. FEMS Microbiol. Lett. 309, 100–104. https://doi.org/10.1111/j.1574-
isolated from rice field soil, and proposal of the archaeal family Methanoregulaceae 6968.2010.02023.x.
fam. nov. within the order Methanomicrobiales. Int. J. Syst. Evol. Microbiol. 62, Yen, H.-W., Brune, D.E., 2007. Anaerobic co-digestion of algal sludge and waste paper to
1389–1395. https://doi.org/10.1099/ijs.0.035048-0. produce methane. Bioresour. Technol. 98, 130–134. https://doi.org/10.1016/j.
Sekiguchi, Y., Imachi, H., Susilorukmi, A., Muramatsu, M., Ohashi, A., Harada, H., Hanada, biortech.2005.11.010.
S., Kamagata, Y., 2006. Tepidanaerobacter syntrophicus gen. nov., sp. nov., an anaero- Yu, D., Kurola, J.M., Lähde, K., Kymäläinen, M., Sinkkonen, A., Romantschuk, M., 2014. Bio-
bic, moderately thermophilic, syntrophic alcohol- and lactate-degrading bacterium gas production and methanogenic archaeal community in mesophilic and thermo-
isolated from thermophilic digested sludges. Int. J. Syst. Evol. Microbiol. 56, philic anaerobic co-digestion processes. J. Environ. Manag. 143, 54–60. https://doi.
1621–1629. https://doi.org/10.1099/ijs.0.64112-0. org/10.1016/j.jenvman.2014.04.025.
Shiratori, H., Sasaya, K., Ohiwa, H., Ikeno, H., Ayame, S., Kataoka, N., Miya, A., Beppu, T., Yuan, X., Wen, B., Ma, X., Zhu, W., Wang, X., Chen, S., Cui, Z., 2014. Enhancing the anaerobic
Ueda, K., 2009. Clostridium clariflavum sp. nov. and Clostridium caenicola sp. nov., digestion of lignocellulose of municipal solid waste using a microbial pretreatment
moderately thermophilic, cellulose-/cellobiose-digesting bacteria isolated from me- method. Bioresour. Technol. 154, 1–9. https://doi.org/10.1016/j.biortech.2013.11.090.
thanogenic sludge. Int. J. Syst. Evol. Microbiol. 59, 1764–1770. https://doi.org/ Zamanzadeh, M., Hagen, L.H., Svensson, K., Linjordet, R., Horn, S.J., 2016. Anaerobic diges-
10.1099/ijs.0.003483-0. tion of food waste – effect of recirculation and temperature on performance and mi-
Valdez-Vazquez, I., Sparling, R., Risbey, D., Rinderknecht-Seijas, N., Poggi-Varaldo, H.M., crobiology. Water Res. 96, 246–254. https://doi.org/10.1016/j.watres.2016.03.058.
2005. Hydrogen generation via anaerobic fermentation of paper mill wastes. Zamorano-López, N., Moñino, P., Borrás, L., Aguado, D., Barat, R., Ferrer, J., Seco, A., 2018.
Bioresour. Technol. 96, 1907–1913. https://doi.org/10.1016/j.biortech.2005.01.036. Influence of food waste addition over microbial communities in an anaerobic mem-
Wasserfallen, A., Nölling, J., Pfister, P., Reeve, J., Conway De Macario, E., 2000. Phylogenetic brane bioreactor plant treating urban wastewater. J. Environ. Manag. 217, 788–796.
analysis of 18 thermophilic Methanobacterium isolates supports the proposals to create https://doi.org/10.1016/j.jenvman.2018.04.018.
a new genus, Methanothermobacter gen. nov., and to reclassify several isolates in three Zhu, C., Zhang, J., Tang, Y., Zhengkai, X., Song, R., 2011. Diversity of methanogenic archaea
species, Methanothermobacter thermautotrophicus comb. nov., Methanothermobacter in a biogas reactor fed with swine feces as the mono-substrate by mcrA analysis.
wolfeii comb. nov., and Methanothermobacter marburgensis sp. nov. Int. J. Syst. Evol. Microbiol. Res. 166, 27–35. https://doi.org/10.1016/j.micres.2010.01.004.
Microbiol. 50, 43–53. https://doi.org/10.1099/00207713-50-1-43.
Wery, N., Moricet, J.M., Cueff, V., Jean, J., Pignet, P., Lesongeur, F., Cambon-Bonavita, M.A.,
Barbier, G., 2001. Caloranaerobacter azorensis gen. nov., sp. nov., an anaerobic