PCR reactions were performed in 96-well plates with an ABI PRISM® 7500 Sequence Detection System (Applied Biosystems) using SYBR® Green to detect dsDNA synthesis. Reactions were done in 25 μl volumes containing 200 nM of each primer, 5 μl cDNA (corresponding to ~3 ng), and 12.5 μl 2× SYBR Green Master Mix Reagent (Applied Biosystems). Aliquots from the same cDNA sample were used with all primer sets in each experiment. Reactions were run using the manufacturer's recommended cycling parameters of 50°C for 2 min, 95°C for 10 min, 40 cycles of 95°C for 15 s, and 60°C for 1 min. No-template controls were included for each primer pair and each PCR reaction was completed in triplicate. Dissociation curves for each amplicon were then analyzed to verify the specificity of each amplification reaction; the dissociation curve was obtained by heating the amplicon from 60°C to 95°C (See additional file 1 : Dissociation curve data).
Data were analyzed using the SDS 1.2.2 software (Applied Biosystems). Expression levels were determined as the number of amplification cycles needed to reach a fixed threshold in the exponential phase of the PCR reaction (Ct). All amplification plots were analyzed with an Rn threshold of 0.2 to obtain Ct values. The PCR efficiency was determined for each gene with LinReg software, which uses absolute fluorescence data captured during the exponential phase of amplification of each reaction [23 (link)]. Results from the SDS and LinReg software were imported into Microsoft Excel for further analyses and to correct for the different PCR efficiencies [24 (link)]. All primer pairs had efficiencies higher than 1.80 with the exception of EF1-α (m) (1.70) and β-tubulin (1.54). Each was run on the full pericarp developmental series for 2003 and 2004. In addition, all primer pairs except those targeting β-tubulin and cyclophilin were run on 2003 mesocarp samples.
In order to evaluate reference gene stability among samples, three statistical approaches were incorporated. In the first approach, Ct difference (Ct max-Ct min) and CV were calculated for each gene throughout each development series tested. During the second approach, Ct values were converted into relative quantities and imported into geNorm v.3.4 software [16 (link)]. Analyses were performed both with and without EF1-α (m) data to evaluate whether co-regulation with EF1-α biased the results, considering that the EF1-α (m) targets EF1-α as well as other paralogs. Finally, in the third approach, the standard deviation was calculated for each mean Ct difference (Ct-mean Ct) [17 (link)]. β-tubulin, cyclophilin, and EF1-α (m) were excluded from this third approach due to their poor performances during earlier analysis.
To score reference genes based on gene stability, a scoring scheme was implemented whereby results were combined from each statistical approach for each sampling year. Given that cyclophilin, β-tubulin and EF1-α (m) were excluded from the final analyses (Table3 ), they could not be included in the final scoring scheme. For all other genes, the results from each statistical approach (CV, M or 2× SD) within a dataset (2003 or 2004 pericarp development series) were distributed and assigned a score value between 1 to 100, an arbitrary scoring range. For example, the CV values for the 2003 samples ranged from 0.98 for SAND to 2.24 for UBQ-L40. SAND was assigned an arbitrary value of 1 and UBQ-L40, an arbitrary value of 100. Then all other genes were assigned values between 1 and 100, scaled based on their relative distribution. Once all scores were derived for each statistical approach, a cumulative score was used to deduce the final standing.
Data were analyzed using the SDS 1.2.2 software (Applied Biosystems). Expression levels were determined as the number of amplification cycles needed to reach a fixed threshold in the exponential phase of the PCR reaction (Ct). All amplification plots were analyzed with an Rn threshold of 0.2 to obtain Ct values. The PCR efficiency was determined for each gene with LinReg software, which uses absolute fluorescence data captured during the exponential phase of amplification of each reaction [23 (link)]. Results from the SDS and LinReg software were imported into Microsoft Excel for further analyses and to correct for the different PCR efficiencies [24 (link)]. All primer pairs had efficiencies higher than 1.80 with the exception of EF1-α (m) (1.70) and β-tubulin (1.54). Each was run on the full pericarp developmental series for 2003 and 2004. In addition, all primer pairs except those targeting β-tubulin and cyclophilin were run on 2003 mesocarp samples.
In order to evaluate reference gene stability among samples, three statistical approaches were incorporated. In the first approach, Ct difference (Ct max-Ct min) and CV were calculated for each gene throughout each development series tested. During the second approach, Ct values were converted into relative quantities and imported into geNorm v.3.4 software [16 (link)]. Analyses were performed both with and without EF1-α (m) data to evaluate whether co-regulation with EF1-α biased the results, considering that the EF1-α (m) targets EF1-α as well as other paralogs. Finally, in the third approach, the standard deviation was calculated for each mean Ct difference (Ct-mean Ct) [17 (link)]. β-tubulin, cyclophilin, and EF1-α (m) were excluded from this third approach due to their poor performances during earlier analysis.
To score reference genes based on gene stability, a scoring scheme was implemented whereby results were combined from each statistical approach for each sampling year. Given that cyclophilin, β-tubulin and EF1-α (m) were excluded from the final analyses (Table
Full text: Click here