TY - JOUR
AU - Cooke, Janice E.K.
AB - Abstract Conifers exhibit a number of constitutive and induced mechanisms to defend against attack by pests and pathogens such as mountain pine beetle (Dendroctonus ponderosae Hopkins) and their fungal associates. Ecological studies have demonstrated that stressed trees are more susceptible to attack by mountain pine beetle than their healthy counterparts. In this study, we tested the hypothesis that water deficit affects constitutive and induced responses of mature lodgepole pine × jack pine hybrids (Pinus contorta Dougl. ex Loud. var. latifolia Engelm. ex S. Wats. × Pinus banksiana Lamb.) to inoculation with the mountain pine beetle fungal associate Grosmannia clavigera (Robinson-Jeffrey and Davidson) Zipfel, de Beer and Wingfield. The degree of stress induced by the imposed water-deficit treatment was sufficient to reduce photosynthesis. Grosmannia clavigera-induced lesions exhibited significantly reduced dimensions in water-deficit trees relative to well-watered trees at 5 weeks after inoculation. Treatment-associated cellular-level changes in secondary phloem were also observed. Quantitative RT-PCR was used to analyze transcript abundance profiles of 18 genes belonging to four families classically associated with biotic and abiotic stress responses: aquaporins (AQPs), dehydration-responsive element binding (DREB), terpene synthases (TPSs) and chitinases (CHIs). Transcript abundance profiles of a TIP2 AQP and a TINY-like DREB decreased significantly in fungus-inoculated trees, but not in response to water deficit. One TPS, Pcb(+)-3-carene synthase, and the Class II CHIs PcbCHI2.1 and PcbCHI2.2 showed increased expression under water-deficit conditions in the absence of fungal inoculation, while another TPS, Pcb(E)-β-farnesene synthase-like, and two CHIs, PcbCHI1.1 and PcbCHI4.1, showed attenuated expression under water-deficit conditions in the presence of fungal inoculation. The effects were observed both locally and systemically. These results demonstrate that both constitutive and induced carbon- and nitrogen-based defenses are affected by water deficit, suggesting potential consequences for mountain pine beetle dynamics, particularly in novel environments. Introduction The current outbreak of mountain pine beetle (MPB; Dendroctonus ponderasae Hopkins) has resulted in the loss of more than 28 million hectares of pine forest in western North America since 1999, including more than 19 million hectares in the western Canadian provinces of British Columbia and Alberta (Government of Alberta 2012, Government of British Columbia 2012, Man 2012). In 2006, the MPB outbreak spread to the lodgepole pine (Pinus contorta Dougl. ex Loud. var. latifolia Engelm. ex S. Wats.) × jack pine (Pinus banksiana Lamb.) hybrid zone in north-central Alberta, and more recently has undergone a host-shift expansion into pure jack pine, a boreal forest species with a range extending east to the Atlantic coast (Cullingham et al. 2011). This unprecedented expansion of MPB into more northerly latitudes, into higher elevation forests and east into jack pine forests—all formerly considered marginal habitats for MPB—is postulated to have been facilitated at least in part by warmer than average temperatures during recent decades (Carroll et al. 2004, Trenberth et al. 2007). Mountain pine beetles overcome host defenses through a mass-attack strategy, i.e., concerted attack of a tree by an aggregation of insects, overcoming the critical threshold of resistance for physiologically stressed trees at lower attack densities than that of healthy and vigorously growing trees (Berryman 1982, Kolb et al. 1998, Wallin and Raffa 2002, Raffa et al. 2005, Boone et al. 2011). It is surmised that larger, healthier trees provide greater resources to MPB and their fungal pathogen associates, enabling greater fecundity and offspring survival (Raffa et al. 2008). Several fungal species belonging to the Ophiostomaceae are associated with MPB, facilitating an MPB attack and subsequent colonization of the host. Grosmannia clavigera (Robinson-Jeffrey and Davidson) Zipfel, de Beer and Wingfield (formerly Ophiostoma clavigerum) is considered to be the most pathogenic of the characterized MPB fungal associates, and even in the absence of the MPB is able to overcome tree defenses (Yamaoka et al. 1995, Lee et al. 2006, Safranyik and Carroll 2006). Grosmannia clavigera grows in the phloem and sapwood of the host tree (Parmeter et al. 1992), inducing water and mineral deficits that decrease the tree's overall vigor, including the ability to synthesize and mobilize secondary compounds. Fungal growth eventually disrupts water transport via tracheid occlusion, resulting in tree mortality (Christiansen et al. 1987, Yamaoka et al. 1995, Huber et al. 2004). Trees defend against attack by pests and pathogens via a suite of constitutive and inducible defense mechanisms that include physical, cellular, biochemical and molecular processes (Keeling and Bohlmann 2006, Van Loon et al. 2006, Eyles et al. 2010). In response to bark invaders such as MPB and G. clavigera, pine species produce a reaction zone called a lesion in both phloem and xylem in an attempt to contain the invading organism (Hudgins et al. 2003, Franceschi et al. 2005). This lesion contains necrotic cells near the point of invasion, as well as cells containing quantities of defensive compounds (Franceschi et al. 2005). Conifer chemical defenses include carbon-based metabolites such as terpenoids and phenolics that create chemical and physical barriers to attackers (Byun-McKay et al. 2003, 2006, Fäldt et al. 2003, Miller et al. 2005, Keeling and Bohlmann 2006, Van Loon et al. 2006). Terpenoids, a major component of the oleoresins that are a cornerstone of the conifer defense response, are synthesized by a complex network of reactions that include steps catalyzed by terpene synthases (TPSs). Terpene synthases can be further classified as mono-, sesqui- and diterpene synthases (Byun-McKay et al. 2006, Zulak and Bohlmann 2010). Mono- and diterpene synthases producing many of the resin terpenoids and volatiles of lodgepole pine and jack pine have recently been described (Hall et al. 2013a, 2013b). Protein, i.e., nitrogen-based, defenses also appear to figure prominently in the conifer defense response (Lippert et al. 2007). Chitinases (CHIs) are a well-characterized class of defense proteins in both angiosperms (Bishop et al. 2000) and conifers (Davis et al. 2002, Fossdal et al. 2005). Chitinases that exhibit chitinolytic activity hydrolyze chitin, a cell wall constituent of many fungi, and may also be involved in generating signal molecules that function as endogenous elicitors of other defense mechanisms (Bishop et al. 2000, Van Loon et al. 2006, Miranda et al. 2007). Members of the CHI gene family also function in non-chitinolytic roles in plants, such as antifreeze proteins and vegetative storage proteins (Islam et al. 2011). Water deficit is an abiotic stress that affects tree defenses (Dunn and Lorio 1993, Franceschi et al. 2005, McDowell et al. 2008). Lombardero et al. (2000) have demonstrated that constitutive oleoresin flow was increased in water-stressed loblolly pine (Pinus taeda L.) relative to unstressed trees, but wound-induced resin flow was decreased. Lusebrink et al. (2011) demonstrated that both constitutive and G. clavigera-induced monoterpene volatile profiles of both lodgepole pine and jack pine seedlings are affected by water deficit. Northern Alberta has experienced periods of drought in the past two decades (Hanesiak et al. 2011) that have had measureable negative impacts on these forests (Chhin et al. 2008, Michaelian et al. 2011). If the drought experienced in the region has adversely impacted constitutive and induced defense responses of pines, the critical threshold for a successful MPB attack may be lowered. Given that MPB population densities are comparatively low at the leading edges of an outbreak, factors that serve to reduce the attack threshold at which MPB successfully overcome host tree defenses could increase the probability of MPB establishment in novel habitats, facilitating continued range expansion (Safranyik and Carroll 2006, Raffa et al. 2008). In this study, we have examined the effect of water limitation on molecular defense responses of mature lodgepole × jack pine hybrid trees in central Alberta to the MPB fungal associate G. clavigera. These hybrids are considered ‘naïve’ hosts, as it is thought that the region has not experienced an MPB outbreak in the past (Safranyik et al. 2010). We hypothesized that water deficit would lead to reduced photosynthesis, concomitantly reducing carbon gain and the availability of photosynthate for allocation to defenses. The degree of water stress imposed by the water limitation treatment was determined by measuring gas exchange parameters and conducting quantitative RT-PCR (qRT-PCR) transcript profiling of genes well known for their roles in mediating drought responses: aquaporin (AQP) water channels and dehydration-responsive element binding (DREB) transcription factors. The defense response of these trees was assessed by measuring lesion length, changes to phloem cellular structure and qRT-PCR transcript profiling of secondary phloem for selected genes encoding TPSs and CHIs. We tested whether all defense-associated genes would respond in the same manner to water-deficit-induced stress, and if not, whether reduced carbon gain would impact the trees' carbon-based defenses more than nitrogen-based defenses. We also examined whether the effect of water deficit on tree defense gene expression was manifested both locally and systemically. The results of the study are discussed within the context of pine–MPB co-evolutionary relationships and ecological plant defense theory. Materials and methods Field study The study was conducted at Chickadee, a naturally regenerated, thinned stand of ∼60-year-old P. contorta × banksiana hybrid pine trees (Pcb) located 37 km northwest of Whitecourt, Alberta, Canada (54°13′N, 116°03′W). The plot size area was 5513 m2. The hybrid status of the trees in the stand was confirmed by microsatellite analysis according to Cullingham et al. (2011). The experiment was conducted from May to August 2009, corresponding to Julian Days 147–226. Meteorological data to describe environmental conditions at the study site were obtained from Environment Canada's Whitecourt weather station (http://www.climate.weatheroffice.gc.ca). Regional meteorological data from Environment Canada were used as the input to derive Climate Moisture Index (CMI) values according to Hogg (1997). The model was implemented within the BioSIM framework (Régnière 1996); 1971–2000 values were computed using 20 random replicates, and used to generate spatially interpolated maps. Water deficit and fungal inoculation treatments A total of 40 trees with no evident history of disease or mechanical damage, of similar height and with a diameter at breast height between 16 and 24 cm were randomly selected. A two-by-two factorial design with two treatment factors, water level and inoculation with G. clavigera, was performed. The 40 experimental trees were randomly assigned to one of four treatments, with 10 trees per treatment: well-watered/uninoculated control, water-deficit/uninoculated control, well-watered/G. clavigera inoculated and water-deficit/G. clavigera inoculated. Well-watered and water-deficit trees were separated by a minimum of 11 m and a maximum of 126 m in the experimental design. Water treatments (well-watered or water deficit) were applied for 6 weeks before the fungal inoculation was applied. Water deficit was induced placing woven laminated polyethylene tarps (12′ × 14′; GH Factory Sales, Surrey, BC, Canada) around the base of each tree extending to the drip line from 27 May until termination of the experiment on 13 August. Well-watered trees were watered every 2 weeks from 12 June to 7 August using 40-gallon arborist's water bags (Treegator®, Rittenhouse, St. Catharines, ON, Canada) placed around the base of each tree. Soil relative water content (RWC) was measured using a time-domain reflectometer (TDR; Tektronix 1502B Cable TDR Cable Tester; Tektronix, Inc., Scarborough, ON, Canada). The relative water content (θ) was determined as a percentage calculated using an equation for organic soils calibrated for sandy soil (Topp et al. 1980, Robinson et al. 2003). Probe rods for measuring the RWC were placed 1-m away from the trunk, at three different depths (30, 60 and 90 cm) in the ground. Trees were inoculated with G. clavigera—considered a moderately virulent fungal pathogen (Six and Wingfield 2011)—on 9 July 2009 according to Rice et al. (2007). The G. clavigera inoculum was derived from cultures initiated from an isolate cultured from a single MPB larva collected in 2007 from a naturally attacked tree near Fox Creek, Alberta, Canada (54°48′N, 116°63′W), close to the study site. It was reasoned that the haplotype represented by this culture would be the same or similar to G. clavigera carried by MPB attacking trees in the region. Molecular analyses confirmed the identity of the isolate as G. clavigera (Roe et al. 2010, 2011). The isolate corresponds to specimen M001-03-03-07-UC04DL09 described in Roe et al. (2010, 2011). A cordless drill was used to make 40 holes of ∼8 mm diameter per tree; holes were positioned 2–3 cm apart and arranged in two rings, one at breast height and the second 40 cm above breast height. Following inoculation, each hole was plugged with a small piece of sterilized dowel rod, and the ring of inoculations was wrapped using Parafilm® strips (Pechiney PM992) to reduce environmental contamination. Mock-inoculated trees were not included in the experimental design due to the logistics of adding additional trees to a study conducted in a natural stand. Approximately 5 weeks after inoculation, phloem was sampled by removing the bark above and below the inoculated area, and peeling the phloem layer away from the remaining outer tissue of the bark. Phloem adjacent to the maximum vertical extension of the lesion (outside lesion), and ∼1 cm inside of the lesion boundary (inside lesion), were collected separately for gene expression. Tissues collected for gene expression analysis were frozen in liquid nitrogen immediately after harvesting and stored at –80 °C until analysis. Eight to 10 trees were sampled for each treatment. Following harvest of material for molecular analyses, lesions were measured from the remaining bolts (wood) on 12 August 2009. Vertical extents of the lesion were delimited as the outermost discolored edge on the sapwood above and below the inoculation point. Gas exchange Stomatal conductance (gs) and photosynthesis (A) were measured between 10:00 and 15:00 h using an open gas exchange system (Li-6400; Li-Cor, Inc, Lincoln, NE, USA). Measurements were performed on three fascicles of fully expanded current year needles. The plants were exposed to red/blue light (6400-02B LED Light Source) of ∼1200 μmol photons m−2 s−1 for 5 min prior to measurements. The cuvette CO2 concentration was set to 400 μmol CO2 mol−1. Leaf temperature in the chamber was controlled between 18 and 20 °C using a Peltier cooling system. Gas exchange was measured three times on each set of needles and the mean of the three measurements used for the analysis. Microscopy Preparation of material for staining was carried out as described in El Kayal et al. (2011). Pieces of phloem ∼3 mm2 were placed directly in fixative (2% v/v glutaraldehyde, 1% caffeine buffered in 0.1 M sodium phosphate buffer; pH 7.2). Following successive ethanol dehydration and infiltration, samples were embedded in JB-4 Plus® (Polysciences, Hatfield, PA, USA) at 55 °C overnight. Four-micrometer sections were stained with Richardson's stain (Richardson et al. 1960) or periodic acid-Schiff's (PAS) reagent for sugar residues with vicinal hydroxyl groups (Franceschi et al. 2000). In the latter procedure, sections were incubated for 30 min in 1% w/v periodic acid at room temperature, washed with distilled water and air-dried. Some sections were incubated in distilled water as non-oxidized control samples. Sections were then incubated in Schiff's reagent in the dark at room temperature for 1 h, rinsed with distilled water, counterstained with hematoxylin, and then dehydrated and mounted. Sequence analysis Expressed sequence tags are publicly accessible at National Center for Biotechnology Information (NCBI) under accession numbers GW738148-GW774480 (P. banksiana) and GT229582-GT270713 (P. contorta). The corresponding transcriptome resources are described in Hall et al. (2013b). Contigs corresponding to putative AQP, DREB, CHI and TPS genes were identified by BLASTX using previously characterized sequences from other species as queries. Sequences used for subsequent analyses are provided in File S1 available as Supplementary Data at Tree Physiology Online. Open reading frames were predicted using Geneious 4.8.2 (Drummond et al. 2010). Multiple sequence alignments of deduced amino acid sequences and neighbor-joining phylogenetic analyses were performed in MEGA v4.0 (Tamura et al. 2007). Dendrograms were constructed using pairwise deletion and 1000 bootstrap replicates. The phylogenetic analyses were used to select genes for qRT-PCR based on position relative to functionally characterized genes from other species. Quantitative RT-PCR Quantitative RT-PCR was carried out as described in El Kayal et al. (2011). Total RNA extractions (∼100 mg tissue per extraction) were carried out according to Pavy et al. (2008), quantified with a NanoDrop 1000 (Thermo Fisher Scientific, Waltham, MA, USA) and assessed with a 2100 Bioanalyzer (Agilent, Mississauga, ON, Canada). Two micrograms of total RNA was treated with DNaseI (Invitrogen; Life Technologies, Burlington, ON, Canada) prior to cDNA synthesis using Superscript II reverse transcriptase, following the manufacturer's protocols (Invitrogen). Primer Express (v3.0, Applied Biosystems, Mississauga, ON, Canada) was used to design primers for qRT-PCR (see Table S1 available as Supplementary Data at Tree Physiology Online). Quantitative RT-PCR reactions (10 μl) consisting of a master mix (0.2 mM dNTPs, 0.3 U Platinum Taq Polymerase (Invitrogen), 0.25× SYBR Green, and 0.1× ROX), 20 ng diluted cDNA and 0.8 μM primers were analyzed as described in El Kayal et al. (2011). Three technical replicates were run for each biological replicate. Transcript abundance was quantified using standard curves generated from serial dilutions of P. contorta PCR products. Eukaryotic translation initiation factor 5A-1 (eIF-5A; accession number KF322083) was used as the reference gene, since transcript abundance corresponding to this gene was not significantly different across treatments (P = 0.34). Statistical data analysis Analysis of variance (ANOVA) was performed using the Proc Mixed procedure of SAS (9.1) (SAS Institute, Cary, NC, USA, 2002). One-way ANOVA was performed to test for variation in lesion length data. Two-way ANOVA was used to test for variation in gas exchange and gene expression under water level and pathogen treatment combinations. A Kenward–Roger option was used to calculate denominator degrees of freedom. Least squared means were separated using the pdiff option, and P ≤ 0.05 was used to determine significances between treatments. To meet the assumptions of normality and homogeneity of variance, log transformation was applied to the gene expression values for all but two genes. Principal component analysis (PCA) was used to assess the relationship between genes with similar expression profiles. The PCA was performed using the Vegan package from R version 2.11.0 (R Development Core Team 2010; http://www.R-project.org). Significant contribution to the principal components was determined using the broken-stick distribution. Results Soil water status in the field experiment The period during which the field experiment was conducted—May to August 2009—was drier than normal in much of Alberta, particularly in the central region of the province where the study was conducted. BioSIM (Régnière 1996) was used to generate spatially interpolated maps of Hogg's CMI, which takes into account both precipitation and potential evapotranspiration (Hogg 1997). Figure 1 illustrates the reduced moisture availability in the region of the study site in 2009 relative to the 1971–2000 mean, particularly for the summer months (May to August) when the experiment took place. The CMI at the Chickadee study site was negative throughout the duration of the study—indicating a moisture deficit—and was lower than the CMI for the 1971–2000 mean (Figure 2a). Precipitation in the study area was lower than the 1971–2000 mean throughout much of the experiment except for July, when near-normal precipitation levels were due mainly to a single weather event (Figure 2). In contrast, monthly temperature averages for 2009 were very similar to those for the 1971–2000 mean, particularly during the experiment (Figure 2). Daily maximum temperatures were at or below normal for much of the duration of the experiment, notably from mid-June to mid-July (Julian Days 170–196) and during data and sample collection in August (Julian Days 222–226). Figure 1. Open in new tabDownload slide A comparison of normal (1971–2000) average CMI values for western Canada with those for 2009, the year during which the study was conducted. The inset maps depict the area of central Alberta in which the study site was located. The CMI was computed according to Hogg (1997), implemented within BioSIM (Régnière 1996) to generate spatially interpolated maps. (a) Annual average CMI for 1971–2000. (b) Annual CMI for 2009. (c) Average CMI for the months of May–August 1971–2000. (d) Average CMI for the months of May–August 2009, the period during which the study was conducted. Negative CMI values signify that mean precipitation is less than the estimated mean potential evapotranspiration, indicative of a predicted moisture deficit (Hogg 1997). Figure 2. Open in new tabDownload slide Comparison of 2009 to 1971–2000 average monthly and daily environmental data for the study site district. Data were obtained from Environment Canada's Whitecourt weather station, ∼37 km from the study site. (a) Climate Moisture Index calculated on a monthly basis for 2009 and the 1971–2000 average. Negative CMI values indicate a predicted moisture deficit. (b) Precipitation expressed as a monthly average, (c) maximum daily temperature expressed as a monthly average, (c) climate moisture index using the average of 1971–2000 and recorded data from 2009 in Whitecourt (AB). (d) Minimum and maximum temperature and precipitation averaged for 1971–2000 and the year of 2009. In a–c, black bars illustrate 2009 data, and white bars illustrate 1971–2000 data. In d, gray diamonds indicate daily temperature minima for 2009, black squares indicate daily temperature maxima for 2009, the gray thick line indicates average daily temperature minima for 1971–2000, the black thick line indicates average daily temperature maxima for 1971–2000, and white triangles indicate 2009 daily precipitation. Soil relative water content values at a depth of 90 cm averaged 17 ± 1.3% and 11 ± 0.6% for well-watered and water-deficit trees (Figure 3), respectively, demonstrating that water availability differed between well-watered and water-deficit treatments (P < 0.0001). These values bracketed those obtained for ambient soil RWC obtained at the study site. Relative water content values were relatively constant in well-watered trees at each of the three measured depths. Rain events did not have substantive effects on the RWC. Figure 3. Open in new tabDownload slide Soil RWC corresponding to well-watered and water-deficit trees over the course of the study. Mean values (±SE) for 20 trees across the field site, measured at several time points, are shown. Measurements were taken at 90 cm below the soil surface. RWC data obtained at 30 and 60 cm below the soil surface are not shown. P < 0.0001 indicates significance at the 5% level between RWC of well-watered and water-deficit treatments. Precipitation over the duration of the experiment, as measured at the Environment Canada Whitecourt station, is also shown. Physiological and anatomical responses to plant water and fungal inoculation status Stomatal conductance and A are key physiological indicators of water limitation stress for isohydric species such as pines (Klein et al. 2011). Mature P. contorta × banksiana hybrids trees subjected to water deficit exhibited significantly reduced A (P = 0.02) and nearly significantly reduced gs (P = 0.09) after 11 weeks of treatment (Figure 4). Inoculation with G. clavigera did not significantly affect either A (P = 0.38) or gs (P = 0.56). Similarly, the interaction terms for both A and gs were not significant (P = 0.89 and P = 0.86, respectively). Figure 4. Open in new tabDownload slide Water-deficit-induced changes to A and gs. (a) Photosynthesis (μmol CO2 m−2 s−1) and (b) gs (mmol m−2 s−1) in mature P. contorta × banksiana trees subjected to well-watered or water-deficit conditions, either as unioculated controls or 5 weeks following G. clavigera inoculation. For each box plot, the upper and lower limits of the box indicate the 25th and 75th percentiles; the line through the box is the median; the capped lines indicate the 10th and 90th percentiles, and the black circle is the mean (n = 8–10). Analysis of variance results are indicated in the top right corner of each panel. At 5 weeks post-inoculation, lesions surrounding each inoculation point were generally of a uniform elongated oval shape extending equally above and below the inoculation point. Lesion lengths were significantly greater in well-watered trees than in water-deficit trees (P = 0.0004; Figure 5). Figure 5. Open in new tabDownload slide Lesion lengths of mature P. contorta × banksiana trees 5 weeks post-inoculation with G. clavigera. Different letters denote statistical differences at P = 0.05. Values are the mean ± SE for 10 trees. Ten lesions were measured per tree. Comparison of phloem anatomy between the four treatments revealed slight swelling of the polyphenolic parenchyma cells (PP cells) in inoculated trees, concomitant with compressed sieve cells (Figure 6). Polyphenolic parenchyma cells exhibited more densely staining phenolic bodies under water deficit compared with well-watered conditions. In addition, more abundant bodies, assumed to be starch grains due to preferential staining with PAS, were apparent around the PP cell periphery in inoculated well-watered trees. The PP cell phenolic bodies were mostly lightly stained with PAS, with a few of the PP cells in the outer layers having darker stained phenolic bodies (see Figure S1 available as Supplementary Data at Tree Physiology Online). Figure 6. Open in new tabDownload slide Effect of water deficit and G. clavigera on the secondary phloem cellular structure of P. contorta × banksiana mature trees. Transverse sections were stained with Richardson's stain. (a) Well-watered, uninoculated control tree. (b) Well-watered tree, 5 weeks following inoculation with G. clavigera. (c) Water deficit, uninoculated tree. (d) Water-deficit tree, 5 weeks following inoculation with G. clavigera. Polyphenolic parenchyma cells are more densely stained in well-watered inoculated trees relative to the other treatments. Arrows indicate collapsed sieve cell layers that are more evident in inoculated trees (b and d), but were also observed in some water-deficit uninoculated trees. PP, polyphenolic parenchyma cells; R, ray parenchyma; S, sieve cells. Scale bar = 100 μm. Water-deficit- and fungal-induced changes in transcript abundance profiles of putative drought-associated genes To determine whether the water stress imposed upon the trees was sensed at the molecular level and what effect G. clavigera inoculation had on this response, we generated transcript abundance profiles corresponding to several genes belonging to the AQP and DREB families. Members of the AQP and DREB families are classically associated with responses to water deficit and other abiotic stresses, and as such represent appropriate reporters for detecting plant responses to water deficit (Ozturk et al. 2002, Agarwal et al. 2006, Almeida-Rodriguez et al. 2010). Since no sequence resources were available for P. contorta × banksiana, we identified AQP- and DREB-like sequences for transcript abundance profiling by mining Sanger sequence data for pure P. contorta and P. banksiana. The P. contorta and P. banksiana sequences were characterized by phylogenetic analyses together with the corresponding Arabidopsis thaliana gene family and selected functionally characterized genes from other species, as well as by identification of conserved motifs typical of AQPs and DREBs (see Figures S2–S5 available as Supplementary Data at Tree Physiology Online). Because P. contorta × banksiana hybrids in this region typically demonstrate a greater proportion of P. contorta ancestry (Cullingham et al. 2012), we elected to choose P. contorta genes as the basis for designing qRT-PCR assays. Four pine AQP genes of the tonoplast intrinsic protein (TIP) and plasma membrane intrinsic protein (PIP) subfamilies were selected for transcript abundance profiling (see Figures S2 and S3 available as Supplementary Data at Tree Physiology Online). Members of these subfamilies have been shown to function as water channels and some are involved in water-deficit responses (Alexandersson et al. 2005, Zhang et al. 2007, Kaldenhoff et al. 2008). The P. contorta deduced amino acid sequences included typical MIP features such as a highly conserved asparagine–proline–alanine (NPA) motif and histidine residues involved in pH sensing (see Figure S2 available as Supplementary Data at Tree Physiology Online; Shelden et al. 2009, Park et al. 2010). The selected P. contorta genes were named PcPIP1;1, PcPIP2;1, PcTIP1;1 and PcTIP4;1 to reflect their placement in the phylogenetic tree and AQP naming conventions (see Figure S2 available as Supplementary Data at Tree Physiology Online). Transcript abundance profiles for the corresponding genes in P. contorta × banksiana hybrids indicated that water deficit had no significant effect on transcript abundance in control samples for any of the four AQPs that were examined. Similarly, there were no significant differences in transcript levels of PIP2;1 TIP1;1 and TIP4;1 (designated PcbPIP2;1, PcbTIP1;1 and PcbTIP4;1 for the hybrids) across water level and fungal inoculation treatments (Figure 7b–d). However, there was a significant decrease in transcript abundance of PcbPIP1;1 in fungus-inoculated samples (P < 0.0001) both inside and outside of the lesion, although the interaction between water and inoculation treatments was nearly significant (P = 0.07; Figure 7a). Figure 7. Open in new tabDownload slide AQP transcript profiling in secondary phloem of mature P. contorta × banksiana trees subjected to water deficit and/or inoculation with G. clavigera for 35 days. Phloem from trees subjected to G. clavigera inoculation were sampled both outside and within the lesion zone for expression analyses. eIF-5A was used as the reference gene. Values represent the mean ± SE, n = 7–10. Analysis of variance results are indicated in the top right corner of each panel. For statistical analyses, data for most genes were log transformed to meet the assumptions of normality and homogeneity of variance. No pine sequences were identified in the Sanger dataset with similarity to the A. thaliana DREB1 or DREB2, the best characterized of the dehydration-responsive DREB family members. Accordingly, five other members of the family were selected for transcript profiling, and named PcTINY-like1, PcTINY-like2, PcDREB-like, PcRAP2.4-like and PcERF61-like to reflect their relative placement on the tree. In silico analyses of these five P. contorta sequences indicated conservation of key amino acids necessary for DREB function (see Figure S5 available as Supplementary Data at Tree Physiology Online), including conserved tryptophan and arginine residues that interact with the GCC box cis-acting promoter element (Allen et al. 1998), an A38 residue important for stabilizing binding of the AP2 domain (Liu et al. 2006), either a serine or alanine at position 15 important for binding of DREB family transcription factors to ethylene-responsive elements (ERE) (Sun et al. 2008), and conservation of valine (V14) and glutamic acid (E19) residues, important for DREB binding to DRE elements (Sakuma et al. 2002). As with the AQPs, none of the DREB family members that were examined showed significant differences in transcript abundance between control well-watered and control water- deficit treatments in mature P. contorta × banksiana (Figure 8). However, both water deficit and fungal inoculation resulted in significantly decreased transcript abundance of PcbTINY-like1 (P = 0.03 and P < 0.0001, respectively), with the interaction term also being significant (P = 0.02; Figure 8a). Fungal inoculation was associated with a significant decrease in transcript levels of PcbTINY-like2 and PcbRAP2-like (P = 0.005 and P = 0.002, respectively), although water treatment did not have a significant effect on transcript abundance for these two genes, and the interaction term was not significant (Figure 8b and c). Conversely, water deficit resulted in significantly increased transcript abundance for PcbERF61-like and PcbDREB-like, compared with well-watered trees (P = 0.02 and P = 0.0008, respectively; Figure 8d and e). Decreases in transcript abundance on PcbERF61-like were observed due to fungal inoculation, although the effect was marginally significant (P = 0.09). A similar trend was observed in PcbDREB-like under well-watered conditions (Figure 8d and e). Figure 8. Open in new tabDownload slide DREB transcript profiling in secondary phloem of mature P. contorta × banksiana trees subjected to water deficit and/or inoculation with G. clavigera for 35 days. Phloem from trees subjected to G. clavigera inoculation were sampled both outside and within the lesion zone for expression analyses. eIF-5A was used as the reference gene. Values represent the mean ± SE, n = 7–10. Analysis of variance results are indicated in the top right corner of each panel. For statistical analyses, data were log transformed to meet the assumptions of normality and homogeneity of variance. Different letters denote statistical difference at P < 0.05. Water-deficit- and fungal-induced changes in transcript abundance profiles of defense-associated genes Chitinases and TPSs are frequently associated with both constitutive and induced defense responses in conifers (e.g., Davis et al. 2002, Fossdal et al. 2005, Ralph et al. 2006, Zulak et al. 2009, Islam et al. 2010). Chitinases have been traditionally grouped into Classes I–VII based on biochemical and structural characterizations (Hamel et al. 1997, Brunner et al. 1998, Kasprzewska 2003), while more recent phylogenetic analyses group members into six clusters designated Clusters 1–6 (Xu et al. 2007). Phylogenetic analyses showed that the four P. contorta CHIs selected for transcript profiling were positioned in Clusters 1 and 3, which contain well-characterized pathogen-response CHIs belonging to Class IV and Classes I/II, respectively (see Figure S6 available as Supplementary Data at Tree Physiology Online). The P. contorta CHIs were named PcCHI1.1, PcCHI2.1, PcCHI2.2 and PcCHI4.1, based on their placement in the phylogenetic tree. Placement within the tree enabled identification of PcCHI4.1 as a Class IV CHI, while placement of the other three CHIs suggested that they could belong to either Class I or II. Multiple alignment of the deduced amino acid sequences was performed to determine the presence of different CHI domains used to distinguish different classes (see Figure S7 available as Supplementary Data at Tree Physiology Online). The presence of a chitin-binding domain, a complete catalytic domain and a carboxy-terminal extension showed that PcCHI1.1 was a Class I CHI. PcCHI2.1 and PcCHI2.2 exhibited a deletion in the catalytic domain but no chitin-binding domain or carboxy-terminal extension, typical of Class II CHIs. PcCHI4.1 exhibited a chitin-binding domain, no carboxy-terminal extension, and two deletions on the catalytic binding domain, providing further evidence that PcCHI4.1 encodes a Class IV CHI (see Table S2 available as Supplementary Data at Tree Physiology Online). Quantitative RT-PCR transcript analysis revealed that most of the CHIs profiled responded to at least one of the two treatment factors (Figure 9). PcbCHI1.1 transcript abundance significantly increased under both water availability and fungal inoculation treatments (P = 0.0006 and P < 0.0001); the interaction term was also significant (P = 0.005; Figure 9a). PcbCHI4.1 showed significant increases in transcript levels due to fungal inoculation, but expression was unaffected by water deficit (P < 0.0001; Figure 9b). PcbCHI2.1 and PcbCHI2.2 both exhibited significantly increased transcript abundance under water deficit (P = 0.001 and P = 0.002, respectively; Figure 9c and d), but while PcbCHI2.1 transcript abundance was marginally significantly greater in fungus-inoculated versus control samples (P = 0.08), fungal inoculation had no significant effect on PcbCHI2.2 transcript abundance (Figure 9c and d). PcbCHI5.1-like did not show any significant changes in transcript abundance (Figure 9e). Figure 9. Open in new tabDownload slide Chitinase transcript profiling in secondary phloem of mature P. contorta × banksiana trees subjected to water deficit and/or inoculation with G. clavigera for 35 days. Phloem from trees subjected to G. clavigera inoculation were sampled both outside and within the lesion zone for expression analyses. eIF-5A was used as the reference gene. Values represent the mean ± SE, n = 7–10. Analysis of variance results are indicated in the top right corner of each panel. For statistical analyses, data were log transformed to meet the assumptions of normality and homogeneity of variance. Different letters denote statistical difference at P < 0.05. Several lodgepole pine TPS genes have recently been cloned and the corresponding enzyme functions characterized (Hall et al. 2013a, 2013b). Four TPS genes were profiled: PcTPS-(+)α pin1, PcTPS-3car1, PcTPS-LAS1 and Pc(E)-β-farnesene synthase-like (see Figures S8 and S9 available as Supplementary Data at Tree Physiology Online; Hall et al. 2013a, 2013b). Except for Pc(E)-β-farnesene synthase-like, substrate specificity consistent with their naming has been demonstrated for each of these TPSs (Hall et al. 2013a, 2013b). The four TPSs displayed diverse transcript profiles in response to water deficit and fungal inoculation (Figure 10). Pcb(E)-β-farnesene synthase-like transcript levels showed significant differences under water availability (P = 0.01) and fungal inoculation (P = 0.02; Figure 10a). Although the interaction was not significant, increases in transcript abundance were observed under well-watered but not water-deficit conditions (Figure 10a). PcbTPS-3car1 also showed significant differences in transcript abundance under water availability (P = 0.01) and fungal inoculation treatments (P = 0.0018), but the interaction term was not significant (Figure 10b). PcbTPS-LAS1and PcbTPS-(+)α pin1 showed no significant differences between any of the treatments (Figure 10c and d). Figure 10. Open in new tabDownload slide Terpene synthases transcript profiling in secondary phloem of mature P. contorta × banksiana trees subjected to water deficit and/or inoculation with G. clavigera for 35 days. Phloem from trees subjected to G. clavigera inoculation were sampled both outside and within the lesion zone for expression analyses. eIF-5A was used as the reference gene. Values represent the mean ± SE, n = 7–10. Analysis of variance results are indicated in the top right corner of each panel. For statistical analyses, data were log transformed to meet the assumptions of normality and homogeneity of variance. Multivariate analyses of transcript profiles Principal component analysis was conducted on the transcript abundance profiles to identify the gene(s) explaining the greatest proportion of the variance between treatments (Figure 11). The first principal component, which was significant, accounted for 67.55% of the variability. Uninoculated treatments, water-deficit inoculated treatments and well-watered inoculated treatments separated into three distinct groups along PC1. The biplot showed that the inoculated/well-watered treatments were associated with increased transcript abundance of Pcbchi1.1, Pcbchi4.1 and Pcb(E)-β-farnesene synthase-like, while the control uninoculated treatments were associated with increased transcript abundance for PcbPIP1;1, although none was statistically significant according to the broken-stick model. Figure 11. Open in new tabDownload slide Principal component analysis to identify relationships between AQP, DREB, TPS and CHI phloem transcript abundance profiles, water availability and G. clavigera inoculation. The percentage of variation in the data explained by each principal component is provided on the graph axes, and principal components that are found to be significant based on the broken-stick distribution are noted with an asterisk. Gene expression profiles are shown as eigenvectors (arrows) and those that exceed the equilibrium circle contribute significantly to the principal components. WW, well-watered uninoculated trees; WD, water-deficit uninoculated trees; WW-I, well-watered fungus-inoculated trees; WD-I, water-deficit fungus-inoculated trees. Sampling of phloem inside or outside of the lesion is additionally indicated. A2, PcbPIP1;1; C1, Pcbchi2.2; C2, Pcbchi4.1; C3, Pcbchi1.1; C5, Pcbchi2.1; T2, Pcb(E)- ± -farnesene synthase-like; T3, Pcb(+)-3-carene synthase. The second principal component accounted for 27.11% of the variability. Water-deficit treatments were separated from well-watered treatments along PC2, which was driven mainly by increased expression of PcbTPS-3car1, Pcbchi2.1-like and Pcbchi2.2-like in water-deficit treatments and increased expression of Pcb(E)-β-farnesene synthase-like in well-watered treatments, although none was considered significant. Discussion Water deficit significantly reduced A in mature P. contorta × banksiana Bark beetles such as MPB are able to overcome the critical threshold of resistance for physiologically stressed trees at lower attack densities than that of healthy and vigorously growing trees (Berryman 1982, Kolb et al. 1998, Wallin and Raffa 2002, Raffa et al. 2005). Our analyses establish that water- deficit conditions were experienced at the study site in 2009, the year in which the present study took place. The degree of water deficit to which the water-limited trees in the study site were exposed was augmented via water exclusion from the root area extending to the tree's drip line, while well-watered treatments were not subjected to this water deficit through active watering. Pine species exhibit a relatively conservative water-use strategy, and exert strong regulation over gs for control of transpiration under water-deficit conditions (Andrews et al. 2012). As in previous studies with pines (e.g., Guérard et al. 2007, Domec et al. 2009, Lusebrink et al. 2011), the degree of water limitation imposed by the combination of environmental conditions and experimental manipulation was sufficient to cause decreased gs and A in the water-deficit P. contorta × banksiana hybrids, although the former was not statistically significant at α = 0.05. Consequently, we can infer that carbon gain was reduced in water-deficit trees relative to well-watered trees. There are at least three possible explanations for why greater differences in gs were not observed. First, there is considerable tree to tree variation in both photosynthesis and gs values, as illustrated by the box plot in Figure 4. This likely reflects genetic variation between biological replicates—the study was conducted in a naturally regenerated stand rather than a plantation—as well as micro-site variation. Second, temperatures were at or below normal during the week in which A and gs were measured. We speculate that under these average temperatures, evapotranspiration may have been insufficient to promote greater stomatal closure. If the water-limiting conditions had been accompanied by high temperatures, as is often the case under field conditions, it is possible that gs would have been further reduced in water-deficit trees. Third, pure jack pine seedlings from Saskatchewan show greater gs under water-deficit conditions than lodgepole pine seedlings from western Alberta (Lusebrink et al. 2011). Jack pine is also more likely to be found on drier sites than lodgepole pine (Cullingham et al. 2012). Accordingly, it is possible that lodgepole × jack pine hybrids are more tolerant to water deficit than to lodgepole pine, although this remains to be tested. Water deficit alters P. contorta × banksiana defense responses to G. clavigera infection Lesion length is a commonly used measure in studies of bark beetle–fungal associate–pine host interactions. In some studies, lesion length is considered an indicator of fungal pathogenicity and a measure of fungal success (e.g., Wong and Berryman 1977, Plattner et al. 2008, Rice and Langor 2008), with longer lesions generally interpreted to indicate greater fungal virulence (e.g., Lu et al. 2010). Lesion length is also considered an indicator of tree defense capacity (Paine and Stephen 1987), with longer lesions reported to be a reflection of increased release of toxic and/or inhibitory substances by trees in a defense response (e.g., Raffa and Berryman 1983, Christiansen et al. 1987, Raffa and Smalley 1995, Witzell and Martín 2010). In this study, inoculation with G. clavigera resulted in significantly larger lesions in well-watered trees compared with water-deficit trees after 5 weeks. Given that lesion length can be interpreted as an enhanced defense response and/or as a proxy for the extent of fungal invasion of the plant tissue, it is possible that reduced lesion length in water-deficit trees resulted from (i) induction of a defense response reduced either in total magnitude or delayed temporally and/or (ii) reduced fungal growth in the tree due to a less hospitable environment. Although not measured in this experiment, PCR detection of G. clavigera in inoculated lodgepole and jack pine seedlings suggest that the extent of G. clavigera growth is similar to the extent of lesion formation over time (A. Arango-Velez and J.E.K. Cooke, unpublished data), implying that these two possible explanations are not mutually exclusive. Further investigations of lesion development dynamics and the growth of fungi within the developing lesion are warranted to better understand how this measure reflects plant defense capacity. Since measures of lesion length cannot be used to unambiguously assess the tree's defense response to inoculation with bark beetle fungal associates, particularly as a single time point measurement (Wallin and Raffa 2001), we turned to investigate the impact of water limitation and G. clavigera inoculation on transcript abundance profiles for P contorta × banksiana genes putatively involved in responses to biotic and abiotic stress. Several genes profiled by qRT-PCR showed significant changes in transcript profiles in response to inoculation with G. clavigera, suggesting a role for these genes in the response of P. contorta × banksiana to attack by this MPB fungal associate. Given that P. contorta × banksiana is beyond the historic range of MPB, it is considered a naïve host to MPB and by extension to G. clavigera (Cudmore et al. 2010, Safranyik et al. 2010). Consequently, it is likely that the molecular defense responses examined here represent more of a generalized response to phloem- and sapwood-invading fungal pathogens than an evolutionarily co-evolved adaptive or specialized response by the tree to G. clavigera. Historical exposure of trees in the P. contorta–P. banksiana complex zone to other Ophiostomoid fungi vectored by secondary bark beetles such as Ips pini, Ips grandicollis and Dendroctonus valens may also be reflected in the response patterns. We observed that a subset of the defense-associated genes—PcbCHI2.1, PcbCHI2.2 and PcbTPS-3car1—were up-regulated under water deficit in the absence of inoculation, while another subset of defense-associated genes—PcbCHI1.1, PcbCHI4.1 and Pcb(E)- ± -farnesene synthase-like—showed attenuated expression upon challenge with G. clavigera under water deficit. These results are consistent with previously published reports in conifers that water limitation augments constitutive defenses but diminishes induced defenses (Lorio et al. 1995, Lombardero et al. 2000). They are also consistent with the findings of Lusebrink et al. (2011), who reported that water limitation resulted in a transient increase and subsequent decrease in total monoterpene levels, in both lodgepole and jack pine seedlings. Here, we use a finer-scale approach to demonstrate that both monoterpene-based and protein-based defenses are impacted by water limitation. The anatomical and cellular-level differences observed between G. clavigera-induced lesions of well-watered and water-deficit trees suggest that additional defense-associated metabolic pathways such as phenolic biosynthesis and cell wall biosynthesis and/or modification are also affected by water limitation. Analysis of PcTPS-3car1 demonstrated that this sequence encodes a functional 3-carene synthase (Hall et al. 2013b). 3-Carene synthases have a well-documented role in defense against herbivores in conifers (Fäldt et al. 2003, Robert et al. 2010, Hall et al. 2011). Interestingly, 3-carene makes up a greater proportion of the monoterpene profile of jack pine phloem than that of lodgepole pine (Lusebrink et al. 2011). Substrate specificity and product profiles for Pcb-(E)-β-farnesene-like have yet to be determined. Bohlmann et al. (1998) reported that a similar sesquiterpene synthase shows an inducible response upon wound treatments in Abies grandis. These TPS synthesize less of the oleoresin than of mono- or di-terpenes and are also involved in the transformation of terpenes into insect hormone mimics (Phillips and Croteau 1999). PcCHI1.1 encodes a Class I CHI, PcbCHI4.1 encodes a Class IV CHI, and PcbCHI2.1 and PcbCHI2.2 exhibit the hallmarks of Class II CHI. The Class I CHI, PcbCHI1.1, appears to contain the carboxy-terminal extension required for vacuolar localization (Neuhaus et al. 1991), while the Class IV CHI, PcbCHI4.1, is most likely localized to the apoplast. In general, apoplastic CHI show activity in early stages of the pathogenic infection to generate elicitors that increase the signal and activate defense mechanisms, while vacuolar CHI are activated later and to greater levels, acting to degrade chitin of penetrating fungal hyphae (Kasprzewska 2003). Chitinases have often been reported in defense responses of conifers to pathogens. For example, Class II and Class IV CHI have been implicated in defense responses of Norway spruce to Heterobasidion annosum (Hietala et al. 2004) and of Douglas-fir to Phellinus sulphurascens (Islam et al. 2010), while Class I, Class II and Class IV CHI have been described in defense of slash pine to Fusarium subglutinans f.sp. pini (Davis et al. 2002). In the cases of Norway spruce–H. annosum and slash pine–F. subglutinans f.sp. pini interactions, there is an early, lower intensity induction of Class II and Class IV CHI in resistant tree genotypes, whereas there is a late, higher intensity induction of the CHI in susceptible tree genotypes. In these studies, the authors concluded that the early response exhibited by resistant genotypes constituted part of an effective containment strategy, while the high magnitude late response of the susceptible genotypes likely represented a delayed and ultimately unsuccessful attempt to contain the invading pathogen (Davis et al. 2002, Hietala et al. 2004). Similarly, Fossdal et al. (2007) demonstrated that drought resulted in an earlier induction of some defense-related genes in Norway spruce seedlings inoculated with Rhizoctonia sp. These studies reinforce the importance of the timing of the defense response. Due to the scale of this field study with mature trees in a natural stand, we were not able to include a temporal component; however, future studies will incorporate a time course to address the rate and magnitude of induction for CHI and other defense-associated genes. Although originally selected as a candidate to assess abiotic (water limitation) stress, PcbTINY1-like was down-regulated in response to G. clavigera challenge, with greater decreases in transcript abundance exhibited by well-watered versus water- deficit trees. PcbTINY1-like exhibits greater similarity to the ERF subfamily of DREBs than to the DREB subfamily. While DREBs are usually induced by abiotic stresses and bind to DRE elements, ERFs are induced by biotic stresses and bind to ERE elements (Ohme-Takagi et al. 2000, Agarwal et al. 2006). AtTINY from A. thaliana binds to both DRE and ERE elements and is induced by abiotic stresses (drought and cold), and defense signaling molecules (ethylene and methyl jasmonate; Sun et al. 2008). Although PcbTINY1-like shows down-regulation rather than up-regulation in response to G. clavigera, the interaction observed between water deficit and G. clavigera inoculation hints that PcbTINY1-like could also play a role in mediating crosstalk between abiotic and biotic stress signaling. Interestingly, most of the G. clavigera-responsive genes showed similar transcript abundance in phloem tissues within lesions and phloem tissues distinct from but surrounding the lesions, suggestive of a systemic response. This is consistent with previous cellular-level analyses demonstrating that production of both phenolic- and terpenoid-based compounds occurs some distance away from the site of attack, and can be transported to the defense zone (Franceschi et al. 2000, 2005, Hudgins et al. 2003). Induction of defenses in the tissues immediately surrounding the lesion may act to further constrain fungal growth. In the context of an MPB attack, this corresponds to the beetle's entry point into the host tree. If the tree is able to rapidly perceive the presence of G. clavigera upon attack by the first arriving MPB, these data suggest that the tree could mount a systemic induction of defense responses along the length of the trunk in time to defend against G. clavigera vectored by MPB that join the mass attack. Such systemic acquired resistance has been demonstrated in Norway spruce inoculated with Ceratocystis polonica, a fungal associate of the bark beetle Ips typographus (Christiansen et al. 1999). If this is true, we would predict that pine populations sharing a co-evolutionary history with MPB would be likely to perceive MPB-vectored G. clavigera more rapidly and more specifically than evolutionarily naïve pine populations, and be able to mount a systemic response in time to render greater protection against a mass attack. Over longer time frames, if the tree were successfully able to defend against an MPB mass attack, it is possible that the tree may acquire systemic induced resistance, a greater level of basal resistance providing greater protection against attacks in subsequent years (Bonello et al. 2006). Systemic induced resistance, also referred to as priming, has been demonstrated in Monterey pine inoculated with the fungus causing pitch canker, Fusarium circinatum (Bonello et al. 1991), and in Norway spruce inoculated with the root and butt rot pathogen Heterobasidion parviporum (Fossdal et al. 2012). How does water deficit alter P. contorta × banksiana defense response profiles? Given that transcript abundance for genes belonging to both the TPS and CHI families of defense-associated genes are up-regulated under water-limiting conditions in the absence of G. clavigera inoculation, but that different genes belonging to the TPS and CHI families show attenuated expression under water-limiting conditions following G. clavigera inoculation, we infer that both carbon- and nitrogen-based defenses are impacted by water deficit. If the major effect of water limitation reduced photoassimilate availability for allocation to carbon-based defenses (Hodges and Lorio 1975, Lorio 1986, Christiansen et al. 1987, Byers 1995, Croisé et al. 2001), then we might have expected to observe a greater decrease in carbon-based defense gene expression than in nitrogen-based defense gene expression. This is not the case in this study. This finding is inconsistent with one of the predictions of the carbon-nutrient balance hypothesis (Bryant et al. 1983), which states that a change in resource conditions that decreases the carbon: nutrient ratio will lead to a decrease in non-nitrogen- (i.e., carbon-) based defenses, and a concomitant increase in nitrogen-based defenses. However, the finding is consistent with those of Guérard et al. (2007), who used radioisotope studies to demonstrate that carbon resources used in local defenses are derived mainly from locally stored carbon reserves rather than newly assimilated carbon. Another effect of water limitation is acceleration of the seasonal transition from active growth to growth cessation, i.e., dormancy (Lorio 1986, Gruber et al. 2010). Although we did not explicitly examine growth cessation in this study, there are lines of evidence to suggest that water limitation promoted an early transition to growth cessation, particularly in G. clavigera-inoculated water-deficit trees. First, we observed that the bark was consistently much more difficult to remove from water- deficit trees than well-watered trees, particularly if trees had been inoculated, suggesting that the number of cell layers making up the cambial zone and recent derivatives was reduced. This could not be examined directly in this study because the cambial zone was destroyed during sampling. However, the compressed appearance of sieve cells in phloem of water-deficit trees—exacerbated in inoculated trees—hints that these trees were progressing toward growth cessation more quickly. Second, PcbPIP1;1 was significantly down-regulated in G. clavigera-inoculated trees, and somewhat down-regulated in water-deficit control trees. A subset of PIPs mediate cell expansion in plant tissues (Maurel et al. 2002), suggesting that expression of PcbPIP1;1 decreased as cell expansion of developing cells was completed and cells proceeded with maturation. Third, the Class II CHIs PcbCHI2.1 and PcbCHI2.2 exhibited up-regulation under water-limited conditions regardless of inoculation treatment, while G. clavigera inoculation had little effect on transcript abundance. Conifers synthesize a diverse array of abiotic and biotic defense-associated proteins, including CHIs, during autumnal acclimation, concomitant with the transition to growth cessation and well in advance of autumn-like environmental conditions (Holliday et al. 2008, Jarzabek et al. 2009, Cooke et al. 2012, Galindo González et al. 2012). Seasonally accumulated CHIs, including those from conifers, have been shown to exhibit chitinolytic activity, antifreeze activity or both (e.g., Zamani et al. 2003, Griffith and Yaish 2004, Stressmann et al. 2004, Jarzabek et al. 2009). It is tempting to speculate that the effects of water deficit on patterns of constitutive defenses (Lombardero et al. 2000, Lusebrink et al. 2011) are due at least in part to contrasting phenologies between well-watered and water-deficit plants. Accelerated phenology of water-deficit plants may augment constitutive defenses as described above, providing some degree of generalized protection against pests and pathogens, while curtailing the capacity to mount inducible defenses. While seasonal effects on key conifer defenses such as oleoresin terpenoids have received relatively little attention, seasonal variation in oleoresin terpenoid composition has been documented in Pinus halepensis and Pinus pinea (Mita et al. 2002). The diminished capacity to mount an induced defense under water-limited conditions may in part be due to reduced photoassimilate availability that is a function of reduced photosynthesis as described above. In trees that are proceeding more quickly to growth cessation, we would also predict that carbon resources would be preferentially allocated in the stem to resource-intensive processes that accompany growth cessation such as cell wall biosynthesis and modification. We would further predict that carbon resources would be preferentially partitioned to perennating parts—particularly roots—for storage during the transition to growth cessation (Scarascia-Mugnozza et al. 1999), further reducing the availability of carbon skeletons for synthesis of defense-associated compounds. The notion that water deficit accelerates the transition from active growth to growth cessation, thereby influencing the array of constitutive and inducible defenses manifested by the tree, is consistent with the growth–differentiation hypothesis (Lorio 1986, Herms and Mattson 1992). The central tenet of the growth–differentiation hypothesis is that an environmental factor that reduces growth (i.e., cell division and cell expansion) more than photosynthesis will result in a shift in allocation priority, so that a greater proportion of carbon resources are allocated to differentiation-related processes. Differentiation-related processes comprise a myriad of cellular, biochemical and molecular changes that occur as cells reach maturation, including physical and chemical defenses (Lorio 1986). Lorio (1986) invoked the growth–differentiation hypothesis to explain how seasonally correlated developmental events and defense processes impacted loblolly pine–southern pine beetle interactions, and how water deficit could impact these interactions via promoting a shift from early to late season development in the stem, and the reductions in resin production that accompany this developmental transition. Here, we extend Lorio's model to the developmental transition from latewood production to cessation of cambial growth, and the numerous changes that occur to defense capacity as a result of this transition. While the evidence presented in this study is circumstantial for the effect of drought being primarily one of accelerating phenology to promote early growth cessation, with concomitant effects on constitutive and inducible defense, we propose this framework as a testable model to form the foundation for future experimentation. Conclusions Here, we identified a suite of genes from P. contorta × banksiana that respond to water deficit, G. clavigera inoculation or to both stressors in mature trees. The identification of G. clavigera-responsive genes provides us with the means to more accurately assess the defense responses of P. contorta × banksiana trees to an MPB attack than the traditionally used measure of lesion length. The expression profiles for this suite of genes demonstrated that, under water deficit, constitutive expression of some defense-associated genes was enhanced under water deficit, while the induction of other defense-associated genes upon G. clavigera inoculation was attenuated. That some genes respond to both stressors is suggestive of molecular crosstalk between water stress and defense responses in pines. The data also indicate that not all components of the constitutive and induced defenses are affected equally, suggesting that the effect of water deficit on the tree's complex defense machinery is not an ‘all-or-none’ response, but rather is quantitative in nature. The results suggest that if MPB and its fungal associates are able to overcome the constitutive defenses of a water-stressed tree, they are more likely to overcome the induced defenses. In recent years, MPB has undergone considerable range expansion across northern Alberta (Cullingham et al. 2011). The pines of this landscape are considered naïve hosts, and as such may not have acquired specialized defenses against MPB and their fungal associates. Consequently, the threshold density for a successful MPB attack may be lower for hosts in these novel habitats. Our results suggest that water deficit could further lower the critical threshold at which MPB successfully attacks pines in the region by reducing key aspects of the induced defense response. Given the recent periods of drought experienced in the region (Hanesiak et al. 2011), this environmental factor could serve as a positive feedback on MPB dynamics and consequently continued range expansion (Raffa et al. 2008). Furthermore, the region is predicted to receive less precipitation in the future under scenarios of climate change (Schindler and Donahaue 2006). If the gene expression pattern and lesion development differences reported here translate into differences in defense capability of novel hosts against an attack by MPB and their fungal associates, further investigation into how water deficit affects pine defense capacity against MPB and its pathogenic fungal associates will be important to ascertain stand susceptibility in the face of this unprecedented MPB epidemic that continues to expand into novel territory in concert with changing climate conditions. Supplementary data Supplementary data for this article are available at Tree Physiology Online. Conflict of interest None declared. Funding Funding for this research has been provided through grants from Genome Canada, The Government of Alberta through Genome Alberta, and The Government of British Columbia through Genome BC (The TRIA Project, http://www.thetriaproject.ca) to J.E.K.C. and B.J.C.; Alberta Innovates Bio Solutions to J.E.K.C. and B.J.C.; and The Canadian Foundation for Innovation to J.E.K.C. Acknowledgments We are grateful to Drs Jörg Bohlmann and Dawn Hall for access to TPS sequences and functional characterization prior to publication, helpful discussions and a critical reading of the manuscript. We acknowledge the British Columbia Genome Sciences Centre for sequencing, and Dr Christopher Keeling for sequence deposition at NCBI. We thank Tim McCready from Millar Western for providing access to the experimental field site at Whitecourt (AB), and personnel at Alberta Environment and Sustainable Resource Development and the Canadian Forest Service Northern Forestry Centre for assistance in selecting the field site and establishing the experiment. We also thank Dr Nadir Erbilgin, Bill Sperling, Jeremiah Bolstad and Andrew Ho (University of Alberta) for their valuable contributions to the field experiment. We thank Troy Locke (University of Alberta Department of Biological Sciences Molecular Biology Service Unit) for assistance with qRT-PCR and Bioanalyzer analyses, and Dr Patrick James (University of Alberta, Université de Montreal) for guidance with the principal component analysis. References Agarwal PK Agarwal P Reddy MK Sopory SK ( 2006 ) Role of DREB transcription factors in abiotic and biotic stress tolerance in plants . Plant Cell Rep 25 : 1263 – 1274 . Google Scholar Crossref Search ADS WorldCat Alexandersson E Fraysse L Sjövall-Larsen S Gustavsson S Fellert M Karlsson M Johanson U Kjellbom P ( 2005 ) Whole gene family expression and drought stress regulation of aquaporins . Plant Mol Biol 59 : 469 – 484 . Google Scholar Crossref Search ADS WorldCat Allen MD Yamasaki M Ohme-Takagi Tateno M Suzuki M ( 1998 ) A novel mode of DNA recognition by a β-sheet revealed by the solution structure of the GCC-box binding domain in complex with DNA . EMBO J 17 : 5484 – 5496 . Google Scholar Crossref Search ADS WorldCat Almeida-Rodriguez AM Cooke JE Yeh F Zwiazek JJ ( 2010 ) Functional characterization of drought-responsive aquaporins in Populus balsamifera and Populus simonii × balsamifera clones with different drought resistance strategies . Physiol Plant 14 : 321 – 333 . Google Scholar Crossref Search ADS WorldCat Andrews SF Flanagan LB Sharp EJ Cai T ( 2012 ) Variation in water potential, hydraulic characteristics and water source use in montane Douglas-fir and lodgepole pine trees in southwestern Alberta and consequences for seasonal changes in photosynthetic capacity . Tree Physiol 32 : 146 – 160 . Google Scholar Crossref Search ADS WorldCat Berryman AA ( 1982 ) Mountain pine beetle outbreaks in Rocky Mountain lodgepole pine forests . J Forest 80 : 410 – 413 . Google Scholar OpenURL Placeholder Text WorldCat Bishop JG Dean AM Mitchell-Olds T ( 2000 ) Rapid evolution in plant chitinases: molecular targets of selection in plant–pathogen coevolution . Proc Natl Acad Sci USA 97 : 5322 – 5327 . Google Scholar Crossref Search ADS WorldCat Bohlmann J Crock J Jetter R Croteau R ( 1998 ) Terpenoid-based defenses in conifers: cDNA cloning, characterization, and functional expression of wound-inducible (E)-alpha-bisabolene synthase from grand fir (Abies grandis) . Proc Natl Acad Sci USA 95 : 6756 – 6761 . Google Scholar Crossref Search ADS WorldCat Bonello P Gordon TR Storer AJ ( 1991 ) Systemic induced resistance in Monterey pine . For Pathol 31 : 99 – 106 . Google Scholar OpenURL Placeholder Text WorldCat Bonello P Gordon TR Herms DA Wood DL Erbilgin N ( 2006 ) Nature and ecological implications of pathogen-induced systemic resistance in conifers: a novel hypothesis . Physiol Mol Plant Pathol 68 : 95 – 104 . Google Scholar Crossref Search ADS WorldCat Boone CK Aukema BH Bohlmann J Carroll AL Raffa KF ( 2011 ) Efficacy of tree defense physiology varies with bark beetle population density: a basis for positive feedback in eruptive species . Can J For Res 41 : 1174 – 1188 . Google Scholar Crossref Search ADS WorldCat Brunner FA Stintzi A Fritig B Legrand M ( 1998 ) Substrate specificities of tobacco chitinases . Plant J 14 : 225 – 234 . Google Scholar Crossref Search ADS WorldCat Bryant JP Chapin FS III Klein DR ( 1983 ) Carbon/nutrient balance of boreal plants in relation to vertebrate herbivory . Oikos 40 : 357 – 368 . Google Scholar Crossref Search ADS WorldCat Byers JA ( 1995 ) Host-tree chemistry affecting colonization in bark beetles . In: Cardé RT Bell WJ (eds) Chemical Ecology of Insects 2 . Chapman and Hall , New York , pp 154 – 213 . Google Scholar Google Preview OpenURL Placeholder Text WorldCat COPAC Byun-McKay A Hunter WL Godard K-A Wang SX Martin DM Bohlmann J Plant AL ( 2003 ) Insect attack and wounding induce traumatic resin duct development and gene expression of (−)-pinene synthase in Sitka spruce . Plant Physiol 133 : 368 – 378 . Google Scholar Crossref Search ADS WorldCat Byun-McKay A Godard KA Toudefallah M Martin DM Alfaro R King J Bohlmann J Plant AL ( 2006 ) Wound-induced terpene synthase gene expression in Sitka spruce that exhibit resistance or susceptibility to attack by the white pine weevil . Plant Physiol 140 : 1009 – 1021 . Google Scholar Crossref Search ADS WorldCat Carroll AL Taylor SW Régnière J Safranyik L ( 2004 ) Effects of climate change on range expansion by the mountain pine beetle in British Columbia . In: Shore TL Brooks JE Stone JE (eds) Proceedings of the Mountain Pine Beetle Symposium: Challenges and Solutions . 30–31 October 2003. Kelowna, BC, Canada . Natural Resources Canada, Canadian Forest Service, Pacific Forestry Centre , Victoria, BC , Information Report BC-X-399, 298 , pp 223 – 232 . Google Scholar Google Preview OpenURL Placeholder Text WorldCat COPAC Chhin S Hogg EH Lieffers VJ Huang S ( 2008 ) Potential effects of climate change on the growth of lodgepole pine across diameter size classes and ecological regions . For Ecol Manag 256 : 1692 – 1703 . Google Scholar Crossref Search ADS WorldCat Christiansen E Waring RH Berryman AA ( 1987 ) Resistance of conifers to bark beetle attack: searching for general relationships . For Ecol Manag 22 : 89 – 106 . Google Scholar Crossref Search ADS WorldCat Christiansen E Krokene P Berryman AA Franceschi VR Krekling T Lieutier F Lönneborg A Solheim H ( 1999 ) Mechanical injury and fungal infection induce acquired resistance in Norway spruce . Tree Physiol 19 : 399 – 403 . Google Scholar Crossref Search ADS WorldCat Cooke JEK Eriksson ME Junttila O ( 2012 ) The dynamic nature of bud dormancy in trees; environmental control and molecular mechanisms . Plant Cell Environ 35 : 1707 – 1728 . Google Scholar Crossref Search ADS WorldCat Croisé L Lieutier F Cochard H Dreyer E ( 2001 ) Effects of drought stress and high density stem inoculations with Leptographium wingfieldii on hydraulic properties of young Scots pine trees . Tree Physiol 21 : 427 – 436 . Google Scholar Crossref Search ADS WorldCat Cudmore TJ Bjorklund N Carroll AL Lindgren BS ( 2010 ) Climate change and range expansion of an aggressive bark beetle: evidence of higher reproductive success in naive host tree populations . J Appl Ecol 47 : 1036 – 1043 . Google Scholar Crossref Search ADS WorldCat Cullingham CI Cooke JEK Dang S Davis CS Cooke BJ Coltman DW ( 2011 ) Mountain pine beetle host-range expansion threatens the boreal forest . Mol Ecol 20 : 2157 – 2171 . Google Scholar Crossref Search ADS WorldCat Cullingham CI James PMA Cooke JEK Coltman DW ( 2012 ) Characterizing the physical and genetic structure of the lodgepole pine × jack pine hybrid zone: mosaic structure and differential introgression . Evol Appl 5 : 879 – 891 . Google Scholar Crossref Search ADS WorldCat Davis JM Wu HG Cooke JEK Reed JM Luce KS Michler CH ( 2002 ) Pathogen challenge, salicylic acid, and jasmonic acid regulate expression of chitinase gene homologs in pine . Mol Plant Microbe Interact 15 : 380 – 387 . Google Scholar Crossref Search ADS WorldCat Domec JC Noorments A King JS Sun G McNulty SG Gavazzi MJ Boggs JL Treasure EA ( 2009 ) Decoupling the influence of leaf and root hydraulic conductances on stomatal conductance and its sensitivity to vapour pressure deficit as soil dries in a drained loblolly pine plantation . Plant Cell Environ 32 : 980 – 991 . Google Scholar Crossref Search ADS WorldCat Drummond AJ Ashton B Buxton S et al. ( 2010 ) Geneious v4.8.2 . http://www.geneious.com (13 November 2013, date last accessed) . OpenURL Placeholder Text WorldCat Dunn JP Lorio PL Jr ( 1993 ) Modified water regimes effect photosynthesis, xylem water potential, cambial growth, and resistance of juvenile Pinus taeda L. to Dendroctonus frontalis (Coleoptera: Scolytidae) . Environ Entomol 22 : 948 – 957 . Google Scholar Crossref Search ADS WorldCat El Kayal W Allen CC Ju CJ Adams E King-Jones S Zaharia LI Abrams SR Cooke JE ( 2011 ) Molecular events of apical bud formation in white spruce, Picea glauca . Plant Cell Environ 34 : 480 – 500 . Google Scholar Crossref Search ADS WorldCat Eyles A Bonello P Ganley E Mohammed C ( 2010 ) Induced resistance to pests and pathogens in trees . New Phytol 185 : 893 – 908 . Google Scholar Crossref Search ADS WorldCat Fäldt J Martin D Miller B Rawat S Bohlmann J ( 2003 ) Traumatic resin defense in Norway spruce (Picea abies): Methyl jasmonate-induced terpene synthase gene expression, and cDNA cloning and functional characterization of (+)-3-carene synthase . Plant Mol Biol 51 : 119 – 133 . Google Scholar Crossref Search ADS WorldCat Fossdal CG Hietala AM Kvaalen H Solheim H ( 2005 ) Changes in host chitinase isoforms in relation to wounding and colonization by Heterobasidion annosum: early and strong defense response in 33-year-old resistant Norway spruce clone . Tree Physiol 26 : 169 – 177 . Google Scholar Crossref Search ADS WorldCat Fossdal CG Nagy NE Johnsen Ø Dalen LS ( 2007 ) Local and systemic stress responses in Norway spruce: similarities in gene expression between a compatible pathogen interaction and drought stress . Physiol Mol Plant Pathol 70 : 161 – 173 . Google Scholar Crossref Search ADS WorldCat Fossdal CG Nagy NE Hietala AM Kvaalen H Slimestad R Woodward S Solheim H ( 2012 ) Indications of heightened constitutive or primed host response affecting the lignin pathway transcripts and phenolics in mature Norway spruce clones . Tree Physiol 32 : 1137 – 1147 . Google Scholar Crossref Search ADS WorldCat Franceschi VR Krokene P Krekling T Christiansen E ( 2000 ) Phloem parenchyma cells are involved in local and distant defense responses to fungal inoculation of bark-beetle attack in Norway spruce (Pinaceae) . Am J Bot 87 : 314 – 326 . Google Scholar Crossref Search ADS WorldCat Franceschi VR Krokene P Christiansen E Krekling T ( 2005 ) Anatomical and chemical defenses of conifer bark against bark beetles and other pests . New Phytol 167 : 353 – 376 . Google Scholar Crossref Search ADS WorldCat Galindo González LMG El Kayal W Ju CJ-T Allen CCG King-Jones S Cooke JEK ( 2012 ) Integrated transcriptomic and proteomic profiling of white spruce stems during the transition from active growth to dormancy . Plant Cell Environ 35 : 682 – 701 . Google Scholar Crossref Search ADS WorldCat Government of Alberta . Why are they here? http://mpb.alberta.ca/Resources/documents/MPB-WhyAreTheyHere.pdf (9 April 2013, date last accessed) . Google Scholar OpenURL Placeholder Text WorldCat Government of British Columbia . Facts about B.C.'s mountain pine beetle . http://www.for.gov.bc.ca/hfp/mountain_pine_beetle/Updated-Beetle-Facts_May2012.pdf (9 April 2013, date last accessed) . Google Scholar OpenURL Placeholder Text WorldCat Griffith M Yaish MWF ( 2004 ) Antifreeze proteins in overwintering plants: a tale of two activities . Trends Plant Sci 9 : 399 – 405 . Google Scholar Crossref Search ADS WorldCat Gruber A Strobl S Veit B Oberhuber W ( 2010 ) Impact of drought on the temporal dynamics of wood formation in Pinus sylvestris . Tree Physiol 30 : 490 – 501 . Google Scholar Crossref Search ADS WorldCat Guérard N Maillard P Bréchet C Lieutier F Dreyer E ( 2007 ) Do trees use reserve or newly assimilated carbon for their defense reactions? A 13C labeling approach with young Scots pines inoculated with a bark-beetle-associated fungus (Ophiostoma brunneo ciliatum) . Ann For Sci 64 : 601 – 608 . Google Scholar Crossref Search ADS WorldCat Hall D Robert JA Keeling CI et al. ( 2011 ) An integrated genomic, proteomic and biochemical analysis of (+)-3-carene biosynthesis in Sitka spruce (Picea sitchensis) genotypes that are resistant or susceptible to white pine weevil . Plant J 65 : 936 – 948 . Google Scholar Crossref Search ADS WorldCat Hall DE Zerbe P Jancsik S Lara Quesada A Dullat H Madilao LL Yuen M Bohlman J ( 2013a) Evolution of conifer diterpene synthases: diterpene resin acid biosynthesis in lodgepole pine and jack pine involves monofunctional and bifunctional diterpene synthases . Plant Physiol 162 : 660 – 616 . Google Scholar OpenURL Placeholder Text WorldCat Hall DE Yuen MMS Jancsik S et al. ( 2013b) Transcriptome resources and functional characterization of monoterpene synthases for two host species of the mountain pine beetle, lodgepole pine (Pinus contorta) and jack pine (Pinus banksiana) . BMC Plant Biol 13 : 80 . Google Scholar Crossref Search ADS WorldCat Hamel F Boivin R Tremblay C Bellemare G ( 1997 ) Structural and evolutionary relationships among chitinases of flowering plants . J Mol Evol 44 : 614 – 624 . Google Scholar Crossref Search ADS WorldCat Hanesiak JM Stewart RE Bonsal BR et al. ( 2011 ) Characterization and summary of the 1999–2005 Canadian Prairie drought . Atmos-Ocean 49 : 421 – 452 . Google Scholar Crossref Search ADS WorldCat Herms DA Mattson WJ ( 1992 ) The dilemma of plants: to grow or defend . Q Rev Biol 67 : 283 – 335 . Google Scholar Crossref Search ADS WorldCat Hietala AM Kvaalen H Schmidt A Jøhnk N Solheim H Fossdal CG ( 2004 ) Temporal and spatial profiles of chitinase expression by Norway spruce in response to bark colonization by Heterobasidion annosum . Appl Environ Microbiol 70 : 3948 – 3953 . Google Scholar Crossref Search ADS WorldCat Hodges JD Lorio PL Jr ( 1975 ) Moisture stress and composition of xylem oleoresin in loblolly pine . For Sci 21 : 283 – 290 . Google Scholar OpenURL Placeholder Text WorldCat Hogg EH ( 1997 ) Temporal scaling of moisture and the forest-grassland boundary in western Canada . Agric For Meteorol 84 : 115 – 122 . Google Scholar Crossref Search ADS WorldCat Holliday JA Ralph SG White R Bohlmann J Aitken SN ( 2008 ) Global monitoring of autumn gene expression within and among phenotypically divergent populations of Sitka spruce (Picea sitchensis) . New Phytol 178 : 103 – 122 . Google Scholar Crossref Search ADS WorldCat Huber DPW Ralph S Bohlmann J ( 2004 ) Genomic hardwiring and phenotypic plasticity of terpenoid-based defenses in conifers . J Chem Ecol 30 : 2399 – 2418 . Google Scholar Crossref Search ADS WorldCat Hudgins JW Christiansen E Franceschi VR ( 2003 ) Methyl jasmonate induce changes mimicking anatomical defenses in diverse members of the Pinaceae . Tree Physiol 23 : 361 – 367 . Google Scholar Crossref Search ADS WorldCat Islam MA Sturrock RN Williams HL Ekramoddoullah KM ( 2010 ) Identification, characterization, and expression analyses of class II and IV chitinase genes from Douglas-fir seedlings infected by Phellinus sulphurascens . Phytopathology 100 : 356 – 366 . Google Scholar Crossref Search ADS WorldCat Islam MA Sturrock RN Ekramoddoullah KM ( 2011 ) Conifer chitinases . Am J Plant Sci Biotech 5 : 22 – 26 . Google Scholar OpenURL Placeholder Text WorldCat Jarzabek M Pukacki PM Nuc K ( 2009 ) Cold-regulated proteins with potent antifreeze and cryoprotective activities in spruces (Picea spp.) . Cryobiology 58 : 268 – 274 . Google Scholar Crossref Search ADS WorldCat Kaldenhoff R Ribas-Carbo M Sans JF Lovisolo C Hechwolf M Uehlein N ( 2008 ) Aquaporins and plant water balance . Plant Cell Environ 31 : 658 – 666 . Google Scholar Crossref Search ADS WorldCat Kasprzewska A ( 2003 ) Plant chitinase—regulation and function . Cell Mol Biol Lett 8 : 809 – 824 . Google Scholar OpenURL Placeholder Text WorldCat Keeling CI Bohlmann J ( 2006 ) Genes, enzymes and chemicals of terpenoid diversity in the constitutive and induced defense of conifers against insects and pathogens . New Phytol 170 : 657 – 675 . Google Scholar Crossref Search ADS WorldCat Klein T Cohen S Yakir D ( 2011 ) Hydraulic adjustments underlying drought resistance of Pinus halepensis . Tree Physiol 31 : 637 – 648 . Google Scholar Crossref Search ADS WorldCat Kolb TE Homberg KM Wagner MR Stone JE ( 1998 ) Regulation of ponderosa pine foliar physiology and insect resistance mechanisms by basal area treatments . Tree Physiol 18 : 375 – 381 . Google Scholar Crossref Search ADS WorldCat Lee S Kim JJ Breuil C ( 2006 ) Fungal diversity associated with the mountain pine beetle, Dendroctonus ponderosae and infested lodgepole pines in British Columbia . Fungal Divers 22 : 91 – 105 . Google Scholar OpenURL Placeholder Text WorldCat Lippert D Chowrira S Ralph SG Zhuang J Aeschliman D Ritland C Rigland K Bohlmann J ( 2007 ) Conifer defense against insects: proteome analysis of Sitka spruce (Picea sitchensis) bark induced by mechanical wounding or feeding by white pine weevils (Pissodes strobi) . Proteomics 7 : 248 – 270 . Google Scholar Crossref Search ADS WorldCat Liu Y Zhao TJ Liu JM Liu Q Yan YB Zhoue HM ( 2006 ) The conserved Ala37 in the ERF/AP2 domain is essential for binding with the DRE element and the GCC box . FEBS Lett 580 : 1303 – 1308 . Google Scholar Crossref Search ADS WorldCat Lombardero MH Ayres MP Lorio PL Jr Ruel JJ ( 2000 ) Environmental effects on constitutive and inducible resin defences of Pinus taeda . Ecol Lett 3 : 329 – 339 . Google Scholar Crossref Search ADS WorldCat Lorio PL Jr ( 1986 ) Growth and differentiation balance: a basis for understanding southern pine beetle Dendroctonus frontalis and tree interactions . For Ecol Manag 14 : 259 – 274 . Google Scholar Crossref Search ADS WorldCat Lorio PL Stephen FM Paine TD ( 1995 ) Environment and ontogeny modify loblolly pine response to induced acute water deficits and bark beetle attack . For Ecol Manag 3 : 97 – l10 . Google Scholar Crossref Search ADS WorldCat Lu M Wingfield MJ Gillette NE Mori SR Sun JH ( 2010 ) Complex interactions among host pines and fungi vectored by an invasive bark beetle . New Phytol 187 : 859 – 866 . Google Scholar Crossref Search ADS WorldCat Lusebrink I Evenden MI Blanchet FG Cooke JEK Erbilgin N ( 2011 ) Effect of water stress and fungal inoculation on monoterpene émission from an historical and a new pine host of the mountain pine beetle . J Chem Ecol 37 : 1013 – 1026 . Google Scholar Crossref Search ADS WorldCat Man G ( 2012 ) Major forest insect and disease conditions in the United States: 2011 update. FS-1000 . United States Department of Agriculture Forest Service , Washington, DC , 50 pp. Google Scholar Maurel CH Javot V Lauvergeat P Gerbeau C Tournaire V Santoni V Heyes J ( 2002 ) Molecular physiology of aquaporins in plants . Int Rev Cytol 215 : 105 – 148 . Google Scholar OpenURL Placeholder Text WorldCat McDowell N Pockman WT Allen CD et al. ( 2008 ) Mechanisms of plant survival and mortality during drought: why do some plants survive while others succumb to drought? New Phytol 178 : 719 – 739 . Google Scholar Crossref Search ADS WorldCat Michaelian M Hogg EH Hall RJ Arsenault E ( 2011 ) Massive mortality of aspen following severe drought along the southern edge of the Canadian boreal forest . Glob Change Biol 17 : 2084 – 2094 . Google Scholar Crossref Search ADS WorldCat Miller B Madilao LL Ralph S Bohlmann J ( 2005 ) Insect-induced conifer defense. White pine weevil and methyl jasmonate induce traumatic resinosis, de novo formed volatile emissions, and accumulation of terpenoid synthase and putative octadecanoid pathway transcripts in Sitka spruce . Plant Physiol 127 : 369 – 382 . Google Scholar Crossref Search ADS WorldCat Miranda M Ralph SG Mellway R White R Heath MC Bohlmann J Constabel CP ( 2007 ) The transcriptional response of hybrid poplar (Populus trichocarpa × P. deltoides) to infection by Melampsora medusae Leaf rust involves induction of flavonoid pathway genes leading to the accumulation of proanthocyanidins . Mol Plant Microbe Interact 20 : 816 – 831 . Google Scholar Crossref Search ADS WorldCat Mita E Tsitsimpikou C Tsiveleka L Petrakis PV Ortiz A Vagias C Roussis V ( 2002 ) Seasonal variation of oleoresin terpenoids from Pinus halepensis and Pinus pinea and host selection of the scale insect Marchalina hellenica (Homoptera, Coccoidea, Margarodidae, Coelostonidiinae) . Holzforschung 56 : 572 – 578 . Google Scholar Crossref Search ADS WorldCat Neuhaus JM Sticher L Meins F Boller T ( 1991 ) A short C-terminal sequence is necessary and sufficient for the targeting of chitinases to the plant vacuole . Proc Natl Acad Sci USA 88 : 10362 – 10366 . Google Scholar Crossref Search ADS WorldCat Ohme-Takagi M Suzuki K Shinshi H ( 2000 ) Regulation of ethylene-induced transcription of defense genes . Plant Cell Physiol 41 : 1187 – 1192 . Google Scholar Crossref Search ADS WorldCat Ozturk ZN Talame V Dehyolos M Michalowski CB Galbraith DW Gozukimizi N Tuberosa R Bohnert HJ ( 2002 ) Monitoring large-scale changes in transcript abundance in drought- and salt-stressed barley . Plant Mol Biol 48 : 551 – 573 . Google Scholar Crossref Search ADS WorldCat Paine TD Stephen FM ( 1987 ) Influence of tree stress and site quality on the induced defense system of loblolly pine . Can J For Res 17 : 569 – 571 . Google Scholar Crossref Search ADS WorldCat Park W Scheffler BE Bauer PJ Campbell BT ( 2010 ) Identification of the family of aquaporin genes and their expression in upland cotton (Gossypium hirsutum L.) . BMC Plant Biol 10 : 142 . Google Scholar Crossref Search ADS WorldCat Parmeter JR Slaughter GW Chen M Wood DL ( 1992 ) Rate and depth of sapwood occlusion following inoculation of pines with bluestain fungi . For Sci 38 : 34 – 44 . Google Scholar OpenURL Placeholder Text WorldCat Pavy N Boyle B Nelson C et al. ( 2008 ) Identification of conserved core xylem gene stets: conifer cDNA microarray development, transcript profiling and computational analyses . New Phytol 180 : 766 – 786 . Google Scholar Crossref Search ADS WorldCat Phillips MA Croteau RB ( 1999 ) Resin-based defenses in conifers . Trends Plant Sci 4 : 184 – 190 . Google Scholar Crossref Search ADS WorldCat Plattner A Kim JJ DiGuistini S Breuil C ( 2008 ) Variation in pathogenicity of a mountain pine beetle-associated blue-stain fungus, Grosmannia clavigera, on young lodgepole pine in British Columbia . Can J Plant Pathol 30 : 457 – 466 . Google Scholar Crossref Search ADS WorldCat Raffa KF Berryman AA ( 1983 ) Physiological aspects of lodgepole pine wound response to a fungal symbiont of the mountain pine beetle . Can Entomol 115 : 723 – 734 . Google Scholar Crossref Search ADS WorldCat Raffa KF Smalley EB ( 1995 ) Interaction of pre-attack and induced monoterpene concentration in host conifer defense against bark beetle–fungal complexes . Oecologia 102 : 285 – 295 . Google Scholar Crossref Search ADS WorldCat Raffa KF Aukema BH Erbilgin N Klepzig KD Wallin KF ( 2005 ) Interactions among conifer terpenoids and bark beetles across multiple levels of scale: an attempt to understand links between population patterns and physiological processes . Recent Adv Phytochem 39 : 79 – 118 . Google Scholar OpenURL Placeholder Text WorldCat Raffa KF Aukema BH Bentz BJ Carroll AL Hicke JA Turner MG Romme WH ( 2008 ) Cross-scale drivers of natural disturbances prone to anthropogenic amplification: the dynamics of bark beetle eruptions . BioScience 58 : 501 – 517 . Google Scholar Crossref Search ADS WorldCat Ralph SG Yueh H Friedmann M et al. ( 2006 ) Conifer defence against insects: microarray gene expression profiling of Sitka spruce (Picea sitchensis) induced by mechanical wounding or feeding by spruce budworms (Choristoneura occidentalis) or white pine weevils (Pissodes strobi) reveals large-scale changes of the host transcriptome . Plant Cell Environ 29 : 1545 – 1570 . Google Scholar Crossref Search ADS WorldCat R Development Core Team ( 2010 ) R: A Language and Environment for Statistical Computing (http://www.R-project.org) . R Foundation for Statistical Computing , Vienna, Austria . Google Scholar Régnière J ( 1996 ) Generalized approach to landscape-wide seasonal forecasting with temperature-driven simulation models . Environ Entomol 25 : 869 – 881 . Google Scholar Crossref Search ADS WorldCat Rice AV Langor DW ( 2008 ) A comparison of the heat pulse velocity and lesion lengths for assessing the relative virulence of mountain pine beetle-associated fungi on jack pine . For Pathol 38 : 257 – 262 . Google Scholar OpenURL Placeholder Text WorldCat Rice AV Thormann MN Langor WD ( 2007 ) Virulence of, and interactions among mountain pine beetle associated blue-stain on two pine species and their hybrids in Alberta . Can J Bot 85 : 316 – 323 . Google Scholar Crossref Search ADS WorldCat Richardson KC Jarrett L Finke EH ( 1960 ) Embedding in epoxy resins for ultrathin sectioning in electron microscopy . Stain Technol 35 : 313 – 323 . Google Scholar Crossref Search ADS WorldCat Robert JA Madilao LL White R Yanchuk A King J Bohlmann J ( 2010 ) Terpenoid metabolite profiling in Sitka spruce identifies association of dehydroabietic acid, (+)-3-carene and terpinolene with resistance against white pine weevil . Botany 88 : 810 – 820 . Google Scholar Crossref Search ADS WorldCat Robinson DA Jones SB Wraith JM Or D Friedman SP ( 2003 ) A review of advances in dielectric and electrical conductivity measurement in soils using time domain reflectometry . Vadose Zone J 2 : 444 – 475 . Google Scholar Crossref Search ADS WorldCat Roe AD Rice AV Bromilow S Cooke JEK Sperling FAH ( 2010 ) Multilocus species identification and fungal DNA barcoding: insights from blue stain fungal symbionts of the mountain pine beetle . Mol Ecol Res 10 : 946 – 959 . Google Scholar Crossref Search ADS WorldCat Roe AD Rice AV Coltman D Cooke JEK Sperling FAH ( 2011 ) Comparative phylogeography, genetic differentiation, and contrasting reproductive modes in three fungal symbionts of a multipartite bark beetle symbiosis . Mol Ecol 20 : 584 – 600 . Google Scholar Crossref Search ADS WorldCat Safranyik L Carroll A ( 2006 ) The biology and epidemiology of the mountain pine beetle in lodgepole pine forests . In: Safranyik L Wilson B (eds) The mountain pine beetle, a synthesis of biology, management, and impacts on lodgepole pine . Natural Resources Canada, Canadian Forest Service, Pacific Forestry Service , Victoria, BC , pp 3 – 66 . Google Scholar Google Preview OpenURL Placeholder Text WorldCat COPAC Safranyik L Carroll AL Régnière J et al. ( 2010 ) Potential for range expansion of mountain pine beetle into the boreal forest of North America . Can Entomol 142 : 415 – 442 . Google Scholar Crossref Search ADS WorldCat Sakuma Y Liu Q Dubouzet JG Abe H Shinozaki K Yamaguchi-Shinozaki K ( 2002 ) DNA-binding specificity of the ERF/AP2 domain of Arabidopsis DREBs, transcription factors involved in dehydration- and cold-inducible gene expression . Biochem Biophys Res Commun 290 : 998 – 1009 . Google Scholar Crossref Search ADS WorldCat Scarascia-Mugnozza GE Hinckley TM Stettler RF Heilman PE Isebrands JG ( 1999 ) Production physiology and morphology of Populus species and their hybrids grown under short rotation. III. Seasonal carbon allocation patterns from branches . Can J For Res 29 : 1419 – 1432 . Google Scholar Crossref Search ADS WorldCat Schindler DW Donahaue WF ( 2006 ) An impending water crisis in Canada's western prairie provinces . Proc Natl Acad Sci USA 103 : 7210 – 7216 . Google Scholar Crossref Search ADS WorldCat Shelden MC Howitt MSM Kaiser BN Tyerman SD ( 2009 ) Identification and functional characterization of aquaporins in the grapevine, Vitis vinifera . Funct Plant Biol 36 : 1065 – 1078 . Google Scholar Crossref Search ADS WorldCat Six DL Wingfield MJ ( 2011 ) The role of phytopathogenicity in bark beetle–fungus symbioses: a challenge to the classic paradigm . Annu Rev Entomol 56 : 255 – 272 . Google Scholar Crossref Search ADS WorldCat Stressmann M Kitao S Griffith M Moresoli C Bravo LA Marangoni AG ( 2004 ) Calcium interacts with antifreeze proteins and chitinase from cold-acclimated winter rye . Plant Physiol 135 : 364 – 376 . Google Scholar Crossref Search ADS WorldCat Sun S Yu J-P Chen F Zhao T-J Fang X-H Li Y-Q Sui SF ( 2008 ) TINY, a dehydration-responsive element (DRE)-binding protein-like transcription factor connecting the DRE and ethylene-responsive element-mediated signaling pathways in Arabidopsis . J Biol Chem 283 : 6261 – 6271 . Google Scholar Crossref Search ADS WorldCat Tamura K Dudley J Nei M Kumar S ( 2007 ) MEGA4: molecular evolutionary genetics analysis (MEGA) software version 4.0 . Mol Biol Evol 24 : 1596 – 1599 . Google Scholar Crossref Search ADS WorldCat Topp GC Davis JL Annan AP ( 1980 ) Electromagnetic determination of soil water content: measurements in coaxial transmission lines . Water Resour Res 16 : 574 – 582 . Google Scholar Crossref Search ADS WorldCat Trenberth KE Jones PD Ambenje P et al. ( 2007 ) Observations: surface and atmospheric climate change . In: Solomon S Qin D Manning M Chen Z Marquis M Averyt KB Tignor M Miller HL (eds) Climate change 2007: the physical science basis. Contribution of Working Group I to the Fourth Assessment Report of the Intergovernmental Panel on Climate Change . Cambridge University Press , Cambridge, UK, and New York, NY, USA , pp 235 – 336 . Google Scholar Google Preview OpenURL Placeholder Text WorldCat COPAC Van Loon LC Rep M Pieterse CMJ ( 2006 ) Significance of inducible defense-related proteins in infected plants . Annu Rev Phytopathol 44 : 135 – 162 . Google Scholar Crossref Search ADS WorldCat Wallin KF Raffa KF ( 2001 ) Effects of folivory on subcortical plant defenses: can defense theories predict interguild processes? Ecology 82 : 1387 – 1400 . Google Scholar Crossref Search ADS WorldCat Wallin KF Raffa KF ( 2002 ) Density-mediated responses of bark beetles to host allelochemicals: a link between individual behavior and population dynamics . Ecol Entomol 27 : 484 – 492 . Google Scholar Crossref Search ADS WorldCat Witzell J Martín JA ( 2010 ) Phenolic metabolites in the resistance of northern forest trees to pathogens—past experiences and future prospects . Can J For Res 38 : 711 – 2727 . Google Scholar OpenURL Placeholder Text WorldCat Wong BL Berryman AA ( 1977 ) Resistance to the fir engraver beetle. 3. Lesion development and containment of infection by resistant Abies grandis inoculated with Trischosporium symbioticum . Can J Bot 55 : 2358 – 2365 . Google Scholar Crossref Search ADS WorldCat Xu F Fan C He Y ( 2007 ) Chitinases in Oryza sativa ssp. japonica and Arabidopsis thaliana . J Genet Genomics 34 : 138 – 150 . Google Scholar Crossref Search ADS WorldCat Yamaoka Y Hiratsuka Y Maruyama PJ ( 1995 ) The ability of Ophiostoma clavigerum to kill mature lodgepole pine trees . Eur J Forest Pathol 25 : 401 – 404 . Google Scholar Crossref Search ADS WorldCat Zamani A Sturrock R Ekramoddoullah AKM Wiseman SB Griffith M ( 2003 ) Endochitinase activity in the apoplastic fluid of Phellinus weirii-infected Douglas-fir and its association with over wintering and antifreeze activity . For Pathol 33 : 299 – 316 . Google Scholar OpenURL Placeholder Text WorldCat Zhang Y Wang Y Jiang L Xu Y Wang Y Lu D Chen F ( 2007 ) Aquaporin JcPIP2 is involved in drought responses in Jatropha curcas . Acta Biochem Biophys Sin 39 : 787 – 794 . Google Scholar Crossref Search ADS WorldCat Zulak KG Bohlmann J ( 2010 ) Terpenoid biosynthesis and specialized vascular cells of conifer defense . J Integr Plant Biol 52 : 86 – 97 . Google Scholar Crossref Search ADS WorldCat Zulak KG Lippert DN Kuzyk M Domanski D Chou T Borchers CH Bohlmann J ( 2009 ) Targeted proteomics using selected reaction monitoring (SRM) reveals the induction of specific terpene synthases in a multi-level study of methyl jasmonate treated Norway spruce (Picea abies) . Plant J 60 : 1015 – 1030 . Google Scholar Crossref Search ADS WorldCat © The Author 2013. Published by Oxford University Press. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/3.0/), which permits unrestricted reuse, distribution, and reproduction in any medium, provided the original work is properly cited. © The Author 2013. Published by Oxford University Press.
TI - Influence of water deficit on the molecular responses of Pinus contorta × Pinus banksiana mature trees to infection by the mountain pine beetle fungal associate, Grosmannia clavigera
JF - Tree Physiology
DO - 10.1093/treephys/tpt101
DA - 2014-11-01
UR - https://www.deepdyve.com/lp/oxford-university-press/influence-of-water-deficit-on-the-molecular-responses-of-pinus-bz9K89YqnL
SP - 1220
EP - 1239
VL - 34
IS - 11
DP - DeepDyve
ER -