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