Fast-tracking development of homozygous transgenic cereal lines using a simple and highly flexible real-time PCR assay
© Mieog et al.; licensee BioMed Central Ltd. 2013
Received: 30 October 2012
Accepted: 23 April 2013
Published: 30 April 2013
A crucial step in the evaluation of newly produced transgenic plants is the selection of homozygous plants. Here we describe an efficient and highly flexible real-time PCR-based method for the development of homozygous lines in plant models with complex (multiple) genomes and/or relatively long generation times (>3 months) using direct copy number determinations.
An existing DNA extraction method was converted into a high-throughput plant leaf DNA extraction procedure yielding DNA suitable for real-time PCR analyses. Highly specific and efficient primer pairs were developed for a bread wheat reference gene (Epsilon Cyclase) and for standard sequence elements in the gene cassette routinely used for cereal transformations (an intron bridge and the Nopaline Synthase terminator). The real-time PCR assay reliably distinguished wheat plants with a single copy of the transgene from individuals with multiple copies or those lacking the transgene. To obtain homozygous lines carrying a unique insertion event as efficiently as possible, T0 plants (plants raised from transformed callus) with a single copy of the transgene were selected and their progeny screened for homozygous plants. Finally, the assay was adapted to work on rice.
The ability to quickly, easily and accurately quantify the construct copy numbers, as provided by the real-time PCR assay, greatly improved the efficiency and reliability of the selection of homozygous transgenic plants in our case study. We were able to select homozygous plants in early generations, avoiding time-consuming methods such as large scale analysis of segregation patterns of descendants and/or Southern blotting. Additionally, the ability to specifically develop homozygous lines carrying a unique insertion event could be important in avoiding gene silencing due to co-suppression, and if needed assist in the selection of lines suitable for future deregulation. The same primer pairs can be used to quantify many different wheat transgenic events because the construct-specific primer pairs are targeted to standard sequence elements of the cereal gene cassettes, making the method widely applicable in wheat GM research. Moreover, because all procedures described here are standardized, the method may easily be adapted to vectors lacking the target regions used here and/or to other plant models.
Genetic modification, the insertion of a DNA construct into a host genome, is an emerging powerful tool for increasing productivity and/or product quality of crop plants . However, the transformation process produces hemizygous plants (transgene is inserted without allelic counterpart) and the construct can be lost in subsequent generations through Mendelian segregation. Therefore, most downstream applications of genetic transformation require the transgene to be homozygous in the host genome. For models such as Arabidopsis thaliana with simple screening options such as marker selection , short generation times (<1 month), and single-to-low copy number transformation methods (through Agrobacterium tumefaciens), selection of homozygous plants can be achieved within respectable timelines through Mendelian segregation studies over few generations . However, this approach is very slow and cumbersome for agricultural models such as wheat, barley, rice and corn.
Most wheat transformations to date have been achieved using particle bombardment  in which, in many cases, multiple copies of the construct are inserted, either randomly throughout genomes or within a single insertion site (linked copies). Multiple unlinked copies cause complex segregation patterns, potentially requiring screening over many generations to establish the zygosity of a line. Linked copies will segregate as a single copy and will therefore not be noticed in Mendelian segregation studies. However, multiple inserted copies (either linked or unlinked) have been implicated in issues with transgene silencing  and are not suitable for lines being considered for deregulation, indicating a need to identify linked copies. Recently, progress has been made on the efficiency of A. tumefaciens mediated transformation in cereals including wheat (reviewed by Sood et al. ), which may limit the number of copies inserted into the genome and therefore reduce (but not avoid) some of these complications. Nonetheless, the multi-copy issues together with the long generation times of most agricultural models (>3 months) highlight the need for a more efficient alternative to Mendelian segregation studies, such as direct copy number determination.
The two most popular methods of directly quantifying copy numbers of a DNA fragment (gene, construct, etc.) are Southern blot and real-time PCR (reviewed by Bubner and Baldwin ). Southern blotting is the original method and is commonly used in order to determine copy numbers of DNA fragments. It involves blotting of enzymatically digested genomic DNA followed by hybridisation of specific labelled DNA probe corresponding to the sequence of interest . As powerful and reliable as the technique can be, this somewhat cumbersome technique displays some important limitations. Southern blot may fail to accurately assess the exact number of copies due to DNA alterations and/or loss of restriction site, or to distinguish copies that are situated closely together (in particular linked copies) . In addition, Southern blots will not usually give information on zygosity, which is needed to select a homozygous plant in the T1 generation. Using this technique, information on zygosity would have to be gathered in the T2 generation by analyzing the segregation ratios of the trait of interest or additional Southern blots. Furthermore, Southern blotting requires each construct to have a specific and labelled probe, making it less suitable to screen sample collections containing multiple constructs.
Real-time PCR is based on the detection of fluorescence produced during the amplification process of the PCR. This fluorescence can be produced by an intercalating dye that fluoresces when bound to the double-stranded DNA (e.g. SYBR green), or by a probe containing both a fluorophore and a quencher (e.g. Taqman) that binds between the primers. During primer extension, the probe is broken down, releasing the fluorophore from the quencher. The result of a reaction is expressed in a CT value which represents the Cycle number where the fluorescence overtakes a pre-set threshold value. The comparative CT method is considered to be the most robust for copy number determination . This technique requires the comparison of the CT value from the gene of interest (GOI) to the CT value of an endogenous reference gene. Benefits of this method are that no standard lines are needed and the DNA concentration is allowed to vary somewhat between samples, but it is crucial that both the target and reference reactions have near-identical amplification efficiencies. In theory, real-time PCR analyses will detect the inserted transgene regardless of where it is located within the genome.
Real-time PCR has been employed to determine transgene copy numbers in plants such as wheat [7, 8], maize , rice , tomato  and sugarcane . There are, however, concerns that real-time PCR is not accurate enough to reliably determine zygosity [7, 13], especially when the number of technical replicates is constrained due to the large number of plants requiring screening. However, we believe that previous results that found real-time PCR unreliable for copy numbers analyses and/or zygosity determination were likely due to suboptimal optimization procedures that can be overcome. In this case, real-time PCR may have been under-utilized in respect to its capacity to assist in the establishment of transgenic homozygous lines, especially in plants with longer generation times.
In this paper, we describe a real-time PCR based method with the required accuracy and precision to be able to distinguish plants with only a single copy of a transgene from homozygous plants and/or plants with multiple insertions. We show how this information can be best used in the development of transgenic homozygous lines hosting a unique insertion event. The method is specifically designed to be highly versatile and all procedures are standardized and described in sufficient detail to facilitate its adaptation to new organisms and/or vectors.
In the first experiment, we used DNA samples from offspring of self-fertilized plants which originated from crosses between a Glucan, Water Dikinase (GWD) RNAi transformed, single-insertion homozygous wheat line (based on observed Mendelian segregation history over 5 generations by endpoint PCR, as described in Ral et al.) and a wild type plant. The aim of this test was to determine if it was possible distinguish between single copy and homozygous plants using relative quantification (see Methods - real-time PCR copy numbers assay).
In a second experiment we followed four wheat plants, biolistically transformed with the RNAi cassette, over up to four generations (T0 to T3) to investigate the effectiveness of our assay in selecting homozygous transgenic lines.
DNA was isolated from leaf samples using a modified method obtained from Ellis et al.. Ca. 1 cm of young leaf sample was put into a 96x1 mL masterblock (Greiner Bio-one) together with one 5 mm glass bead (Sigma-Aldrich). The block was sealed with a sealing mat (Thermo Scientific) and frozen for at least 1 hour at −80°C. Next, the frozen block was placed in a Retch MM300 ball mill and the tissue was pulverized for 1 min at maximum frequency (30/s). Three hundred μL of extraction buffer (0.1 M Tris–HCl pH8, 0.05 M EDTA, 1.25% SDS, pre-warmed to 60°C) was added to each well after which the sealing mat was replaced. The block was placed in a custom-made metal clamp to prevent leaking, after which it was incubated for 1–2 h in a water bath at 60°C, with regular mixing by inverting. The block was cooled on ice-water, 150 μL of 6 M ammonium acetate was added and the block was sealed with PCR film. Samples were mixed by inverting 6–8 times and incubated in ice-water for 10 min. Next, the block was centrifuged for 20 min at 4,000 g (4°C) after which 250 μL of the supernatant was transferred to a new masterblock. One hundred and fifty μL of ice-cold 2-propanol was added, the block was sealed with PCR film and mixed by inverting 6–8 times and incubated in ice-water for 10 min. This step was followed by centrifuging for 20 min at 4000 g (4°C) after which the supernatant was poured off and the block dried on tissue paper. 150 μL of 70% ethanol was added, the centrifuge step was repeated and after the supernatant was poured off again the samples were allowed to dry for 5 min at room temperature. Finally, the precipitated DNA was re-dissolved in 400 μL 0.05 M Tris pH9 overnight at 4°C. This method routinely yields 10–30 ng/μL genomic DNA for wheat.
Endpoint PCR screening
PCR and real-time PCR primers
Primer sequence (5′ – 3′)
Product size (bp)
Endogenous gene; ca. 980
qPCR Rint4-9 F
qPCR Rint4-9 R
qPCR NOS F
qPCR NOS R
qPCR EC ABD F
qPCR EC ABD R
qPCR EC A F
qPCR EC A R
qPCR Rice SPS F
qPCR Rice SPS R
Real-time PCR primer development
All primers were generated using Primer3 software (http://frodo.wi.mit.edu/primer3/) using the following modifications from the default: Product size = 50 – 150 bp; Primer Tm = 59 – 65°C; Product Tm = 80 – 88°C; concentration of divalent cations = 3 mM, concentration of dNTPs = 0.2 mM; Table of thermodynamic parameters = SantaLucia 1998; Salt correction formula = Santalucia 1998.
The trialling of primer pairs consisted of three steps. First, primer pair specificity was determined by using two positive and two negative controls (using the same samples for all primer pairs), selecting for good melt curves (single peak around the expected temperature), whereas no-template controls were used to select primer pairs with no or only low primer dimer formation. If several primer pairs qualified, a pre-selection for efficiency was made by selecting for lowest CT values. The most promising primer pairs were used in sample dilution series using DNA of a single-copy GWD RNAi line (40 ng DNA/μL – 0.1 ng DNA/μL in 3-fold steps). To establish the dynamic range, the highest DNA concentrations of the dilution series are removed when they exhibit higher-than-expected CT values, an indicator of inhibited reactions. Those with the lowest DNA concentrations are removed when they show large variation between replicates, an indicator of the lower detection limit of the real-time PCR technique. The DNA concentration range left after reactions on both ends have been discarded where necessary constitutes the dynamic range and is used in a regression analysis to calculate R2 and efficiency values. Lastly, a test sample collection was assayed to determine assay accuracy and precision.
Real-time PCR copy numbers assay
Reactions were run in a MyIQ real-time PCR machine (BioRAD). A typical reaction consisted of: 10 μl Sensimix SYBR green with Fluorescein, 5 μL primer mix (1.4 μM of both forward and reverse primer), 5 μl DNA template (1 – 6 ng/uL). A standard 2-step protocol for Sensimix was followed: enzyme activation 10 min at 95°C, followed by 40 cycles of 95°C for 15 s, 60°C for 45 s. After each run, a melt-curve was produced in 1°C increments starting at 60°C to check for primer dimer formation. The fluorescent threshold was set at 150 Relative Fluorescent Units (RFU) in all runs. Reactions were run in triplicate. If the standard deviation within a triplicate was above 0.4 of a CT, the replicates were screened by eye for obvious outliers, which were removed leaving a duplicate with a standard deviation ≤ 0.4. Samples with a CT standard deviation value above 0.4 with no outliers were removed from further analyses, which occurred rarely. All copy numbers were calculated with the ΔCT method using a reference gene and calibrator sample (essentially the same as the ΔΔCT method  but calculated in a different order) on a home-made spreadsheet for the calculations. In short, the “raw” copy numbers were calculated with the formula: raw copy number = 2(CT(reference)-CT(Gene-Of-Interest)) × (number of reference gene copies). Next, raw values were divided by the raw value of the calibrator (a known single-copy sample included in all runs) to obtain the (corrected) copy numbers. The advantages of this variant over the standard ΔΔCT method were that (i) in our experimental procedure we found this method more intuitive to calculate, (ii) the raw copy numbers gave additional verification on how well the primer pairs were matched in efficiency (the closer the raw values are to the corrected values, the better the match), and (iii) when single copy samples were identified in a run, dividing the raw values by the mean of the raw values of the single copy samples (including the calibrator) enhanced the precision of copy number determination when compared to just using the raw value of the calibrator.
Primer pair selection and sample dilution
During sample dilutions, it was noted that samples with a DNA concentrations > 13 ng/μL had CT values higher than expected. Similarly, samples with concentrations < 0.3 ng/μL had increased variability between technical replicates. Samples of the dilution series outside the 0.3 – 13 ng/μL range were therefore not included in the regression analyses. Because of this, experimental samples were diluted to stay within a 1–6 ng/ μL range.
A 1 : 3.3 ratio between EC A and EC ABD was obtained. This is close to the expected 1 : 3 ratio and well within the error margin for real-time PCR when comparing two different primer pairs without a reference sample (i.e. a raw copy number). This result further indicates that our primer design and selection process produced well-matched primer pairs.
Discrimination between hemizygous and homozygous plants
EC A/EC ABD
The real-time PCR results of experiment 2 were verified in two ways. First, we screened 24 offspring of a hemizygous and a homozygous plant by endpoint PCR. Seven offspring of the single copy plant were negative, whereas all offspring from the homozygous plant were positive, confirming the real-time PCR results. Second, we screened 4 offspring of 3 selected homozygous plants with real-time PCR, confirming that they all possessed the same, even number of copies of the transgene.
The assay described here provided a powerful tool in the homozygous transgenic line development process. Ideally, the assay is used to select a single copy transgenic T0 plant as the assay is most precise in detecting a single copy, and in these lines there is no potential for gene-silencing due to multi-copy insertions. Screening 6 plants of the next generation was generally enough to identify a homozygous plant (expected frequency of 0.25) that could be used as founder of a transgenic homozygous line. Low copy number T0 plants can be screened in the next generation for a single copy segregant, which can then be used as the founding plant for a homozygous line carrying a single insertion. In T0 plants with multiple copies, back-crossings may be necessary to achieve a single-copy generation (not tested), although this will not work if the plants carry linked copies. A plant with putative linked copies can be made homozygous with a fair degree of certainty, however there remains uncertainty in these analyses about the exact copy numbers which is not present for single insertion homozygous lines (discussed below).
The required high precision was reliably achieved by following a few simple rules: (i) select primers that are (near-)100% efficient (by designing and testing several primer pairs per locus under standardized PCR conditions and selecting the one that is (near-)100% efficient); (ii) use a high-quality commercial qPCR master mix and (self-made) primer pre-mixes to minimize differences between runs by limiting the number of the pipetting steps; (iii) use good quality DNA at the right concentration (within the dynamic range). This approach consistently led to CT values with standard deviations of ≤ 0.4, as is required for copy number and zygosity determinations .
An important aspect of the identified ca. 30% error margin in these analyses is that there is very little uncertainty about the exact copy number of a single copy plant, but this uncertainty increases incrementally with increasing copy numbers. For instance, a result of 6 copies should be interpreted as having between 4 and 8 copies. This error margin may be reduced by running more technical replicates (not tested); however, for our analyses we deemed this not necessary as the main aim was to select plants with a unique insertion as efficiently as possible. To reflect this error margin, copy number results should be reported as integers. Additionally, to minimize the possibility of false positives, several offspring of an identified homozygous plant should be screened to verify copy number stability.
Bubner and Baldwin  showed that most previous studies have used Taqman probes for copy number analyses (one exception being ), suggesting that this is because these probes are superior since they are not affected by primer dimer formation. However, we considered that SYBR green has an important advantage over Taqman probes for copy number analyses and the focus on Taqman probes may be one of the reasons why previous studies have struggled to get to the required accuracy for reliable zygosity screening. In our experience, the use of SYBR green facilitates the selection of primer pairs that can amplify at near-100% efficiency. The much lower cost per primer pair makes it more feasible to order several primer pairs per target. Also, the chance of finding a suitable DNA sequence for your primers is higher with SYBR green (especially important when working within a small region) because the product can be anywhere between 60 and 200 bp and does not need to be suitable for a probe. Rather than optimising the PCR reaction to the primers, the PCR conditions are standardised and only primer pairs are selected that amplify with near-100% efficiency under the standard conditions and do not produce interfering primer dimers.
The one benefit of Taqman probes over SYBR in copy number determinations is the ability they provide for multiplexing, which is a powerful advantage for high-throughput screening. However, our experience is that the process of selecting homozygous lines is not a high-throughput procedure per se, but rather benefits from the flexibility to screen for many different constructs and easy adaptation procedures to new targets, as is the case for our assay. Any wheat construct containing one of the two standard sequence elements can be screened using the primers provided, and the example of rice indicates that the assay can easily be adapted to other plants.
Because bread wheat is hexaploid, we were able to make genome non-specific and genome specific primer pairs for the reference gene. Either pair could be used as a reference, but it is important to know how many copies the reference primer pair amplifies, because this will affect the raw transgene copy number estimate. Other organisms such as rice and barley do not have this complication and primer design is therefore somewhat easier.
Fitzgerald et al. have used genome-specific real-time PCR in wheat to detect Null-plants in large sample collections (plants with the gene of interest missing in one of the three genomes). We propose that our method can assist in the selection of double (and triple) null lines through crosses of these single null plants. With our method it should be possible to select the plants that are null for 1 genome and hemizygous for the other, which could be of great benefit especially if double nulls are elusive (e.g. due to negative selection pressure). A similar advantage can be gained when crossing multiple constructs into a single plant line (using construct-specific primer pairs). Finally, our assay could be useful in investigating endogenous gene copy number variation for low copy number genes.
Here we present a real-time PCR based method to fast-track the development of transgenic homozygous lines in commercially important plant species with relatively long generation times, such as wheat and rice. Standardized procedures for DNA isolation, primer design, primer selection and real-time PCR resulted in a method that is able to reliably distinguish single copy from low to high copy number plants, which can be done as early as the T0 generation. Zygosity determinations can be performed on their descendants to select homozygous lines carrying a single insertion, avoiding potential issues with transgene silencing and complex deregulation procedures for commercial use associated with plants with multiple insertions. We specifically designed the method to be relatively quick, affordable and versatile so it can be used with many different cereal constructs and may easily be adapted to different vectors and plant species. SYBR green is considered to have several advantages over Taqman probes for copy number analyses, including benefits in primer design, flexibility and cost-effectiveness. The method greatly assisted in the development of homozygous lines and promises a significant saving in time and collection management efforts.
We like to acknowledge Marianne Bloemsma for her assistance with the rice copy numbers analyses and proof-reading the manuscript, and Alex Whan for his assistance with the figures preparation. We furthermore like to thank Andrew Bowerman for the informative and critical discussions on real-time PCR, and the three reviewers for their useful comments on the manuscript. We acknowledge the Grains Research and Development Corporation of Australia for funding.
- Lemaux PG: Genetically engineered plants and foods: A scientist’s analysis of the issues (Part I). Annu Rev Plant Biol. 2008, 59: 771-812.PubMedView ArticleGoogle Scholar
- Bubner B, Baldwin IT: Use of real-time PCR for determining copy number and zygosity in transgenic plants. Plant Cell Rep. 2004, 23 (5): 263-271.PubMedView ArticleGoogle Scholar
- Butaye KMJ, Goderis IJWM, Wouters PFJ, Pues JMTG, Delaure SL, Broekaert WF, Depicker A, Cammue BPA, De Bolle MFC: Stable high-level transgene expression in Arabidopsis thaliana using gene silencing mutants and matrix attachment regions. Plant J. 2004, 39 (3): 440-449.PubMedView ArticleGoogle Scholar
- Sood P, Bhattacharya A, Sood A: Problems and possibilities of monocot transformation. Biol Plantarum. 2011, 55 (1): 1-15.View ArticleGoogle Scholar
- Tang W, Newton RJ, Weidner DA: Genetic transformation and gene silencing mediated by multiple copies of a transgene in eastern white pine. J Exp Bot. 2007, 58 (3): 545-554.PubMedView ArticleGoogle Scholar
- Bhat SR, Srinivasan S: Molecular and genetic analyses of transgenic plants: Considerations and approaches. Plant Sci. 2002, 163 (4): 673-681.View ArticleGoogle Scholar
- Gadaleta A, Giancaspro A, Cardone MF, Blanco A: Real-time PCR for the detection of precise transgene copy number in durum wheat. Cell Mol Biol Lett. 2011, 16 (4): 652-668.PubMedView ArticleGoogle Scholar
- Li ZW, Hansen JL, Liu Y, Zemetra RS, Berger PH: Using real-time PCR to determine transgene copy number in wheat. Plant Mol Biol Rep. 2004, 22 (2): 179-188.View ArticleGoogle Scholar
- Shou HX, Frame BR, Whitham SA, Wang K: Assessment of transgenic maize events produced by particle bombardment or Agrobacterium-mediated transformation. Mol Breeding. 2004, 13 (2): 201-208.View ArticleGoogle Scholar
- Yang LT, Ding JY, Zhang CM, Jia JW, Weng HB, Liu WX, Zhang DB: Estimating the copy number of transgenes in transformed rice by real-time quantitative PCR. Plant Cell Rep. 2005, 23 (10–11): 759-763.PubMedView ArticleGoogle Scholar
- German MA, Kandel-Kfir M, Swarzberg D, Matsevitz T, Granot D: A rapid method for the analysis of zygosity in transgenic plants. Plant Sci. 2003, 164 (2): 183-187.View ArticleGoogle Scholar
- Casu RE, Selivanova A, Perroux JM: High-throughput assessment of transgene copy number in sugarcane using real-time quantitative PCR. Plant Cell Rep. 2012, 31 (1): 167-177.PubMedView ArticleGoogle Scholar
- Bubner B, Gase K, Baldwin IT: Two-fold differences are the detection limit for determining transgene copy numbers in plants by real-time PCR. BMC Biotechnol. 2004, 4: 14.PubMedPubMed CentralView ArticleGoogle Scholar
- Ral JP, Bowerman AF, Li Z, Sirault X, Furbank R, Pritchard JR, Bloemsma M, Cavanagh CR, Howitt CA, Morell MK: Down-regulation of Glucan, Water-Dikinase activity in wheat endosperm increases vegetative biomass and yield. Plant Biotechnol J. 2012, 10 (7): 871-882.PubMedView ArticleGoogle Scholar
- Ellis MH, Rebetzke GJ, Azanza F, Richards RA, Spielmeyer W: Molecular mapping of gibberellin-responsive dwarfing genes in bread wheat. Theor Appl Genet. 2005, 111 (3): 423-430.PubMedView ArticleGoogle Scholar
- Cunningham FX, Gantt E: Genes and enzymes of carotenoid biosynthesis in plants. Annu Rev Plant Phys. 1998, 49: 557-583.View ArticleGoogle Scholar
- Jiang LX, Yang LT, Zhang HB, Guo JC, Mazzara M, Van den Eede G, Zhang DB: International collaborative study of the endogenous reference gene, sucrose phosphate synthase (SPS), used for qualitative and quantitative analysis of genetically modified rice. J Agr Food Chem. 2009, 57 (9): 3525-3532.View ArticleGoogle Scholar
- Fitzgerald TL, Kazan K, Li ZY, Morell MK, Manners JM: A high-throughput method for the detection of homoeologous gene deletions in hexaploid wheat. BMC Plant Biol. 2010, 10: 264.PubMedPubMed CentralView ArticleGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.