Detection and identification of transgenic events by next generation sequencing combined with enrichment technologies

An integrated strategy involving high-throughput sequencing to characterize an unknown GM wheat event in Canada

Glyphosate-resistant wheat plants were discovered in southern Alberta in 2017, representing an unauthorized GM release in Canada. The Canadian Food Inspection Agency undertook a series of experiments to characterize and identify this unknown GM wheat, as well as to develop and validate construct-specific and event-specific qPCR assays. Results of PCR-based assays and Sanger sequencing indicated the presence of CaMV 35S promoter (p35S), Rice Actin 1 intron (RactInt1), CP4-EPSPS gene and nopaline synthase terminator (tNOS) elements in the unknown GM wheat. Genome walking and bead capture strategies, combined with high-throughput sequencing, were used to identify the 5' and 3' wheat junctions and the subsequent mapping of the insert to chromosome 3B of the wheat genome. A probable transformation vector, pMON25497, was recognized, and further testing identified the unknown GM wheat as MON71200 event, one of two events obtained with the pMON25497 vector. The two construct-specific assays targeted the junctions of the RactInt1 and the CP4-EPSPS elements and the CP4-EPSPS and tNOS elements, while the event-specific assay was located at the 3' junction into the wheat genome. Both construct-specific and event-specific assays had limits of detection of 0.10% of MON71200 in a seed pool. As expected, the two construct-specific assays cross-reacted with other wheat and corn events containing the same elements in the same order. No cross-reactivity was observed for the event-specific assay. The integrated strategy employed in this study can serve as a model for other cases when facing similar challenges involving unknown GM events.

Reproducibility of next-generation-sequencing-based analysis of a CRISPR/Cas9 genome edited oil seed rape

Next-generation-sequencing (NGS) becomes increasingly important for laboratories tasked with the detection of genetically modified organisms (GMOs) in food, feed and seeds. Its implementation into standardized workflows demands reliable intra- and inter-laboratory reproducibility. Here, we analyze the reproducibility of short- and long-read targeted NGS and long-read whole genome sequencing (WGS) data between three independent laboratories. Replicate samples were submitted for sequencing and comparatively analyzed. The targeted-NGS-samples consisted of oil seed rape (OSR) sampled from a commodity shipment spiked with a genome edited (GE) OSR and the WGS-samples consisted of leaf material from the GMOs' parental line. All laboratories delivered highly reproducible high-quality targeted NGS data with little variation. The detection of GMO-related sequences works well regardless of the facility, while the mapping to the complex genome is superior using long read data. Long read WGS is currently not suitable for routine use in enforcement laboratories, due to a large inter-laboratory variation.

A Machine Learning Method for Genome Engineering Design Tool Attribution

As the ability to engineer biological systems improves with increasingly advanced technology, the risk of accidental or intentional release of a dangerous genetically modified organism becomes greater. It is important that authorities can carry out attribution for the source of a genetically modified biological agent release. In the absence of evidence that ties a release directly to the individuals responsible, attribution can be carried out in part by discovering the in silico tools used to design the engineered genetic components, which can leave a signature in the DNA of the organism. Previous attribution methods have focused on identifying the laboratory of origin of an engineered organism using machine learning on plasmid signatures. The next logical step is to address attribution using signatures from the tools that are used to create the engineered modifications. A random forest classifier was developed that discriminates between design tools used to optimize coding regions for incorporation into the genome of another organism. To this end, tens of thousands of genes were optimized with 4 different codon optimization methods and relevant features from these sequences were generated for a machine learning classifier. This method achieves more than 97% accuracy in predicting which tools were used to design codon optimized genes for expression in other organisms. The methods presented here lay the groundwork for the creation of effective organism engineering attribution techniques. Such methods can act both as deterrents for future attempts at creating dangerous organisms as well as tools for forensic science.

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

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 GMO 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.

Development and performance evaluation of whole-genome sequencing with paired-end and mate-pair strategies in molecular characterization of GM crops: One GM rice 114-7-2 line as an example

Basic data for the safety assessment of transgenic line involves the molecular characterization of the integration site of exogenous DNA, flanking sequences, copy number, and unintended plasmid backbone residues. However, performing a full molecular characterization remains challenging, especially for GMOs that possess complex exogenous DNA integrations. We established two whole-genome sequencing strategies: paired-end and mate-pair, to characterize the exogenous DNA integration of a human serum albumin gene into rice line 114-7-2, and evaluated the performance of these two strategies in the molecular characterization of transgenic line. The results showed the existence of two exogenous DNA insertion loci (Chr 01 and Chr 04) and their corresponding flanking sequences, five copies of the exogenous rHSA gene, and the presence of unintended residual plasmid backbone sequences. However, the WGS-MP strategy demonstrated higher efficiency, lower cost, and lower background noise compared with the WGS-PE analysis, especially for identification of the exogenous DNA integration site.

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.

LIFE-Seq: a universal Large Integrated DNA Fragment Enrichment Sequencing strategy for deciphering the transgene integration of genetically modified organisms

Molecular characterization of genetically modified organisms (GMOs) yields basic information on exogenous DNA integration, including integration sites, entire inserted sequences and structures, flanking sequences and copy number, providing key data for biosafety assessment. However, there are few effective methods for deciphering transgene integration, especially for large DNA fragment integration with complex rearrangement, inversion and tandem repeats. Herein, we developed a universal Large Integrated DNA Fragments Enrichment strategy combined with PacBio Sequencing (LIFE-Seq) for deciphering transgene integration in GMOs. Universal tilling DNA probes targeting transgenic elements and exogenous genes facilitate specific enrichment of large inserted DNA fragments associated with transgenes from plant genomes, followed by PacBio sequencing. LIFE-Seq were evaluated using six GM events and four crop species. Target DNA fragments averaging ~6275 bp were enriched and sequenced, generating ~26 352 high fidelity reads for each sample. Transgene integration structures were determined with high repeatability and sensitivity. Compared with next-generation whole-genome sequencing, LIFE-Seq achieved better data integrity and accuracy, greater universality and lower cost, especially for transgenic crops with complex inserted DNA structures. LIFE-Seq could be applied in molecular characterization of transgenic crops and animals, and complex DNA structure analysis in genetics research.

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.

Event-specific quantitative polymerase chain reaction methods for detection of double-herbicide-resistant genetically modified corn MON 87419 based on the 3'-junction of the insertion site

MON 87419 was one of the new transgenic corn events developed in US with the trait of herbicide resistance to both dicamba and glyphosate. To monitor unintended release of genetically modified organism in the future, as well as to meet GM-labeling requirements, it is requisite to develop a reliable method for the detection and quantification of MON 87419, an event-specific primer pair was designed to amplify the 3'-junction site between the endogenous genome sequence and the transferred DNA of GM event MON 87419, amplicons of desired size were produced by qualitative polymerase chain reaction (PCR) assay. For the validation of this quantitative method, the mixed samples containing 10%, 1%, and 0.1% MON 87419 ingredient were quantified. The precisions were expressed as relative standard deviations, deviated by 7.87%, 12.94%, and 19.98%, respectively. These results clearly demonstrate that the PCR methods we developed herein can be used for event-specific quantitative testing of the double-herbicide-resistant corn MON 87419.

Comparison of qPCR and Metabarcoding Methods as Tools for the Detection of Airborne Inoculum of Forest Fungal Pathogens

Forest diseases caused by invasive fungal pathogens are becoming more common, sometimes with dramatic consequences to forest ecosystems. The development of early detection systems is necessary for efficient surveillance and to mitigate the impact of invasive pathogens. Windborne spores are an important pathway for introduction of fungal pathogens into new areas; the design of spore trapping devices adapted to forests, capable of collecting different types of spores, and aligned with development of efficient molecular methods for detection of the pathogen, should help forest managers anticipate new disease outbreaks. Two types of Rotorod samplers were evaluated for the collection of airborne inoculum of forest fungal pathogens with a range of spore sizes in five forest types. Detection was by specific quantitative PCR (qPCR) and by high-throughput sequencing (HTS) of amplified internal transcribed spacer sequences using a new bioinformatic pipeline, FungiSearch, developed for diagnostic purposes. Validation of the pipeline was conducted on mock communities of 10 fungal species belonging to different taxa. Although the sensitivity of the new HTS pipeline was lower than the specific qPCR, it was able to detect a wide variety of fungal pathogens. FungiSearch is easy to use, and the reference database is updatable, making the tool suitable for rapid identification of new pathogens. This new approach combining spore trapping and HTS detection is promising as a diagnostic tool for invasive fungal pathogens.

Distribution and degradation of DNA from non-genetically and genetically modified soybean (Roundup Ready): Impact of soybean protein concentrate and soybean protein isolate preparation

To improve genetically modified product labelling legislation and promote the development of genetically modified foods, the mass variations of genomic DNA and length distributions of DNA fragments in non-genetically and genetically modified soybean (Roundup Ready) and the variations in transgenic contents during soybean protein concentrate (SPC) and soybean protein isolate (SPI) preparation were monitored. The material masses throughout the process conformed to the law of mass conservation, and amounts of DNA recovered decreased with SPC and SPI preparation. The successive steps of ethanol extraction led to a decrease in the size distribution of the recovered DNA. For the LECTIN and CP4 EPSPS targets investigated, longer fragments (>800 bp) were more affected than smaller fragments (<200 bp), and both targets degraded slowly upon degradation into small fragments. DNA distribution and degradation thereby affected GMO quantification. After preparation procedures, the transgenic contents of SPC and SPI products were higher than that of raw soybean.

Strategies to produce T-DNA free CRISPRed fruit trees via Agrobacterium tumefaciens stable gene transfer

Genome editing via CRISPR/Cas9 is a powerful technology, which has been widely applied to improve traits in cereals, vegetables and even fruit trees. For the delivery of CRISPR/Cas9 components into dicotyledonous plants, Agrobacterium tumefaciens mediated gene transfer is still the prevalent method, although editing is often accompanied by the integration of the bacterial T-DNA into the host genome. We assessed two approaches in order to achieve T-DNA excision from the plant genome, minimizing the extent of foreign DNA left behind. The first is based on the Flp/FRT system and the second on Cas9 and synthetic cleavage target sites (CTS) close to T-DNA borders, which are recognized by the sgRNA. Several grapevine and apple lines, transformed with a panel of CRISPR/SpCas9 binary vectors, were regenerated and characterized for T-DNA copy number and for the rate of targeted editing. As detected by an optimized NGS-based sequencing method, trimming at T-DNA borders occurred in 100% of the lines, impairing in most cases the excision. Another observation was the leakage activity of Cas9 which produced pierced and therefore non-functional CTS. Deletions of genomic DNA and presence of filler DNA were also noticed at the junctions between T-DNA and genomic DNA. This study proved that many factors must be considered for designing efficient binary vectors capable of minimizing the presence of exogenous DNA in CRISPRed fruit trees.

Targeted MinION sequencing of transgenes

The presence of genetically modified organisms (GMO) is commonly assessed using real-time PCR methods targeting the most common transgenic elements found in GMOs. Once the presence of GM material has been established using these screening methods, GMOs are further identified using a battery of real-time PCR methods, each being specific of one GM event and usually targeting the junction of the plant genome and of the transgenic DNA insert. If, using these specific methods, no GMO could be identified, the presence of an unauthorized GMO is suspected. In this context, the aim of this work was to develop a fast and simple method to obtain the sequence of the transgene and of its junction with plant DNA, with the presence of a screening sequence as only prior knowledge. An unauthorized GM petunia, recently found on the French market, was used as template during the development of this new molecular tool. The innovative proposed protocol is based on the circularization of fragmented DNA followed by the amplification of the transgene and of its flanking regions using long-range inverse PCR. Sequencing was performed using the Oxford Nanopore MinION technology and a bioinformatic pipeline was developed.