Potential link between biotic defense activation and recalcitrance to induction of somatic embryogenesis in shoot primordia from adult trees of white spruce (Picea glauca)
© Rutledge et al.; licensee BioMed Central Ltd. 2013
Received: 11 March 2013
Accepted: 8 August 2013
Published: 12 August 2013
Among the many commercial opportunities afforded by somatic embryogenesis (SE), it is the ability to clonally propagate individual plants with rare or elite traits that has some of the most significant implications. This is particularly true for many long-lived species, such as conifers, but whose long generation times pose substantive challenges, including increased recalcitrance for SE as plants age. Identification of a clonal line of somatic embryo-derived trees whose shoot primordia have remained responsive to SE induction for over a decade, provided a unique opportunity to examine the molecular aspects underpinning SE within shoot tissues of adult white spruce trees.
Microarray analysis was used to conduct transcriptome-wide expression profiling of shoot explants taken from this responsive genotype following one week of SE induction, which when compared with that of a nonresponsive genotype, led to the identification of four of the most differentially expressed genes within each genotype. Using absolute qPCR to expand the analysis to three weeks of induction revealed that differential expression of all eight candidate genes was maintained to the end of the induction treatment, albeit to differing degrees. Most striking was that both the magnitude and duration of candidate gene expression within the nonresponsive genotype was indicative of an intense physiological response. Examining their putative identities further revealed that all four encoded for proteins with similarity to angiosperm proteins known to play prominent roles in biotic defense, and that their high-level induction over an extended period is consistent with activation of a biotic defense response. In contrast, the more temperate response within the responsive genotype, including induction of a conifer-specific dehydrin, is more consistent with elicitation of an adaptive stress response.
While additional evidence is required to definitively establish an association between SE responsiveness and a specific physiological response, these results suggest that biotic defense activation may be antagonistic, likely related to the massive transcriptional and metabolic reprogramming that it elicits. A major issue for future work will be to determine how and if suppressing biotic defense activation could be used to promote a physiological state more conducive to SE induction.
KeywordsConifer Gene expression profiling Microarray analysis Absolute qPCR LRE qPCR Clonal propagation
Plant somatic embryogenesis (SE) has become an established biotechnology within the horticulture, agriculture and forest industries, providing the capability for commercial-scale production of clonal seedlings [1–3]. However, the efficiency of inducing embryogenic tissue formation (SE induction) continues to be problematic, particularly in woody species such as conifers. For example, although zygotic embryos from a few species belonging to the Pinaceae family are highly responsive, many other conifer species are either completely nonresponsive or produce efficiencies too low to be commercially viable. An even more prominent issue is the recalcitrance of tissues from adult trees, which, if overcome, would allow unlimited propagation of individual trees with elite characteristics .
While judicious manipulation of induction media has found success in improving SE induction efficiency from zygotic embryos, particularly for pines , successful application to vegetative tissues has to date been marginal, at best [4, 5]. In addition, although many physiological and genetic factors impacting SE induction have been documented for angiosperms [6–8], lack of an effective experimental system has impeded efforts to identify even the most fundamental aspects underpinning SE induction within vegetative tissues of conifers.
In an attempt to address this deficiency, experiments initiated over a decade ago targeted somatic embryo-derived white spruce trees with the expectation that they would have a greater propensity for SE induction than trees grown from seed. This led to the identification of a clonal line of white spruce that produced shoot buds that have remained responsive to SE induction even after reaching sexual maturity . Combined with advances in conifer genomics [10–12], this presented an unprecedented opportunity to explore the molecular aspects of SE induction within shoot primordia of adult spruce trees.
Using a recently constructed conifer 32 K oligo-probe microarray , transcriptome-wide expression profiling led to the identification of four of the most differentially expressed genes within this and a nonresponsive genotype at day 7 of induction. Expanding the analysis to day 21 using absolute qPCR revealed substantive differences in the expression dynamics of these candidate genes. Most evident was that both the magnitude and duration of candidate gene expression were greater within the nonresponsive genotype, which is indicative of an intense physiological response to the induction treatment that may be antagonistic to SE induction. Examination of their putative identities further revealed that this intense response may be a result of biotic defense elicitation, whereas the moderate response of the responsive genotype is suggestive of an adaptive response.
Induction of somatic embryogenesis within primordial shoots
Microarray analysis and selection of candidate genes
To further investigate how genotype and the induction treatment interacted, a two-way ANOVA analysis was performed (Figure 2C). This revealed that 8433 targets were differentially expressed across all combinations, with about 90% responding solely in relation to the SE induction treatment. Furthermore, about 10% differed between the two genotypes, and about 3% showed a significant interaction between genotype and treatment. It is important to note that all targets showing a genotype effect also showed response to the SE induction, in which 37% showed a genotype X treatment interaction.
Greatest differential expression within G6 and G12 at day 7 of induction based on microarray analysis
G6 / G12
G12 / G6
G6 candidate genes showing the largest fold differences relative to G12 at day 7 of induction
G6 / G12
G6 / G12
Class III Peroxidase
G12 candidate genes showing the largest fold differences relative to G6 at day 7 of induction
G12 / G6
G12 / G6
Day 7/ 0
Class III Peroxidase
Cell wall invertase
With respect to changes in expression over time, comparing day 7 with day 0 revealed that all but one of the candidate genes increased significantly within both genotypes, with no example of a reduction in gene expression in the apposing genotype. This indicates that differential expression at day 7 was due to higher levels of activation within the originating genotype (Tables 2 and 3). Also notable is that the expression of all but one of the G12 candidates increased to greater levels within the G12 explants than that of the G6 candidates within the G6 explants, suggesting that a major distinguishing characteristic of the nonresponsive G12 genotype is higher levels of candidate gene activation.
Concordance of microarray analysis with absolute qPCR
Conducting absolute quantification greatly increased the resolution of the analysis, in addition to allowing the expression of any gene to be directly compared with that of any other gene, within and between multiple samples. This was accomplished using a method developed by our group called LRE qPCR that greatly simplifies absolute quantification, in large part by abrogating the need to construct target-specific standard curves [16–18].
As is described in the Methods section, expression analysis of nine reference genes revealed that within the five biological replicates taken for each time point used in the microarray analysis, the average variance was found to be about ±20% (intra-sample group variance), which is in part indicative of the analytical precision that can be achieved with LRE qPCR . Furthermore, when their average expression level was compared across the four sample groups, six of nine references generated inter-sample group variances below ±20%, reflective of a remarkably low level of biological variability within this experimental system. This in turn circumvented the need to conduct reference gene normalization as is commonly practiced for qPCR-based gene expression analysis, particularly for those employing relative quantification [19–22].
Average transcript quantities of the candidate genes within the four sample groups used in the microarray analysis
Fold differences in candidate gene expression as determined by microarray and absolute qPCR quantification
Profiling the dynamics of candidate gene expression
The dynamics of DHN1 expression was similar in nature to PgPrx52, peaking at day 7 followed by a progressive 3-fold reduction by day 21 within the G6 explants, indicative of an early, transient-like activation (Figure 4C). However, DHN1 expression was not only apparent within the G12 explants, but progressively increased up to day 15, suggesting that activation of this G6 candidate gene is much less genotype-specific. While differential expression of the proline-rich candidate was maintained up to day 21, both genotypes generated similar expression dynamics, again reflective of modest, if any, genotypic specificity (Figure 4D).
In summary, qPCR expression profiling confirmed the efficacy of candidate gene selection based on microarray analysis at day 7 of induction, in that all of the candidate genes maintained differential expression within their originating genotype throughout the entire induction treatment, although to varying degrees. In addition, the overall biological variability was sufficiently low to reveal strong trends in gene expression dynamics, and to illustrate the utility of absolute quantification, which, among other attributes, provides the ability to assess the magnitude of expression of individual genes. Also consistent with that predicted by the microarray analysis was the high level of induction of all four G12 candidates within the G12 explants, which was maintained well beyond day 7. This suggests that a major distinguishing characteristic of these nonresponsive explants is an intense physiological response to the SE induction treatment.
Recalcitrance of plant explants to many types of tissue culture manipulation, including SE induction, has long been an impediment to clonal propagation of individual plants with elite characteristics, a capability that has significant commercial implications. These also include rare traits, such as insect or pathogen resistance, for which rapid propagation and dissemination could have important ecological implications. This is particularly relevant to forest trees, whose long generation times pose substantive challenges to traditional propagation approaches such as rooting of cuttings, and to genetic improvement programs based on conventional breeding. These limitations are further exacerbated by the fact that many elite traits become most evident in mature plants, a time at which many woody species, including conifers, become recalcitrant to tissue culture manipulation . The identification of a clonal line of adult white spruce trees that produce shoot primordia responsive to SE induction thus presented a unique opportunity to address the molecular aspects of SE induction, with the expectation that this may be broadly applicable to other plant species.
Stress-response as a likely determinant for SE induction responsiveness
As is often the case for SE-induction research, a major presumption entering into this study was that genes selectively expressed within the responsive G6 bud explants would become the primary targets for investigation. Indeed, expression profiling revealed an unusual glutamine-threonine repeat protein (Figure 3B) that was specifically expressed throughout the entire SE induction treatment (Figure 4A). Although such a dramatic differential expression implies a potentially important role in induction responsiveness, this protein appears to be conifer-specific with no known function, excluding the likelihood that this protein plays a universal role in somatic embryogenesis. This is also the case for a proline-rich protein (Figure 3C), which showed marginal, albeit persistent, differential expression in the G6 bud explants up to day 21 (Figure 4D); however, this protein also appears to be conifer-specific with no known function. Although the two remaining candidates encode for proteins with putative functions (dehydrin and peroxidase; Figure 3A and Additional file 3, respectively), the greatest differential expression occurred during early stages of the induction treatment (Figure 4B, C), long before embryogenic tissues begin to emerge, which again does not greatly support a role, at least directly, in determining SE induction responsiveness.
In contrast to this temperate response, the intense candidate gene activation within the nonresponsive G12 bud explants was not only found to persist into the late stages of the induction treatment, but also reached very high levels (maxima 60,000-600,000 transcripts per 10 ng RNA) (Figure 5). While the limited use of absolute quantification has to date provided an inadequate context for universally comparing gene expression levels, EF1α does provide some perspective as to the magnitude of candidate gene expression. Utilized as a reference gene, EF1α expression averaged around 700,000 transcripts per 10 ng RNA (see Methods), which is the most highly expressed gene encountered in this study. This is consistent with it being one of the most highly expressed genes in developing conifer shoots, as based on EST clone frequency and microarray analysis (data not shown). It is thus the magnitude of candidate gene expression within the nonresponsive G12 explants that suggests the physiological response of a primordial shoot explant, particularly during the early stages of SE induction, could determine its ability to become responsive to SE induction.
While the subject area of stress physiology is vast, it is apparent that a number of general principles described for angiosperms could provide insights into the physiological response of conifer bud explants, particularly in relation to abiotic and biotic defense responses. For example, a transient oxidative burst of low intensity is predicted to occur immediately after an explant is introduced into culture, which is indicative of early oxidative signaling generated during both abiotic stress (including wounding) and biotic defense responses [23–25]. A major distinction, however, is that elicitation of a biotic defense response is associated with a second prolonged oxidative burst of high intensity, which in turn triggers massive transcriptional and metabolic reprogramming, including high level induction of defense protein expression, slowing of growth, and, in the most extreme cases, induced cell death [23, 24, 26–29]. This possess the question as to whether the intense response generated by the G12 explants is reflective of a biotic defense response that makes them physiologically recalcitrant to SE induction, in contrast to an adaptive stress response within the G6 explants that generated a physiological state conducive to SE induction.
Although such generalizations provide an attractive model, it should be stressed that angiosperm defense responses have been found to be highly dynamic and complex processes that involve cross-talk between signaling networks regulated by salicylate, jasmonates, and ethylene, in combination with other plant hormones [27–32]. It is therefore difficult to draw specific parallels to conifer bud explants without direct supporting evidence. Nevertheless, examining the putative functions of the proteins encoded by the four G12 candidate genes (an apoplastic class III peroxidase, a cell wall invertase, and two closely related extracellular serine protease inhibitors) provides support for the contention that the SE induction treatment elicited a biotic defense response within these nonresponsive explants.
Apoplastic class III peroxidases play a prominent role in biotic defense activation
Class III peroxidases have been found to generate apoplastic hydrogen peroxide that acts as a signal for biotic defense elicitation [27, 33–36]. Direct demonstration of this was recently reported for maize, in which U. maydis (corn smut) leaf infection was found to induce the expression of a single class III peroxidase gene (POX12), and that virulence of this pathogen is dependent on its excretion of a peroxidase inhibitor peptide (Pep 1) [37, 38]. Furthermore, induced gene silencing of POX 12 was found to restore virulence even in the absence of Pep 1, providing direct evidence that POX12 activity is essential for mounting a successful biotic defense response against this pathogen . The progressive, high level activation of a class III peroxidase within the G12 explants (Figure 5A) is thus consistent with a role in biotic defense elicitation.
Additional, albeit tentative, support for such a contention is that this G12 candidate is most similar in amino acid sequence to the Arabidopsis class III peroxidase, AtPrx21 (Additional file 3), which belongs to an unusual evolutionary branch of plant peroxidases . Increased expression of AtPrx21 produced by wounding and microbial attack has led to the suggestion that it has a protective role against pathogens . A more direct demonstration of a role in biotic defense comes from the fact that overexpression of AtPrx21 in Arabidopsis produces resistance to Botrytis cinerea. Thus, currently available data suggests that persistent, high level expression of an apoplastic class III peroxidase gene plays a central role in biotic defense activation, a role that the G12 PgPrx21 may play in conifers.
Conversely, the transient nature of the G6 class III peroxidase PgPrx52 gene activation (Figure 4B) could be reflective of an adaptive stress response in which the initial oxidative burst dissipates during the first few days of induction treatment, restoring cellular redox homeostasis [41–43]. This is similar to the activation, although slightly later, of the G6 dehydrin DHN1 gene (Figure 4C), which is a conifer-specific dehydrin that has been reported to play a role in bud dormancy and overwintering in Norway spruce . This too supports the contention that the G6 explants elicited an adaptive stress response, in that dehydrins have long been recognized as playing a fundamental role in adapting to environmental stresses . Thus, although speculative, it could be argued that an adaptive stress response may be an important determinant for establishing SE induction responsiveness.
Induction of a cell wall invertase
Activation of cell wall invertases, which catalyze the hydrolysis of sucrose to glucose and fructose within the apoplast, have also been found to play a prominent role in biotic defense by providing the large quantities of energy required for mounting an intense metabolic response [45–49]. RNAi inhibition of expression has provided direct support for a central role of a cell wall invertase in defense response elicitation in tobacco . Indeed, the prominence of cell wall invertase induction has led to a proposal that it be classified as a pathogenesis-related protein , which is a large class of defense proteins that become highly expressed following biotic defense activation .
Amino acid sequence comparison to Arabidopsis revealed that the G12 invertase is most similar to AtcwINV1 (also called AtβFruct1) (Additional file 3), which is a member of a small gene family within Arabidopsis . A number of studies have reported induction of AtcwINV1 expression following pathogen infection of Arabidopsis leaves [52–54], in addition to induction following Agrobacterium infection of Arabidopsis cell cultures . Importantly, selective induction of AtcwINV1 following wounding also suggests a general role in adaptive stress responses . This is consistent with the progressive increase in expression of PgcwINV1 observed within the G6 explants. However, this is contrasted by a rapid activation within the G12 explants between day 3 and 7, reaching an average of about twice that found within the G6 explants, a level that is maintained up to day 21 of the induction treatment (Figure 5B). This could be reflective of a more intense metabolic response of the G12 explants during the earliest stages of the induction treatment, a presumption consistent with the high levels of expression observed for the two remaining G12 candidate genes.
Activation of two I20 serine protease inhibitors
The two most highly expressed G12 candidate genes were found to encode for closely related protease inhibitors belonging to an unassigned subclass of the MEROPS I20 family of serine protease inhibitors [57, 58]. This subclass has a number of distinguishing features, including the presence of a transient peptide, an unusually small mature protein (less than 60 amino acids), and the presence of eight highly conserved cysteines (Figure 3D). In fact, these features have led to their classification into the superfamily of small cysteine-rich peptides (CRP), a very large family of secreted peptides composed of several hundred genes within Arabidopsis. Initially founded on structural similarities with defensins, which is an ancient form of antimicrobial peptide , a large number of CRPs have been shown to play a role in biotic defense . Another distinctive characteristic of this subclass of I20 protease inhibitors is the occurrence of closely related homologues throughout the Angiospermae, principally as a single gene . In fact, this high level of conservation led Hartl et al.  to suggest that in addition to a probable role in biotic defense, they could have an essential role in plant physiology. The presence of two highly conserved homologs within conifers provides support for such a possibility.
A direct link of the angiosperm homolog to biotic defense activation comes from expression analysis of U. maydis Δpep1-infected maize leaves, which elicits a massive biotic defense response, as was described earlier. Among the 220 genes found to be induced, the maize I20 homolog (EF406275 in Figure 3D) was the fourth most highly induced, increasing by 166-fold within 24 h post-inoculation (Additional file 1: Table S1 in ). Similarly, the Arabidopsis I20 homolog (At1G72060) was the fifth most induced gene following infection by Trichoderma harzianum, reaching a 2.79-fold increase within 24 h . Application of biotic defense elicitors to 10-day-old Arabidopsis seedlings was also found to induce this I20 gene, with the bacterial flagellin peptide Flg22 generating 5.08- and 7.45-fold increases after 1 and 3 h, respectively (Additional file 1: Table S1 in ). Oxidative stress induction of the maize homolog was also demonstrated in maize leaves (referred to as serpin), reaching 2.91- and 7.56-fold increases 16 h after application of H2O2 and methyl viologen, respectively . However, a role in abiotic defense is unlikely, as ectopic expression of the Arabidopsis homolog in transgenic seedlings dramatically reduced resistance to oxidative, osmotic and salt stress .
The induction kinetics of the two conifer I20 homologs within the G12 explants is thus consistent with biotic defense activation as observed in angiosperms, with both genes reaching close to 200,000 transcripts per 10 ng RNA by day 7, which, in comparison with day 0, roughly represents a 2000-fold increase. Another notable feature was a progressive increase in expression throughout the entire induction treatment, reaching about 600,000 transcripts per 10 ng RNA by day 21 (Figure 5C and D). Thus, while the precise biochemical function of this unusual class of serine protease inhibitors remains to be determined, the high level of amino acid sequence conservation, combined with their expression dynamics, provide support to the supposition that these two conifer protease inhibitors play a role in conifer biotic defense, similar to that observed in angiosperms.
A paradigm shift towards physiological processes that may antagonize SE induction
While it has long been recognized that the physiological state of an explant can be a major determinant for responsiveness to SE induction, very little is understood about the underlying mechanisms. Furthermore, the vast majority of research efforts have historically focused on defining factors that promote SE induction, for example, through the activation of SE-promoting genes. This study expands this paradigm, suggesting that nonresponsiveness of an explant is not necessarily due to an innate lack of SE promoting activity, but that biotic defense activation could potentially be a dominant antagonist.
A number of physiological aspects of biotic defense elicitation provide general support for such a presumption. For example, activation of biotic defense through exogenous application of various elicitors has been shown to dramatically reduce plant growth, a result of redirecting metabolic energy from growth to defense [65, 66]. This includes the action of jasmonates (JAs), the central regulators of biotic defense elicitation associated with wounding, which have been shown to directly mediate a switch from growth to the production of biocidal compounds, cell-wall remodeling, and defense protein expression . Indeed, growth repression by JAs has been directly linked to inhibition of cell-cycle progression , in addition to directly antagonizing the growth promoting activity of gibberellic acid , thus providing evidence that JAs could be direct antagonists of SE induction, based in part on the assumption that cell division is necessary for embryogenic tissue formation.
Another line of supporting evidence, albeit indirect, comes from proteomic studies that have reported a correlation between expression of biotic defense proteins, primarily pathogenesis-related proteins, and a lack of embryogenic competency of tissues in culture [68–71]. Another notable observation related to the quantitative nature of proteomic analysis, is the magnitude of defense protein expression within these nonembryogenic tissues, often being the most prominent proteins in the analysis. While it is difficult to draw a direct comparison, this is consistent with the intense activation of the G12 candidate genes within the G12 explants, and with the supposition that redirecting metabolic resources towards the production of such large quantities of defense proteins could itself be antagonistic to the formation of embryogenic tissues.
SE induction within leaf explants from the model legume Medicago truncatula, an experimental system similar to the bud explants used in this study, has also been used to directly compare responsive and nonresponsive genotypes . Proteomic analysis revealed large physiological differences as reflected by high levels of protein accumulation, some of which were identified as stress proteins, within a nonresponsive line during the first week of the SE induction treatment. Although an association with biotic defense elicitation was not evident from the data presented, this led the authors to suggest that a hyperresponse to the stress produced by the induction treatment could be related to a lack of responsiveness . This is a scenario similar to the intense physiological response of the G12 explants, and reminiscent of the high levels of defense protein accumulation within nonembryogenic tissues, as discussed above.
Finally, although this study provides evidence for a link between biotic defense elicitation and recalcitrance to SE induction, it has also generated hypotheses that could not be directly addressed within the experimental design. First, only one nonresponsive genotype was analyzed, so the question as to whether a similar physiological response occurs in other nonresponsive genotypes remains unanswered. Second, the hypothesis of a direct association between biotic defense elicitation and suppression of SE induction remains to be tested. Thus, it is possible that the differences in physiological response observed within the G6 and G12 explants is a result of a genotypic difference unrelated to SE responsiveness. Third, it is unknown whether SE could be induced in the G12 explants even if biotic defense elicitation were to be mitigated. Additional work is thus required before a definitive understanding of the broad applicability and implications of these findings can be achieved.
The central conclusion of this study is that the physiological response of conifer bud explants, particularly in relation to elicitation of a defense response, could be an important determinant of SE induction responsiveness. Although definitive demonstration that biotic defense activation is antagonistic to SE induction requires additional evidence, many general characteristics, such as the dramatic metabolic and transcriptional reprogramming associated with its elicitation, support a role contributing to the recalcitrance of explants to SE induction. In addition, it opens new avenues of investigation into the mechanisms regulating the activation and intensity of defense responses within explants placed into culture, along with the prospective of developing methods that could be used to suppress them, with the expectation that this could generate a physiological state more conducive not only to SE induction, but potentially to other types of tissue culture manipulation.
Primordial shoot collection and somatic embryogenesis induction
Shoot buds were collected on May 4 and 6, 2009, from the second and third whorls of branches of 9-year-old Picea glauca (white spruce) trees that were generated from somatic embryos as previously described . These consisted of a responsive (893-6: G6) and nonresponsive (893-12: G12) genotype, from which a total of 700 shoot buds were collected from several clonal trees per genotype. SE induction was conducted as previously described . Briefly, the buds were disinfected, primordial shoots were excised and cut longitudinally into two or four equal parts, and the explants were placed onto semi-solid MLV-S medium containing 9.5 μM 2,4-dichlorophenoxyacetic acid and 4.5 μM benzyl adenine. Replicate explant samples were collected after 3, 7, 15 and 21 days of culture, frozen in liquid nitrogen and stored at −80°C. All remaining explants were cultured for 16 weeks, during which formation of embryonal masses was verified under a stereomicroscope and recorded.
Two rounds of RNA extraction were conducted, referred to as sample series 1 and 2. For sample series 1, which was used in the microarray analysis, a CTAB-LiCl precipitation protocol  was used to extract RNA from five biological replicates of both genotypes, consisting of approximately 80 mg of buds that were either placed into liquid nitrogen immediately after collection in the field (day 0), or after one week of SE induction treatment (day 7). For qPCR analysis, aliquots of these RNA samples were DNase treated before cDNA production, as described in the reverse transcription section.
For sample series 2, three replicate collections were taken at day 3, 15 and 21 of induction, with each replicate consisting of approximately 80 mg of fresh mass. These were placed into a 2-ml Sarstedt conical microtube containing a single 5 mm stainless steel bead, frozen in liquid nitrogen, and stored at −80°C. However, it should be noted that subsequent work has revealed that smaller amounts (less than 50 mg) can significantly improve both the quality and quantity of RNA recovered for some sample types. The tubes were transferred into an adapter set that was prechilled at −80°C for a minimum of two hours and transported in a cooler containing a few inches of liquid nitrogen in order to prevent the samples from thawing. The tissues were disrupted twice for 45 s at 26 Hz using the Qiagen Tissuelyzer II bead mill. The adapter set was returned to the cooler, and each tube was removed one at a time. In each tube 550 μl of lysis buffer (4 M guanidine isothiocyanate, 0.2 M sodium acetate, pH 5.0, 25 mM EDTA, 2.5% (wt/vol) PVP-40) was added . The tubes were vortexed at high speed and incubated for 2 min at 56°C, during which one or two more vortexing steps were conducted. Samples were then centrifuged briefly to remove cell debris and 450 μl were removed for RNA extraction using the Qiagen RNEasy plant mini kit (Cat. # 74904). RNA extractions were performed using a Qiacube DNA/RNA purification robot (Qiagen), which included an on column DNase treatment (Qiagen RNase free DNase, Cat. # 79254). However, variable quantities of genomic DNA were subsequently detected by qPCR, such that a second DNase treatment was necessary (see the reverse transcription section for details). Following this DNase treatment, RNA was quantified using a Nanodrop 1000 Spectrophotometer (Thermo Scientific), and RNA integrity was assessed using the Agilent 2100 Bioanalyzer (Agilent Technologies), which generated RIN values of 8.2-9.8.
Microarray experiment and analysis were conducted as previously described , with minor modifications. Briefly, 1 μg of each total RNA sample was amplified using the Amino Allyl MessageAmpII aRNA Amplification Kit (Life Technologies, Carlsbad, CA, USA), fragmented and quantified, and 5 μg of amplified RNA was labeled with AlexaFluor 647 (Life Technologies). Prehybridization of the oligonucleotide arrays, hybridization of the labeled samples to the slides, slide washing and drying were performed on HS400Pro hybridization stations (Tecan Group Ltd., Männedorf, Switzerland). Slide scanning and feature extraction were done on a ScanArray Express scanner (PerkinElmer, Waltham, MA, USA) and QuantArray v3.0 (PerkinElmer), respectively. Each array image was analyzed in two sections (top and bottom half), and both sections were fused into one file in Excel using an in-house macro. All of the arrays (20) were kept for further analysis. The experiment was analyzed as a one-color design with four groups of five arrays (two time points and two genotypes, 5 biological replicates). Quality control and data processing, namely background and buffer subtraction, aQuantile normalization and correction for multiple testing (Benjamini-Hotchberg), were done in R version 2.8.1 . All microarray data have been deposited in NCBI’s Gene Expression Omnibus and are accessible through GEO Series accession number GSE46977. Genes were selected as being differentially expressed on the basis of their adjusted p-values (p-value < 0.05). Analysis of variance and t-tests of selected genes were conducted using Flexarray v1.6 .
Candidate target selection and primer design
Candidate targets were selected based on fold differences in expression between the two genotypes at day 7 of SE induction with no consideration of target gene identity (Tables 2, 3; Additional file 1). A secondary consideration was similar expression within the two genotypes at the point of bud collection (day 0), which led to the rejection of one G12 and six G6 targets. The probe sequence from each candidate target was blasted against the NCBI Picea EST libraries and nucleotide collection databases, from which nucleotide alignments were constructed for primer selection (Additional files 4 and 5), which was based primarily on positioning the amplicon close to the stop codon in order to minimize variance due to partial reverse transcription.
The central parameter used for primer design was determining a length sufficient to generate a predicted Tm of 70°C using the Integrated DNA Technologies online OligoAnalyzer. This was based on calibrating the program by arbitrarily adjusting the Mg++ parameter concentration to 50 mM, such that the predicted Tm of CAL1 F1 and CAL1 R1 primers (Additional file 5) reached 70°C, an approach found successful for other oligo design programs. Following calibration, candidate primers were then designed by simply adjusting their length until the predicted Tm just exceeded 70°C. Extensive self-complementary primers or those complementary to the apposing primer were rejected, as were any primer pair that generated non-specific products in a no template control amplification or that generated amplification efficiencies <99%. Amplicon size was restricted to 80–200 bp; however, extensive efforts to predict primer performance, such as analyzing the secondary structure of the primers or of the resulting amplicon, were unsuccessful. It was therefore necessary to test multiple primer pair combinations for some targets, rejecting those that generated profile collapse or extensive plateau drifting as indicated by LRE analysis (see  for details about these anomalies).
Before conducting reverse transcription, genomic DNA contamination was quantified by amplifying 20 ng samples of raw RNA. This revealed that many samples contained small amounts of gDNA contamination (10–100 genomes) so that all RNA samples were DNase-treated using the Ambion Turbo DNA-free DNase kit (Cat. # AM1907), which was found to reduce gDNA contamination to undetectable levels.
Reverse transcription was conducted in 20 μl reactions containing 50 ng/μl RNA, 25 ng/μl oligo dT primer (Invitrogen, Cat.# 18418–012), 5 U/μl Superscript II (Invitrogen, Cat. # 18064–014) using the manufacturer’s supplied buffer, and incubated at 42°C for 50 min, followed by the addition of 180 μl 10 mM Tris, pH 8.0 to generate a final concentration of 5 ng RNA per μl. Testing revealed that the commonly practiced addition of RNase inhibitor and DTT provided no benefit. Note also that we have found that RNase H treatment produces extensive scattering of replicate profiles for some cDNA targets, which can greatly reduce quantitative accuracy. As previously reported, three replicate reverse transcriptase reactions using an identical RNA sample, generated an average variance of about ±12% for three reference gene targets , indicating that this method is highly repeatable, consistent with the small variances in reference gene expression observed in this study (see below).
PCR amplification was conducted with an Applied Biosystems 7500 Fast qPCR system (normal ramping), QuantiTect enzyme formulation (Qiagen, Cat. # 204145) in a 10 μl reaction volume containing 500 nM of primers, an aliquot of reverse transcriptase reaction equivalent to 5 ng RNA, using 96 well BrightWhite plates (Primerdesign, BW-Fast) sealed with MicroAmp film (Applied Biosystems, Cat # 4311971). The cycling regime consisted of a 15 min activation at 95°C, followed by 50 cycles of 95°C for 10 s, and 65°C for 120 s. Amplicon Tm was determined for each amplification reaction by melt curve analysis (65 to 90°C) conducted at the end of each run. Raw fluorescence readings were imported into the LRE Analyzer (Version 0.8.7), an open source Java program available for download [79, 80]. The program contains an extensive help dataset that describes how the program functions, along with various guidelines for setting up and testing the performance of qPCR assays.
All primer pairs produced a single prominent amplicon peak following melt curve analysis, and generated amplification efficiencies that ranged from 102 to 107%. The propensity for QuantiTect to generate amplification efficiencies greater than 100%, a systemic bias that generates moderate underestimates of target quantity, was compensated by fixing the amplification efficiency to 100%, as described in the LRE Analyzer help. The LRE Analyzer databases are provided in Additional file 6, along with the amplicon and optical calibration databases used in this study. A summary of the individual quantitative determinations is provided in Additional file 7.
Reference gene analysis
Expression of nine reference genes was used to determine the levels of biological variability, in addition to serving as internal quality controls for assessing the technical variance associated with sample preparation and LRE qPCR analysis. Two reference genes were taken from microarray analysis of Sitka spruce apical shoots (hypothetical protein (HP) and peroxisomal targeting signal receptor (PTSR)) , with the remaining seven being conifer homologs to Arabidopsis reference genes also identified from microarray analysis . Primer sequences along with UniGene accession numbers are provided in Additional file 5.
The authors thank A. Séguin for critically reviewing the manuscript and I. Lamarre for editorial assistance. Funding for this work was provided by Natural Resources Canada (RGR, KK) and from Genome Canada, Genome Quebec, Genome BC for the Arborea, and SMarTForests Projects (JM), and the Discovery grant program of the Natural Sciences and Engineering Council (JM).
- Preil W: Plant tissue culture, 100 years since Gottlieb Haberlandt. Micropropagation of ornamental plants. Edited by: Laimer M, Rücker W, Wien W. 2003, New York: Springer, 115-133.Google Scholar
- Nehra NS, Becwar MR, Rottmann WH, Pearson L, Chowdhury K, Chang S, Dayton Wilde H, Kodrzycki RJ, Zhang C, Gause KC, Parks DW, Hinchee M: Vitro Cell Dev Biol-Plant. Forest biotechnology: innovative methods, emerging opportunities 41st edition. 2005, 701-717.Google Scholar
- Klimaszewska K, Trontin JF, Becwar M, Devillard C, Park YS, Lelu-Walter MA: Recent progress in somatic embryogenesis of four Pinus spp. Tree For Sci Biotechnol. 2007, 1: 11-25.Google Scholar
- Bonga JM, Klimaszewska K, Aderkas P: Recalcitrance in clonal propagation, in particular of conifers. Plant Cell Tissue and Organ Cult. 2010, 100: 241-254. 10.1007/s11240-009-9647-2.View ArticleGoogle Scholar
- Park SY, Klimaszewska K, Park JY, Mansfield SD: Lodgepole pine: the first evidence of seed-based somatic embryogenesis and the expression of embryogenesis marker genes in shoot bud cultures of adult trees. Tree Physiol. 2010, 30: 1469-1478. 10.1093/treephys/tpq081.PubMedView ArticleGoogle Scholar
- Verdeil JL, Alemanno L, Niemenak N, Tranbarger TJ: Pluripotent versus totipotent plant stem cells: dependence versus autonomy?. Trends Plant Sci. 2007, 12: 245-252. 10.1016/j.tplants.2007.04.002.PubMedView ArticleGoogle Scholar
- Zhang H, Ogas J: An epigenetic perspective on developmental regulation of seed genes. Mol Plant. 2009, 2: 610-627. 10.1093/mp/ssp027.PubMedView ArticleGoogle Scholar
- Grafi G, Florentin A, Ransbotyn V, Morgenstern Y: The stem cell state in plant development and in response to stress. Front Plant Sci. 2011, 2: 1-10.View ArticleGoogle Scholar
- Klimaszewska K, Overton C, Stewart D, Rutledge RG: Initiation of somatic embryos and regeneration of plants from primordial shoots of 10-year-old somatic white spruce and expression profiles of 11 genes followed during the tissue culture process. Planta. 2011, 233: 635-647. 10.1007/s00425-010-1325-4.PubMedView ArticleGoogle Scholar
- Rigault P, Boyle B, Lepage P, Cooke JEK, Bousquet J, MacKay JJ: A white spruce gene catalog for conifer genome analyses. Plant Physiol. 2011, 157: 14-28. 10.1104/pp.111.179663.PubMedPubMed CentralView ArticleGoogle Scholar
- MacKay JJ, Dean JFD, Plomion C, Peterson DG, Cánovas FM, Pavy N, Ingvarsson PK, Savolainen O, Guevara MA, Fluch S, Vinceti B, Abarca D, Díaz-Sala C, Cervera MT: Towards decoding the conifer giga-genome. Plant Mol Biol. 2012, 80 (6): 555-569. 10.1007/s11103-012-9961-7.PubMedView ArticleGoogle Scholar
- Raherison ES, Rigault P, Caron S, Poulin PL, Boyle B, Verta JP, Giguère I, Bomal C, Bohlmann J, MacKay JJ: Transcriptome profiling in conifers and the PiceaGenExpress database show patterns of diversification within gene families and interspecific conservation in vascular gene expression. BMC Genomics. 2012, 13: 434-10.1186/1471-2164-13-434.PubMedPubMed CentralView ArticleGoogle Scholar
- Yakovlev IA, Asante DKA, Gunnar C, Jouni F, Olavi P: Dehydrins expression related to timing of bud burst in Norway spruce late flushing family. Planta. 2008, 228: 459-472. 10.1007/s00425-008-0750-0.PubMedView ArticleGoogle Scholar
- Petersen TN, Brunak S, Von Heijne G, Nielsen H: SignalP 4.0: discriminating signal peptides from transmembrane regions. Nat Meth. 2011, 8: 785-786. 10.1038/nmeth.1701.View ArticleGoogle Scholar
- Silverstein KT, Moskal W, Wu HC, Underwood B, Graham M, Town CD, VandenBosch K: Small cysteine-rich peptides resembling antimicrobial peptides have been under-predicted in plants. Plant J. 2007, 51: 262-280. 10.1111/j.1365-313X.2007.03136.x.PubMedView ArticleGoogle Scholar
- Rutledge RG, Stewart D: A kinetic-based sigmoidal model for the polymerase chain reaction and its application to high-capacity absolute quantitative real-time PCR. BMC Biotechnol. 2008, 8: 47-10.1186/1472-6750-8-47.PubMedPubMed CentralView ArticleGoogle Scholar
- Rutledge RG, Stewart D: Assessing the performance capabilities of LRE-based assays for absolute quantitative real-time PCR. PLoS ONE. 2010, 5: e9731-10.1371/journal.pone.0009731.PubMedPubMed CentralView ArticleGoogle Scholar
- Rutledge RG: A Java program for LRE-based real-time qPCR that enables large-scale absolute quantification. PLoS ONE. 2011, 6: e17636-10.1371/journal.pone.0017636.PubMedPubMed CentralView ArticleGoogle Scholar
- Vandesompele J, Kubista M, Pfaffl M: Real-time PCR: current technology and applications. Reference gene validation software for improved normalization 4th edition. Edited by: Logan J, Edwards K, Saunders NB. 2009, Norfolk, UK: Caister Academic Press, 47-64.Google Scholar
- Schmittgen TD, Livak KJ: Analyzing real-time PCR data by the comparative CT method. Nat Protocols. 2008, 3: 1101-1108. 10.1038/nprot.2008.73.PubMedView ArticleGoogle Scholar
- Huggett JF, Dheda K, Bustin SA: Real-time RT-PCR normalisation; strategies and considerations. Genes Immun. 2005, 6: 279-284. 10.1038/sj.gene.6364190.PubMedView ArticleGoogle Scholar
- Gutierrez L, Mauriat M, Guénin S, Pelloux J, Lefebvre JF, Louvet R, Rusterucci C, Moritz T, Guerineau F, Bellini C, Van Wuytswinkel O: The lack of a systematic validation of reference genes: a serious pitfall undervalued in reverse transcription-polymerase chain reaction (RT-PCR) analysis in plants. Plant Biotechnol J. 2008, 6: 609-618. 10.1111/j.1467-7652.2008.00346.x.PubMedView ArticleGoogle Scholar
- Lamb C, Dixon RA: The oxidative burst in plant disease resistance. Annu Rev Plant Physiol Plant Mol Biol. 1997, 48: 251-275. 10.1146/annurev.arplant.48.1.251.PubMedView ArticleGoogle Scholar
- Alvarez M, Pennell RI, Meijer P-j, Ishikawa A, Dixon RA, Lamb C: Reactive oxygen intermediates mediate a systemic signal network in the establishment of plant immunity. Cell. 1998, 92: 773-784. 10.1016/S0092-8674(00)81405-1.PubMedView ArticleGoogle Scholar
- Marino D, Dunand C, Puppo A, Pauly N: A burst of plant NADPH oxidases. Trends Plant Sci. 2012, 17: 9-15. 10.1016/j.tplants.2011.10.001.PubMedView ArticleGoogle Scholar
- Cheong YH, Chang HS, Gupta R, Wang X, Zhu T, Luan S: Transcriptional profiling reveals novel interactions between wounding, pathogen, abiotic stress, and hormonal responses in Arabidopsis. Plant Physiol. 2002, 129: 661-677. 10.1104/pp.002857.PubMedPubMed CentralView ArticleGoogle Scholar
- Torres MA, Jones JDG, Dangl JL: Reactive oxygen species signaling in response to pathogens. Plant Physiol. 2006, 141: 373-378. 10.1104/pp.106.079467.PubMedPubMed CentralView ArticleGoogle Scholar
- Pauwels L, Inzé D, Goossens A: Jasmonate-inducible gene: what does it mean?. Trends Plant Sci. 2009, 14: 87-91. 10.1016/j.tplants.2008.11.005.PubMedView ArticleGoogle Scholar
- Chaouch S, Queval G, Noctor G: AtRbohF is a crucial modulator of defence-associated metabolism and a key actor in the interplay between intracellular oxidative stress and pathogenesis responses in Arabidopsis. Plant J. 2012, 69: 613-627. 10.1111/j.1365-313X.2011.04816.x.PubMedView ArticleGoogle Scholar
- Robert-Seilaniantz A, Grant MR, Jones JDG: Hormone crosstalk in plant disease and defense: more than just jasmonate-salicylate antagonism. Ann Rev Phytopathol. 2011, 49: 317-343. 10.1146/annurev-phyto-073009-114447.View ArticleGoogle Scholar
- Kazan K, Manners JM: JAZ repressors and the orchestration of phytohormone crosstalk. Trends Plant Sci. 2012, 17: 22-31. 10.1016/j.tplants.2011.10.006.PubMedView ArticleGoogle Scholar
- Thaler JS, Humphrey PT, Whiteman NK: Evolution of jasmonate and salicylate signal crosstalk. Trends Plant Sci. 2012, 17: 260-270. 10.1016/j.tplants.2012.02.010.PubMedView ArticleGoogle Scholar
- Apel K, Hirt H: Reactive oxygen species: metabolism, oxidative stress, and signal transduction. Ann Rev Plant Biol. 2004, 55: 373-399. 10.1146/annurev.arplant.55.031903.141701.View ArticleGoogle Scholar
- Passardi F, Cosio C, Penel C, Dunand C: Peroxidases have more functions than a Swiss army knife. Plant Cell Rep. 2005, 24: 255-265. 10.1007/s00299-005-0972-6.PubMedView ArticleGoogle Scholar
- Almagro L, Gómez Ros LV, Belchi-Navarro S, Bru R, Ros Barceló A, Pedreño M: Class III peroxidases in plant defence reactions. J Exp Bot. 2009, 60: 377-390. 10.1093/jxb/ern277.PubMedView ArticleGoogle Scholar
- Cosio C, Dunand C: Specific functions of individual class III peroxidase genes. J Exp Bot. 2009, 60: 391-408. 10.1093/jxb/ern318.PubMedView ArticleGoogle Scholar
- Doehlemann G, Linde KVD, Aßmann D, Schwammbach D, Mohanty A, Jackson D, Kahmann R: Pep1, a secreted effector protein of Ustilago maydis, is required for successful invasion of plant cells. PLoS Pathog. 2009, 5: e1000290-10.1371/journal.ppat.1000290.PubMedPubMed CentralView ArticleGoogle Scholar
- Hemetsberger C, Herrberger C, Zechmann B, Hillmer M, Doehlemann G: The Ustilago maydis effector Pep1 suppresses plant immunity by inhibition of host peroxidase activity. PLoS Pathog. 2012, 8: e1002684-10.1371/journal.ppat.1002684.PubMedPubMed CentralView ArticleGoogle Scholar
- Kjaersgård IV, Jespersen HM, Rasmussen SK, Welinder KG: Sequence and RT-PCR expression analysis of two peroxidases from Arabidopsis thaliana belonging to a novel evolutionary branch of plant peroxidases. Plant Mol Biol. 1997, 33: 699-708. 10.1023/A:1005707813801.PubMedView ArticleGoogle Scholar
- Chassot C, Nawrath C, Métraux JP: Cuticular defects lead to full immunity to a major plant pathogen. Plant J. 2007, 49: 972-980. 10.1111/j.1365-313X.2006.03017.x.PubMedView ArticleGoogle Scholar
- Potters G, Horemans N, Jansen MAK: The cellular redox state in plant stress biology–a charging concept. Plant Physiol Biochem. 2010, 48: 292-300. 10.1016/j.plaphy.2009.12.007.PubMedView ArticleGoogle Scholar
- Mittler R, Vanderauwera S, Suzuki N, Miller G, Tognetti VB, Vandepoele K, Gollery M, Shulaev V, Van Breusegem F: ROS signaling: the new wave?. Trends Plant Sci. 2011, 16: 300-309. 10.1016/j.tplants.2011.03.007.PubMedView ArticleGoogle Scholar
- Foyer CH, Noctor G: Ascorbate and glutathione: the heart of the redox hub. Plant Physiol. 2011, 155: 2-18. 10.1104/pp.110.167569.PubMedPubMed CentralView ArticleGoogle Scholar
- Hanin M, Brini F, Ebel C, Toda Y, Takeda S, Masmoudi K: Plant dehydrins and stress tolerance: versatile proteins for complex mechanisms. Plant Signal Behav. 2011, 6: 1503-1509. 10.4161/psb.6.10.17088.PubMedPubMed CentralView ArticleGoogle Scholar
- Roitsch T, Balibrea ME, Hofmann M, Proels R, Sinha AK: Extracellular invertase: key metabolic enzyme and PR protein. J Exp Bot. 2003, 54: 513-524. 10.1093/jxb/erg050.PubMedView ArticleGoogle Scholar
- Roitsch T, González MC: Function and regulation of plant invertases: sweet sensations. Trends Plant Sci. 2004, 9: 606-613. 10.1016/j.tplants.2004.10.009.PubMedView ArticleGoogle Scholar
- Fotopoulos V: Plant invertases: structure, function and regulation. J Biol Res. 2005, 4: 127-137.Google Scholar
- Berger S, Sinha AK, Roitsch T: Plant physiology meets phytopathology: plant primary metabolism and plant-pathogen interactions. J Exp Bot. 2007, 58: 4019-4026. 10.1093/jxb/erm298.PubMedView ArticleGoogle Scholar
- Bolton MD: Primary metabolism and plant defense–fuel for the fire. Mol Plant-Microbe Interact. 2009, 22: 487-497. 10.1094/MPMI-22-5-0487.PubMedView ArticleGoogle Scholar
- Essmann J, Schmitz-Thom I, Schön H, Sonnewald S, Weis E, Scharte J: RNA interference-mediated repression of cell wall invertase impairs defense in source leaves of tobacco. Plant Physiol. 2008, 147: 1288-1299. 10.1104/pp.108.121418.PubMedPubMed CentralView ArticleGoogle Scholar
- Van Loon LC, Rep M, Pieterse CMJ: Significance of inducible defense-related proteins in infected plants. Annu Rev Phytopathol. 2006, 44: 135-162. 10.1146/annurev.phyto.44.070505.143425.PubMedView ArticleGoogle Scholar
- Chou HM, Bundock N, Rolfe SA, Scholes JD: Infection of Arabidopsis thaliana leaves with Albugo candida (white blister rust) causes a reprogramming of host metabolism. Mol Plant Pathol. 2000, 1: 99-113. 10.1046/j.1364-3703.2000.00013.x.PubMedView ArticleGoogle Scholar
- Fotopoulos V, Gilbert MJ, Pittman JK, Marvier AC, Buchanan AJ, Sauer N, Hall JL, Williams LE: The monosaccharide transporter gene, AtSTP4, and the cell-wall invertase, Atbetafruct1, are induced in Arabidopsis during infection with the fungal biotroph Erysiphe cichoracearum. Plant Physiol. 2003, 132: 821-829. 10.1104/pp.103.021428.PubMedPubMed CentralView ArticleGoogle Scholar
- Bonfig KB, Schreiber U, Gabler A, Roitsch T, Berger S: Infection with virulent and avirulent P. syringae strains differentially affects photosynthesis and sink metabolism in Arabidopsis leaves. Planta. 2006, 225: 1-12. 10.1007/s00425-006-0303-3.PubMedView ArticleGoogle Scholar
- Ditt RF, Kerr KF, De Figueiredo P, Delrow J, Comai L, Nester EW: The Arabidopsis thaliana transcriptome in response to Agrobacterium tumefaciens. Mol Plant-Microbe Interact. 2006, 19: 665-681. 10.1094/MPMI-19-0665.PubMedView ArticleGoogle Scholar
- Quilliam RS, Swarbrick PJ, Scholes JD, Rolfe SA: Imaging photosynthesis in wounded leaves of Arabidopsis thaliana. J Exp Bot. 2006, 57: 55-69.PubMedView ArticleGoogle Scholar
- Rawlings ND, Tolle DP, Barrett AJ: Evolutionary families of peptidase inhibitors. Biochem J. 2004, 378: 705-716. 10.1042/BJ20031825.PubMedPubMed CentralView ArticleGoogle Scholar
- Rawlings ND, Barrett AJ, Bateman A: MEROPS: the database of proteolytic enzymes, their substrates and inhibitors. Nucleic Acids Res. 2012, 40: D343-D350. 10.1093/nar/gkr987.PubMedPubMed CentralView ArticleGoogle Scholar
- Silverstein KAT, Graham MA, Paape TD, VandenBosch KA: Genome organization of more than 300 defensin-like genes in Arabidopsis. Plant Physiol. 2005, 138: 600-610. 10.1104/pp.105.060079.PubMedPubMed CentralView ArticleGoogle Scholar
- Hartl M, Giri AP, Kaur H, Baldwin IT: The multiple functions of plant serine protease inhibitors: defense against herbivores and beyond. Plant Signal Behav. 2011, 6: 1009-1011. 10.4161/psb.6.7.15504.PubMedPubMed CentralView ArticleGoogle Scholar
- Morán-Diez E, Rubio B, Domínguez S, Hermosa R, Monte E, Nicolás C: Transcriptomic response of Arabidopsis thaliana after 24 h incubation with the biocontrol fungus Trichoderma harzianum. J Plant Physiol. 2012, 169: 614-620. 10.1016/j.jplph.2011.12.016.PubMedView ArticleGoogle Scholar
- Denoux C, Galletti R, Mammarella N, Gopalan S, Werck D, De Lorenzo G, Ferrari S, Ausubel FM, Dewdney J: Activation of defense response pathways by OGs and Flg22 elicitors in Arabidopsis seedlings. Mol Plant. 2008, 1: 423-445. 10.1093/mp/ssn019.PubMedPubMed CentralView ArticleGoogle Scholar
- Prins A, Mukubi JM, Pellny TK, Verrier PJ, Beyene G, Lopes MS, Emami K, Treumann A, Lelarge-Trouverie C, Noctor G, Kunert KJ, Kerchev P, Foyer CH: Acclimation to high CO2 in maize is related to water status and dependent on leaf rank. Plant Cell Environ. 2011, 34: 314-331. 10.1111/j.1365-3040.2010.02245.x.PubMedView ArticleGoogle Scholar
- Luhua S, Ciftci-Yilmaz S, Harper J, Cushman J, Mittler R: Enhanced tolerance to oxidative stress in transgenic Arabidopsis plants expressing proteins of unknown function. Plant Physiol. 2008, 148: 280-292. 10.1104/pp.108.124875.PubMedPubMed CentralView ArticleGoogle Scholar
- Van Hulten M, Pelser M, Van Loon LC, Pieterse CMJ, Ton J: Costs and benefits of priming for defense in Arabidopsis. Proc Natl Acad Sci USA. 2006, 103: 5602-5607. 10.1073/pnas.0510213103.PubMedPubMed CentralView ArticleGoogle Scholar
- Canet JV, Dobón A, Ibáñez F, Perales L, Tornero P: Resistance and biomass in Arabidopsis: a new model for salicylic acid perception. Plant Biotechnol J. 2010, 8: 126-141. 10.1111/j.1467-7652.2009.00468.x.PubMedView ArticleGoogle Scholar
- Pauwels L, Morreel K, De Witte E, Lammertyn F, Van Montagu M, Boerjan W, Inzé D, Goossens A: Mapping methyl jasmonate-mediated transcriptional reprogramming of metabolism and cell cycle progression in cultured Arabidopsis cells. Proc Natl Acad Sci U S A. 2008, 105: 1380-1385. 10.1073/pnas.0711203105.PubMedPubMed CentralView ArticleGoogle Scholar
- Marsoni M, Bracale M, Espen L, Prinsi B, Negri AS, Vannini C: Proteomic analysis of somatic embryogenesis in Vitis vinifera. Plant Cell Rep. 2008, 27: 347-356. 10.1007/s00299-007-0438-0.PubMedView ArticleGoogle Scholar
- Zhang J, Ma H, Chen S, Ji M, Perl A, Kovacs L, Chen S: Stress response proteins’ differential expression in embryogenic and non-embryogenic callus of Vitis vinifera L. cv. Cabernet Sauvignon–A proteomic approach. Plant Sci. 2009, 177: 103-113. 10.1016/j.plantsci.2009.04.003.View ArticleGoogle Scholar
- Sharifi G, Ebrahimzadeh H, Ghareyazie B, Gharechahi J, Vatankhah E: Identification of differentially accumulated proteins associated with embryogenic and non-embryogenic calli in saffron (Crocus sativus L.). Proteome Sci. 2012, 10: 3-10.1186/1477-5956-10-3.PubMedPubMed CentralView ArticleGoogle Scholar
- Correia S, Vinhas R, Manadas B, Lourenço AS, Veríssimo P, Canhoto JM: Comparative proteomic analysis of auxin-induced embryogenic and nonembryogenic tissues of the solanaceous tree Cyphomandra betacea (Tamarillo). J Proteome Res. 2012, 11: 1666-1675. 10.1021/pr200856w.PubMedView ArticleGoogle Scholar
- Rose RJ, Nolan KE: Vitro Cell Dev Biol. Invited review: Genetic regulation of somatic embryogenesis with particular reference to Arabidopsis thaliana and Medicago truncatula 42edition. 2006, 473-481.Google Scholar
- Almeida AM, Parreira JR, Santos R, Duque AS, Francisco R, Tomé DFA, Ricardo CP, Coelho AV, Fevereiro P: A proteomics study of the induction of somatic embryogenesis in Medicago truncatula using 2DE and MALDI-TOF/TOF. Physiol Plant. 2012, 146: 236-249. 10.1111/j.1399-3054.2012.01633.x.PubMedView ArticleGoogle Scholar
- Pavy N, Boyle B, Nelson C, Paule C, Giguère I, Caron S, Parsons LS, Dallaire N, Bedon F, Bérubé H, Cooke J, MacKay JJ: Identification of conserved core xylem gene sets: conifer cDNA microarray development, transcript profiling and computational analyses. New Phytol. 2008, 180: 766-786. 10.1111/j.1469-8137.2008.02615.x.PubMedView ArticleGoogle Scholar
- MacKenzie DJ, McLean MA, Mukerji S, Green M: Improved RNA extraction from woody plants for the detection of viral pathogens by reverse transcription-polymerase chain reaction. Plant Dis. 1997, 81: 222-226. 10.1094/PDIS.19184.108.40.206.View ArticleGoogle Scholar
- R Core Team: A language and environment for statistical computing. http://www.R-project.org,
- FlexArray: A statistical data analysis software for gene expression microarrays. McGill University and Genome Quebec Innovation Centre, http://www.gqinnovationcenter.com/services/bioinformatics/flexarray/index.aspx?l=e,
- Rutledge RG, Stewart D: Critical evaluation of methods used to determine amplification efficiency refutes the exponential character of real-time PCR. BMC Mol Biol. 2008, 9: 96-10.1186/1471-2199-9-96.PubMedPubMed CentralView ArticleGoogle Scholar
- LRE qPCR Website. http://sites.google.com/site/lreqpcr/,
- LRE qPCR Analyzer. http://code.google.com/p/lreqpcr/,
- Friedmann M, Ralph SG, Aeschliman D, Zhuang J, Ritland K, Ellis BE, Bohlmann J, Douglas CJ: Microarray gene expression profiling of developmental transitions in Sitka spruce (Picea sitchensis) apical shoots. J Exp Bot. 2007, 58: 593-614.PubMedView ArticleGoogle Scholar
- Czechowski T, Stitt M, Altmann T: 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
- Hruz T, Wyss M, Docquier M, Pfaffl MW, Masanetz S, Borghi L, Verbrugghe P, Kalaydjieva L, Bleuler S, Laule O, Descombes P, Gruissem W, Zimmermann P: RefGenes: identification of reliable and condition specific reference genes for RT-qPCR data normalization. BMC Genomics. 2011, 12: 156-10.1186/1471-2164-12-156.PubMedPubMed CentralView ArticleGoogle Scholar