Subscriber access provided by UNIV OF DURHAM
New Analytical Methods
Development and Inter-laboratories Validation of Event-Specific Quantitative Real-Time PCR Method for Genetically Modified Rice G6H1 Event Litao Yang, Yu Yang, Wujun Jin, Xiujie Zhang, Xiaying Li, Yuhua Wu, Jun Li, and Liang Li J. Agric. Food Chem., Just Accepted Manuscript • DOI: 10.1021/acs.jafc.8b01519 • Publication Date (Web): 09 Jul 2018 Downloaded from http://pubs.acs.org on July 9, 2018
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 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 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.
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 29
Journal of Agricultural and Food Chemistry
1
Development and Inter-laboratories Validation of Event-Specific
2
Quantitative Real-Time PCR Method for Genetically Modified Rice
3
G6H1 Event
4 5
Litao Yang1*, Yu Yang1, Wujun Jin2, Xiujie Zhang3, Xiaying Li3, Yuhua Wu4, Jun Li4,
6
Liang Li2*
7
1
8
Organisms, School of Life Sciences and Biotechnology, Shanghai Jiao Tong University,
9
Shanghai 200240, China
National Center for the Molecular Characterization of Genetically Modified
10
2
11
100081, China
12
3
13
Republic of China. Beijing 100025, China.
14
4
15
Institute, Chinese Academy of Agricultural Sciences, No. 2 Xudong 2nd Road, Wuhan
16
430062, People’s Republic of China
Biotechnology Research Institute, Chinese Academy of Agricultural Sciences, Beijing
Development Center of Science and Technology, Ministry of Agriculture of People’s
Key Laboratory of Oil Crop Biology of the Ministry of Agriculture, Oil Crops Research
17 18
*
19
[email protected] (Liang Li).
To whom correspondence should be addressed:
[email protected] (Litao Yang) or
ACS Paragon Plus Environment
Journal of Agricultural and Food Chemistry
20
Abstract
21
The transgenic rice G6H1 was a new event with the traits of herbicide-tolerance and
22
insect-resistant. Herein, we developed one event-specific real-time PCR method with
23
high specificity and sensitivity for G6H1 event quantitative analysis, and validated its
24
performance on practical samples quantification through a collaborative ring trial. A
25
total of eight laboratories participated in this validation and quantified three blind
26
G6H1 powder samples including DNA extraction and real-time PCR analysis. The
27
statistical analyzed results from returned data confirmed its high PCR efficiency, and
28
good linearity, trueness and precision, indicating that the developed G6H1 real-time
29
PCR assay were accurate, reliable, and comparable for G6H1 identification and
30
quantification.
31
Keywords: Collaborative ring trial; transgenic rice G6H1; Real-time PCR
ACS Paragon Plus Environment
Page 2 of 29
Page 3 of 29
Journal of Agricultural and Food Chemistry
32
Introduction
33
Rice is a staple food enriched with protein, vitamins and minerals, and which is
34
widely planted in the world, especially in Southeast Asia and East Asia. However,
35
some problems about rice planting that exist and disturb long are infestation of
36
vermin and the use of herbicide, which bring huge damage to rice harvest and
37
environment pollution.1 In recent years, DNA recombinant technique is widely
38
applied in agriculture and enables continuous increase of genetically modified
39
organisms (GMOs). Update to the end of 2016, 185.1 million hectares of GM crops
40
have been planted in global area, more than 450 genetically modified (GM) events
41
from 29 crop species have been developed and approved for commercialization.2 In
42
GMOs, exogenous genes are integrated into recipient genome, and which can create
43
one or more traits, such as insect-resistant and herbicide-tolerant.3 With great
44
progress in rice genetic modification, more and more GM rice events with different
45
traits have been developed.1, 4-5 Eight GM rice events have been approved for
46
commercialization in several countries, such as PWC16, CL121/CL141/CFX51 and
47
IMINTA-1/IMINTA-4 in Canada, LLRICE06/LLRICE62 in Australia, Canada, Colombia,
48
Mexico, and the USA. LLRICE601 in Colombia, and Huahui No.1/Bt Shanyou 63 in
49
China, Tarom molaii +Cry1Ab in Iran and 7 crp#10 in Japan.2, 6-7 In China, the
50
transgenic rice research with the purpose of high insect-resistance,
51
herbicide-tolerance, yields and quality is one hot topic because rice is one primary
52
alimentary crop.1 Several transgenic rice events with the trait of insect-resistant have
53
been developed in China, such as TT51-1, Kemingdao, KF6, T1C-19, and T2A-1, etc. 1, 5
54
GM rice G6H1 developed by Zhejiang University is a novel event, which was
55
produced by inserting a modified Cry1Ab/Vip3H gene encoding insect-resistant HJC-1
ACS Paragon Plus Environment
Journal of Agricultural and Food Chemistry
56
protein and a G6-EPSPS gene encoding herbicide-tolerant protein. 8-10 The G6H1 rice
57
showed very high insect-resistance and herbicide tolerance in the phases of field test
58
and environment, and which are in the pipeline of applying the safety production
59
certificate.8
60
However, with the large-scale applications of transgenic technology in modern
61
agricultural production, several controversial issues are being discussed including
62
food safety, environment risk and ethical concerns.3-4 Several countries and areas
63
have issued labeling regulations to protect the authority of consumers.11 For the
64
implement of GMO labeling regulations, rapid and accurate detection methodologies
65
of GM content have been developed.11 DNA-based PCR techniques were mainly used
66
in GMO identification and quantification, such as PCR,12 real-time PCR,13-16 and digital
67
PCR, etc.17-19 Real-time PCR employing the TaqMan probe technique is the golden
68
standard for GMO analysis because of high sensitivity, accuracy, and practicability. A
69
lot of real-time PCR assays were developed and validated to detect and quantify GM
70
contents in raw materials and processed food/feed samples.20-25 Also, the
71
collaborative ring trial validation is one prerequisite step for proposing a method to
72
be used in routine lab analysis and included as one national or ISO standard.26-27
73
According to the requests of GMO regulations in China and EU, the validated
74
detection method for one GM event should be supplied before its application for
75
commercialization.11 In this study, one event-specific real-time PCR method for G6H1
76
rice inspection and monitoring were developed based on the 3’ junction sequence
77
between exogenous DNA and rice DNA, and well validated employing three blind
78
samples by a collaborative ring trial.
ACS Paragon Plus Environment
Page 4 of 29
Page 5 of 29
Journal of Agricultural and Food Chemistry
79
Material and Methods
80
Materials and DNA extraction
81
Homozygous seeds of G6H1 rice and its recipient material Xiushui110 were grown
82
and supplied by Zhejiang Academy of Agricultural Science (ZJAAS). The seed samples
83
of GM rice TT51-1, T1C-19, Kemingdao, and KF6 were supplied by their developers.
84
The mixed seed powder samples of GM maize (MON810, MON863, and Bt176), GM
85
cotton (MON531 and MON88913), GM soybean (GTS40-3-2 and MON89788), and
86
GM canola (Ms1Rsf1, Ms3Rf83, and GT73) were prepared by our laboratory. Three
87
G6H1 matrix-based certified reference material (CRM) samples with the GM
88
mass/mass ratio of 4.983%, 0.997%, and 0.499% were produced by Shanghai Jiao
89
Tong University (SJTU), China and used as blind samples in collaborative ring trial.28
90
All plant genomic DNAs were extracted firstly by CTAB extraction according to the
91
standard of ISO 21571:2005 and purified using DNeasy 96 Plant Kit (Qiagen, Hilden,
92
Germany).29 The quality and concentration of purified DNAs were measured and
93
estimated using 1.0% agarose gel electrophoresis and the ultraviolet spectrometric
94
method with a NanoDrop 1000 UV/vis spectrophotometer (NanoDrop, Wilmington,
95
DE, USA).
96 97
Oligonucleotide primers and probes
98
The event-specific primers (G6H1-F/G6H1-R) and probe (qG6H1-p) were designed
99
based on the 3’ end junction sequence of transgene insertion for G6H1 content
100
quantification. Previously reported primer pair of SPS-F/SPS-R and SPS-P probe were
101
used for quantification of rice endogenous gene saccharose phosphate syntheatase
102
(SPS).21 All primers and probes were synthesized by Invitrogen Co., Ltd. Shanghai,
ACS Paragon Plus Environment
Journal of Agricultural and Food Chemistry
103
China), and the sequence information was listed in Table 1 and Figure 1.
104 105
Quantitative Real-time PCR
106
Real-time PCR assays for the G6H1 event and SPS gene were optimized with a 25-uL
107
reaction volume, containing 5 μL genomic DNA and 20 μL Real-time PCR Reaction
108
Master Mix. The Real-time PCR Reaction Master Mix was consisted of the following
109
reagents: 1 X PCR buffer, 200 μM dNTPs, 400 nM primers, 200 nM TaqMan probes,
110
1.25 U Taq DNA polymerase (Huirui Biotech company, Shanghai, China). All real-time
111
PCR reactions were carried out with the following program: 95°C for 10min, and then
112
45 cycles of 94 °C for 15s and 60 °C for 60s. The fluorescent signal was monitored
113
during each PCR cycle at the annealing step. Each quantitative PCR was repeated
114
three times and each time with three replicate reactions.
115 116
Collaborative Ring Trial
117
The collaborative ring trial was organized by Chinese Academy of Agricultural
118
Sciences (CCAAS). A total of 8 laboratories were invited and participated in the ring
119
trial. Each laboratory received a kit for the validation study, including one tube of
120
homozygous G6H1 genomic DNAs named with G1 (50 ng/μL, 100 μL), three bottle of
121
G6H1 CRM powder samples (coded with S1, S2, and S3, respectively), one tube of 0.1
122
X TE dilution buffer (5 mL in volume) and pre-mixed Real-time PCR Reaction Master
123
Mix of G6H1 and SPS (2 tubes for each assay, and each tube with the volume of 2.0
124
mL). Kits were shipped in dry ice to all laboratories. All participating laboratories also
125
received a detailed protocol to follow and an excel file for reporting results. All
126
laboratories were requested to report any deviation from the protocol which may
ACS Paragon Plus Environment
Page 6 of 29
Page 7 of 29
Journal of Agricultural and Food Chemistry
127
have occurred during the preparation and execution of the real-time PCR
128
experiments. The DNA extraction kit was not supplied, DNeasy 96 Plant Kit (Qiagen,
129
Hilden, Germany) were recommend but no limited to each lab for blind samples DNA
130
extraction. For each lab, the G6H1 genomic DNAs of G1 was requested to dilute to a
131
series of ten-fold dilutions with the concentrations of 10.0, 1.0, 0.1, 0.01, and 0.001
132
ng/μL (named from R1 to R5) using 0.1 X TE dilution buffer, and which were used to
133
construct standard curves of G6H1 and SPS assays. For blind samples quantification,
134
each lab was requested to randomly sample three sub-samples from each sample for
135
DNA extraction and further quantitative analysis. The GM G6H1 content of each blind
136
sample was determined by the average values of three sub-samples. The extracted
137
DNAs of blind samples should be diluted to the final concentration of 5 ng/μL for
138
real-time PCR amplification. During each PCR reaction, the negative and blank control
139
reactions should be performed.
140 141
Data Analysis
142
In this ring trial, all participants were requested to return their results sheet within
143
two months after receiving the trial package. Based on the returned results, PCR
144
efficiency, linearity of regression, precision, and trueness were analyzed and
145
calculated according to ISO 5725-2 and definition of minimum performance
146
requirements for analytical methods of GMO testing reported by European Network
147
of GMO Laboratories (ENGL).27, 30 The identified outlying results and/or laboratories
148
were excluded from the analysis, in accordance with ISO 5725-2.30 The measurement
149
uncertainty of G6H1 method was estimated according to the guidance document
150
issued by the European Commission Joint research Center, Institute for Reference
ACS Paragon Plus Environment
Journal of Agricultural and Food Chemistry
151
Materials and Measurements (JRC-IRMM).31 Statistical analysis of result differences
152
among participating laboratories, such as least significant difference (LSD) test or the
153
Q test, were performed using SPSS 13.0 software.
ACS Paragon Plus Environment
Page 8 of 29
Page 9 of 29
Journal of Agricultural and Food Chemistry
154
Results and Discussion
155
Establishment of G6H1 event-specific real-time PCR method
156
For developing the event-specific method for G6H1, four sets of primers and probe
157
were designed and pre-evaluated according to their specificity and sensitivity. The
158
pre-evaluated results showed that the real-time PCR assay with primer pair (G6H1-F/
159
G6H1-R) and probe (G6H1-P) presented one 90-bp amplicon in length could be used
160
for further establishing the event-specific method of G6H1 (Data not shown). The
161
sequences and locations of primer pair (G6H1-F/ G6H1-R) and probe (G6H1-P)with
162
were shown in Figure 1.
163
Specificity test. The specificity of G6H1 real-time PCR method was tested among
164
different GM rice events (TT51-1, T1C-19, Kemingdao, and KF6) and GM plant species
165
(maize, soybean, canola, and cotton). The amplified results showed that only positive
166
fluorescent signal and traditional amplified curves were observed in the reactions
167
with G6H1 DNA template, no fluorescent signal and amplified curves were observed
168
in the reactions of mixed GM maize of MON810, MON863, and Bt176, mixed GM
169
cotton MON531 and MON88913, mixed GM soybean of GTS40-3-2, and MON89788,
170
mixed GM canola of Ms1Rsf1, Ms3Rf83, and GT73, non-GM rice Xiushui 110, and no
171
template control (NTC) (As shown in Figure 2a). The results revealed the high
172
specificity of developed G6H1 event-specific assay among different GM events.
173
Construction of Standard Curve. To quantify the GM content of unknown samples,
174
the standard curve of G6H1 real-time PCR assay was constructed employing five
175
concentrations of diluted G6H1 genomic DNA (50.0, 5.0, 0.5, 0.05, and 0.005
176
ng/reaction). High PCR efficiency of 100.95% was obtained (Figure 2b). The linearity
177
between the logarithm of DNA copy numbers and Ct values was evaluated by R2
ACS Paragon Plus Environment
Journal of Agricultural and Food Chemistry
178
values. The R2 value of 0.9999 was obtained (Figure 2b). The high PCR efficiencies
179
and good linearity indicated that the standard curve is suitable for further
180
quantitative measurement of GM rice with a large dynamic range.
181
Sensitivity evaluation. To determine the lowest amount of the initial template DNA
182
copies that could be reliably detected and quantified with 95% confidence level (i.e.,
183
absolute limit of detection, LOD; limit of quantification, LOQ), four concentrations of
184
genomic DNA samples (100, 50, 10, and 5 copies/reaction) were prepared and tested.
185
Each sample were tested with three repeats, and each repeat with 20 parallel
186
reactions. The expected positive fluorescent signals were obtained from all four
187
concentrations in all replicated reactions. However, a high relative SD (RSD) value
188
(49.81%) of genomic DNA amount was observed when using five haploid genomic
189
DNA copies as template (Table 2). Therefore, the LOD value of five haploid genomic
190
DNA copies was determined. The bias values of quantified results increased as the
191
DNA amounts decrease, and a large bias of 38.54% was observed in the reactions of
192
five haploid genomic DNA copies. To be quantified reliably, the approximately 10
193
initial copies were required, and the LOQ was determined with 10 copies of haploid
194
genome.
195 196
Validation of G6H1 event-specific real-time PCR method
197
In order to transfer the developed G6H1 method to other labs, one collaborative ring
198
trial of the G6H1 real-time PCR assay was designed and organized according to the
199
ISO 5725:1994. 26, 29 Since the SPS real-time PCR assay had been already validated for
200
rice genomic DNA detection and quantification,21 we focused more on the validation
201
of G6H1 real-time PCR assay in this ring trial, such as the PCR efficiency, linearity,
ACS Paragon Plus Environment
Page 10 of 29
Page 11 of 29
Journal of Agricultural and Food Chemistry
202
precision, trueness, and measurement uncertainty.
203
PCR efficiency and linearity. All participants were requested to prepare the calibrate
204
DNA solutions (R1-R5) by gradient dilution from sample G1, and run the quantitative
205
real-time PCR. Each reaction should be performed with three replicates and three
206
repeats. The returned Ct values of three repeats were plotted against the log
207
transformation of the copy number of the diluted solutions to obtain the standard
208
curves, respectively (Figure 3). The PCR efficiency, linearity of regression, and slope
209
were calculated from the standard curve to examine the applicability of quantitative
210
real-time PCR assays.27 All the data showed high PCR efficiency, good linearity, and
211
acceptable slope value in the curves for the G6H1 and SPS real-time PCR assays
212
(Table 3). The key parameters of PCR efficiency, slope and the regression coefficients
213
(R2) for the G6H1 and SPS assays from 8 laboratories were all in the acceptable range
214
based on the requirements proposed by the European Network of GMO Laboratories
215
(ENGL) guidelines.27 The linearity from each lab was higher than 0.9826, and the
216
mean square regression correlations (R2) were 0.9989 and 0.9991 for both G6H1 and
217
SPS real-time PCR assays, respectively. Meanwhile, the values of slope for the G6H1
218
assay and the SPS assay from all labs were ranged from -3.24 to -3.58. The average
219
values of slope for G6H1 assay (-3.43) and SPS assay (-3.40) were very similar,
220
indicating no difference existed between these two assays in quantification even the
221
target amplicons were completely different. The PCR efficiency was calculated using
222
the following equation E= (10(-1/slope) -1) X 100%. As is shown in Table 3, the mean
223
PCR efficiency was 95.55% for the G6H1 assay and 97.13% for the SPS assay. The
224
results demonstrated good linearity and high PCR efficiency of G6H1 and SPS
225
real-time PCR assays.
ACS Paragon Plus Environment
Journal of Agricultural and Food Chemistry
226
Quantification of blind samples. In order to evaluate the accuracy of the G6H1
227
real-time PCR method, three parameters including the trueness, precision, and
228
measure uncertainty were calculated through the verification. Three G6H1 CRMs
229
with the property values of 4.983%, 0.997%, and 0.499% (mass/mass ratio) were
230
used as blind samples with the code of S1, S2, and S3 in this validation, respectively.
231
In each participated laboratory, the S1-S3 blind samples was quantified according to
232
the constructed standard curves, and the G6H1 content were calculated based on the
233
formula of GM%= (G6H1 gene copy number/ SPS gene copy number) X100%. Each
234
lab reported three values of G6H1 contents for each sample, and which were
235
quantified from three sub-samples individually. The statistical data about GM content
236
of three samples was exhibited in Table 4. The calculated values of GM content were
237
checked with Dixon’s test and Grubbs’ test to find out remarkable outliers according
238
to IS0 5725.29 No outliers was discovered from both the test.
239
Trueness. The calculated GM G6H1 content of each sample was presented in Table 5.
240
No significant outliers were excluded from all the returned data. The average test
241
values of GM content for S1, S2, and S3 were 5.03%, 1.09%, and 0.53%, respectively.
242
The bias between the quantified values and true values of three samples from all
243
eight laboratories were ranged from 0.74% to 20.24% and shown in Figure 4. The
244
average bias between the quantified values and true values of three samples among
245
eight laboratories was within the range from 0.74% to 9.33%. In terms of real-time
246
PCR method acceptance standard enacted by ENGL (within±25%), the G6H1
247
event-specific method showed an ideal and reliable trueness.
248
Precision. In accordance with ISO5725-2, the precision of a test method is mainly
249
evaluated according to its reproducibility and repeatability. Both repeatability and
ACS Paragon Plus Environment
Page 12 of 29
Page 13 of 29
Journal of Agricultural and Food Chemistry
250
reproducibility are usually reported as relative standard deviation. For instance, the
251
relative repeatability (RSDr) describing the intra-laboratory variation, and relative
252
reproducibility (RSDR) describing the inter-laboratory variation. The repeatability and
253
reproducibility of the quantification result of each sample were estimated and listed
254
in Table 4. The RSDr values ranged from 3.91% to 6.33% and RSDR values ranged from
255
5.68% to 7.89% in blind samples analysis. All of these values were below the method
256
acceptance criteria value of 25%, indicating that the G6H1 real-time PCR method was
257
very reliable and repeatable even by different operators and equipment in different
258
labs.
259
Measurement uncertainty of quantified results. The measurement uncertainty (MU)
260
of G6H1 method in blind samples quantification was also evaluated. Based on the
261
guideline of measurement uncertainty for GMO testing laboratories, 30 the MU is
262
calculated using the formulas of = + ( × ) and = 2 × . Herein,
263
uc is the standard uncertainty, u0 is absolute standard uncertainty, RSU is relative
264
standard uncertainty, c is the mean concentration of tested sample, and U is
265
expanded uncertainty. The measurement uncertainty of the tested results is usually
266
reported as an expanded uncertainty (U) under the 95% confidence level. The u0 and
267
RSU was calculated by plotting the reproducibility standard deviation values (SDR)
268
against the mean concentration (c) of the blind samples in the ring trial (Figure 5). u0
269
is same with the absolute value of intercept of the linear regression (u0= 0.000044),
270
RSU value is equal to the slope of the linear regression (RSU=0.056728). For the three
271
tested blind samples, the expanded uncertainty (U) were calculated with the value of
272
0.57%, 0.12% and 0.06% for S1, S2 and S3 sample, respectively. As a result, the
273
measurement concentrations of the three samples were (5.03±0.57)% for S1,
ACS Paragon Plus Environment
Journal of Agricultural and Food Chemistry
274
(1.09±0.12)% for S2, and (0.53±0.06)% for S3, respectively.
275 276
Conclusion
277
In summary, the G6H1 event-specific real-time PCR method was successfully
278
established for accurate quantification of G6H1 contents with high specificity. The
279
LOD and LOQ was determined with the value of 5 and 10 haploid genomic DNA
280
copies, respectively. For the purpose of standardization, the developed G6H1
281
event-specific real-time PCR method was well validated through a collaborative ring
282
trial of 8 laboratories. The PCR efficiency, wide quantitative dynamic range, and good
283
PCR performance in blind samples quantification indicated that the developed
284
real-time PCR method is suitable for application in identification and quantification of
285
GM rice event G6H1 and its derived products, and could be used in routine lab
286
analysis in the future.
ACS Paragon Plus Environment
Page 14 of 29
Page 15 of 29
Journal of Agricultural and Food Chemistry
287
Acknowledgements
288
We heartily thank the following collaborators for their participation in this ring trial:
289
GMOs and Derived Products Inspection and Supervision Center (Tianjin), Ministry of
290
Agriculture, P.R. China ; Inspection Test Center for Environmental Safety of Transgenic
291
Crops, (Jinan), Ministry of Agriculture, P.R. China; Inspection Test Center for
292
Environmental Safety of Transgenic Crops, (Guangzhou), China; Inspection Test
293
Center for Environmental Safety of Transgenic Crops, (Anyang), China; Chinese
294
Academy of Inspection and Quarantine, China; National Institute of Metrology, China;
295
Shanghai Entry-Exit Inspection and Quarantine Bureau, China; Cereal Quality
296
Supervision and Testing Center, (Haerbin), Ministry of Agriculture, P.R. China. This
297
work was supported by the National Transgenic Plant Special Fund
298
(2016ZX08012-003 and 2016ZX08012-002), and the National Natural Science
299
Foundation of China (31471670).
ACS Paragon Plus Environment
Journal of Agricultural and Food Chemistry
300
Reference
301
1. Lu, B. R.; Qiang, F.; Shen, Z. Commercialization of transgenic rice in China: potential
302 303 304
environmental biosafety issues. Syst. Biodivers. 2008, 16(5), 426-436. 2. James, C. Global status of commercialized biotech/GM crops: 2016. ISAAA Briefs, 2016, No. 52.
305
3. Fraiture, M. A.; Herman, P.; De, Loose. M.; Debode, F.; Roosens, N. H. How can we
306
better detect unauthorized GMOs in food and feed chains? Trends Biotechnol.
307
2017, 35(6), 508-517.
308 309
4. Bajaj, S.; Mohanty, A. Recent advances in rice biotechnology towards genetically superior transgenic rice. Plant Biotechnol. J. 2005, 3, 275-307.
310
5. Lu, C. The first approved transgenic rice in China. GM Crops. 2010, 1(3), 113-115.
311
6. www.isaaa.org/gmapprovaldatabase/crop/default.asp?CropID=17&Crop=Rice.
312
7. http://www.cera-gmc.org/GMCropDatabase.
313
8. Chen, Y.; Tian, JC.; Shen, ZC.; Peng, YF.; Hu, C.; Guo, YY. Transgenic rice plants
314
expressing a fused protein of Cry1Ab/Vip3H has resistance to rice stem borers
315
under laboratory and field conditions. J. Econ. Entomol. 2010, 103(4), 1444-1453.
316
9. Tian, JC.; Liu, ZC.; Chen, M.; Chen, Y.; Chen, XX.; Peng, YF. Laboratory and field
317
assessments of prey-mediated effects of transgenic Bt rice on Ummeliata
318
insecticeps (Araneida: Linyphiidae). Environ. Entomol. 2010, 39(4), 1369-1377.
319
10. Lu, ZB.; Han, NS.; Tian, JC.; Peng, YF.; Cui, HU.; Guo, YY. Transgenic
320
cry1Ab/vip3H+epsps rice with insect and herbicide resistance acted no adverse
321
impacts on the population growth of a non-target herbivore, the white-backed
322
planthopper, under laboratory and field conditions. J. Integr. Agric. 2014, 13(12),
323
2678-2689.
ACS Paragon Plus Environment
Page 16 of 29
Page 17 of 29
Journal of Agricultural and Food Chemistry
324 325
11. Zhang, DB.; Guo, JC. The development and standardization of testing methods for GMO and their derived products. J. Integr. Plant Biol. 2011, 3, 539–551.
326
12. Yang, C.; Zhang, DB.; Yang, LT. Development of event-specific PCR detection
327
methods for genetically modified tomato Huafan No. 1. J. Sci. Food Agric. 2013,
328
93(3), 652-660.
329
13. Jiang, Y.; Yang, H.; Quan, S.; Liu, YN.; Shen, P.; Yang, LT. Development of certified
330
matrix-based reference material of genetically modified rice event TT51-1 for
331
real-time PCR quantification. Anal. Bioanal .Chem. 2015, 407(22), 6731-6739.
332
14. Yang, LT.; Guo, JC.; Pan, AH.; Zhang, HB.; Zhang, KW.; Wang, ZM.; Zhang, DB.
333
Event-specific quantitative detection of nine genetically modified maizes using
334
one novel standard reference molecule. J. Agric. Food Chem. 2007, 55(1), 15-24.
335
15. Li, ZQ.; Li, X.; Wang, CH.; Song, GW.; Pi, LQ.; Zheng, L.; Zhang, DB.; Yang, LT. One
336
novel multiple-target plasmid reference molecule targeting eight genetically
337
modified canola events for genetically modified canola detection. J. Agric. Food
338
Chem. 2017, 65(38), 8489-8500
339
16. Pi, LQ.; Li, X.; Cao, YW.; Wang, CH.; Pan, LW.; Yang, LT. Development and
340
application of a multi-targeting reference plasmid as calibrator for analysis of five
341
genetically modified soybean events. Anal. Bioanal. Chem. 2015, 407, 2877-2886.
342
17. Morisset, D.; Štebih, D.; Milavec, M.; Gruden, K.; Žel, J. Quantitative analysis of
343
food and feed samples with droplet digital PCR. PLOS One. 2013, 8(5), e62583.
344
18. Dobnik, D.; Spilsberg, B.; Košir, A. B.; Holst-Jensen, A.; Žel, J. Multiplex
345
quantification of 12 European Union authorized genetically modified maize lines
346
with droplet digital polymerase chain reaction. Anal. Chem. 2015, 87(16),
347
8218-8226.
ACS Paragon Plus Environment
Journal of Agricultural and Food Chemistry
348
19. Košir, A. B.; Spilsberg, B.; Holst-Jensen, A.; Žel, J.; Dobnik, D. Development and
349
inter-laboratory assessment of droplet digital PCR assays for multiplex
350
quantification of 15 genetically modified soybean lines. Sci. Rep. 2017, 7(1), 8601.
351
20. Meng, YN.; Liu, X.; Wang, S.; Zhang, DB.; Yang, LT. Applicability of plasmid
352
calibrant pTC1507 in quantification of TC1507 maize: an interlaboratory study. J.
353
Agric. Food Chem. 2012, 60(1), 23-28.
354
21. Jiang, LX.; Yang, LT.; Zhang, HB.; Guo, JC.; Mazzara, M.; Eede, G. V. D. International
355
collaborative study of the endogenous reference gene, sucrose phosphate
356
synthase (SPS), used for qualitative and quantitative analysis of genetically
357
modified rice. J. Agric. Food Chem. 2009, 57(9), 3525-3532.
358
22. Wu, YH.; Yang, LT.; Cao, YL.; Song, GW.; Shen, P.; Zhang, DB; Wu, G. Collaborative
359
validation of an event-specific quantitative real-time PCR method for genetically
360
modified rice event TT51-1 detection. J. Agric. Food Chem. 2013, 61(25),
361
5953-5960.
362
23. Wei, JJ.; Li, FW.; Guo, JC.; Li, X.; Xu, JF.; Wu, G.; Zhang, DB.; Yang, LT. Collaborative
363
ring trial of the papaya endogenous reference gene and its polymerase chain
364
reaction assays for genetically modified organism analysis. J. Agric. Food Chem.
365
2013, 61(47), 11363-11370.
366
24. Jacchia, S.; Nardini, E.; Bassani, N.; Savini, C.; Shim, J. H.; Trijatmiko, K.
367
International ring trial for the validation of an event-specific golden rice 2
368
quantitative real-time polymerase chain reaction method. J. Agric. Food Chem.
369
2015, 63 (20), 4954-4965.
370
25. Wei, JJ.; Le, HY.; Pan, AH.; Xu, JF.; Li, FW.; Li, X.; Quan, S.; Guo, JC.; Yang, LT.
371
Collaborative trial for the validation of event-specific PCR detection methods of
ACS Paragon Plus Environment
Page 18 of 29
Page 19 of 29
Journal of Agricultural and Food Chemistry
372 373 374 375
genetically modified papaya Huanong No.1. Food Chem. 2016, 194, 20-25. 26. Horwitz, W. Protocol for the design, conduct and interpretation of method performance studies. Pure Appl. Chem. 1995, 67, 331-343. 27. Definition of minimum performance requirements for analytical methods of GMO
376
testing-European
Network
of
GMO
Laboratories
(ENGL).
377
http://gmo-crl.jrc.ec.europa.eu/doc/Min_Perf_Requirements_Analytical_method
378
s.pdf.
379
28. Yang, Y.; Li, L.; Yang, H.; Li, XY.; Zhang, XJ.; Xu, JF.; Zhang, DB.; Jin, WJ.; Yang, LT.
380
Development of Certified Matrix-Based Reference Material as a Calibrator for
381
Genetically Modified Rice G6H1 Analysis. J. Agric. Food Chem. 2018. 66 (14),
382
3708-3715.
383
29. International Organization for Standardization. ISO 21571:2005, Foodstuffs-
384
Methods of analysis for the detection of genetically modified organisms and
385
derived products-Nucleic acid extraction. 2005.
386
30. Institution, B. S. Accuracy (trueness and precision) of measurement methods and
387
results -- part 2: basic method for the determination of repeatability and
388
reproducibility of a standard measurement method. 1994.
389 390
31. Trapman, S.; Burns, M.; Broll, H.; Macarthur, R.; Wood, R.; Zel, J. Guidance document on measurement uncertainty for GMO testing laboratories. 2007.
ACS Paragon Plus Environment
Journal of Agricultural and Food Chemistry
391
Figure captions
392
Figure 1. The 3’ end event-specific sequence of GM rice event G6H1. The italic capital
393
letters were the exogenous DNA sequence of G6H1, and the capital letters were rice
394
genome DNA sequence. The designed real-time PCR primers were labeled with
395
underline, and the probe was labeled with red color and square frame.
396 397
Figure 2. The specificity test and construction of standard curve of G6H1
398
event-specific real-time PCR method. a) The real-time PCR amplified curves from
399
different GM rice events, GM maize, GM soybean, GM cotton, and GM canola
400
samples, etc. b) Real-time PCR amplified curves and constructed standard curve of
401
G6H1 method.
402 403
Figure 3. Standard curves of G6H1 event-specific and rice SPS real-time PCR assays
404
employing serially diluted solutions of R1-R5 in each laboratory. For each laboratory,
405
three standard curves were obtained from three individual repeat for each real-time
406
PCR assay.
407 408
Figure 4. Relative deviation from the true value of three blind samples for all
409
participating laboratories.
410 411
Figure 5. Linear regression produced by plotting mean measurement concentration
412
(c) against reproducibility standard deviation (SDR).
ACS Paragon Plus Environment
Page 20 of 29
Page 21 of 29
Journal of Agricultural and Food Chemistry
413
Tables
414
Table 1. Primers and Probes used for G6H1 real-time PCR method Amplicon size PCR assay
Name
Sequence (5’-3’) (bp)
SPS
G6H1
SPS-F
TTGCGCCTGAACGGATAT
SPS-R
CGGTTGATCTTTTCGGGATG
SPS-P
FAM-TCCGAGCCGTCCGTGCGTC-BHQ
G6H1-F
AAGCGTCAATTTGTTTACACC
G6H1-R
TCGATCTGCTGCAGCTTG
G6H1-P
FAM-ATGGATGTATCGCCACCAGCACC-BHQ1
415
ACS Paragon Plus Environment
81
90
Journal of Agricultural and Food Chemistry
416
Page 22 of 29
Table 2. The LOD and LOQ evaluation of G6H1 real-time PCR method using diluted calibrators Mean Ct values
DNA amount
Amount (copies) Bias (%)
SD
RSD
102.12
2.12
6.43
6.29
42.78
46.84
6.31
4.26
9.09
9.41
7.68
8.65
13.47
0.89
10.34
1.68
4.68
3.07
38.54
1.53
49.81
(copies)
Rep 1
Rep 2
Rep 3
Rep 1
Rep 2
Rep 3
Mean
100.00
33.59
33.48
33.41
95.46
103.08
108.24
50.00
34.48
34.61
34.74
51.29
46.84
10.00
36.98
36.91
37.20
8.96
5.00
38.25
39.38
37.91
3.69
ACS Paragon Plus Environment
Page 23 of 29
Journal of Agricultural and Food Chemistry
417
Table 3. Slope, PCR efficiency, and R2 values of the standard curves Lab
A
B
C
D
E
F
G
H
Rep
G6H1
SPS
slope
PCR efficiency (%)
Linearity (R )
slope
PCR efficiency (%)
Linearity (R2)
1
-3.41
96.45
0.9998
-3.42
96.06
0.9999
2
-3.40
96.84
0.9995
-3.36
98.44
0.9989
3
-3.46
94.54
0.9993
-3.40
96.84
0.9999
1
-3.45
94.92
1.0000
-3.36
98.44
1.0000
2
-3.44
95.30
0.9999
-3.29
101.35
0.9997
3
-3.45
94.92
1.0000
-3.56
90.94
0.9958
1
-3.49
93.43
0.9998
-3.31
100.50
0.9999
2
-3.39
97.24
0.9999
-3.35
98.84
0.9998
3
-3.41
96.45
0.9994
-3.24
103.54
0.9992
1
-3.41
96.45
0.9999
-3.42
96.06
0.9995
2
-3.41
96.45
0.9996
-3.40
96.84
0.9999
3
-3.38
97.63
0.9999
-3.45
94.92
0.9995
1
-3.48
93.80
0.9988
-3.38
97.63
0.9998
2
-3.47
94.17
0.9997
-3.38
97.63
0.9997
3
-3.35
98.84
0.9999
-3.38
97.63
0.9998
1
-3.58
90.25
0.9993
-3.55
91.29
0.9997
2
-3.51
92.71
0.9987
-3.54
91.64
0.9994
3
-3.48
93.80
0.9996
-3.45
94.92
0.9980
1
-3.33
99.66
0.9826
-3.40
96.84
0.9995
2
-3.48
93.80
1.000
-3.25
103.09
1.0000
3
-3.36
98.44
0.9994
-3.40
96.84
0.9920
1
-3.42
96.06
0.9998
-3.42
96.06
0.9998
2
-3.40
96.84
0.9994
-3.37
98.03
0.9988
3
-3.47
94.17
0.9993
-3.40
96.84
0.9999
-3.43
95.55
0.9989
-3.40
97.13
0.9991
Mean
2
418
ACS Paragon Plus Environment
Journal of Agricultural and Food Chemistry
419
Page 24 of 29
Table 4. Summary of validation results for the G6H1 event-Specific method Expected value Blind samples S1 (4.983%)
S2 (0.997%)
S3 (0.499%)
Labs returning results
8
8
8
Samples per lab
3
3
3
Number of outliers
0
0
0
Reason for exclusion
/
/
/
Mean value
5.02%
1.09%
0.53%
Repeatability SD
0.20%
0.04%
0.03%
Repeatability RSD
3.91%
3.45%
6.33%
Reproducibility SD
0.29%
0.06%
0.04%
Reproducibility RSD
5.68%
5.79%
7.89%
Bias (absolute value)
0.04%
0.09%
0.03%
Bias
0.74%
9.33%
6.21%
420
ACS Paragon Plus Environment
Page 25 of 29
421
Journal of Agricultural and Food Chemistry
Table 5. Determined GM% values of three blind samples from eight laboratories S1 (4.983%)
S2 (0.997%)
S3 (0.499%)
Lab Sub sample 1
Sub sample 2
Sub sample 3
Sub sample 1
Sub sample 2
Sub sample 3
Sub sample 1
Sub sample 2
Sub sample 3
A
4.72%
4.75%
4.56%
1.14%
1.23%
1.24%
0.57%
0.56%
0.49%
B
5.27%
4.90%
4.59%
1.02%
1.03%
1.05%
0.47%
0.44%
0.49%
C
5.36%
5.43%
5.81%
1.15%
1.13%
1.05%
0.61%
0.56%
0.52%
D
5.11%
4.57%
4.65%
1.11%
1.18%
1.15%
0.51%
0.58%
0.48%
E
5.04%
4.75%
4.61%
1.13%
1.08%
1.07%
0.55%
0.61%
0.62%
F
5.02%
5.26%
4.97%
1.01%
0.98%
1.04%
0.60%
0.57%
0.56%
G
5.25%
5.26%
5.13%
1.14%
1.06%
1.03%
0.47%
0.52%
0.51%
H
5.20%
5.05%
5.33%
1.07%
1.02%
1.06%
0.53%
0.52%
0.49%
ACS Paragon Plus Environment
Journal of Agricultural and Food Chemistry
422
Figures
423
Figure 1
424 425 426 427
Figure 2.
428 429
ACS Paragon Plus Environment
Page 26 of 29
Page 27 of 29
Journal of Agricultural and Food Chemistry
430
Figure 3.
431 432
ACS Paragon Plus Environment
Journal of Agricultural and Food Chemistry
433
Figure 4
434 435 436
Figure 5
437
ACS Paragon Plus Environment
Page 28 of 29
Page 29 of 29
Journal of Agricultural and Food Chemistry
438
TOC
439
ACS Paragon Plus Environment