Development of a Gene-Based Soybean-Origin Discrimination Method Using Allele-Specific Polymerase Chain Reaction.

A low soybean self-sufficiency rate in South Korea has caused a high import dependence and considerable price variation between domestic and foreign soybeans, causing the false labeling of foreign soybeans as domestic. Conventional soybean origin discrimination methods prevent a single-grain analysis and rely on the presence or absence of several compounds or concentration differences. This limits the origin discrimination of mixed samples, demonstrating the need for a method that analyzes individual grains. Therefore, we developed a method for origin discrimination using genetic analysis. The whole-genome sequencing data of the Williams 82 reference cultivar and 15 soybean varieties cultivated in South Korea were analyzed to identify the dense variation blocks (dVBs) with a high single-nucleotide polymorphism density. The PCR primers were prepared and validated for the insertion-deletion (InDel) sequences of the dVBs to discriminate each soybean variety. Our method effectively discriminated domestic and foreign soybean varieties, eliminating their false labeling.

Dry-milled flour rice 'Seolgaeng' harbors a mutated fructose-6-phosphate 2-kinase/fructose-2,6-bisphosphatase2.

'Seolgaeng', an opaque-endosperm rice (Oryza sativa) mutant, is used to prepare high-quality dry-milled rice flour. The mutation causing its opaque-endosperm phenotype was unknown. Map-based cloning identified a missense mutation in the gene FRUCTOSE-6-PHOSPHATE 2-KINASE/FRUCTOSE-2,6-BISPHOSPHATASE 2 (OsF2KP2) in Seolgaeng. Transfer DNA insertion and clustered regularly interspaced short palindromic repeat (CRISPR)/CRISPR-associated nuclease 9 (Cas9)-induced f2kp2 mutants exhibited opaque endosperm. Rice harbors another F2KP gene, OsF2KP1. CRISPR/Cas9-induced double mutants of OsF2KP1 and OsF2KP2 (f2kp-d) possessed more opaque endosperm compared to f2kp2 single mutants, whereas the endosperm of the f2kp1 single mutant was normal. Grain hardness and damaged starch content were significantly reduced in f2kp2 mutants compared to the wild type and f2kp1. Amylose content was lower than normal in f2kp2 mutants but not f2kp1. Grain hardness and amylose content were much lower in f2kp-d than in f2kp2. Starch polymerization analysis revealed altered amylopectin structure in f2kp2 and f2kp-d mutants. F2KP activity was lower in f2kp2 and much lower in the double mutants when compared to the wild types, but f2kp1 showed no significant difference. In coleoptiles, hypoxia induced OsF2KP2 expression but downregulated OsF2KP1. These results suggest that OsF2KP2 functions as the main F2KP isoform in endosperm experiencing hypoxia, but OsF2KP1 may partially compensate for the absence of OsF2KP2. We propose that F2KP has a crucial role in inorganic pyrophosphate-utilizing energy metabolism for starch biosynthesis in rice endosperm.

Over-Expression of Dehydroascorbate Reductase Improves Salt Tolerance, Environmental Adaptability and Productivity in Oryza sativa.

Abiotic stress induces reactive oxygen species (ROS) generation in plants, and high ROS levels can cause partial or severe oxidative damage to cellular components that regulate the redox status. Here, we developed salt-tolerant transgenic rice plants that overexpressed the dehydroascorbate reductase gene (OsDHAR1) under the control of a stress-inducible sweet potato promoter (SWPA2). OsDHAR1-expressing transgenic plants exhibited improved environmental adaptability compared to wild-type plants, owing to enhanced ascorbate levels, redox homeostasis, photosynthetic ability, and membrane stability through cross-activation of ascorbate-glutathione cycle enzymes under paddy-field conditions, which enhanced various agronomic traits, including root development, panicle number, spikelet number per panicle, and total grain yield. dhar2-knockdown plants were susceptible to salt stress, and owing to poor seed maturation, exhibited reduced biomass (root growth) and grain yield under paddy field conditions. Microarray revealed that transgenic plants highly expressed genes associated with cell growth, plant growth, leaf senescence, root development, ROS and heavy metal detoxification systems, lipid metabolism, isoflavone and ascorbate recycling, and photosynthesis. We identified the genetic source of functional genomics-based molecular breeding in crop plants and provided new insights into the physiological processes underlying environmental adaptability, which will enable improvement of stress tolerance and crop species productivity in response to climate change.

Hd1 Allele Types and Their Associations with Major Agronomic Traits in Korean Rice Cultivars.

Optimizing flowering time in crop plants is critical for maximizing yield and quality under target environments. While there is a wide range of heading date variation in Korean rice cultivars, the underlying gene mechanisms are unclear. Here, we sequenced the protein coding regions of Hd1, the major rice heading date gene, from 293 Korean rice cultivars and investigated the associations between Hd1 allele types and major agronomic traits under four different environments. There were four functional Hd1 and five nonfunctional hd1 alleles distributed among the 293 Korean rice cultivars. The effects of the Hd1 allele types were highly significant for days to heading in all four environments, explaining 51.4-65.8% of the phenotypic variation. On average, cultivars carrying nonfunctional hd1 headed 13.7 days earlier than those carrying functional Hd1. While the Hd1 allele types exhibited highly significant effects on culm length and protein content under all four environments, the differences between cultivars carrying Hd1 and hd1 were minimal. The effects of the Hd1 allele types on amylose content were significant in only one of the four environments. Our results provide useful information for fine-tuning rice heading dates by utilizing different Hd1 alleles in rice breeding programs.

Improved stress tolerance and productivity in transgenic rice plants constitutively expressing the Oryza sativa glutathione synthetase OsGS under paddy field conditions.

Reactive oxygen species, which increase under various environmental stresses, have deleterious effects on plants. An important antioxidant, glutathione, is used to detoxify reactive oxygen species in plant cells and is mainly produced by two enzymes: gamma-glutamylcysteine synthetase (γ-ECS) and glutathione synthetase (GS). To evaluate the functional roles of the glutathione synthetase gene (OsGS) in rice, we generated four independent transgenic rice plants (TG1-TG4) that overexpressed OsGS under the control of the constitutively expressed OsCc1 promoter. When grown under natural paddy field conditions, the TG rice plants exhibited greater growth development, higher chlorophyll content, and higher GSH/GSSH ratios than control wild-type (WT) rice plants. Subsequently, the TG rice plants enhanced redox homeostasis by preventing hydroperoxide-mediated membrane damage, which improved their adaptation to environmental stresses. As a result, TG rice plants improved rice grain yield and total biomass following increases in panicle number and number of spikelets per panicle, despite differences in climate during the cultivation periods of 2014 and 2015. Overall, our results indicate that OsGS overexpression improved redox homeostasis by enhancing the glutathione pool, which resulted in greater tolerance to environmental stresses in the paddy fields.

Barcode System for Genetic Identification of Soybean [Glycine max (L.) Merrill] Cultivars Using InDel Markers Specific to Dense Variation Blocks.

For genetic identification of soybean [Glycine max (L.) Merrill] cultivars, insertions/deletions (InDel) markers have been preferred currently because they are easy to use, co-dominant and relatively abundant. Despite their biological importance, the investigation of InDels with proven quality and reproducibility has been limited. In this study, we described soybean barcode system approach based on InDel makers, each of which is specific to a dense variation block (dVB) with non-random recombination due to many variations. Firstly, 2,274 VBs were mined by analyzing whole genome data in six soybean cultivars (Backun, Sinpaldal 2, Shingi, Daepoong, Hwangkeum, and Williams 82) for transferability to dVB-specific InDel markers. Secondly, 73,327 putative InDels in the dVB regions were identified for the development of soybean barcode system. Among them, 202 dVB-specific InDels from all soybean cultivars were selected by gel electrophoresis, which were converted as 2D barcode types according to comparing amplicon polymorphisms in the five cultivars to the reference cultivar. Finally, the polymorphism of the markers were assessed in 147 soybean cultivars, and the soybean barcode system that allows a clear distinction among soybean cultivars is also detailed. In addition, the changing of the dVBs in a chromosomal level can be quickly identified due to investigation of the reshuffling pattern of the soybean cultivars with 27 maker sets. Especially, a backcross-inbred offspring, "Singang" and a recurrent parent, "Sowon" were identified by using the 27 InDel markers. These results indicate that the soybean barcode system enables not only the minimal use of molecular markers but also comparing the data from different sources due to no need of exploiting allele binning in new varieties.

Map-based Cloning and Characterization of the BPH18 Gene from Wild Rice Conferring Resistance to Brown Planthopper (BPH) Insect Pest.

Brown planthopper (BPH) is a phloem sap-sucking insect pest of rice which causes severe yield loss. We cloned the BPH18 gene from the BPH-resistant introgression line derived from the wild rice species Oryza australiensis. Map-based cloning and complementation test revealed that the BPH18 encodes CC-NBS-NBS-LRR protein. BPH18 has two NBS domains, unlike the typical NBS-LRR proteins. The BPH18 promoter::GUS transgenic plants exhibited strong GUS expression in the vascular bundles of the leaf sheath, especially in phloem cells where the BPH attacks. The BPH18 proteins were widely localized to the endo-membranes in a cell, including the endoplasmic reticulum, Golgi apparatus, trans-Golgi network, and prevacuolar compartments, suggesting that BPH18 may recognize the BPH invasion at endo-membranes in phloem cells. Whole genome sequencing of the near-isogenic lines (NILs), NIL-BPH18 and NIL-BPH26, revealed that BPH18 located at the same locus of BPH26. However, these two genes have remarkable sequence differences and the independent NILs showed differential BPH resistance with different expression patterns of plant defense-related genes, indicating that BPH18 and BPH26 are functionally different alleles. These findings would facilitate elucidation of the molecular mechanism of BPH resistance and the identified novel alleles to fast track breeding BPH resistant rice cultivars.

Comparative evaluation of nutritional compositions between transgenic rice harboring the CaMsrB2 gene and the conventional counterpart.

As a part of a safety assessment of new transgenic crops, compositional equivalence studies between transgenic crops with non-transgenic comparators are almost universally required. This study was conducted to compare nutritional profiles of proximate composition, and fatty acid, amino acid, mineral, and vitamin contents, and anti-nutrients, between transgenic drought-tolerant Agb0103 rice harboring the pepper methionine sulfoxide reductase B2 gene CaMsrB2 and the parental rice cultivar, 'Ilmi' as a non-transgenic control. Both transgenic and non-transgenic rice were grown and harvested in 2 different locations. Proximate compositions of moisture, starch, protein, lipid, and ash content of Agb0103 rice were similar to parental non-transgenic rice. There were no differences between transgenic and non-transgenic rice with respect to the whole nutritional composition, except for minor locality differences for a few nutritional components. Agb0103 rice with improved resistance to drought is nutritionally equivalent to the parental rice cultivar.

A correlation between tocopherol content and antioxidant activity in seeds and germinating seeds of soybean cultivars.

BACKGROUND: Tocopherols are crucial lipid-soluble antioxidants and essential nutrients. There is increasing interest in the biofortification of crops with vitamin E for reducing micronutrient malnutrition. However, relatively little is known about the development of soybean cultivars with high levels of tocopherol through combined breeding. RESULT: Tocopherol contents of seeds and germinating seeds of 28 Korean soybean cultivars were analyzed and evaluated for health-promoting activities. Total tocopherol concentrations ranged from 203.9 to 503.1 µg g⁻¹ in seeds and from 20.1 to 230.1 µg g⁻¹ in germinating seeds. The traditional landraces of HaNagari (HN, 503.1 µg g⁻¹), Orialtae (OL, 486.6 µg g⁻¹), SuMoktae (SM, 476.5 µg g⁻¹) and SoRitae (SR, 475.5 µg g⁻¹) showed high levels of tocopherol content. The contents of the four isomers of tocopherol in seeds and germinating seeds were correlated with lipid peroxidation. The γ- and δ-tocopherol contents in seeds were related to 1,1-diphenyl-2-picrylhydrazyl free radical scavenging activity (0.434; P < 0.01 and 0.373; P < 0.05). CONCLUSION: Total tocopherol content was higher in soybean landraces as compared with modern cultivars developed by cross-breeding. These results suggest that soybean breeding is necessary to increase tocopherol levels.

Variation block-based genomics method for crop plants.

BACKGROUND: In contrast with wild species, cultivated crop genomes consist of reshuffled recombination blocks, which occurred by crossing and selection processes. Accordingly, recombination block-based genomics analysis can be an effective approach for the screening of target loci for agricultural traits. RESULTS: We propose the variation block method, which is a three-step process for recombination block detection and comparison. The first step is to detect variations by comparing the short-read DNA sequences of the cultivar to the reference genome of the target crop. Next, sequence blocks with variation patterns are examined and defined. The boundaries between the variation-containing sequence blocks are regarded as recombination sites. All the assumed recombination sites in the cultivar set are used to split the genomes, and the resulting sequence regions are termed variation blocks. Finally, the genomes are compared using the variation blocks. The variation block method identified recurring recombination blocks accurately and successfully represented block-level diversities in the publicly available genomes of 31 soybean and 23 rice accessions. The practicality of this approach was demonstrated by the identification of a putative locus determining soybean hilum color. CONCLUSIONS: We suggest that the variation block method is an efficient genomics method for the recombination block-level comparison of crop genomes. We expect that this method will facilitate the development of crop genomics by bringing genomics technologies to the field of crop breeding.

A pepper MSRB2 gene confers drought tolerance in rice through the protection of chloroplast-targeted genes.

BACKGROUND: The perturbation of the steady state of reactive oxygen species (ROS) due to biotic and abiotic stresses in a plant could lead to protein denaturation through the modification of amino acid residues, including the oxidation of methionine residues. Methionine sulfoxide reductases (MSRs) catalyze the reduction of methionine sulfoxide back to the methionine residue. To assess the role of this enzyme, we generated transgenic rice using a pepper CaMSRB2 gene under the control of the rice Rab21 (responsive to ABA protein 21) promoter with/without a selection marker, the bar gene. RESULTS: A drought resistance test on transgenic plants showed that CaMSRB2 confers drought tolerance to rice, as evidenced by less oxidative stress symptoms and a strengthened PSII quantum yield under stress conditions, and increased survival rate and chlorophyll index after the re-watering. The results from immunoblotting using a methionine sulfoxide antibody and nano-LC-MS/MS spectrometry suggest that porphobilinogen deaminase (PBGD), which is involved in chlorophyll synthesis, is a putative target of CaMSRB2. The oxidized methionine content of PBGD expressed in E. coli increased in the presence of H2O2, and the Met-95 and Met-227 residues of PBGD were reduced by CaMSRB2 in the presence of dithiothreitol (DTT). An expression profiling analysis of the overexpression lines also suggested that photosystems are less severely affected by drought stress. CONCLUSIONS: Our results indicate that CaMSRB2 might play an important functional role in chloroplasts for conferring drought stress tolerance in rice.

Co-expression of monodehydroascorbate reductase and dehydroascorbate reductase from Brassica rapa effectively confers tolerance to freezing-induced oxidative stress.

Plants are exposed to various environmental stresses and have therefore developed antioxidant enzymes and molecules to protect their cellular components against toxicity derived from reactive oxygen species (ROS). Ascorbate is a very important antioxidant molecule in plants, and monodehydroascorbate reductase (MDHAR; EC 1.6.5.4) and dehydroascorbate reductase (DHAR; EC 1.8.5.1) are essential to regeneration of ascorbate for maintenance of ROS scavenging ability. The MDHAR and DHAR genes from Brassica rapa were cloned, transgenic plants overexpressing either BrMDHAR and BrDHAR were established, and then, each transgenic plant was hybridized to examine the effects of co-expression of both genes conferring tolerance to freezing. Transgenic plants co-overexpressing BrMDHAR and BrDHAR showed activated expression of relative antioxidant enzymes, and enhanced levels of glutathione and phenolics under freezing condition. Then, these alteration caused by co-expression led to alleviated redox status and lipid peroxidation and consequently conferred improved tolerance against severe freezing stress compared to transgenic plants overexpressing single gene. The results of this study suggested that although each expression of BrMDHAR or BrDHAR was available to according tolerance to freezing, the simultaneous expression of two genes generated synergistic effects conferring improved tolerance more effectively even severe freezing.

Functional implication of ß-carotene hydroxylases in soybean nodulation.

Legume-Rhizobium spp. symbiosis requires signaling between the symbiotic partners and differential expression of plant genes during nodule development. Previously, we cloned a gene encoding a putative ß-carotene hydroxylase (GmBCH1) from soybean (Glycine max) whose expression increased during nodulation with Bradyrhizobium japonicum. In this work, we extended our study to three GmBCHs to examine their possible role(s) in nodule development, as they were additionally identified as nodule specific, along with the completion of the soybean genome. In situ hybridization revealed the expression of three GmBCHs (GmBCH1, GmBCH2, and GmBCH3) in the infected cells of root nodules, and their enzymatic activities were confirmed by functional assays in Escherichia coli. Localization of GmBCHs by transfecting Arabidopsis (Arabidopsis thaliana) protoplasts with green fluorescent protein fusions and by electron microscopic immunogold detection in soybean nodules indicated that GmBCH2 and GmBCH3 were present in plastids, while GmBCH1 appeared to be cytosolic. RNA interference of the GmBCHs severely impaired nitrogen fixation as well as nodule development. Surprisingly, we failed to detect zeaxanthin, a product of GmBCH, or any other carotenoids in nodules. Therefore, we examined the possibility that most of the carotenoids in nodules are converted or cleaved to other compounds. We detected the expression of some carotenoid cleavage dioxygenases (GmCCDs) in wild-type nodules and also a reduced amount of zeaxanthin in GmCCD8-expressing E. coli, suggesting cleavage of the carotenoid. In view of these findings, we propose that carotenoids such as zeaxanthin synthesized in root nodules are cleaved by GmCCDs, and we discuss the possible roles of the carotenoid cleavage products in nodulation.

Homologous expression of cytosolic dehydroascorbate reductase increases grain yield and biomass under paddy field conditions in transgenic rice (Oryza sativa L. japonica).

Dehydroascorbate reductase (DHAR, EC 1.8.5.1) maintains redox pools of ascorbate (AsA) by recycling oxidized AsA to reduced AsA. To investigate whether DHAR affects rice yield under normal environmental conditions, cDNA-encoding DHAR (OsDHAR1) was isolated from rice and used to develop OsDHAR1-overexpressing transgenic rice plants, under the regulation of a maize ubiquitin promoter. Incorporation and expression of the transgene in transgenic rice plants was confirmed by genomic polymerase chain reaction (PCR), semi-quantitative reverse transcription PCR (RT-PCR), western blot, and enzyme activity. The expression levels were at least twofold higher in transgenic (TG) rice plants than in control wild-type (WT) rice plants. In addition, OsDHAR1-overexpression in seven-independent homologous transgenic plants, as compared to WT plants, increased photosynthetic capacity and antioxidant enzyme activities under paddy field conditions, which led to an improved AsA pool and redox homeostasis. Furthermore, OsDHAR1 overexpression significantly improved grain yield and biomass due to the increase of culm and root weights and to enhance panicle and spikelet numbers in the same seven independent TG rice plants during the farming season (2010 and 2011) in South Korea. The OsDHAR protein contained the redox-active site (Cys20), as well as the conserved GSH-binding region, GSH-binding motif, glutathione-S-transferase (GST) N-terminal domain, C-terminal domain interface, and GST C-terminal domain. Therefore, our results indicate that OsDHAR1 overexpression, capable of functioning in AsA recycling, and protein folding increases environmental adaptation to paddy field conditions by the improving AsA pool and redox homeostasis, which enhances rice grain yield and biomass.

Homologous expression of γ-glutamylcysteine synthetase increases grain yield and tolerance of transgenic rice plants to environmental stresses.

Various environmental stresses induce reactive oxygen species (ROS), causing deleterious effects on plant cells. Glutathione (GSH), a critical antioxidant, is used to combat ROS. GSH is produced by γ-glutamylcysteine synthetase (γ-ECS) and glutathione synthetase (GS). To evaluate the functional roles of the Oryza sativa L. Japonica cv. Ilmi ECS (OsECS) gene, we generated transgenic rice plants overexpressing OsECS under the control of an inducible promoter (Rab21). When grown under saline conditions (100mM) for 4 weeks, 2-independent transgenic (TGR1 and TGR2) rice plants remained bright green in comparison to control wild-type (WT) rice plants. TGR1 and TGR2 rice plants also showed a higher GSH/GSSG ratio than did WT rice plants in the presence of 100mM NaCl, which led to enhanced redox homeostasis. TGR1 and TGR2 rice plants also showed lower ion leakage and higher chlorophyll-fluorescence when exposed to 10µM methyl viologen (MV). Furthermore, the TGR1 and TGR2 rice seeds had approximately 1.5-fold higher germination rates in the presence of 200mM salt. Under paddy field conditions, OsECS-overexpression in transgenic rice plants increased rice grain yield (TGW) and improved biomass. Overall, our results show that OsECS overexpression in transgenic rice increases tolerance and germination rate in the presence of abiotic stress by improving redox homeostasis via an enhanced GSH pool. Our findings suggest that increases in grain yield by OsECS overexpression could improve crop yields under natural environmental conditions.

FSTVAL: a new web tool to validate bulk flanking sequence tags.

BACKGROUND: Information about a transgene locus is one of the major concerns in transgenic research because expression of the transgene or a gene interrupted by the integration event could be affected. Thus, the flanking sequences obtained from transgenic plants need to be analyzed in terms of genomic context, such as genic and intergenic regions. This process may consist of several steps: 1) elimination of a vector sequence from the flanking sequence, 2) finding the locations in the target genome, and 3) statistics of the integration sites. These steps could be automated for flanking sequences from several dozens of transgenic plants generated in an ordinary targeted gene expression strategy. It would be indispensable in a genome-wide mutagenesis screen using T-DNA or transposons because these projects often generate several thousands of transgenic lines and just as many loci of the transgene among the transgenic plants. RESULTS: We present an open access web tool, flanking sequence tags validator (FSTVAL), to manage bulk flanking sequence tags (FSTs). FSTVAL automatically evaluates the FSTs and finds the best mapping positions of the FST against a known genome sequence. The statistics, in terms of genic and intergenic regions, are presented as a table, a distribution map, and a frequency graph along the chromosomes. Currently, 17 plant genome sequences, including Arabidopsis thaliana, Oryza sativa, and Glycine max, are available as reference genomes. We evaluated the utility and accuracy of the tool with 5,144 rice FSTs. The whole process, from uploading the sequences to generating tables of insertions, required a few minutes, with less than 4 clicks in the web environment. CONCLUSIONS: Run for 1 year and tested over 1,000 times, we have confirmed FSTVAL efficiently handles bulk FSTs. FSTVAL is freely available without login at http://bioinfo.mju.ac.kr/fstval/.

MicroRNA-34a modulates cytoskeletal dynamics through regulating RhoA/Rac1 cross-talk in chondroblasts.

MicroRNAs (miRNAs) have been implicated in various cellular processes, such as cell fate determination, cell death, and tumorigenesis. In the present study, we investigated the role of miRNA-34a (miR-34a) in the reorganization of the actin cytoskeleton, which is essential for chondrocyte differentiation. miRNA arrays to identify genes that appeared to be up-regulated or down-regulated during chondrogenesis were applied with chondrogenic progenitors treated with JNK inhibitor. PNA-based antisense oligonucleotides and miRNA precursor were used for investigation of the functional roles of miR-34a. We found that, in chick chondroprogenitors treated with JNK inhibitor, which suppresses chondrogenic differentiation, the expression levels of miR-34a and RhoA1 are up-regulated through modulation of Rac1 expression. Blockade of miR-34a via the use of PNA-based antisense oligonucleotides was associated with decreased protein expression of RhoA (a known modulator of stress fiber expression), down-regulation of stress fibers, up-regulation of Rac1, and recovery of protein level of type II collagen. miR-34a regulates RhoA/Rac1 cross-talk and negatively modulates reorganization of the actin cytoskeleton, which is one of the essential processes for establishing chondrocyte-specific morphology.

Suppression of NS3 and MP is important for the stable inheritance of RNAi-mediated rice stripe virus (RSV) resistance obtained by targeting the fully complementary RSV-CP gene.

Rice stripe virus (RSV) is a viral disease that seriously impacts rice production in East Asia, most notably in Korea, China, and Japan. Highly RSV-resistant transgenic japonica rice plants were generated using a dsRNAi construct designed to silence the entire sequence region of the RSV-CP gene. Transgenic rice plants were inoculated with a population of viruliferous insects, small brown planthoppers (SBPH), and their resistance was evaluated using ELISA and an infection rate assay. A correlation between the expression of the RSV-CP homologous small RNAs and the RSV resistance of the transgenic rice lines was discovered. These plants were also analyzed by comparing the expression pattern of invading viral genes, small RNA production and the stable transmission of the RSV resistance trait to the T3 generation. Furthermore, the agronomic trait was stably transmitted to the T4 generation of transgenic plants.