Use of next generation sequencing data to develop a qPCR method for specific detection of EU-unauthorized genetically modified Bacillus subtilis overproducing riboflavin
BMC Biotechnology volume 15, Article number: 103 (2015)
Recently, the presence of an unauthorized genetically modified (GM) Bacillus subtilis bacterium overproducing vitamin B2 in a feed additive was notified by the Rapid Alert System for Food and Feed (RASFF). This has demonstrated that a contamination by a GM micro-organism (GMM) may occur in feed additives and has confronted for the first time,the enforcement laboratories with this type of RASFF. As no sequence information of this GMM nor any specific detection or identification method was available, Next GenerationSequencing (NGS) was used to generate sequence information. However, NGS data analysis often requires appropriate tools, involving bioinformatics expertise which is not alwayspresent in the average enforcement laboratory. This hampers the use of this technology to rapidly obtain critical sequence information in order to be able to develop a specific qPCRdetection method.
Data generated by NGS were exploited using a simple BLAST approach. A TaqMan® qPCR method was developed and tested on isolated bacterial strains and on the feed additive directly.
In this study, a very simple strategy based on the common BLAST tools that can be used by any enforcement lab without profound bioinformatics expertise, was successfully used toanalyse the B. subtilis data generated by NGS. The results were used to design and assess a new TaqMan® qPCR method, specifically detecting this GM vitamin B2 overproducing bacterium. The method complies with EU critical performance parameters for specificity, sensitivity, PCR efficiency and repeatability. The VitB2-UGM method also could detect the B. subtilis strain in genomic DNA extracted from the feed additive, without prior culturing step.
The proposed method, provides a crucial tool for specifically and rapidly identifying this unauthorized GM bacterium in food and feed additives by enforcement laboratories. Moreover, this work can be seen as a case study to substantiate how the use of NGS data can offer an added value to easily gain access to sequence information needed to develop qPCR methods to detect unknown andunauthorized GMO in food and feed.
Riboflavin (vitamin B2) is necessary for basic cell metabolism but is not synthetized by higher animals and needs to be ingested. Therefore, the vitamin is commercialized for nutritional use in the fortification of various food and feed products. It is also commonly used as a colorant in food such as ice cream, cheeses, mayonnaises and sauces, and as a medical identification aid, because of its intense yellow color [1–3].
At the present time, fermentation using microbiological processes is widely applied to commercially produce riboflavin [4, 5] and genetically modified (GM) Bacillus (B.) subtilis overproducing riboflavin is one of the main organism currently used for riboflavin’s industrial production [2, 4, 6–8]. One common genetic method used to enhance riboflavin production derived from fermentation by GM B. subtilis, is to clone the riboflavin biosynthesis operon (ribG, ribB, ribA and ribH) from B. subtilis under the control of a strong constitutive promoter into (integrative) plasmids (e.g. pUC19, pBR322) [9–11].
In the framework of EU legislation [12, 13], companies wanting to market in the EU a specific additive produced by Genetically Modified Microorganisms (GMMs), like vitamin B2, need to submit an application that will be evaluated by the European Food Safety Authority (EFSA). This will result in a scientific opinion concerning the safety and the efficacy of the product EFSA . According to the EFSA guidance, food and feed additives produced by GMMs intended for human and animal consumption must be pure. This means that both GMMs (either alive or killed) and newly introduced genes should have been removed from the final product . This is especially important for safety evaluation as the amount of recombinant residual host cell DNA as well as the risk of gene transfer need to be assessed. Based on this information, EFSA will be able to deliver its scientific opinion concerning the product, including whether or not the production strain or its recombinant DNA was detected in the final product (for example see ).
According to EU legislation [16, 17], labelling is not required for vitamins and additives produced by a GMM because the producer strain (either alive or killed) and its components, including DNA, are no longer present in the commercialized product as the final product should be carefully purified, not containing any residual GMM material. Therefore, unlike for plant GMOs, the companies do not have to provide an event-specific method to trace the GMO in food and feed products on the EU market and there is no control by enforcement laboratories, for the correctness of the labelling. It is assumed that the company bringing the GMM-derived riboflavin to the market has verified the absence of the GMM or its recombinant DNA in the final product, so that only “pure” product can be found on the EU market. However, very recently, it was demonstrated that GMM-contamination in feed additives could occur as a viable GM -B. subtilis was found in an imported lot of vitamin B2 feed additive placed on the EU market. This has led in September 2014 to a notification for unauthorized GMO (UGM) in the European Rapid Alert System for Food and Feed (RASFF 2014–1249 ).
Confronted with this RASFF, as no official method is available for detecting this GMM in food and feed, two methods previously developed for other purposes have been proposed and used by enforcement laboratories controlling the GMO-content in food and feed at the EU market. The first one allows identifying the presence of B. subtilis strains and includes a PCR amplification of the 16S rRNA gene followed by its sequencing . This method is not especially convenient and is more laborious. Indeed, it requires sequencing in addition to the PCR analysis and this is not commonly used by the routine laboratories in GMO analysis. The second method allows detecting the presence of recombinant DNA. It is a real-time PCR method that was used as an internal amplification control targeting the plasmid pUC19, commonly used to construct genetically modified B. subtilis. This method was originally developed in duplex with a qPCR method for the detection of pathogenic Yersinia enterocolitica in food samples . It targets a standard vector construct (i.e. pUC19), and consequently has a large spectrum of targets commonly used in many GMM constructs. Moreover, traces of such kind of recombinant DNA could be found in mixes containing optimized Taq polymerase enzyme used for PCR. Therefore, these 2 methods both provide ambiguous results in the detection of GM-B. subtilis overproducing riboflavin. There is a need to develop a faster, specific method targeting this particular GM-strain to be put at the disposition of the enforcement laboratories.
In order to develop such a specific method, the GM-B. subtilis overproducing riboflavin was isolated from three different samples of imported vitamin B2 feed additive and collected by the French competent authorities in the framework of the RASFF 2014–1249. Next generation sequencing (NGS) was performed on the genomic DNA extracted from one of the isolates. The sequencing reads were de novo assembled and genome annotation was performed on the contig sequences . However, NGS data analysis often requires appropriate tools, involving bioinformatics expertise which is not always present in the average enforcement laboratory. This hampers the use of this technology to easily obtain critical sequence information in order to be able to develop a specific qPCR detection method.
In this study, the DNA sequences previously generated by the NGS approach  were used to develop a TaqMan® qPCR method targeting the junction between the B. subtilis riboflavin operon and the vector used to construct this GM-strain. Different performance criteria of the developed qPCR method such as specificity, sensitivity, PCR efficiency and repeatability were evaluated according to the GMO guidelines [22, 23]. Finally, the qPCR method was assessed directly on the gDNA extracts from the three different samples of imported vitamin B2 feed additive and collected by the French competent authorities. The method proposed in the present study provides a crucial tool for identifying specifically and rapidly the GM-B. subtilis overproducing riboflavin firstly detected in the imported vitamin B2 feed additive. Moreover, it illustrates that NGS data can be used very easily, without specific bioinformatics expertise, in order to obtain crucial information to develop specific TaqMan® qPCR methods, to detect unknown and unauthorized GMO in food and feed.
Results and Discussion
Sequence analysis to identify the junction and design of the GM-specific qPCR assay
The first step in the development of the specific qPCR assay consisted of identifying the junction of the GM-insert into the Bacillus subtilis GM-strains extracted from vitamin B2 samples 2014–3557 [Genbank: JYFL01000000] . As elaborated in Methods, a BLAST study of the contigs of the previously published sequence of the GM-Bacillus subtilis 2014–3557 isolated from vitamin B2 80 % [Genbank:JYFL01000000.1]  was used to identify the contigs containing the riboflavin biosynthesis operon ribGBAH, with a focus on locating the ribA gene as it has been previously reported that additional expression of this gene encoding the rate limiting enzyme in riboflavin synthesis increases riboflavin synthesis even more as compared to only overexpression of ribGBAH [7, 11, 24]. The overlapping contigs Contig0022 [Genbank:JYFL01000022.1] and Contig0016 [Genbank:JYFL01000016.1] were retrieved as containing genes having 100 % similarity with B. subtilis subsp. subtilis str. 168 genes ribH, ribBA, ribE and ribD, involved in riboflavin biosynthesis [7, 25], as it could also be deduced from the annotation of the GM-B. subtilis genome sequence [Genbank:JYFL01000000]  (Fig. 1). Using the publicly available genome sequence of B. subtilis subsp. subtilis str. 168 [Genbank:CP010052.1] and a BLAST analysis, Contig00019 was identified as the other contig adjacent to Contig00022 (Fig. 1). However, a gap of 37 basepairs between Contig00019 and Contig00022 was observed when comparing these contig sequences to the B. subtilis subsp. subtilis str. 168 genome sequence [Genbank:CP010052.1] . PCR and sequence analysis was used to confirm that the genomic regions present on both contigs are indeed adjacent in the GM-Bacillus genome (see Additional file 1: Figure S1 and Additional file 2: Figure S2). Further investigation of Contig00019 and Contig00016 revealed the presence of a region with 100 % similarity to the pSM19035 [Genbank:AY357120.1] or pBT233 [Genbank:NG_034603.1] plasmid (Fig. 1). It needs to be further investigated whether these plasmid sequences are integrated in the GM-Bacillus genome. Nevertheless, this vector has probably been used to construct the GM-strain. Therefore, the junctions between the non-naturally present vector sequences and the Bacillus region containing the riboflavin biosynthesis genes (GM-cassette) were to be targeted by the event-specific TaqMan® qPCR assay. The GM-cassette junctions were found around position 1454 in the contig0019 and around position 5887 in the contig0016 (Fig. 1). Primer pairs and probes have been designed at these positions and tested on the GM-Bacillus subtilis 2014–3557 (data not shown). One of the assays targeting the junction on the Contig0019 was kept for further evaluation (Table 1 and Fig. 1). This assay, called the VitB2-UGM qPCR assay, was then evaluated for its specificity.
Determination of VitB2-UGM qPCR assay’s specificity
The specificity test allows testing the inclusivity and exclusivity of the VitB2-UGM qPCR assay. The experimental design involved 3 Bacillus subtilis GM-strains overproducing riboflavin that were extracted from three samples of imported vitamin B2 80 %, 51 non-target strains representing 28 species belonging to 19 genera, a CTD (DNA diluent) and an NTC (Table 2). The non-target microorganisms to test the exclusivity were chosen among taxonomically closely related (e.g. other Bacillus species) and not closely related (pathogenic or not) bacteria that can be present in the environment and in food and feed matrices . In feed/food additives, normally there should not be any bacteria present, as elaborated above. However, some species also known to be used as GMM were included in the specificity test (e.g. E. coli, B. licheniformis). The VitB2-UGM qPCR assay amplified 100 % (3/3) of the tested GM-B. subtilis strains overproducing riboflavin and none (0/51) of the non-target samples nor the No Template Control (NTC) nor the diluent DNA (CTD) (Table 2). Thus, the detection of the GM-Bacillus subtilis 2014–3557 using the VitB2-UGM qPCR assay is 100 % specific for the GM-strain. This was somehow expected as the VitB2 UGM qPCR assay is a construct specific assay (designed on the junction between the between inserted endogenous riboflavin biosynthesis genes (“GM-cassette”) and non-naturally present sequences. The presence of this junction is not expected in wild-type bacteria, no matter what their origin or phylogenetic relationship is. The sequence of the amplicon was determined and corresponded to the expected one (data not shown).
Determination of VitB2-UGM qPCR assay’s dynamic range and PCR efficiency
The dynamic range of a qPCR assay is the range of concentrations where the assay performs linearly. This was assessed for the VitB2-UGM qPCR assay by the analysis in duplicate of a serial dilution of gDNA (10,000 to 0.01 theoretical genomic copies) of the GM-Bacillus subtilis strain 2014–3557 overproducing riboflavin . In addition, this analysis allowed for the assessment of the coefficient of determination (R2) and the PCR efficiency (E). The coefficient of determination (R2) is an indicator of how well the data fit the linear regression curve. The VitB2-UGM qPCR assay performed linearly between 1 and 10,000 genomic copies as its R2-value, i.e. 0.99, exceeds the required 0.98 [22, 23] (Fig. 2). From the dynamic range analysis, the PCR efficiency (E) was calculated. The VitB2-UGM qPCR assay displayed a PCR efficiency of 105.4 %, thereby corresponding to the accepted limits fixed for qPCR qualitative method used to detect GMO of 80–120 % as well as to the ones used to quantify GMO of 90–110 % [22, 23].
Determination of VitB2-UGM qPCR assay sensitivity and repeatability
The sensitivity test was performed to determine the LOD of the VitB2-UGM qPCR assay. The LOD is defined as the lowest concentration of an analyte giving a positive result with a probability of 95 % [22, 23]. The LOD of the VitB2-UGM qPCR assay was determined to be between 1 and 2 copies (Table 3) complying with the requirement “below 25 genomic copies” [23, 27]. The r-values of the Cq-values ranged between 1.3 to 3.2 Cq. The RSDr-values of the Cq-values were between 1.4 to 3.3 % (Table 3). There is no limit fixed for these criteria for qualitative qPCR methods .
Test of the VitB2-UGM qPCR assay on vitamin B2 80 % feed additive samples
To verify if the developed VitB2-UGM qPCR assay targeting the GM-Bacillus subtilis 2014–3557 overproducing riboflavin is able to detect this GM strain in real-life samples, the developed qPCR assay has been tested on gDNA extracts from the three samples of imported vitamin B2 80 % feed additive and containing these GM-Bacillus strains.
The DNA extracted from each of the samples gave a clear signal when using the VitB2-UGM qPCR assay at a Cq-value around 14 for 10 times dilutions and around 17.5 for 100 times dilutions (Table 4).
This demonstrates that the VitB2-UGM qPCR assay can be used to detect the GM-Bacillus strain in real-life samples.
As no method was available to specifically detect the presence of an unauthorized GMO consisting of a GM Bacillus subtilis bacterium overproducing vitamin B2 found in a feed additive (RASFF 2014–1249), in the present study a new TaqMan® qPCR method specifically detecting this GM vitamin B2 overproducing bacterium was designed. To be able to develop the VitB2-UGM qPCR assay, critical sequence information was retrieved and used from massive sequence data previously obtained by an NGS approach . To analyse these data, a very simple strategy based on the common BLAST tools, that can be used by any enforcement lab without profound bioinformatics expertise, was successfully followed. This allowed to design a specific qPCR method targeting the junction between inserted endogenous riboflavin biosynthesis genes (“GM-cassette”) and non-naturally present sequences (vector).
To guarantee the high-standard quality of the VitB2-UGM method, acceptance parameters like specificity, sensitivity, PCR efficiency and repeatability commonly used to validate a qualitative qPCR method for the detection of GMO were evaluated [22, 23, 27]. The method complies with the acceptance criteria for all of the assessed parameters. In addition, to further assess the applicability of the newly developed method not only for isolated strains but also for feed products, the VitB2-UGM method was successfully tested on gDNA immediately extracted from the feed additive from where the GM-strains were isolated. This confirms its effective use in routine analysis for reliable detection by enforcement laboratories of the GM- B. subtilis bacterium overproducing vitamin B2 (RASFF 2014–1249). In the future, this specific method can be used to screen feed and food products regarding this specific GM-bacterium overproducing vitamin B2. However, this method might have a broader application, as vitamin B2 is used in a wide range of products present at both the food and feed market . Moreover, this vitamin is not a unique example of additive produced by GMM. Indeed, GMM are used to produce other vitamins, additives and processing agents for the food industry . Companies are responsible for ensuring that their food and feed put at the market are safe, pure, and that it complies with legislation on food and feed additives and regulation on reducing or eliminating human health risks caused by possible contaminants. However, in view of the RASFF demonstrating the presence of an unauthorized GM (UGM) contaminant in feed, the intensive use of GMM to synthetize food and feed additives, as well the possible import of these products in EU, the question is raised if a more systematic survey of the EU market for GMM food and feed contaminants in additives is not necessary.
Hereto, our study underlines that when this NGS strategy is applied on unknown and unauthorized GMM, it delivers a massive amount of information that can be easily used by enforcement laboratories without specific bioinformatics expertise for the development of specific qPCR methods targeting these UGM in the food and feed market.
The bacterial strains used in this study are listed in Table 2. A panel of DNA extracted from 51 bacterial strains (24 Bacillus subtilis strains, 7 others species of the Bacillus genus (B.) and 20 strains from 18 other bacterial genera, all obtained from National Reference Centres and Laboratories and C. Nannan from UCL, Belgium) and 3 GM-B. subtilis strains has been tested. The 3 GM-B. subtilis strains (2014–3557, 2014–3558 and 2014–3559) were isolated from three samples of imported vitamin B2 80 % feed additive powder and analysed by the French GMO-Laboratory “Service commun des Laboratoires” in the framework of the RASFF 2014–1249. No ethical approval was required for this study as no humans, human data or animals were involved.
Bacterial growth conditions, DNA extraction quantification and sequencing
For the specificity test, DNA of the Bacillus strains (including the GM ones) was extracted by boiling as previously described . For the other strains, DNA was extracted using the “Gram-Negative or Gram-Positive Bacteria” protocol of the DNeasy Blood and Tissue Kit (Qiagen Benelux–B.V., KJ Venlo, the Netherlands) from the pellet of 2 ml of overnight cultures of each bacterial strain grown in Brain-Heart Infusion (BHI) broth at the adequate temperature and oxygen condition. For the sensitivity test, the DNA from the GM-B. subtilis 2014–3557 strain was extracted from the pellet of 84 ml of overnight culture under the conditions specified above, using the Genomic-tip 100/G (Qiagen, Benelux–B.V., KJ Venlo, the Netherlands). All kits were used according to the manufacturer’s recommendations. The DNA concentration was determined using a Nanodrop® 2000 device. The DNA quality was verified via the A260/A280 and A260/A230 ratio and on agarose gel (1 %). Sequencing was performed on a ABI3130xl Genetic Analyzer (Life Technologies, Gent, Belgium) following the manufacturer’s instructions.
Sequence analysis to identify the contig containing the junction
To identify which contig of the GM-B. subtilis 2014–3557 genome sequence [Genbank:JYFL00000000] ) contains the natural riboflavin biosynthesis genes (“GM-cassette”), a megablast of each of the 36 contigs was performed. Contigs matching with riboflavin biosynthesis genes were further analysed by megablast to determine the position of the junction between inserted natural riboflavin biosynthesis genes (“GM-cassette”) and the non-naturally present sequence (vector). These sequences were then used to develop the GM-B. subtilis 2014–3557 overproducing riboflavin-specific qPCR assay.
PCR analysis to confirm the order of the contigs Contig00019 and Contig00016
The PCR was performed in a final reaction volume of 25 μl containing 1× DreamTaq Green Master Mix (Thermo Fisher Scientific), 400 nM of forward (Scaf-19-F3-seq, 5′ TTTCGGTACTCAATCAGCTTTTC 3′) and reverse primer (Scaf-22-R-seq, 5′ CAAGGTGCTTCCTCCTTTAAT 3′), and 5 μl DNA template (corresponding to 106 genomic copies of the GM-B. subtilis 2014–3557 strain extracted with the Genomic-tip 100/G kit (Qiagen Benelux – B.V., KJ Venlo, the Netherlands). The following protocol was used in a thermal cycler: 2 min at 95 °C, 35 cycles of 30 s at 95 °C, 30 s at 55 °C and 1 min at 72 °C, 10 min at 72 °C. The PCR product (600 bp) was visualized by agarose gel electrophoresis (1 % agar) with ethidium bromide staining. The PCR product was subsequently purified with the Wizard PCR clean-up kit (Promega Benelux B.V., Leiden, the Netherlands). This DNA was used for bidirectional Sanger sequencing (on a ABI3130XL Genetic Analyzer instrument, according to the manufacturer’s instructions), using the primers Scaf-19-F3-seq and Scaf-22-R-seq. The obtained sequences were aligned to the corresponding one of the B. subtilis subsp. subtilis str. 168 genome sequence [Genbank:CP010052.1] using DNASTAR Lasergene software.
Development and in silico assessment of primer pairs
The qPCR assay developed makes use of an hydrolysis probe, also known as TaqMan® probe. The design of the primers and probe for the TaqMan® qPCR assay to detect the GM-junction was done using the “Primer 3” program [29, 30] with the “product size range” parameter set at “60 to 120 bp” and “primer size” optimal set at “22 bases” and a part of the DNA sequence of contig-19 (described in [Genbank: JYFL01000000], and containing the B. subtilis riboflavin biosynthesis genes (“GM-cassette”) and the non-naturally present sequences (vector) as input target sequence. Subsequently, the selectivity of the designed primer pairs were in silico tested. This test was performed using the “wprimersearch” software available on wEMBOSS Open Source Software package [31–33]. This software mimics the PCR amplification with the tested primers and using as template the DNA sequences present in a collection of bacterial DNA sequences of non-GM bacterial genomes of 248 strains, representing 122 species belonging to 72 genera, retrieved from the NCBI public database . Only primer pairs that did not yield an in silico amplicon using the DNA sequences of this collection of non-GM bacteria strains were retained for the following steps. As a last verification step, a megablast of the amplicon was performed to verify that no similar sequence exists [35–37].
VitB2-UGM qPCR assay
All qPCR assays were performed on an Applied Biosystems 7300 Real-Time PCR System (Applied Biosystems, Life Technologies, Gent, Belgium) using MicroAmp® Optical 96-Well Reaction Plates closed with MicroAmp® Optical 8-Cap Strips (Applied Biosystems, Life Technologies, Gent, Belgium). The reaction was performed in a final volume of 25 μl containing 5 μl of the appropriate template (104 copies of gDNA of each of the tested strains (Table 2) for the specificity test or a serial dilution of gDNA of GM-B. subtilis strain 2014–3557 overproducing riboflavin  for the sensitivity test), 1X TaqMan® PCR Mastermix (Diagenode, Liège, Belgium), 250 nM of each primer and 100nM of the probe (Table 1, i.e. VitB2-UGM qPCR assay). The estimated bacterial genomic copy number was calculated according to the genome size of each targeted bacteria (information available in the NCBI database; for Bacillus strains, the previously reported genome size of B. subtilis of 4, 214, 810 bp was used) using the formula previously presented . The following thermal program was applied: A single cycle of DNA polymerase activation for 10 min at 95 °C followed by 45 amplification cycles of 15 sec at 95 °C (denaturing step) and 1 min at 60 °C (annealing-extension step). The fluorescent reporter signal was normalized against the internal reference dye (ROX) signal and the threshold limit setting was performed in automatic mode, according to the ABI Sequence Detection Software version 1.4 (Applied Biosystems, Life Technologies, Gent, Belgium), unless manual adjustment was considered necessary. “No Template” Controls (NTC) using DNase and RNase free water (Acros, Geel, Belgium) and DNA diluent control using calf thymus DNA (CTD, Invitrogen) at 4 ng/μl were included in each assay to assess respectively primer dimer formation and non-specific amplification. Sanger sequencing on a ABI3130xl Genetic Analyzer (Life Technologies, Gent, Belgium) according to the manufacturer’s instructions, was used to confirm the obtained amplicon.
The qPCR assay gives one result, i.e. the quantification cycle (Cq) value which represents the fractional cycle at which the PCR amplification reaches the threshold level for the reaction .
Dynamic range and calculation of the PCR efficiency
The dynamic range was assessed for the VitB2-UGM qPCR assay by the analysis in duplicate of a serial dilution in a carrier DNA background (4 ng/μL Calf Thymus DNA (CTD) (Invitrogen, Life Technologies, Gent, Belgium)) of gDNA (10,000 to 0.01 theoretical genomic copies) of the GM-Bacillus subtilis strain 2014–3557 overproducing riboflavin . The carrier DNA avoids the improper dilution due to low concentration of gDNA. The PCR efficiency was calculated according to the formula previously reported .
Sensitivity test and repeatability calculation
The GM-B. subtilis strain 2014–3557 overproducing riboflavin  was used as target. To determine the LOD, a range between 10 and 0.1 theoretical genomic copies was tested (i.e. 10, 5, 2, 1 and 0.1). Each dilution was tested in six replicates. Moreover, the analysis was performed three times independently, under repeatable conditions, resulting in 18 repeats for each dilution point. In order to verify that the DNA concentration used to calculate the LOD was not overestimated, the number of positive signals was recorded over the 18 replicates. The “0.1 copy” was included to verify that the concentration of DNA used to calculate the LOD was approximately correct.
The repeatability of the assay can be evaluated thanks to the independent tests performed, i.e. with the same protocol, with the same samples, by the same operator using the same apparatus within a short interval of time . The repeatability limit (r) and the relative standard deviation of repeatability (RSDr) were calculated according to ISO 16140 . The RSDr and r-values of the Cq-values were calculated at each dilution point.
Extraction of genomic DNA from vitamin B2 samples
Genomic DNA (gDNA) was extracted from feed additive samples using a CTAB-based method adapted from Dellaporta et al. . Four volumes of CTAB extraction buffer (NaCl 1.4 M, EDTA 0.02 M, Tris–HCl 0.1 M, CTAB 2 %), supplemented with Ribonuclase A (at a final concentration of 0.1 mg/ml) were added to 250 mg of vitamin B2 powder, mixed regularly and incubated for 30 min at 65 °C. Next, Proteinase K (at a final concentration of 0.4 mg/ml) was added, mixed regularly and incubated for 30 min at 65 °C. After centrifugation (10 min at 12,000 g), the supernatant was collected and 1 volume of chloroform was added and vortexed for 30 sec. After mixing and centrifugation (15 min at 12,000 g), the upper phase was collected. 0.6 volumes of isopropanol and 0.4 volumes of ammonium acetate 10 M were added, gently mixed by inversion, incubated for 30 min on ice and centrifuged (15 min at 12,000 g). The supernatant was discarded and the pellet was washed with 1 ml of 70 % ethanol and subsequently centrifuged (15 min at 12,000 g). The pellet was dried for a few minutes at 37 °C. Finally, the pellet was re-suspended in 350 μl of TE (10 mM Tris-0.2 mM EDTA pH 8.0) pre-warmed at 50 °C for 2 min. gDNA extracts were then stored at +4 °C. For the qPCR analysis, 1/10 and 1/100 dilutions were made in TE (10 mM Tris-0.2 mM EDTA pH 8.0) and analysed in duplicate with the VitB2-UGM qPCR assay as described above.
Availability of supporting data
The data sets supporting the results of this article are included within the article and its additional files:
Additional file 1: Figure S1: Sequence of B. subtilis subsp. subtilis str. 168 genome [Genbank:CP010052.1] corresponding to the region containing part of Contig00019 and Contig00022, and the “37 bp gap”.
Additional file 2: Figure S2: Alignment of sequences of PCR fragment containing the “37 bp gap” between Contig00019 and Contig00022 and the B. subtilis subsp. subtilis str. 168 genome sequence [Genbank:CP010052.1].
Genetically modified micro-organism
Rapid alert system for food and feed
Next generation sequencing
Genetically modified organism
European food safety authority
No template control
European network of GMO laboratories
Unauthorized genetically modified organism
Calf Thymus DNA
- R2 :
Coefficient of determination
GMO Compass web site. http://www.gmo-compass.org/eng/home/. Accessed 18 March 2015.
Wang Z, Chen T, Ma X, Shen Z, Zhao X. Enhancement of riboflavin production with Bacillus subtilis by expression and site-directed mutagenesis of zwf and gnd gene from Corynebacterium glutamicum. Bioresour Technol. 2011;102:3934–40.
Burgess C, O’connell-Motherway M, Sybesma W, Hugenholtz J, van Sinderen D. Riboflavin production in Lactococcus lactis: potential for in situ production of vitamin-enriched foods. Appl Environ Microbiol. 2004;70:5769–77.
Stahmann KP, Revuelta JL, Seulberger H. Three biotechnical processes using Ashbya gossypii, Candida famata, or Bacillus subtilis compete with chemical riboflavin production. Appl Microbiol Biotechnol. 2000;53:509–16.
Lin Z, Xu Z, Li Y, Wang Z, Chen T, Zhao X. Metabolic engineering of Escherichia coli for the production of riboflavin. Microb Cell Fact. 2014;13:104.
Hao T, Han B, Ma H, Fu J, Wang H, Wang Z, et al. In silico metabolic engineering of Bacillus subtilis for improved production of riboflavin, Egl-237, (R, R)-2,3-butanediol and isobutanol. Mol Biosyst. 2013;9:2034–44.
Hümbelin M, Griesser V, Keller T, Schurter W, Haiker M, Hohmann HP, et al. GTP cyclohydrolase II and 3,4-dihydroxy-2-butanone 4-phosphate synthase are rate-limiting enzymes in riboflavin synthesis of an industrial Bacillus subtilis strain used for riboflavin production. J Ind Microbiol Biotech. 1999;22:1–7.
Schallmey M, Singh A, Ward OP. Developments in the use of Bacillus species for industrial production. Can J Microbiol. 2004;50:1–17.
Apeldoorn MEV, Van Steeg H, Speijers GJA. Riboflavin derived by fermentation with genetically modified Bacillus subtilis. 42 ed. WHO Food Additives Series: World Health Organization, Geneva and IPCS - International Programme on Chemical Safety. 1999. Available: http://www.inchem.org/documents/jecfa/jecmono/v042je05.htm
Härtl B, Wehrl W, Wiegert T, Homuth G, Schumann W. Development of a new integration site within the Bacillus subtilis chromosome and construction of compatible expression cassettes. J Bacteriol. 2001;183:2696–9.
Duan Y, Chen T, Chen X, Wang J, Zhao X. Enhanced riboflavin production by expressing heterologous riboflavin operon from B. cereus ATCC14579 in Bacillus subtilis. Chin J Chem Eng. 2010;18:129–36.
Commission of the European Communities. Regulation (EC) No 1333/2008 of the European Parliament and of the Council of 16 December 2008 on food additives. 2008.
Commission of the European Communities. Regulation (EC) No 1331/2008 of the European Parliament and of the Council of 16 December 2008 establishing a common authorisation procedure for food additives, food enzymes and food flavourings. 2008.
EFSA panel on genetically modified organisms (GMO). Scientific opinion on guidance on the risk assessment of genetically modified microorganisms and their products intended for food and feed use. 9(6):2193 ed. EFSA Journal. 2011:1–54.
EFSA Panel on Additives and Products or Substances used in Animal Feed (FEEDAP). Scientific opinion on the safety and efficacy of vitamin B2 (80 %) as riboflavin produced by Bacillus subtilis for all animal species, based on a dossier submitted by VITAC EEIG. 12(1):3531 ed. 2014:1–2.
Commission of the European Communities. Regulation (EC) No 1830/2003 of the European Parliament and of the Council of 22 September 2003 concerning the traceability and labelling of genetically modified organisms and the traceability of food and feed products produced from genetically modified organisms and amending Directive 2001/18/EC. 2003.
Commission of the European Communities. Regulation (EC) No 1829/2003 of the European Parliament and of the Council of 22 September 2003 on genetically modified food and feed. 2003.
Rapid Alert System for Food and Feed (RASFF). Notification details-2014.1249 unauthorised genetically modified (Bacillus subtilis) bacteria in vitamin B2 from China, via Germany. 2014.
Amtliche Sammlung von Untersuchungsverfahren nach 28b GenTG. Amplifizierung von Teilsequenzen des bakteriellen 16S-rRNA-Gens zur Gattungs- und Speziesidentifizierung. G 21.40-1 ed. 2010:1–5.
Mäde D, Reiting R, Strauch E, Ketteritzsch K, Wicke A. A real-time PCR for detection of pathogenic Yersinia enterocolitica in food combined with an universal internal amplification control system. J Verbr Lebensm. 2008;3:141–51.
Barbau-Piednoir E, De Keersmaecker SCJ, Wuyts V, Gau C, Pirovano W, Costessi A, et al. Genome sequence of the EU-unauthorized genetically modified Bacillus subtilis strain 2014–3557 overproducing riboflavin, isolated from an imported lot of Vitamin B2 80 % feed additive. Genome Announc. 2015;3:e00214–5.
European Network of GMO Laboratories (ENGL). Definition of Minimum Performance Requirements for Analytical Methods of GMO Testing. 2008:1–8.
European Network of GMO Laboratories (ENGL). Definition of minimum performance requirements for analytical methods of GMO testing. JRC Technical Report. 2015:1–24.
Perkins JB, Sloma A, Hermann T, Theriault K, Zachgo E, Erdenberger T, et al. Genetic engineering of Bacillus subtilis for the commercial production of riboflavin. J Ind Microbio Biotech. 1999;22:8–18.
Lim SH, Choi JS, Park EY. Microbial production of riboflavin using riboflavin overproducers, Ashbya gossypii, Bacillus subtilis, and Candida famate: an overview. Biotechnol Bioprocess Eng. 2001;6:75–88.
ISO: International Organization for Standardization. ISO 22118:2011-Microbiology of food and animal feeding stuffs - Polymerase chain reaction (PCR) for the detection and quantification of food-borne pathogens - Performance characteristics of molecular detection methods. 2011.
Broeders S, Huber I, Grohmann L, Berben G, Taverniers I, Mazzara M, et al. Guidelines for validation of qualitative real-time PCR methods. Trends Food Sci Tech. 2014;37:115–26.
Lin A, Kase JA, Moore MM, Son I, Tran N, Clotilde LM, et al. Multilaboratory validation of a Luminex microbead-based suspension array for the identification of the 11 most clinically relevant Shiga toxin-producing Escherichia coli O serogroups. J Food Prot. 2013;76:867–70.
Primer 3. http://frodo.wi.mit.edu/primer3/. Accessed January 2015.
Rozen S, Skaletsky H. Primer 3 on the WWW for general users and for biologist programmers. Methods Mol Biol. 2000;132:365–86.
wEMBOSS Open Source Software package. http://emboss.bioinformatics.nl/cgi-bin/emboss/primersearch.
Rice P, Longden I, Bleasby A. EMBOSS: the European Molecular Biology Open Software Suite. Trends Genet. 2000;16:276–7.
Sarachu M, Colet M. wEMBOSS: a web interface for EMBOSS. Bioinformatics. 2005;21:540–1.
NCBI public database. http://www.ncbi.nlm.nih.gov/sites/entrez.
Madden TL, Tatusov RL, Zhang J. Applications of network BLAST server. Methods Enzymol. 1996;266:131–41.
Altschul SF, Gish W, Miller W, Myers EW, Lipman DJ. Basic local alignment search tool. J Mol Biol. 1990;215:403–10.
Barbau-Piednoir E, Botteldoorn N, Yde M, Mahillon J, Roosens NH. Development and validation of qualitative SYBR®Green real-time PCR for detection and discrimination of Listeria spp. and Listeria monocytogenes. Appl Microbiol Biotechnol. 2013;97:4021–37.
Bustin SA. Absolute quantification of mRNA using real-time reverse transcription polymerase chain reaction assays. J Mol Endocrinol. 2000;25:169–93.
Rutledge RG, Cote C. Mathematics of quantitative kinetic PCR and the application of standard curves. Nucleic Acids Res. 2003;31, e93.
International Organization for Standardization (ISO). ISO 3534–1:1993 Statistics-Vocabulary and Symbols-part1. 1993.
ISO: International Organization for Standardization. ISO 16140:2003-Microbiology of food and animal feeding stuffs - Protocol for the validation of alternative methods. Edited by International Organization for Standardization (ISO). 2003.
Dellaporta S, Wood J, Hicks J. A plant DNA minipreparation: Version II. Plant Mol Biol Rep. 1983;1:19–21.
The authors would like to greatly thank Catherine Nannan from the Université Catholique de Louvain (Belgium) for providing the non-GM Bacillus subtilis strains. This work was supported by the Belgian Federal Public Service: Public Heath, Food Chain Safety and Environment.
The authors declare that they have no competing interests.
EBP performed the sequence analysis, designed the qPCR method, extracted the DNA, carried out the specificity test, analysed the data, made the tables and drafted part of the manuscript. SDK carried out part of the sequence analysis, participated in the study design and in the discussion of the data, made the Fig. 1, drafted part of the manuscript, prepared the manuscript for submission and assisted in writing of the rebuttal and revision of the manuscript after reviewing. MD carried out the sensitivity and dynamic range test. CG and PP isolated the GM-Bacillus strains and performed the qPCR on real life samples, and provided comments on the draft of the manuscript. NR conceived of the study, participated in its design and coordination, participated in the sequence analysis and in the discussion of the data, drafted part of the manuscript, and responded to the reviewers and revised the manuscript after review. All authors read and approved the final manuscript.
Sequence of B. subtilis subsp. subtilis str. 168 genome sequence [Genbank:CP010052.1] corresponding to the region containing part of Contig00019 and Contig00022. There is a gap of 37 basepairs (uncoloured region in the figure) between Contig00019 (pink region in figure) and Contig00022 (yellow region in figure), when aligning the obtained contigs to B. subtilis subsp. subtilis str. 168 genome sequence [Genbank:CP010052.1]. PCR and sequence analysis were used to confirm that the genomic regions present on both contigs are indeed adjacent in the GM-Bacillus genome. Hereto, primers Scaf-19-F3-seq (positioned on Contig00019, indicated in green) and Scaf-22-R-seq (positioned on Contig00022, indicated in blue) were used to amplify the flanking regions of the junction between Contig00019 and Contig00022. Subsequently, the obtained PCR fragment was sequenced. (PDF 101 kb)
Alignment of sequences of PCR fragment containing the “37 bp gap” between Contig00019 and Contig00022 and the B. subtilis subsp. subtilis str. 168 genome sequence [Genbank:CP010052.1]. Alignment of the obtained sequences (“57-scaf19-F3” and “57-scaf22-R”) of the PCR fragment containing the junction between Contig00019 and Contig00022 to the corresponding region in the B. subtilis subsp. subtilis str. 168 genome sequence [Genbank:CP010052.1] (region taken as displayed in Additional file 1: Figure S1). The nucleotides at the beginning of the obtained sequence are to be ignored, as these are of lower quality, which is known for Sanger sequencing with BigDye Terminator v3.1. The position corresponding to the “37 bp gap” has been indicated with a red box. It is clear that both contigs are adjacent in the GM-Bacillus genome, and that the sequence region corresponding to the “gap” was not present in the reads/contigs obtained through the the NGS analysis. (PDF 534 kb)
About this article
Cite this article
Barbau-piednoir, E., De Keersmaecker, S.C.J., Delvoye, M. et al. Use of next generation sequencing data to develop a qPCR method for specific detection of EU-unauthorized genetically modified Bacillus subtilis overproducing riboflavin. BMC Biotechnol 15, 103 (2015). https://doi.org/10.1186/s12896-015-0216-y
- Unauthorized GM-Bacillus subtilis
- Vitamin B2