Effect of Substrate on Identification of Microbial Communities in Poultry

Aug 22, 2016 - ... abundance of microorganisms associated with the decomposition of poultry carcasses. Understanding the potential composting driver c...
0 downloads 0 Views 4MB Size
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