# Targeted High-Throughput Sequencing Enables the Detection of Single Nucleotide Variations in CRISPR/Cas9 Gene-Edited Organisms

Marie-Alice Fraiture, Jolien D'aes, Emmanuel Guiderdoni, Anne-Cécile

Meunier, Thomas Delcourt, Stefan Hoffman, Els Vandermassen, Sigrid C J de
Keersmaecker, Kevin Vanneste, Nancy H. C. Roosens

► To cite this version:

Marie-Alice Fraiture, Jolien D'aes, Emmanuel Guiderdoni, Anne-Cécile Meunier, Thomas Delcourt, et al.. Targeted High-Throughput Sequencing Enables the Detection of Single Nucleotide Variations in CRISPR/Cas9 Gene-Edited Organisms. Foods, 2023, 12 (3), pp.455. 10.3390/foods12030455. hal-04021719

HAL Id: hal-04021719

https://hal.inrae.fr/hal-04021719

Submitted on 9 Mar 2023

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
Article

# Targeted High-Throughput Sequencing Enables the Detection of Single Nucleotide Variations in CRISPR/Cas9 Gene-Edited Organisms

Marie-Alice Fraiture <fnref n="1" />,† ID, Jolien D'aes <fnref n="1" />,† ID, Emmanuel Guiderdoni <fnref n="2" />,<fnref n="3" /> ID, Anne-Cécile Meunier <fnref n="2" />,<fnref n="3" />, Thomas Delcourt <fnref n="1" />, Stefan Hoffman <fnref n="1" />, Els Vandermassen <fnref n="1" />, Sigrid C. J. De Keersmaecker <fnref n="1" /> ID, Kevin Vanneste <fnref n="1" />,‡ and Nancy H. C. Roosens <fnref n="1" />,*,‡ ID

<figure><img src="image_2.png" /></figure>

check for
updates

**Citation:** Fraiture, M.-A.; D'aes, J.; Guiderdoni, E.; Meunier, A.-C.; Delcourt, T.; Hoffman, S.; Vandermassen, E.; De Keersmaecker, S.C.J.; Vanneste, K.; Roosens, N.H.C. Targeted High-Throughput Sequencing Enables the Detection of Single Nucleotide Variations in CRISPR/Cas9 Gene-Edited Organisms. *Foods* **2023**, *12*, 455. https://doi.org/10.3390/foods12030455

Academic Editor: Diego Bonatto

Received: 16 December 2022

Revised: 9 January 2023

Accepted: 16 January 2023

Published: 18 January 2023

<figure><img src="image_3.png" /></figure>

**Copyright:** © 2023 by the authors.
Licensee MDPI, Basel, Switzerland.
This article is an open access article
distributed under the terms and
conditions of the Creative Commons
Attribution (CC BY) license (https://
creativecommons.org/licenses/by/
4.0/).

**Abstract:** Similar to genetically modified organisms (GMOs) produced by classical genetic engineering, gene-edited (GE) organisms and their derived food/feed products commercialized on the European Union market fall within the scope of European Union Directive 2001/18/EC. Consequently, their control in the food/feed chain by GMO enforcement laboratories is required by the competent authorities to guarantee food/feed safety and traceability (2003/1829/EC; 2003/1830/EC). However, their detection is potentially challenging at both the analytical and interpretation levels since this requires methodological approaches that can target and detect a specific single nucleotide variation (SNV) introduced into a GE organism. In this study, we propose a targeted high-throughput sequencing approach, including (i) a prior PCR-based enrichment step to amplify regions of interest, (ii) a sequencing step, and (iii) a data analysis methodology to identify SNVs of interest. To investigate if the performance of this targeted high-throughput sequencing approach is compatible with the performance criteria used in the GEO detection field, several samples containing different percentages of a GE rice line carrying a single adenosine insertion in OsMADS26 were prepared and analyzed. The SNV of interest in samples containing the GE rice line could successfully be detected, both at high and low percentages. No impact related to food processing or to the presence of other crop species was observed. The present proof-of-concept study has allowed us to deliver the first experimental-based evidence indicating that the proposed targeted high-throughput sequencing approach may constitute, in the future, a specific and sensitive tool to support the safety and traceability of the food/feed chain regarding GE plants carrying SNVs.

**Keywords:** CRISPR/Cas9; single nucleotide variations; transgenic rice; detection; PCR-based enrichment; high-throughput sequencing; food and feed chain; genetically modified organism

## 1. Introduction

Gene-editing techniques such as CRISPR/Cas9 have revolutionized the genetic manipulations of organisms, allowing the introduction of a precise modification (substitution, insertion, or deletion) in the genome of an organism of interest. Such modifications include both large nucleotide sequence variations as well as single nucleotide variations (SNVs) [1–7]. Given the endless possibilities offered by these techniques, the commercialization of such gene-edited (GE) organisms by biotech companies is expected to increase in the next few years. Currently, a GE high-GABA (gamma-aminobutyric acid) tomato has been commercialized in Japan, while two GE crops, a herbicide-tolerant canola and a

<footnote n="1">¹ Sciensano, Transversal Activities in Applied Genomics (TAG), rue Juliette Wytsman 14, 1050 Brussels, Belgium</footnote>

<footnote n="2">² CIRAD, UMR AGAP Institut, F-34398 Montpellier, France</footnote>

<footnote n="3">³ UMR AGAP Institut, Université de Montpellier, CIRAD, INRAE, Institut Agro, F-34398 Montpellier, France</footnote>

<footnote>* Correspondence: nancy.roosens@sciensano.be</footnote>

<footnote>† These authors contributed equally to this work.</footnote>

<footnote>‡ These authors contributed equally to this work.</footnote>
soybean with a modified oil composition, have been commercialized in the United States of America [1-7].

Since 2018, the commercialization of GE organisms and their derived food/feed products on the European Union market has fallen within the scope of European Union Directive 2001/18/EC, similar to genetically modified organisms (GMOs) produced by classical genetic engineering. To guarantee the traceability of such products in the food/feed chain, control by the competent authorities using enforcement laboratories is therefore requested [2,4,7-12]. For GE organisms carrying a large nucleotide sequence variation, the development of detection methods is highly similar to the detection of GMOs produced by classical genetic engineering. However, for GE organisms carrying an SNV of only one or a few nucleotides, their detection and discrimination from the parental lines are not trivial, both at the analytical and interpretation levels. At the analytical level, the design of a detection method restricted to a particular region carrying such an SNV of interest can be challenging. At the interpretation level, successful detection of an SNV does not automatically prove the use of gene-editing techniques because an SNV of only one or few nucleotides could also have occurred naturally or through random mutagenesis techniques [2,4,7,12].

Among the available detection methods, the potential use of PCR-based methods, such as real-time PCR and digital PCR technologies that are widely used by GMO enforcement laboratories for the detection of GMOs produced by classical genetic engineering, has already previously explored and has been assessed as suitable to analytically detect crops carrying an SNV introduced by gene editing [3,13-16]. However, without information demonstrating that the SNV is only present in a particular GE organism, the presence of this GE organism cannot be unambiguously deduced from a positive PCR signal associated with this SNV [2,3,13-17]. Therefore, in addition to detecting the SNV(s) introduced by gene-editing techniques, the detection of additional targets is required to identify the specific line that was GE. These targets can include, for example, certain nucleotide variations associated specifically with a particular genotype or cultivar as well as potential off-target mutations introduced through the gene-editing process. Unintended changes to the genome can indeed occur using gene-editing techniques, such as CRISPR/Cas9, at sites with high sequence similarity to the on-target sequence, although gene-editing processes are generally designed to minimize them and their frequency is usually low [2,18,19]. However, the design and development of real-time PCR and digital PCR methods to simultaneously detect multiple targets of interest can quickly become fastidious, as multiplexing with these technologies is usually limited to only a few targets [20-22]. Therefore, such PCR-based technologies, although widely used in GMO enforcement laboratories for the detection of classical GMOs, do not necessarily represent the most appropriate strategy in this context.

A targeted high-throughput sequencing strategy can overcome these bottlenecks as it offers both sequence information and a high potential for multiplexing by simultaneously targeting multiple samples as well as several SNVs of interest [2,18,23-25]. This sequencing strategy requires a prior enrichment step of the principal targets using probe-based capture hybridization or PCR-based amplification. The latter enrichment approach is usually considered the most efficient and presents the advantage of avoiding potential difficulties related to the design of probes specific to SNVs. Targeted high-throughput sequencing strategies have previously been applied in different fields, including the genotyping of plants based on key SNVs [2,18,23-25]. The potential ability of high-throughput sequencing, including targeted strategies, to be used for the detection of GMOs produced by classical genetic engineering has also been previously illustrated [26-39]. However, in contrast to PCR-based methods (e.g., conventional PCR, real-time PCR, digital droplet PCR (ddPCR)), high-throughput sequencing-based methods still require further harmonization and standardization to be used in the context of GMO control. Nowadays, in the GMO detection field, no minimum quality performance criteria are established for sequencing-based methods [40]. In addition, to our knowledge, no targeted high-throughput sequencing strategy has currently yet been experimentally investigated for the detection of SNVs carried by
GE organisms, leading to a lack of information about the performance of the methodology and the potential compatibility of such methodology with the European GMO regulations, including the threshold of 0.9% for the labeling of authorized GMOs and the technical zero-tolerance threshold of 0.1% for the low-level presence of GMOs in feed. In the case of non-European countries, for example, in Japan, a threshold of 5% is established for GMO labeling [9,41,42].

In this proof-of-concept study, for the first time, a targeted high-throughput sequencing strategy has been developed to detect, in multiple samples in parallel, an SNV of interest carried by a GE plant, and the performance of this approach is assessed. A GE rice line carrying a single adenosine insertion in *OsMADS26* was used. *OsMADS26*, encoding a MADS-box transcription factor, was previously reported as putatively involved in rice stress responses [16]. Several samples mimicking various contaminations of this GE rice line in food products were prepared to assess the feasibility and performance of the proposed targeted high-throughput sequencing approach (e.g., specificity, sensitivity, applicability). The generated sequencing data were analyzed using an in-house pipeline developed to identify SNVs, including both single polymorphism nucleotides and indels, as well as to estimate their frequency. The potential added value of the proposed targeted high-throughput sequencing to support GMO enforcement laboratories in the control of GE organisms in the food/feed chain is discussed in light of the experimental-based results obtained in this proof-of-concept study.

## 2. Materials and Methods

### 2.1. Plant Materials

Rice (*Oryza sativa* L. *Nipponbare*) seeds from a GE line and its parental line were used in this study. The GE rice line, carrying a single adenosine insertion in *OsMADS26* (locus: Os08g02070) encoding a MADS-box transcription factor, was generated using CRISPR/Cas9 and is currently not commercialized on the market. The homozygous adenosine insertion was integrated into the coding region close to the start codon between genomic positions 679646 and 679647 on chromosome VIII. This frameshift mutation inactivates *OsMADS26*, consequently putatively increasing biotic resistance and biotic stress tolerance (Rice SNP-Seek Database; Meunier et al., unpublished) [16,43]. Using rice seeds from the GE line and its parental line, rice noodles were prepared in-house as described previously [16]. DNA from homogenous powders of ground rice seeds and rice noodles was extracted using a CTAB-based procedure [44] in combination with the Genomic-tip20/G kit (QIAGEN, Hilden, Germany), as described previously [45,46].

Certified reference materials of wild-type (WT) maize (*Zea mays*) and WT soybean (*Glycine max*) were collected from the Institute for Reference Materials and Measurements (IRMM, Geel, Belgium). DNA from these crop species was obtained as described previously [47].

Using the Nanodrop® 2000 (ThermoFisher, Wilmington, DE, USA) device, the DNA concentration was measured by spectrophotometry. DNA purity was considered as complying with the acceptance criteria (A260/A280 ration ~1.8 and A260/A230 ratio ~2.0–2.2).

### 2.2. Sample Preparation

To assess the performance of the proposed targeted high-throughput sequencing approach, several samples were prepared. For samples n°1–11 (Table 1), DNA from rice seeds of the GE rice line and its parental line were mixed to create samples containing 100%, 99.9%, 99.1%, 95%, 90%, 50%, 10%, 5%, 0.9%, 0.1%, or 0% of the GE rice line, ranging from ~14,000 to 0 estimated haploid genome copies [16] (Table S1).
<table><caption>**Table 1.** Observed frequency of the CRISPR/Cas9 adenosine insertion, specific to the GE rice line in rice seed samples (n°1–11). The samples were composed of DNA extracted from rice seeds of the GE rice line and its parental rice line. The samples contained 100%, 99.9%, 99.1%, 95%, 90%, 50%, 10%, 5%, 0.9%, 0.1%, or 0% of the GE rice line. For each sample, the copy number associated with the GE rice line, previously measured using a 2-plex ddPCR method [16], is indicated (Table S1). The targeted high-throughput sequencing approach was performed in quadruplicate for each sample, allowing us to subsequently determine the allele frequency of the CRISPR/Cas9 adenosine insertion specific to the GE rice line. The associated depth of coverage of the reference at the position of this specific indel (bp 679,646) is indicated.</caption><thead><tr><th rowspan="3">Sample n°</th><th colspan="3">GE Rice Line Content</th><th colspan="3">Allele Frequency ~</th><th colspan="2">Depth of Coverage ~</th></tr><tr><th colspan="2">Percentage</th><th>Detection</th><th colspan="2">Percentage</th><th>StDev</th><th>Value</th><th>StDev</th></tr><tr><th>100</th><th>13,733.3</th><th>+</th><th>99.92</th><th>0.03</th><th></th><th>119,592</th><th>29,499</th></tr></thead><tbody><tr><td>1</td><td>99.9</td><td>13,346.7</td><td>+</td><td>99.86</td><td>0.03</td><td></td><td>97,038</td><td>25,007</td></tr><tr><td>2</td><td>99.1</td><td>13,273.3</td><td>+</td><td>99.73</td><td>0.02</td><td></td><td>70,873</td><td>12,339</td></tr><tr><td>3</td><td>95</td><td>12,113.3</td><td>+</td><td>97.97</td><td>0.15</td><td></td><td>101,581</td><td>63,687</td></tr><tr><td>4</td><td>90</td><td>10,753.3</td><td>+</td><td>92.06</td><td>0.24</td><td></td><td>75,787</td><td>36,739</td></tr><tr><td>5</td><td>50</td><td>5640.0</td><td>+</td><td>53.80</td><td>0.57</td><td></td><td>71,791</td><td>19,099</td></tr><tr><td>6</td><td>10</td><td>1153.3</td><td>+</td><td>8.01</td><td>0.29</td><td></td><td>86,841</td><td>10,531</td></tr><tr><td>7</td><td>5</td><td>528.0</td><td>+</td><td>4.00</td><td>0.15</td><td></td><td>98,283</td><td>40,588</td></tr><tr><td>8</td><td>0.9</td><td>85.3</td><td>+</td><td>0.66</td><td>0.01</td><td></td><td>101,464</td><td>33,345</td></tr><tr><td>9</td><td>0.1</td><td>10.8</td><td>+</td><td>0.13*</td><td>0.02</td><td></td><td>126,724</td><td>42,631</td></tr><tr><td>10</td><td>0</td><td>0</td><td>+</td><td>0.00</td><td>0.00</td><td></td><td>158,436</td><td>99,002</td></tr></tbody></table>

~ For each sample, the average value of the 4 replicates and the associated standard deviation (StDev) are indicated.

* The CRISPR/Cas9 adenosine insertion was observed in 3 of the 4 replicates.

For samples n°12–17 (Table 2), DNA from in-house-prepared rice noodles of the GE rice line and its parental line were mixed to obtain samples containing 100%, 99.9%, 99.1%, 0.9%, 0.1%, or 0% of the GE rice line, ranging from ~14,000 to 0 estimated haploid genome copies [16] (Table S1).

<table><caption>**Table 2.** Observed frequency of the CRISPR/Cas9 adenosine insertion, specific to the GE rice line in rice noodle samples (n°12–17). The samples were composed of DNA extracted from rice noodles of the GE rice line and its parental rice line. The samples contained 100%, 99.9%, 99.1%, 0.9%, 0.1%, or 0% of the GE rice line. For each sample, the copy number associated with the GE rice line, measured using a 2-plex ddPCR method [16], is indicated (Table S1). The targeted high-throughput sequencing approach was performed in quadruplicate for each sample, allowing us to subsequently determine the allele frequency of the CRISPR/Cas9 adenosine insertion specific to the GE rice line. The associated depth of coverage of the reference at the position of this specific indel (bp 679,646) is indicated.</caption><thead><tr><th rowspan="3">Sample n°</th><th colspan="3">GE Rice Line Content</th><th colspan="3">Allele Frequency ~</th><th colspan="2">Depth of Coverage ~</th></tr><tr><th colspan="2">Percentage</th><th>Detection</th><th colspan="2">Percentage</th><th>StDev</th><th>Value</th><th>StDev</th></tr><tr><th>100</th><th>13,610.0</th><th>+</th><th>99.85</th><th>0.06</th><th></th><th>120,711</th><th>60,941</th></tr></thead><tbody><tr><td>12</td><td>99.9</td><td>13,596.4</td><td>+</td><td>99.68</td><td>0.15</td><td></td><td>47,261</td><td>27,646</td></tr><tr><td>13</td><td>99.1</td><td>13,487.0</td><td>+</td><td>99.05</td><td>0.16</td><td></td><td>69,972</td><td>26,283</td></tr><tr><td>14</td><td>0.9</td><td>122.5</td><td>+</td><td>1.03</td><td>0.03</td><td></td><td>62,408</td><td>16,122</td></tr><tr><td>15</td><td>0.1</td><td>13.6</td><td>+</td><td>0.24</td><td>0.04</td><td></td><td>68,359</td><td>8826</td></tr><tr><td>16</td><td>0</td><td>0</td><td>+</td><td>0.00</td><td>0.00</td><td></td><td>67,906</td><td>37,437</td></tr></tbody></table>

~ For each sample, the average value of the 4 replicates and the associated standard deviation (StDev) are indicated.
For samples n°18–19 (Table 3), DNA from rice seeds of the GE rice line and its parental line, both at 14 estimated haploid genome copies, were mixed with DNA from WT maize or WT soybean at 13,972 estimated haploid genome copies [16] (Table S1).

**Table 3.** Observed frequency of the CRISPR/Cas9 adenosine insertion, specific to the GE rice line in crop mixture samples (n°18–19). (A) The samples were composed of 0.1% of the GE rice line (corresponding to 14 estimated haploid genome copies), 0.1% of the parental rice line (corresponding to 14 estimated haploid genome copies), and 99.8% of WT maize (sample n°12) or WT soybean (sample n°13) (corresponding to 13,972 estimated haploid genome copies). (B) For each sample, the copy number associated with the GE rice line, previously measured using a 2-plex ddPCR method [16], is indicated (Table S1). Additionally, regarding the rice ingredient, the rice ratio of the GE rice line is specified. The targeted high-throughput sequencing approach was performed in quadruplicate for each sample, allowing us to subsequently determine the allele frequency of the CRISPR/Cas9 adenosine insertion specific to the GE rice line. The associated depth of coverage of the reference at the position of this specific indel (bp 679,646) is indicated.

(A)

Sample Composition

Sample n°18

0.1 % GE rice line (14 estimated haploid genome copies)

0.1 % parental rice line (14 estimated haploid genome copies)

99.8 % WT maize (113,972 estimated haploid genome copies)

Sample n°19

0.1 % GE rice line (14 estimated haploid genome copies)

0.1 % parental rice line (14 estimated haploid genome copies)

99.8 % WT soybean (113,972 estimated haploid genome copies)

(B)

Crop Mixture Sample Description

CRISPR/Cas9 Adenosine Insertion

<table><thead><tr><th rowspan="2">Sample n°</th><th colspan="3">GE Rice Line Content</th><th colspan="2">Allele Frequency ~</th><th colspan="2">Depth of Coverage ~</th></tr><tr><th>Copy Number</th><th>Rice Ratio *</th><th>Detection</th><th>Percentage</th><th>StDev</th><th>Value</th><th>StDev</th></tr></thead><tbody><tr><td>18</td><td>10.6</td><td>50</td><td>+</td><td>40.46</td><td>10.78</td><td>533</td><td>50</td></tr><tr><td>19</td><td>9.3</td><td>50</td><td>+</td><td>48.52</td><td>18.12</td><td>23,124</td><td>3181</td></tr></tbody></table>

~ For each sample, the average value of the 4 replicates and the associated standard deviation (StDev) are indicated.

* The rice ratio of the GE rice line is determined according to the total amount of rice (ingredient) present in the investigated samples.

The calculation of the estimated haploid genome copy number was based on the size of the rice (0.5 pg), maize (2.6 pg), and soybean (1.16 pg) haploid genomes [48,49].

### 2.3. *ddPCR Assays*

A recently developed 2-plex ddPCR method targeting the GE rice line was applied in triplicate to all samples (n°1–19), as described previously [16]. For each ddPCR assay, an NTC (no template control) was included. On this basis, the copy number of the GE rice line was measured (Table S1).

### 2.4. *Conventional PCR Assays*

A standard 25 µL reaction volume was used, containing 1x KAPA HiFi HotStart ReadyMix (Roche, Brussels, Belgium), 900 nM of each primer (Eurogentec, Liège, Belgium), and 5 µL of DNA. The primer pair (OsMADS26-F: GACAGGAGGAGAGGAGGAAGA; OsMADS26-R: AAGGTGACCTGACGGTGAAC) used in this study was previously designed to cover 113 bp, for the GE rice line, and 112 bp, for the parental rice line, of the OsMADS26 region carrying the CRISPR/Cas9 SNV (Fraiture et al., 2022). In this study, these primers were supplemented by Illumina® overhang adapter sequences (OsMADS26-F: TCGTCGGCAGCGTCAGATGTGTATAAGAGACAG-[GACAGGAGGAGAGGAGGA AGA]; OsMADS26-R: GTCTCGTGGGCTCGGAGATGTGTATAAGAGACAG-[AAGGTGA
CCTGACGGTGAAC]) in order to be compatible with the subsequent Illumina® MiSeq sequencing run (Illumina support instructions). The PCR assays were performed on a T100TM Thermal Cycler (Bio-Rad, Temse, Belgium). The PCR program consisted of a single cycle at 95 °C for 3 min (initial denaturation), 30 cycles at 95 °C for 30 s (denaturation), at 60 °C for 30 s (annealing), and at 72 °C for 30 s (extension), and a single cycle at 72 °C for 5 min (final extension). The final PCR products were visualized by electrophoresis using the Tapestation 4200 device with the associated D1000 Screen Tape and reagents (Agilent, Machelen, Belgium) (Figure S1). Each sample (n°1–19) and the included NTC (no template control) were tested in quadruplicate.

## 2.5. *Library Preparation, Sequencing, and Data Analysis*

Following the manufacturer's instructions, each PCR product was purified using Agencourt® Ampure® XP (Beckman Coulter, Danvers, MA, USA), and amplicon sequencing libraries were prepared according to Illumina's instructions (Illumina, San Diego, CA, USA) (Illumina support instructions). Then, 15% PhiX was added to the library pool. Each of the 80 PCR products, generated from the 20 different samples (n°1–20), was individually barcoded and then pooled for sequencing on the same run. Sequencing was carried out on an Illumina® MiSeq system using V3 chemistry, obtaining 250 bp paired-end reads.

The quality of raw sequencing data was evaluated using FastQC 0.11.5 (www.bioinfor-matics.babraham.ac.uk/projects/fastqc, accessed on 14 August 2021) with default settings (Figure S2). The frequency of the GE rice line in the samples was estimated with an in-house pipeline, consisting of the alignment of the amplicon sequencing reads to a rice reference genome, followed by variant calling to detect variants and estimate their allelic frequency. The reference for the alignment, the reference genome of the *Oryza sativa* Japonica Group (RefSeq GCF_001433935.1), was indexed with Samtools 1.9 [50] and Bowtie2 2.3.4.3 [51], while Picard 2.8.14 was used to generate a dictionary of the indexed reference FASTA file. The raw sequencing data were pre-processed using Trimmomatic 0.38 [52] with the following settings: ILLUMINACLIP:NexteraPE-PE.fa:2:30:10, LEADING:10, TRAILING:10, SLIDINGWINDOW:4:20, MINLEN: 40, after which the trimmed reads were aligned to the reference with Bowtie2 with the "-end-to-end" and "-very-sensitive" settings. The resulting alignments in SAM format were converted to BAM format with Samtools and sorted using Picard [53] with the option "SORT_ORDER = coordinate", followed by assignment of the reads to a new read group with Picard, with the flags "LB", "PL", "PU", and "SM" set to the arbitrary placeholder value "test." The resulting BAM files were indexed using Samtools and used as input for the creation of a target intervals file and indel realignment using GATK 3.7 [54], with "-maxReads" set to 1,000,000. Samtools was used to index the generated BAM files, followed by the insertion of indel qualities into the BAM files with LoFreq 2.1.3.1 [55]. Next, variants were called with the LoFreq package using the options "-call-indels", "-min-bq 35", "-min-mq 42", and "-sig 0.001". Additionally, the output calls were limited to positions 679,589–679,659 of chromosome 8 of the reference (NC_029263.1), which represents the amplicon sequence targeted by the PCR assay described in Section 2.3, excluding the primer sequences. No off-target mutation was observed in the target amplicons. The resulting VCF file was filtered with LoFreq, setting the minimal allelic frequency to 0.01 %. All samples were analyzed in this study with a cut-off for the allelic frequency at 0.1%.

# 3. Results

## 3.1. *Development of a Workflow for Targeted High-Throughput Sequencing*

The possibility of using a targeted high-throughput sequencing approach for the detection of an SNV in CRISPR/Cas9 plants has been investigated in this study. As a proof-of-concept, a GE rice line carrying a homozygous single adenosine insertion in *OsMADS26* and its parental rice line were used. Using CRISPR/Cas9, this SNV was introduced in the coding region close to the start codon, resulting in a gene inactivation that was expected to increase the biotic resistance and biotic stress tolerance of the GE rice line (Meunier et al.,
unpublished) [16,56]. In this study, a targeted high-throughput sequencing approach has been investigated for the detection of the single nucleotide insertion carried by the GE rice line. This sequencing approach is composed of three main successive steps (Figure 1).

<figure><img src="image_4.png" /><figcaption>**Figure 1.** Schematic workflow of the targeted high-throughput sequencing approach performed in this study. A total of 20 different samples were analyzed using the proposed targeted high-throughput</figcaption></figure>
sequencing approach, including 19 samples containing various percentages of the GE rice line as well as 1 NTC. The proposed targeted high-throughput sequencing approach is composed of three successive steps. First, a PCR-based enrichment step is applied to the whole DNA extracted from a given sample. The used primer pair was previously designed to amplify the OsMADS26 region carrying an SNV of interest (insertion or deletion of a single adenosine), resulting in a PCR amplicon of 113 bp for the GE rice line and a PCR amplicon of 112 bp for its parental rice line [16]. For each of these 20 samples, the PCR assay was performed in quadruplicate, resulting in a total of 80 final PCR products. Second, each PCR product is individually coupled to a unique barcode and then pooled for high-throughput sequencing. Third, raw sequencing reads are analyzed using an in-house pipeline based on LoFreq, developed to identify SNVs as well as to estimate their frequency.

First, a PCR-based enrichment step is applied to the whole DNA extracted from a given food matrix. A primer pair (OsMADS26-F and OsMADS26-R), previously designed to amplify the OsMADS26 region carrying the SNV of interest, was used (Fraiture et al., 2022) to create two different PCR amplicons: one with a length of 113 bp for the GE rice line and one with a length of 112 bp for the parental rice line (Figure 1). The specificity of this primer pair has been demonstrated previously [16]. Among the entire NCBI nucleotide (nr/nt) collection, only hits with sequences belonging to rice species were observed for these PCR amplicons [16]. Moreover, at the experimental level, PCR amplification was only observed with rice materials, both using WT and transgenic rice lines, and no PCR amplification was observed for non-targeted DNA, including various animal, microbial, and plant materials [16]. Second, each PCR product from each sample is individually coupled to a unique barcode during the sequencing library preparation steps, and all barcoded PCR products are subsequently pooled for high-throughput sequencing. The multiplexing of several samples, as illustrated in this study, can also be used to target additional SNVs of interest if these are known and/or available (e.g., species/genotype-specific sites and off-targeting signatures). Finally, based on the associated unique barcode, all raw sequencing data are demultiplexed. For each PCR product, the raw sequencing reads are analyzed using an in-house pipeline that compares each nucleotide from the generated amplicons to a rice reference genome, allowing us to (i) identify SNVs, including both single nucleotide polymorphisms and small indels, and (ii) estimate the frequency of the SNVs of interest. This pipeline was adapted from a previously developed pipeline to detect low-frequency variants in SARS-CoV-2 [57]. The pipeline input consists of the demultiplexed raw sequencing data, which are first processed to filter out low-quality reads and trim low-quality sequences from the reads. The resulting processed reads of the target amplicon are then aligned to the corresponding target region in the reference. To determine the frequency at which the SNV of interest is present in the aligned reads, variants are called for the target region with very strict quality cut-offs to retain only high-confidence variant calls and filtering to remove variants present at a frequency below 0.1%, reflecting the required threshold for the detection of GMOs in food/feed samples.

In this study, the multiplexing and performance of the proposed targeted high-throughput sequencing approach to detect the CRISPR/Cas9 adenosine insertion found in the GE rice line are explored. To this end, as schematized in Figure 1, several samples were prepared, containing different percentages of the GE rice line and its parental line, and subsequently pooled to be processed in the same sequencing run.

## 3.2. Assessment of Sensitivity

To assess whether the proposed targeted high-throughput sequencing approach is able to detect GMO contaminations at trace levels, its sensitivity was investigated using samples n°1–11. For these mixture samples, DNA from rice seeds of the GE rice line and its parental line (DNA extraction yield: ~6 × 10³ ng/g of rice) were mixed to obtain samples containing 100%, 99.9%, 99.1%, 95%, 90%, 50%, 10%, 5%, 0.9%, 0.1%, or 0% of the GE rice line (Table 1). For each sample, the copy number associated with the GE rice line was
measured by ddPCR (Tables 1 and S1). These mixture samples were specifically composed to mimic different contamination levels of a GE plant in its WT background.

The presence of the single variation of interest was not detected in sample n°11, which was exclusively composed of the parental rice line, whereas this SNV was observed in all samples containing the GE rice line (n°1–10), both at high and low percentages (Table 1). This SNV was detected in samples containing as low as 0.9% for all four replicates of the GE rice line, corresponding to approximately 85 haploid genome copies. The SNV of interest was also observed for the majority of replicates in sample n°10, containing 0.1% of the GE rice line, corresponding to approximately 11 haploid genome copies. In addition, the observed frequency of the detected single variation of interest was close to the expected value for all tested samples.

Based on these results, using the proposed targeted high-throughput sequencing approach, with an observed depth of coverage varying between 70,873 and 158,436 for the analyzed samples (n°1–11), the detection of the SNV of interest was possible even at low contamination levels (as low as 11 haploid genome copies and 0.1% of transgenic material per ingredient), which is crucial for GMO enforcement laboratories. The proposed targeted high-throughput sequencing approach was also able to deal with samples containing less than 25 copies of the target, which is one of the minimum performance criteria for GMO detection methods [40]. However, although the proposed targeted high-throughput sequencing approach is compatible with the sensitivity performance criteria used in the GMO detection field, such performance criteria are currently only intended for PCR-based detection methods (e.g., conventional PCR, real-time PCR, ddPCR), and no criteria have currently been established for sequencing-based detection methods. In addition, despite the initial objective of assessing the sensitivity of the proposed targeted high-throughput sequencing approach, the generated results also allowed us to provide preliminary data confirming the possible detection of heterozygous mutations. For example, sample n°6 may also be used to mimic a GE rice line carrying a heterozygous mutation (Table 1).

### 3.3. Assessment of Applicability

The applicability of the proposed targeted high-throughput sequencing approach was investigated using two different types of samples.

Firstly, samples n°12–17 were prepared to mimic different contamination levels of the GE rice line in its parental rice line within a processed food matrix, such as rice noodles. Food processing, defined as any physical, chemical, or mechanical food manipulations from the raw material to the final product, may induce DNA damage, such as the fragmentation of high molecular weight DNA. In order to mimic processed samples with strong DNA degradation, rice noodle materials were used in this study [58]. These rice noodle samples contained 100%, 99.9%, 99.1%, 0.9%, 0.1%, or 0% of the GE rice line (Table 2). For each sample, the copy number associated with the GE rice line was measured by ddPCR (Tables 2 and S1). The depth of coverage of the processed rice samples varied between 47,261 and 120,711. The SNV of interest was detected in all rice noodle samples containing the GE rice line, both at high and low percentages. The SNV was observed in all four replicates, with a frequency of the GE rice line as low as 0.1%, corresponding to approximately 14 haploid genome copies. Regarding the observed frequency of the SNV of interest, the values were close to the expected ones. In addition, the results obtained from the rice noodle samples (n°12–17) were generally similar to the results observed for the unprocessed rice seed samples (n°1–11), supporting that the proposed targeted high-throughput sequencing approach does not appear to be impacted by the DNA degradation of the analyzed samples.

Secondly, samples n°18–19 were prepared to mimic low-level contamination of both the GE rice line and its parental line in another crop species, such as maize or soybean (Table 3). More precisely, these samples contained 0.1% of the GE rice line (corresponding to 14 estimated haploid genome copies), 0.1% of the parental rice line (corresponding to 14 estimated haploid genome copies), and 99.8% of WT maize or WT soybean (corresponding to 13,972 estimated haploid genome copies). For each sample, the copy number
associated with the GE rice line was measured by ddPCR (Tables 3 and S1). Although the observed depth of coverage of these crop mixture samples was lower than for the other analyzed samples, the presence of the SNV of interest was detected in the four replicates, with an estimated frequency close to the expected value. These results, therefore, similarly support that the proposed targeted high-throughput sequencing approach is not affected by the presence of a high percentage of untargeted materials from other plant species.

## 4. Discussion

Despite their increasing use to strengthen the safety and traceability of the food/feed chain at several levels [59], sequencing-based methods are currently not widely used by GMO enforcement laboratories. Further harmonization and standardization are indeed still required to be used in the context of GMO control. The proposed targeted high-throughput sequencing approach developed in this proof-of-concept study allowed us to generate data in order to contribute to the assessment of minimum quality performance criteria for the detection of GE plants using sequencing-based methods. To our knowledge, the present study is the first available study using a targeted high-throughput sequencing approach to detect SNVs of interest in a GE organism. As a proof-of-concept, a CRISPR/Cas9 GE rice line carrying a single adenosine insertion, currently not commercialized on the market, was used. Using several samples containing different percentages of the GE rice line, the feasibility and performance of this targeted high-throughput sequencing approach were investigated, for the first time, as a tool to enable GMO control, guaranteeing the safety and traceability of the food/feed chain. The targeted high-throughput sequencing approach was able to detect the SNV of interest, even at the low level of contamination by the GE rice line of 0.1%. Moreover, although more conditions will need to be tested, this targeted high-throughput sequencing approach does not seem to be affected by food processing and is compatible with a low contamination level of the GE rice line in another plant species. Furthermore, although this targeted high-throughput sequencing approach was used as a qualitative detection method, the observed frequency of the single nucleotide insertion of interest was generally close to the amount of the GE rice line contained in each tested sample. Consequently, the proposed targeted high-throughput sequencing approach was demonstrated to be a promising tool to detect specific SNVs in GE organisms as well as to estimate their frequency.

One of the main advantages of the proposed targeted high-throughput sequencing approach is its high-multiplexing capacity, as illustrated in this study through the simultaneous sequencing of several different samples. This high-multiplexing capacity could be used to detect, in parallel, additional variations introduced by different routes (e.g., natural occurrences, random mutagenesis, or gene editing) as well as potential off-targeting sites. Key targets could also be selected in order to identify the specific line used to introduce modifications through gene-editing techniques. Although it will not completely solve the bottlenecks related to the unambiguous discrimination of one GE line from all other lines, high-throughput sequencing represents, however, non-negligible help in solving this problem. In this context, precise and detailed sequence information on the genomic background is necessary to identify additional key targets, such as SNVs specifically associated with a particular genotype, cultivar, or variety. However, this approach requires access to appropriate databases with high-quality sequence data in order to compare a given sample to reference genome sequences, allowing us to identify specific single mutation points and genetic elements within a particular genetic background. Ideally, such database(s) will need to represent the entire diversity of all genotypes, cultivars, and varieties for each species of interest to aid the specific identification of a particular GE organism. Nonetheless, such crucial information is currently not available. Moreover, even if the specific GE line could be identified, other genetic modifications may be naturally introduced during vegetative propagation or during crossing with other cultivars. Therefore, it would be very difficult to trace genotypes associated with the offspring of GE plants if the sequences from these offspring are not included in reference databases [2, 17, 18, 60]. This represents a major
bottleneck for the detection of GE organisms by enforcement laboratories, making the data analysis even more complex. This problem is increased with GE plants carrying SNVs by reason of divergences in international legislation to determine if a GE organism falls under the scope of GMO legislation or not [2,61,62]. Consequently, the traceability and associated sequence information will be even more complex to obtain from GE organisms considered as GMOs in Europe but not in other countries.

Although also subject to the aforementioned constraints regarding database reference information, shotgun metagenomics constitutes an interesting alternative to the targeted high-throughput sequencing approach because it requires no prior enrichment and it allows us to sequence an entire sample without any prior knowledge. However, no investigation has currently been performed to detect GM plants in the food/feed chain [63].

## 5. Conclusions

For the first time, the present proof-of-concept study has delivered experimental-based results on the performance of a targeted high-throughput sequencing approach to detect a specific SNV introduced into a CRISPR/Cas9 GE plant. These results highlight the potential performance compatibility of the proposed targeted high-throughput sequencing approach with the detection of GE plants in the GMO field. However, although this approach is highly promising, its current implementation in the context of GMO control remains challenging. Therefore, it will be essential for future case studies using this approach to generate sufficient data for further performance assessments. In addition, it is important to explore how this approach will be able to deal with major technical and analytical bottlenecks associated, for example, with the large size and complexity of plant genomes or with complex food/feed products. Moreover, if a GE organism is commercialized on the European market, a full validation process of the proposed sequencing approach, including transferability and robustness assessments, will be initially necessary before its use for control by enforcement laboratories.

**Supplementary Materials:** The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/foods12030455/s1, Table S1: Copy number associated with the GE rice line in samples n°1-19.; Figure S1: Visualization of final PCR products from samples n°1-19 and NTC (no template control); Figure S2: FastQC of the raw sequencing data. Reference [16] is cited in the supplementary materials.

**Author Contributions:** Conceptualization, M.-A.F. and N.H.C.R.; methodology, M.-A.F. and S.C.J. D.K.; formal analysis, M.-A.F., J.D., S.H., E.V., S.C.J.D.K. and K.V.; investigation, M.-A.F., J.D. and K.V.; writing—original draft preparation, M.-A.F. and J.D.; writing—review and editing, M.-A.F., J.D., E.G., A.-C.M., T.D., S.H., E.V., S.C.J.D.K., K.V. and N.H.C.R.; supervision, N.H.C.R.; funding acquisition, N.H.C.R. All authors have read and agreed to the published version of the manuscript.

**Funding:** The research that yielded these results was funded by the Belgian Federal Public Service of Health, Food Chain Safety and Environment through the contract GENEDIT (RF 20/6342).

**Data Availability Statement:** The original contributions presented in the study are included in the article/supplementary material. Raw sequencing data are deposited in the European Nucleotide Archive under study accession number PRJEB58426. Further inquiries can be directed to the corresponding author.

**Acknowledgments:** The sequencing was performed using the Transversal Activities in Applied Genomics Service at Sciensano. The GE rice line was produced in the frame of the ANR-11-BTBR-0001_GENIUS project.

**Conflicts of Interest:** The authors declare no conflict of interest.

## References

1. ISAAA. Available online: https://www.isaaa.org/kc/cropbiotechupdate/article/default.asp?ID=19024 (accessed on 23 October 2022).
2. Grohmann, L.; Keilwagen, J.; Duensing, N.; Dagand, E.; Hartung, F.; Wilhelm, R.; Bendiek, J.; Sprink, T. Detection and Identification of Genome Editing in Plants: Challenges and Opportunities. *Front. Plant Sci.* **2019**, *10*, 236. [CrossRef] [PubMed]

3. Chhalliyil, P.; Ilves, H.; Kazakov, S.A.; Howard, S.J.; Johnston, B.H.; Fagan, J.A. Real-Time Quantitative PCR Method Specific for Detection and Quantification of the First Commercialized Genome-Edited Plant. *Foods* **2020**, *9*, 1245. [CrossRef] [PubMed]

4. Menz, J.; Modrzejewski, D.; Hartung, F.; Wilhelm, R.; Sprink, T. Genome Edited Crops Touch the Market: A View on the Global Development and Regulatory Environment. *Front. Plant Sci.* **2020**, *11*, 586027. [CrossRef]

5. Tyagi, S.; Kumar, R.; Das, A.; Won, S.Y.; Shukla, P. CRISPR-Cas9 system: A genome-editing tool with endless possibilities. *J. Biotechnol.* **2022**, *10*, 36–53. [CrossRef] [PubMed]

6. Zhang, D.; Hussain, A.; Manghwar, H.; Xie, K.; Xie, S.; Zhao, S.; Larkin, R.M.; Qing, P.; Jin, S.; Ding, F. Genome editing with the CRISPR-Cas system: An art, ethics and global regulatory perspective. *Plant Biotechnol. J.* **2020**, *18*, 1651–1669. [CrossRef]

7. Ribarits, A.; Eckerstorfer, M.; Simon, S.; Stepanek, W. Genome-Edited Plants: Opportunities and Challenges for an Anticipatory Detection and Identification Framework. *Foods* **2021**, *10*, 430. [CrossRef]

8. Directive 2001/18/EC of the European Parliament and of the Council of 12 March 2001 on the deliberate release into the environment of genetically modified organisms and repealing Council Directive 90/220/EEC. *Off. J. Eur. Commun.* **2001**, L106, 1–38.

9. Regulation (EC) No 1829/2003 of the European Parliament and of the Council of 22 September 2003 on genetically modified food and feed. *Off. J. Eur. Union* **2003**, L268, 1–23.

10. 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. *Off. J. Eur. Union* **2003**, L268, 24–28.

11. European Court of Justice C-528/16-Judgement of 25 July 2018 on New Mutagenesis Techniques. Available online: http://curia.europa.eu/juris/document/document.jsf?text=&docid=204387&pageIndex=0&doclang=EN&mode=1st&dir=&occ=first&part=1&cid=138460 (accessed on 23 October 2022).

12. ENGL. Evaluation of the Scientific Publication "A Real-Time Quantitative PCR Method Specific for Detection and Quantification of the First Commercialized Genome-Edited Plant" P. Chhalliyil et al. in: *Foods* **2020**, *9*, 1245. Available online: https://gmo-crl.jrc.ec.europa.eu/ENGL/docs/ENGL%20Evaluation%20of%20the%20scientific%20publication%2002-10-2020.pdf (accessed on 22 October 2022). [CrossRef]

13. Peng, C.; Wang, H.; Xu, X.; Wang, X.; Chen, X.; Wei, W.; Lai, Y.; Liu, G.; Godwin, I.D.; Li, J.; et al. High throughput detection and screening of plants modified by gene editing using quantitative real-time polymerase chain reaction. *Plant J.* **2018**, *95*, 557–567. [CrossRef]

14. Peng, C.; Zheng, M.; Ding, L.; Chen, X.; Wang, X.; Feng, X.; Wang, J.; Xu, J. Accurate Detection and Evaluation of the Gene-Editing Frequency in Plants Using Droplet Digital PCR. *Front. Plant Sci.* **2020**, *11*, 610790. [CrossRef] [PubMed]

15. Zhang, H.; Li, J.; Zhao, S.; Yan, X.; Si, N.; Gao, H.; Li, Y.; Zhai, S.; Xiao, F.; Wu, G.; et al. An Editing-Site-Specific PCR Method for Detection and Quantification of CAO1-Edited Rice. *Foods* **2021**, *10*, 1209. [CrossRef] [PubMed]

16. Fraiture, M.A.; Guiderdoni, E.; Meunier, A.C.; Papazova, N.; Roosens, N.H.C. ddPCR strategy to detect a gene-edited plant carrying a single variation point: Technical feasibility and interpretation issues. *Food Control* **2022**, *137*, 108904. [CrossRef]

17. ENGL. Detection of Food and Feed Plant Products Obtained by New Mutagenesis Techniques. Available online: https://gmo-crl.jrc.ec.europa.eu/doc/JRC116289-GE-report-ENGL.pdf (accessed on 23 October 2022).

18. Shillito, R.D.; Whitt, S.; Ross, M.; Ghavami, F.; De Vleesschauwer, D.; D'Halluin, K.; Van Hoecke, A.; Meulewaeter, F. Detection of genome edits in plants—From editing to seed. *In Vitro Cell. Dev. Biol. Plant* **2021**, *57*, 595–608. [CrossRef]

19. Sturme, M.H.J.; van der Berg, J.P.; Bouwman, L.M.S.; De Schrijver, A.; de Maagd, R.A.; Kleter, G.A.; Battaglia-de Wilde, E. Occurrence and Nature of Off-Target Modifications by CRISPR-Cas Genome Editing in Plants. *ACS Agric. Sci. Technol.* **2022**, *2*, 192–201. [CrossRef]

20. Fraiture, M.A.; Herman, P.; Taverniers, I.; De Loose, M.; Deforce, D.; Roosens, N.H.C. Current and new approaches in gmo detection: Challenges and solutions. *Biomed. Res. Int.* **2015**, *2015*, 392872. [CrossRef]

21. Whale, A.S.; Huggett, J.F.; Tzonev, S. Fundamentals of multiplexing with digital PCR. *Biomol. Detect. Quantif.* **2016**, *10*, 15–23. [CrossRef]

22. Grohmann, L.; Barbante, A.; Eriksson, R.; Gatto, F.; Georgieva, T.; Huber, I.; Hulin, J.; Köppel, R.; Marchesi, U.; Marmin, L.; et al. Guidance Document on Multiplex Real-Time PCR Methods; EUR 30708 EN; Publications Office of the European Union: Luxembourg, 2021.

23. Shirasawa, K.; Kuwata, C.; Watanabe, M.; Fukami, M.; Hirakawa, H.; Isobe, S. Target Amplicon Sequencing for Genotyping Genome-Wide Single Nucleotide Polymorphisms Identified by Whole-Genome Resequencing in Peanut. *Plant Genome* **2016**, *9*, 3. [CrossRef]

24. Onda, Y.; Takahagi, K.; Shimizu, M.; Inoue, K.; Mochida, K. Multiplex PCR Targeted Amplicon Sequencing (MTA-Seq): Simple, Flexible, and Versatile SNP Genotyping by Highly Multiplexed PCR Amplicon Sequencing. *Front. Plant Sci.* **2018**, *9*, 201. [CrossRef]

25. Jo, J.; Kim, Y.; Kim, G.W.; Kwon, J.K.; Kang, B.C. Development of a Panel of Genotyping-in-Thousands by Sequencing in Capsicum. *Front. Plant Sci.* **2021**, *12*, 769473. [CrossRef] [PubMed]
26. Kovalic, D.; Garnaat, C.; Guo, L.; Yan, Y.; Groat, J.; Silvanovich, A.; Ralston, L.; Huang, M.; Tian, Q.; Christian, A.; et al. The use of next generation sequencing and junction sequence analysis bioinformatics to achieve molecular characterization of crops improved through modern biotechnology. *Plant Genome* **2022**, *5*, 149–163. [CrossRef]

27. Wahler, D.; Schäuser, L.; Bendiek, J.; Grohmann, L. Next generation sequencing as a tool for detailed molecular characterisation of genomic insertions and flanking regions in genetically modified plants: A pilot study using a rice event unauthorised in the EU. *Food Anal. Methods* **2013**, *6*, 1718–1727. [CrossRef]

28. Liang, C.; van Dijk, J.P.; Scholtens, I.M.; Staats, M.; Prins, T.W.; Voorhuijzen, M.M.; da Silva, A.M.; Arisi, A.C.; den Dunnen, J.T.; Kok, E.J. Detecting authorized and unauthorized genetically modified organisms containing vip3A by real-time PCR and next-generation sequencing. *ABC* **2014**, *406*, 2603–2611. [CrossRef] [PubMed]

29. Holst-Jensen, A.; Spilsberg, B.; Arulandhu, A.J.; Kok, E.; Shi, J.; Zel, J. Application of whole genome shotgun sequencing for detection and characterization of genetically modified organisms and derived products. *ABC* **2016**, *408*, 4595–4614. [CrossRef] [PubMed]

30. Willems, S.; Fraiture, M.A.; Deforce, D.; De Keersmaecker, S.; De Loose, M.; Ruttink, T.; Van Nieuwerburgh, F.; Roosens, N. Statistical framework for detection of genetically modified organisms based on Next Generation Sequencing. *Food Chem.* **2016**, *192*, 788–798. [CrossRef] [PubMed]

31. Fraiture, M.A.; Herman, P.; Papazova, N.; De Loose, M.; Ruttink, T.; Roosens, N.H.C. An integrated strategy combining DNA walking and NGS to detect GMO. *Food Chem.* **2017**, *232*, 351–358. [CrossRef]

32. Fraiture, M.A.; Herman, P.; De Loose, M.; Debode, F.; Roosens, N.H.C. How can we better detect unauthorized GMO in the food and feed chain. *Trends Biotechnol.* **2017**, *35*, 508–517. [CrossRef]

33. Bogožalec Košir, A.; Arulandhu, A.J.; Voorhuijzen, M.M.; Xiao, H.; Hagelaar, R.; Staats, M.; Costessi, A.; Žel, J.; Kok, E.J.; van Dijk, J.P. ALF: A strategy for identification of unauthorized GMOs in complex mixtures by a GW-NGS method and dedicated bioinformatics analysis. *Sci. Rep.* **2017**, *7*, 14155. [CrossRef]

34. Fraiture, M.A.; Saltykova, A.; Hoffman, S.; Winand, R.; Deforce, D.; Vanneste, K.; De Keersmaecker, S.C.J.; Roosens, N.H.C. Nanopore sequencing technology: A new route for the fast detection of unauthorized GMO. *Sci. Rep.* **2018**, *8*, 7903. [CrossRef]

35. Debode, F.; Hulin, J.; Charloteaux, B.; Coppieters, W.; Hanikenne, M.; Karim, L.; Berben, G. Detection and identification of transgenic events by next generation sequencing combined with enrichment technologies. *Sci. Rep.* **2019**, *9*, 15595. [CrossRef] [PubMed]

36. Fraiture, M.A.; Ujhelyi, G.; Ovesná, J.; Van Geel, D.; De Keersmaecker, S.C.J.; Saltykova, A.; Papazova, N.; Roosens, N.H.C. MinION sequencing technology to characterize unauthorized GM petunia plants circulating on the European Union market. *Sci. Rep.* **2019**, *9*, 7141. [CrossRef] [PubMed]

37. Boutigny, A.L.; Fioriti, F.; Rolland, M. Targeted MinION sequencing of transgenes. *Sci. Rep.* **2020**, *10*, 15144. [CrossRef] [PubMed]

38. Fraiture, M.A.; Papazova, N.; Vanneste, K.; De Keersmaecker, S.C.J.; Roosens, N.H.C. *GMO Detection and Identification Using Next-generation Sequencing in DNA Techniques to Verify Food Authenticity: Applications in Food Fraud*; Burns, M., Foster, L., Walker, M., Eds.; Royal Society of Chemistry Publishing: Cambridge, UK, 2020; pp. 96–106.

39. Saltykova, S.; Van Braekel, J.; Papazova, N.; Fraiture, M.A.; Deforce, D.; Vanneste, K.; De Keersmaecker, S.C.J.; Roosens, N.H.C. Detection and identification of authorized and unauthorized GMOs using high-throughput sequencing with the support of a sequence-based GMO database. *Food Chem. Mol. Sci.* **2022**, *4*, 100096. [CrossRef] [PubMed]

40. ENGL. Definition of Minimum Performance Requirements for Analytical Methods of GMO Testing. Available online: https://gmo-crl.jrc.ec.europa.eu/guidance-documents (accessed on 2 November 2022).

41. Regulation (EC) No 619/2011 of 24 June 2011 laying down the methods of sampling and analysis for the official control of feed as regards presence of genetically modified material for which an authorisation procedure is pending or the authorisation of which has expired. *Off. J. Eur. Union* **2011**, L166, 9–15.

42. Takabatake, R.; Egi, E.; Soga, K.; Narushima, J.; Yoshiwa, S.; Shibata, N.; Nakamura, K.; Kondo, K.; Kishine, M.; Mano, J.; et al. Development and Interlaboratory Validation of a Novel Reproducible Qualitative Method for GM Soybeans Using Comparative Cq-Based Analysis for the Revised Non-GMO Labeling System in Japan. *Anal. Chem.* **2022**, *94*, 13447–13454. [CrossRef] [PubMed]

43. Mansueto, L.; Fuentes, R.R.; Borja, F.N.; Detras, J.; Abriol-Santos, J.M.; Chebotarov, D.; Sanciangco, M.; Palis, K.; Copetti, D.; Poliakov, A.; et al. Rice SNP-seek database update: New SNPs, indels, and queries. *Nucleic Acids Res.* **2017**, *45*, D1075–D1081. [CrossRef]

44. ISO 21571:2005. Foodstuffs—Methods of Analysis for the Detection of Genetically Modified Organisms and Derived Products—Nucleic acid Extraction. ISO: Genève, Switzerland, 2015.

45. EURL. Sampling and DNA extraction of cotton seeds. Report from the Validation of the "CTAB/Genomic-tip 20" method for DNA extraction from ground cotton seeds. Available online: http://gmo-crl.jrc.ec.europa.eu/summaries/281-3006%20Cotton_D_NAExtr.pdf (accessed on 7 October 2022).

46. Fraiture, M.A.; Herman, P.; Taverniers, I.; De Loose, M.; Deforce, D.; Roosens, N.H.C. An innovative and integrated approach based on DNA walking to identify unauthorised GMOs. *Food Chem.* **2014**, *147*, 60–69. [CrossRef]

47. Broeders, S.; Fraiture, M.A.; Vandermassen, E.; Delvoye, M.; Barbau-Piednoir, E.; Lievens, A.; Roosens, N.H.C. New qualitative trait-specific SYBR®Green qPCR methods to expand the panel of GMO screening methods used in the CoSYPS. *Eur. Food Res. Technol.* **2015**, *241*, 275–287. [CrossRef]
48. Arumuganathan, K.; Earle, E.D. Nuclear DNA content of some important plant species. *Plant Mol. Biol. Rep.* **1991**, *9*, 208–218.
[CrossRef]

49. Leitch, I.J.; Johnston, E.; Pellicer, J.; Hidalgo, O.; Bennett, M.D. The Plant DNA C-values database (release 7.1): An updated online repository of plant genome size data for comparative studies. *New Phytol.* **2019**, *226*, 301–305.

50. Li, H.; Handsaker, B.; Wysoker, A.; Fennell, T.; Ruan, J.; Homer, N.; Marth, G.; Abecasis, G.; Durbin, R. The Sequence Alignment/Map Format and SAMtools. *Bioinformatics* **2009**, *25*, 2078–2079. [CrossRef] [PubMed]

51. Langmead, B.; Salzberg, S.L. Fast gapped-read alignment with Bowtie 2. *Nat. Methods* **2013**, *9*, 357–359. [CrossRef] [PubMed]

52. Bolger, A.M.; Lohse, M.; Usadel, B. Trimmomatic: A flexible trimmer for Illumina sequence data. *Bioinformatics* **2014**, *30*, 2114–2120.
[CrossRef] [PubMed]

53. Picard Toolkit. Broad Institute, GitHub Repository. 2019. Available online: https://broadinstitute.github.io/picard/ (accessed on 12 March 2022).

54. McKenna, A.; Hanna, M.; Banks, E.; Sivachenko, A.; Cibulskis, K.; Kernytsky, A.; Garimella, K.; Altshuler, D.; Gabriel, S.; Daly, M.; et al. The Genome Analysis Toolkit: A MapReduce framework for analyzing next-generation DNA sequencing data. *Genome Res.* **2014**, *20*, 1297–1303. [CrossRef]

55. Wilm, A.; Aw, P.P.K.; Bertrand, D.; Yeo, G.H.T.; Ong, S.H.; Wong, C.H.; Khor, C.C.; Petric, R.; Hibberd, M.L.; Nagarajan, N. LoFreq: A sequence-quality aware, ultra-sensitive variant caller for uncovering cell-population heterogeneity from high-throughput sequencing datasets. *Nucleic Acids Res.* **2012**, *40*, 11189–11201. [CrossRef]

56. Khong, N.G.; Pati, P.K.; Richaud, F.; Parizot, B.; Bidzinski, P.; Mai, C.D.; Bes, M.; Bourrié, I.; Meynard, D.; Beeckman, T.; et al. OsMADS26 negatively regulates resistance to pathogens and drought tolerance in rice. *Plant Physiol.* **2015**, *169*, 2935–2949.
[CrossRef]

57. Van Poelvoorde, L.A.E.; Delcourt, T.; Coucke, W.; Herman, P.; De Keersmaecker, S.C.J.; Saelens, X.; Roosens, N.H.C.; Vanneste, K. Strategy and performance evaluation of low-frequency variant calling for SARS-CoV-2 using targeted deep Illumina sequencing. *Front. Microbiol.* **2021**, *12*, 747458. [CrossRef]

58. Fraiture, M.A.; Herman, P.; Taverniers, I.; De Loose, M.; Van Nieuwerburgh, F.; Deforce, D.; Roosens, N.H.C. Validation of a sensitive DNA walking strategy to characterise unauthorised GMOs using model food matrices mimicking common rice products. *Food Chem.* **2015**, *173*, 1259–1265. [CrossRef]

59. Jagadeesan, B.; Gerner-Smidt, P.; Allard, M.W.; Leuillet, S.; Winkler, A.; Xiao, Y.; Chaffron, S.; Van Der Vossen, J.; Tang, S.; Katase, M.; et al. The use of next generation sequencing for improving food safety: Translation into practice. *Food Microbiol. 2019*, *79*, 96–115. [CrossRef]

60. Zhu, H.; Misel, L.; Graham, M.; Robinson, M.L.; Liang, C. CT-Finder: A Web Service for CRISPR Optimal Target Prediction and Visualization. *Sci Rep.* **2016**, *6*, 25516. [CrossRef] [PubMed]

61. Entine, J.; Felipe, M.S.S.; Groenewald, J.H.; Kershen, D.L.; Lema, M.; McHughen, A.; Nepomuceno, A.L.; Ohsawa, R.; Ordonio, R.L.; Parrott, W.A.; et al. Regulatory approaches for genome edited agricultural plants in select countries and jurisdictions around the world. *Transgenic Res.* **2021**, *30*, 551–584. [CrossRef] [PubMed]

62. Zimny, T.; Sowa, S. Potential effects of asymmetric legal classification of gene edited plant products in international trade, from the perspective of the EU. *EFB Bioecon. J.* **2021**, *1*, 100016. [CrossRef]

63. D'aes, J.; Fraiture, M.A.; Bogaerts, B.; De Keersmaecker, S.C.J.; Roosens, N.H.C.; Vanneste, K. Metagenomic Characterization of Multiple Genetically Modified Bacillus Contaminations in Commercial Microbial Fermentation Products. *Life* **2022**, *12*, 1971. [CrossRef] [PubMed]

**Disclaimer/Publisher's Note:** The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.