Subscriber access provided by CORNELL UNIVERSITY LIBRARY
Article
Effect of Substrate on Identification of Microbial Communities in Poultry Carcass Composting and Microorganisms Associated with Poultry Carcass Decomposition Jie Wang, Xueqing Du, Yitao Zhang, Ting Li, and Xindi Liao J. Agric. Food Chem., Just Accepted Manuscript • DOI: 10.1021/acs.jafc.6b02442 • Publication Date (Web): 22 Aug 2016 Downloaded from http://pubs.acs.org on August 27, 2016
Just Accepted “Just Accepted” manuscripts have been peer-reviewed and accepted for publication. They are posted online prior to technical editing, formatting for publication and author proofing. The American Chemical Society provides “Just Accepted” as a free service to the research community to expedite the dissemination of scientific material as soon as possible after acceptance. “Just Accepted” manuscripts appear in full in PDF format accompanied by an HTML abstract. “Just Accepted” manuscripts have been fully peer reviewed, but should not be considered the official version of record. They are accessible to all readers and citable by the Digital Object Identifier (DOI®). “Just Accepted” is an optional service offered to authors. Therefore, the “Just Accepted” Web site may not include all articles that will be published in the journal. After a manuscript is technically edited and formatted, it will be removed from the “Just Accepted” Web site and published as an ASAP article. Note that technical editing may introduce minor changes to the manuscript text and/or graphics which could affect content, and all legal disclaimers and ethical guidelines that apply to the journal pertain. ACS cannot be held responsible for errors or consequences arising from the use of information contained in these “Just Accepted” manuscripts.
Journal of Agricultural and Food Chemistry is published by the American Chemical Society. 1155 Sixteenth Street N.W., Washington, DC 20036 Published by American Chemical Society. Copyright © American Chemical Society. However, no copyright claim is made to original U.S. Government works, or works produced by employees of any Commonwealth realm Crown government in the course of their duties.
Page 1 of 59
Journal of Agricultural and Food Chemistry
Effect of Substrate on Identification of Microbial Communities in Poultry Carcass Composting and Microorganisms Associated with Poultry Carcass Decomposition Jie Wang‡, Xueqing Du†, Yitao Zhang†, Ting Li†, and Xindi Liao*† ‡ Department of Bioengineering, College of Food Science, South China Agricultural University, Guangzhou, 510642, PR China; † Department of Animal Production, College of Animal Science, South China Agricultural University, Guangzhou, 510642, PR China
*Corresponding Author. Phone: +86 020 85280279. Fax: +86 020 85280740. E-mail:
[email protected] (Xindi Liao).
1 ACS Paragon Plus Environment
Journal of Agricultural and Food Chemistry
2 ACS Paragon Plus Environment
Page 2 of 59
Page 3 of 59
Journal of Agricultural and Food Chemistry
1
ABSTRACT
2
Three composting systems, which consisted of different ratio of chicken manure,
3
sawdust and poultry carcasses, were conducted to investigate the effect of substrate on
4
the identification of microbial communities and microorganisms associated with poultry
5
carcasses decomposition by characterizing of microbial communities and
6
physico-chemical parameters. The physico-chemical and Miseq Illumina sequencing
7
results showed the composition of substrate had a significant effect on the identification
8
and metabolic capabilities of microbial communities in decomposting process. Poultry
9
carcasses might be the potential driver for the identification of bacterial communities in
10
poultry carcass composting, while initial C/N ratio may mainly contribute to the
11
different diversity of fungal communities and the similar dominant microbial
12
communities in treatments. Poultry carcasses and initial C/N ratio could respectively
13
affect the composition and abundance of microorganisms associated with the
14
decomposition of poultry carcasses. Understanding of the potential composting driver
15
could allow development of efficient carcasses degradation system.
16
KEYWORDS: Poultry carcass composting; Substrate; Identification; Microbial
17
community; Functional microorganisms
18 19 20 21 22
3 ACS Paragon Plus Environment
Journal of Agricultural and Food Chemistry
23 24 25
26
As poultry and livestock production operations around the world have become larger
27
and concentrated, the intensive production in one area inevitably results in the
28
production of high number of animal carcasses. In China, for example, approximately
29
one million tons poultry carcasses are produced per year in recent years1. In case of an
30
emergency situation such as disease outbreak or ventilation failure, timely and
31
effectively handling carcasses has become a significant issue.
32
INTRODUCTION
Composting is recognized as an attractive option for disposal of animal carcasses in
33
many regions of the world because it is relatively inexpensive as the supplies and
34
equipment for the process are usually available on farm and could provide biological
35
conversion of carcasses into usable fertilizer products with pathogen destruction 2-4.
36
However, one issue which may curtail the widespread adoption of compost for carcass
37
disposal is the low efficiency with regards to the rate of carcass breakdown.
38
Consequently, improvements on disposal efficiency of animal carcasses are being
39
sought by researches, with several advancements in composting supplies and equipment
40
5,6
41
. Since composting is a microbial consortia-mediated bioprocessing of once-living
42
materials, the degradation efficiency relies greatly on the microbial activity involved
43
through the process, in which the complex organic materials undergo decomposition
44
jointly by microbial communities 7-9. And it is generally accepted that microorganisms
45
tend to act not alone but in association with others and efficient degradation occurs by 4 ACS Paragon Plus Environment
Page 4 of 59
Page 5 of 59
Journal of Agricultural and Food Chemistry
46
the highly collaborative action of microorganisms 10,11. Therefore, the characterization
47
of microorganisms that relate to animal carcass composting could be the starting point
48
for novel biotechnological applications related to the efficient conversion of organic
49
materials.
50
It was reported that the physicochemical and nutritional properties of the raw material
51
being composted are mainly responsible for the identification of the microbial
52
communities, which could affect the metabolic capabilities for transforming organic
53
substrates and then the subsequent microbial functional succession associated to specific
54
processes 12-14. In carcasses composting system, the organic matter is composed of
55
animal carcasses and composting feedstock. Accordingly, the microbial communities
56
involved in the decomposition of organic matter should be respectively affected by
57
animal carcasses and composting feedstock. And the initial composition of animal
58
carcasses and feedstock in substrate may exert their influence by altering the microbial
59
community assembly and function during composting. So, it is important for
60
establishing efficient composting system to character the relationship between the raw
61
material being composted and the microbial diversity or succession and the change in
62
the compost’s physicochemical properties in composting process. In this regard, little
63
information has been reported in the literature regarding the microorganisms or their
64
metabolic activity associated with the decomposition of animal carcasses and
65
composting feedstock at the same time, especially if we consider the potential driver for
66
identification and metabolic capabilities of microbial communities and the functional
67
microorganisms contributing to carcass decomposition 4,6,15.
68 69
In the present study, a deeper exploration of poultry carcass composting was conducted by establishing different composting systems, which consisted of different 5 ACS Paragon Plus Environment
Journal of Agricultural and Food Chemistry
70
ratio of chicken manure, sawdust and poultry carcasses. The aims of this research were
71
focused on: 1) how does the raw material being composted affect the diversity and
72
composition of bacterial-fungal consortia in composting system; 2) analyzing the effect
73
of substrate types on the metabolic capabilities of microbial communities by
74
investigating the physico-chemical properties; 3) identifying the potential functional
75
microorganisms associated with poultry carcasses decomposition. Investigation the role
76
of substrate and their importance in controlling microbial community assembly,
77
biomass turnover and nutrient cycling could allow development of efficient carcasses
78
degradation.
79
80
MATERIALS AND METHODS Composting schedule. Composting was conducted in the special composting bins
81
with mixing device (Figure 1, 60 cm × 60 cm × 80 cm, patent number:
82
ZL201320021724.5) from August to September at the Zengcheng experiment station of
83
South China Agricultural University, Guangzhou, China. Three substrate materials, fresh
84
poultry manure, poultry carcasses weighing 1-1.5kg (18 week-old laying hen) and
85
sawdust were collected from a farm in Zhongshan and a wood processing factory in
86
Zengcheng, respectively.
87
Three treatments with three replicates were carried out in the composting bins,
88
namely group A, group B and group C, each of which has the same volume. Group A
89
was a treatment without the addition of poultry carcass, in which a thorough mixture of
90
poultry manure and sawdust at a ratio of 5:1 (manure weight: sawdust weight) was
91
conducted to achieve the initial C/N of 15. Group B was the treatment that consisted of
92
12 poultry carcasses and the same ratio of poultry manure and sawdust as group A.
6 ACS Paragon Plus Environment
Page 6 of 59
Page 7 of 59
Journal of Agricultural and Food Chemistry
93
Group C was the treatment that contained a thorough mixture of poultry manure and
94
sawdust at a ratio of 4:1 (weight) with the initial C/N of 22 and 12 poultry carcasses. In
95
group B and group C, the compost constituents from bottom to top layer were
96
respectively 25 cm manure-sawdust mixture, six poultry carcasses (every chicken was
97
cut into four pieces and uniformly placed ), 10 cm manure-sawdust mixture, another six
98
poultry carcasses (the same treatment as mentioned above), and 20 cm manure-sawdust
99
mixture.
100
Tissues obtained from the same part of chicken carcasses were cut into 50 g pieces
101
and sealed in nylon bag (100 µm pore size). Then every two nylon bag were
102
respectively suspended adjacent to the carcasses layers in each compost bin of group A
103
and group B to test the degradation rate.
104
Three nanotechnology microporous trachea used for ventilation were fixed on the
105
bottom of each composting bin. Forced aeration was applied to the compost piles, and
106
the ventilation rate was 1.1 m3 h-1bin-1. Fresh air was introduced into each bin through
107
an air pump, which was connected with a flow meter and nanotechnology microporous
108
trachea through tubes (Figure 1). The air supply for each treatment was controlled by a
109
timer to operate for 40 min and then turn off for 20 min. The air moved upward through
110
the compost pile and was exhausted through a outlet above on one side of each bin. An
111
sucking pump inlet was connected with the outlet, and the outlet was connected with a
112
shunt device to divide the exhaust air into two parts. One part gas was collected by a
113
storage bag volume of three M3, and the other part was released into the outside. Pile
114
turning was started from the third day and was conducted before every sampling by
115
turning in two directions for three minutes. A schematic representation of the
116
experimental setup is shown in Figure 1. 7 ACS Paragon Plus Environment
Journal of Agricultural and Food Chemistry
117
Compost sample collection and analysis. Samples from day 0, day 1, day 3, day 6,
118
day 9, day 12, day 16, day 20, day 24, day 28, and day 32 of compost were collected
119
and analyzed. One kilogram of samples collected at different depths and directions were
120
mixed. The samples were divided into two sections and stored at -20 °C and -80 °C for
121
further analysis, which were respectively used for chemical parameters analysis and
122
Illumina sequencing.
123
Physical and chemical characterization consisted of analyses for temperature, pH,
124
moisture content, total nitrogen (TN), organic matter, C:N ratio of non purgeable
125
organic C and total dissolved N, ammonium nitrogen (NH4+-N) and nitrate nitrogen
126
(NO3--N). Ambient temperature and pile core temperature were measured real-time by a
127
PT100 temperature sensor (Platinum Resistance Thermometers, Shanghai Feilong,
128
China) during the composting. Moisture content was measured by drying in an oven at
129
105 °C until no change in dry weight was observed. The pH value was determined on a
130
1:10 compost: water (w/v) suspension using a Five pH meter (Mettler Toledo,
131
Switzerland) after 1 h equilibrium with shaking. The TN, C:N ratio, organic matter,
132
ammonium nitrogen and nitrate nitrogen were respectively determined according to the
133
methods of Chowdhury et al. 16.
134
Gas sample collection and analysis. Gas generated from each composting bin was
135
collected with storage bag for gap sampling. Gap was sampled from day 0, day 1, day 3,
136
day 6, day 9, day 12, day 16, day 20, day 24, day 28, and day 32 after mixing in storage
137
bag and fitted with 1L foil bag. Gap samples were used to determinate the concentration
138
of CO2, CH4 and NH3.
139
The concentration of CO2 and CH4 were analyzed by gas-chromatography with a 8 ACS Paragon Plus Environment
Page 8 of 59
Page 9 of 59
Journal of Agricultural and Food Chemistry
140
7890 A GC-System (Agilent, Santa Clara, CA, USA) connected to an automatic
141
headspace-sample injection syste (Dani HSS 86.50 Dani, Cologno Monzese, Italy).
142
Helium served as the carrier gas. Concentration of CO2 was determined by HP-PLOT/Q
143
column (19095P-Q04, 30 m×0.53 mm×40 µm). Temperatures of the injector oven,
144
column oven and detector were respectively 200, 30 and 200 °C, and the flow rate for
145
nitrogen, H2 and air were 35, 40 and 400 ml/min, respectively. Peaks were identified by
146
comparison with a know concentration of pure CO2. Concentration of CH4 was
147
determined by HP-PLOT/Q column (19091P-Q04, 30 m×0.32 mm×20 µm).
148
Temperatures of the injector oven, column oven and detector were 200, 50 and 250 °C,
149
respectively, and the flow rates for nitrogen, H2 and air were 40, 40 and 400 ml/min,
150
respectively. Peaks were identified by comparison with a know concentration of pure
151
CH4.
152
The concentration of NH3 were measured by PN2000 ammonia gas detector (PNLE
153
Technology Co., Ltd., China). The following formula was used for the calculation of
154
gas emissions:
155 156
F = (C1-C2) V Where, F = Gas emissions (mg), C1 = the concentration of gas in the outlet port (mg /
157
m3), C2 = the concentration of gas in the inlet (mg / m3), V = the total flow of outlet
158
(m³).
159
Microbial community sampling and DNA extraction. The Miseq Illumina
160
sequencing technology was used to explore the microbial ecosystems of three
161
treatments. Microbial communities were respectively sampled from the compost at the 0,
162
1st, 3rd, 6th, 9th, 16th, 24th and 32th day based on the observation of biochemical 9 ACS Paragon Plus Environment
Journal of Agricultural and Food Chemistry
163
changes in composing system. Samples were collected at different depths and directions
164
and mixed. The total microbial DNA of triplicate compost samples were extracted by
165
EZNATM Soil DNA extraction kit produced by OMEGA company. DNA products of
166
the three replicate samples were mixed in equal concentrations and aliquoted for
167
analysis by MiSeq high-throughput sequencing.
168
Illumina sequencing and quality filtering. PCR primers were designed to amplify a
169
459-bp region across the V3-V4 region of the 16S rRNA gene and a 200-600 bp
170
spanning ITS1 region of the 18S rRNA gene, respectively. The 16S rRNA gene was
171
amplified with primers 341F (GTACTCCTACGGGAGGCAGCA) and R806
172
(GTGGACTACHVGGGTWTCTAAT), as previously described 17. 18S primers were
173
designed as F (CTTGGTCATTTAGAGGAAGTAA) and R
174
(R:TGCGTTCTTCATCGATGC) based on Courtin et al.18and Ihrmark et al.19.
175
The amplification conditions were as follows: 94 °C for 5 min, 30 cycles of PCR
176
(94 °C for 30 s, 56 °C for 30 s, and 72 °C for 30s ), and a final elongation step of 72 °C
177
for 10 min. The samples were kept in 4 °C. Libraries were constructed and quantitated
178
using a KAPA Hyper Prep Kit (Illumina). Equimolar amounts of each of the 96 libraries
179
were pooled and submitted for sequencing on a MiSeq Sequencer using a 250-bp
180
paired-end protocol at the Honor tech company of Beijing.
181
The sample containing indexed ampllicons was loaded onto the Miseq reagent, and
182
then onto the instrument along with the flow cell. Automated cluster generation and
183
paried-end sequencing with dual index reads was performed. Total run time for this
184
2×300 bp run was 40 hours. The primary date analysis were performed directly on
185
Illumina MiSeq platform. Removing low quality data reads according to the Q20, 90% 10 ACS Paragon Plus Environment
Page 10 of 59
Page 11 of 59
Journal of Agricultural and Food Chemistry
186
of the filter standard. The quality score of most reads were above 20. After splicing, 16S
187
rRNA sequences less than 200 bp and 18S rRNA sequences less than 100 bp in length
188
were removed using mothur software, where the control condition is minlength = 100,
189
maxhomop=10 and maxambig=0. The date were used for further information analysis.
190
16S rRNA and 18S rRNA gene Sequence Analysis. The 16S rRNA dataset were
191
processed and analyzed using QIIME software (v1.8.0). After optimizing, sequences
192
with ≥ 97% similarities were selected and grouped into the Operational Taxonomic
193
Units (OTU) for species classification. Alpha diversity was estimated on the abundance
194
information of OTUs in each sample by calculating the Shannon index, Chao1 index,
195
Phylogenetic diversity and observed number of species.
196
The mothur software (v1.8.0) was used to process the 18S rRNA dataset obtained
197
from Illumine sequencing. The optimized sequences were grouped into OTUs based on
198
a 97 % sequence similarity. Community diversity indices, including Shannon, Chao1,
199
ace, Coverage and Simpson, were calculated to estimate the diversity within each
200
treatment.
201
Beta diversity patterns were analyzed to compare the community diversity among
202
treatments by performing UniFrac-based Principal Coordinates Analysis (PCoA), and
203
the heatmaps were drawn by R (v3.1.1) software. Taxonomic classifications at the
204
phylum and genus or species level of each OTU were done using RDP classifier and
205
BLAST algorithms against the Greengenes (16S rRNA), UNITE and GenBank (ITS1)
206
databases. The possible functional profiles of microbial consortia in three treatments
207
were speculated by analyzing the OTU overlap among samples from Venn diagrams
208
drawn by R (v3.1.1) software. 11 ACS Paragon Plus Environment
Journal of Agricultural and Food Chemistry
Statistical analysis. All the physicochemical properties of samples were carried out
209 210
in triplicates, and the results were expressed as mean values and standard errors. A
211
one-way ANOVA analysis was used to examine the effects of substrate on
212
physico-chemical parameters and microbial communities. Significant differences
213
between the treatments were compared by Tukey test at p ˂ 0.05 using the general
214
linear model procedures of the SPSS v. 18.0 programs.
215
216
RESULTS AND DISCUSSION Mass loss. Mass loss represented the degradation degree of material composted to
217
some extent 20. As observed in Table 1, the highest loss rates of mass was found for
218
group C, followed by group B and group A. Poultry carcasses in group B and group C
219
started to degrade rapidly within the first 10 days, and no whole carcasses were evident
220
after 20 days. At the end of composting, only some solid material was present in the
221
form of bones, and the loss rates of poultry carcasses in these two treatments were about
222
98%. The chicken degradation rates of group B and group C were lower than the loss
223
rates of poultry carcasses, which might be ascribed to the inadequate contact of chicken
224
inside the nylon bags with substrate. Degradation rates of chicken in group C were
225
higher than those of group B, especially after 16 days of composting (Figure 2A). These
226
results indicated that the C:N ratio may exerted more influence on substrate
227
decomposition, but poultry carcasses may promote the decomposition of substrate,
228
although no significant difference was observed between treatments (p > 0.05).
229
Temperature, pH and moisture content. Temperature course in three treatments fit
230
the typical evolution of composting processes operated with turning treatments (Figure
231
2B). The thermophilic phase (> 55 °C) started within 12 h and lasted for about 16 days 12 ACS Paragon Plus Environment
Page 12 of 59
Page 13 of 59
Journal of Agricultural and Food Chemistry
232
in all treatments. The peak temperature of group B and group C was obtained at 10 days,
233
but group A was on the 7 th day. In the first 16 days of composting, average temperature
234
of group A was significant lower (p < 0.05) than those of group B and group C. The
235
subsequent cooling stage of group B registered values slightly higher than group A and
236
group C.
237
A increasing tendency in pH curves was observed in three treatments except for the
238
slight decreasing during the 3th to 6th days in group B and group C (Figure 2C). During
239
6 to 16 days, pH in group A was significantly higher (p < 0.05) than the other two
240
treatments. Thereafter, no significant difference (p > 0.05) was observed among three
241
treatments. Moisture content in three treatments followed similar decreasing profile
242
(Figure 2D), with differences observed in decreasing amplitude on trail day, where
243
group A exhibited significantly lower (p < 0.05) moisture contents than group B and C.
244
Rapid decomposition of organic matter could generate more heat, and high water
245
retention of substrate also could increase the thermophilic temperature of composting 12.
246
Accordingly, we observed that samples from group B and C had higher degradation
247
rates, temperature profile and moisture content. The results indicated that initial
248
composition of raw material, and especially the combined addition of carcasses, resulted
249
in the differences on the factors of the surroundings.
250
Organic matter and C:N ratio. It was found in other composting systems that rapid
251
biodegradation of organic matter could result in the disappearance of organic carbon 21.
252
In the present study, group C had the highest initial organic matter content. During 6 to
253
20 days, the highest organic matter content was found for group B, followed by group C
254
and group A, respectively (Figure 2E). The highest initial value and a decreasing
255
tendency was observed on C:N ratio of group C (Figure 2F). The C:N ratio of group A 13 ACS Paragon Plus Environment
Journal of Agricultural and Food Chemistry
256
and group B presented a slightly increasing and then significantly decreasing tendency.
257
At the last stage, group B had the lowest C:N ratio, while group A had the highest C:N
258
ratio. In addition, the standard errors on day 3 and day 6 were higher in Figure 2E and
259
Figure 2F, which might result from the partially nonuniform mixture of poultry manure
260
and sawdust.
261
TN content, ammonium nitrogen and nitrate nitrogen content. The change of TN
262
content was shown in Figure 3A, where a decreasing tendency was observed on group A,
263
which differed from the increasing first and then decreasing profile of group B and
264
group C. According to Moorhead and Sinsabaugh 22, nitrogen seems to accelerate
265
decomposition due to a higher microbial nutrient demand. The initial TN content in
266
group C was lowest, but in later composting stages, group C exhibited the highest
267
degradation rates and TN content.
268
Previous findings showed that high NH4 + concentrations generally promote N losses
269
via volatilization 8,23. Ammonium nitrogen (NH4+-N) content in three treatments all
270
increased significantly during composting (Figure 3B). Accordingly, group A had the
271
highest NH4+-N content and the lowest TN content. Three treatments had similar initial
272
value and increasing profile of nitrate nitrogen (NO3--N) (Figure 3C). The highest and
273
lowest NO3--N contents were observed in group A and group B, respectively. Nutrients
274
in compost system could affect the conversion of ammonia to nitrate and the retain of N
275
in compost 24. Compared to group A, treatments with addition of carcasses had lower
276
NH4 + and NO3 -N concentrations by the end of composting, indicating a much more
277
dynamic N transformation during composting.
14 ACS Paragon Plus Environment
Page 14 of 59
Page 15 of 59
Journal of Agricultural and Food Chemistry
278
Gaseous emissions. A rapid increase of initial daily NH3 emission volume was
279
registered in three treatments on the second day (Figure 3D). Subsequently, the daily
280
NH3 emission volume decreased in all treatments. The highest total NH3 emission
281
volume was observed in group A, while the lowest was in group C. The most important
282
factors that influence the emission of NH3 during composting are pH value and C:N
283
ratio 25, which supported by our results, where the highest total NH3 emission volume
284
was accompanied by high pH value and low C:N ratio in group A.
285
The CO2 emissions of three treatments had similar change profile and fluctuated over
286
time during composting (Figure 3E). Group C exhibited highest daily and total CO2
287
emission, followed by group B and group A until the last stage. The change of CO2
288
emission was found to be consistent with the changes of C:N ratio, which ranked in
289
order of group C, group B and group A.
290
The daily CH4 emission significantly decreased (P < 0.05) on the first day, and then
291
slight increased during day 6 to day 12, followed by a slowly decreasing trend (Figure
292
3F). Production of CH4 occurs through the degradation of lipids, carbohydrates and
293
proteins in anaerobic conditions and reduces in aerobic condition 26. We speculated that
294
central region of compost group was under anaerobic conditions due to the lack of
295
oxygen at the beginning, and thus resulted in the highest CH4 production. With the
296
supply of fresh air, CH4 production significant decreased. In addition, addition of
297
carcasses in group B and group C resulted in larger gap among substrate material, which
298
might be responsible for the less CH4 production.
299
Analysis of the community diversities of three microbial consortia. After sequence
300
quality filtering and removal of failed samples with low numbers of sequences,
15 ACS Paragon Plus Environment
Journal of Agricultural and Food Chemistry
301
the Miseq Illumina sequence dataset included 24 samples and 642242 16S rRNA
302
sequences that generated 60212 OTUs based on a 97 % sequence similarity . The 18S
303
rRNA dataset included 23 samples and 2574729 18S rRNA sequences that generated
304
9884 OTUs.
305
Alpha diversity. Alpha diversity analysis showed that both richness and diversity
306
estimates of the 16S rRNA and 18S rRNA sequences were variable across sample days.
307
16S rRNA diversity indices showed that richness and diversity of samples from group B
308
and group C exhibited the similar change pattern, where two increasing peaks were
309
observed during the sampling days 1- 6 and days 16-24 (Table S1), which was different
310
from the steady increasing trend of group A. With respect to fungal community,
311
diversity indices of samples from group B presented similar change profile with those of
312
group A, but different from those of group C (Table S2).
313
Beta diversity. The abundance-unweighted UniFrac analysis of bacterial and fungal
314
communities (Figure 4A, C) revealed the substrate-specific clusters, and the separation
315
of samples at all sample sites in three treatments was significant (P < 0.05), reflecting
316
the effect of substrate composition on the identification of microbial communities.
317
However, the abundance-weighted analysis (Figure 4B, D) results showed that
318
microbial communities clusters formed in a substrate-independent pattern, where three
319
treatments presented similar bacterial and fungal diversity at most sample sites. We
320
speculated that the similar microbial diversity may be ascribed to the decomposition of
321
poultry manure and sawdust commonly contained in three treatments. In addition, it was
322
found that weighted UniFrac distances of bacterial communities between samples from
323
group A and group B during 0-16 days were lower than between samples from group A
324
and group C or from group B and group C. While the shorter distances of fungal 16 ACS Paragon Plus Environment
Page 16 of 59
Page 17 of 59
Journal of Agricultural and Food Chemistry
325
communities were observed between samples from group B and group C when
326
compared with the samples from group A and group B or from group A and group C.
327
In conjunction with the analysis results of alpha diversity, we speculated that the
328
presence of carcasses might be the potential contributing factor to the different change
329
patterns of bacterial diversity between group A and the other two groups, and the initial
330
C/N ratio might contribute to the diversity difference of fungal community among
331
treatments, supporting previous reports that raw material being composting could affect
332
the richness and diversity of microbial community 27,28.
333
Analysis of community structure and composition of three dominant microbial
334
consortia. Three treatments displayed the involvement of similar phyla of bacteria and
335
fungi, but differed in relative abundance at all sampled days (Figure 5A, C) . For
336
bacterial community, the Firmicutes, members of which were considered to play an
337
important role in degrading organic materials in composting process, were dominant in
338
these treatments. Proteobacteria was the next most abundant phylum, followed by
339
Actinobacteria and Bacteroidetes. The similar structure phylum has been observed in
340
many bacterial diversity surveys of composting, including carcasses decomposition 7,15.
341
For fugal community, Ascomycota phyla were predominant during 0-16 days, thereafter
342
the relative abundance decreased sharply. At the last stage of 24 to 32 days, the relative
343
abundance of Basidiomycota and Glomeromycota significantly increased in samples
344
from group A and group B. The fungal community structure was different from the
345
dominant fungal community observed in a mouse model system 29.
346 347
A closer look of bacterial community at genus level (Figure 5B) revealed that beginning at the initial day (0), three treatments were populated by several genera,
17 ACS Paragon Plus Environment
Journal of Agricultural and Food Chemistry
348
including Corynebacterium, Virgibacillus, Bacteroides and Gallicola. Corynebacterium
349
and Gallicola were significantly more abundant in group A and group B than in group C,
350
whereas Virgibacillus and Bacteroides were significantly more abundant in group C
351
than in the other two group. According to previous reports, some species from
352
Corynebacterium genus are generally recognized as compost inhabitants associated with
353
the transformation of nitrogen transformation in composting 30. Virgibacillus members
354
were found to be more abundant in cellulose decomposing 31. These data corroborated
355
the evidences analyzed above of chemical properties, where group A and B had higher
356
initial TN content and exhibited more efficient N transformation, while group C
357
contained more cellulase because of the high content of sawdust. During 6-16 days,
358
similar communities of microbes dominated within each treatments, but the relative
359
abundance and change pattern were different. Among these dominant genera, members
360
of Corynebacterium and Lactobacillus could break down lipids and complex
361
carbohydrates associated with animal carcass tissue 32,33. Caldicoprobacter is generally
362
recognized as bacterial community associated with keratin decomposition of carcasses
363
34
364
in group A and highest in group B. Group C contained higher abundances of
365
Tepidimicrobium and Lactobacillus. In addition, it was found that group C and group B
366
had the same dominant taxa and similar change profiles of Caldicoprobacter and
367
Lactobacillus, which may contribute to decomposition of keratin and lipids, suggesting
368
that this same functional community may generalize across composting systems with
369
the induction of poultry carcasses. After 16 days composting, the prominent genera
370
sharply decreased in three treatments.
371
. The abundances of Corynebacterium and Caldicoprobacter were significantly lowest
As shown in Figure 5D, two key stages were observed before and after the sampling 18 ACS Paragon Plus Environment
Page 18 of 59
Page 19 of 59
Journal of Agricultural and Food Chemistry
372
day 16 on change profile of fungal community at specie level. During 0-16 days, fungal
373
community at all sampled sites were dominated by Aspergillus, Trichocomaceae and
374
Penicillium at specie level, with other minor fungal populations of Ophiostoma and
375
Lasiodiplodia in three treatments, of which species of Aspergillus and Penicillium were
376
thought to play a major role in the decomposition of cellulose and lipids, respecitvely
377
35-37
378
community at all sites became dominated by a nematode, O. tipulae, in the family
379
Rhabditidae and fungal community composition did not change significantly over time
380
in mouse model system 29. After 16 days composting, Several genera, such as Gigaspora,
381
Pleurotus and Cladosporium, differently dominated in samples. Pleurotus and
382
Cladosporium are known to be common members of the fresh poultry manure or
383
sawdust 38. Carcasses decomposition might be responsible for the shifts because of the
384
breakdown of tissue, which release an ammonia-rich, high nutrient fluid that alters both
385
the pH and nutrient content of the substrate 39. Accordingly the predictable spike in pH
386
(Figure 1B) and total nitrogen and the rapid degradation rates of chicken were observed
387
on the samples collected at the 16 days from group B and group C. In conjunction with
388
the analysis results of bacterial community, we speculated that the initial C/N ratio
389
might the potential contributor for the identification of dominant microbial
390
communities.
391
. These data differed from the findings of Metcalf et al., who found that the fungal
Predicted (potential) functional profiles of microbial consortia. The Venn
392
diagrams represented the shared bacterial and fungal OTUs between treatments as well
393
as the unique OTUs in each treatment (Figure 6). The OTUs shared by three treatments
394
were observed throughout composting. The common in three treatments was the
395
comprising of poultry manure and sawdust in substrate, indicating that the shared 19 ACS Paragon Plus Environment
Journal of Agricultural and Food Chemistry
396
microbial communities might be the feedstock-specific microbiota. Most of the OTUs
397
unique to group A and B and the OTUs unique to group A and C sharply decreased after
398
6 days composting. Group A contained the same ratio of poultry manure and sawdust as
399
group B, while group B had the same amounts of poultry carcasses with group C.
400
Results illustrated that these shared microbiota may originally present in substrate, and
401
decreased because of the poor adaption to severe environment. The OTUs unique to
402
each treatment were found to be generally abundant in the first six days or during the
403
last composting stage. According to previous reports, some microorganism are
404
generally considered to play roles on the adaptation to the composting conditions 9.
405
Thus, we speculated that microbial consortia unique to each treatment might consist of
406
microbiota originally presenting in the substrate or coming from the surroundings and
407
those adapting to the unique composting condition in each treatment changed with
408
composting.
409
Group B and C shared 1158 bacterial OTUs and 109 fungal OTUs, most of which
410
exhibited higher relative abundance during 3-24 days, which was in accordance with the
411
stage of poultry carcasses decomposition. And the dominant microbiota of these shared
412
OTUs belonged to carcasses-specific microbial phyla, indicating that the microbial
413
communities may contribute to the decomposition of poultry carcasses. In addition, the
414
relative abundances of these shared OTUs were different between group B and group C.
415
These results suggested that poultry carcasses and initial C/N ratio could respectively
416
affect the composition and abundance of microorganisms associated with the
417
decomposition of poultry carcasses.
418
In conclusion, this is the first report on microorganisms and their metabolic activity
419
associated with the decomposition of animal carcasses and composting feedstock at the 20 ACS Paragon Plus Environment
Page 20 of 59
Page 21 of 59
Journal of Agricultural and Food Chemistry
420
same time. It is important for establishment of efficient composting system to character
421
the relationship between the raw material being composted and the microbial
422
community and metabolic activity in composting process. Further investigation is still
423
necessary to determine the optimum growth conditions of the functional microbiota
424
associated with poultry carcass decomposition.
425
426
Funding
427
This work was funded by the earmarked fund for Modern Agro-industry Technology
428
Research System, China (Grant Nos. CARS-41); and the National Special Fund for
429
Agro-scientific Research in the Public Interest, China (Grant Nos. 201303091).
430
Notes
431
The authors declare no competing financial interest.
432
Supporting Information Available
433
Analysis of microbial alpha-diversity indices was found in Supporting Information. This
434
material is available free of charge via the Internet at http://pubs.acs.org.
435
436
(1) National Bureau of Statistics of the People’s Republic of China. Agriculture,
437
forestry, animal husbandry and fishery statistics. In China Rural Statistical Yearbook,
438
China Statistics Press: Beijing. China, 2015; 174.
439
(2) Akdeniz, N.; Koziel, J.A.; Ahn, H.K.; Glanville, T.D.; Crawford, B.P.; Raman, D.R.
440
Air sampling and analysis method for volatile organic compounds (VOCs) related to
AUTHOR IMFORMATION
REFERENCES
21 ACS Paragon Plus Environment
Journal of Agricultural and Food Chemistry
441
field-scale mortality composting operations. J. Agric. Food Chem. 2009, 57: 5658-5664.
442
(3) Fongaro, G.; Viancelli, A.; Magri, M.E.; Elmahdy, E.M.; Biesus, L.L.; Kich, J.D.;
443
Kunz, A.; Barardi, C.R.M. Utility of specific biomarkers to assess safety of swine
444
manure for biofertilizing purposes. Sci. Total Envrion. 2014, 1: 479-480, 277-283.
445
(4) McLaughlin, M.R.; Brooks, J.P.; Adeli, A.; Miles, D.M., 2015. Using broiler litter
446
and swine manure lagoon effluent in sawdust-based swine mortality composts: Effects
447
on nutrients, bacteria, and gaseous emissions. Sci. Total Environ. 2015, 532: 265-280.
448
(5) Glanville, T.D.; Ahn, H.K.; Richard, T.L.; Harmon, J.D.; Reynolds, D.L.; Akinc, S.
449
Effect of envelope material on biosecurity during emergency bovine mortality
450
composting. Bioresour. Technol. 2013, 130: 543-551.
451
(6) Glanville, T.D.; Ahn, H.; Akdeniz, N.; Crawford, B.P.; Koziel, J.A. Performance of
452
a plastic-wrapped composting system for biosecure emergency disposal of
453
disease-related swine mortalities. Waste Manage. 2016, 48: 483-491.
454
(7) Pechal, J.L.; Crippen, T.L.; Tarone, A.M.; Lewis, A.J.; Tomberlin, J.K.; Benbow,
455
M.E. Microbial community functional change during vertebrate carrion decomposition.
456
PLoS One 2013, 8: e79035.
457
(8) Brandstatter, C.; Keiblinger, K.; Wanek, W.; Zechmeister-Boltenstern, S. A closeup
458
study of early beech litter decomposition: potential drivers and microbial interactions on
459
a changing substrate. Plant Soil. 2013, 371: 139-154.
460
(9) Zuroff, T.R.; Curtis, W.R. Developing symbiotic consortia for lignocellulosic
461
biofuel production. Appl. Microbiol. Biotechnol. 2012, 93: 1423-1435.
22 ACS Paragon Plus Environment
Page 22 of 59
Page 23 of 59
Journal of Agricultural and Food Chemistry
462
(10) Bokulich, N.A.; Mills, D.A. Facility-specific "house" microbiome drives microbial
463
landscapes of artisan cheesemaking plants. Appl. Environ. Microb. 2013, 79:
464
5214-5223.
465
(11) Embree, M.; Nagarajan, H.; Movahedi, N.; Chitsaz, H.; Zengler, K. Single-cell
466
genome and metatranscriptome sequencing reveal metabolic interactions of an
467
alkane-degrading methanogenic community. ISME J. 2014, 8: 757-767.
468
(12) Lyons, G.A.; Sharma, H.S.S.; Kilpatrick, M.; Cheung, L.; Moore, S. Monitoring of
469
changes in substrate characteristics during mushroom compost production. J. Agric.
470
Food Chem. 2006, 54: 4658-4667.
471
(13) Bonilla, N.; Vida, C.; Martinez-Alonso, M.; Landa, B.B.; Gaju, N.; Cazorla, F.M.;
472
de Vicente, A. Organic amendments to avocado crops induce suppressiveness and
473
influence the composition and activity of soil microbial communities. Appl. Environ.
474
Microb. 2015, 81: 3405-3418.
475
(14) Zhang, L.; Jia, Y.; Zhang, X.; Feng, X.; Wu, J.; Wang, L.; Chen, G. Wheat straw:
476
An inefficient substrate for rapid natural lignocellulosic composting. Bioresour. Technol.
477
2016, 209: 402-406.
478
(15) Tkachuk, V.L.; Krause, D.O.; Knox, N.C.; Hamm, A.C.; Zvomuya, F.; Ominski,
479
K.H.; McAllister, T.A. Targeted 16S rRNA high-throughput sequencing to characterize
480
microbial communities during composting of livestock mortalities. J. Appl. Microbiol.
481
2014, 116: 1181-1194.
482
(16) Chowdhury, M.A.; Neergaard, A.D.; Jensen, L.S. Composting of solids separated
483
from anaerobically digested animal manure: Effect of different bulking agents and
23 ACS Paragon Plus Environment
Journal of Agricultural and Food Chemistry
484
mixing ratios on emisssions of greenhouse gases and ammonia. Biosyst. Eng. 2014, 124:
485
63-77.
486
(17) Takahashi, S.; Tomita, J.; Nishioka, K.; Hisada, T.; Nishijima, M. Development of
487
a Prokaryotic Universal Primer for Simultaneous Analysis of Bacteria and Archaea
488
Using Next-Generation Sequencing. PLoS One 2014, 9: e1055928.
489
(18) Courtin, F.; Stone, W.B.; Risatti, G.; Gilbert, K.; Kruiningen, H.J.V. Pathologic
490
findings and liver elements in hibernating bats with white-nose syndrome. Vet. Pathol.
491
2010, 47: 214-219.
492
(19) Ihrmark, K.; Bodeker, I.T.; Cruz-Martinez, K.; Friberg, H.; Kubartova, A.;
493
Schenck, J.; Strid, Y.; S-tenlid, J.; Brandstrom-Durling, M.; Clemmensen, K.E.; Lindahl,
494
B.D. New primers to amplify the fungal ITS2 region--evaluation by 454-sequencing of
495
artificial and natural communities. FEMS Microbiol. Ecol. 2012, 82: 666-677.
496
(20) Jacob, M.; Viedenz, K.; Polle, A.; Thomas, F.M. Leaf litter decomposition in
497
temperate deciduous forest stands with a decreasing fraction of beech (Fagus sylvatica).
498
Oecologia 2010, 164: 1083-1094.
499
(21) Cleveland, C.C.; Neff, J.C.; Townsend, A.R.; Hood, E. Composition, dynamics,
500
and fate of leached dissolved organic matter in terrestrial ecosystems: results from a
501
decomposition experiment. Ecosystems 2004, 7: 175-285.
502
(22) Moorhead, D.L.; Sinsabaugh, R.L. A theoretical model of litter decay and
503
microbial interaction. Ecol. Monogr. 2006, 76: 151-174.
504
(23) Chen, Y.N.; Zhou, W; Li, Y.P.; Zhang, J.C.; Zeng, G.M.; Huang, A.Z.; Huang, J.X.
505
Nitrite reductase genes as functional markers to investigate diversity of denitrifying 24 ACS Paragon Plus Environment
Page 24 of 59
Page 25 of 59
Journal of Agricultural and Food Chemistry
506
bacteria during agricultural waste composting. Appl Microbiol Biotechnol, 2014, 98:
507
4233-4243.
508
(24) Güsewell, S.; Gessner, M.O. N: P ratios influence litter decomposition and
509
colonization by fungi and bacteria in microcosms. Funct. Ecol. 2009, 23: 211-219.
510
(25) Jiang, T.; Schuchardt, F.; Li, G.; Guo, R.; Zhao, Y. Effect of C/N ratio, aeration
511
rate and moisture content on ammonia and greenhouse gas emission during the
512
composting. J. Environ. Sci. 2011, 23: 1754-1760.
513
(26) Parravicini, V.; Svardal, K.; Hornek, R.; Kroiss, H. Aeration of anaerobically
514
digested sewage sludge for COD and nitrogen removal: optimization at largescale.
515
Water Sci. Technol. 2008, 57: 257-264.
516
(27) Jiménez, D.J.; Dini-Andreote, F.; Elsas, J.D.V. Metataxonomic profiling and
517
prediction of functional behaviour of wheat straw degrading microbial
518
consortia. Biotechnol. Biofuels 2014, 7: 1-18.
519
(28) Lopez-Gonzalez, J.A.; Vargas-Garcia Mdel, C.; Lopez, M.J.; Suarez-Estrella, F.;
520
Jurado, M.D.M.; Moreno, J. Biodiversity and succession of mycobiota associated to
521
agricultural lignocellulosic waste-based composting. Bioresour. Technol. 2015, 187:
522
305-313.
523
(29) Metcalf, J.L.; Wegener Parfrey, L.; Gonzalez, A.; Lauber, C.L.; Knights, D.;
524
Ackermann, G.; Hum-phrey, G.C.; Gebert, M.J.; Van Treuren, W.; Berg-Lyons, D.;
525
Keepers, K.; Guo, Y.; Bullar-d, J.; Fierer, N.; Carter, D.O.; Knight, R. A microbial
526
clock provides an accurate estimate of the postmortem interval in a mouse model
527
system. Elife 2013, 2: e01104.
25 ACS Paragon Plus Environment
Journal of Agricultural and Food Chemistry
528
(30) Schröder, J.; Maus, I.; Trost, E.; Tauch, A. Complete genome sequence of
529
Corynebacterium variabile DSM 44702 isolated from the surface of smear-ripened
530
cheeses and insights into cheese ripening and flavor generation. BMC genomics 2011,
531
12: 545.
532
(31) Chen, C.Y.; Mei, H.C.; Cheng, CY.; Lin, J.H.; Chung, Y.C. Enhancing the
533
conversion of organic waste into biofertilizer with thermophilic bacteria. Environ. Eng.
534
Sci. 2012, 29: 726-730.
535
(32) Wei,Z.; Xi, B.; Zhao, Y.; Wang, S.; Liu, H.; Jiang, Y. Effect of inoculating
536
microbes in municipal solid waste composting on characteristics of humic
537
acid. Chemosphere 2007, 68: 368-374.
538
(33) Yang, S.; Xie, J.; Hu, N.; Liu, Y.; Zhang, J.; Ye, X.; Liu, Z. Bioconversion of
539
gibberellin fermentation residue into feed supplement and organic fertilizer employing
540
housefly (musca domestica l.) assisted by corynebacterium variabile. PloS One 2015, 10:
541
e0110809.
542
(34) Bouacem, K.; Bouanane-Darenfed, A.; Jaouadi, N.Z.; Joseph, M.; Hacene, H.;
543
Ollivier, B.; Fardeau, M.L.; Bejar, S.; Jaouadi, B. Novel serine keratinase from
544
caldicoprobacter algeriensis exhibiting outstanding hide dehairing abilities. Int. J. Biol.
545
Macromol. 2016, 86: 321-328.
546
(35) Liu, J.; Wan, L.I.; Xiu-Hong, X.U.; Hong-Tao, L.I. Effect of
547
cellulose-decomposing strain on microbial community of cow manure compost.
548
Environ. Sci. Technol. 2011, 32: 3073-3081.
549
(36) Sun, W.; Yu, G.; Louie, T.; Liu, T.; Zhu, C.; Xue, G.; Gao, P. From mesophilic to
26 ACS Paragon Plus Environment
Page 26 of 59
Page 27 of 59
Journal of Agricultural and Food Chemistry
550
thermophilic digestion: the transitions of anaerobic bacterial, archaeal, and fungal
551
community structures in sludge and manure samples. Appl. Microbiol. Biotechnol. 2015,
552
99: 10271-10282.
553
(37) Chen, Y.N.; Huang, J.X.; Li, Y.P.; Zeng, G.M.; Zhang, J.C.; Huang, A.Z.; Zhang, J.;
554
Ma, S.; Tan, X.B.; Xu, W.; Zhou, W. Study of the rice straw biodegradation in mixed
555
culture of Trichoderma viride and Aspergillus niger by GC-MS and FTIR. Environ Sci
556
Pollut Res. 2015, 22: 9807-9815.
557
(38) Cho, K.M.; Kwon, E.J.; Kim, S.K.; Kambiranda, D.M.; Math, R.K.; Lee, Y.H.;
558
Kim, J.; Yun, H.D.; Kim, H. Fungal diversity in composting process of pig manure and
559
mushroom cultural waste based on partial sequence of large subunit rRNA. J. Microbiol.
560
Biotechnol. 2009, 19: 743-748.
561
(39) Gwyther, C.L.; Williams, A.P.; Golyshin, P.N.; Edwards-Jones, G.; Jones, D.L.
562
The environmental and biosecurity characteristics of livestock carcass disposal methods:
563
A review. Waste Manage. 2011, 31: 767-778.
564
Figure captions
565
Figure 1. Schematic representation of the experimental setup (M=poultry manure,
566
S=sawdust)
567
Figure 2. Chicken degradation rate, temperature, pH value, moisture, organic matter
568
and C/N ratio at all sampled sites in three treatments (A, B, C, D, E and F). Different
569
lowercase letters indicate significant difference between days in each treatment (p