Next-generation sequencing as a tool for the molecular characterisation and risk assessment of genetically modified plants: Added value or not?

Technical Note on the quality of DNA sequencing for the molecular characterisation of genetically modified plants

As part of the risk assessment (RA) requirements for genetically modified (GM) plants, according to Regulation (EU) No 503/2013 and the EFSA guidance on the RA of food and feed from GM plants (EFSA GMO Panel 2011), applicants need to perform a molecular characterisation of the DNA sequences inserted in the GM plant genome. This Technical Note to the applicants puts together requirements and recommendations for the quality assessment of the methodology, analysis and reporting when DNA sequencing is used for the molecular characterisation of GM plants. In particular, it applies to the use of Sanger sequencing and next-generation sequencing for the characterisation of the inserted genetic material and its flanking regions at each insertion site, the determination of the copy number of all detectable inserts and the analysis of the genetic stability of the inserts. This updated document replaces the EFSA 2018 Technical Note and reflects the current knowledge in scientific-technical methods for generating and verifying, in a standardised manner, DNA sequencing data in the context of RA of GM plants. It does not take into consideration the verification and validation of the detection method which remains under the remit of the Joint Research Centre (JRC).

Development and evaluation of triple gene transgenic cotton lines expressing three genes (Cry1Ac-Cry2Ab-EPSPS) for lepidopteran insect pests and herbicide tolerance

Cotton is an international agricultural commodity and the main cash crop of Pakistan of which quality and quantity are subject to various whims of nature. Climate change, insect pest complex, and weeds are reducing its productivity. Here, we have developed triple gene cotton containing EPSPS gene along with two Bt toxin genes Cry1Ac and Cry2Ab using a strategy where all three genes are cloned in the same T-DNA, followed by successful cotton transformation via Agrobacterium-mediated transformation. This strategy has been developed to help cotton breeders in developing new cultivars by incorporating these genes into the non-transgenic or single Bt (Cry1Ac) gene cotton background where all three genes will inherit together. The expression of all three proteins was confirmed through immunostrips and was quantified through enzyme-linked immunosorbent assay (ELISA). The spatio-temporal expression of Bt protein in different parts of triple gene NIBGE cotton plants was determined. Maximum expression was found in leaves followed by seeds and boll rinds. Insect bioassays with cotton bollworms (Helicoverpa armigera), armyworms (Spodoptera litura), and pink bollworms (Pectinophora gossypiella) showed more than 90% mortality. The best performing line (NIBGE-E2) on the basis of spatiotemporal expression, glyphosate assays, and insect mortality data, was used for event characterization by using the genome sequencing approach. The event was successfully characterized and named NIBGE 20-01. A diagnostics test based on event-specific PCR was developed and its ability to distinguish NIBGE 20-01 event from other commercial transgenic cotton events was confirmed. To confirm stable expression of all three proteins in the field conditions, homozygous transgenic lines were grown in the field and the expression was confirmed through immunostrip assays. It was found that all three genes are expressed under field conditions. To show that all three genes are inherited together upon crossing with local elite cotton lines, the F₁ generation was grown under glasshouse and field conditions. The expression of all three genes was confirmed under field conditions. Our results showed that transgenic cotton with three genes cloned in the same T-DNA can express all genes and can be conveniently transferred into elite cotton lines through a single cross.

Transcriptomic Analysis of Mature Transgenic Poplar Expressing the Transcription Factor JERF36 Gene in Two Different Environments

During the last several decades, a number of transgenic or genetically modified tree varieties with enhanced characteristics and new traits have been produced. These trees have become associated with generally unsubstantiated concerns over health and environmental safety. We conducted transcriptome sequencing of transgenic Populus alba × P. berolinensis expressing the transcription factor JERF36 gene (ABJ01) and the non-transgenic progenitor line (9#) to compare the transcriptional changes in the apical buds. We found that 0.77% and 1.31% of the total expressed genes were significant differentially expressed in ABJ01 at the Daqing and Qiqihar sites, respectively. Among them, 30%–50% of the DEGs contained cis-elements recognized by JERF36. Approximately 5% of the total number of expressed genes showed significant differential expression between Daqing and Qiqihar in both ABJ01 and 9#. 10 DEGs resulting from foreign gene introduction, 394 DEGs that resulted solely from the environmental differences, and 47 DEGs that resulted from the combination of foreign gene introduction and the environment were identified. The number of DEGs resulting from environmental factors was significantly greater than that resulting from foreign gene introduction, and the combined effect of the environmental effects with foreign gene introduction was significantly greater than resulting from the introduction of JERF36 alone. GO and KEGG annotation showed that the DEGs mainly participate in the photosynthesis, oxidative phosphorylation, plant hormone signaling, ribosome, endocytosis, and plant-pathogen interaction pathways, which play important roles in the responses to biotic and abiotic stresses ins plant. To enhance its adaptability to salt-alkali stress, the transgenic poplar line may regulate the expression of genes that participate in the photosynthesis, oxidative phosphorylation, MAPK, and plant hormone signaling pathways. The crosstalk between biotic and abiotic stress responses by plant hormones may improve the ability of both transgenic and non-transgenic poplars to defend against pathogens. The results of our study provide a basis for further studies on the molecular mechanisms behind improved stress resistance and the unexpected effects of transgenic gene expression in poplars, which will be significant for improving the biosafety evaluation of transgenic trees and accelerating the breeding of new varieties of forest trees resistant to environmental stresses.

GmoDetector: An accurate and efficient GMO identification approach and its applications

The rapid increase of genetically modified organisms (GMOs) entering the food and feed markets, and the contamination of donor (micro)organisms of transgenic elements make it more challenging for the existing GMO detection. In this study, we developed a high-throughput and contamination-removal GMO detection approach named as GmoDetector. GmoDetector targeted 64 common transgenic elements and 76 GMO-specific events collected from 251 singular GM events, and combined with next generation sequencing (NGS) and target enrichment technology to detect various GMOs. As a result, GmoDetector was able to exclude the donor (micro)organism contamination, and detect the authorized and unauthorized GMOs (UGMOs) in any forms of food or feed, such as processed or unprocessed. The sensitivity of GmoDetector is as low as 0.1% (GMO content), which has met the GMO labeling threshold for all countries. Therefore, GmoDetector is a robust tool for accurate and efficient detection of the authorized and UGMOs.

Influence of transgenesis on genome variability in cucumber lines with a thaumatin II gene

The development of new plant varieties by genetic modification aims at improving their features or introducing new qualities. However, concerns about the unintended effects of transgenes and negative environmental impact of genetically modified plants are an obstacle for the use of these plants in crops. To analyze the impact of transgenesis on plant genomes, we analyze three cucumber transgenic lines with an introduced thaumatin II gene. After genomes sequencing, we analyzed the transgene insertion site and performed variant prediction. As a result, we obtained similar number of variants for all analyzed lines (average of 4307 polymorphisms), with high abundance in one region of chromosome 4. According to SnpEff analysis, the presence of genomic variants generally does not influence the genome functionality, as less than 2% of polymorphisms have high impact. Moreover, analysis indicates that these changes were more likely induced by in vitro culture than by the transgenesis itself. The insertion site analysis shows that the region of transgene integration could cause changes in gene expression, by gene disruption or loss of promoter region continuity.

SUPPLEMENTARY INFORMATION

The online version contains supplementary material available at 10.1007/s12298-021-00990-8.

Molecular characterization for food safety assessment of a genetically modified late blight resistant potato: an unusual case

Standard food safety assessments of genetically modified crops require a thorough molecular characterization of the novel DNA as inserted into the plant that is intended for commercialization, as well as a comparison of agronomic and nutritional characteristics of the genetically modified to the non-modified counterpart. These characterization data are used to identify any unintended changes in the inserted DNA or in the modified plant that would require assessment for safety in addition to the assessment of the intended modification. An unusual case of an unintended effect discovered from the molecular characterization of a genetically modified late blight resistant potato developed for growing in Bangladesh and Indonesia is presented here. Not only was a significant portion of the plasmid vector backbone DNA inserted into the plant along with the intended insertion of an R-gene for late blight resistance, but the inserted DNA was split into two separate fragments and inserted into two separate chromosomes. One fragment carries the R-gene and the other fragment carries the NPTII selectable marker gene and the plasmid backbone DNA. The implications of this for the food safety assessment of this late blight resistant potato are considered.

Development of event-specific detection method for identification of insect resistant NIBGE-1601 cotton harboring double gene Cry1Ac-Cry2Ab construct

Bt cotton expressing Cry1Ac is being cultivated in Pakistan. It has been observed that pink bollworm may have developed resistance against single Bt gene (Cry1Ac). For durable resistance, insect resistant NIBGE-1601 cotton harboring double gene Cry1Ac-Cry2Ab construct was developed. There was a need to characterize NIBGE-1601 event for intellectual property rights protection. The Presence of NIBGE Cry1Ac and NIBGE Cry2Ab genes was checked in NIBGE-1601 cotton plants through PCR, while there was no amplification using primers specific for Monsanto events (MON531, MON15985, MON1445). Using genome walking technology, NIBGE-601 event has been characterized. Event-specific primers of NIBGE-1601 were designed and evaluated to differentiate it from other cotton events mentioned above. NIBGE-1601 event detection primers are highly specific, therefore, can detect NIBGE 1601 event at different conditions using single or multiplex PCR. In the qualitative PCR, using NIBGE-1601 event specific primers, 0.05 ng was the limit of detection for NIBGE-1601double gene cotton genomic DNA. Thus event characterization and development of event-specific diagnostics will help in breeding new cotton varieties resistant to cotton bollworms.

Rapid and Detailed Characterization of Transgene Insertion Sites in Genetically Modified Plants via Nanopore Sequencing

Molecular characterization of genetically modified plants can provide crucial information for the development of detection and identification methods, to comply with traceability, and labeling requirements prior to commercialization. Detailed description of the genetic modification was previously a challenging step in the safety assessment, since it required the use of laborious and time-consuming techniques. In this study an accurate, simple, and fast method was developed for molecular characterization of genetically modified (GM) plants, following a user-friendly workflow for researchers with limited bioinformatic capabilities. Three GM events from a diverse array of crop species-perennial ryegrass, white clover, and canola-were used to test the approach that exploits long-read sequencing by the MinION device, from Oxford Nanopore Technologies. The method delivered a higher degree of resolution of the transgenic events within the host genome than has previously been possible with the standard Illumina short-range sequencing strategies. The flanking sequences, copy number, and presence of backbone sequences, and overall transgene insertion structure were determined for each of the plant genomes, with the additional identification of moderate-sized secondary insertions that would have previously been missed. The proposed workflow takes only about 1 week from DNA extraction to analyzed result, and the method will complement the existing approaches for molecular characterization of GM plants, since it makes the process faster, simpler, and more cost-effective.

Evaluation on reprogramed biological processes in transgenic maize varieties using transcriptomics and metabolomics

Genetic engineering (GM) has great potential to improve maize productivity, but rises some concerns on unintended effects, and equivalent as their comparators. There are some limitations through targeted analysis to detect the UE in genetically modified organisms in many previous studies. We here reported a case-study on the effects of introducing herbicides and insect resistance (HIR) gene cassette on molecular profiling (transcripts and metabolites) in a popular maize variety Zhengdan958 (ZD958) in China. We found that introducing HIR gene cassette bring a limited numbers of differential abundant genes (DAGs) or differential abundant metabolites (DAMs) between transgenic events and non-transgenic control. In contrast, averaged 10 times more DAGs and DAMs were observed when performed comparison under different growing environments in three different ecological regions of China than the numbers induced by gene effects. Major biological pathways relating to stress response or signaling transduction could explain somehow the effects of growing environments. We further compared two transgenic events mediated ZD958 (GM-ZD958) with either transgenic parent GM-Z58, and other genetic background nonGM-Z58, nonGM-ZD958, and Chang7-2. We found that the numbers of DAGs and DAMs between GM-ZD958 and its one parent maize variety, Z58 or GM-Z58 is equivalent, but not Chang7-2. These findings suggest that greater effects due to different genetic background on altered molecular profiling than gene modification itself. This study provides a case evidence indicating marginal effects of gene pleiotropic effects, and environmental effects should be emphasized.

Elucidation of genomic organizations of transgenic soybean plants through de novo genome assembly with short paired-end reads

Elucidation of the genomic organizations of transgene insertion sites is essential for the genetic studies of transgenic plants. Herein, we establish an analysis pipeline that identifies the transgene insertion sites as well as the presence of vector backbones, through de novo genome assembly with high-throughput sequencing data in two transgenic soybean lines, AtYUCCA6-#5 and 35S-UGT72E3/2-#7. Sequencing data of approximately 28× and 29× genome coverages for each line generated by high-throughput sequencing were de novo assembled. The databases generated from the de novo assembled sequences were used to search contigs that contained putative insertion sites and their flanking sequences (integration sites) of transgene fragments using transgenic vector sequences as queries. The predicted integration site sequences, which are located at three annotated genes that might regulate plant development or confer disease resistance, were then confirmed by local alignment against the soybean reference genome and PCR amplification. As results, we revealed the precise transgene-flanking sequences and sequence rearrangements at insertion sites in both the transgenic lines, as well as the aberrant insertion of a transgene fragment. Consequently, relative to experimental or enrichment technologies, our approach is straightforward and time-effective, providing an alternative method for the identification of insertion sites in transgenic plants.

Defining the 'HoneySweet' insertion event utilizing NextGen sequencing and a de novo genome assembly of plum (Prunus domestica)

'HoneySweet' plum (Prunus domestica) is resistant to Plum pox potyvirus, through an RNAi-triggered mechanism. Determining the precise nature of the transgene insertion event has been complicated due to the hexaploid genome of plum. DNA blots previously indicated an unintended hairpin arrangement of the Plum pox potyvirus coat protein gene as well as a multicopy insertion event. To confirm the transgene arrangement of the insertion event, 'HoneySweet' DNA was subjected to whole genome sequencing using Illumina short-read technology. Results indicated two different insertion events, one containing seven partial copies flanked by putative plum DNA sequence and a second with the predicted inverted repeat of the coat protein gene driven by a double 35S promoter on each side, flanked by plum DNA. To determine the locations of the two transgene insertions, a phased plum genome assembly was developed from the commercial plum 'Improved French'. A subset of the scaffolds (2447) that were >10 kb in length and representing, >95% of the genome were annotated and used for alignment against the 'HoneySweet' transgene reads. Four of eight matching scaffolds spanned both insertion sites ranging from 157,704 to 654,883 bp apart, however we were unable to identify which scaffold(s) represented the actual location of the insertion sites due to potential sequence differences between the two plum cultivars. Regardless, there was no evidence of any gene(s) being interrupted as a result of the insertions. Furthermore, RNA-seq data verified that the insertions created no new transcriptional units and no dramatic expression changes of neighboring genes.

Detection of 30 bp DNA fragments with a sensitive modified Southern blot analysis

To evaluate crops generated by new breeding techniques, it is important to confirm the removal of recombinant DNAs (rDNAs) derived from foreign genes including unintentionally introduced short rDNA(s). We attempted to develop a sensitive detection method for such short rDNAs using Southern blot analysis and performed a model study targeting single-copy endogenous genes in plants. To increase the detection sensitivity, the general protocol for Southern blot analysis was modified. In the model study, we used endogenous-gene-targeting probes in which complementary sequences were serially replaced by dummy sequences, and detected complementary sequences as well as 30 bp. We further evaluated the sensitivity using short rDNAs derived from GM sequences as pseudoinsertions, and the results demonstrated that rDNA-insertions as small as 30 bp could be detected. The results suggested that unintentionally introduced rDNA-insertions were 30 bp or more in length could be detected by the Southern blot analysis.

Sequencing of two transgenic early-flowering poplar lines confirmed vector-free single-locus T-DNA integration

Next-generation sequencing (NGS) approaches are attractive alternatives to the PCR-based characterisation of genetically modified plants for safety assessment and labelling since NGS is highly sensitive to the detection of T-DNA inserts as well as vector backbone sequences in transgenic plants. In this study, two independent transgenic male Populus tremula lines, T193-2 and T195-1, both carrying the FLOWERING LOCUS T gene from Arabidopsis thaliana under control of a heat-inducible promoter (pHSP::AtFT) and the non-transgenic control clone W52, were further characterised by NGS and third-generation sequencing. The results support previous findings that the T-DNA was hemizygously inserted in one genomic locus of each line. However, the T-DNA insertions consist of conglomerations of one or two T-DNA copies together with a small T-DNA fragment without AtFT parts. Based on NGS data, no additional T-DNA splinters or vector backbone sequences could be identified in the genome of the two transgenic lines. Seedlings derived from crosses between the pHSP::AtFT transgenic male parents and female wild type plants are therefore expected to be T-DNA splinter or vector backbone free. Thus, PCR analyses amplifying a partial T-DNA fragment with AtFT-specific primers are sufficient to determine whether the seedlings are transgenic or not. An analysis of 72 second generation-seedlings clearly showed that about 50% of them still reveal the presence of the T-DNA, confirming data already published. To prove if unanticipated genomic changes were induced by T-DNA integration, extended future studies using long-range sequencing technologies are required once a suitable chromosome-level P. tremula reference genome sequence is available.

Identification of genomic insertion and flanking sequences of the transgenic drought-tolerant maize line "SbSNAC1-382" using the single-molecule real-time (SMRT) sequencing method

Safety assessment of genetically modified (GM) crops is crucial at the product-development phase before GM crops are placed on the market. Determining characteristics of sequences flanking exogenous insertion sequences is essential for the safety assessment and marketing of transgenic crops. In this study, we used genome walking and whole-genome sequencing (WGS) to identify the flanking sequence characteristics of the SbSNAC1 transgenic drought-tolerant maize line "SbSNAC1-382", but both of the two methods failed. Then, we constructed a genomic fosmid library of the transgenic maize line, which contained 4.18×105 clones with an average insertion fragment of 35 kb, covering 5.85 times the maize genome. Subsequently, three positive clones were screened by pairs of specific primers, and one of the three positive clones was sequenced by using single-molecule real-time (SMRT) sequencing technology. More than 1.95 Gb sequence data (~105× coverage) for the sequenced clone were generated. The junction reads mapped to the boundaries of T-DNA, and the flanking sequences in the transgenic line were identified by comparing all sequencing reads with the maize reference genome and the sequence of the transgenic vector. Furthermore, the putative insertion loci and flanking sequences were confirmed by PCR amplification and Sanger sequencing. The results indicated that two copies of the exogenous T-DNA fragments were inserted at the same genomic site, and the exogenous T-DNA fragments were integrated at the position of Chromosome 5 from 177155650 to 177155696 in the transgenic line 382. In this study, we demonstrated the successful application of the SMRT technology for the characterization of genomic insertion and flanking sequences.

Foreign DNA detection by high-throughput sequencing to regulate genome-edited agricultural products

Although the advent of several new breeding techniques (NBTs) is revolutionizing agricultural production processes, technical information necessary for their regulation is yet to be provided. Here, we show that high-throughput DNA sequencing is effective for the detection of unintended remaining foreign DNA segments in genome-edited rice. A simple k-mer detection method is presented and validated through a series of computer simulations and real data analyses. The data show that a short foreign DNA segment of 20 nucleotides can be detected and the probability that the segment is overlooked is 10-3 or less if the average sequencing depth is 30 or more, while the number of false hits is less than 1 on average. This method was applied to real sequencing data, and the presence and absence of an external DNA segment were successfully proven. Additionally, our in-depth analyses also identified some weaknesses in current DNA sequencing technologies. Hence, for a rigorous safety assessment, the combination of k-mer detection and another method, such as Southern blot assay, is recommended. The results presented in this study will lay the foundation for the regulation of NBT products, where foreign DNA is utilized during their generation.

Identification of T-DNA Insertion Site and Flanking Sequence of a Genetically Modified Maize Event IE09S034 Using Next-Generation Sequencing Technology

Molecular characteristics including information of insertion site, flanking sequence, and copy numbers are the base for the safety assessment and subsequent monitoring of genetically modified organisms (GMOs), which has to be revealed thoroughly in a case-by-case manner. Although both polymerase chain reaction (PCR)-based and next-generation sequencing (NGS)-based approaches are proven to be effective in the molecular characterization of most of GM events, they often fail to work with GM maize events, mainly due to the genome complexity. In this study, by using NGS, we successfully identified the 3' end T-DNA insertion site and flanking sequence of a GM maize event IE09S034, which were confirmed by PCR amplification and Sanger sequencing. Notably, insertions of unintended exogenous elements were revealed in this event although the single copy of target exogenous genes was also confirmed by digital PCR. The output of this study provides novel and important genetic evidence for the safety assessment and monitoring of GM maize event IE09S034.

Safety Assessment of Genetically Modified Feed: Is There Any Difference From Food?

Food security is one of major concerns for the growing global population. Modern agricultural biotechnologies, such as genetic modification, are a possible solution through enabling an increase of production, more efficient use of natural resources, and reduced environmental impacts. However, new crop varieties with altered genetic materials may be subjected to safety assessments to fulfil the regulatory requirements, prior to marketing. The aim of the assessment is to evaluate the impact of products from the new crop variety on human, animal, and the environmental health. Although, many studies on the risk assessment of genetically modified (GM) food have been published, little consideration to GM feedstuff has been given, despite that between 70 to 90% of all GM crops and their biomass are used as animal feed. In addition, in some GM plants such as forages that are only used for animal feeds, the assessment of the genetic modification may be of relevance only to livestock feeding. In this article, the regulatory framework of GM crops intended for animal feed is reviewed using the available information on GM food as the baseline. Although, the majority of techniques used for the safety assessment of GM food can be used in GM feed, many plant parts used for livestock feeding are inedible to humans. Therefore, the concentration of novel proteins in different plant tissues and level of exposure to GM feedstuff in the diet of target animals should be considered. A further development of specific methodologies for the assessment of GM crops intended for animal consumption is required, in order to provide a more accurate and standardized assessment to the GM feed safety.

A high-throughput and ultrasensitive identification methodology for unauthorized GMP component based on suspension array and logical calculator

To solve the problem of the unauthorized GMP components within import and export goods, the LI-US (Logic Identification of unauthorized GMP content by Universal-primer Suspension-array) system, which takes advantage of suspension array and logic calculator, was developed in the present study. Seventeen signal input channels have been optimized and validated in our research to ensure the multiplex practicality of the LI-US system. Three LI-US logic gates, including a YES gate, an OR gate and an AND gate, were designed as different detection strategies for GMP identification. The feasibility and specificity of the LI-US system were validated in the present study. Combining the optimization and evaluation of the signal input procedure, the sensitivity of this LI-US system reached 0.05% of the GMP mass concentration. The practicability evaluation of LI-US demonstrated its application within different substrates and varieties. In conclusion, the LI-US system was developed with extremely high specificity, sensitivity and practicability among different substrates and varieties, which could meet the demands of unauthorized GMP contents for both import and export goods.

Unintended effects of transgenic rice revealed by transcriptome and metabolism

Genetically modified (GM) organisms have been developed for decades. However, unintended effects are the main concerns of safety assessment that needs to be carefully investigated. Here, eight varieties of GM rice that were developed in China were selected to assess the unintended effects through transcriptome and metabolism. There are 2892-8758 differentially expressed genes (DEGs) and 7-50 metabolites at significant level between GM varieties and their isogenic counterparts, which were far fewer than that between traditional rice varieties. The function enrichment analysis showed altered transcription in stress-related pathway and starch and sucrose metabolism. DEGs shared among eight GM samples constitute less than 1% of the genes in the genome, and none of them is reported more than four times. The insertion effect on the nearby gene expression and the associated metabolism is only restricted to 50 genes. All the results provide a comprehensive analysis of unintended effects and indication of difference in Chinese transgenic rice based on their backgrounds, transformation, and insertion elements.

Detection and Identification of Genome Editing in Plants: Challenges and Opportunities

Conventional genetic engineering techniques generate modifications in the genome via stable integration of DNA elements which do not occur naturally in this combination. Therefore, the resulting organisms and (most) products thereof can unambiguously be identified with event-specific PCR-based methods targeting the insertion site. New breeding techniques such as genome editing diversify the toolbox to generate genetic variability in plants. Several of these techniques can introduce single nucleotide changes without integrating foreign DNA and thereby generate organisms with intended phenotypes. Consequently, such organisms and products thereof might be indistinguishable from naturally occurring or conventionally bred counterparts with established analytical tools. The modifications can entirely resemble random mutations regardless of being spontaneous or induced chemically or via irradiation. Therefore, if an identification of these organisms or products thereof is demanded, a new challenge will arise for (official) seed, food, and feed testing laboratories and enforcement institutions. For detailed consideration, we distinguish between the detection of sequence alterations - regardless of their origin - the identification of the process that generated a specific modification and the identification of a genotype, i.e., an organism produced by genome editing carrying a specific genetic alteration in a known background. This article briefly reviews the existing and upcoming detection and identification strategies (including the use of bioinformatics and statistical approaches) in particular for plants developed with genome editing techniques.

Technical Note on the quality of DNA sequencing for the molecular characterisation of genetically modified plants

As part of the risk assessment (RA) requirements for genetically modified (GM) plants, according to Regulation (EU) No 503/2013 and the EFSA guidance on the RA of food and feed from GM plants (EFSA GMO Panel, 2011), applicants need to perform a molecular characterisation of the DNA sequences inserted in the GM plant genome. The European Commission has mandated EFSA to develop a technical note to the applicants on, and checking of, the quality of the methodology, analysis and reporting covering complete sequencing of the insert and flanking regions, insertion site analysis of the GM event, and generational stability and integrity. This Technical Note puts together requirements and recommendations for when DNA sequencing is part of the molecular characterisation of GM plants, in particular for the characterisation of the inserted genetic material at each insertion site and flanking regions, the determination of the copy number of all detectable inserts, and the analysis of the genetic stability of the inserts, when addressed by Sanger sequencing or NGS. This document reflects the current knowledge in scientific-technical methods for generating and verifying, in a standardised manner, DNA sequencing data in the context of RA of GM plants. From 1 October 2018, this Technical Note will replace the JRC guideline of 2016 (updated April 2017) related to the verification and quality assessment of the sequencing of the insert(s) and flanking regions. It does not take into consideration the verification and validation of the detection method which remains under the remit of the JRC.

The importance of genotype identity, genetic heterogeneity, and bioinformatic handling for properly assessing genomic variation in transgenic plants

BACKGROUND

The advent of -omics technologies has enabled the resolution of fine molecular differences among individuals within a species. DNA sequence variations, such as single nucleotide polymorphisms or small deletions, can be tabulated for many kinds of genotype comparisons. However, experimental designs and analytical approaches are replete with ways to overestimate the level of variation present within a given sample. Analytical pipelines that do not apply proper thresholds nor assess reproducibility among samples are susceptible to calling false-positive variants. Furthermore, issues with sample genotype identity or failing to account for heterogeneity in reference genotypes may lead to misinterpretations of standing variants as polymorphisms derived de novo.

RESULTS

A recent publication that featured the analysis of RNA-sequencing data in three transgenic soybean event series appeared to overestimate the number of sequence variants identified in plants that were exposed to a tissue culture based transformation process. We reanalyzed these data with a stringent set of criteria and demonstrate three different factors that lead to variant overestimation, including issues related to the genetic identity of the background genotype, unaccounted genetic heterogeneity in the reference genome, and insufficient bioinformatics filtering.

CONCLUSIONS

This study serves as a cautionary tale to users of genomic and transcriptomic data that wish to assess the molecular variation attributable to tissue culture and transformation processes. Moreover, accounting for the factors that lead to sequence variant overestimation is equally applicable to samples derived from other germplasm sources, including chemical or irradiation mutagenesis and genome engineering (e.g., CRISPR) processes.

Ultrasensitive Single Fluorescence-Labeled Probe-Mediated Single Universal Primer-Multiplex-Droplet Digital Polymerase Chain Reaction for High-Throughput Genetically Modified Organism Screening

As genetically modified (GM) technology develops and genetically modified organisms (GMOs) become more available, GMOs face increasing regulations and pressure to adhere to strict labeling guidelines. A singleplex detection method cannot perform the high-throughput analysis necessary for optimal GMO detection. Combining the advantages of multiplex detection and droplet digital polymerase chain reaction (ddPCR), a single universal primer-multiplex-ddPCR (SUP-M-ddPCR) strategy was proposed for accurate broad-spectrum screening and quantification. The SUP increases efficiency of the primers in PCR and plays an important role in establishing a high-throughput, multiplex detection method. Emerging ddPCR technology has been used for accurate quantification of nucleic acid molecules without a standard curve. Using maize as a reference point, four heterologous sequences ( 35S, NOS, NPTII, and PAT) were selected to evaluate the feasibility and applicability of this strategy. Surprisingly, these four genes cover more than 93% of the transgenic maize lines and serve as preliminary screening sequences. All screening probes were labeled with FAM fluorescence, which allows the signals from the samples with GMO content and those without to be easily differentiated. This fiveplex screening method is a new development in GMO screening. Utilizing an optimal amplification assay, the specificity, limit of detection (LOD), and limit of quantitation (LOQ) were validated. The LOD and LOQ of this GMO screening method were 0.1% and 0.01%, respectively, with a relative standard deviation (RSD) < 25%. This method could serve as an important tool for the detection of GM maize from different processed, commercially available products. Further, this screening method could be applied to other fields that require reliable and sensitive detection of DNA targets.

A bioinformatics approach for identifying transgene insertion sites using whole genome sequencing data

BACKGROUND

Genetically modified crops (GM crops) have been developed to improve the agricultural traits of modern crop cultivars. Safety assessments of GM crops are of paramount importance in research at developmental stages and before releasing transgenic plants into the marketplace. Sequencing technology is developing rapidly, with higher output and labor efficiencies, and will eventually replace existing methods for the molecular characterization of genetically modified organisms.

METHODS

To detect the transgenic insertion locations in the three GM rice gnomes, Illumina sequencing reads are mapped and classified to the rice genome and plasmid sequence. The both mapped reads are classified to characterize the junction site between plant and transgene sequence by sequence alignment.

RESULTS

Herein, we present a next generation sequencing (NGS)-based molecular characterization method, using transgenic rice plants SNU-Bt9-5, SNU-Bt9-30, and SNU-Bt9-109. Specifically, using bioinformatics tools, we detected the precise insertion locations and copy numbers of transfer DNA, genetic rearrangements, and the absence of backbone sequences, which were equivalent to results obtained from Southern blot analyses.

CONCLUSION

NGS methods have been suggested as an effective means of characterizing and detecting transgenic insertion locations in genomes. Our results demonstrate the use of a combination of NGS technology and bioinformatics approaches that offers cost- and time-effective methods for assessing the safety of transgenic plants.

Safety, Security, and Policy Considerations for Plant Genome Editing

Genome editing with engineered nucleases (GEEN) is increasingly used as a tool for gene discovery and trait development in crops through generation of targeted changes in endogenous genes. The development of the CRISPR-Cas9 system (clustered regularly interspaced short palindromic repeats with associated Cas9 protein), in particular, has enabled widespread use of genome editing. Research to date has not comprehensively addressed genome-editing specificity and off-target mismatches that may result in unintended changes within plant genomes or the potential for gene drive initiation. Governance and regulatory considerations for bioengineered crops derived from using GEEN will require greater clarity as to target specificity, the potential for mismatched edits, unanticipated downstream effects of off-target mutations, and assurance that genome reagents do not occur in finished products. Since governance and regulatory decision making involves robust standards of evidence extending from the laboratory to the postcommercial marketplace, developers of genome-edited crops must anticipate significant engagement and investment to address questions of regulators and civil society.

Molecular characterization of an unauthorized genetically modified Bacillus subtilis production strain identified in a vitamin B2 feed additive

•

Genetically modified Bacillus subtilis identified in a vitamin B₂ product.

•

Whole genome sequencing runs are performed for characterization of the isolated strain.

•

Complex modifications of the genome are identified.

•

Four putative recombinant plasmids are characterized.

•

Real-time PCR methods are developed and available for testing vitamin B₂ products.

Many food and feed additives result from fermentation of genetically modified (GM) microorganisms. For vitamin B2 (riboflavin), GM Bacillus subtilis production strains have been developed and are often used. The presence of neither the GM strain nor its recombinant DNA is allowed for fermentation products placed on the EU market as food or feed additive. A vitamin B₂ product (80% feed grade) imported from China was analysed. Viable B. subtilis cells were identified and DNAs of two bacterial isolates (LHL and LGL) were subjected to three whole genome sequencing (WGS) runs with different devices (MiSeq, 454 or HiSeq system). WGS data revealed the integration of a chloramphenicol resistance gene, the deletion of the endogenous riboflavin (rib) operon and presence of four putative plasmids harbouring rib operons. Event- and construct-specific real-time PCR methods for detection of the GM strain and its putative plasmids in food and feed products have been developed.

Re-sequencing transgenic plants revealed rearrangements at T-DNA inserts, and integration of a short T-DNA fragment, but no increase of small mutations elsewhere

KEY MESSAGE

Transformation resulted in deletions and translocations at T-DNA inserts, but not in genome-wide small mutations. A tiny T-DNA splinter was detected that probably would remain undetected by conventional techniques. We investigated to which extent Agrobacterium tumefaciens-mediated transformation is mutagenic, on top of inserting T-DNA. To prevent mutations due to in vitro propagation, we applied floral dip transformation of Arabidopsis thaliana. We re-sequenced the genomes of five primary transformants, and compared these to genomic sequences derived from a pool of four wild-type plants. By genome-wide comparisons, we identified ten small mutations in the genomes of the five transgenic plants, not correlated to the positions or number of T-DNA inserts. This mutation frequency is within the range of spontaneous mutations occurring during seed propagation in A. thaliana, as determined earlier. In addition, we detected small as well as large deletions specifically at the T-DNA insert sites. Furthermore, we detected partial T-DNA inserts, one of these a tiny 50-bp fragment originating from a central part of the T-DNA construct used, inserted into the plant genome without flanking other T-DNA. Because of its small size, we named this fragment a T-DNA splinter. As far as we know this is the first report of such a small T-DNA fragment insert in absence of any T-DNA border sequence. Finally, we found evidence for translocations from other chromosomes, flanking T-DNA inserts. In this study, we showed that next-generation sequencing (NGS) is a highly sensitive approach to detect T-DNA inserts in transgenic plants.

Achieving Plant CRISPR Targeting that Limits Off-Target Effects

The CRISPR-Cas9 system (clustered regularly interspaced short palindromic repeats with associated Cas9 protein) has been used to generate targeted changes for direct modification of endogenous genes in an increasing number of plant species; but development of plant genome editing has not yet fully considered potential off-target mismatches that may lead to unintended changes within the genome. Assessing the specificity of CRISPR-Cas9 for increasing editing efficiency as well as the potential for unanticipated downstream effects from off-target mutations is an important regulatory consideration for agricultural applications. Increasing genome-editing specificity entails developing improved design methods that better predict the prevalence of off-target mutations as a function of genome composition and design of the engineered ribonucleoprotein (RNP). Early results from CRISPR-Cas9 genome editing in plant systems indicate that the incidence of off-target mutation frequencies is quite low; however, by analyzing CRISPR-edited plant lines and improving both computational tools and reagent design, it may be possible to further decrease unanticipated effects at potential mismatch sites within the genome. This will provide assurance that CRISPR-Cas9 reagents can be designed and targeted with a high degree of specificity. Improved and experimentally validated design tools for discriminating target and potential off-target positions that incorporate consideration of the designed nuclease fidelity and selectivity will help to increase confidence for regulatory decision making for genome-edited plants.

Identification of Genomic Insertion and Flanking Sequence of G2-EPSPS and GAT Transgenes in Soybean Using Whole Genome Sequencing Method

Molecular characterization of sequence flanking exogenous fragment insertion is essential for safety assessment and labeling of genetically modified organism (GMO). In this study, the T-DNA insertion sites and flanking sequences were identified in two newly developed transgenic glyphosate-tolerant soybeans GE-J16 and ZH10-6 based on whole genome sequencing (WGS) method. More than 22.4 Gb sequence data (∼21 × coverage) for each line was generated on Illumina HiSeq 2500 platform. The junction reads mapped to boundaries of T-DNA and flanking sequences in these two events were identified by comparing all sequencing reads with soybean reference genome and sequence of transgenic vector. The putative insertion loci and flanking sequences were further confirmed by PCR amplification, Sanger sequencing, and co-segregation analysis. All these analyses supported that exogenous T-DNA fragments were integrated in positions of Chr19: 50543767-50543792 and Chr17: 7980527-7980541 in these two transgenic lines. Identification of genomic insertion sites of G2-EPSPS and GAT transgenes will facilitate the utilization of their glyphosate-tolerant traits in soybean breeding program. These results also demonstrated that WGS was a cost-effective and rapid method for identifying sites of T-DNA insertions and flanking sequences in soybean.

Synthetic strategies for plant signalling studies: molecular toolbox and orthogonal platforms

Plants deploy a wide array of signalling networks integrating environmental cues with growth, defence and developmental responses. The high level of complexity, redundancy and connection between several pathways hampers a comprehensive understanding of involved functional and regulatory mechanisms. The implementation of synthetic biology approaches is revolutionizing experimental biology in prokaryotes, yeasts and animal systems and can likewise contribute to a new era in plant biology. This review gives an overview on synthetic biology approaches for the development and implementation of synthetic molecular tools and techniques to interrogate, understand and control signalling events in plants, ranging from strategies for the targeted manipulation of plant genomes up to the spatiotemporally resolved control of gene expression using optogenetic approaches. We also describe strategies based on the partial reconstruction of signalling pathways in orthogonal platforms, like yeast, animal and in vitro systems. This allows a targeted analysis of individual signalling hubs devoid of interconnectivity with endogenous interacting components. Implementation of the interdisciplinary synthetic biology tools and strategies is not exempt of challenges and hardships but simultaneously most rewarding in terms of the advances in basic and applied research. As witnessed in other areas, these original theoretical-experimental avenues will lead to a breakthrough in the ability to study and comprehend plant signalling networks.

Application of whole genome shotgun sequencing for detection and characterization of genetically modified organisms and derived products

The emergence of high-throughput, massive or next-generation sequencing technologies has created a completely new foundation for molecular analyses. Various selective enrichment processes are commonly applied to facilitate detection of predefined (known) targets. Such approaches, however, inevitably introduce a bias and are prone to miss unknown targets. Here we review the application of high-throughput sequencing technologies and the preparation of fit-for-purpose whole genome shotgun sequencing libraries for the detection and characterization of genetically modified and derived products. The potential impact of these new sequencing technologies for the characterization, breeding selection, risk assessment, and traceability of genetically modified organisms and genetically modified products is yet to be fully acknowledged. The published literature is reviewed, and the prospects for future developments and use of the new sequencing technologies for these purposes are discussed.

Consumer acceptance of food crops developed by genome editing

One of the major problems regarding consumer acceptance of genetically modified organisms (GMOs) is the possibility that their transgenes could have adverse effects on the environment and/or human health. Genome editing, represented by the CRISPR/Cas9 system, can efficiently achieve transgene-free gene modifications and is anticipated to generate a wide spectrum of plants. However, the public attitude against GMOs suggests that people will initially be unlikely to accept these plants. We herein explored the bottlenecks of consumer acceptance of transgene-free food crops developed by genome editing and made some recommendations. People should not pursue a zero-risk bias regarding such crops. Developers are encouraged to produce cultivars with a trait that would satisfy consumer needs. Moreover, they should carefully investigate off-target mutations in resultant plants and initially refrain from agricultural use of multiplex genome editing for better risk-benefit communication. The government must consider their regulatory status and establish appropriate regulations if necessary. The government also should foster communication between the public and developers. If people are informed of the benefits of genome editing-mediated plant breeding and trust in the relevant regulations, and if careful risk-benefit communication and sincere considerations for the right to know approach are guaranteed, then such transgene-free crops could gradually be integrated into society.

Current trends in Bt crops and their fate on associated microbial community dynamics: a review

Cry protein expressing insect-resistant trait is mostly deployed to control major devastating pests and minimize reliance on the conventional pesticides. However, the ethical and environmental issues are the major constraints in their acceptance, and consequently, the cultivation of genetically modified (GM) crops has invited intense debate. Since root exudates of Bacillus thuringiensis (Bt) crops harbor the insecticidal protein, there is a growing concern about the release and accumulation of soil-adsorbed Cry proteins and their impact on non-target microorganisms and soil microbial processes. This review pertains to reports from the laboratory studies and field trials to assess the Bt toxin proteins in soil microbes and the processes determining the soil quality in conjunction with the existing hypothesis and molecular approaches to elucidate the risk posed by the GM crops. Ecological perturbations hinder the risk aspect of soil microbiota in response to GM crops. Therefore, extensive research based on in vivo and interpretation of results using high-throughput techniques such as NGS on risk assessment are imperative to evaluate the impact of Bt crops to resolve the controversy related to their commercialization. But more studies are needed on the risk associated with stacked traits. Such studies would strengthen our knowledge about the plant-microbe interactions.

Molecular Characterization of Transgenic Events Using Next Generation Sequencing Approach

Demand for the commercial use of genetically modified (GM) crops has been increasing in light of the projected growth of world population to nine billion by 2050. A prerequisite of paramount importance for regulatory submissions is the rigorous safety assessment of GM crops. One of the components of safety assessment is molecular characterization at DNA level which helps to determine the copy number, integrity and stability of a transgene; characterize the integration site within a host genome; and confirm the absence of vector DNA. Historically, molecular characterization has been carried out using Southern blot analysis coupled with Sanger sequencing. While this is a robust approach to characterize the transgenic crops, it is both time- and resource-consuming. The emergence of next-generation sequencing (NGS) technologies has provided highly sensitive and cost- and labor-effective alternative for molecular characterization compared to traditional Southern blot analysis. Herein, we have demonstrated the successful application of both whole genome sequencing and target capture sequencing approaches for the characterization of single and stacked transgenic events and compared the results and inferences with traditional method with respect to key criteria required for regulatory submissions.

Efficiency to Discovery Transgenic Loci in GM Rice Using Next Generation Sequencing Whole Genome Re-sequencing

Molecular characterization technology in genetically modified organisms, in addition to how transgenic biotechnologies are developed now require full transparency to assess the risk to living modified and non-modified organisms. Next generation sequencing (NGS) methodology is suggested as an effective means in genome characterization and detection of transgenic insertion locations. In the present study, we applied NGS to insert transgenic loci, specifically the epidermal growth factor (EGF) in genetically modified rice cells. A total of 29.3 Gb (~72× coverage) was sequenced with a 2 × 150 bp paired end method by Illumina HiSeq2500, which was consecutively mapped to the rice genome and T-vector sequence. The compatible pairs of reads were successfully mapped to 10 loci on the rice chromosome and vector sequences were validated to the insertion location by polymerase chain reaction (PCR) amplification. The EGF transgenic site was confirmed only on chromosome 4 by PCR. Results of this study demonstrated the success of NGS data to characterize the rice genome. Bioinformatics analyses must be developed in association with NGS data to identify highly accurate transgenic sites.