Association genetics of coastal Douglas fir (Pseudotsuga menziesii var. menziesii, Pinaceae). I. Cold-hardiness related traits

Adaptation to cold is one of the greatest challenges to forest trees. This process is highly synchronized with environmental cues relating to photoperiod and temperature. Here, we use a candidate gene-based approach to search for genetic associations between 384 single-nucleotide polymorphism (SNP) markers from 117 candidate genes and 21 cold-hardiness related traits. A general linear model approach, including population structure estimates as covariates, was implemented for each marker-trait pair. We discovered 30 highly significant genetic associations [false discovery rate (FDR) Q < 0.10] across 12 candidate genes and 10 of the 21 traits. We also detected a set of 7 markers that had elevated levels of differentiation between sampling sites situated across the Cascade crest in northeastern Washington. Marker effects were small (r(2) < 0.05) and within the range of those published previously for forest trees. The derived SNP allele, as measured by a comparison to a recently diverged sister species, typically affected the phenotype in a way consistent with cold hardiness. The majority of markers were characterized as having largely nonadditive modes of gene action, especially underdominance in the case of cold-tolerance related phenotypes. We place these results in the context of trade-offs between the abilities to grow longer and to avoid fall cold damage, as well as putative epigenetic effects. These associations provide insight into the genetic components of complex traits in coastal Douglas fir, as well as highlight the need for landscape genetic approaches to the detection of adaptive genetic diversity.

Evolution of class III homeodomain-leucine zipper genes in streptophytes

Land plants underwent tremendous evolutionary change following the divergence of the ancestral lineage from algal relatives. Several important developmental innovations appeared as the embryophyte clade diversified, leading to the appearance of new organs and tissue types. To understand how these changes came about, we need to identify the fundamental genetic developmental programs that are responsible for growth, patterning, and differentiation and describe how these programs were modified and elaborated through time to produce novel morphologies. Class III homeodomain-leucine zipper (class III HD-Zip) genes, identified in the model plant Arabidopsis thaliana, provide good candidates for basic land plant patterning genes. We show that these genes may have evolved in a common ancestor of land plants and their algal sister group and that the gene family has diversified as land plant lineages have diversified. Phylogenetic analysis, expression data from nonflowering lineages, and evidence from Arabidopsis and other flowering plants indicate that class III HD-Zip genes acquired new functions in sporophyte apical growth, vascular patterning and differentiation, and leaf development. Modification of expression patterns that accompanied diversification of class III HD-Zip genes likely played an important role in the evolution of land plant form.

Architecture of the wood-wide web: Rhizopogon spp. genets link multiple Douglas-fir cohorts

*The role of mycorrhizal networks in forest dynamics is poorly understood because of the elusiveness of their spatial structure. We mapped the belowground distribution of the fungi Rhizopogon vesiculosus and Rhizopogon vinicolor and interior Douglas-fir trees (Pseudotsuga menziesii var. glauca) to determine the architecture of a mycorrhizal network in a multi-aged old-growth forest. *Rhizopogon spp. mycorrhizas were collected within a 30 x 30 m plot. Trees and fungal genets were identified using multi-locus microsatellite DNA analysis. Tree genotypes from mycorrhizas were matched to reference trees aboveground. Two trees were considered linked if they shared the same fungal genet(s). *The two Rhizopogon species each formed 13-14 genets, each colonizing up to 19 trees in the plot. Rhizopogon vesiculosus genets were larger, occurred at greater depths, and linked more trees than genets of R. vinicolor. Multiple tree cohorts were linked, with young saplings established within the mycorrhizal network of Douglas-fir veterans. A strong positive relationship was found between tree size and connectivity, resulting in a scale-free network architecture with small-world properties. *This mycorrhizal network architecture suggests an efficient and robust network, where large trees play a foundational role in facilitating conspecific regeneration and stabilizing the ecosystem.

Nucleotide diversity and linkage disequilibrium in cold-hardiness- and wood quality-related candidate genes in Douglas fir

Nuclear sequence variation and linkage disequilibrium (LD) were studied in 15 cold-hardiness- and 3 wood quality-related candidate genes in Douglas fir [Pseudotsuga menziesii (Mirb.) Franco]. This set of genes was selected on the basis of its function in other plants and collocation with cold-hardiness-related quantitative trait loci (QTL). The single-nucleotide polymorphism (SNP) discovery panel represented 24 different trees from six regions in Washington and Oregon plus parents of a segregating population used in the QTL study. The frequency of SNPs was one SNP per 46 bp across coding and noncoding regions on average. Haplotype and nucleotide diversities were also moderately high with H(d) = 0.827 +/- 0.043 and pi = 0.00655 +/- 0.00082 on average, respectively. The nonsynonymous (replacement) nucleotide substitutions were almost five times less frequent than synonymous ones and substitutions in noncoding regions. LD decayed relatively slowly but steadily within genes. Haploblock analysis was used to define haplotype tag SNPs (htSNPs). These data will help to select SNPs for association mapping, which is already in progress.

Genecology of Douglas fir in western Oregon and Washington

BACKGROUND AND AIMS

Genecological knowledge is important for understanding evolutionary processes and for managing genetic resources. Previous studies of coastal Douglas fir (Pseudotsuga menziesii var. menziesii) have been inconclusive with respect to geographical patterns of variation, due in part to limited sample intensity and geographical and climatic representation. This study describes and maps patterns of genetic variation in adaptive traits in coastal Douglas fir in western Oregon and Washington, USA.

METHODS

Traits of growth, phenology and partitioning were measured in seedlings of 1338 parents from 1048 locations grown in common gardens. Relations between traits and environments of seed sources were explored using regressions and canonical correlation analysis. Maps of genetic variation as related to the environment were developed using a geographical information system (GIS).

KEY RESULTS

Populations differed considerably for adaptive traits, in particular for bud phenology and emergence. Variation in bud-set, emergence and growth was strongly related to elevation and cool-season temperatures. Variation in bud-burst and partitioning to stem diameter versus height was related to latitude and summer drought. Seedlings from the east side of the Washington Cascades were considerably smaller, set bud later and burst bud earlier than populations from the west side.

CONCLUSIONS

Winter temperatures and frost dates are of overriding importance to the adaptation of Douglas fir to Pacific Northwest environments. Summer drought is of less importance. Maps generated using canonical correlation analysis and GIS allow easy visualization of a complex array of traits as related to a complex array of environments. The composite traits derived from canonical correlation analysis show two different patterns of variation associated with different gradients of cool-season temperatures and summer drought. The difference in growth and phenology between the westside and eastside Washington Cascades is hypothesized to be a consequence of the presence of interior variety (P. menziessii var. glauca) on the eastside.

A SNP resource for Douglas-fir: de novo transcriptome assembly and SNP detection and validation

BACKGROUND

Douglas-fir (Pseudotsuga menziesii), one of the most economically and ecologically important tree species in the world, also has one of the largest tree breeding programs. Although the coastal and interior varieties of Douglas-fir (vars. menziesii and glauca) are native to North America, the coastal variety is also widely planted for timber production in Europe, New Zealand, Australia, and Chile. Our main goal was to develop a SNP resource large enough to facilitate genomic selection in Douglas-fir breeding programs. To accomplish this, we developed a 454-based reference transcriptome for coastal Douglas-fir, annotated and evaluated the quality of the reference, identified putative SNPs, and then validated a sample of those SNPs using the Illumina Infinium genotyping platform.

RESULTS

We assembled a reference transcriptome consisting of 25,002 isogroups (unique gene models) and 102,623 singletons from 2.76 million 454 and Sanger cDNA sequences from coastal Douglas-fir. We identified 278,979 unique SNPs by mapping the 454 and Sanger sequences to the reference, and by mapping four datasets of Illumina cDNA sequences from multiple seed sources, genotypes, and tissues. The Illumina datasets represented coastal Douglas-fir (64.00 and 13.41 million reads), interior Douglas-fir (80.45 million reads), and a Yakima population similar to interior Douglas-fir (8.99 million reads). We assayed 8067 SNPs on 260 trees using an Illumina Infinium SNP genotyping array. Of these SNPs, 5847 (72.5%) were called successfully and were polymorphic.

CONCLUSIONS

Based on our validation efficiency, our SNP database may contain as many as ~200,000 true SNPs, and as many as ~69,000 SNPs that could be genotyped at ~20,000 gene loci using an Infinium II array-more SNPs than are needed to use genomic selection in tree breeding programs. Ultimately, these genomic resources will enhance Douglas-fir breeding and allow us to better understand landscape-scale patterns of genetic variation and potential responses to climate change.

The Mutualist Laccaria bicolor Expresses a Core Gene Regulon During the Colonization of Diverse Host Plants and a Variable Regulon to Counteract Host-Specific Defenses

The coordinated transcriptomic responses of both mutualistic ectomycorrhizal (ECM) fungi and their hosts during the establishment of symbiosis are not well-understood. This study characterizes the transcriptomic alterations of the ECM fungus Laccaria bicolor during different colonization stages on two hosts (Populus trichocarpa and Pseudotsuga menziesii) and compares this to the transcriptomic variations of P. trichocarpa across the same time-points. A large number of L. bicolor genes (≥ 8,000) were significantly regulated at the transcriptional level in at least one stage of colonization. From our data, we identify 1,249 genes that we hypothesize is the 'core' gene regulon necessary for the mutualistic interaction between L. bicolor and its host plants. We further identify a group of 1,210 genes that are regulated in a host-specific manner. This variable regulon encodes a number of genes coding for proteases and xenobiotic efflux transporters that we hypothesize act to counter chemical-based defenses simultaneously activated at the transcriptomic level in P. trichocarpa. The transcriptional response of the host plant P. trichocarpa consisted of differential waves of gene regulation related to signaling perception and transduction, defense response, and the induction of nutrient transfer in P. trichocarpa tissues. This study, therefore, gives fresh insight into the shifting transcriptomic landscape in both the colonizing fungus and its host and the different strategies employed by both partners in orchestrating a mutualistic interaction.

Climate-related genetic variation in drought-resistance of Douglas-fir (Pseudotsuga menziesii)

There is a general assumption that intraspecific populations originating from relatively arid climates will be better adapted to cope with the expected increase in drought from climate change. For ecologically and economically important species, more comprehensive, genecological studies that utilize large distributions of populations and direct measures of traits associated with drought-resistance are needed to empirically support this assumption because of the implications for the natural or assisted regeneration of species. We conducted a space-for-time substitution, common garden experiment with 35 populations of coast Douglas-fir (Pseudotsuga menziesii var. menziesii) growing at three test sites with distinct summer temperature and precipitation (referred to as 'cool/moist', 'moderate', or 'warm/dry') to test the hypotheses that (i) there is large genetic variation among populations and regions in traits associated with drought-resistance, (ii) the patterns of genetic variation are related to the native source-climate of each population, in particular with summer temperature and precipitation, (iii) the differences among populations and relationships with climate are stronger at the warm/dry test site owing to greater expression of drought-resistance traits (i.e., a genotype × environment interaction). During midsummer 2012, we measured the rate of water loss after stomatal closure (transpiration(min)), water deficit (% below turgid saturation), and specific leaf area (SLA, cm(2) g(-1)) on new growth of sapling branches. There was significant genetic variation in all plant traits, with populations originating from warmer and drier climates having greater drought-resistance (i.e., lower transpiration(min), water deficit and SLA), but these trends were most clearly expressed only at the warm/dry test site. Contrary to expectations, populations from cooler climates also had greater drought-resistance across all test sites. Multiple regression analysis indicated that Douglas-fir populations from regions with relatively cool winters and arid summers may be most adapted to cope with drought conditions that are expected in the future.

Genomic selection of juvenile height across a single-generational gap in Douglas-fir

Here, we perform cross-generational GS analysis on coastal Douglas-fir (Pseudotsuga menziesii), reflecting trans-generational selective breeding application. A total of 1321 trees, representing 37 full-sib F1 families from 3 environments in British Columbia, Canada, were used as the training population for (1) EBVs (estimated breeding values) of juvenile height (HTJ) in the F1 generation predicting genomic EBVs of HTJ of 136 individuals in the F2 generation, (2) deregressed EBVs of F1 HTJ predicting deregressed genomic EBVs of F2 HTJ, (3) F1 mature height (HT35) predicting HTJ EBVs in F2, and (4) deregressed F1 HT35 predicting genomic deregressed HTJ EBVs in F2. Ridge regression best linear unbiased predictor (RR-BLUP), generalized ridge regression (GRR), and Bayes-B GS methods were used and compared to pedigree-based (ABLUP) predictions. GS accuracies for scenarios 1 (0.92, 0.91, and 0.91) and 3 (0.57, 0.56, and 0.58) were similar to their ABLUP counterparts (0.92 and 0.60, respectively) (using RR-BLUP, GRR, and Bayes-B). Results using deregressed values fell dramatically for both scenarios 2 and 4 which approached zero in many cases. Cross-generational GS validation of juvenile height in Douglas-fir produced predictive accuracies almost as high as that of ABLUP. Without capturing LD, GS cannot surpass the prediction of ABLUP. Here we tracked pedigree relatedness between training and validation sets. More markers or improved distribution of markers are required to capture LD in Douglas-fir. This is essential for accurate forward selection among siblings as markers that track pedigree are of little use for forward selection of individuals within controlled pollinated families.

Characterization of four terpene synthase cDNAs from methyl jasmonate-induced Douglas-fir, Pseudotsuga menziesii

Numerous terpenoid compounds are present in copious amounts in the oleoresin produced by conifers, especially following exposure to insect or fungal pests. CDNA clones for many terpene synthases responsible for the biosynthesis of these defense compounds have been recovered from several conifer species. Here, the use of three terpene synthase sequences as heterologous probes for the discovery of related terpene synthase genes in Douglas-fir, Pseudotsuga menziesii (Mirbel) Franco (Pinaceae), is reported. Four full-length terpene synthase cDNAs were recovered from a methyl jasmonate-induced Douglas-fir bark and shoot cDNA library. These clones encode two multi-product monoterpene synthases [a (-)-alpha-pinene/(-)-camphene synthase and a terpinolene synthase] and two single-product sesquiterpene synthases [an (E)-beta-farnesene synthase and a (E)-gamma-bisabolene synthase].

Impact of climate change on cold hardiness of Douglas-fir (Pseudotsuga menziesii): environmental and genetic considerations

The success of conifers over much of the world's terrestrial surface is largely attributable to their tolerance to cold stress (i.e., cold hardiness). Due to an increase in climate variability, climate change may reduce conifer cold hardiness, which in turn could impact ecosystem functioning and productivity in conifer-dominated forests. The expression of cold hardiness is a product of environmental cues (E), genetic differentiation (G), and their interaction (G × E), although few studies have considered all components together. To better understand and manage for the impacts of climate change on conifer cold hardiness, we conducted a common garden experiment replicated in three test environments (cool, moderate, and warm) using 35 populations of coast Douglas-fir (Pseudotsuga menziesii var. menziesii) to test the hypotheses: (i) cool-temperature cues in fall are necessary to trigger cold hardening, (ii) there is large genetic variation among populations in cold hardiness that can be predicted from seed-source climate variables, (iii) observed differences among populations in cold hardiness in situ are dependent on effective environmental cues, and (iv) movement of seed sources from warmer to cooler climates will increase risk to cold injury. During fall 2012, we visually assessed cold damage of bud, needle, and stem tissues following artificial freeze tests. Cool-temperature cues (e.g., degree hours below 2 °C) at the test sites were associated with cold hardening, which were minimal at the moderate test site owing to mild fall temperatures. Populations differed 3-fold in cold hardiness, with winter minimum temperatures and fall frost dates as strong seed-source climate predictors of cold hardiness, and with summer temperatures and aridity as secondary predictors. Seed-source movement resulted in only modest increases in cold damage. Our findings indicate that increased fall temperatures delay cold hardening, warmer/drier summers confer a degree of cold hardiness, and seed-source movement from warmer to cooler climates may be a viable option for adapting coniferous forest to future climate.

Ethylene in induced conifer defense: cDNA cloning, protein expression, and cellular and subcellular localization of 1-aminocyclopropane-1-carboxylate oxidase in resin duct and phenolic parenchyma cells

Members of the Pinaceae family have complex chemical defense strategies. Conifer defenses associated with specialized cell types of the bark involve constitutive and inducible accumulation of phenolic compounds in polyphenolic phloem parenchyma cells and oleoresin terpenoids in resin ducts. These defenses can protect trees against insect herbivory and fungal colonization. The phytohormone ethylene has been shown to induce the same anatomical and cellular defense responses that occur following insect feeding, mechanical wounding, or fungal inoculation in Douglas fir (Pseudotsuga menziesii) stems (Hudgins and Franceschi in Plant Physiol 135:2134-2149, 2004). However, very little is known about the genes involved in ethylene formation in conifer defense or about the temporal and spatial patterns of their protein expression. The enzyme 1-aminocyclopropane-1-carboxylate oxidase (ACO) catalyzes the final step in ethylene biosynthesis. We cloned full-length and near full-length ACO cDNAs from three conifer species, Sitka spruce (Picea sitchensis), white spruce (P. glauca), and Douglas fir, each with high similarity to Arabidopsis thaliana ACO proteins. Using an Arabidopsis anti-ACO antibody we determined that ACO is constitutively expressed in Douglas fir stem tissues and is up-regulated by mechanical wounding, consistent with the wound-induced increase of ethylene levels. Immunolocalization showed cytosolic ACO is predominantly present in specialized cell types of the wound-induced bark, specifically in epithelial cells of terpenoid-producing cortical resin ducts, in polyphenolic phloem parenchyma cells, and in ray parenchyma cells.

Functional analysis of a Douglas-fir metallothionein-like gene promoter: transient assays in zygotic and somatic embryos and stable transformation in transgenic tobacco

Douglas-fir (Pseudotsuga menziesii [Mirb] Franco) metallothionein (PmMT) cDNA encodes a novel cysteine- and serine-rich MT, indicating a new subtype or prototype MT from which other plant MTs may have evolved. A genomic library of Douglas-fir was screened using MT cDNA probes, and genomic sequences that mediate tissue-specific, temporal as well as inducible expression of the embryo-specific MT-gene were analyzed. The promoter region of the PmMT genomic clone (gPmMT) contained a hexameric G-box, two putative ethylene-responsive elements and an inverted repeat of a motif similar to the core metal regulatory element. Interestingly, comparison of the upstream region of Douglas-fir gPm2S1 and gPmMTa genes revealed a conserved motif, CATTATTGA, not found in any known angiosperm gene promoter. Chimeric gene constructs containing a series of deletions in the gPmMTa promoter fused to the uidA reporter gene were assayed in Douglas-fir and transgenic tobacco (Nicotiana tabacum L.). Transient-expression assays in Douglas-fir megagametophyte and zygotic embryos indicated that the sequence -190 to +88 of gPmMTa was sufficient to drive the expression of the reporter gene and that the 225-bp fragment (-677 to -453) contained sequences necessary for high-level expression. In transgenic tobacco seedlings the beta-glucuronidase activity was localized in the vacuolar tissue and proliferating tissue of the auxiliary buds and stem elongation zone. The gPmMTa promoter was not active in the seeds of transgenic tobacco or in the roots of seedlings up to 3 weeks old. Detailed studies of transient expression and stable transformation provided important information on evolutionary conservation as well as novel features found in the conifer promoter. This is the first report of an MT-like gene promoter from conifers.

Metabolite profiling of Douglas-fir (Pseudotsuga menziesii) field trials reveals strong environmental and weak genetic variation

The primary objective of this study was to assess metabolomics for its capacity to discern biological variation among 10 full-sib families of a Douglas-fir tree breeding population, replicated on two sites. The differential accumulation of small metabolites in developing xylem was examined through metabolite profiles (139 metabolites common to 181 individual trees) generated by gas chromatography mass spectrometry and a series of statistical analyses that incorporated family, site, and tree growth and quantitative phenotypic wood traits (wood density, microfibril angle, wood chemistry and fiber morphology). Multivariate discriminant, canonical discriminant and factor analyses and broad-sense heritabilities revealed that metabolic and phenotypic traits alike were strongly related to site, while similar associations relating to genetic (family) structure were weak in comparison. Canonical correlation analysis subsequently identified correlations between specific phenotypic traits (i.e. tree growth, fibre morphology and wood chemistry) and metabolic traits (i.e. carbohydrate and lignin biosynthetic metabolites), demonstrating a coherent relationship between genetics, metabolism, environmental and phenotypic expression in wood-forming tissue. The association between cambial metabolites and tree phenotype, as revealed by metabolite profiling, demonstrates the value of metabolomics for systems biology approaches to understanding tree growth and secondary cell wall biosynthesis in plants.

Aminocyclopropane carboxylic acid synthase is a regulated step in ethylene-dependent induced conifer defense. Full-length cDNA cloning of a multigene family, differential constitutive, and wound- and insect-induced expression, and cellular and subcellular localization in spruce and Douglas fir

In conifer stems, formation of chemical defenses against insects or pathogens involves specialized anatomical structures of the phloem and xylem. Oleoresin terpenoids are formed in resin duct epithelial cells and phenolics accumulate in polyphenolic parenchyma cells. Ethylene signaling has been implicated in the induction of these chemical defenses. Recently, we reported the cloning of 1-aminocyclopropane-1-carboxylic acid oxidase (ACO) from spruce (Picea spp.) and Douglas fir (Pseudotsuga menziesii). ACO protein was constitutively expressed in Douglas fir and only weakly induced upon wounding. We now cloned seven full-length and one near full-length cDNA representing four distinct 1-aminocyclopropane-1-carboxylic acid synthases (ACS; ACS1, ACS2, ACS3, and ACS4) from spruce and Douglas fir. Cloning of ACS has not previously been reported for any gymnosperm. Using gene-specific, quantitative real-time polymerase chain reaction, we measured constitutive expression for the four ACS genes and the single-copy ACO gene in various tissues of Sitka spruce (Picea sitchensis) and in white spruce (Picea glauca) somatic embryos. ACO and ACS4 were ubiquitously expressed at high levels; ACS1 was predominantly expressed in developing embryos and ACS2 and ACS3 were expressed only at very low levels. Insect attack or mechanical wounding caused strong induction of ACS2 and ACS3 in Sitka spruce bark, a moderate increase in ACO transcripts, but had no effect on ACS1 and ACS4. ACS protein was also strongly induced following mechanical wounding in Douglas fir and was highly abundant in resin duct epithelial cells and polyphenolic parenchyma cells. These results suggest that ACS, but not ACO, is a regulated step in ethylene-induced conifer defense.

Growth phenology of coast Douglas-fir seed sources planted in diverse environments

The timing of periodic life cycle events in plants (phenology) is an important factor determining how species and populations will react to climate change. We evaluated annual patterns of basal-area and height growth of coast Douglas-fir (Pseudotusga menziesii var. menziesii (Mirb.) Franco) seedlings from four seed sources that were planted in four diverse environments as part of the Douglas-fir Seed-Source Movement Trial. Stem diameters and heights were measured periodically during the 2010 growing season on 16 open-pollinated families at each study installation. Stem diameters were measured on a subset of trees with electronic dendrometers during the 2010 and 2011 growing seasons. Trees from the four seed sources differed in phenology metrics that described the timing of basal-area and height-growth initiation, growth cessation and growth rates. Differences in the height-growth metrics were generally larger than differences in the basal-area growth metrics and differences among installations were larger than differences among seed sources, highlighting the importance of environmental signals on growth phenology. Variations in the height- and basal-area growth metrics were correlated with different aspects of the seed-source environments: precipitation in the case of height growth and minimum temperature in the case of basal-area growth. The detailed dendrometer measurements revealed differences in growth patterns between seed sources during distinct periods in the growing season. Our results indicate that multiple aspects of growth phenology should be considered along with other traits when evaluating adaptation of populations to future climates.

Dynamic changes in concentrations of auxin, cytokinin, ABA and selected metabolites in multiple genotypes of Douglas-fir (Pseudotsuga menziesii) during a growing season

Changes in concentrations of several endogenous phytohormones and metabolites were analyzed in the long shoots of nine genotypes of coastal Douglas-fir (Pseudotsuga menziesii (Mirb.) Franco var. menziesii) at five developmental stages: (1) closed buds, (2) flushing buds, (3) rapidly elongating shoots, (4) growing shoots and (5) near full-length shoots during one growing season. When averaged across genotypes, indole-3-acetic acid (IAA) concentration was high at stages 1 and 3. The only pattern that correlated with cone productivity was the one that was unique to IAA, in which high concentrations at stages 3 and 4 were found in all genotypes with high female cone productivity. Concentrations of isopentenyl adenosine (iPA) decreased and zeatin riboside (ZR) concentrations increased as the buds initiated and differentiated; ZR was 30 and 28 ng g(-1) dry weight (DW) at stages 1 and 4, respectively, before increasing to 166 ng g(-1) DW at stage 5. Isopentenyl adenosine peaked at 92 ng g(-1) DW at stage 2 and declined to low concentrations at stages 4 and 5. Zeatin-O-glucoside was 30 ng g(-1) DW at stage 1, declined at stages 2 and 3 and increased at stages 4 and 5. High abscisic acid (ABA) concentrations were positively correlated with rapid shoot elongation (stages 1 and 2), but as growth slowed and terminated, ABA concentrations decreased. Abscisic acid was 7 microg g(-1) DW at stage 1, increased to 13 microg g(-1) DW at stage 2 and then declined. The glucosyl ester (GE) of ABA decreased rapidly in early summer, and increased inversely with an increase in ABA. Between stages 1 and 2, ABA-GE decreased from 10 to 0.2 microg g(-1) DW and then increased. Of the ABA catabolites studied, 7'-hydroxy-ABA was about 2 microg g(-1) DW at stage 1, declined at stages 2 and 3 and increased at stages 4 and 5; phaseic acid concentrations were low at all stages, whereas dihydrophaseic acid was detected only at stages 4 and 5.

Clair JB, Crepeau MW, Salzberg SL, Langley CH, Allen B, Neale DB. Dissecting the Polygenic Basis of Cold Adaptation Using Genome-Wide Association of Traits and Environmental Data in Douglas-fir

Understanding the genomic and environmental basis of cold adaptation is key to understand how plants survive and adapt to different environmental conditions across their natural range. Univariate and multivariate genome-wide association (GWAS) and genotype-environment association (GEA) analyses were used to test associations among genome-wide SNPs obtained from whole-genome resequencing, measures of growth, phenology, emergence, cold hardiness, and range-wide environmental variation in coastal Douglas-fir (Pseudotsuga menziesii). Results suggest a complex genomic architecture of cold adaptation, in which traits are either highly polygenic or controlled by both large and small effect genes. Newly discovered associations for cold adaptation in Douglas-fir included 130 genes involved in many important biological functions such as primary and secondary metabolism, growth and reproductive development, transcription regulation, stress and signaling, and DNA processes. These genes were related to growth, phenology and cold hardiness and strongly depend on variation in environmental variables such degree days below 0c, precipitation, elevation and distance from the coast. This study is a step forward in our understanding of the complex interconnection between environment and genomics and their role in cold-associated trait variation in boreal tree species, providing a baseline for the species' predictions under climate change.

Photoperiod cues and patterns of genetic variation limit phenological responses to climate change in warm parts of species' range: Modeling diameter-growth cessation in coast Douglas-fir

The phenology of diameter-growth cessation in trees will likely play a key role in mediating species and ecosystem responses to climate change. A common expectation is that warming will delay cessation, but the environmental and genetic influences on this process are poorly understood. We modeled the effects of temperature, photoperiod, and seed-source climate on diameter-growth-cessation timing in coast Douglas-fir (an ecologically and economically vital tree) using high-frequency growth measurements across broad environmental gradients for a range of genotypes from different seed sources. Our model suggests that cool temperatures or short photoperiods can induce cessation in autumn. At cool locations (high latitude and elevation), cessation seems to be induced primarily by low temperatures in early autumn (under relatively long photoperiods), so warming will likely delay cessation and extend the growing season. But at warm locations (low latitude or elevation), cessation seems to be induced primarily by short photoperiods later in autumn, so warming will likely lead to only slight extensions of the growing season, reflecting photoperiod limitations on phenological shifts. Trees from seed sources experiencing frequent frosts in autumn or early winter tended to cease growth earlier in the autumn, potentially as an adaptation to avoid frost. Thus, gene flow into populations in warm locations with little frost will likely have limited potential to delay mean cessation dates because these populations already cease growth relatively late. In addition, data from an abnormal heat wave suggested that very high temperatures during long photoperiods in early summer might also induce cessation. Climate change could make these conditions more common in warm locations, leading to much earlier cessation. Thus, photoperiod cues, patterns of genetic variation, and summer heat waves could limit the capacity of coast Douglas-fir to extend its growing season in response to climate change in the warm parts of its range.

The Douglas-fir BiP promoter is functional in Arabidopsis and responds to wounding

A DNA sequence representing the promoter region of the Douglas-fir (Pseudotsuga menziesii [Mirb.] Franco) luminal binding protein PmBiP (PmBiPPro1) was isolated using inverse polymerase chain reaction (iPCR). Transient expression analysis of PmBiPPro1 fused to the beta-glucuronidase (GUS) reporter gene demonstrated that this promoter is functional in germinating Douglas-fir embryos. Transgenic Arabidopsis plants containing PmBiPPro1:GUS reporter gene constructs revealed strong staining associated with actively dividing/expanding cells and secretory tissues in developing seedlings. Wounding of cotyledons resulted in an increase in local staining associated with cells surrounding the wound site. Deletion analysis showed that elements necessary for basal-level expression reside within a -261 to +16 bp region, although upstream elements are necessary for higher-level expression in germinating Douglas-fir embryos, developing Arabidopsis seedlings and wounded cotyledons. Correlation of the observed expression pattern with the known function of BiP suggests that pathways controlling expression are highly conserved between angiosperms and gymnosperms.

A transcriptomic resource for Douglas-fir seed development and analysis of transcription during late megagametophyte development

Douglas-fir transcriptomics. Douglas-fir (Pseudotsuga menziesii (Mirb.) Franco) is economically important with extensive breeding programs and seed trade. However, the molecular genetics of its seed development are largely unknown. We developed a transcriptome resource covering key developmental stages of megagametophytes over time: prefertilization, fertilization, embryogenesis, and early, unfertilized abortion. RNA sequencing reads were assembled de novo into 105,505 predicted high-confidence transcripts derived from 34,521 predicted genes. Expression levels were estimated based on alignment of the original reads to the reference. Megagametophytes express a distinct set of genes compared to those of vegetative tissues. Transcripts related to signaling, protein turnover, and RNA biogenesis have lower expression values in vegetative tissues, whereas cell wall remodeling, solute transport, and seed storage protein transcripts have higher expression values in megagametophytes. Seed storage protein transcripts become very abundant in both pollinated and unpollinated megagametophytes over time, even in aborting ovules. However, the absence of protein storage bodies in unfertilized megagametophytes suggests extensive posttranscriptional mechanisms that either inhibit storage protein translation or their aggregation into protein bodies. This novel transcriptome resource provides a foundation for further important insights into conifer seed development.

Spatiotemporal activities of Douglas-fir BiP Pro1 promoter in transgenic potato

The PmBiPPro1 promoter of the luminal binding protein (BiP) from Douglas-fir is fully functional in transgenic potato, responsive to wounding, and has high transcriptional activity in tubers. A predefined pattern and level of transgene expression targeted to specific tissues or organs and at a particular developmental stage is a pre-requisite for the successful development of plants with desired traits. Here, we evaluated the transcriptional activity of the PmBiPPro1 promoter of the luminal binding protein (BiP) from Douglas-fir, by expressing reporter β-D-glucuronidase (GUS) gene constructs containing three different PmBiPPro1 promoter versions (2258 bp, 1259 bp, and 278 bp) in transgenic potato. In conifers, this promoter regulates the endoplasmic reticulum (ER) molecular chaperon of the HSP70 stress-related protein family and is essential for proper functioning of the ER. Stable expression analysis demonstrated that two of three PmBiPPro1 promoter versions (PmBiPPro1-1 and PmBiPPro1-3) were fully functional in the heterologous host, exhibited high transcriptional activities in the leaves of unstressed potatoes, and were responsive to wounding. Deletion analysis showed that the positive cis-active regulatory elements necessary for higher level expression resided within the - 1243 to - 261 region, whereas negative cis-active elements encompassed nucleotides - 2242 to - 1243. Histochemical staining revealed high level of GUS activities in tissues associated with a high rate of cell division and secretory activities. Most importantly, the PmBiPPro1 promoters, especially the full-length version, had activity in microtubers at a level that was much higher than in any other potato organ or tissue. The - 2242 to - 1243 bp region likely contains important cis element(s) that interact with tuber-specific transcription factors required for promoter activation in the storage organs. The organ-specific activity of the PmBiPPro1 promoters may be useful for targeted expression of heterologous molecules in potato tubers.

Effects of inbreeding on coastal Douglas fir growth and yield in operational plantations: a model-based approach

In advanced generation seed orchards, tradeoffs exist between genetic gain obtained by selecting the best related individuals for seed orchard populations, and potential losses due to subsequent inbreeding between these individuals. Although inbreeding depression for growth rate is strong in most forest tree species at the individual tree level, the effect of a small proportion of inbreds in seed lots on final stand yield may be less important. The effects of inbreeding on wood production of mature stands cannot be assessed empirically in the short term, thus such effects were simulated for coastal Douglas fir [ Pseudotsuga menziesii var. menziesii (Mirb.) Franco] using an individual-tree growth and yield model TASS (Tree and Stand Simulator). The simulations were based on seed set, nursery culling rates, and 10-year-old field test performance for trees resulting from crosses between unrelated individuals and for inbred trees produced through mating between half-sibs, full-sibs, parents and offspring and self-pollination. Results indicate that inclusion of a small proportion of related clones in seed orchards will have relatively low impacts on stand yields due to low probability of related individuals mating, lower probability of producing acceptable seedlings from related matings than from unrelated matings, and a greater probability of competition-induced mortality for slower growing inbred individuals than for outcrossed trees. Thus, competition reduces the losses expected due to inbreeding depression at harvest, particularly on better sites with higher planting densities and longer rotations. Slightly higher breeding values for related clones than unrelated clones would offset or exceed the effects of inbreeding resulting from related matings. Concerns regarding the maintenance of genetic diversity are more likely to limit inclusion of related clones in orchards than inbreeding depression for final stand yield.

Evaluation of resistance to Rhabdocline needlecast in Douglas fir variety Shuswap, with quantitative polymerase chain reaction

A quantitative polymerase chain reaction assay was developed that could detect DNA of Rhabdocline pseudotsugae and R. oblonga among DNA of Douglas fir needles to a limit as low as three copies of target DNA. Differential infection rates of two varieties (seed sources) of Douglas fir interplanted in a field were studied in relation to staggered bud breaks. Infection of Douglas fir var. San Isabel corresponded to ascospore release times for Rhabdocline spp., whereas infection of var. Shuswap Lake did not occur throughout the spore release period during 2 years of study, despite abundant inoculum and adequate moisture during bud break. Rhabdocline spp. DNA was never detected in Shuswap Lake and disease symptoms were not observed in any year. We provide evidence that Shuswap Lake is resistant and probably immune to Rhabdocline spp. infection and Rhabdocline needlecast under Michigan conditions.

Spatially heterogeneous selection and inter-varietal differentiation maintain population structure and local adaptation in a widespread conifer

BACKGROUND

Douglas-fir (Pseudotsuga menziesii [Mirb.] Franco) plays a critical role in the ecology and economy of Western North America. This conifer species comprises two distinct varieties: the coastal variety (var. menziesii) along the Pacific coast, and the interior variety (var. glauca) spanning the Rocky Mountains into Mexico, with instances of inter-varietal hybridization in Washington and British Columbia. Recent investigations have focused on assessing environmental pressures shaping Douglas-fir's genomic variation for a better understanding of its evolutionary and adaptive responses. Here, we characterize range-wide population structure, estimate inter-varietal hybridization levels, identify candidate loci for climate adaptation, and forecast shifts in species and variety distribution under future climates.

RESULTS

Using a custom SNP-array, we genotyped 540 trees revealing four distinct clusters with asymmetric admixture patterns in the hybridization zone. Higher genetic diversity observed in coastal and hybrid populations contrasts with lower diversity in inland populations of the southern Rockies and Mexico, exhibiting a significant isolation by distance pattern, with less marked but still significant isolation by environment. For both varieties, we identified candidate loci associated with local adaptation, with hundreds of genes linked to processes such as stimulus response, reactions to chemical compounds, and metabolic functions. Ecological niche modeling revealed contrasting potential distribution shifts among the varieties in the coming decades, with interior populations projected to lose habitat and become more vulnerable, while coastal populations are expected to gain suitable areas.

CONCLUSIONS

Overall, our findings provide crucial insights into the population structure and adaptive potential of Douglas-fir, with the coastal variety being the most likely to preserve its evolutionary path throughout the present century, which carry implications for the conservation and management of this species across their range.