Correlating gene expression to physiological parameters and environmental conditions during cold acclimation of Pinus sylvestris, identification of molecular markers using cDNA microarrays

Comparative transcriptome analyses reveal two distinct transcriptional modules associated with pollen shedding time in pine

Background

Seasonal flowering time is an ecologically and economically important trait in temperate trees. Previous studies have shown that temperature in many tree species plays a pivotal role in regulating flowering time. However, genetic control of flowering time is not synchronised in different individual trees under comparable temperature conditions, the underlying molecular mechanism is mainly to be investigated.

Results

In the present study, we analysed the transcript abundance in male cones and needles from six early pollen-shedding trees (EPs) and six neighbouring late pollen-shedding trees (LPs) in Pinus tabuliformis at three consecutive time points in early spring. We found that the EPs and LPs had distinct preferred transcriptional modules in their male cones and, interestingly, the expression pattern was also consistently maintained in needles even during the winter dormancy period. Additionally, the preferred pattern in EPs was also adopted by other fast-growing tissues, such as elongating new shoots. Enhancement of nucleic acid synthesis and stress resistance pathways under cold conditions can facilitate rapid growth and maintain higher transcriptional activity.

Conclusions

During the cold winter and early spring seasons, the EPs were more sensitive to relatively warmer temperatures and showed higher transcriptomic activity than the LPs, indicating that EPs required less heat accumulation for pollen shedding than LPs. These results provided a transcriptomic-wide understanding of the temporal regulation of pollen shedding in pines.

Combined transcriptomic and metabolomic analyses uncover rearranged gene expression and metabolite metabolism in tobacco during cold acclimation

Cold temperatures often severely restrict the growth, distribution and productivity of plants. The freezing tolerance of plants from temperate climates can be improved by undergoing periods of cold acclimation (CA). Tobacco is an important economic plant and is sensitive to cold stress. However, the dynamic changes and regulatory mechanisms of gene expression and metabolic processes during CA remain largely unknown. In this study, we performed RNA sequencing and metabolomic profiling analyses to identify the genes and metabolites specifically expressed during CA. Our transcriptomic data revealed 6905 differentially expressed genes (DEGs) during CA. Functional annotation and enrichment analyses revealed that the DEGs were involved mainly in signal transduction, carbohydrate metabolism and phenylpropanoid biosynthesis. Moreover, a total of 35 significantly changed metabolites were identified during CA via an LC-MS platform. Many protective metabolites, such as amino acids, carbohydrates, tricarboxylic acid (TCA) cycle intermediates and phenylpropanoid-related substances, were identified during CA. The gene-metabolite network extensively outlined the biological processes associated with the utilization of sugars, activation of amino acid metabolism, TCA cycle and phenylpropanoid biosynthesis in tobacco under CA. The results of our present study provide a comprehensive view of signal transduction and regulation, gene expression and dynamic changes in metabolites during CA.

Morpho-Physiological and Genomic Evaluation of Juglans Species Reveals Regional Maladaptation to Cold Stress

Climate change may have unpredictable effects on the cold hardiness of woody species planted outside of their range of origin. Extreme undulations in temperatures may exacerbate susceptibility to cold stress, thereby interfering with productivity and ecosystem functioning. Juglans L. and their naturally occurring interspecific F1 hybrids, are distributed natively across many temperate regions, and J. regia has been extensively introduced. Cold hardiness, an environmental and genetic factor yet to be evaluated in many native and introduced Juglans species, may be a limiting factor under future climate change and following species introductions. We evaluated cold hardiness of native North American and Eastern Asian Juglans along with J. regia genotypes using field data from the Midwestern United States (Indiana), controlled freezing tests, and genome sequencing with close assessment of Juglans cold hardy genes. Many Juglans species previously screened for cold-hardiness were genotypes derived from the Midwest, California, and Europe. In 2014, despite general climate adaptation, Midwestern winter temperatures of -30°C killed J. regia originating from California; however, naturalized Midwestern J. regia survived and displayed low damage. Hybridization of J. regia with black walnut (J. nigra) and butternut (J. cinerea) produced F1s displaying greater cold tolerance than pure J. regia. Cold hardiness and growth are variable in Midwestern J. regia compared to native Juglans, East Asian Juglans, and F1 hybrids. Phylogeny analyses revealed that J. cinerea sorted with East Asian species using the nuclear genome but with North American species using the organellar genome. Investigation of selected cold hardy genes revealed that J. regia was distinct from other species and exhibited less genetic diversity than native Juglans species Average whole genome heterozygosity and Tajima's D for cold hardy genes was low within J. regia samples and significantly higher for hybrid as well as J. nigra. We confirmed that molecular and morpho-physiological data were highly correlated and thus can be used effectively to characterize cold hardiness in Juglans species. We conclude that the genetic diversity within local J. regia populations is low and additional germplasm is needed for development of more regionally adapted J. regia varieties.

Expansion of the dehydrin gene family in the Pinaceae is associated with considerable structural diversity and drought-responsive expression

Temperatures are expected to increase over the next century in all terrestrial biomes and particularly in boreal forests, where drought-induced mortality has been predicted to rise. Genomics research is helping to develop hypotheses regarding the molecular basis of drought tolerance and recent work proposed that the osmo-protecting dehydrin proteins have undergone a clade-specific expansion in the Pinaceae, a major group of conifer trees. The objectives of this study were to identify all of the putative members of the gene family, trace their evolutionary origin, examine their structural diversity and test for drought-responsive expression. We identified 41 complete dehydrin coding sequences in Picea glauca, which is four times more than most angiosperms studied to date, and more than in pines. Phylogenetic reconstructions indicated that the family has undergone an expansion in conifers, with parallel evolution implicating the sporadic resurgence of certain amino acid sequence motifs, and a major duplication giving rise to a clade specific to the Pinaceae. A variety of plant dehydrin structures were identified with variable numbers of the A-, E-, S- and K-segments and an N-terminal (N1) amino acid motif including assemblages specific to conifers. The expression of several of the spruce dehydrins was tissue preferential under non-stressful conditions or responded to water stress after 7-18 days without watering, reflecting changes in osmotic potential. We found that dehydrins with N1 K2 and N1 AESK2 sequences were the most responsive to the lack of water. Together, the family expansion, drought-responsive expression and structural diversification involving loss and gain of amino acid motifs suggests that subfunctionalization has driven the diversification seen among dehydrin gene duplicates. Our findings clearly indicate that dehydrins represent a large family of candidate genes for drought tolerance in spruces and in other Pinaceae that may underpin adaptability in spatially and temporally variable environments.

The cold responsive mechanism of the paper mulberry: decreased photosynthesis capacity and increased starch accumulation

Background

Most studies on the paper mulberry are mainly focused on the medicated and pharmacology, fiber quality, leaves feed development, little is known about its mechanism of adaptability to abiotic stress. Physiological measurement, transcriptomics and proteomic analysis were employed to understand its response to cold stress in this study.

Methods

The second to fourth fully expanded leaves from up to down were harvested at different stress time points forthe transmission electron microscope (TEM) observation. Physiological characteristics measurement included the relative electrolyte leakage (REL), SOD activity assay, soluble sugar content, and Chlorophyll fluorescence parameter measurement. For screening of differentially expressed genes, the expression level of every transcript in each sample was calculated by quantifying the number of Illumina reads. To identify the differentially expressed protein, leaves of plants under 0, 6, 12, 24, 48 and 72 h cold stress wereharvested for proteomic analysis. Finally, real time PCR was used to verify the DEG results of the RNA-seq and the proteomics data.

Results

Results showed that at the beginning of cold stress, respiratory metabolism was decreased and the transportation and hydrolysis of photosynthetic products was inhibited, leading to an accumulation of starch in the chloroplasts. Total of 5800 unigenes and 38 proteins were affected, including the repressed expression of photosynthesis and the enhanced expression in signal transduction, stress defense pathway as well as secondary metabolism. Although the transcriptional level of a large number of genes has been restored after 12 h, sustained cold stress brought more serious injury to the leaf cells, including the sharp rise of the relative electrolyte leakage, the declined Fv/Fm value, swelled chloroplast and the disintegrated membrane system.

Conclusion

The starch accumulation and the photoinhibition might be the main adaptive mechanism of the paper mulberry responded to cold stress. Most of important, enhancing the transport and hydrolysis of photosynthetic products could be the potential targets for improving the cold tolerance of the paper mulberry.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-2047-6) contains supplementary material, which is available to authorized users.

Nucleotide polymorphisms in a pine ortholog of the Arabidopsis degrading enzyme cellulase KORRIGAN are associated with early growth performance in Pinus pinaster

We have carried out a candidate-gene-based association genetic study in Pinus pinaster Aiton and evaluated the predictive performance for genetic merit gain of the most significantly associated genes and single nucleotide polymorphisms (SNPs). We used a second generation 384-SNP array enriched with candidate genes for growth and wood properties to genotype mother trees collected in 20 natural populations covering most of the European distribution of the species. Phenotypic data for total height, polycyclism, root-collar diameter and biomass were obtained from a replicated provenance-progeny trial located in two sites with contrasting environments (Atlantic vs Mediterranean climate). General linear models identified strong associations between growth traits (total height and polycyclism) and four SNPs from the korrigan candidate gene, after multiple testing corrections using false discovery rate. The combined genomic breeding value predictions assessed for the four associated korrigan SNPs by ridge regression-best linear unbiased prediction (RR-BLUP) and cross-validation accounted for up to 8 and 15% of the phenotypic variance for height and polycyclic growth, respectively, and did not improve adding SNPs from other growth-related candidate genes. For root-collar diameter and total biomass, they accounted for 1.6 and 1.1% of the phenotypic variance, respectively, but increased to 15 and 4.1% when other SNPs from lp3.1, lp3.3 and cad were included in RR-BLUP models. These results point towards a desirable integration of candidate-gene studies as a means to pre-select relevant markers, and aid genomic selection in maritime pine breeding programs.

Comparative transcriptomics of a complex of four European pine species

Background

Pinus sylvestris, P. mugo, P. uliginosa and P. uncinata are closely related but phenotypically and ecologically very distinct European pine species providing an excellent study system for analysis of the genetic basis of adaptive variation and speciation. For comparative genomic analysis of the species, transcriptome sequence was generated for 17 samples collected across the European distribution range using Illumina paired-end sequencing technology.

Results

De novo transcriptome assembly of a reference sample of P. sylvestris contained 40968 unigenes, of which fewer than 0.5% were identified as putative retrotransposon sequences. Based on gene annotation approaches, 19659 contigs were identified and assigned to unique genes covering a broad range of gene ontology categories. About 80% of the reads from each sample were successfully mapped to the reference transcriptome of P. sylvestris. Single nucleotide polymorphisms were identified in 22041-24096 of the unigenes providing a set of ~220-262 k SNPs identified for each species. Very similar levels of nucleotide polymorphism were observed across species (π=0.0044-0.0053) and highest pairwise nucleotide divergence (0.006) was found between P. mugo and P. sylvestris at a common set of unigenes.

Conclusions

The study provides whole transcriptome sequence and a large set of SNPs to advance population and association genetic studies in pines. Our study demonstrates that transcriptome sequencing can be a very useful approach for development of novel genomic resources in species with large and complex genomes.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1401-z) contains supplementary material, which is available to authorized users.

Clock genes and diurnal transcriptome dynamics in summer and winter in the gymnosperm Japanese cedar (Cryptomeria japonica (L.f.) D.Don)

BACKGROUND

The circadian clock and diurnal dynamics of the transcriptome are presumed to play important roles in the regulation of physiological, biological and developmental processes synchronized with diurnal and annual cycles of plant environments. However, little is known about the circadian clock and its regulation in gymnosperms, including conifers. Here we present the diurnal transcriptome dynamics of Japanese cedar (Cryptomeria japonica (L.f.) D.Don) in both active (summer) and dormant (winter) periods.

RESULTS

Microarray analysis revealed significant differences in transcripts between summer and winter, and diurnal transcriptome dynamics only in the summer. About 7.7% of unique genes (556 out of 7,254) on the microarray were periodically expressed in summer. Expression patterns of some genes, especially light-related genes, did not show significant oscillation in Japanese cedar, thus differing from those reported in angiosperms. Gene network analysis of the microarray data revealed a network associated with the putative core clock genes (CjLHYa, CjLHYb, CjTOC1, CjGI and CjZTL), which were also isolated, indicating their importance in the diurnal regulation of the transcriptome.

CONCLUSION

This study revealed the existence of core clock genes and diurnal rhythms of the transcriptome in summer in Japanese cedar. Dampening of diurnal rhythms in winter indicated seasonal change in the rhythms according to environmental conditions. The data also revealed genes that showed different expression patterns compared to angiosperms, suggesting a unique gene regulatory network in conifers. This study provides fundamental data to understand transcriptional regulatory mechanisms in conifers.

Molecular and physiological responses to abiotic stress in forest trees and their relevance to tree improvement

Abiotic stresses, such as drought, salinity and cold, are the major environmental stresses that adversely affect tree growth and, thus, forest productivity, and play a major role in determining the geographic distribution of tree species. Tree responses and tolerance to abiotic stress are complex biological processes that are best analyzed at a systems level using genetic, genomic, metabolomic and phenomic approaches. This will expedite the dissection of stress-sensing and signaling networks to further support efficient genetic improvement programs. Enormous genetic diversity for stress tolerance exists within some forest-tree species, and due to advances in sequencing technologies the molecular genetic basis for this diversity has been rapidly unfolding in recent years. In addition, the use of emerging phenotyping technologies extends the suite of traits that can be measured and will provide us with a better understanding of stress tolerance. The elucidation of abiotic stress-tolerance mechanisms will allow for effective pyramiding of multiple tolerances in a single tree through genetic engineering. Here we review recent progress in the dissection of the molecular basis of abiotic stress tolerance in forest trees, with special emphasis on Populus, Pinus, Picea, Eucalyptus and Quercus spp. We also outline practices that will enable the deployment of trees engineered for abiotic stress tolerance to land owners. Finally, recommendations for future work are discussed.

Climatic adaptation and ecological divergence between two closely related pine species in Southeast China

Climate is one of the most important drivers for adaptive evolution in forest trees. Climatic selection contributes greatly to local adaptation and intraspecific differentiation, but this kind of selection could also have promoted interspecific divergence through ecological speciation. To test this hypothesis, we examined intra- and interspecific genetic variation at 25 climate-related candidate genes and 12 reference loci in two closely related pine species, Pinus massoniana Lamb. and Pinus hwangshanensis Hisa, using population genetic and landscape genetic approaches. These two species occur in Southeast China but have contrasting ecological preferences in terms of several environmental variables, notably altitude, although hybrids form where their distributions overlap. One or more robust tests detected signals of recent and/or ancient selection at two-thirds (17) of the 25 candidate genes, at varying evolutionary timescales, but only three of the 12 reference loci. The signals of recent selection were species specific, but signals of ancient selection were mostly shared by the two species likely because of the shared evolutionary history. FST outlier analysis identified six SNPs in five climate-related candidate genes under divergent selection between the two species. In addition, a total of 24 candidate SNPs representing nine candidate genes showed significant correlation with altitudinal divergence in the two species based on the covariance matrix of population history derived from reference SNPs. Genetic differentiation between these two species was higher at the candidate genes than at the reference loci. Moreover, analysis using the isolation-with-migration model indicated that gene flow between the species has been more restricted for climate-related candidate genes than the reference loci, in both directions. Taken together, our results suggest that species-specific and divergent climatic selection at the candidate genes might have counteracted interspecific gene flow and played a key role in the ecological divergence of these two closely related pine species.

Conservation and divergence of gene expression plasticity following c. 140 million years of evolution in lodgepole pine (Pinus contorta) and interior spruce (Picea glauca×Picea engelmannii)

Species respond to environmental stress through a combination of genetic adaptation and phenotypic plasticity, both of which may be important for survival in the face of climatic change. By characterizing the molecular basis of plastic responses and comparing patterns among species, it is possible to identify how such traits evolve. Here, we used de novo transcriptome assembly and RNAseq to explore how patterns of gene expression differ in response to temperature, moisture, and light regime treatments in lodgepole pine (Pinus contorta) and interior spruce (a natural hybrid population of Picea glauca and Picea engelmannii). We found wide evidence for an effect of treatment on expression within each species, with 6413 and 11,658 differentially expressed genes identified in spruce and pine, respectively. Comparing patterns of expression among these species, we found that 74% of all orthologs with differential expression had a pattern that was conserved in both species, despite 140 million yr of evolution. We also found that the specific treatments driving expression patterns differed between genes with conserved versus diverged patterns of expression. We conclude that natural selection has probably played a role in shaping plastic responses to environment in these species.

Decoding the massive genome of loblolly pine using haploid DNA and novel assembly strategies

Background

The size and complexity of conifer genomes has, until now, prevented full genome sequencing and assembly. The large research community and economic importance of loblolly pine, Pinus taeda L., made it an early candidate for reference sequence determination.

Results

We develop a novel strategy to sequence the genome of loblolly pine that combines unique aspects of pine reproductive biology and genome assembly methodology. We use a whole genome shotgun approach relying primarily on next generation sequence generated from a single haploid seed megagametophyte from a loblolly pine tree, 20-1010, that has been used in industrial forest tree breeding. The resulting sequence and assembly was used to generate a draft genome spanning 23.2 Gbp and containing 20.1 Gbp with an N50 scaffold size of 66.9 kbp, making it a significant improvement over available conifer genomes. The long scaffold lengths allow the annotation of 50,172 gene models with intron lengths averaging over 2.7 kbp and sometimes exceeding 100 kbp in length. Analysis of orthologous gene sets identifies gene families that may be unique to conifers. We further characterize and expand the existing repeat library based on the de novo analysis of the repetitive content, estimated to encompass 82% of the genome.

Conclusions

In addition to its value as a resource for researchers and breeders, the loblolly pine genome sequence and assembly reported here demonstrates a novel approach to sequencing the large and complex genomes of this important group of plants that can now be widely applied.

Evidence for extensive heterotrophic metabolism, antioxidant action, and associated regulatory events during winter hardening in Sitka spruce

BACKGROUND

Cold acclimation in woody perennials is a metabolically intensive process, but coincides with environmental conditions that are not conducive to the generation of energy through photosynthesis. While the negative effects of low temperatures on the photosynthetic apparatus during winter have been well studied, less is known about how this is reflected at the level of gene and metabolite expression, nor how the plant generates primary metabolites needed for adaptive processes during autumn.

RESULTS

The MapMan tool revealed enrichment of the expression of genes related to mitochondrial function, antioxidant and associated regulatory activity, while changes in metabolite levels over the time course were consistent with the gene expression patterns observed. Genes related to thylakoid function were down-regulated as expected, with the exception of plastid targeted specific antioxidant gene products such as thylakoid-bound ascorbate peroxidase, components of the reactive oxygen species scavenging cycle, and the plastid terminal oxidase. In contrast, the conventional and alternative mitochondrial electron transport chains, the tricarboxylic acid cycle, and redox-associated proteins providing reactive oxygen species scavenging generated by electron transport chains functioning at low temperatures were all active.

CONCLUSIONS

A regulatory mechanism linking thylakoid-bound ascorbate peroxidase action with "chloroplast dormancy" is proposed. Most importantly, the energy and substrates required for the substantial metabolic remodeling that is a hallmark of freezing acclimation could be provided by heterotrophic metabolism.

Molecular physiology and breeding at the crossroads of cold hardiness improvement

Alfalfa (Medicago sativa L.) is a major forage legume grown extensively worldwide with important agronomic and environmental attributes. Insufficient cold hardiness is a major impediment to its reliable production in northern climates. Improvement of freezing tolerance using conventional breeding approaches is slowed by the quantitative nature of inheritance and strong interactions with the environment. The development of gene-based markers would facilitate the identification of genotypes with superior stress tolerance. Successive cycles of recurrent selection were applied using an indoor screening method to develop populations with significantly higher tolerance to freezing (TF). Bulk segregant analysis of heterogeneous TF populations identified DNA variations that are progressively enriched in frequency in response to selection. Polymorphisms resulting from intragenic variations within a dehydrin gene were identified and could potentially lead to the development of robust selection tools. Our results illustrate the benefits of feedback interactions between germplasm development programs and molecular physiology for a deeper understanding of the molecular and genetic bases of cold hardiness.

Novel conserved segments are associated with differential expression patterns for Pinaceae dehydrins

Dehydrins are thought to play an essential role in the response, acclimation and tolerance to different abiotic stresses, such as cold and drought. These proteins have been classified into five groups according to the presence of conserved and repeated motifs in their amino acid sequence. Due to their putative functions in the response to stress, dehydrins have been often used as candidate genes in studies on population variability and local adaptation to environmental conditions. However, little is still known regarding the differential role played by such groups or the mechanism underlying their function. Based on the sequences corresponding to dehydrins available in public databases we have isolated eight different dehydrins from cDNA of Pinus pinaster. We have obtained also their genomic sequences and identified their intron/exon structure. Quantitative RT-PCR analysis of their expression pattern in needles, stems and roots during a severe and prolonged drought stress, similar to the ones trees must face in nature, is also reported. Additionally, we have identified two amino acid motifs highly conserved and repeated in Pinaceae dehydrins and absent in angiosperms, presumably related to the divergent expression profiles observed.

The dynamic nature of bud dormancy in trees: environmental control and molecular mechanisms

In tree species native to temperate and boreal regions, the activity-dormancy cycle is an important adaptive trait both for survival and growth. We discuss recent research on mechanisms controlling the overlapping developmental processes that define the activity-dormancy cycle, including cessation of apical growth, bud development, induction, maintenance and release of dormancy, and bud burst. The cycle involves an extensive reconfiguration of metabolism. Environmental control of the activity-dormancy cycle is based on perception of photoperiodic and temperature signals, reflecting adaptation to prevailing climatic conditions. Several molecular actors for control of growth cessation have been identified, with the CO/FT regulatory network and circadian clock having important coordinating roles in control of growth and dormancy. Other candidate regulators of bud set, dormancy and bud burst have been identified, such as dormancy-associated MADS-box factors, but their exact roles remain to be discovered. Epigenetic mechanisms also appear to factor in control of the activity-dormancy cycle. Despite evidence for gibberellins as negative regulators in growth cessation, and ABA and ethylene in bud formation, understanding of the roles that plant growth regulators play in controlling the activity-dormancy cycle is still very fragmentary. Finally, some of the challenges for further research in bud dormancy are discussed.

Mining and visualization of microarray and metabolomic data reveal extensive cell wall remodeling during winter hardening in Sitka spruce (Picea sitchensis)

Microarray gene expression profiling is a powerful technique to understand complex developmental processes, but making biologically meaningful inferences from such studies has always been challenging. We previously reported a microarray study of the freezing acclimation period in Sitka spruce (Picea sitchensis) in which a large number of candidate genes for climatic adaptation were identified. In the current paper, we apply additional systems biology tools to these data to further probe changes in the levels of genes and metabolites and activities of associated pathways that regulate this complex developmental transition. One aspect of this adaptive process that is not well understood is the role of the cell wall. Our data suggest coordinated metabolic and signaling responses leading to cell wall remodeling. Co-expression of genes encoding proteins associated with biosynthesis of structural and non-structural cell wall carbohydrates was observed, which may be regulated by ethylene signaling components. At the same time, numerous genes, whose products are putatively localized to the endomembrane system and involved in both the synthesis and trafficking of cell wall carbohydrates, were up-regulated. Taken together, these results suggest a link between ethylene signaling and biosynthesis, and targeting of cell wall related gene products during the period of winter hardening. Automated Layout Pipeline for Inferred NEtworks (ALPINE), an in-house plugin for the Cytoscape visualization environment that utilizes the existing GeneMANIA and Mosaic plugins, together with the use of visualization tools, provided images of proposed signaling processes that became active over the time course of winter hardening, particularly at later time points in the process. The resulting visualizations have the potential to reveal novel, hypothesis-generating, gene association patterns in the context of targeted subcellular location.

High genetic diversity at the extreme range edge: nucleotide variation at nuclear loci in Scots pine (Pinus sylvestris L.) in Scotland

Nucleotide polymorphism at 12 nuclear loci was studied in Scots pine populations across an environmental gradient in Scotland, to evaluate the impacts of demographic history and selection on genetic diversity. At eight loci, diversity patterns were compared between Scottish and continental European populations. At these loci, a similar level of diversity (θ(sil)= ~0.01) was found in Scottish vs mainland European populations, contrary to expectations for recent colonization, however, less rapid decay of linkage disequilibrium was observed in the former (ρ=0.0086±0.0009, ρ=0.0245±0.0022, respectively). Scottish populations also showed a deficit of rare nucleotide variants (multi-locus Tajima's D=0.316 vs D=-0.379) and differed significantly from mainland populations in allelic frequency and/or haplotype structure at several loci. Within Scotland, western populations showed slightly reduced nucleotide diversity (π(tot)=0.0068) compared with those from the south and east (0.0079 and 0.0083, respectively) and about three times higher recombination to diversity ratio (ρ/θ=0.71 vs 0.15 and 0.18, respectively). By comparison with results from coalescent simulations, the observed allelic frequency spectrum in the western populations was compatible with a relatively recent bottleneck (0.00175 × 4N(e) generations) that reduced the population to about 2% of the present size. However, heterogeneity in the allelic frequency distribution among geographical regions in Scotland suggests that subsequent admixture of populations with different demographic histories may also have played a role.

Molecular events of apical bud formation in white spruce, Picea glauca

Bud formation is an adaptive trait that temperate forest trees have acquired to facilitate seasonal synchronization. We have characterized transcriptome-level changes that occur during bud formation of white spruce [Picea glauca (Moench) Voss], a primarily determinate species in which preformed stem units contained within the apical bud constitute most of next season's growth. Microarray analysis identified 4460 differentially expressed sequences in shoot tips during short day-induced bud formation. Cluster analysis revealed distinct temporal patterns of expression, and functional classification of genes in these clusters implied molecular processes that coincide with anatomical changes occurring in the developing bud. Comparing expression profiles in developing buds under long day and short day conditions identified possible photoperiod-responsive genes that may not be essential for bud development. Several genes putatively associated with hormone signalling were identified, and hormone quantification revealed distinct profiles for abscisic acid (ABA), cytokinins, auxin and their metabolites that can be related to morphological changes to the bud. Comparison of gene expression profiles during bud formation in different tissues revealed 108 genes that are differentially expressed only in developing buds and show greater transcript abundance in developing buds than other tissues. These findings provide a temporal roadmap of bud formation in white spruce.

Forest tree genomics: growing resources and applications

Over the past two decades, research in forest tree genomics has lagged behind that of model and agricultural systems. However, genomic research in forest trees is poised to enter into an important and productive phase owing to the advent of next-generation sequencing technologies, the enormous genetic diversity in forest trees and the need to mitigate the effects of climate change. Research on long-lived woody perennials is extending our molecular knowledge of complex life histories and adaptations to the environment - enriching a field that has traditionally drawn biological inference from a few short-lived herbaceous species.

The environmental contribution to gene expression profiles

Microarray analysis provides a bridge between the molecular genetic analysis of model organisms in laboratory settings and studies of physiology, development, and adaptation in the wild. By sampling species across a range of environments, it is possible to gain a broad picture of the genomic response to environmental perturbation. Incorporating estimates of genetic relationships into study designs will facilitate genomic analysis of environmental plasticity by aiding the identification of major regulatory loci in natural populations.

Global monitoring of autumn gene expression within and among phenotypically divergent populations of Sitka spruce (Picea sitchensis)

Cold acclimation in conifers is a complex process, the timing and extent of which reflects local adaptation and varies widely along latitudinal gradients for many temperate and boreal tree species. Despite their ecological and economic importance, little is known about the global changes in gene expression that accompany autumn cold acclimation in conifers. Using three populations of Sitka spruce (Picea sitchensis) spanning the species range, and a Picea cDNA microarray with 21,840 unique elements, within- and among-population gene expression was monitored during the autumn. Microarray data were validated for selected genes using real-time PCR. Similar numbers of genes were significantly twofold upregulated (1257) and downregulated (967) between late summer and early winter. Among those upregulated were dehydrins, pathogenesis-related/antifreeze genes, carbohydrate and lipid metabolism genes, and genes involved in signal transduction and transcriptional regulation. Among-population microarray hybridizations at early and late autumn time points revealed substantial variation in the autumn transcriptome, some of which may reflect local adaptation. These results demonstrate the complexity of cold acclimation in conifers, highlight similarities and differences to cold tolerance in annual plants, and provide a solid foundation for functional and genetic studies of this important adaptive process.