Plant breeding involving genetic engineering does not result in unacceptable unintended effects in rice relative to conventional cross-breeding

Genetically modified soybean lines exhibit less transcriptomic variation compared to natural varieties

Genetically modified (GM) soybeans provide a huge amount of food for human consumption and animal feed. However, the possibility of unexpected effects of transgenesis has increased food safety concerns. High-throughput sequencing profiling provides a potential approach to directly evaluate unintended effects caused by foreign genes. In this study, we performed transcriptomic analyses to evaluate differentially expressed genes (DEGs) in individual soybean tissues, including cotyledon (C), germ (G), hypocotyl (H), and radicle (R), instead of using the whole seed, from four GM and three non-GM soybean lines. A total of 3,351 DEGs were identified among the three non-GM soybean lines. When the GM lines were compared with their non-GM parents, 1,836 to 4,551 DEGs were identified. Furthermore, Gene Ontology (GO) analysis of the DEGs showed more abundant categories of GO items (199) among non-GM lines than between GM lines and the non-GM natural varieties (166). Results of Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis showed that most KEGG pathways were the same for the two types of comparisons. The study successfully employed RNA sequencing to assess the differences in gene expression among four tissues of seven soybean varieties, and the results suggest that transgenes do not induce massive transcriptomic alterations in transgenic soybeans compared with those that exist among natural varieties. This work offers empirical evidence to investigate the genomic-level disparities induced by genetic modification in soybeans, specifically focusing on seed tissues.

Coronatine-Induced Maize Defense against Gibberella Stalk Rot by Activating Antioxidants and Phytohormone Signaling

One of the most destructive diseases, Gibberella stalk rot (GSR), caused by Fusarium graminearum, reduces maize yields significantly. An induced resistance response is a potent and cost-effective plant defense against pathogen attack. The functional counterpart of JAs, coronatine (COR), has attracted a lot of interest recently due to its ability to control plant growth and stimulate secondary metabolism. Although several studies have focused on COR as a plant immune elicitor to improve plant resistance to pathogens, the effectiveness and underlying mechanisms of the suppressive ability against COR to F. graminearum in maize have been limited. We investigated the potential physiological and molecular mechanisms of COR in modulating maize resistance to F. graminearum. COR treatment strongly enhanced disease resistance and promoted stomatal closure with H2O2 accumulation, and 10 μg/mL was confirmed as the best concentration. COR treatment increased defense-related enzyme activity and decreased the malondialdehyde content with enhanced antioxidant enzyme activity. To identify candidate resistance genes and gain insight into the molecular mechanism of GSR resistance associated with COR, we integrated transcriptomic and metabolomic data to systemically explore the defense mechanisms of COR, and multiple hub genes were pinpointed using weighted gene correlation network analysis (WGCNA). We discovered 6 significant modules containing 10 candidate genes: WRKY transcription factor (LOC100279570), calcium-binding protein (LOC100382070), NBR1-like protein (LOC100275089), amino acid permease (LOC100382244), glutathione S-transferase (LOC541830), HXXXD-type acyl-transferase (LOC100191608), prolin-rich extensin-like receptor protein kinase (LOC100501564), AP2-like ethylene-responsive transcription factor (LOC100384380), basic leucine zipper (LOC100275351), and glycosyltransferase (LOC606486), which are highly correlated with the jasmonic acid-ethylene signaling pathway and antioxidants. In addition, a core set of metabolites, including alpha-linolenic acid metabolism and flavonoids biosynthesis linked to the hub genes, were identified. Taken together, our research revealed differentially expressed key genes and metabolites, as well as co-expression networks, associated with COR treatment of maize stems after F. graminearum infection. In addition, COR-treated maize had higher JA (JA-Ile and Me-JA) levels. We postulated that COR plays a positive role in maize resistance to F. graminearum by regulating antioxidant levels and the JA signaling pathway, and the flavonoid biosynthesis pathway is also involved in the resistance response against GSR.

Effect of T-DNA Integration on Growth of Transgenic Populus × euramericana cv. Neva Underlying Field Stands

Multigene cotransformation has been widely used in the study of genetic improvement in crops and trees. However, little is known about the unintended effects and causes of multigene cotransformation in poplars. To gain insight into the unintended effects of T-DNA integration during multigene cotransformation in field stands, here, three lines (A1-A3) of Populus × euramericana cv. Neva (PEN) carrying Cry1Ac-Cry3A-BADH genes and three lines (B1-B3) of PEN carrying Cry1Ac-Cry3A-NTHK1 genes were used as research objects, with non-transgenic PEN as the control. Experimental stands were established at three common gardens in three locations and next generation sequencing (NGS) was used to identify the insertion sites of exogenous genes in six transgenic lines. We compared the growth data of the transgenic and control lines for four consecutive years. The results demonstrated that the tree height and diameter at breast height (DBH) of transgenic lines were significantly lower than those of the control, and the adaptability of transgenic lines in different locations varied significantly. The genotype and the experimental environment showed an interaction effect. A total of seven insertion sites were detected in the six transgenic lines, with B3 having a double-site insertion and the other lines having single copies. There are four insertion sites in the gene region and three insertion sites in the intergenic region. Analysis of the bases near the insertion sites showed that AT content was higher than the average chromosome content in four of the seven insertion sites within 1000 bp. Transcriptome analysis suggested that the differential expression of genes related to plant hormone transduction and lignin synthesis might be responsible for the slow development of plant height and DBH in transgenic lines. This study provides an integrated analysis of the unintended effects of transgenic poplar, which will benefit the safety assessment and reasonable application of genetically modified trees.

Effects of stacking breeding on the methylome and transcriptome profile of transgenic rice with glyphosate tolerance

MAIN CONCLUSION

Transcriptomics and methylomics were used to identify the potential effects resulting from GM rice breeding stacks, which provided scientific data for the safety assessment strategy of stacked GM crops in China. Gene interaction is one of the main concerns for stacked genetically modified crop safety. With the development of technology, the combination of omics and bioinformatics has become a useful tool to evaluate the unintended effects of genetically modified crops. In this study, transcriptomics and methylomics were used as molecular profiling techniques to identify the potential effects of stack through breeding. Stacked transgenic rice En-12 × Ec-26 was used as material, which was obtained through hybridization using parents En-12 and Ec-26, in which the foreign protein can form functional EPSPS protein by intein-mediated trans-splitting. Differentially methylated region (DMR) analysis showed that the effect of stacking breeding on methylation was less than that of genetic transformation at the methylome level. Differentially expressed gene (DEG) analysis showed that the DEGs between En-12 × Ec-26 and its parents were far fewer than those between transgenic rice and Zhonghua 11 (ZH11), and no unintended new genes were found in En-12 × Ec-26. Statistical analysis of gene expression and methylation involved in shikimic acid metabolism showed that there was no difference in gene expression, although there were 16 and 10 DMR genes between En-12 × Ec-26 and its parents (En and Ec) in methylation, respectively. The results indicated that the effect of stacking breeding on gene expression and DNA methylation was less than the effect of genetic transformation. This study provides scientific data supporting safety assessments of stacked GM crops in China.

Analysis of the Genetic Stability of Insect and Herbicide Resistance Genes in Transgenic Rice Lines: A Laboratory and Field Experiment

A lack of stability in the expression of Bacillus thuringiensis genes (CRY) and the dialaninophosphate resistance gene (BAR) in transgenic rice plants can lead to the loss of important characters. The genetic stability of transgenic expression in high-generation lines is thus critically important for ensuring the success of molecular breeding efforts. Here, we studied the genetic stability of resistance to insect pests and herbicides in transgenic rice lines at the molecular and phenotypic levels in a pesticide-free environment. Southern blot analysis, real-time polymerase chain reaction, and enzyme-linked immunosorbent assays revealed high stability in the copy numbers and expression levels of CRY1C, CRY2A, and BAR in transgenic lines across different generations, and gene expression levels were highly correlated with protein expression levels. The insecticide resistance of the transgenic rice lines was high. The larval mortality of Chilo suppressalis was 50.25% to 68.36% higher in transgenic lines than in non-transgenic control lines. Percent dead hearts and percent white spikelets were 16.66% to 22.15% and 27.07% to 33.47% lower in transgenic lines than in non-transgenic control lines, respectively. The herbicide resistance of the transgenic rice lines was also high. The bud length and root length ranged were 2.53 cm to 4.20 cm and 0.28 cm to 0.73 cm higher in transgenic lines than in non-transgenic control lines in the budding stage, respectively. Following application of the herbicide Basta, the chlorophyll content of the transgenic lines began to recover 2 d later in the seedling and tillering stages and 3 d later in the booting and heading stages, by contrast, the chlorophyll content of the non-transgenic lines did not recover and continued to decrease. These findings revealed high genetic stability of the resistance to insect pests and herbicides across several generations of transgenic rice regardless of the genetic background.

Comparison of the Phenotypic Performance, Molecular Diversity, and Proteomics in Transgenic Rice

The extent of molecular diversity and differentially expressed proteins (DEPs) in transgenic lines provide valuable information to understand the phenotypic performance of transgenic crops compared with their parents. Here, we compared the differences in the phenotypic variation of twelve agronomic and end-use quality traits, the extent of microsatellite diversity, and DEPs of a recurrent parent line with three transgenic rice restorer lines carrying either CRY1C gene on chromosome 11 or CRY2A gene on chromosome 12 or both genes. The three transgenic lines had significantly smaller stem borer infestation than the recurrent parent without showing significant differences among most agronomic traits, yield components, and end-use quality traits. Using 512 microsatellite markers, the three transgenic lines inherited 2.9-4.3% of the Minghui 63 donor genome and 96.3-97.1% of the CH891 recurrent parent genome. As compared with the recurrent parent, the number of upregulated and down-regulated proteins in the three transgenic lines varied from 169 to 239 and from 131 to 199, respectively. Most DEPs were associated with the secondary metabolites biosynthesis transport and catabolism, carbohydrate transport and metabolism, post-translational modification, and signal transduction mechanisms. Although several differentially expressed proteins were observed between transgenic rice and its recurrent parent, the differences may not have been associated with grain yield and most other phenotypic traits in transgenic rice.

Proteomic and Metabolomic Evaluation of Insect- and Herbicide-Resistant Maize Seeds

Label-free quantitative proteomic (LFQ) and widely targeted metabolomic analyses were applied in the safety evaluation of three genetically modified (GM) maize varieties, BBL, BFL-1, and BFL-2, in addition to their corresponding non-GM parent maize. A total of 76, 40, and 25 differentially expressed proteins (DEPs) were screened out in BBL, BFL-1, and BFL-2, respectively, and their abundance compared was with that in their non-GM parents. Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis showed that most of the DEPs participate in biosynthesis of secondary metabolites, biosynthesis of amino acids, and metabolic pathways. Metabolomic analyses revealed 145, 178, and 88 differentially accumulated metabolites (DAMs) in the BBL/ZH58, BFL-1/ZH58, and BFL-2/ZH58×CH72 comparisons, respectively. KEGG pathway enrichment analysis showed that most of the DAMs are involved in biosynthesis of amino acids, and in arginine and proline metabolism. Three co-DEPs and 11 co-DAMs were identified in the seeds of these GM maize lines. The proteomic profiling of seeds showed that the GM maize varieties were not dramatically different from their non-GM control. Similarly, the metabolomic profiling of seeds showed no dramatic changes in the GM/non-GM maize varieties compared with the GM/GM and non-GM/non-GM maize varieties. The genetic background of the transgenic maize was found to have some influence on its proteomic and metabolomic profiles.

No adverse dietary effect of a cisgenic fire blight resistant apple line on the non-target arthropods Drosophila melanogaster and Folsomia candida

Genetic modification of apple cultivars through cisgenesis can introduce traits, such as disease resistance from wild relatives, quickly and without crossing. This approach was used to generate the cisgenic apple line C44.4.146, a 'Gala Galaxy' carrying the fire blight resistance gene FB_MR5. In contrast to traditionally bred apple cultivars, genetically modified (GM) plants need to undergo a regulatory risk assessment considering unintended effects before approval for commercial release. To determine potential unintended effects of C44.4.146, we assessed major leaf components and effects on the fitness of the decomposers Drosophila melanogaster (fruit fly) and Folsomia candida (collembolan), which were fed a diet amended with powdered apple leaf material. Leaf material of 'Gala Galaxy', several natural 'Gala' mutants, and the unrelated apple cultivar 'Ladina' were used for comparison. The genetic modification did not alter major leaf components and did not adversely affect survival, growth, or fecundity of the two decomposers. Consistent with previous studies with other GM crops, the differences between conventionally bred cultivars were greater than between the GM line and its non-GM wild type. These data provide a baseline for future risk assessments.

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.

Modified Gene Editing Systems: Diverse Bioengineering Tools and Crop Improvement

Gene-editing systems have emerged as bioengineering tools in recent years. Classical gene-editing systems include zinc-finger nucleases (ZFNs), transcription activator-like effector nucleases (TALENs), and clustered regularly interspaced short palindromic repeats (CRISPR) with CRISPR-associated protein 9 (Cas9) (CRISPR/Cas9), and these tools allow specific sequences to be targeted and edited. Various modified gene-editing systems have been established based on classical gene-editing systems. Base editors (BEs) can accurately carry out base substitution on target sequences, while prime editors (PEs) can replace or insert sequences. CRISPR systems targeting mitochondrial genomes and RNA have also been explored and established. Multiple gene-editing techniques based on CRISPR/Cas9 have been established and applied to genome engineering. Modified gene-editing systems also make transgene-free plants more readily available. In this review, we discuss the modifications made to gene-editing systems in recent years and summarize the capabilities, deficiencies, and applications of these modified gene-editing systems. Finally, we discuss the future developmental direction and challenges of modified gene-editing systems.

RNA Interference-Based Genetic Engineering Maize Resistant to Apolygus lucorum Does Not Manifest Unpredictable Unintended Effects Relative to Conventional Breeding: Short Interfering RNA, Transcriptome, and Metabolome Analysis

The use of omics techniques to analyze the differences between genetic engineering organisms and their parents can identify unintended effects and explore whether such unintended effects will have negative consequences. In order to evaluate whether genetic engineering will cause changes in crops beyond the changes introduced by conventional plant breeding, we compared the extent of transcriptome and metabolome modification in the leaves of three lines developed by RNA interference (RNAi)-based genetic engineering and three lines developed by conventional breeding. The results showed that both types of plant breeding methods can manifest changes at the short interfering RNA (siRNA), transcriptomic, and metabolic levels. Relative expression analysis of potential off-target gene revealed that there was no broad gene decline in the three RNAi-based genetic engineering lines. We found that the number of DEGs and DAMs between RNAi-based genetic engineering lines and the parental line was less than that between conventional breeding lines. These unique DEGs and DAMs between RNAi-based genetic engineering lines and the parental lines were not enriched in detrimental metabolic pathways. The results suggest that RNAi-based genetic engineering do not cause unintended effects beyond those found in conventional breeding in maize.

Cooperative herbivory between two important pests of rice

Normally, when different species of herbivorous arthropods feed on the same plant this leads to fitness-reducing competition. We found this to be different for two of Asia's most destructive rice pests, the brown planthopper and the rice striped stem borer. Both insects directly and indirectly benefit from jointly attacking the same host plant. Double infestation improved host plant quality, particularly for the stemborer because the planthopper fully suppresses caterpillar-induced production of proteinase inhibitors. It also reduced the risk of egg parasitism, due to diminished parasitoid attraction. Females of both pests have adapted their oviposition behaviour accordingly. Their strong preference for plants infested by the other species even overrides their avoidance of plants already attacked by conspecifics. This cooperation between herbivores is telling of adaptations resulting from the evolution of plant-insect interactions, and points out mechanistic vulnerabilities that can be targeted to control these major pests.

Metabolic Analysis Reveals Cry1C Gene Transformation Does Not Affect the Sensitivity of Rice to Rice Dwarf Virus

Metabolomics is beginning to be used for assessing unintended changes in genetically modified (GM) crops. To investigate whether Cry1C gene transformation would induce metabolic changes in rice plants, and whether the metabolic changes would pose potential risks when Cry1C rice plants are exposed to rice dwarf virus (RDV), the metabolic profiles of Cry1C rice T1C-19 and its non-Bt parental rice MH63 under RDV-free and RDV-infected status were analyzed using gas chromatography–mass spectrometry (GC-MS). Compared to MH63 rice, slice difference was detected in T1C-19 under RDV-free conditions (less than 3%), while much more metabolites showed significant response to RDV infection in T1C-19 (15.6%) and in MH63 (5.0%). Pathway analysis showed biosynthesis of lysine, valine, leucine, and isoleucine may be affected by RDV infection in T1C-19. No significant difference in the contents of free amino acids (AAs) was found between T1C-19 and MH63 rice, and the free AA contents of the two rice plants showed similar responses to RDV infection. Furthermore, no significant differences of the RDV infection rates between T1C-19 and MH63 were detected. Our results showed the Cry1C gene transformation did not affect the sensitivity of rice to RDV, indicating Cry1C rice would not aggravate the epidemic and dispersal of RDV.

Compositional equivalence assessment of insect-resistant genetically modified rice using multiple statistical analyses

The safety of transgenic Bt rice containing bacteria-derived mCry1Ac gene from Bacillus thuringiensis (Bt) was assessed by conducting field trials at two locations for two consecutive years in South Korea, using the near-isogenic line comparator rice cultivar ('Ilmi', non-Bt rice) and four commercial cultivars as references. Compositional analyses included measurement of proximates, minerals, amino acids, fatty acids, vitamins, and antinutrients. Significant differences between Bt rice and non-Bt rice were detected; however, all differences were within the reference range. The statistical analyses, including analysis of % variability, analysis of similarities (ANOISM), similarity percentage (SIMPER) analysis, and permutational multivariate analysis of variance (PERMANOVA) were performed to study factors contributing to compositional variability. The multivariate analyses revealed that environmental factors more influenced rice components' variability than by genetic factors. This approach was shown to be a powerful method to provide meaningful evaluations between Bt rice and its comparators. In this study, Bt rice was proved to be compositionally equivalent to conventional rice varieties through multiple statistical methods.

Integrated proteomics and metabolomics analysis of transgenic and gene-stacked maize line seeds

Unintended effects of genetically modified (GM) crops may pose safety issues. Omics techniques provide researchers with useful tools to assess such unintended effects. Proteomics and metabolomics analyses were performed for three GM maize varieties, 2A-7, CC-2, and 2A-7×CC-2 stacked transgenic maize, and the corresponding non-GM parent Zheng58.Proteomics revealed 120, 271 and 135 maize differentially expressed proteins (DEPs) in the 2A-7/Zheng58, CC-2/Zheng58 and 2A-7×CC-2/Zheng58 comparisons, respectively. Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis showed that most DEPs participated in metabolic pathways and the biosynthesis of secondary metabolite. Metabolomics revealed 179, 135 and 131 differentially accumulated metabolites (DAMs) in the 2A-7/Zheng58, CC-2/Zheng58 and 2A-7×CC-2/Zheng58 comparisons, respectively. Based on KEGG enrichment analysis, most DAMs are involved in the biosynthesis of secondary metabolite and metabolic pathways. According to integrated proteomics and metabolomics analysis, the introduction of exogenous EPSPS did not affect the expression levels of six other enzymes or the abundance of seven metabolites involved in the shikimic acid pathway in CC-2 and 2A-7×CC-2 seeds. Six co-DEPs annotated by integrated proteomics and metabolomics pathway analysis were further analyzed by qRT-PCR.This study successfully employed integrated proteomic and metabolomic technology to assess unintended changes in maize varieties. The results suggest that GM and gene stacking do not cause significantly unintended effects.

Transcriptomic and Metabolomic Responses of Rice Plants to Cnaphalocrocis medinalis Caterpillar Infestation

Simple Summary

The transcriptomic and metabolomic differences in rice leaves after infestation by the rice leaf folder Cnaphalocrocis medinalis were investigated for better understanding of the mechanisms of rice defenses against this species. The results suggest that C. medinalis infestation can induce rapid and precise defense responses involved in many primary and secondary metabolic processes in rice leaves, and the jasmonic acid (JA)-dependent signaling pathway plays vital roles in the response of rice plants to this pest species. These results provide comprehensive insights into the defense system of rice to the rice leaf folder and may facilitate the development of insect-resistant rice varieties by identifying molecular targets for selection.

Abstract

Interactions between plants and insect herbivores are important determinants of plant productivity in cultivated and natural agricultural fields. The rice leaf folder (Cnaphalocrocis medinalis) causes tremendous damage to rice production in Asian countries. However, little information is available about how rice plants defend themselves against this destructive pest at molecular and biochemical levels. Here, we observed the transcriptomic and metabolomic differences in rice leaves after 0, 1, 6, 12, and 24 h of being fed by C. medinalis using RNA sequencing and metabolome profiling. Transcriptional analyses showed that gene expression responds rapidly to leaf folder infestation, with the most significant transcriptional changes occurring within 6 h after the initiation of feeding. Metabolite abundance changed more slowly than gene expression. Gene Ontology and Kyoto Encyclopedia of Genes and Genomes enrichment analyses indicated that the rice transcriptional response to infestation involved genes encoding protein kinases, transcription factors, biosynthesis of secondary metabolites, photosynthesis, and phytohormone signaling. Moreover, the jasmonic acid-dependent signaling pathway triggered by leaf folder herbivory played a vital role in rice defense against this pest. Taken together, our results provide comprehensive insights into the defense system of rice to this species and may inform the development of insect-resistant rice varieties.