Subscriber access provided by UNIVERSITY OF ADELAIDE LIBRARIES
Article
One novel multiple-target plasmid reference molecule targeting eight genetically modified canola events for GM canola detection Zhuqing Li, Xiang Li, Canhua Wang, Guiwen Song, Liqun Pi, Lan Zheng, Dabing Zhang, and Litao Yang J. Agric. Food Chem., Just Accepted Manuscript • DOI: 10.1021/acs.jafc.7b02453 • Publication Date (Web): 28 Aug 2017 Downloaded from http://pubs.acs.org on August 30, 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 44
Journal of Agricultural and Food Chemistry
1
One novel multiple-target plasmid reference molecule
2
targeting eight genetically modified canola events for GM
3
canola detection
4 5
Zhuqing Li1,†, Xiang Li2,†, Canhua Wang1, Guiwen Song3, Liqun Pi1, Lan Zheng1, Dabing
6
Zhang1,4, Litao Yang1*.
7 8
1
9
Organisms, School of Life Sciences and Biotechnology, Shanghai Jiao Tong University,
National Center for the Molecular Characterization of Genetically Modified
10
Shanghai 200240, China.
11
2
12
1208 Minsheng Road, Shanghai 200135, P. R. China
13
3
14
Republic of China. Beijing 100025, China.
15
4
16
Collaborative Innovation Center of Regional Modern Agriculture and Environmental
17
Protection, Huaian 223300, China
Shanghai Entry-Exit Inspection and Quarantine Bureau, GMO Detection Laboratory,
Development Center of Science and Technology, Ministry of Agriculture of People’s
Key Laboratory of Crop Marker-Assisted Breeding of Huaian Municipality, Jiangsu
18 19
†These two authors contributed equally to this work.
20
*
21
34207174; Email:
[email protected].
To whom correspondence should be addressed: Tel.: +86 21 34205073; Fax: +86 21
22 1
ACS Paragon Plus Environment
Journal of Agricultural and Food Chemistry
Page 2 of 44
23
Abstract
24
Multiple-target plasmid DNA reference materials have been generated and utilized as
25
good substitutes of matrix-based reference materials in the analysis of genetically
26
modified organisms (GMOs). Herein, we report the construction of one
27
multiple-target plasmid reference molecule, pCAN, which harbors eight GM canola
28
event-specific sequences (RF1, RF2, MS1, MS8, Topas 19/2, Oxy235, RT73, and T45)
29
and a partial sequence of the canola endogenous reference gene PEP. The
30
applicability of this plasmid reference material in qualitative and quantitative PCR
31
assays of the eight GM canola events was evaluated, including the analysis of
32
specificity, limit of detection (LOD), limit of quantification (LOQ), and performance of
33
pCAN in the analysis of various canola samples, etc. The LODs are 15 copies for RF2,
34
MS1, and RT73 assays using pCAN as the calibrator, and 10 genome copies for the
35
other events. The LOQ in each event-specific real-time PCR assay is 20 copies. In
36
quantitative real-time PCR analysis, the PCR efficiencies of all event-specific and PEP
37
assays are between 91% and 97 %, and the squared regression coefficients (R2) are all
38
higher than 0.99. The quantification bias values varied from 0.47% to 20.68% with
39
relative standard deviation (RSD) from 1.06% to 24.61% in the quantification of
40
simulated samples. Furthermore, 10 practical canola samples sampled from
41
imported shipments in the port of Shanghai, China were analyzed employing pCAN
42
as the calibrator, and the results were comparable with those assays using
43
commercial certified materials as calibrator. Concluding from these results, we
44
believe that this newly developed pCAN plasmid is one good candidate for being a 2
ACS Paragon Plus Environment
Page 3 of 44
Journal of Agricultural and Food Chemistry
45
plasmid DNA reference material in the detection and quantification of the eight GM
46
canola events in routine analysis.
47
Keywords: Multiple-target reference plasmid; Genetically Modified canola;
48
Qualitative PCR; Quantitative real-time PCR.
49
3
ACS Paragon Plus Environment
Journal of Agricultural and Food Chemistry
Page 4 of 44
50
Introduction
51
Canola (Brassica napus L.), one of the staple oil-crops, is a significant source of edible
52
oil and plant proteins in the world. It is mainly planted in China, India, and Canada
53
and holds an important status in the international agricultural product trade. GM
54
canola is the fourth commercialized GM crops. In 2016, the planted area of GM
55
canola has reached about 8.6 million hectares in global area, such as USA, Canada,
56
and Argentina, etc. 1 To date, a total of 43 GM Argentine canola and Polish canola
57
events, including the gene stacked events, have been approved for commercialization
58
and planting. In China, seven GM canola lines (MS1×RF1, MS8×RF3, MS1×RF2, RT73,
59
Oxy235, T45, and Topas 19/2) including nine GM events have been approved for
60
usage as raw and processed materials since 2004.2 China now imports more than one
61
million tons of canola each year, and 60 percent are GM canola. According to the
62
GMO labeling regulation issued by Chinese government in 2001, products derived
63
from GM canola should be labeled with a zero threshold. 3 To implement the labeling
64
regulation, pre-validated event-specific detection methods and reference materials
65
for each GM line need to be well developed. Currently, several techniques have been
66
developed for GMO analysis in the past decade, such as PCR 4-5, real-time PCR 6,
67
digital PCR 7-8, microarray 9, loop-mediated isothermal amplification (LAMP) 10-11, and
68
next generation sequencing 12. Regardless of which detection method to choose,
69
reference materials (RMs) which are used as calibrator are essential for quality
70
insurance to deliver accurate GM contents.
71
Reference materials play key roles in the development and evaluation of new 4
ACS Paragon Plus Environment
Page 5 of 44
Journal of Agricultural and Food Chemistry
72
methods, calibration of instruments, and routine analysis in many fields, such as
73
clinical diagnosis, food safety testing, chemical analysis, and GMO detection, etc. 13 In
74
GMO analysis, matrix-based RMs, genomic DNA RMs, and plasmid DNA RMs are the
75
mainly used RMs. Among them, matrix-based RMs and genomic DNA RMs have the
76
advantage of being most closely to practical samples. However, the production of
77
matrix-based or genomic DNA RM was limited by the availability of GM materials
78
with high purity and clear genotype. 2 Currently, commercial matrix-based and
79
genomic DNA certified reference materials (CRMs) are mainly produced by Reference
80
Materials and Measurements, European Commission (IRMM-EU) and the American
81
Oil Chemists’ Society (AOCS).
82
In 2002, Kuribara et al., developed plasmid reference molecules for GM soybean and
83
maize detection, and validated their applicability in GM soybean and maize
84
quantification through a collaborative ring trial. 14 Since then, more and more
85
plasmid reference molecules have begun to be used as calibrator in GMO detection.
86
The advantages and disadvantages of plasmid reference molecules have been well
87
discussed in several previous publications 14-17. Their main advantages are: i) plasmid
88
reference molecules can be easily constructed and produced with high quantity and
89
quality; ii) multiple GM sequences from different GM events can be recombined into
90
one plasmid, and one plasmid can be used as calibrator for multiple GM events; iii)
91
the production cost is quite low. The disadvantages include easy contamination,
92
limitation of specific DNA fragments and corresponding methods, and different
93
amplification efficiency in PCR reactions compared with using genomic DNA. Despite 5
ACS Paragon Plus Environment
Journal of Agricultural and Food Chemistry
Page 6 of 44
94
of these disadvantages, plasmid reference molecules have been widely accepted in
95
GMO analysis. Up to date, several publications have reported the development and
96
application of new plasmid reference molecules in GMO analysis. For example, Yang
97
et al. developed a multiple-target plasmid reference molecule for quantification of
98
nine GM maize events. 18 Zhang et al. constructed a plasmid reference material for
99
multiplex quantitative analysis of GM soybean GTS 40-3-2. 19 Lievens et al. validated
100
pJANUS-02-001 as a valid plasmid calibrator for GM soybean GTS 40-3-2. 20 Meng et
101
al. developed and validated a plasmid reference material for GM maize TC1507. 21
102
Kim et al. developed a simple plasmid reference for GM wheat MON71800. 22 Pi et al.
103
reported a plasmid reference material targeting five GM soybean events. 23 Also,
104
IRMM-EU developed and certified four plasmid DNA CRMs for GM maize events
105
MON810 (ERM-AD413), NK603 (ERM-AD415), 356043 (ERM- AD425), and 98140
106
(ERM-AD427).
107
For GM canola detection, the qualitative PCR 24-25, real-time PCR 26, and Microarray 27
108
methods have been developed and validated employing different GM canola RMs.
109
Up to date, a total of nine CRMs have been developed by AOCS and are commercially
110
available. All of the nine CRMs are matrix-based CRM or genomic DNA CRM,
111
including one seed CRM for event GT73/RT73, one seed powder CRM of MON88302,
112
and seven leaf genomic DNA CRMs for GM events MS1, Rf1, Rf2, Rf3, Topas 19/2, T45,
113
and MS8. However, few plasmid CRM or multiple-target CRM for GM canola events
114
detection were developed except for GT73 and Oxy235 events. 24, 28 In this study, we
115
developed one multiple-target plasmid CRM, pCAN, which contains the 6
ACS Paragon Plus Environment
Page 7 of 44
Journal of Agricultural and Food Chemistry
116
event-specific DNA fragments of eight GM canola events and a partial sequence of
117
canola endogenous reference gene PEP. The applicability of pCAN as a calibrator in
118
qualitative PCR and quantitative real-time PCR analysis of GM canola events were
119
also evaluated.
7
ACS Paragon Plus Environment
Journal of Agricultural and Food Chemistry
Page 8 of 44
120
Material and Methods
121
Materials and DNA Extraction
122
Reference materials of GM maize NK603 (ERM-BF415f) were purchased from
123
Institute for Reference Materials and Measurements (IRMM). Leaf tissue DNA or
124
seed powder of GM canola event RT73 (0304-B2), Ms1 (0711-A2), Rf1 (0711-B2), Rf2
125
(0711-C2), Topas 19/2 (0711-D3), T45 (0208-A5), Ms8 (0306-F6) and Rf3 (0306-G5)
126
were purchased from American Oil Chemists’ Society (AOCS). Seeds or leaves of GM
127
canola event Oxy235 and GM rice TT51-1 were kindly provided by the developers.
128
Non-GM canola seeds were purchased from a local market in Shanghai, China and
129
validated to be free of GM contents in our laboratory. Ten practical canola seed
130
samples were sampled from the port of Shanghai, China and supplied by Shanghai
131
Entry-exit Inspection and Quarantine Bureau. A Plant DNA Mini-Prep Kit (Ruifeng
132
Agro-tech Co. Ltd., Shanghai, China) was used to extract and purify genomic DNA
133
after the plant materials were ground in liquid nitrogen. A Plasmid Mini Kit (Watson
134
Biotechnologies, Inc., Shanghai, China) was used for isolating and purifying plasmid
135
DNA. The concentration and purity of each extracted genomic DNA or plasmid DNA
136
was evaluated according to the absorbance at 230nm, 260nm and 280 nm
137
determined using a NanoDrop 1000 spectrophotometer (NanoDrop Technologies,
138
LLC, Wilmington, DE, USA). The integrity of extracted genomic DNAs and plasmid
139
DNAs were further checked by 0.8% agarose gel electrophoresis. The copy number in
140
each genomic DNA or plasmid DNA sample was calculated according to the
141
determined DNA quantity and genome size 29. 8
ACS Paragon Plus Environment
Page 9 of 44
Journal of Agricultural and Food Chemistry
142
Oligonucleotide Primers and Probes
143
Primers for cloning the event-specific DNA fragments were designed based on the
144
event-specific sequence of each event using Primer Premier 5.0 software (PREMIER
145
Biosoft International, CA). Primers for qualitative PCR analysis were adopted from
146
Chinese national standards for GM canola detection. Primers and probes for
147
quantitative real-time PCR assays of the eight events were designed using Beacon
148
Designer 7.2 software (PREMIER Biosoft International, Palo Alto,CA, USA). For GM
149
canola events RF1, RF2, MS1, MS8, Oxy235 and T45, the 5' end flanking sequence
150
were used as target fragments. For GM canola events Topas 19/2 and GT73/RT73,
151
the 3'-end event-specific sequence was amplified in real-time PCR analysis. The
152
primer pair of Q-PEP-1F/2R and probe Q-PEP-P were designed and used to amplify
153
canola endogenous reference gene PEP. The primers and probes for real-time PCR
154
assays were listed in Table 1. All primers and probes were synthesized by Invitrogen
155
Co., Ltd. (shanghai China).
156
Construction of the multiple-target plasmid pCAN
157
To construct the multiple-target pCAN by molecular cloning, the pEASY-T3 Cloning
158
Vector (Beijing TransGen Biotech Co., Ltd. China) was selected as the basic vector,
159
and nine DNA fragments (eight event-specific sequences of RF1, RF2, MS1, MS8,
160
Topas19/2, Oxy235, RT73, and T45 as well as the partial sequence of canola PEP gene)
161
were inserted into the pEASY-T3 vector by conventional molecular cloning
162
techniques, such as overlapping PCR, restriction enzyme digestion and ligation. The
163
sequences of eight canola events and PEP gene were from NCBI database, including 9
ACS Paragon Plus Environment
Journal of Agricultural and Food Chemistry
Page 10 of 44
164
RF1 with Accession number EU090199.1, RF2 with EU090197.1, MS1 with
165
EU090198.1, MS8 with EU020107.1, Topas 19/2 with EU124676.1, Oxy235 with
166
EU099578.1, RT73 with HJ231561.1, T45 with FJ154954.1 and PEP gene with
167
D13987.1.
168
The eight event-specific DNA fragments and PEP gene fragment were first amplified
169
individually from templates by conventional PCR, and the corresponding primers are
170
listed in Supplementary Table S1. Next, three integrated DNA fragments of
171
MS1-RT73 named F1 (MS1+RT73), F2 (Oxy235+ MS8), and F3 (T45+Topas+RF2+RF1)
172
were obtained by overlapping PCR. Firstly, fragment F1 was cloned into the pEASY-T3
173
vector by T-A ligation, then fragment F2 was integrated by endogenous restriction
174
enzyme Pst I and Bgl II, and finally fragment F3 was ligated after enzyme Bgl II and
175
Xho I digestion. The sequence of the new integrated vector was confirmed by
176
sequencing. Then the new big pEASY-T3 vector was digested by enzyme Bgl II and
177
Xba I, and then ligated with the PEP gene fragment after digestion by enzyme BamH1
178
and Sal I. The finally assembled pCAN was confirmed twice by DNA sequencing
179
(Shanghai Meiji Biological Science and Technology Co. Ltd, Shanghai, China).
180
Qualitative PCR assays
181
Qualitative PCR assays were performed to evaluate the specificity and sensitivity
182
when pCAN was used as templates in PCR reactions. All PCR reactions were carried
183
out using a Verity Thermal Cycler (Applied Biosystems, Carlsbad, CA, USA) in a 25 μl
184
reaction volume, including 1×PCR buffer, 100 nM dNTP, 200 nM of each primer, 1.5
185
units of Taq DNA polymerase (TaKaRa Biotechnology Co., Ltd.,Dalian, China), 2 μl of 10
ACS Paragon Plus Environment
Page 11 of 44
Journal of Agricultural and Food Chemistry
186
DNA template, and ddH2O. The PCR amplification program was: 95°C for 7 min; 35
187
cycles of 95 °C for 30 s, 58 °C for 30 s, 72 °C for 30 s; a final extension step of 72 °C for
188
7 min. PCR products were examined by 2% agarose gel electrophoresis with Gelred
189
staining. Each PCR test was repeated three times and each time with three parallel
190
replicates.
191
Quantitative real-time PCR assays
192
The real-time PCR reactions were carried out using a Prism ABI 7900 detection
193
system (Applied Biosystems, Carlsbad, CA, USA) in a reaction volume of 25 μL, which
194
contained 1 × TaqMan Universal PCR Master Mix (Applied Biosystems, Carlsbad, CA,
195
USA), 400 nM of each primer, 200nM of the probe, 5 μL DNA template, and ddH2O to
196
fill up the final volume. The Real-time PCR program was: pre-denaturation at 95 °C
197
for 10 min, followed by 45 cycles of 95 °C for 15 s and 60 °C for 60 s. The fluorescent
198
signals were monitored and analyzed in the extension step using the Sequence
199
Detector Software (SDS) V2.4. Each real-time PCR reaction was repeated three times,
200
and each time with three parallel reactions.
201
Quantification of GM canola content using pCAN as a calibrator
202
The performance of pCAN as a calibrator in real-time PCR analysis was evaluated
203
using three mixed DNA samples and ten practical canola samples. Since the
204
discrepancy of PCR efficiency between using plasmid and genomic DNA templates
205
might influence the ratio of exogenous to endogenous genome copy number, a
206
conversion factor (Cf) is often used. 14, 21-23, 28 In general, a Cf value is defined as the
207
genome copy number ratio of the event-specific DNA to endogenous reference gene 11
ACS Paragon Plus Environment
Journal of Agricultural and Food Chemistry
Page 12 of 44
208
in a homogenous genomic DNA sample using the plasmid DNA CRM as a calibrator,
209
and calculated using the formula Cf =
210
determine the Cf values of pCAN, homogenous genomic DNAs of each GM event
211
were prepared at three concentrations (10.0, 5.0, 1.0 ng/μl). The Cf value for each
212
event is defined as the average of three Cf values determined in the three DNAs at
213
different concentrations. For Cf determination, each real-time PCR reaction was
214
repeated three times, and each time with three parallel replicates. After the Cf values
215
were determined, the GM content in each canola sample was calculated using the
216
formula GM amount (%) = × +.
. To
' × )**
217
12
ACS Paragon Plus Environment
Page 13 of 44
Journal of Agricultural and Food Chemistry
218
Results and Discussion
219
Construction of the reference Plasmid pCAN
220
We have developed a novel reference plasmid pCAN containing eight event-specific
221
GM canola DNA sequences and one canola endogenous reference gene PEP (Figure
222
1). The length of pCAN is 7972 bp including 3253 bp of the eight event-specific GM
223
canola and endogenous PEP sequences. The full length pCAN DNA sequence was
224
determined by DNA sequencing and the sequences of the eight event-specific GM
225
canola and endogenous PEP sequences were confirmed. The confirmed sequences of
226
the inserted GM event and PEP DNA fragments are same with our design and
227
expectation (As shown in Supplementary File 1).
228
Qualitative PCR assays using pCAN as calibrator
229
Specificity. For a given GM event, an event-specific DNA sequence is a unique
230
sequence specific to that GM event, so the specificity of an event-specific PCR assay
231
should be quite high. For a multiple-target reference plasmid, however, since
232
multiple event-specific DNA sequences are within close proximity in the same vector,
233
cross interferences might happen and unexpected amplicons might be produced in
234
an event-specific PCR assay using the reference plasmid as template. Therefore, the
235
specificity of using pCAN in each event-specific PCR assay should be tested
236
individually. We performed qualitative PCR assays of eight GM canola events (RF1,
237
RF2, MS1, MS8, Topas 19/2, Oxy235, RT73, and T45) as well as endogenous gene PEP,
238
and the specificity of each PCR assay when using pCAN as the template was
239
evaluated. In all nine PCR assays, only expected amplicons corresponding to their 13
ACS Paragon Plus Environment
Journal of Agricultural and Food Chemistry
Page 14 of 44
240
specific targets were obtained, and no unexpected amplicon was observed (Figure 2).
241
For instance, a 200 bp amplicon was expected and observed in RF1 event PCR assay
242
when pCAN plasmid and RF1 genomic DNA were used as templates, and no other
243
amplified DNA fragment was observed. The same held true for PCR assays of the
244
other GM canola events (RF2, MS1, MS8, Topas 19/2, Oxy235, RT73, and T45), as well
245
as assays of non-GM canola, GM maize NK603, GM rice TT51-1, and GM soybean
246
MON89788 (Figure 2A and 2J). In PEP assay, the expected 248 bp PEP DNA fragment
247
was specifically amplified using pCAN and each GM canola genomic DNAs as
248
template; no PEP amplification was observed when GM maize, soybean, and rice
249
genomic DNAs were used as template (Figure 2I and 2J). Also, the high specificity
250
was confirmed using the real-time PCR assays with Sybr Green I dyer (As shown in
251
supplementary Figure S1). These results indicate that pCAN has high specificity for
252
each of the nine targets contained in the plasmid, and it could be used as calibrator
253
in PCR assays of the eight GM canola events.
254
Sensitivity. The sensitivities of using pCAN as calibrator in the nine PCR assays were
255
evaluated using serially diluted pCAN DNA solutions. These serially diluted pCAN
256
DNA solutions had concentrations of 5000, 1000, 250, 125, 100, 50, 25, 15, and 10
257
copies per reaction, respectively. The test results showed that the sensitivities of
258
using pCAN in the nine PCR assays could be as low as 15 genome copies per reaction
259
except for the assays of events RF2, MS1, and GT73, which have sensitivities of 25
260
copies per reaction (Figure 3). The sensitivity of 25 genome copies is approximately
261
equal to 0.05% GM DNA in 100 ng canola genomic DNA, which is accepted in the 14
ACS Paragon Plus Environment
Page 15 of 44
Journal of Agricultural and Food Chemistry
262
requirements of GMO labeling regulations around the world. 2, 30
263
Applicability of pCAN as calibrator in quantitative real-time PCR analysis
264
In order to evaluate the applicability of pCAN as calibrator in quantitative real-time
265
PCR assays of the eight GM canola events, primary real-time PCR parameters were
266
tested according to the guideline of real-time PCR method for GMOs, such as PCR
267
efficiency, linearity of standard curve, repeatability and reproducibility, limit of
268
detection (LOD), limit of quantification (LOQ), and conversion factor (Cf), etc. 30
269
PCR efficiencies and linearity of standard curves. Standard curves of real-time PCR
270
assays of the 10-fold serially diluted pCAN DNA solutions (2×105, 2×104, 2×103, 2×102,
271
and 20 copies/μL) were generated, and PCR efficiencies and linearity of the standard
272
curves were calculated. The PCR efficiencies of the eight event-specific GM canola
273
assays and endogenous PEP gene assay ranged from 91% to 97%, which were within
274
the range of 90–110 % recommended by the European Commission (Figure 4). 30 The
275
linearity of standard curves of the nine real-time PCR assays varied from 0.9987 to
276
1.0000 (Figure 4), indicating excellent linearity between the DNA copy numbers in
277
the template and corresponding fluorescent intensities (Ct). The high PCR efficiency
278
and good linearity of standard curves indicate that pCAN could be used as calibrator
279
in real-time PCR assays of the eight GM canola events and the endogenous PEP gene.
280
Repeatability and reproducibility. Repeatability and reproducibility of the nine
281
real-time PCR assays using pCAN as calibrator were evaluated by calculating the
282
relative standard deviation (RSD) of the nine repeats of one PCR assay (three times
283
and each time with three parallel replicates). To determine repeatability, the nine 15
ACS Paragon Plus Environment
Journal of Agricultural and Food Chemistry
Page 16 of 44
284
PCR reactions were performed by the same operator. For reproducibility, the nine
285
PCR reactions were performed by three different operators (each operator
286
performed PCR assay once with three parallel replicate reactions). As shown in Table
287
2, the standard deviation of repeatability (SDr) ranged from 0.01 to 0.41, and the
288
relative standard deviation of repeatability (RSDr) ranged from 0.05% to 1.60%. The
289
standard deviation of reproducibility (SDR) ranged from 0.04 to 0.39, and the relative
290
standard deviation of reproducibility (RSDR) ranged from 0.13% to 1.79%. The
291
acceptable RSDr and RSDR by ENGL guideline were ≤25%, 30 and our RSDr and RSDR
292
values are well below 25%, indicating pCAN has good repeatability and
293
reproducibility when employed as calibrator for the nine real-time PCR assays.
294
LOD and LOQ.
295
In real-time PCR analysis, the LOD and LOQ are the lowest genome copy number of
296
the target DNA that can be reliably detected and accurately quantified with a
297
confidence interval of more than 95%. To determine the LODs and LOQs of pCAN in
298
each assay, pCAN DNA solutions with concentrations of 4, 2 and 1 copy/μl) were
299
prepared as PCR templates, and 20 real-time PCR reactions for each concentration
300
was performed in parallel. For the nine real-time PCR assays (eight event-specific
301
assays and PEP assay), good PCR amplification curves and reliable Ct values were
302
observed in all 20 parallel reactions when 20 or 10 copies of pCAN DNA was used as
303
template, indicating that the LODs are as low as 10 genome copies for all 9 assays.
304
When 20 genome copies of pCAN DNA was used as template, the RSD values of the
305
20 parallel reactions ranged from 9.38% to 22.53%, and the quantification biases 16
ACS Paragon Plus Environment
Page 17 of 44
Journal of Agricultural and Food Chemistry
306
were between 1.34% and 14.36% (Table 3). However, when 10 genome copies of
307
pCAN DNA was used as template, the RSD values were higher than 25% (27.60% to
308
49.12%), even though the quantification biases for RF2, MS8, and Topas 19/2 assays
309
were lower than 25% (Table 3). Concluding from the abovementioned results, the
310
LOQs of the nine real-time PCR assays using pCAN plasmid DNA as templates are
311
determined to be as low as 20 genome copies.
312
Conversion Factor. The Cf value for pCAN in each event-specific GM canola assay was
313
determined using homogenous genomic DNAs of each event at three different
314
concentrations (10.0, 5.0, and 1.0 ng/μl). The results were listed in Table 4. Cf values
315
for RF1, RF2, MS1, MS8, Topas 19/2, Oxy235, RT73 and T45 event-specific
316
quantitative PCR assays were 4.63, 2.93, 0.87, 1.29, 1.31, 1.93, 2.04, and 2.04,
317
respectively; the RSD of these Cfs are ranged from 3.00% to 16.35%, far below 35%.
318
In theory, Cf values should be consistent with the value of 1.0 among different
319
event-specific assays if the endogenous reference gene is single copy in haploid
320
genome and there is no difference in PCR efficiency between using genomic DNA and
321
plasmid DNA as template. 14 Because of the complex genome of canola, the haploid
322
copy number of currently used endogenous reference genes, acetyl-CoA carboxylase
323
gene (BnACCg8), phosphoenolpyruvate carboxylase (PEP), oleoyl hydrolase gene
324
(FatA), high-mobility-group protein I/Y gene (HMG-I/Y) and cruciferin A gene (CruA),
325
were two copies at least, and the ideal Cf value should be 2.0. Our results showed
326
that the Cf values of MS8, Topas 19/2, Oxy235, RT73 and T45 were around 2.0, RF1,
327
RF2, and MS1 have slightly large variation, with the lowest Cf for MS1 at 0.87 and the 17
ACS Paragon Plus Environment
Journal of Agricultural and Food Chemistry
Page 18 of 44
328
highest cf value for RF1 at 4.63.
329
To assess the discrepancy between experimental and theoretical Cf values, we
330
constructed PCR standard curves (As shown in Figure 4 and Table S2) using genomic
331
DNA templates of each of the eight GM canola events at different concentrations
332
(50.0, 25.0, 5.0, 2.5, 0.5 ng/reaction). We found that the deviations of slopes and
333
intercepts of plasmid and genomic DNA standard curves (Table S2) might cause the
334
difference between the experimental and theoretical Cf values, especially for the
335
deviation of intercepts.31 For instance, in RF1 event and endogenous PEP gene assays,
336
standard curves between reactions using genomic DNA and pCAN plasmid DNA
337
templates had similar slopes (-3.4471 and -3.4466) but different intercepts of 39.583
338
and 37.052, the ⊿intercept of 2.531 was the mainly origin caused the high Cf value.
339
In RF1 event and endogenous PEP gene assays, the standard curves between
340
reactions using genomic DNA and pCAN plasmid DNA templates were almost
341
coincidnece resulted ideal Cf value. Also, one previous report on the comparison of
342
five canola endogenous reference genes indicated that the significant allelic variation
343
might cause the Ct values variation among different cultivars, such as HMG-I/Y and
344
BnACCg8. 32 The allelic variation of PEP gene among the eight GM canola events
345
might result differences of Cf values. Therefore, Cf values should be used to
346
normalize the ratio of exogenous to endogenous genome copies in quantitative
347
analysis of the eight GM canola events when using reference plasmid pCAN as
348
calibrator.
349
Performance of pCAN as a calibrator in quantitative analysis of canola samples 18
ACS Paragon Plus Environment
Page 19 of 44
Journal of Agricultural and Food Chemistry
350
To test the performance of pCAN as a reference calibrator in quantitative real-time
351
PCR analysis of GM canola, two sets of canola samples were tested. One set was
352
three simulated DNA samples (G1, G2 and G3) prepared by mixing GM canola
353
genomic DNA with non-GM canola genomic DNA at different ratios. The three
354
simulated samples are named G1, G2 and G3, which has 5.0%, 3.0%, and 1.0% GM
355
canola genomic DNA, respectively. The other set was ten canola seed samples
356
(named from R1 to R10) sampled from the shipped canola seeds to Shanghai port
357
and was kindly provided by Shanghai Entry-exit Inspection and Quarantine Bureau.
358
The accuracy and precision of the test results were used to evaluate the performance
359
of pCAN in real-time PCR assays of the canola samples. The accuracy of a test result
360
was measured as degree of bias, which was calculated using the formula of
361
|) -
|
-
× 100%. The precision of a test result was evaluated using the
362
standard deviation (SD) and relative standard deviation (RSD).
363
As shown in Table 5 and Figure 5, the mean of the tested values of the GM content in
364
each mixed canola genomic DNA sample using pCAN as the reference calibrator was
365
close to the actual value after applying corresponding Cf. The quantification bias of
366
all samples are ranged from 0.47% to 20.68%, which were less than the acceptance
367
criterion of