Methodology

A plastome primer set for comprehensive quantitative real time RT-PCR analysis of Zea mays: a starter primer set for other Poaceae species

Richard M Sharpe , Sade N Dunn and A Bruce Cahoon

Department of Biology, Middle Tennessee State University, Murfreesboro, TN, USA

author email corresponding author email

Plant Methods 2008, 4:14doi:10.1186/1746-4811-4-14

The electronic version of this article is the complete one and can be found online at: http://www.plantmethods.com/content/4/1/14

Received:	1 February 2008
Accepted:	2 June 2008
Published:	2 June 2008

© 2008 Sharpe et al; licensee 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.

Abstract

Background

Quantitative Real Time RT-PCR (q2(RT)PCR) is a maturing technique which gives researchers the ability to quantify and compare very small amounts of nucleic acids. Primer design and optimization is an essential yet time consuming aspect of using q2(RT)PCR. In this paper we describe the design and empirical optimization of primers to amplify and quantify plastid RNAs from Zea mays that are robust enough to use with other closely related species.

Results

Primers were designed and successfully optimized for 57 of the 104 reported genes in the maize plastome plus two nuclear genes. All 59 primer pairs produced single amplicons after end-point reverse transcriptase polymerase chain reactions (RT-PCR) as visualized on agarose gels and subsequently verified by q2(RT)PCR. Primer pairs were divided into several categories based on the optimization requirements or the uniqueness of the target gene. An in silico test suggested the majority of the primer sets should work with other members of the Poaceae family. An in vitro test of the primer set on two unsequenced species (Panicum virgatum and Miscanthus sinensis) supported this assumption by successfully producing single amplicons for each primer pair.

Conclusion

Due to the highly conserved chloroplast genome in plant families it is possible to utilize primer pairs designed against one genomic sequence to detect the presence and abundance of plastid genes or transcripts from genomes that have yet to be sequenced. Analysis of steady state transcription of vital system genes is a necessary requirement to comprehensively elucidate gene expression in any organism. The primer pairs reported in this paper were designed for q2(RT)PCR of maize chloroplast genes but should be useful for other members of the Poaceae family. Both in silico and in vitro data are presented to support this assumption.

Background

Chloroplasts are semi-autonomous organelles believed to have developed from free-living photosynthetic bacteria [1,2] They are members of a diverse and flexible family of organelles called plastids that are responsible for photosynthesis plus other essential biosynthetic pathways and cellular functions. Plastids have maintained a small remnant genome with a species-specific number of genes mostly involved in photosynthesis and gene expression. The full function and development of all the plastid types, however, requires thousands of nuclear encoded gene products.

Maize is an agriculturally important monocot grass that has served as a genetic model system for decades [3], is the focus of a major genome project [4,5], and is especially valuable for the study of chloroplast biology [6]. Other members of the grass family (Poaceae) also offer unique opportunities to study differential plastid gene expression. The presence of closely related species with either C3 or C4 photosynthetic capabilities enables comparison of dimorphic C4 chloroplast development to monomorphic C3 species. In addition, grasses such as switchgrass (Panicum virgatum) and miscanthus (Miscanthus sinensis) are the subjects of an increased focus on cellulosic ethanol production [7].

To date, most Poaceae chloroplast gene expression studies have concentrated on proteomic or physiologic assays or focused on post-transcriptional modification and regulation of transcripts, [8-10]. While these approaches elucidate end point development of transcribed genes, they rarely illuminate the transcription activity of a particular gene. q2(RT)PCR is a maturing tool sensitive enough to detect the existence of small amounts of nucleic acid [11]. This allows for in-depth, comprehensive investigations into transcript abundance and offers a useful tool to help elucidate the relationship between transcription, translation, and expression.

q2(RT)PCR is dependent upon the amplification and quantification of a single amplicon that makes primer design and amplification conditions key to the success of an experiment. Substantial time and resources may be spent in the design, testing, and subsequent reworking of primers for optimal efficiency [12].

In this paper we describe the development and empirical optimization of primer pairs to amplify each rRNA and mRNA from maize plastids. Primer sequences, optimal annealing temperatures, and extension times are reported. In addition, each primer set was tested, in silico, against published plastome sequences and in vitro against switchgrass and miscanthus transcripts. Using the conditions optimized for maize, all primer pairs successfully produced a single amplicon for these two grass species.

Results and Discussion

One of the challenging and time consuming aspects of q2(RT)-PCR is the design and optimization of primer pairs which yield single amplicons. The aim of this study was to design and optimize a comprehensive set of plastid specific primers for q2(RT)PCR specific enough to yield robust steady state transcript data from maize yet flexible enough to recognize transcripts from multiple members of the family Poaceae. To accomplish this, primers were designed against highly homologous protein coding regions from the maize plastome, as well as rRNA genes. Primers were designed to produce amplicons 75–150 bp in length, have similar annealing temperatures, and were carefully evaluated for favorable melting temperatures to insure a lack of intra-molecular folding. Each primer set was qualitatively examined and optimized with endpoint RT-PCR as seen in Figure 1(A–E). These optimal conditions were then utilized in q2(RT)-PCR for each primer pair to verify the amplification of single amplicons via melt curves and their ability to produce quantitative data as seen in figure 1(F) and 1(G). The complete set of primer pairs yielding single amplicons and their optimized conditions is listed in Additional File 1.

Throughout the design and optimization process, we found it useful to categorize the primers based on their ease of use and optimal amplification conditions. For example, primer pairs that produced a single amplicon over a relatively large range of annealing temperatures and elongation times as shown in figure 1A were designated category 1 primers. Category 2 primer pairs produced multiple products at non-optimal annealing temperatures or elongation times and are shown in figure 1B,1C,1D, and 1E; category 3 primer pairs produced multiple amplicons due to the presence of multiple organellar genomes or primers having a high affinity to unknown additional nucleotide sequences. Category 1 primer pairs consistently were the most forgiving of strict protocol adherence. Elongation and annealing temperatures could be adjusted with no visible endpoint differences. Category 2 primer pairs required stricter elongation times and annealing temperatures for optimal single product results. Subcategory 2a includes primer pairs where an elongation time greater than 30 seconds resulted in multiple erroneous amplicons. Subcategory 2b includes primer pairs requiring specific annealing temperatures either to generate the most product as seen in Figure 1D or to produce only one product as shown in Figure 1E. One lesson learned is that some primer pairs initially in the 2b subcategory could be moved to category one if the primers were re-designed such that their individual melting temperatures were more similar. Figure 1B illustrates inconsistent and nonspecific product bands from the first pair of primers designed for psaB that had individual primer melting point temperatures that differed by 4°C. The subsequent set of primers designed for psaB have more similar melting temperatures and optimized end point results can be seen in Figure 1C. The third category of primer pairs consist of gene sequences where a highly similar homolog exists either in the nucleus or mitochondria (4.5S, 16S, atpE, rpl2, and 23S) and\or genes where multiple attempts to design primer pairs produced unacceptable products (rps12 and rps7). Two category 3 examples, 4.5S and 16S, are shown in Figure 2. For 4.5S, an end-point reaction utilizing suboptimal primers appeared to produce a single band but a melt peak curve clearly showed three very similar but distinct products. A subsequent redesign of the 3' primer for 4.5S and 23S resulted in single amplicons verified by their respective melt curve signatures and are reported as the 3' primers in Additional File 1. 16S, on the other hand, should also produce multiple bands according to screens of the maize genomes but neither end-point nor q2(RT)PCR can distinguish separate products. We could not successfully design some category 3 primer pairs. They either repeatedly produced multiple amplicons or could not be verified as amplified from the chloroplast genome. These primer pairs have been excluded from Additional File 1 or do not have annealing temperatures associated with them.

Some plastid genes with intra-cellular homologs could be specifically amplified when primers were designed to take advantage of small non-similar regions within the genes. This was achieved by comparing the homology between the genes and designing primers to recognize unique regions within chloroplast genes. Primers were then compared to other competing homologs, if the primer homology was less than 50% and there were no 3' matching nucleotides, then the design process was completed and the primers synthesized. The similarity estimation was based on the BLAST [13] percent identity score and not the E value which can be inflated when sequence sizes are small. One exception to this design strategy was the atpA gene primers which have >50% homology and should anneal to mitochondrial homologous regions. The possible mitochondrial amplicon however was predicted to be 125,630 base pairs long. An experimental end-point check of the atpA primers amplified the desired 77 base pair chloroplast amplicon and the single product was verified by a q2(RT)PCR melt curve. Multi-genomic genes exhibiting high homology but whose primer annealing sites and subsequent amplicons are plastid specific when analyzed in this manner are reported in Additional File 1 as category 1 or category 2 genes. Multi-genomic genes with homologies below 50% and whose primer pairs were determined to amplify only chloroplast genomic material are also reported in Additional File 1 as category 1 or category 2 genes.

The homology between plastomes from closely related plant species is well documented [14-16]. This homology has been previously utilized to design PCR primer pairs that would produce amplicons across genera [17,18]. Based on this, we hypothesized that our primer pairs should amplify similar regions of the same genes from other members of the grass family (Poaceae). To test this, an in silico comparison of the primers to grasses with completed plastomes was conducted. The possibility of a successful reaction was judged based on overall percent homology and the presence of mismatched bases at the 3' end of the primers. This comparison suggested that primer homology and subsequent amplification success would closely follow established taxonomic designations (Figure 3 and Additional File 2). Members of the PACCAD clade (S. officinarum and S. bicolor), to which maize belongs, have the greatest possibility of successful amplifications for all primer pairs while there is a noticeable drop in predicted success among members of the BEP clade (H. vulgare, T. aestivum, L. perrene). Arabidopsis was included as an outgroup, with only 22% of the primers predicted to successfully anneal to the intended target gene.

The in silico results were tested in vitro with two species in the PACCAD clade with unsequenced (or unpublished) plastomes – switchgrass (Panicum virgatum) and miscanthus (Miscanthus. sinensis). As predicted, we found that the primers and their empirically determined optimal conditions worked well across these genera and produced reliable quantitative data [Figure 4]. Representative data from the chloroplast Large Single Copy area (rpoC1 and matK), Small Single Copy area (psaC), and the Inverted Repeat (16S) areas are illustrated in Figure 5. As shown, these primer sets and amplification conditions eliminated the need for additional primer design for these species.

As seen in Figure 5, there are obvious differences in expression amounts and amplification slopes between maize (A), switchgrass (B) and miscanthus (C). Although equal amounts of total RNA were added for each round there was no attempt to normalize the data to account for any differences that may exist between species. It also must be stressed that the primer pairs have been optimized for the maize plastome and although all were checked and verified to produce single amplicons for both switchgrass and miscanthus, the reported annealing temperatures were not optimized to ensure the greatest signal was detected for these two species. Individual reagents, specific types of equipment and reaction volumes may require some user-specific modifications for optimal performance for species specific work.

Conclusion

Primers were designed and empirically optimized to quantify plastid transcripts from Zea mays. Due to the highly conserved nature of plastomes within the family Poaceae, this set of primers works equally as well for the majority of targeted genes of other species in the Poaceae family, particularly members of the PACCAD clade. This primer set and the conditions published in this report should reduce the time and cost associated with primer design and testing for future plastid q2(RT)PCR experiments.

Methods

q2(RT)PCR Primer Design and validation

q2(RT)PCR primer design and validation procedural flow can be seen in Figure 6. Primer pairs were designed from the NCBI published Zea mays complete chloroplast genome sequence accession numbers [GenBank: X86563 and GenBank: NC_001666], NAPD-Malic enzyme accession number [GenBank: J05130], and NAPD-Malate dehydrogenase accession number [GenBank: AY105634] [19]. Gene sequences were obtained from each coding sequence (CDS) listing, rRNA listing, or accession number sequence. The coding sequence was then analyzed using the BLAST [13] program for homology to mitochondrial and known nuclear maize genomes. Primers were designed with an attempt to utilize dissimilar sequences between the genomes using the online software primer3 [20]. Suggested primers were checked for inhibitory secondary structures using mFold [21] and edited until the recommended primer produced a positive delta G at 45°C. This ensured selected primer pairs would possess a null or positive delta G energy signature at reverse transcriptase annealing temperatures and amplify chloroplast specific gene sequences. Primer pairs were then optimized to produce a single amplicon of predicted length by end-point temperature gradient PCR with maize whole leaf tissue RNA template and visualized on a 3% TAE agarose gel run alongside a 2-Log ladder from New England Biolabs (Ipswitch, MA). Optimization of primer pairs consisted of running eight temperature gradient dependent reactions and visualizing them on a 3% agarose gel. The gradient spanned a 10°C range, ± 5° from the average melting temperatures (Tm) of the primer pairs as reported by the supply company. Reaction elongation times were set at two minutes in the initial primer set testing phase. Several of the initial primer pairs produced amplicons of greater length than the predicted length. Elongation times for these and subsequent primer pairs were reduced to 30 seconds to better match calculated elongation time [22] for amplicons of the desired size. Amplicons from switchgrass and miscanthus were then verified as the only product via end point RT-PCR visualized on 3% TAE agarose gel and then subsequently verified by the q2(RT)PCR melt curve for a single amplicon with BioRad's IQ™5 Optical System software, Hercules, CA,. All primers were obtained from Integrated DNA Technologies, Coralville, IA.

RNA material preparation

Maize (inbred line T42) stage 2 leaves were dissected into basal, mid, and tip sectional tissue as described in Cahoon, et al [23] and mixed into a single sample. Switchgrass stem and leaf and miscanthus leaves were harvested from mature demonstration specimens maintained by the MTSU greenhouse. Tissue samples were cut into approximately 2 mm by 2 mm sections, placed in liquid nitrogen, and ground with ceramic mortars and pestles. RNA was extracted with Trizol^® according to the manufacturer's suggested protocol, Invitrogen (Carlsbad, CA). Each extraction was checked for a DNA/RNA concentration utilizing a Spectronic^® Genesys™ spectrophotometer (Lincoln, NE). Ambion's Turbo DNase-free™ or DNase-free™, (Applied Biosystems, Foster City, CA) was used according to manufacturers suggested protocol to remove contaminating DNA. DNase treatments were repeated until each RNA extraction was visually confirmed to be DNA free by end point PCR reactions.

Gene Expression Verification

Each gene primer pair was verified utilizing BioRad's (Hercules, CA) iCycler and IQ™5 Multi-color Real-Time PCR Detection System. 25 μl reactions were prepared with BioRad iScript™ One-Step RT-PCR kit with SYBER^®Green utilizing approximately 200 ng/ul cleaned RNA per reaction for each of the three tissue samples. Amplicon presence and purity was confirmed by analysis of the IQ™5 Optical System's melt curve, melt curve peak, and PCR amplification software.

Abbreviations

RT-PCR: reverse transcriptase PCR; q2(RT)PCR: quantitative real time reverse transcriptase PCR

Competing interests

The authors declare that they have no competing interests.

Authors' contributions

RMS designed the primers, performed the q2(RT)PCR reactions, interpreted reaction results and prepared the manuscript, SND processed and performed the q2(RT)PCR for the majority of the switchgrass samples, ABC was involved in the conception and design of the project, revised and edited the manuscript, and is the PI of the laboratory.

Acknowledgements

Work was funded by an internal research enhancement program awarded to ABC. The authors wish to thank Larry Sizemore of the MTSU greenhouse facility for providing switchgrass and miscanthus samples, Drs. Rebecca Seipelt and Matthew Elrod-Erickson for their critical reading, and members of "Lab 122" for their helpful discussions and insight. Dai Gu contributed by designing the rpl14 and rps16 primer pairs.

References

1. Vothknecht UC, Westhoff P: Biogenesis and origin of thylakoid membranes.

   Biochim Biophys Acta 2001, 1541:91-101. PubMed Abstract | Publisher Full Text

2. De Las Rivas J, Lozano JJ, Ortiz AR: Comparative analysis of chloroplast genomes: functional annotation, genome-based phylogeny, and deduced evolutionary patterns.

   Genome Res 2002, 12:567-583. PubMed Abstract | Publisher Full Text | PubMed Central Full Text

3. Brendel V, Kurtz S, Walbot V: Comparative genomics of Arabidopsis and maize: prospects and limitations.

   Genome Biol 2002, 3:REVIEWS1005. PubMed Abstract | BioMed Central Full Text | PubMed Central Full Text

4. Chandler VL, Brendel V: The Maize Genome Sequencing Project.

   Plant Physiol 2002, 130:1594-1597. PubMed Abstract | Publisher Full Text | PubMed Central Full Text

5. Dong Q, Roy L, Freeling M, Walbot V, Brendel V: ZmDB, an integrated database for maize genome research.

   Nucleic Acids Res 2003, 31:244-247. PubMed Abstract | Publisher Full Text | PubMed Central Full Text

6. Stern DB, Hanson MR, Barkan A: Genetics and genomics of chloroplast biogenesis: maize as a model system.

   Trends Plant Sci 2004, 9:293-301. PubMed Abstract | Publisher Full Text

7. Morrow WR, Griffin WM, Matthews HS: Modeling switchgrass derived cellulosic ethanol distribution in the United States.

   Environ Sci Technol 2006, 40:2877-2886. PubMed Abstract | Publisher Full Text

8. Majeran W, Cai Y, Sun Q, van Wijk KJ: Functional differentiation of bundle sheath and mesophyll maize chloroplasts determined by comparative proteomics.

   Plant Cell 2005, 17:3111-3140. PubMed Abstract | Publisher Full Text | PubMed Central Full Text

9. Darie CC, Biniossek ML, Winter V, Mutschler B, Haehnel W: Isolation and structural characterization of the Ndh complex from mesophyll and bundle sheath chloroplasts of Zea mays.

   Febs J 2005, 272:2705-2716. PubMed Abstract | Publisher Full Text

10. Romanowska E, Drozak A, Pokorska B, Shiell BJ, Michalski WP: Organization and activity of photosystems in the mesophyll and bundle sheath chloroplasts of maize.

    J Plant Physiol 2006, 163:607-618. PubMed Abstract | Publisher Full Text

11. Ramakers C, Ruijter JM, Deprez RH, Moorman AF: Assumption-free analysis of quantitative real-time polymerase chain reaction (PCR) data.

    Neurosci Lett 2003, 339:62-66. PubMed Abstract | Publisher Full Text

12. Brunner AM, Yakovlev IA, Strauss SH: Validating internal controls for quantitative plant gene expression studies.

    BMC Plant Biol 2004, 4:14. PubMed Abstract | BioMed Central Full Text | PubMed Central Full Text

13. Altschul SF, Gish W, Miller W, Myers EW, Lipman DJ: Basic local alignment search tool.

    J Mol Biol 1990, 215:403-410. PubMed Abstract | Publisher Full Text

14. Demesure B, Sodzi N, Petit RJ: A set of universal primers for amplification of polymorphic non-coding regions of mitochondrial and chloroplast DNA in plants.

    Mol Ecol 1995, 4:129-131. PubMed Abstract | Publisher Full Text

15. Dumolin-Lapegue S, Pemonge MH, Petit RJ: An enlarged set of consensus primers for the study of organelle DNA in plants.

    Mol Ecol 1997, 6:393-397. PubMed Abstract | Publisher Full Text

16. Hamilton MB: Four primer pairs for the amplification of chloroplast intergenic regions with intraspecific variation.

    Mol Ecol 1999, 8:521-523. PubMed Abstract | Publisher Full Text

17. Dhingra A, Folta KM: ASAP: amplification, sequencing & annotation of plastomes.

    BMC Genomics 2005, 6:176. PubMed Abstract | BioMed Central Full Text | PubMed Central Full Text

18. Heinze B: A database of PCR primers for the chloroplast genomes of higher plants.

    Plant Methods 2007, 3:4. PubMed Abstract | BioMed Central Full Text | PubMed Central Full Text

19. NCBI: The National Center for Biotechnology Information. [http://www.ncbi.nlm.nih.gov/blast] webcite

20. Rozen S, Skaletsky HJ: Primer3 on the WWW for general users and for biologist programmers. In Methods Mol Biol. Edited by: Krawetz S, Misener S. Humana Press, Totowa, NJ; 2000:365-386. PubMed Abstract

21. Zuker M: Mfold web server for nucleic acid folding and hybridization prediction.

    Nucleic Acids Res 2003, 31:3406-3415. PubMed Abstract | Publisher Full Text | PubMed Central Full Text

22. Sambrook J, Russell DW: Molecular cloning: a laboratory manual. 3rd edition. Cold Spring Harbor Laboratory Press, New York; 2001.

23. Cahoon AB, Harris FM, Stern DB: Analysis of developing maize plastids reveals two mRNA stability classes correlating with RNA polymerase type.

    EMBO Reports 2004, 5:801-806. PubMed Abstract | Publisher Full Text | PubMed Central Full Text

24. Cahoon AB, Takacs EM, Sharpe RM, Stern DB: Nuclear, Chloroplast, and Mitochondrial Transcript Abundance Along a Maize Leaf Developmental Gradient.

    Plant Molecular Biology 2008, 66:33-46. PubMed Abstract | Publisher Full Text

Have something to say? Post a comment on this article!