Modulation of the Fecal Microbiota in Sprague ... - ACS Publications

Modulation of the Fecal Microbiota in Sprague-Dawley Rats Using Genetically. Modified and Isogenic Corn Lines. Penggao Li †, ‡, Chun Yang†, ‡,...
0 downloads 0 Views 3MB Size
Subscriber access provided by EDINBURGH UNIVERSITY LIBRARY | @ http://www.lib.ed.ac.uk

Article

Modulation of the fecal microbiota in Sprague-Dawley rats using genetically modified and isogenic corn lines Penggao Li, Chun Yang, Rong Yue, Yaping Zhen, Qin Zhuo, Jianhua Piao, Xiaoguang Yang, and Rong Xiao J. Agric. Food Chem., Just Accepted Manuscript • DOI: 10.1021/acs.jafc.7b05285 • Publication Date (Web): 21 Dec 2017 Downloaded from http://pubs.acs.org on December 24, 2017

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 34

Journal of Agricultural and Food Chemistry

Modulation of the Fecal Microbiota in Sprague-Dawley Rats Using Genetically Modified and Isogenic Corn Lines

Penggao Li †, ‡, Chun Yang †, ‡, Rong Yue§, Yaping Zhen∥, Qin Zhuo⊥,*, Jianhua Piao⊥, Xiaoguang Yang⊥, Rong Xiao †, ‡,*



School of Public Health, Capital Medical University, Beijing 100069, People’s Republic of China ‡

Beijing Key Laboratory of Environmental Toxicology, Beijing 100069, People’s Republic of China §

Yuncheng Central Hospital, Yuncheng, Shanxi 044000, People’s Republic of China



Youanmen Clinical Detection Center, Capital Medical University, Beijing 100069, People’s Republic of China ⊥

Key Laboratory of Trace Element Nutrition NHFPC, Institute for Nutrition and Health, Chinese Center for Disease Control and Prevention, Beijing 100050, People’s Republic of China

*Correspondence: Q Zhuo, Tel: +86-10-66237240. Fax: +86-10-67711813. E-mail: [email protected] or [email protected], or R Xiao, Telephone/Fax: +86-10-83911512. E-mail: [email protected].

2

ACS Paragon Plus Environment

Journal of Agricultural and Food Chemistry

1

ABSTRACT:

2

This study investigated the composition and proportions of fecal microbiota in

3

Sprague-Dawley rats after consuming two genetically modified (GM) corn lines in

4

comparison with the isogenic corn and the AIN93G standard feed for 10 weeks using

5

bar-coded 16S rRNA gene sequencing. As a result, GM corn didn’t significantly alter

6

the overall health and alpha-diversity of fecal microbiota. Fecal microbiota structures

7

could be separated into non-corn and corn but not non-GM and GM corn subgroups.

8

Both non-GM and GM corn caused the increase in bacterial populations related to

9

carbohydrates utilization, such as Lactobacillus, Barnesiella and Bifidobacterium, and

10

the reduction in potentially pathogenic populations, such as Tannerella and

11

Moraxellaceae. In conclusion, similar effects on the fecal microbiota were observed

12

after consuming a GM- and non-GM-corn-based diet for long periods. Further studies

13

are warranted to elucidate the functional relevance of the changes in the proportions

14

of bacterial populations in these diets.

15

KEYWORDS: gut microbiota; genetically modified corn; diet; feces; rat

3

ACS Paragon Plus Environment

Page 2 of 34

Page 3 of 34

Journal of Agricultural and Food Chemistry

16

INTRODUCTION

17

Corn is a principal crop worldwide and is first among the grain crops in terms of

18

production. It is widely used, directly and indirectly, as human food and animal feed

19

in most countries. It is also an important raw material for various industries. As a

20

result, the demand for this grain is increasing. Over the past two decades, to improve

21

the production and quality of corn, genetic engineering with desired target genes has

22

been extensively employed to produce transgenic corn cultivars with improved traits,

23

and a great number of genetically modified (GM) corn lines have been

24

commercialized. However, concerns over the safety of GM food consumption still

25

exists.1 Although a scientific consensus that currently available food derived from

26

GM crops poses no greater risk to human health than conventional food, 2-3 each GM

27

food and its safety must be tested on a case-by-case basis before introduced to the

28

market.4

29

In recent years, it has been recognized that the gut microbiota plays an important

30

role in maintaining the health of the host. Changes in the structure of the bacterial

31

communities within the gastrointestinal tract have been implicated in many conditions,

32

such as inflammatory bowel diseases, type 2 diabetes and brain abnormalities.5

33

However, the composition of the gut microbiota is influenced by a myriad of

34

exogenous factors, such as diet, antibiotic use and exercise.6-8 Among these

35

exogenous factors, a long-term dietary pattern may play a significant role in shaping

36

the gut microbial community.9-11 Thus, several studies have investigated the influence 4

ACS Paragon Plus Environment

Journal of Agricultural and Food Chemistry

37

of GM corn on gut microbiota in different animal models and have provided valuable

38

insights into the effects of GM corn on microbial communities.12-17 However, in these

39

studies, the amount of the GM corn added to the animal feed was usually less than

40

70%. Most studies only compared one GM corn cultivar with its non-GM counterpart,

41

and some studies did not contain a blank control that consumed the standard feed.

42

Additionally, many studies investigated only the short-term effects of consumption,

43

which may not be informative because it may take six weeks of exposure for the

44

intestinal microbiota to adapt to the feed structure.18

45

Thus, in the present study, two GM corn cultivars, one herbicide-resistant cultivar

46

carrying the endolpyruvylshikimate-3-phosphate synthase (EPSPS) and the

47

phosphinothricin acetyl transferase (PAT) genes and one insect- and

48

herbicide-resistant cultivar carrying the insecticidal Cry1Ab gene from Bacillus

49

thuringiensis and EPSPS genes, in comparison with their common isogenic parental

50

line and the standard AIN93G rodent chow, were investigated simultaneously to

51

observe changes in the gut microbiota after 10 weeks of dietary intervention in

52

Sprague-Dawley (SD) rats. The dietary pattern in early life is vital for shaping the gut

53

microbial structure of an individual.19,20 Thus, weaning animals were used in the

54

present study. Moreover, because the main objective was to observe the effects of the

55

consumption of GM corn on the gut microbiota, as long as the nutritional needs of the

56

animals were met, the percentage of the corn in their diet was as high as possible.

57

After careful determination, calculation and comparison with the AIN93G standard

58

feed, approximately 75% of the feed composition was replaced by corn flour without 5

ACS Paragon Plus Environment

Page 4 of 34

Page 5 of 34

Journal of Agricultural and Food Chemistry

59

significantly changing the energy and nutrient composition of the feed. Thus, the

60

percentages of the three corn flours in their respective feeds were all greater than 75%

61

in the present study.

62

6

ACS Paragon Plus Environment

Journal of Agricultural and Food Chemistry

63

64

MATERIALS AND METHODS

Animals.

65

All animal work was performed according to the guidelines of the Institutional

66

Animal Care and Use Committee and a local veterinary office. The experimental

67

protocol was approved by, and a license (No. 20141011) obtained from, the

68

Experimental Animal Ethics Committee of the Institute for Nutrition and Food Safety,

69

Chinese Center for Disease Control and Prevention. All animal experiments were

70

strictly supervised to ensure that the animals were treated humanely and their needs

71

for habitat, food type and food supply mode, as well as their olfactory, entertainment,

72

training and social needs were met according to the relevant guidelines. Forty (20

73

female and 20 male) 4-week-old weaning Sprague-Dawley (SD) rats with body

74

weights of 60–80 g were purchased from Beijing HFK Bioscience Co., Ltd (Beijing,

75

China) and housed in a specific pathogen-free animal experiment center at the

76

Chinese Center for Disease Control and Prevention. The rats were housed at

77

24 °C±2 °C with 65%±5% humidity on a 12h light/dark cycle with free access to feed

78

and water. After one week of adaption, they were randomly and evenly assigned to

79

four groups to receive different feeds for 10 weeks. Five animals of the same group

80

and same sex were housed in one regular rat cage when they were young and then

81

separated to two cages when they grow older. All animals were monitored daily for

82

abnormalities and mortality. The abnormality signs include, but were not limited to,

83

changes in skin, fur, eyes and mucous membranes, as well as the occurrence of 7

ACS Paragon Plus Environment

Page 6 of 34

Page 7 of 34

Journal of Agricultural and Food Chemistry

84

secretions and excretions, and autonomic activities (e.g. unusual respiratory patterns).

85

Other abnormalities were behavioral, such as changes in gait, posture and response to

86

handling, as well as the presence of clonic or tonic movements, stereotypy (e.g.

87

excessive grooming and repetitive circling) or bizarre behaviors (e.g. self-mutilation

88

and walking backwards). The body weight of each animal was recorded once a week.

89

Average food consumption for each animal was determined twice a week by

90

subtracting the amount remaining from the initial amount provided for each cage and

91

then dividing the days and the number of the animals in the cage. After 10 weeks of

92

dietary intervention, fresh feces of each animal were collected in sterile tubes and

93

stored at −80 °C immediately.

94

Feeds.

95

The herbicide-resistant corn carrying the EPSPS and PAT genes (C0010.3.7), the

96

insect- and herbicide-resistant corn carrying the Cry1Ab and EPSPS genes

97

(C0030.3.5), and their parental non-GM counterpart (DBN318) were used in the

98

current work. The expression of the EPSPS protein imparts the tolerance to the

99

herbicide glyphosate. The expression of the PAT gene imparts the tolerance to the

100

herbicide phosphinotricin (glufosinate). The C0030.3.5 corn has been engineered to

101

express the insecticidal Cry1Ab protein from B. thuringiensis to resist crop damage

102

caused by insects. It was also transformed with the EPSPS gene to provide

103

herbicide-resistant capabilities. The two GM corn lines were produced by Beijing

104

DBN Biotech Co., Ltd. (Beijing China) using Agrobacterium-mediated 8

ACS Paragon Plus Environment

Journal of Agricultural and Food Chemistry

105

transformations. Both the GM corn lines and their parental corn were simultaneously

106

cultivated in adjoining plots in the experimental field of Yitong Manchu Autonomous

107

County, Siping City, Jilin Province, P. R. China. They were planted and raised under

108

identical conditions, and were harvested and stored in the same way. The presence of

109

the Cry1Ab, EPSPS and PAT transformation cassettes was confirmed by PCR using

110

standard protocols.

111

The nutritional compositions of the corn lines (supplementary Table S1) were

112

determined and compared with that of the AIN93G standard feed to determine which

113

ingredient was needed to make all of the feeds used in the current experiment meet the

114

specifications of the AIN93G standard feed, which was used as a negative control in

115

the present study. The formulas of the corn-based feeds are provided in supplementary

116

Table S2. The percentages of the flour powders from DBN318, C0010.3.7 and

117

C0030.3.5 corn in the feeds were 75.0%, 76.5% and 76.2%, respectively. After

118

mixing, all feeds were pelleted by the Beijing HFK Bioscience company (Beijing,

119

China), and then vacuum-packed and sterilized by gamma-irradiation. The final

120

compositions of the feeds were determined again to ensure that they met the

121

specifications of AIN93G as shown in supplementary Table S3.

122

Fecal DNA Extraction and PCR Amplification.

123

Genomic DNA was extracted from 0.18 to 0.22 g of stool using the E.Z.N.ATM

124

Mag-Bind Soil DNA Kit (Omega Bio-tek, Norcross, GA, USA) according to

125

manufacturer’s protocol. The V3–V4 regions of the bacterial 16S ribosomal RNA 9

ACS Paragon Plus Environment

Page 8 of 34

Page 9 of 34

Journal of Agricultural and Food Chemistry

126

genes were amplified by PCR (94 °C for 3 min, followed by 5 cycles at 94 °C for 30 s,

127

45 °C for 20 s and 65 °C for 30 s and another 20 cycles at 94 °C for 20 s, 55 °C for 20

128

s and 72 °C for 30 s, with a final extension at 72 °C for 10 min, using primers 341F

129

5′-CCCTACACGACGCTCTTCCGATCTG-barcode-CCTACGGGNGGCWGCAG-3′

130

and 805R

131

5′-GACTGGAGTTCCTTGGCACCCGAGAATTCCAGACTACHVGGGTATCTAAT

132

CC-3′. An 8-bp barcode sequence unique to each sample was attached to primers for

133

multiplexing. PCRs were performed in triplicate in 50-µL reaction mixtures

134

containing 5 µL of 10× PCR buffer, 0.5 µL of 10 mM dNTPs, 10 ng of genomic DNA,

135

0.5 µL of each primer (50 µM) and 0.5 µL of Plantium Taq polymerase (5 U/µL). The

136

PCR products were separated with 1.2% agarose gel electrophoresis and bands of the

137

desired size (> 600 bp) were purified using the SanPrep DNA Gel Extraction Kit

138

(Sangon Biotech, Shanghai, China) according to the manufacturer’s instructions.

139

Illumina MiSeq Sequencing.

140

Samples were sent to the Sangon Biotech Co., Ltd. (Shanghai, China) for

141

sequencing. The DNA concentration of each PCR product was determined using

142

Qubit2.0 DNA (Life Tech, CA, USA) before sequencing. Equimolar purified products

143

were pooled and paired-end sequenced (2 × 300) on an Illumina MiSeqPE300

144

platform according to the standard protocol at Sangon Biotech.

145

Bioinformatics Analysis.

10

ACS Paragon Plus Environment

Journal of Agricultural and Food Chemistry

146

Demultiplexing and quality-filtering of the raw fastq files were performed using

147

QIIME (version 1.8.0).21 An operational taxonomic unit (OTU)-based method was

148

performed to analyze where sequences were split into bins based on taxonomy. They

149

were clustered into each bin using a cutoff point of 0.05. UPARSE was used to cluster

150

OTUs with a 97% similarity cutoff level, and UCHIME was applied to identify and

151

remove chimeric sequences. The phylogenetic affiliation analysis of each 16S rRNA

152

gene sequence was introduced by RDP Classifier against the SILVA (SSU115) 16S

153

rRNA database with a confidence threshold of 70%.22 An alpha-diversity analysis was

154

applied to analyze the complexity of the species’ diversity through several indices,

155

including Richness, Shannon, Coverage and Simpson. Estimators of community

156

richness, evenness and diversity were calculated based on OTUs (97% similarity).

157

The overall microbiota’ dissimilarities among all of the samples were accessed using

158

the weighted principal coordinates analysis by the R package (R 3.0.2). The heatmap

159

construction and clustering analysis were also performed using the R package. Data

160

are expressed as mean ± standard deviation. A statistical analysis of differences

161

among groups was performed by one-way ANOVA using SPSS Statistics 18.0 (IBM,

162

New York, USA), and P < 0.05 was considered statistically significant.

11

ACS Paragon Plus Environment

Page 10 of 34

Page 11 of 34

Journal of Agricultural and Food Chemistry

163

RESULTS

164

Health and Behavioral Features.

165

During the experiment, rats in all of the groups were in good general health with no

166

abnormalities in skin, fur, eyes or mucous membranes, no abnormities in secretion,

167

excretion or autonomic activities, no changes in gait, posture or response to handling,

168

and no other abnormal behaviors, such as clonic or tonic movements, excessive

169

grooming, repetitive circling, self-mutilation, or walking backwards. Figure 1 shows

170

that the daily feed intakes of the rats eating corn-based feeds were relatively higher

171

than those of rats eating the standard chow along the intervention; these differences

172

were more evident in females than in males. Generally, the rats eating the non-GM

173

corn had the highest intakes and the rats eating the standard chow had the lowest

174

intakes, with the rats eating the two GM corn lines between them. In spite of this, the

175

body weights of the rats from the four feeding groups were not significantly different

176

from each other throughout the experimental period, indicating that their growth was

177

not affected by the differences in the feeds.

178

Microbial Richness and Diversity Differences Owing to Diet.

179

The Illumina MiSeq PE300 high-throughput 16S rRNA sequencing provided 361,351

180

usable sequence reads, with a mean length of 413.2 ± 10.3 bp, that were obtained

181

from 20 fecal samples from female rats, corresponding to 686 OTUs, with an average

182

of 329 ± 61 per biological sample at a similarity level of 97%. Each sample was 12

ACS Paragon Plus Environment

Journal of Agricultural and Food Chemistry

183

covered by an average of 18,068 reads. For male rats, 441,844 usable sequence reads,

184

with a mean length of 413.1 ± 12.3 bp, were obtained from 20 fecal samples,

185

corresponding to 634 OTUs, with an average of 333 ± 39 per biological sample at a

186

similarity level of 97%. Each sample was covered by an average of 22,092 reads. The

187

raw sequencing data have been uploaded to the NCBI SRA database (accession No.

188

SRX3299851). The coverage indices for all of the samples were above 0.99,

189

indicating that the depth of sequencing was sufficient. The flattening of the rarefaction

190

curve based on the values of species richness (OTU number) and Shannon indices

191

(supplementary Figure S1) indicated that our data volume covered all species in the

192

fecal samples’ communities. Figure 2 shows that the alpha-diversity indices, including

193

the species richness, the Shannon and the Simpson indices, were slightly higher in the

194

corn-fed groups but were not significantly different among the four dietary groups (P >

195

0.05). This indicates that the microbial richness and diversity in the feces of the

196

GM-corn-fed rats did not significantly change in comparison with the non-GM-corn

197

and the standard-chow-fed rats.

198

Effects of Dietary GM and non-GM Corn Intake on Gut Microbial Community

199

Structure.

200

The weighted principal coordinate analysis plot in Figure 3A shows that in female

201

animals, 43.4%, 23.6% and 8.2% of variation could be explained by PC1, PC2 and

202

PC3, respectively, and most samples from corn-fed rats had a negative PC1 value,

203

while most samples from the standard feed group had a positive PC1 value. Figure 3B 13

ACS Paragon Plus Environment

Page 12 of 34

Page 13 of 34

Journal of Agricultural and Food Chemistry

204

shows that in male animals, 48.0%, 18.4% and 9.8% of variation could be explained

205

by PC1, PC2 and PC3, respectively, and most samples from corn-fed rats had a

206

negative PC2 value, while most samples from the standard feed group had a positive

207

PC2 value. Samples from the GM- and non-GM-corn-fed rats were interspersed

208

among each other with no clear distribution pattern to separate them from each other,

209

suggesting that the microbiota structures of the rats consuming the GM corns were not

210

significantly different from that of the non-GM-corn-fed rats. The similarities among

211

the fecal samples were further analyzed by the hierarchically clustered heatmap

212

analysis (Figure 3, right panel), which showed that samples from the GM- and

213

non-GM-corn-fed rats could not be clearly clustered into separate subgroups and that

214

the composition and relative abundances of the predominant gut bacterial genera were

215

very similar among individuals from the different feeding groups.

216

Fecal Bacterial Taxonomic Compositions in Rats Consuming Different Feeds.

217

Although it was hard to clearly distinguish the microbiota community structures of the

218

GM- and non-GM-corn-fed rats, there were differences in the relative abundance

219

levels of the fecal bacteria. An outline of the relative abundance levels of the bacterial

220

populations at different taxonomic levels are presented in supplementary Table S4–S8

221

to show the influence of the different diets. As shown in Figure 4, distribution

222

histograms of the taxonomic composition of each sample from the respective groups

223

at the phylum and genus level were constructed. 9 and 10 different bacterial phyla

224

were identified in male and female rats, respectively. In accordance with previous 14

ACS Paragon Plus Environment

Journal of Agricultural and Food Chemistry

225

reports, Firmicutes, Bacteroidetes, Actinobacteria, Proteobacteria and

226

Verrucomicrobia were the major phyla in all of the samples. Lactobacillus,

227

Barnesiella, Ruminococcus, Bacteroides, Clostridium XI and Clostridium XlVa were

228

the predominant genera in all of the samples. Figure 4 also shows that the different

229

bacterial genera were more evenly distributed in the feces of the rats eating the

230

standard feed than in those eating a corn-based diet. The dominant bacteria appeared

231

to become more predominant after dietary intervention with corn, indicating that

232

consuming a large amount of corn for a long period of time might facilitate the growth

233

of the already dominant bacterial populations in the gut. And, compared with the

234

standard feed, the similarities of the gut microbiota patterns among the three different

235

groups of corn-fed rats could be partly attributed to similar changes occurring in the

236

abundance levels of the main gut bacteria in the three groups of rats fed corn. The

237

Firmicutes/Bacteroidetes ratio also showed a tendency to decrease in both sexes for

238

all of the corn-fed groups, although differences were not statistically significant owing

239

to the high inter-individual variances (supplementary Tables S4).

240

The relative proportion of phyla and genera as depending on the types of the diets

241

differed between males and females. As shown in Figure 5, in female animals, the

242

abundance levels of the major phyla were not significantly different among the four

243

feeding groups. However, in male animals, in comparison with the standard feed

244

group, the proportion of Firmicutes in the C0030.3.5 group decreased significantly

245

(59.59 vs. 80.25%, P < 0.05), and the proportions of Verrucomicrobia and Candidatus

246

Saccharibacteria increased significantly in the DBN318 group (5.74% and 1.85% vs. 15

ACS Paragon Plus Environment

Page 14 of 34

Page 15 of 34

Journal of Agricultural and Food Chemistry

247

0.32% and 0.19%, respectively, P < 0.05).

248

At the genus level, 119 and 120 different bacterial populations were identified in

249

the female and male animals, respectively. The proportions of 29 and 28 genera were

250

significantly different among the four dietary groups in female and male rats,

251

respectively (supplementary Tables S8). Compared with the standard feed,

252

Lactobacillus, the largest genus in all of the fecal samples, showed an increase after

253

corn-feeding (Figure 5), which was significant in the C0030.3.5 group of female rats

254

(36.13 vs. 12.02%, P