Development of Marker-Free Insect-Resistant Indica Rice by Agrobacterium tumefaciens-Mediated Co-transformation

Agrobacterium-mediated co-transformation is an efficient strategy to generate marker-free transgenic plants. In this study, the vectors pMF-2A^∗ containing a synthetic cry2A^∗ gene driven by maize ubiquitin promoter and pCAMBIA1301 harboring hygromycin phosphotransferase gene (hpt) were introduced into Minghui86 (Oryza sativa L. ssp. indica), an elite indica restorer line. Two independent transformants containing both the cry2A^∗ gene and hpt gene were regenerated. Several homozygous marker-free transgenic progenies were derived from family 2AH2, and three of them were selected for further insect bioassay in the laboratory and field. Insect-resistance assays revealed that all the three transgenic lines were highly resistant to striped stem borer (Chilo suppressalis), yellow stem borer (Tryporyza incertulas) and rice leaf folder (Cnaphalocrocis medinalis). The measurement of Cry2A protein concentration showed that Cry2A protein was stably expressed in leaves and stems of homozygous transgenic lines and their hybrids. The yields of the marker-free homozygous transgenic lines and their hybrids were not significantly different from those of their corresponding controls. Furthermore, the results of flanking sequence isolation showed that the T-DNA in line 8-30 was integrated into the intergenic region of chromosome 2 (between Os02g43680 and Os02g43690). These results indicate that the marker-free transgenic rice has the potential for commercial production.

Introduction

Rice is one of the most important staple food crops in the world. Insect damage can cause severe rice yield loss up to 10% every year (Chen et al., 2009). Rice leaf folder, striped stem borer (SSB), and yellow stem borer (YSB) are three major pests in rice production. Traditionally, chemical pesticides are used to control these pests, which cause a series of problems, such as environmental pollution, pesticide residue and increase of production cost. Bt proteins have been widely used as a bioinsecticide and the Bt genes have been introduced into rice, corn, cotton, soybean, and other crops to control pests, such as Lepidoptera pests. Transgenic crops such as corn and cotton that express Bt proteins for commercial application have brought great economic benefits (Clive, 2012).

Currently, Agrobacterium tumefaciens-mediated genetic transformation is widely applied to the development of transgenic crops, and antibiotic or herbicide-resistant selectable marker genes (SMGs) are used to select transformed cells in the transformation process. SMGs are redundant after transformation, and are not conducive to continuous transformation and safety. Thus, a series of marker-free strategies were developed, such as biolistic method (Tu et al., 2000; Shiva Prakash et al., 2009) and A. tumefaciens-mediated co-transformation, site-specific recombination (Li et al., 2007; Qiu et al., 2010; Sengupta et al., 2010; Yu et al., 2013; Woo et al., 2015b), multi-auto-transformation vector system (Endo et al., 2002; Khan et al., 2006; Scaramelli et al., 2009), transposon-mediated method (Gao et al., 2015), intrachromosomal recombination (Zubko et al., 2000), and PCR-based method (de Vetten et al., 2003; Woo et al., 2015a). There are three approaches for the co-transformation: one plasmid harboring two T-DNAs in one Agrobacterium strain (Breitler et al., 2004; Matheka et al., 2013); two T-DNAs located on two plasmids in the same Agrobacterium strain (Parkhi et al., 2005; Sripriya et al., 2008); two T-DNAs in separate Agrobacterium strains (Dutt et al., 2012). The approach of ‘two T-DNAs in separate Agrobacterium strains’ was adopted in this study, because the marker-free progeny has the target T-DNA sequence but does not include other redundant sequences.

Similar to the application of chemical pesticides, wide cultivation of Bt crops also can cause the insects to evolve resistance to Bt genes. The cry1A genes such as cry1Aa, cry1Ab, and cry1Ac are used to control Lepidoptera pest in transgenic rice at present (Nayak et al., 1997; Cheng et al., 1998; Ye et al., 2001; Khanna and Raina, 2002). It has been reported that Cry1Aa, Cry1Ab, and Cry1Ac have cross-resistance because of the common binding sites (Escriche et al., 1997). Thus, these proteins cannot be applied to develop Bt gene-stacking rice. Studies focused on transgenic insect-resistant rice with other types of Bt genes are rather limited so far, and there have been only a few studies about marker-free transgenic insect-resistant rice (Tu et al., 2000; Kumar et al., 2010). Several previous studies have shown that the Cry2A and Cry1A proteins have different binding sites in brush border membrane vesicles (BBMVs) and exhibit no cross-resistance (Lee et al., 1997; Karim and Dean, 2000; Alcantara et al., 2004; Gouffon et al., 2011). Therefore, cry2A gene can be applied to develop marker-free Bt transgenic rice. Then, the marker-free cry2A transgenic rice can be used as a resistance source for developing gene-pyramiding marker-free transgenic rice, which can delay the evolution of resistance of insects to Bt rice. In this study, the cry2A^∗ gene, which was synthesized by rice codon optimization with a 69.45% sequence homology with the original cry2Aa gene (Chen et al., 2005), was introduced into an elite indica restorer line Minghui86 (MH86) with the purpose of developing marker-free cry2A transgenic rice.

Materials and Methods

Constructs and Genetic Transformation

pCAMBIA1300 (provided by the Centre for the Application of Molecular Biology in International Agriculture, Australia) was digested by BstXI and XhoI to remove the CaMV35s promoter and hygromycin phosphotransferase (hpt) gene, and then the cry2A^∗ gene driven by ubiquitin promoter (Christensen et al., 1992) was cloned into the multiple cloning sites to obtain the transformation vector pMF-2A^∗. The T-DNA region is displayed in Supplementary Figure S1. The pMF-2A^∗ and pCAMBIA1301 were introduced into A. tumefaciens EHA105, respectively. The two strains were used to infect an elite indica restorer line MH86 together at the ratio of 3:1. The callus culture and genetic transformation procedures were based on the method described by Lin and Zhang (Lin and Zhang, 2005).

Detection of PCR and Southern Blotting

The cry2A^∗ and hpt genes were detected by PCR and Southern blotting. The primers for the amplification of cry2A^∗ were Cry2A-F (5′-CGTGTCAATGCTGACCTGAT-3′) and Cry2A-R (5′-GATGCCGGACAGGATGTAGT-3′). The primers for the amplification of hpt were hpt-F (5′-AGAATCTCGTGCTTTCAGCTTCGA-3′) and hpt-R (5′-TCAAGACCAATGCGGAGCATATAC-3′). The PCR mixture contained 25 ng of template DNA, 2 μl of 10× buffer, 1.5 μl of 2 mM dNTP, 0.2 μl of each primer (10 μM), 0.2 μl of rtaq. The PCR reaction for the amplification of cry2A^∗ was implemented at 94°C for 5 min; 32 cycles of 94°C for 30 s, 57°C for 30 s, 72°C for 30 s; 72°C for 7 min. The PCR reaction to amplify hpt was implemented at 94°C for 5 min; 32 cycles of 94°C for 30 s, 58°C for 30 s, 72°C for 30 s; 72°C for 7 min.

The total genomic DNA of the transformants was isolated by CTAB method for Southern blotting (Murray and Thompson, 1980). There are HindIII sites in the T-DNAs of pMF2A^∗ (Supplementary Figure S1) and pCAMBIA1301, and there is no HindIII site between the probe (cry2A^∗ or hpt) and LB of the T-DNA. Therefore, the DNA samples (8 μg) were digested by HindIII. The plasmids pMF-2A^∗ and pCAMBIA1301 were also digested by HindIII and 0.5 ng of digested plasmids was loaded as positive control, respectively. The digested DNA samples were separated by 1% TAE agarose gel and transferred to nylon membrane. Each sample was hybridized by cry2A^∗ probe (a PCR-amplified fragment of cry2A^∗ by Cry2A-F and Cry2A-R) and hpt probe (a PCR-amplified fragment of hpt by hpt-F and hpt-R). The hybridization procedures were carried out as described by Southern (Southern, 1975).

Selection of Marker-Free Homozygous Cry2A^∗ Lines with Single Copy

The 160 T₁ generation transgenic plants derived from the T₀ transgenic plant with single copy of both cry2A^∗ and hpt were planted. The T₁ transgenic lines containing cry2A^∗ and no hpt were selected by PCR analysis, and were harvested individually. Then, the harvested seeds of these T₁ transgenic lines were sown in seedling bed and cry2A^∗ gene was detected in the seedlings (T₂ generation) through PCR analysis. If the cry2A^∗ gene did not segregate in T₂ generation population derived from the same T₁ plant, the T₁ plant was considered as a cry2A^∗ homozygous line. Three homozygous lines were used for further assays and the presence of cry2A^∗ and absence of hpt were confirmed by southern blotting.

Measurement of Cry2A^∗ Protein Content

In 2016, the concentrations of Cry2A^∗ protein in leaves (at tillering, heading and filling stage), stems (at tillering, heading and filling stage) and endosperms (at filling stage) of three transgenic homozygous lines and their hybrids were measured using the enzyme-linked immune sorbent assay (ELISA) kit AP005 from EnviroLogix Inc. (Portland, ME, USA). These hybrids were obtained by crossing the homozygous lines with II-32A, an elite indica cytoplasmic male sterility (CMS) line, respectively. Fresh leaf samples of approximately 20 mg were collected and homogenized by grinding in 500 μl of extraction buffer, and fresh stem and endosperm samples of approximately 40 mg were collected and homogenized by grinding in 1 ml of extraction buffer. All of the sample extracts were diluted in appropriate proportions. To quantify the Cry2A protein content, the Cry2Aa Calibrators (1, 5, and 10 ppb) provided by the manufacturer were used to create a standard curve. The Cry2A protein immunoassay procedures followed the protocol of the manufacturer. The optical density values were measured with a microplate reader at 450 nm wavelength. The Cry2A^∗ protein concentration (ppb) of each diluted sample was obtained from the standard curve. The final content of Cry2A^∗ protein was calculated according to the following equation: Cry2A^∗ protein content (μg/g) = [concentration value from the standard curve (ppb) × dilution multiple × extraction buffer volume (ml)]/fresh weight of each sample (mg). To calculate the percentage of Cry2A^∗ protein, the Bradford assay (Bradford, 1976) was used to quantify the total protein concentration. One-way analysis of variance and the LSD test were used to analyze the significant difference in the mean values between different stages of transgenic lines.

Insect Assay in the Laboratory

The resistance of the transgenic lines against SSB and YSB was tested in the laboratory. The SSB were fed in our lab, and the eggs of yellow rice borers were collected in the field and hatched in laboratory. The resistance of the three transgenic homozygous lines (MH86 as the susceptible control) and their hybrids (II Youming86 produced by crossing MH86 with II-32A as the susceptible control) was tested. The tested rice stems were collected from field at elongation stage and dissected into 8 cm in length. One glass tube (12 cm in length and 2 cm in diameter), which contained three dissected stems, was set as a biological repeat. Fifteen first-instar larvae of SSB of YSB were placed into the tube and sealed with cotton and black gauze. Three biological repeats were performed for each tested rice line. All of the glass tubes were incubated at 28°C and 80% relative humidity. After 6 days the larvae mortality was investigated. One-way analysis of variance and the LSD test were used to compare the significant difference in the mean values between transgenic lines and corresponding controls.

Insect Assay in the Field

Insect resistance of the three transgenic homozygous lines in the field was evaluated at Huazhong Agricultural University experimental field in 2011, Wuhan, China. The transformation receptor MH86 was planted as susceptible control. The seeds of all tested materials were sown on rice seedling bed in June 6, 2011 and the seedlings were transplanted to paddy field in July 2. The seedlings of the four rice materials were distributed according to a randomized block design with three replications. Each replication consisted of 24 plants of each material in two rows, with a planting distance of 20 cm within a row and 26.7 cm between rows. The experimental plots were bordered by three rows of non-transgenic indica rice plants as protection. The experimental field was not sprayed with chemical pesticides during the entire growth period.

To evaluate the resistance of transgenic lines against stem borers in the field, each rice plant was infested with about 15 first-instar larvae of YSB at the booting stage. Rice dead hearts caused by stem borers was investigated after two weeks of artificial infestation. Meanwhile, leaf folds caused by natural infestation of rice leaf folder were counted. The insect resistance of transgenic hybrids and II Youming86 was evaluated in the field in 2013 using the similar method in 2011 except for that the planting distance was 17 cm within a row and 26.7 cm between rows. Because the eggs of YSB were not collected in this year, the insecticidal activity of the hybrids was evaluated under natural infestation in the experimental paddy. Non-repeated two-factor analysis of variance and LSD test were used to compare the significant difference in the mean values between transgenic lines and corresponding controls.

Evaluation of Agronomic Performance in the Field

The design of experiments for evaluating the agronomic performance of the transgenic homozygous lines and their hybrids were conducted as described above in the insect assay in the field in 2011 and 2013, respectively. Chemical insecticides were sprayed during the growth period to protect the crop. The plant height of 20 plants in the middle of each plot was measured in the field at maturity, and these plants were harvested to record the number of panicles per plant, panicle length, number of grains per panicle, 1000-grain weight, seed setting rate and yield per plant. Non-repeated two-factor analysis of variance and LSD test were also used to compare the significant difference in the mean values between transgenic lines and corresponding controls.

Isolation of Flanking Sequences

Thermal asymmetric interlaced polymerase chain reaction (TAIL PCR) was adopted for isolating the flanking sequences of the T-DNA integration site in the transgenic plant. Three nested specific primers (SP1, SP2, and SP3) near the T-DNA right border (RB) and degenerate primer (AD8) were designed (Supplementary Table S1). The PCR mixture and reaction programs were based on the method described by Liu and Chen (2007). The tertiary amplification product was cloned to pEasy-T3 vector and sequenced. The flanking sequence of T-DNA RB flanking region was obtained by analysis of the sequencing results, and then the sequence was used for blastn in the bioinformatic website¹ in order to determine the T-DNA integration site in rice genome. Meanwhile, the putative flanking sequence of T-DNA left border (LB) flanking region was obtained through querying in bioinformatic website² according to the blastn results. The genome primer 2AH2-tF based on the putative flanking sequence of LB and primer Tail-L1 in vector pMF-2A^∗ were designed (Supplementary Table S1), and the PCR product of the two primers was sequenced to confirm the flanking sequence of LB.

Results

Genetic Transformation and DNA Assay

Only two positive co-transformation plants, which were named as 2AH1 and 2AH2, respectively, were obtained in T₀ generation by co-transformation of MH86 callus. The copy number of transgenic plants was confirmed by southern blotting. The results showed that both cry2A^∗ and hpt genes were single-copy in the 2AH2 family (Supplementary Figure S2). Three marker-free homozygous lines containing cry2A^∗ (named as 7-61, 8-30, and 8-62) derived from 2AH2 were selected by PCR and southern blot analysis (Figure 1).

FIGURE 1

FIGURE 1. Southern blotting of T₂ generation marker-free transgenic homozygous lines. (A) cry2A^∗ probe, (B) hpt probe. M DNA Marker, N MH86 control, P positive plasmid control, 1 line 7-61, 2 line 8-30, 3 line 8-62.

Cry2A^∗ Concentrations of Transgenic Homozygous Lines and Hybrids

In 2016, the Cry2A protein contents and the percentage of Cry2A protein in soluble protein of transgenic homozygous lines were measured at tillering, heading and filling stage (Figures 2A-C; Supplementary Table S2). In leaves, the concentrations of Cry2A protein ranged from 5.01 μg g^-1 to 26.79 μg g^-1 fresh weight, and the percentage of Cry2A protein ranged from 0.016 to 0.099% in soluble protein. In stems, the Cry2A protein concentrations ranged from 0.70 μg g^-1 to 2.35 μg g^-1 fresh weight, and the percentage of Cry2A protein ranged from 0.018 to 0.039% in soluble protein. In endosperms, the contents of Cry2A protein ranged from 3.34 μg g^-1 to 4.09 μg g^-1 fresh weight, and the percentage of Cry2A protein ranged from 0.026 to 0.031% in soluble protein. These results indicated that the Cry2A protein was stably expressed in the three transgenic homozygous lines at each growth stage.

FIGURE 2

FIGURE 2. Cry2A protein concentrations in three homozygous lines and their hybrids. (A-C) homozygous lines, (D-F) hybrids. a, b, c: significant difference at P < 0.05. Error bars indicate SD (n = 3).

To evaluate the commercial potential of the three transgenic homozygous lines, three hybrids, namely H_7-61, H_8-30, and H_8-62, were produced from the crossing of II-32A with the transgenic homozygous lines, respectively. The Cry2A protein content and its proportion in soluble protein in the Bt hybrids were also measured in 2016. In the hybrids, the concentrations of Cry2A protein in leaves ranged from 1.79 μg g^-1 to 14.56 μg g^-1 fresh weight, and the proportion of Cry2A protein ranged from 0.007 to 0.049% in soluble protein. In stems, the Cry2A protein concentrations ranged from 0.27 μg g^-1 to 0.95 μg g^-1 fresh weight, and the proportion of Cry2A protein ranged from 0.008 to 0.017% in soluble protein. In endosperms, the contents of Cry2A protein ranged from 1.56 μg g^-1 to 1.92 μg g^-1 fresh weight, and the proportion of Cry2A protein ranged from 0.012 to 0.013% in soluble protein (Figures 2D-F and Supplementary Table S3). These results suggested that the Cry2A contents in the hybrids were lower than in the homozygous lines.

Insect-Resistance of Transgenic Rice in the Laboratory

In the laboratory assay of the insecticidal activity of the transgenic homozygous lines, the larvae mortality of SSB or YSB fed on non-transgenic control MH86 stems was approximately 22% after 6 days, and the surviving larvae showed normal development. In addition, the MH86 stems were severely damaged and plenty of frass was observed. The larvae mortality of SSB fed on transgenic rice stems was about 50%, and the surviving larvae showed developmental retardation: their body length was significantly shorter than that of survival from non-transgenic control MH86. Besides, there was a smaller amount of frass (Figures 3A,C). All YSB larvae fed on transgenic rice stem were killed after 6 days, and the stems suffered from slight damage and little frass (Figures 3A,D). These results revealed the high resistance of transgenic homozygous lines against YSB.

FIGURE 3

FIGURE 3. Insect-resistance assay of transgenic materials in the laboratory. Larvae mortality of SSB and YSB fed on three transgenic homozygous lines (A); larvae mortality of SSB fed on three transgenic hybrids (B); resistance of transgenic homozygous lines against first-instar SSB (C) and YSB (D); resistance of transgenic hybrids against first-instar SSB (E). a, b, c: significant difference at P < 0.05. Error bars indicate SD (n = 3).

In the insect-resistance assay of the hybrids, the stem cuttings of II Youming86 showed similar level of damage as the stems of MH86, and the mortality of SSB larvae was about 25%. Besides, the surviving larvae developed normally. However, the three transgenic hybrids could kill about 50% of the SSB larvae, and the development of the surviving larvae was retarded (Figures 3B,E). The results indicated that the Cry2A protein content of the transgenic hybrids was lower than that of transgenic homozygous lines, but the hybrid stems showed equivalent insecticidal activity with their parents.

Insect-Resistance of Transgenic Rice in the Field

In the insect-resistance assay of the transgenic homozygous lines in the field in 2011, the leaves of non-transgenic control MH86 showed serious folds caused by rice leaf folders. Averagely two leaves were recorded as being folded per tiller and the percentage of tillers suffering from damage per plant was up to 99%. However, the percentage of tillers damaged in the three transgenic homozygous lines ranged from 4.78 to 8.25% and the leaves were only slightly damaged (0.06-0.11 leaves per tiller). Very few white heads were observed in the transgenic homozygous lines in contrast to the white heads in non-transgenic control MH86 per plant (11.84%) (Table 1; Figure 4A).

TABLE 1

TABLE 1. Resistance of transgenic homozygous lines against rice leaf folder and stem borers in the field under artificial infestation combined with natural infestation in 2011.

FIGURE 4

FIGURE 4. Insect resistance performance of transgenic rice in the field. (A) Resistance of transgenic homozygous lines against rice leaf folder and rice stem borer. CK: MH86, TL: transgenic line. Red arrow indicates white head. (B) Resistance of transgenic hybrids against rice stem borer. CK: II Youming86, TL: transgenic hybrid. Red arrow indicates dead hearts.

The resistance of the transgenic hybrids against target insects in the field was evaluated under natural infestation in 2013. The damage caused by rice leaf folders was not severe in this year. 15.46% of plant tillers suffered from damage in susceptible control II Youming86, and the transgenic hybrids showed only slight damage (Table 2). Some dead hearts caused by rice stem borers were observed in II Youming86 (Figure 4B), while almost no damage was observed in the transgenic hybrids. These results revealed that the transgenic rice was highly resistant to rice leaf folder and rice stem borer in the field.

TABLE 2

TABLE 2. Resistance performance of transgenic hybrids in the field after natural infestation in 2013. Values are the means ± SD (n = 3).

Agronomic Performance of Transgenic Rice in the Field

Some agronomic traits of the three transgenic homozygous lines showed different degrees of variation. The panicle length and 1000-grain weight were significantly lower than those of non-transgenic control MH86 (P < 0.01), but the numbers of grains per panicle of line 7-61 and line 8-30 were significantly larger than that of MH86 (P < 0.05). All the transgenic homozygous lines showed no change in yield compared with wide-type MH86 (Table 3). The three transgenic hybrids showed little difference with non-transgenic control II Youming86 in agronomic performance, especially the hybrid of line 8-30 (Table 4).

TABLE 3

TABLE 3. Agronomic performance of transgenic homozygous lines in the field.

TABLE 4

TABLE 4. Agronomic performance of transgenic hybrids in the field.

Flanking Sequences of T-DNA Integration Site in Rice Genome

The sequence analysis results of TAIL-PCR and subsequent routine PCR revealed that one complete T-DNA copy containing cry2A* gene was inserted into rice genome chromosome 2 in line 8-30. The information of the genes nearby T-DNA integration site was obtained by querying the flanking sequence in the bioinformatics website³. The results indicated that the T-DNA LB and transcription initiation site of gene Os02g43680 had a distance of about 1.5 kb, and the T-DNA RB was 2.99 kb away from the transcription initiation site of gene Os02g43690 (Figure 5A); besides, the T-DNA insertion led to a 7-bp deletion of rice genome (ctcgcgc) (Figure 5B).

FIGURE 5

FIGURE 5. Information of T-DNA integration site in transgenic line 8-30. (A) T-DNA inserted into rice chromosome 2 between gene Os02g43680 and gene Os02g43690. (B) Nucleotide sequences of T-DNA integration site. The two flanking sequences of T-DNA are displayed in lowercase letter. The sequence in bracket is the deletion sequence of genome. The T-DNA sequence is shown in capital letter. Bold capital letters represent the truncated T-DNA left border (LB). T-DNA right border (RB) is missing in line 8-30.

Discussion

In this study, a synthetic cry2A* gene was introduced into an elite indica restorer MH86 via Agrobacterium-mediated co-transformation. Subsequently, three single-copy marker-free transgenic homozygous lines were selected by PCR analysis and Southern blotting. The three transgenic lines not only exhibited high resistance against rice leaf folder and rice stem borer (Figures 3 and 4), but also showed no significant difference in yield compared with their corresponding controls under the spraying of chemical insecticides during the entire growth period (Tables 3 and 4). In the laboratory insect assay, the larvae mortality of SSB was lower than that of YSB. This result was consistent with the previous study in our laboratory (Yang et al., 2011), which indicated that SSB is less susceptible to Cry2A protein than YSB. The difference in susceptibility is mainly due to the difference in Cry2A protein binding sites and binding affinities in the BBMVs of the insect midguts (Lee et al., 1997; Alcantara et al., 2004). It has been reported that the LC₅₀ of Cry2A protein in the artificial diet for SSB was significant higher than that for YSB (Lee et al., 1997), which also suggests that Cry2A protein is more toxic to YSB than to SSB.

As we know, the presence of SMGs in transgenic plants can cause a variety of problems. Herbicide resistance genes, which can endow the transgenic plants with herbicide tolerance, are often used as selectable markers in transformation. However, the majority of SMGs are useless after the generation of transgenic plants, and their presence will increase the metabolism burden. More importantly, the presence of SMGs in transgenic plants also brings about safety concerns, particularly in rice, which is directly consumed by humans. In the past two decades, a number of marker-free plant transformation strategies have been reported and commonly used, including the Agrobacterium-mediated co-transformation and site-specific recombination. In this study, the approach ‘two T-DNAs in separate Agrobacterium strains’ was adopted. Two Agrobacterium strains, respectively, harboring the vector with SMG (hpt gene) and the vector with insect-resistant gene cry2A* were used for genetic transformation, and marker-free transgenic rice was developed successfully, which might address the consumers’ concerns about food safety of transgenic rice.

With the large-scale cultivation of Bt transgenic crops, the insect populations that evolve resistance to Bt proteins have been found in the field. The evolution of insect resistance is mainly due to the mutation of cadherin, which binds with Bt protein in the midgut of insect. Gene pyramiding is an effective strategy to deal with insect resistance. The cry1A genes, including cry1Ab, cry1Ac, and the fusion gene cry1Ab/Ac, are most commonly used in transgenic rice so far, especially in marker-free transgenic rice. The proteins encoded by cry1A genes have common binding sites in the BBMVs of insect midgut, and they have cross-resistance between each other. Therefore, other types of Bt genes should be combined with cry1A genes to develop gene-stacking Bt transgenic rice. In this study, the marker-free transgenic rice expressing Cry2A protein, which showed different binding sites with Cry1A protein in the midgut of SSB or YSB, can be used to cross with marker-free Cry1A transgenic rice for developing marker-free transgenic rice with double Bt genes to delay the evolution of insect resistance. The marker-free cry1Ab/Ac transgenic indica rice TT51 developed previously by our laboratory has been issued with security certificate again by Chinese government in 2014, indicating the great commercialization potential of marker-free Bt transgenic rice. Thus, it can be speculated that the marker-free cry2A* transgenic rice developed in this study has the similar potential for commercialization.

The results of flanking sequence isolation revealed that T-DNA containing cry2A* was integrated into the intergenic region of chromosome 2 (Figure 5). Although the transgenic homozygous lines were significantly different from non-transgenic wide-type MH86 in several agronomic traits (Table 3), the transgenic hybrids showed little difference from the non-transgenic control II Youming86 in agronomic performance (Table 4). These results suggest that the variations of agronomic traits in transgenic homozygous lines might be caused by the somatic mutation in the tissue culture. However, most of these variations were recovered in the hybrids. Yang et al. (2011) reported the similar results when they recorded the agronomic traits of the hybrids derived from four single Bt-gene transgenic rice lines and 10 two Bt-gene transgenic rice lines crossing with an indica CMS Zhenshan97A, respectively. These results suggest that the consecutive backcross of the transgenic homozygous line 8-30 with the wide-type MH86 can eliminate the somatic variations (Bregitzer et al., 2008; Yang et al., 2012).

Author Contributions

FL, FZ, HC, and YL designed the experiments. FL performed the experiments, collected, and analyzed the data. FL and YL wrote the manuscript.

Conflict of Interest Statement

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Acknowledgments

This research was funded by the National Program of Transgenic Variety Development of China (2016ZX08001-001) and the National High Technology Research and Development Program of China (863 Program). And we thank Prof. Zuoxiong Liu (Huazhong Agricultural University) for language improvement.

Supplementary Material

The Supplementary Material for this article can be found online at: http://journal.frontiersin.org/article/10.3389/fpls.2016.01608/full#supplementary-material

Footnotes

1. ^ http://rice.plantbiology.msu.edu/
2. ^ http://rice.plantbiology.msu.edu/cgi-bin/gbrowse/rice/#search
3. ^ http://rice.plantbiology.msu.edu/cgi-bin/gbrowse/rice/\#search

References

Alcantara, E. P., Aguda, R. M., Curtiss, A., Dean, D. H., and Cohen, M. B. (2004). Bacillus thuringiensis δ-endotoxin binding to brush border membrane vesicles of rice stem borers. Arch. Insect Biochem. Physiol. 55, 169–177. doi: 10.1002/arch.10128

PubMed Abstract | CrossRef Full Text | Google Scholar

Bradford, M. M. (1976). A rapid and sensitive method for the quantitation of microgram quantities of protein utilizing the principle of protein-dye binding. Anal. Biochem. 72, 248–254. doi: 10.1016/0003-2697(76)90527-3

CrossRef Full Text | Google Scholar

Bregitzer, P., Dahleen, L. S., Neate, S., Schwarz, P., and Manoharan, M. (2008). A single backcross effectively eliminates agronomic and quality alterations caused by somaclonal variation in transgenic barley. Crop Sci. 48, 471–479. doi: 10.2135/cropsci2007.06.0370

CrossRef Full Text | Google Scholar

Breitler, J. C., Meynard, D., Van Boxtel, J., Royer, M., Bonnot, F., Cambillau, L., et al. (2004). A novel two T-DNA binary vector allows efficient generation of marker-free transgenic plants in three elite cultivars of rice (Oryza sativa L.). Transgenic Res. 13, 271–287. doi: 10.1023/B:TRAG.0000034626.22918.0a

PubMed Abstract | CrossRef Full Text | Google Scholar

Chen, H., Lin, Y., and Zhang, Q. (2009). Review and prospect of transgenic rice research. Chin. Sci. Bull. 54, 4049–4068. doi: 10.1007/s11434-009-0645-x

CrossRef Full Text | Google Scholar

Chen, H., Tang, W., Xu, C., Li, X., Lin, Y., and Zhang, Q. (2005). Transgenic indica rice plants harboring a synthetic cry2A* gene of Bacillus thuringiensis exhibit enhanced resistance against lepidopteran rice pests. Theor. Appl. Genet. 111, 1330–1337. doi: 10.1007/s00122-005-0062-8

PubMed Abstract | CrossRef Full Text | Google Scholar

Cheng, X., Sardana, R., Kaplan, H., and Altosaar, I. (1998). Agrobacterium-transformed rice plants expressing synthetic cryIA (b) and cryIA (c) genes are highly toxic to striped stem borer and yellow stem borer. Proc. Natl. Acad. Sci. U.S.A. 95, 2767–2772. doi: 10.1073/pnas.95.6.2767

PubMed Abstract | CrossRef Full Text | Google Scholar

Christensen, A. H., Sharrock, R. A., and Quail, P. H. (1992). Maize polyubiquitin genes: structure, thermal perturbation of expression and transcript splicing, and promoter activity following transfer to protoplasts by electroporation. Plant Mol. Biol. 18, 675–689. doi: 10.1007/BF00020010

PubMed Abstract | CrossRef Full Text | Google Scholar

Clive, J. (2012). Global Status of Commercialized Biotech/GM Crops: 2012. Ithaca, NY: ISAAA.

Google Scholar

de Vetten, N., Wolters, A. M., Raemakers, K., Van Der Meer, I., Ter Stege, R., Heeres, E., et al. (2003). A transformation method for obtaining marker-free plants of a cross-pollinating and vegetatively propagated crop. Nat. Biotechnol. 21, 439–442. doi: 10.1038/nbt801

PubMed Abstract | CrossRef Full Text | Google Scholar

Dutt, M., Li, Z. T., Dhekney, S. A., and Gray, D. J. (2012). Co-transformation of grapevine somatic embryos to produce transgenic plants free of marker genes. Methods Mol. Biol. 847, 201–213. doi: 10.1007/978-1-61779-558-9_17

PubMed Abstract | CrossRef Full Text | Google Scholar

Endo, S., Sugita, K., Sakai, M., Tanaka, H., and Ebinuma, H. (2002). Single-step transformation for generating marker-free transgenic rice using the ipt-type MAT vector system. Plant J. 30, 115–122. doi: 10.1046/j.1365-313X.2002.01272.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Escriche, B., Ferré, J., and Silva, F. J. (1997). Occurrence of a common binding site in Mamestra brassicae, Phthorimaea operculella, and Spodoptera exigua for the insecticidal crystal proteins CryIA from Bacillus thuringiensis. Insect Biochem. Mol. Biol. 27, 651–656.

PubMed Abstract | Google Scholar

Gao, X., Zhou, J., Li, J., Zou, X., Zhao, J., Li, Q., et al. (2015). Efficient generation of marker-free transgenic rice plants using an improved transposon-mediated transgene reintegration strategy. Plant Physiol. 167, 11–24. doi: 10.1104/pp.114.246173

PubMed Abstract | CrossRef Full Text | Google Scholar

Gouffon, C., Van Vliet, A., Van Rie, J., Jansens, S., and Jurat-Fuentes, J. (2011). Binding sites for Bacillus thuringiensis Cry2Ae toxin on heliothine brush border membrane vesicles are not shared with Cry1A, Cry1F, or Vip3A toxin. Appl. Environ. Microbiol. 77, 3182–3188. doi: 10.1128/AEM.02791-10

PubMed Abstract | CrossRef Full Text | Google Scholar

Karim, S., and Dean, D. (2000). Toxicity and receptor binding properties of Bacillus thuringiensisδ-endotoxins to the midgut brush border membrane vesicles of the rice leaf folders, Cnaphalocrocis medinalis and Marasmia patnalis. Curr. Microbiol. 41, 276–283. doi: 10.1007/s002840010134

PubMed Abstract | CrossRef Full Text | Google Scholar

Khan, R. S., Chin, D. P., Nakamura, I., and Mii, M. (2006). Production of marker-free transgenic Nierembergia caerulea using MAT vector system. Plant Cell Rep. 25, 914–919. doi: 10.1007/s00299-006-0125-6

PubMed Abstract | CrossRef Full Text | Google Scholar

Khanna, H. K., and Raina, S. K. (2002). Elite Indica transgenic rice plants expressing modified Cry1Ac endotoxin of Bacillus thuringiensis show enhanced resistance to yellow stem borer (Scirpophaga incertulas). Transgenic Res. 11, 411–423. doi: 10.1023/A:1016378606189

PubMed Abstract | CrossRef Full Text | Google Scholar

Kumar, S., Arul, L., and Talwar, D. (2010). Generation of marker-free Bt transgenic indica rice and evaluation of its yellow stem borer resistance. J. Appl. Genet. 51, 243–257. doi: 10.1007/BF03208854

PubMed Abstract | CrossRef Full Text | Google Scholar

Lee, M. K., Aguda, R., Cohen, M., Gould, F., and Dean, D. (1997). Determination of binding of Bacillus thuringiensis (delta)-endotoxin receptors to rice stem borer midguts. Appl. Environ. Microbiol. 63, 1453–1459.

PubMed Abstract | Google Scholar

Li, Z., Xing, A., Moon, B. P., Burgoyne, S. A., Guida, A. D., Liang, H., et al. (2007). A Cre/loxP-mediated self-activating gene excision system to produce marker gene free transgenic soybean plants. Plant Mol. Biol. 65, 329–341. doi: 10.1007/s11103-007-9223-2

PubMed Abstract | CrossRef Full Text | Google Scholar

Lin, Y. J., and Zhang, Q. (2005). Optimising the tissue culture conditions for high efficiency transformation of indica rice. Plant Cell Rep. 23, 540–547. doi: 10.1007/s00299-004-0843-6

PubMed Abstract | CrossRef Full Text | Google Scholar

Liu, Y.-G., and Chen, Y. (2007). High-efficiency thermal asymmetric interlaced PCR for amplification of unknown flanking sequences. BioTechniques 43, 649–656. doi: 10.2144/000112601

PubMed Abstract | CrossRef Full Text | Google Scholar

Matheka, J. M., Anami, S., Gethi, J., Omer, R. A., Alakonya, A., Machuka, J., et al. (2013). A new double right border binary vector for producing marker-free transgenic plants. BMC Res. Notes 6:448. doi: 10.1186/1756-0500-6-448

PubMed Abstract | CrossRef Full Text | Google Scholar

Murray, M. G., and Thompson, W. F. (1980). Rapid isolation of high molecular weight plant DNA. Nucleic Acids Res. 8, 4321–4325. doi: 10.1093/nar/8.19.4321

CrossRef Full Text | Google Scholar

Nayak, P., Basu, D., Das, S., Basu, A., Ghosh, D., Ramakrishnan, N. A., et al. (1997). Transgenic elite indica rice plants expressing CryIAc ∂-endotoxin of Bacillus thuringiensis are resistant against yellow stem borer (Scirpophaga incertulas). Proc. Natl. Acad. Sci. U.S.A. 94, 2111–2116. doi: 10.1073/pnas.94.6.2111

PubMed Abstract | CrossRef Full Text | Google Scholar

Parkhi, V., Rai, M., Tan, J., Oliva, N., Rehana, S., Bandyopadhyay, A., et al. (2005). Molecular characterization of marker-free transgenic lines of indica rice that accumulate carotenoids in seed endosperm. Mol. Genet. Genomics 274, 325–336. doi: 10.1007/s00438-005-0030-7

PubMed Abstract | CrossRef Full Text | Google Scholar

Qiu, C., Sangha, J. S., Song, F., Zhou, Z., Yin, A., Gu, K., et al. (2010). Production of marker-free transgenic rice expressing tissue-specific Bt gene. Plant Cell Rep. 29, 1097–1107. doi: 10.1007/s00299-010-0893-x

PubMed Abstract | CrossRef Full Text | Google Scholar

Scaramelli, L., Balestrazzi, A., Bonadei, M., Piano, E., Carbonera, D., and Confalonieri, M. (2009). Production of transgenic barrel medic (Medicago truncatula Gaernt.) using the ipt-type MAT vector system and impairment of Recombinase-mediated excision events. Plant Cell Rep. 28, 197–211. doi: 10.1007/s00299-008-0634-6

PubMed Abstract | CrossRef Full Text | Google Scholar

Sengupta, S., Chakraborti, D., Mondal, H. A., and Das, S. (2010). Selectable antibiotic resistance marker gene-free transgenic rice harbouring the garlic leaf lectin gene exhibits resistance to sap-sucking planthoppers. Plant Cell Rep. 29, 261–271. doi: 10.1007/s00299-010-0819-7

PubMed Abstract | CrossRef Full Text | Google Scholar

Shiva Prakash, N., Bhojaraja, R., Shivbachan, S. K., Hari Priya, G. G., Nagraj, T. K., Prasad, V., et al. (2009). Marker-free transgenic corn plant production through co-bombardment. Plant Cell Rep. 28, 1655–1668. doi: 10.1007/s00299-009-0765-4

PubMed Abstract | CrossRef Full Text | Google Scholar

Southern, E. M. (1975). Detection of specific sequences among DNA fragments separated by gel electrophoresis. J. Mol. Biol. 98, 503–517. doi: 10.1016/S0022-2836(75)80083-0

CrossRef Full Text | Google Scholar

Sripriya, R., Raghupathy, V., and Veluthambi, K. (2008). Generation of selectable marker-free sheath blight resistant transgenic rice plants by efficient co-transformation of a cointegrate vector T-DNA and a binary vector T-DNA in one Agrobacterium tumefaciens strain. Plant Cell Rep. 27, 1635–1644. doi: 10.1007/s00299-008-0586-x

PubMed Abstract | CrossRef Full Text | Google Scholar

Tu, J., Zhang, G., Datta, K., Xu, C., He, Y., Zhang, Q., et al. (2000). Field performance of transgenic elite commercial hybrid rice expressing Bacillus thuringiensis δ-endotoxin. Nat. Biotechnol. 18, 1101–1104. doi: 10.1038/80310

PubMed Abstract | CrossRef Full Text | Google Scholar

Woo, H.-J., Lee, S. B., Qin, Y., Lim, M.-H., Lee, J.-H., Shin, K.-S., et al. (2015a). Generation and molecular characterization of marker-free Bt transgenic rice plants by selectable marker-less transformation. Plant Biotechnol. Rep. 9, 351–360. doi: 10.1007/s11816-015-0371-6

CrossRef Full Text | Google Scholar

Woo, H. J., Qin, Y., Park, S. Y., Park, S. K., Cho, Y. G., Shin, K. S., et al. (2015b). Development of selectable marker-free transgenic rice plants with enhanced seed tocopherol content through FLP/FRT-mediated spontaneous auto-excision. PLoS ONE 10:e0132667. doi: 10.1371/journal.pone.0132667

PubMed Abstract | CrossRef Full Text | Google Scholar

Yang, Z., Chen, H., Tang, W., Hua, H., and Lin, Y. (2011). Development and characterisation of transgenic rice expressing two Bacillus thuringiensis genes. Pest Manag. Sci. 67, 414–422. doi: 10.1002/ps.2079

PubMed Abstract | CrossRef Full Text | Google Scholar

Yang, Z., Chen, H., Tang, W., and Lin, Y.-J. (2012). Effect of successive backcrossing on eliminating somaclonal variation caused by Agrobacterium-mediated transformation in rice. Acta Agron. Sin. 38, 814–819. doi: 10.3724/SP.J.1006.2012.00814

CrossRef Full Text | Google Scholar

Ye, G. Y., Shu, Q. Y., Yao, H. W., Cui, H. R., Cheng, X. Y., Hu, C., et al. (2001). Field evaluation of resistance of transgenic rice containing a synthetic cry1Ab gene from Bacillus thuringiensis Berliner to two stem borers. J. Econ. Entomol. 94, 271–276. doi: 10.1603/0022-0493-94.1.271

PubMed Abstract | CrossRef Full Text | Google Scholar

Yu, Y., Wang, Y., Tong, Q., Liu, X., Su, F., Quan, F., et al. (2013). A site-specific recombinase-based method to produce antibiotic selectable marker free transgenic cattle. PLoS ONE 8:e62457. doi: 10.1371/journal.pone.0062457

PubMed Abstract | CrossRef Full Text | Google Scholar

Zubko, E., Scutt, C., and Meyer, P. (2000). Intrachromosomal recombination between attP regions as a tool to remove selectable marker genes from tobacco transgenes. Nat. Biotechnol. 18, 442–445. doi: 10.1038/74515

PubMed Abstract | CrossRef Full Text | Google Scholar