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

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

Analysis of the Candidate Genes and Underlying Molecular Mechanism of P198, an RNAi-Related Dwarf and Sterile Line

The genome-wide long hairpin RNA interference (lhRNAi) library is an important resource for plant gene function research. Molecularly characterizing lhRNAi mutant lines is crucial for identifying candidate genes associated with corresponding phenotypes. In this study, a dwarf and sterile line named P198 was screened from the Brassica napus (B. napus) RNAi library. Three different methods confirmed that eight copies of T-DNA are present in the P198 genome. However, only four insertion positions were identified in three chromosomes using fusion primer and nested integrated polymerase chain reaction. Therefore, the T-DNA insertion sites and copy number were further investigated using Oxford Nanopore Technologies (ONT) sequencing, and it was found that at least seven copies of T-DNA were inserted into three insertion sites. Based on the obtained T-DNA insertion sites and hairpin RNA (hpRNA) cassette sequences, three candidate genes related to the P198 phenotype were identified. Furthermore, the potential differentially expressed genes and pathways involved in the dwarfism and sterility phenotype of P198 were investigated by RNA-seq. These results demonstrate the advantage of applying ONT sequencing to investigate the molecular characteristics of transgenic lines and expand our understanding of the complex molecular mechanism of dwarfism and male sterility in B. napus.

Triple-transgenic soybean in conjunction with glyphosate drive patterns in the rhizosphere microbial community assembly

Plant roots continuously influence the rhizosphere, which also serves as a recruitment site for microorganisms with desirable functions. The development of genetically engineered (GE) crop varieties has offered unparalleled yield advantages. However, in-depth research on the effects of GE crops on the rhizosphere microbiome is currently insufficient. We used a triple-transgenic soybean cultivar (JD606) that is resistant to insects, glyphosate, and drought, along with its control, ZP661, and JD606 treated with glyphosate (JD606G). Using 16S and ITS rDNA sequencing, their effects on the taxonomy and function of the bacterial and fungal communities in the rhizosphere, surrounding, and bulk soil compartment niches were determined. Alpha diversity demonstrated a strong influence of JD606 and JD606G on bacterial Shannon diversity. Both treatments significantly altered the soil's pH and nitrogen content. Beta diversity identified the soil compartment niche as a key factor with a significant probability of influencing the bacterial and fungal communities associated with soybeans. Further analysis showed that the rhizosphere effect had a considerable impact on bacterial communities in JD606 and JD606G soils but not on fungal communities. Microbacterium, Bradyrhizobium, and Chryseobacterium were found as key rhizobacterial nodes. In addition, the LEfSe analysis identified biomarker taxa with plant-beneficial attributes, demonstrating rhizosphere-driven microbial recruitment. FUNGuild, Bugbase, and FAPROTAX functional predictions showed that ZP661 soils had more plant pathogen-associated microbes, while JD606 and JD606G soils had more stress-tolerance, nitrogen, and carbon cycle-related microbes. Bacterial rhizosphere networks had more intricate topologies than fungal networks. Furthermore, correlation analysis revealed that the bacteria and fungi with higher abundances exhibited varying degrees of positive and negative correlations. Our findings shed new light on the niche partitioning of bacterial and fungal communities in soil. It also indicates that following triple-transgenic soybean cultivation and glyphosate application, plant roots recruit microbes with beneficial taxonomic and functional traits in the rhizosphere.

Identification of T-DNA structure and insertion site in transgenic crops using targeted capture sequencing

The commercialization of GE crops requires a rigorous safety assessment, which includes a precise DNA level characterization of inserted T-DNA. In the past, several strategies have been developed for identifying T-DNA insertion sites including, Southern blot and different PCR-based methods. However, these methods are often challenging to scale up for screening of dozens of transgenic events and for crops with complex genomes, like potato. Here, we report using target capture sequencing (TCS) to characterize the T-DNA structure and insertion sites of 34 transgenic events in potato. This T-DNA is an 18 kb fragment between left and right borders and carries three resistance (R) genes (RB, Rpi-blb2 and Rpi-vnt1.1 genes) that result in complete resistance to late blight disease. Using TCS, we obtained a high sequence read coverage within the T-DNA and junction regions. We identified the T-DNA breakpoints on either ends for 85% of the transgenic events. About 74% of the transgenic events had their T-DNA with 3R gene sequences intact. The flanking sequences of the T-DNA were from the potato genome for half of the transgenic events, and about a third (11) of the transgenic events have a single T-DNA insertion mapped into the potato genome, of which five events do not interrupt an existing potato gene. The TCS results were confirmed using PCR and Sanger sequencing for 6 of the best transgenic events representing 20% of the transgenic events suitable for regulatory approval. These results demonstrate the wide applicability of TCS for the precise T-DNA insertion characterization in transgenic crops.

Discrepancies in rhizobacterial assembly caused by glyphosate application and herbicide-tolerant soybean Co-expressing GAT and EPSPS

There are concerns that the innovation of genetically modified herbicide-tolerant (GMHT) plants, as well as the application of herbicide to such GMHT plants, could have an impact on ecological interactions and unintentionally harm non-targeted organisms. Consequently, we intend to use full-length 16 S rDNA amplicon sequencing to examine changes in the bacterial community in the rhizosphere of GMHT soybean (Z106) harboring 5-enolpyruvylshikimate-3-phosphate synthase and Glyphosate N-acetyltransferase genes and GMHT soybean treated with glyphosate (Z106G). Glyphosate application significantly impacted bacterial alpha diversity (species richness, and Shannon diversity). Permutational multivariate analysis of variance of beta diversity demonstrated that soil compartments and growth stages had a substantial impact on soybean rhizobacterial communities (soil compartments, growth stages, P = 0.001). Community composition revealed that Z106G soils were abundant in Taibaiella and Arthrobacter pascens at maturity, while Chryseobacterium joostei and Stenotrophomonas maltophilia predominated in Z106 soils during flowering. Nitrogen-fixing and phosphate-solubilizing microbes were found in higher proportions in the rhizosphere than in bulk soil, with Sinorhizobium being more abundant in Z106 and Bacillus and Stenotrophomonas being more prevalent in Z106G rhizosphere soils. Collectively, our findings suggest glyphosate application and glyphosate-tolerant soybean as potential regulators of soybean rhizobacterial composition.

Molecular Characterization and Efficacy Evaluation of Transgenic Maize Harboring cry2Ab-vip3A-cp4epsps for Insect Resistance and Herbicide Tolerance

Insect infestation and weed interference have a seriously negative impact on the growth, yield, and grain quality of maize. In this study, transgenic maize plants harboring three exogenous genes, cry2Ab, vip3A, and cp4epsps, that were constructed into a single T-DNA were developed for protection against insects and weeds. The transgene integration sites on the chromosomes in two transgenic maize events, CVC-1 and CVC-2, were determined using whole genome sequencing and specific PCR detection. As revealed by laboratory insect bioassays, these two transgenic events exhibited strong insecticidal toxicity against three major species of Lepidoptera insects, including Mythimna separata, Helicoverpa armigera, and Spodoptera frugiperda, with mortality rates exceeding 96%, 100%, and 100%, respectively, after six days of infestation. In addition, CVC-1 exhibited a high tolerance to glyphosate under field conditions. The successful expressions of cry2Ab, vip3A, and cp4epsps in various tissues at different developmental stages of CVC-1 were validated at the transcriptional and translational levels using quantitative real-time reverse transcription PCR (qRT-PCR) and enzyme-linked immunosorbent assay (ELISA), respectively. These findings demonstrated that the transgenic maize CVC-1 developed using this triple gene construct has excellent insect resistance and herbicide tolerance, which may provide a valuable germplasm resource and data support for future maize breeding of insect and weed control.

Review of CRISPR/Cas Systems on Detection of Nucleotide Sequences

Nowadays, with the rapid development of biotechnology, the CRISPR/Cas technology in particular has produced many new traits and products. Therefore, rapid and high-resolution detection methods for biotechnology products are urgently needed, which is extremely important for safety regulation. Recently, in addition to being gene editing tools, CRISPR/Cas systems have also been used in detection of various targets. CRISPR/Cas systems can be successfully used to detect nucleic acids, proteins, metal ions and others in combination with a variety of technologies, with great application prospects in the future. However, there are still some challenges need to be addressed. In this review, we will list some detection methods of genetically modified (GM) crops, gene-edited crops and single-nucleotide polymorphisms (SNPs) based on CRISPR/Cas systems, hoping to bring some inspiration or ideas to readers.

The host niches of soybean rather than genetic modification or glyphosate application drive the assembly of root-associated microbial communities

Plant roots significantly influence soil microbial diversity, and soil microorganisms play significant roles in both natural and agricultural ecosystems. Although the genetically modified (GM) crops with enhanced insect and herbicide resistance are thought to have unmatched yield and stress resistance advantages, thorough and in-depth case studies still need to be carried out in a real-world setting due to the potential effects of GM plants on soil microbial communities. In this study, three treatments were used: a recipient soybean variety Jack, a triple transgenic soybean line JD321, and the glyphosate-treated JD321 (JD321G). Three sampling stages (flowering, seed filling and maturing), as well as three host niches of soybean rhizosphere [intact roots (RT), rhizospheric soil (RS) and surrounding soil (SS)] were established. In comparison to Jack, the rhizospheric soil of JD321G had higher urease activity and lower nitrite reductase at the flowering stage. Different treatments and different sampling stages existed no significant effects on the compositions of microbial communities at different taxonomic levels. However, at the genus level, the relative abundance of three plant growth-promoting fungal genera (i.e. Mortierella, Chaetomium and Pseudombrophila) increased while endophytic bacteria Chryseobacterium and pathogenic bacteria Streptomyces decreased from the inside to the outside of the roots (i.e. RT → RS → SS). Moreover, two bacterial genera, Bradyrhizobium and Ensifer were more abundant in RT than in RS and SS, as well as three species, Agrobacterium radiobacter, Ensifer fredii and Ensifer meliloti, which are closely related to nitrogen-fixation. Furthermore, five clusters of orthologous groups (COGs) associated to nitrogen-fixation genes were higher in RT than in RS, whereas only one COG annotated as dinitrogenase iron-molybdenum cofactor biosynthesis protein was lower. Overall, the results imply that the rhizosphere host niches throughout the soil-plant continuum largely control the composition and function of the root-associated microbiome of triple transgenic soybean.

Comparative Transcriptome Profiling Reveals Potential Candidate Genes, Transcription Factors, and Biosynthetic Pathways for Phosphite Response in Potato (Solanum tuberosum L.)

The study was conducted with C31 and C80 genotypes of the potato (Solanum tuberosum L.), which are tolerant and susceptible to phosphite (Phi, H₂PO₃), respectively. To decipher the molecular mechanisms underlying tolerance and susceptibility to Phi in the potato, RNA sequencing was used to study the global transcriptional patterns of the two genotypes. Media were prepared with 0.25 and 0.50 mM Phi, No-phosphorus (P), and 1.25 mM (phosphate, Pi as control). The values of fragments per kilobase of exon per million mapped fragments of the samples were also subjected to a principal component analysis, grouping the biological replicates of each sample. Using stringent criteria, a minimum of 819 differential (DEGs) were detected in both C80-Phi-0.25_vs_C80-Phi-0.50 (comprising 517 upregulated and 302 downregulated) and C80-Phi-0.50_vs_C80-Phi-0.25 (comprising 302 upregulated and 517 downregulated) and a maximum of 5214 DEGs in both C31-Con_vs_C31-Phi-0.25 (comprising 1947 upregulated and 3267 downregulated) and C31-Phi-0.25_vs_C31-Con (comprising 3267 upregulated and 1947 downregulated). DEGs related to the ribosome, plant hormone signal transduction, photosynthesis, and plant–pathogen interaction performed important functions under Phi stress, as shown by the Kyoto Encyclopedia of Genes and Genomes annotation. The expressions of transcription factors increased significantly in C31 compared with C80. For example, the expressions of Soltu.DM.01G047240, Soltu.DM.08G015900, Soltu.DM.06G012130, and Soltu.DM.08G012710 increased under P deficiency conditions (Phi-0.25, Phi-0.50, and No-P) relative to the control (P sufficiency) in C31. This study adds to the growing body of transcriptome data on Phi stress and provides important clues to the Phi tolerance response of the C31 genotype.

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.

A paired-end whole-genome sequencing approach enables comprehensive characterization of transgene integration in rice

Efficient, accurate molecular characterization of genetically modified (GM) organisms is challenging, especially for those transgenic events transferred with genes/elements of recipient species. Herein, we decipher the comprehensive molecular characterization of one novel GM rice event G281 which was transferred with native promoters and an RNA interference (RNAi) expression cassette using paired-end whole genome sequencing (PE-WGS) and modified TranSeq approach. Our results show that transgenes integrate at rice chromosome 3 locus 16,439,674 included a 36 bp deletion of rice genomic DNA, and the whole integration contains two copies of the complete transfer DNA (T-DNA) in a head-to-head arrangement. No unintended insertion or backbone sequence of the transformed plasmid is observed at the whole genome level. Molecular characterization of the G281 event will assist risk assessment and application for a commercial license. In addition, we speculate that our approach could be further used for identifying the transgene integration of cisgenesis/intragenesis crops since both ends of T-DNA in G281 rice were from native gene or elements which is similar with that of cisgenesis/intrasgenesis. Our results from the in silico mimicking cisgenesis event confirm that the mimic rice Gt1 gene insertion and its flanking sequences are successfully identified, demonstrating the applicability of PE-WGS for molecular characterization of cisgenesis/intragenesis crops.

Coupling paired-end whole-genome sequencing with droplet digital PCR enabled precise identification of a transgene insertion in the genetically modified rice event G281 on chromosome 3 and the potential for exploring the native gene integration.

Generation and Identification of the Number of Copies of Exogenous Genes and the T-DNA Insertion Site in SCN-Resistance Transformation Event ZHs1-2

Soybean cyst nematode (SCN, Heterodera glycines Ichinohe) causes an estimated economic loss of about USD 3 billion each year in soybean (Glycine max L.) production worldwide. Overexpression of resistance genes against SCN provides a powerful approach to develop SCN resistance cultivars in soybean. The clarification of molecular characterization in transformation events is a prerequisite for ecological risk assessment, food safety, and commercial release of genetically modified crops. Here, we generated transgenic events harboring the BCN (beet cyst nematode) resistance Hs1pro-1 gene using the Agrobacterium-mediated method in soybean, evaluated their resistance to SCN infection, and clarified the molecular characterization of one of the transformation events. Five independent and stable inheritable transformation events were generated by an Agrobacterium-mediated transformation method. SCN resistance tests showed the average number of developed females per plant and female index (FI) in T4 ZHs1-1, ZHs1-2, ZHs1-3, ZHs1-4, and ZHs1-5 transformation events were significantly lower than that in the nontransgenic control. Among these, the ZHs1-2 transformation event had the lowest number of developed females per plant and FI. Southern hybridization showed the exogenous target Hs1pro-1 gene was inserted in one copy and the Bar gene was inserted two copies in the ZHs1-2 transformation event. The exogenous T-DNA fragment was integrated in the reverse position of Chr02: 5351566-5231578 (mainly the Bar gene expression cassette) and in the forward position of Chr03: 17083358-17083400 (intact T-DNA, including Hs1pro-1 and Bar gene expression cassette) using a whole genome sequencing method (WGS). The results of WGS method and Southern hybridization were consistent. All the functional elements of exogenous T-DNA fragments were verified by PCR using specific primer pairs in the T5 and T6 ZHs1-2 transformation events. These results demonstrated that the overexpression of Hs1pro-1 gene enhanced SCN resistance, and provide an important reference for the biosafety assessment and the labeling detection in transformation event ZHs1-2.

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.

Activation of the VQ Motif-Containing Protein Gene VQ28 Compromised Nonhost Resistance of Arabidopsis thaliana to Phytophthora Pathogens

Nonhost resistance refers to resistance of a plant species to all genetic variants of a non-adapted pathogen. Such resistance has the potential to become broad-spectrum and durable crop disease resistance. We previously employed Arabidopsis thaliana and a forward genetics approach to identify plant mutants susceptible to the nonhost pathogen Phytophthora sojae, which resulted in identification of the T-DNA insertion mutant esp1 (enhanced susceptibility to Phytophthora). In this study, we report the identification of VQ motif-containing protein 28 (VQ28), whose expression was highly up-regulated in the mutant esp1. Stable transgenic A. thaliana plants constitutively overexpressing VQ28 compromised nonhost resistance (NHR) against P. sojae and P. infestans, and supported increased infection of P. parasitica. Transcriptomic analysis showed that overexpression of VQ28 resulted in six differentially expressed genes (DEGs) that are involved in the response to abscisic acid (ABA). High performance liquid chromatography-mass spectrometry (HPLC-MS) detection showed that the contents of endogenous ABA, salicylic acid (SA), and jasmonate (JA) were enriched in VQ28 overexpression lines. These findings suggest that overexpression of VQ28 may lead to an imbalance in plant hormone homeostasis. Furthermore, transient overexpression of VQ28 in Nicotiana benthamiana rendered plants more susceptible to Phytophthora pathogens. Deletion mutant analysis showed that the C-terminus and VQ-motif were essential for plant susceptibility. Taken together, our results suggest that VQ28 negatively regulates plant NHR to Phytophthora pathogens.

Using Combined Methods of Genetic Mapping and Nanopore-Based Sequencing Technology to Analyze the Insertion Positions of G10evo-EPSPS and Cry1Ab/Cry2Aj Transgenes in Maize

The insertion position of the exogenous fragment sequence in a genetically modified organism (GMO) is important for the safety assessment and labeling of GMOs. SK12-5 is a newly developed transgenic maize line transformed with two trait genes [i.e., G10evo-5-enolpyrul-shikimate-3-phosphate synthase (EPSPS) and Cry1Ab/Cry2Aj] that was recently approved for commercial use in China. In this study, we tried to determine the insertion position of the exogenous fragment for SK12-5. The transgene-host left border and right border integration junctions were obtained from SK12-5 genomic DNA by using the thermal asymmetric interlaced polymerase chain reaction (TAIL-PCR) and next-generation Illumina sequencing technology. However, a Basic Local Alignment Search Tool (BLAST) analysis revealed that the flanking sequences in the maize genome are unspecific and that the insertion position is located in a repetitive sequence area in the maize genome. To locate the fine-scale insertion position in SK12-5, we combined the methods of genetic mapping and nanopore-based sequencing technology. From a classical bulked-segregant analysis (BSA), the insertion position in SK12-5 was mapped onto Bin9.03 of chromosome 9 between the simple sequence repeat (SSR) markers umc2337 and umc1743 (26,822,048-100,724,531 bp). The nanopore sequencing results uncovered 10 reads for which one end was mapped onto the vector and the other end was mapped onto the maize genome. These observations indicated that the exogenous T-DNA fragments were putatively integrated at the position from 82,329,568 to 82,379,296 bp of chromosome 9 in the transgenic maize SK12-5. This study is helpful for the safety assessment of the novel transgenic maize SK12-5 and shows that the combined method of genetic mapping and the nanopore-based sequencing technology will be a useful approach for identifying the insertion positions of transgenic sequences in other GM plants with relatively large and complex genomes.

Measurement of Transgene Copy Number in Plants Using Droplet Digital PCR

Transgenic plants are produced both to investigate gene function and to confer desirable traits into crops. Transgene copy number is known to influence expression levels, and consequently, phenotypes. Similarly, knowledge of transgene zygosity is desirable for making quantitative assessments of phenotype and tracking the inheritance of transgenes in progeny generations. Since the first transgenic plants were produced, several methods for determining copy number have been applied, including Southern blotting, quantitative real-time PCR, and more recently, sequencing methods; however, each method has specific disadvantages, compromising throughput, accuracy, or expense. Digital PCR (dPCR) divides reactions into partitions, converting the exponential, analogue nature of PCR into a linear, digital signal that allows the frequency of occurrence of specific sequences to be accurately estimated. Confidence increases with the number of partitions; therefore, the availability of emulsion technologies that enable reactions to be divided into tens of thousands of nanodroplets allows accurate determination of copy number in what has become known as digital droplet PCR (ddPCR). ddPCR offers similar benefits of low costs and scalability as other PCR techniques but with superior accuracy and reliability. Graphic abstract: Digital PCR (dPCR) divides reactions into partitions, converting the exponential, analogue nature of PCR into a linear, digital signal that allows the frequency of transgene copy number to be accurately assessed.

Whole-genome resequencing using next-generation and Nanopore sequencing for molecular characterization of T-DNA integration in transgenic poplar 741

BACKGROUND

The molecular characterization information of T-DNA integration is not only required by public risk assessors and regulators, but is also closely related to the expression of exogenous and endogenous genes. At present, with the development of sequencing technology, whole-genome resequencing has become an attractive approach to identify unknown genetically modified events and characterise T-DNA integration events.

RESULTS

In this study, we performed genome resequencing of Pb29, a transgenic high-resistance poplar 741 line that has been commercialized, using next-generation and Nanopore sequencing. The results revealed that there are two T-DNA insertion sites, located at 9,283,905-9,283,937 bp on chromosome 3 (Chr03) and 10,868,777-10,868,803 bp on Chr10. The accuracy of the T-DNA insertion locations and directions was verified using polymerase chain reaction amplification. Through sequence alignment, different degrees of base deletions were detected on the T-DNA left and right border sequences, and in the flanking sequences of the insertion sites. An unknown fragment was inserted between the Chr03 insertion site and the right flanking sequence, but the Pb29 genome did not undergo chromosomal rearrangement. It is worth noting that we did not detect the API gene in the Pb29 genome, indicating that Pb29 is a transgenic line containing only the BtCry1AC gene. On Chr03, the insertion of T-DNA disrupted a gene encoding TAF12 protein, but the transcriptional abundance of this gene did not change significantly in the leaves of Pb29. Additionally, except for the gene located closest to the T-DNA integration site, the expression levels of four other neighboring genes did not change significantly in the leaves of Pb29.

CONCLUSIONS

This study provides molecular characterization information of T-DNA integration in transgenic poplar 741 line Pb29, which contribute to safety supervision and further extensive commercial planting of transgenic poplar 741.

Efficient identification of genomic insertions and flanking regions through whole-genome sequencing in three transgenic soybean events

Genomic insertions and flanking regions of transgenes in host genomes constitute a critical component of precise molecular characterization and event-specific detection, which are required in the development and assessment for regulatory approval of genetically modified (GM) crops. Previously, we reported three transgenic soybean events harboring the inverted repeats of the soybean mosaic virus NIb (nuclear inclusion b) gene, exhibiting significantly enhanced resistance to multiple Potyvirus strains. To facilitate safety assessment and event-specific detection, we identified the transgene insertion sites and flanking sequences of the events L120, L122, and L123 using whole-genome sequencing. More than 14.48 Gb sequence data (13 × coverage) were generated using the Illumina HiSeq Xten platform for each event. The sequence reads corresponding to boundaries of inserted T-DNA, and associated native flanking sequences were identified by bioinformatic comparison with the soybean reference genome (Wm82.a2.v1) and the transformation vector sequence. The results indicated that two T-DNA insertions occurred in L120, on Chr07 and Chr13, while L122 and L123 showed single insertions, on Chr02 and Chr06, respectively. Based on the flanking sequences of the inserted T-DNA, the event-specific detection for each event was established using specific PCR primers, and PCR amplification followed by sequencing of PCR products further confirmed the putative insertion loci and flanking regions in the transgenic lines. Our results demonstrate the efficacy and robustness of whole-genome sequencing in identifying the genomic insertions and flanking regions in GM crops. Moreover, the characterization of insertion loci and the establishment of event-specific detection will facilitate the application and development of broad-spectrum virus-resistant transgenic soybean cultivars.

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.

Whole-Genome Sequencing: An Effective Strategy for Insertion Information Analysis of Foreign Genes in Transgenic Plants

Molecular characterization is a key step in the risk assessment of genetically modified organisms (GMOs) for regulatory approval. Herein, we describe a method for analyzing copy number, insertion loci, and flanking sequences through whole-genome sequencing (WGS) and bioinformatics. Comprehensive molecular characterization of G2-6 transgenic rice was performed using this pipeline. The results showed that one copy of the foreign gene was inserted into rice chromosome 8. There was no vector backbone insertion but an unexpected insertion and DNA rearrangement at the 3' end of the T-DNA. We also obtained the 5' and 3' flanking sequences of the T-DNA. Our results suggested that the use of a combination of WGS and bioinformatics is an effective strategy for the molecular characterization of GMOs.

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.

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.

Impact of a G2-EPSPS & GAT Dual Transgenic Glyphosate-Resistant Soybean Line on the Soil Microbial Community under Field Conditions Affected by Glyphosate Application

In the past thirty years, the biosafety of the aboveground part of crops, including horizontal gene transferal through pollen dispersal and hybridization, has been the focus of research; however, microbial communities in the underground part are attracting increasing attention. In the present study, the soybean root-associated bacterial communities of the G2-EPSPS plus GAT transgenic soybean line Z106, its recipient variety ZH10, and Z106 treated with glyphosate (Z106J) were compared at the seedling, flowering, and seed filling stages by high-throughput sequencing of the V4 hypervariable regions of 16S rRNA gene amplicons using Illumina MiSeq. The results obtained showed no significant differences in the alpha/beta diversities of root-associated bacterial communities at the three stages among ZH10, Z106, and Z106J under field growth conditions; however, the relative abundance of four main nitrogen-fixing bacterial genera significantly differed among ZH10, Z106, and Z106J. Ternary plot results indicated that in the root compartment, the proportional contributions of rhizobial nitrogen-fixing Ensifer fredii and Bradyrhizobium elkanii, which exhibit an extremely broad nodulation host range, markedly differed among the three treatments at the three stages. Thus, the present results indicate that transgenic G2-EPSPS and GAT soybean may induce different changes in functional bacterial species in soil, such as E. fredii and B. elkanii, from ZH10, which were compensated for/enriched at the flowering and seed filling stages, respectively, to some extent through as of yet unknown mechanisms by transgenic soybean treated with glyphosate.

Generation of Epichloë Strains Expressing Fluorescent Proteins Suitable for Studying Host-Endophyte Interactions and Characterisation of a T-DNA Integration Event

Methods for the identification and localisation of endophytic fungi are required to study the establishment, development, and progression of host-symbiont interactions, as visible reactions or disease symptoms are generally absent from host plants. Fluorescent proteins have proved valuable as reporter gene products, allowing non-invasive detection in living cells. This study reports the introduction of genes for two fluorescent proteins, green fluorescent protein (GFP) and red fluorescent protein, DsRed, into the genomes of two distinct perennial ryegrass (Lolium perenne L.)-associated Epichloë endophyte strains using A. tumefaciens-mediated transformation. Comprehensive characterisation of reporter gene-containing endophyte strains was performed using molecular genetic, phenotypic, and bioinformatic tools. A combination of long read and short read sequencing of a selected transformant identified a single complex T-DNA insert of 35,530 bp containing multiple T-DNAs linked together. This approach allowed for comprehensive characterisation of T-DNA integration to single-base resolution, while revealing the unanticipated nature of T-DNA integration in the transformant analysed. These reporter gene endophyte strains were able to establish and maintain stable symbiotum with the host. In addition, the same endophyte strain labelled with two different fluorescent proteins were able to cohabit the same plant. This knowledge can be used to provide the basis to develop strategies to gain new insights into the host-endophyte interaction through independent and simultaneous monitoring in planta throughout its life cycle in greater detail.

Mapping of transgenic alleles in soybean using a nanopore-based sequencing strategy

Transgenic technology was developed to introduce transgenes into various organisms to validate gene function and add genetic variations >40 years ago. However, the identification of the transgene insertion position is still challenging in organisms with complex genomes. Here, we report a nanopore-based method to map the insertion position of a Ds transposable element originating in maize in the soybean genome. In this method, an oligo probe is used to capture the DNA fragments containing the Ds element from pooled DNA samples of transgenic soybean plants. The Ds element-enriched DNAs are then sequenced using the MinION-based platform of Nanopore. This method allowed us to rapidly map the Ds insertion positions in 51 transgenic soybean lines through a single sequencing run. This strategy is high throughput, convenient, reliable, and cost-efficient. The transgenic allele mapping protocol can be easily translated to other eukaryotes with complex genomes.

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.

Development and Event-specific Detection of Transgenic Glyphosate-resistant Rice Expressing the G2-EPSPS Gene

Glyphosate is a widely used herbicide, due to its broad spectrum, low cost, low toxicity, high efficiency, and non-selective characteristics. Rice farmers rarely use glyphosate as a herbicide, because the crop is sensitive to this chemical. The development of transgenic glyphosate-tolerant rice could greatly improve the economics of rice production. Here, we transformed the Pseudomonas fluorescens G2 5-enolpyruvyl shikimate-3-phosphate synthase (EPSPS) gene G2-EPSPS, which conferred tolerance to glyphosate herbicide into a widely used japonica rice cultivar, Zhonghua 11 (ZH11), to develop two highly glyphosate-tolerant transgenic rice lines, G2-6 and G2-7, with one exogenous gene integration. Seed germination tests and glyphosate-tolerance assays of plants grown in a greenhouse showed that the two transgenic lines could greatly improve glyphosate-tolerance compared with the wild-type; The glyphosate-tolerance field test indicated that both transgenic lines could grow at concentrations of 20,000 ppm glyphosate, which is more than 20-times the recommended concentration in the field. Isolation of the flanking sequence of transgenic rice G2-6 indicated that the 5'-terminal of T-DNA was inserted into chromosome 8 of the rice genome. An event-specific PCR test system was established and the limit of detection of the primers reached five copies. Overall, the G2-EPSPS gene significantly improved glyphosate-tolerance in transgenic rice; furthermore, it is a useful candidate gene for the future development of commercial transgenic rice.