Selection and validation of reference genes for gene expression analysis in apomictic and sexual Cenchrus ciliaris
© Simon et al.; licensee BioMed Central Ltd. 2013
Received: 24 April 2013
Accepted: 25 September 2013
Published: 2 October 2013
Apomixis is a naturally occurring asexual mode of seed reproduction resulting in offspring genetically identical to the maternal plant. Identifying differential gene expression patterns between apomictic and sexual plants is valuable to help deconstruct the trait. Quantitative RT-PCR (qRT-PCR) is a popular method for analyzing gene expression. Normalizing gene expression data using proper reference genes which show stable expression under investigated conditions is critical in qRT-PCR analysis. We used qRT-PCR to validate expression and stability of six potential reference genes (EF1alpha, EIF4A, UBCE, GAPDH, ACT2 and TUBA) in vegetative and reproductive tissues of B-2S and B-12-9 accessions of C. ciliaris.
Among tissue types evaluated, EF1alpha showed the highest level of expression while TUBA showed the lowest. When all tissue types were evaluated and compared between genotypes, EIF4A was the most stable reference gene. Gene expression stability for specific ovary stages of B-2S and B-12-9 was also determined. Except for TUBA, all other tested reference genes could be used for any stage-specific ovary tissue normalization, irrespective of the mode of reproduction.
Our gene expression stability assay using six reference genes, in sexual and apomictic accessions of C. ciliaris, suggests that EIF4A is the most stable gene across all tissue types analyzed. All other tested reference genes, with the exception of TUBA, could be used for gene expression comparison studies between sexual and apomictic ovaries over multiple developmental stages. This reference gene validation data in C. ciliaris will serve as an important base for future apomixis-related transcriptome data validation.
KeywordsApomixis Apospory Reference genes qRT-PCR Buffelgrass Cenchrus ciliaris
Sexual and asexual modes of reproduction occur in flowering plants. The sexual reproduction pathway is highly regulated and results in the production of seed via fusion of male and female gametes. Apomixis is a natural form of asexual reproduction through seeds. There are two basic forms of apomixis, sporophytic and gametophytic. Sporophytic apomixis (or adventitious embryony) results in embryos that develop directly from nongenerative cells of the ovule. Gametophytic apomixis involves three steps: formation of a meiotically unreduced egg cell (apomeiosis), parthenogenetic development of this egg cell without fertilization, and formation of functional endosperm with (pseudogamy) or without (autonomous) fertilization of the central cell of the ovule [1, 2]. Apomeiosis in gametophytic apomixis consists of two forms: apospory and diplospory. In apospory, the embryo sacs develop from nucellar cells, while in diplospory the generative cell undergoes mitosis to form an embryo sac. In both apospory and diplospory, a chromosomally unreduced cell gives rise to the megagametophyte, in which the unreduced egg cell parthenogenetically develops into an embryo that is genetically identical to the maternal plant [1, 2]. Hence the successful introgression of apomixis into crop plants would greatly facilitate maintenance of hybrid vigor over successive generations and also reduce costs associated with hybrid seed production . Apomixis is a common mode of asexual reproduction in numerous families but is most frequent in the eudicot families Rosaceae and Asteraceae and in the monocot family Poaceae [4, 5]. The most common hypothesis behind apomictic reproduction is that it is evolved from the sexual pathway, possibly by deregulation in the timing of expression of sexual reproductive genes . The global regulatory effects of polyploidy and hybridity have also been proposed to be the possible triggers for conversion of sexual to apomictic forms of reproduction [1, 7]. Genetic data from several species now indicate that at least two genes, one for apomeiosis and one for parthenogenesis , are required for apomixis, although lack of recombination in some species manifests as monogenic inheritance of the trait.
Cenchrus ciliaris (buffelgrass, syn. Pennisetum ciliare) is found in tropical and southern Africa, the Mediterranean region, India and Pakistan. Most accessions of C. ciliaris are tetraploid with 2n = 4 × = 36 . Apospory in buffelgrass is conferred by an apospory-specific genomic region (ASGR) that is hemizygous in nature and recalcitrant for recombination [9–14]. C. ciliaris displays the 4-nucleate Panicum-type aposporous embryo sac pheno-type [15, 16]. A sexual genotype discovered in C. ciliaris was crossed as the female parent to apomictic genotypes producing progeny in three phenotypic classes, obligate sexual, obligate apomictic and facultative apomicts [14, 17, 18]. Regions of the ASGR in C. ciliaris have been studied via partial sequencing of ASGR-linked BAC (Bacterial Artificial Chromosome) clones, leading to the identification of various genes with putative transcription factor or signaling-related functions .
Differential gene expression studies between sexual and apomictic plants enables comparison between apomictic and sexual pathways. In Boechera, advances in cell isolation methods, along with next generation sequencing technology, have allowed global comparisons of gene expression patterns between sexual and apomictic reproductive tissues [20–22]. Quantitative real-time, reverse-transcriptase PCR (qRT-PCR) is one of the most extensively used methods to analyze and validate transcript expression profiles among different species, treatments or developmental stages due to its high sensitivity, specificity and broad quantification range [20, 23]. However accurate normalization steps should be followed to obtain reliable quantification of gene expression levels via qRT-PCR. The purpose of normalization is to correct any non-biological variability during the experimental procedure [24, 25]. Among the various normalization approaches, the use of internal controls or reference genes has become the method of choice [26, 27]. The success of using reference genes for proper normalization in qRT-PCR is highly dependent on the choice of the appropriate reference gene: its expression should be relatively constant across tissues and should not be significantly altered by the experimental conditions [28, 29]. Since reference genes often are housekeeping genes required for cellular survival, it is assumed that they are stably expressed across all tissues and/or treatments which often is not the case [24, 30]. Different statistical software like GeNorm  and NormFinder  are available to test expression stability of reference genes. The most commonly used reference genes for various qRT-PCR analyses are actin, glyceraldehyde-3-phosphate dehydrogenase, ribosomal genes, cyclophilin, elongation factor 1 alpha, adenine phosphoribosyl transferase and tubulin [33–38].
With respect to testing reference genes for apomictic gene expression analysis, two plants, Boechera and Brachiaria brizantha, have been used to validate reference genes in their respective apomictic and sexual accessions in different tissues and developmental stages [20, 39]. In Boechera, EF1alpha (Elongation factor 1 alpha) and ACT2 (Actin 2), were among the stable genes detected  and in Brachiaria brizantha, UBCE (ubiquitin conjugating enzyme), EIF4A (eukaryotic initiation factor 4A) and EF1 alpha were the most stable genes in both apomictic and sexual plants . TUBA (tubulin alpha) has been used as an internal control gene to normalize the qRT-PCR data for Pennisetum glaucum interspecific hybrids . In the present paper, we selected six reference genes – EF1alpha, EIF4A, ACT2, UBCE (chosen based on their stability in Boechera and Brachiaria brizantha), GAPDH (glyceraldehyde-3-phosphate dehydrogenase, the most commonly used reference gene in different systems), and TUBA (reported as unstable in Boechera and Brachiaria brizantha[20, 39] but used as an internal control for Pennisetum glaucum interspecific hybrids ). All the six reference genes were tested for gene expression stability in both sexual and apomictic C. ciliaris in multiple tissues but specifically in reproductive tissues encompassing four developmental stages for ovary and three developmental stages for anthers.
Plant materials and sample collection for qRT-PCR
Two genotypes of C. ciliaris (buffelgrass syn. Pennisetum ciliare) used in this study (obligate sexual, B-2S, and obligate aposporous, B-12-9) were maintained in the greenhouse with a temperature ranging from 24°C to 30°C. The plants were maintained by vegetative propagation as described in Roche et al..
Leaves and roots were collected from greenhouse-grown B-2S and B-12-9 plants, frozen in liquid nitrogen and stored at −80°C prior to total RNA extraction. To collect reproductive tissues, young inflorescences from B-12-9 and B-2S grown in the greenhouse were covered with pollination bags to prevent cross-pollination during the flowering season. Four developmental stages were collected for ovary samples based on anther developmental stages: stage I, early premeiotic; stage II, tetrad; stage III, DOP (the head contained one or two florets with anthers emerged, although tissue was collected from florets in which anthers had not yet emerged); and stage IV, DOP + 5 days (pollinated and seed were developing). Three different stages were collected for anthers: stage I, early premeiotic; stage II, tetrad; stage III, DOP (tissue was collected from anthers that had not yet emerged). The developmental stages I and II were determined by carbol fuchsin staining of anther squashes . Staged inflorescence segments were initially stored in RNALater (Ambion) at −20°C and later ovary and anther samples were dissected and stored in RLT lysis buffer (Qiagen) at −80°C. At least two biological replicates were collected for all samples.
Total RNA extraction and cDNA synthesis
Total RNA was extracted using the RNeasy plant mini kit (Qiagen) and then subjected to DNAase treatment (Turbo DNA free DNase, Ambion). The DNase treated RNA was concentrated using RNeasy mini elute columns (Qiagen). The final total RNA samples were checked by a Nanodrop spectrofluorometer (260/280, 260/230 ratios). RNA with 260/280 in the range of 1.9-2.0 and 260/230 > 2 were checked for RNA integrity using an Agilent Bioanalyzer. Only RNA with a RIN score > 7 was used for the experiment. First strand cDNA synthesis used 5 μg of total RNA with oligo-dT and superscript II enzyme (Invitrogen), according to the manufacturer’s instructions. The generated cDNAs were quantitated in triplicate using ribogreen , diluted to 3 ng/μl and stored in aliquots to avoid freeze thaw cycles.
Sequencing of conserved reference gene domains from B-2S and B-12-9
Candidate reference gene description and details of the first set of primer sequences designed from the conserved regions of the specific reference genes
Amplicon size (bp)
Initiation factor 4A
Primer design and real-time PCR assay
A set of real-time primers was designed from the aligned region of each of the six reference genes (EF1alpha, EIF4A, UBCE, GAPDH, ACT2 and TUBA) using Primer3, (v. 0.4.0) software  and oligo analysis was performed using IDT oligoanalyzer tool (http://www.idtdna.com) with care taken to ensure that the reference gene real-time primer sequences were from the conserved regions of B-2S and B-12-9 and did not encompass any SNPs (single nucleotide polymorphisms). Real-time RT-PCR amplification reactions were performed using SYBR Green detection chemistry and run on 384-well plates with the Light Cycler 4.8 (Roche Applied Science). To estimate PCR efficiency of each reference gene primer pair, standard curves were generated using cDNA samples derived from different organs (leaf, root, anther and ovary) for both B-2S and B-12-9 genotypes. Samples were run in a 2-fold serial dilution range across 5 points with initial concentration of template starting at 0.6 ng. The corresponding PCR efficiencies and the error values were determined by the Light Cycler 4.8 software. All samples to be analyzed were run in the same plate in triplicate. Reactions were prepared in a total volume of 10 μl containing: 2 μl of template, 1 μl of each reference gene primer pair (optimized concentration: 300 nM), 5 μl of 2x FastStart SYBR Green Master (Roche Applied Science) and 2 μl of sterile water. Each PCR analysis also included reverse transcription negative control (RT minus, without reverse transcriptase) to confirm the absence of genomic DNA and a no-template negative control to check for primer-dimer and contamination. Uracil-N-Glycosylase (1 μl per 100 μl of reaction mixture) was also included in the reaction mixture in order to avoid PCR product carryover contamination. The PCR cycling conditions were set as follows: pre-incubation at 37°C for 10 min to activate Uracil-N-Glycosylase, an initial denaturation step of 95°C for 10 min to inactivate Uracil-N-Glycosylase and activate the FastStart Taq DNA polymerase, followed by 45 cycles of denaturation at 95°C for 10s, annealing at 60°C for 10s and extension at 72°C for 10s. The amplification process was followed by a melting curve analysis, ranging from 65°C to 90°C, increasing temperature in steps of 0.2°C every 10s. Cp values were automatically determined using the Light Cycler 4.8 software (absolute quantitation via second-derivative method). For every cDNA sample, the mean expression level and standard deviation for each set replicate was calculated (cut-off value for standard deviation was kept as 0.3; only samples below this cut-off were considered). To confirm the reproducibility of the assay and to reconfirm reference gene stability, the experiment was repeated for B-2S and B-12-9 ovary samples. Two biological and three technical replicates were included.
Gene expression stability analysis
For ranking reference genes based on their stability, GenEx (version 4.1.7, MultiD Analyses) which includes both GeNorm and NormFinder software was used. The Microsoft Excel file with raw expression Cp values for each tested gene in the 18 different samples generated with the Light Cycler 4.8 software was transferred into GenEx software.
Primer sequences used for qRT-PCR analysis
Amplicon size (bp)
PCR Efficiency ± SD1
1.90 ± 0.03
1.90 ± 0.05
1.85 ± 0.03
1.87 ± 0.03
1.95 ± 0.04
2.05 ± 0.05
In order to minimize bias introduced by the reference gene validation approach, two programs with different algorithms, GeNorm and NormFinder, were used for measuring gene expression stability. GeNorm software assumes that none of the tested genes being analyzed are co-regulated . The stability measure provided by GeNorm (M-value) is the mean pairwise variation between a gene and all other tested candidates, and hence a pair of highly co-regulated genes could be eliminated during the selection if they show high inter-sample variability . Genes with the lowest M-value are considered most stable. M-values below 0.5 indicate good measure of stability [39, 46]. We also tested the NormFinder software since it is less sensitive to co-regulation . NormFinder relies on a ‘model-based approach’; it determines expression stability of candidate reference genes by comparing the expression variation between and within groups and then combines both results in a stability value for each tested gene . For NormFinder, the gene with the lowest standard deviation (SD) will be top ranked [28, 32]. NormFinder software was found to be more robust with larger sample sizes . Differences in gene ranking between different stability software have been reported .
Tissue specific reference gene stability analyses of B-2S (sex) and B-12-9 (apo) using GeNorm and NormFinder
GeNorm - M values
NormFinder_Gene ranking (SD values)
Ovary (early/tetrad/DOP/DOP + 5)
Ovary (second independent assay) reference gene stability analyses of B-2S (sex) and B-12-9 (apo) using GeNorm and NormFinder
GeNorm - M values
NormFinder_Gene ranking (SD values)
Ovary (early/tetrad/DOP/DOP + 5)
Our data present the first in-depth analysis of reference gene stability validation in C. ciliaris, a tetraploid and important model for understanding apomixis in grasses. Using both apomictic and sexual accessions of C. ciliaris and different organs/developmental stages, we were able to understand the reference gene stability specific to each organ, stage, or mode of reproduction. Our analysis used six reference genes based on their previous use in other plant species. These genes were analyzed in 18 samples, from both apomictic and sexual plants. All- tissue analysis using GeNorm and NormFinder found EIF4A as the most stable gene. There was tight correlation between analyses with GeNorm and NormFinder, as genes detected stable by GeNorm occupied top gene rankings in NormFinder analysis.
Detailed analyses of ovary tissue in two independent assays, confirmed stability of genes in specific ovary stages that could be used in both B-2S and B-12-9, irrespective of the mode of reproduction. Based on our results, all tested reference genes were found to be stable (as per GeNorm, M-value below 0.5) for early to tetrad stage ovary development between the two species, while ACT2, UBCE, GAPDH and/or EIF4A would be most suitable for ovary stages up to DOP and DOP + 5. These reference gene expression and stability analyses will provide an important guideline for our future study involving apomixis-related gene expression studies via qRT-PCR in C. ciliaris or other related grasses.
This work was supported by a cooperative research agreement with Pioneer Hi-Bred.
- Grimanelli D, Leblanc O, Perotti E, Grossniklaus U: Developmental genetics of gametophytic apomixis. Trends Genet. 2001, 17: 597-604. 10.1016/S0168-9525(01)02454-4.PubMedView ArticleGoogle Scholar
- Ozias-Akins P, Van-Dijk PJ: Mendelian genetics of apomixis in plants. Annu Rev Genet. 2007, 11: 509-537.View ArticleGoogle Scholar
- Spillane C, Curtis MD, Grossniklaus U: Apomixis technology development–virgin births in farmers’ fields?. Nat Biotechnol. 2004, 22: 687-691. 10.1038/nbt976.PubMedView ArticleGoogle Scholar
- Asker SE, Jerling L: Apomixis in Plants. 1992, Boca Raton, FL: CRC PressGoogle Scholar
- Richards AJ: Agamospery. Plant Breeding Systems. Edited by: Richards AJ. 1986, Boston: George Allen & Unwin, 403-509.Google Scholar
- Koltunow AM, Grossniklaus U: Apomixis: A developmental perspective. Annu Rev Plant Biol. 2003, 54: 547-574. 10.1146/annurev.arplant.54.110901.160842.PubMedView ArticleGoogle Scholar
- Carman JG: Asynchronous expression of duplicate genes in angiosperms may cause apomixis, bispory, tetraspory, and polyembryony. Biol J Linnean Soc. 1997, 61: 51-94. 10.1111/j.1095-8312.1997.tb01778.x.View ArticleGoogle Scholar
- Hignight KW, Bashaw EC, Hussey MA: Cytological and morphological diversity of native apomictic buffelgrass, Pennisetum ciliare (L.) Link. Bot Gaz. 1991, 152: 214-218. 10.1086/337882.View ArticleGoogle Scholar
- Akiyama Y, Conner JA, Goel S, Morishige DT, Mullet JE, DeBarry J, Yang L, Bennetzen JL, Klein PE, Ozias-Akins P: High-resolution physical mapping in Pennisetum squamulatum reveals extensive chromosomal heteromorphism of the genomic region associated with apomixis. Plant Physiol. 2004, 134: 1733-1741. 10.1104/pp.103.033969.PubMedPubMed CentralView ArticleGoogle Scholar
- Akiyama Y, Hanna WW, Ozias-Akins P: High-resolution physical mapping reveals that the apospory-specific genomic region (ASGR) in Cenchrus ciliaris is located on a heterochromatic and hemizygous region of a single chromosome. Theor Appl Genet. 2005, 111: 1042-1051. 10.1007/s00122-005-0020-5.PubMedView ArticleGoogle Scholar
- Goel S, Chen Z, Conner JA, Akiyama Y, Hanna WW, Ozias-Akins P: Physical evidence that a single hemizygous chromosomal region is sufficient to confer aposporous embryo sac formation in Pennisetum squamulatum and Cenchrus ciliaris. Genetics. 2003, 163: 1069-1082.PubMedPubMed CentralGoogle Scholar
- Goel S, Chen Z, Akiyama Y, Conner JA, Basu M, Gualtieri G, Hanna WW, Ozias-Akins P: Comparative physical mapping of the apospory-specific genomic region in two apomictic grasses: Pennisetum squamulatum and Cenchrus ciliaris. Genetics. 2006, 173: 389-400. 10.1534/genetics.105.054429.PubMedPubMed CentralView ArticleGoogle Scholar
- Ozias-Akins P, Roche D, Hanna WW: Tight clustering and hemizygosity of apomixis-linked molecular markers in Pennisetum squamulatum implies genetic control of apospory by a divergent locus which may have no allelic form in sexual genotypes. Proc Natl Acad Sci U S A. 1998, 95: 5127-5132. 10.1073/pnas.95.9.5127.PubMedPubMed CentralView ArticleGoogle Scholar
- Roche D, Cong P, Chen ZB, Hanna WW, Gustine DL, Sherwood RT, Ozias-Akins P: An apospory-specific genomic region is conserved between buffelgrass (Cenchrus ciliaris L.) and Pennisetum squamulatum Fresen. Plant J. 1999, 19: 203-208. 10.1046/j.1365-313X.1999.00514.x.PubMedView ArticleGoogle Scholar
- Nogler GA: Gametophytic apomixis. Embryology of Angiosperms. Edited by: Johri BM. 1984, Berlin, Heidelberg, New York, Tokyo: Springer, 475-518.View ArticleGoogle Scholar
- Vielle J-P, Burson BL, Bashaw EC, Hussey MA: Early fertilization events in the sexual and aposporous egg apparatus of Pennisetum ciliare (L.) Link. Plant J. 1995, 8: 309-316. 10.1046/j.1365-313X.1995.08020309.x.View ArticleGoogle Scholar
- Bashaw EC: Apomixis and sexuality in buffelgrass. Crop Sci. 1962, 2: 412-415. 10.2135/cropsci1962.0011183X000200050015x.View ArticleGoogle Scholar
- Gustine DL, Sherwood RT, Huff DR: Apospory-linked molecular markers in buffelgrass. Crop Sci. 1997, 37: 947-951. 10.2135/cropsci1997.0011183X003700030040x.View ArticleGoogle Scholar
- Conner JA, Goel S, Gunawan G, Cordonnier-Pratt MM, Johnson VE, Liang C, Wang H, Pratt LH, Mullet JE, Debarry J, Yang L, Bennetzen JL, Klein PE, Ozias-Akins P: Sequence analysis of bacterial artificial chromosome clones from the apospory-specific genomic region of Pennisetum and Cenchrus. Plant Physiol. 2008, 147: 1396-1411. 10.1104/pp.108.119081.PubMedPubMed CentralView ArticleGoogle Scholar
- Pellino M, Sharbel TF, Mau M, Amiteye S, Corral JM: Selection of reference genes for quantitative real-time PCR expression studies of microdissected reproductive tissues in apomictic and sexual Boechera. BMC Res Notes. 2011, 4: 303-10.1186/1756-0500-4-303.PubMedPubMed CentralView ArticleGoogle Scholar
- Sharbel TF, Voigt ML, Corral JM, Thiel T, Varshney A, Kumlehn J, Vogel H, Rotter B: Molecular signatures of apomictic and sexual ovules in the Boechera holboellii complex. Plant J. 2009, 58: 870-882. 10.1111/j.1365-313X.2009.03826.x.PubMedView ArticleGoogle Scholar
- Sharbel TF, Voigt ML, Corral JM, Galla G, Kumlehn J, Klukas C, Schreiber F, Vogel H, Rotter B: Apomictic and sexual ovules of Boechera display heterochronic global gene expression patterns. Plant Cell. 2010, 22: 655-10.1105/tpc.109.072223.PubMedPubMed CentralView ArticleGoogle Scholar
- Hong SY, Seo PJ, Yang MS, Xiang F, Park CM: Exploring valid reference genes for gene expression studies in Brachypodium distachyon by real-time PCR. BMC Plant Biol. 2008, 8: 112-10.1186/1471-2229-8-112.PubMedPubMed CentralView ArticleGoogle Scholar
- Pfaffl MW: Quantification strategies in real time PCR. A-Z of quantitative PCR. Edited by: Bustin SA. 2004, La Jolla, CA: International University Line, 1-20.Google Scholar
- Pfaffl MW, Tichopad A, Prgomet C, Neuvians TP: Determination of stable housekeeping genes, differentially regulated target genes and sample integrity: BestKeeper–Excel-based tool using pair-wise correlations. Biotechnol Lett. 2004, 26: 509-515.PubMedView ArticleGoogle Scholar
- Nolan T, Hands RE, Bustin SA: Quantification of mRNA using real-time RT-PCR. Nat Protoc. 2006, 1: 1559-1582. 10.1038/nprot.2006.236.PubMedView ArticleGoogle Scholar
- VanGuilder HD, Vrana KE, Freeman WM: Twenty-five years of quantitative PCR for gene expression analysis. Biotechniques. 2008, 44: 619-626.PubMedView ArticleGoogle Scholar
- Expósito-Rodríguez M, Borges AA, Borges-Pérez A, Pérez JA: Selection of internal control genes for quantitative real-time RT-PCR studies during tomato development process. BMC Plant Biol. 2008, 8: 131-10.1186/1471-2229-8-131.PubMedPubMed CentralView ArticleGoogle Scholar
- Huggett J, Dheda K, Bustin SA: Normalization. Real-Time PCR. Edited by: Dorak MT. 2006, New York: BIOS Advanced Methods, 83-91.Google Scholar
- Thellin O, Zorzi W, Lakaye B, De-Borman B, Coumans B, Hennen G, Grisar T, Igout A, Heinen E: Housekeeping genes as internal standards: use and limits. J Biotechnol. 1999, 75: 291-295. 10.1016/S0168-1656(99)00163-7.PubMedView ArticleGoogle Scholar
- Vandesompele J, De-Preter K, Pattyn F, Poppe B, Van-Roy N, De-Paepe A, Speleman F: Accurate normalization of real-time quantitative RT-PCR data by geometric averaging of multiple internal control genes. Genome Biol. 2002, 3: 0034.1-0034.11.View ArticleGoogle Scholar
- Andersen CL, Jensen JL, Orntoft TF: Normalization of real-time quantitative reverse transcription-PCR data: a model-based variance estimation approach to identify genes suited for normalization, applied to bladder and colon cancer data sets. Cancer Res. 2004, 64: 5245-5250. 10.1158/0008-5472.CAN-04-0496.PubMedView ArticleGoogle Scholar
- Czechowski T, Stitt M, Altmann T, Udvardi MK, Scheible WR: Genome-wide identification and testing of superior reference genes for transcript normalization in Arabidopsis. Plant Physiol. 2005, 139: 5-17. 10.1104/pp.105.063743.PubMedPubMed CentralView ArticleGoogle Scholar
- Dean JD, Goodwin PH, Hsiang T: Comparison of relative RT-PCR and northern blot analyses to measure expression of β-1,3-glucanase in Nicotiana benthamiana infected with Colletotrichum destructivum. Plant Mol Biol Rep. 2002, 20: 347-356. 10.1007/BF02772122.View ArticleGoogle Scholar
- Orsel M, Krapp A, Daniel-Vedele F: Analysis of the NRT2 nitrate transporter family in Arabidopsis: Structure and gene expression. Plant Physiol. 2002, 129: 886-888. 10.1104/pp.005280.PubMedPubMed CentralView ArticleGoogle Scholar
- Ozturk ZN, Talamé V, Deyholos M, Michalowski CB, Galbraith DW, Gozukirmizi N, Tuberosa R, Bohnert HJ: Monitoring large-scale changes in transcript abundance in drought- and salt-stressed barley. Plant Mol Biol. 2002, 48: 551-573. 10.1023/A:1014875215580.View ArticleGoogle Scholar
- Stürzenbaum SR, Kille P: Control genes in quantitative molecular biological techniques: the variability of invariance. Comp Biochem Physiol B. 2007, 130: 281-289.View ArticleGoogle Scholar
- Thomas C, Meyer D, Wolff M, Himber C, Alioua M, Steinmetz A: Molecular characterization and spatial expression of the sunflower ABP1 gene. Plant Mol Biol. 2003, 52: 1025-1036. 10.1023/A:1025482432486.PubMedView ArticleGoogle Scholar
- Silveira ED, Alves-Ferreira M, Guimarães LA, Da-Silva FR, Carneiro VTC: Selection of reference genes for quantitative real-time PCR expression studies in the apomictic and sexual grass Brachiaria brizantha. BMC Plant Biol. 2009, 9: 84-10.1186/1471-2229-9-84.PubMedPubMed CentralView ArticleGoogle Scholar
- Sahu PP, Gupta S, Malaviya DR, Roy AK, Kaushal P, Prasad M: Transcriptome analysis of differentially expressed genes during embryo sac development in apomeiotic non-parthenogenetic interspecific hybrid of Pennisetum glaucum. Mol Biotechnol. 2012, 51: 262-271. 10.1007/s12033-011-9464-9.PubMedView ArticleGoogle Scholar
- Zeng Y, Conner J, Ozias-Akins P: Identification of ovule transcripts from the Apospory-Specific Genomic Region (ASGR)-carrier chromosome. BMC Genomics. 2011, 12: 206-10.1186/1471-2164-12-206.PubMedPubMed CentralView ArticleGoogle Scholar
- Libus J, Štorchová H: Quantification of cDNA generated by reverse transcription of total RNA provides a simple alternative tool for quantitative RT-PCR normalization. Biotechniques. 2006, 41: 156-164. 10.2144/000112232.PubMedView ArticleGoogle Scholar
- Rozen S, Skaletsky HJ: Primer3 on the WWW for general users and for biologist programmers. Bioinformatics Methods and Protocols: Methods in Molecular Biology. Edited by: Krawetz S, Misener S. 2000, Totowa, NJ: Humana Press, 365-386.Google Scholar
- Zhao SFR: Comprehensive algorithm for quantitative real-time polymerase chain reaction. J Comp Biol. 2005, 12: 1047-1064. 10.1089/cmb.2005.12.1047.View ArticleGoogle Scholar
- Bustin SA, Nolan T: Pitfalls of quantitative real-time reverse-transcription polymerase chain reaction. J Biomol Tech. 2004, 15: 155-156.PubMedPubMed CentralGoogle Scholar
- Yan J, Yuan F, Long G, Qin L, Deng Z: Selection of reference genes for quantitative real-time RT-PCR analysis in citrus. Mol Biol Rep. 2012, 39: 1831-1838. 10.1007/s11033-011-0925-9.PubMedView ArticleGoogle Scholar
- Demidenko NV, Logacheva MD, Penin AA: Selection and validation of reference genes for quantitative real-time PCR in buckwheat (Fagopyrum esculentum) based on transcriptome sequence data. PLoS One. 2011, 6: e19434-10.1371/journal.pone.0019434.PubMedPubMed CentralView ArticleGoogle Scholar
- Mehta R, Birerdinc A, Hossain N, Afendy A, Chandhoke V, Younossi Z, Baranova A: Validation of endogenous reference genes for qRT-PCR analysis of human visceral adipose samples. BMC Mol Biol. 2010, 11: 39-10.1186/1471-2199-11-39.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.