Next Article in Journal
Comparative Transcriptome Analyses of Gene Response to Different Light Conditions of Camellia oleifera Leaf Using Illumina and Single-Molecule Real-Time-Based RNA-Sequencing
Next Article in Special Issue
Complete Chloroplast Genome of Michelia shiluensis and a Comparative Analysis with Four Magnoliaceae Species
Previous Article in Journal
Whole-Genome DNA Methylation Associated With Differentially Expressed Genes Regulated Anthocyanin Biosynthesis Within Flower Color Chimera of Ornamental Tree Prunus mume
Previous Article in Special Issue
Bark Features for Identifying Resonance Spruce Standing Timber
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Forests
Volume 11
Issue 1
10.3390/f11010089
forests-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Constitutive and Cold Acclimation-Regulated Protein Expression Profiles of Scots Pine Seedlings Reveal Potential for Adaptive Capacity of Geographically Distant Populations

by
Danas Baniulis
1,*,
Monika Sirgėdienė
2,
Perttu Haimi
1,
Inga Tamošiūnė
1 and
Darius Danusevičius
2
1
Lithuanian Research Centre for Agriculture and Forestry, Babtai, 54333 Kaunas, Lithuania
2
Vytautas Magnus University Agriculture Academy, Institute of Forest Biology and Silviculture, Akademija, 53361 Kaunas, Lithuania
*
Author to whom correspondence should be addressed.
Forests 2020, 11(1), 89; https://doi.org/10.3390/f11010089
Submission received: 5 December 2019 / Revised: 6 January 2020 / Accepted: 8 January 2020 / Published: 10 January 2020
(This article belongs to the Special Issue Genetic and Phenotypic Variation in Tree Crops Biodiversity)
Browse Figures
Review Reports Versions Notes

Abstract

:
Geographically distant Scots pine (Pinus sylvestris L.) populations are adapted to specific photoperiods and temperature gradients, and markedly vary in the timing of growth patterns and adaptive traits. To understand the variability of adaptive capacity within species, molecular mechanisms that govern the physiological aspects of phenotypic plasticity should be addressed. Protein expression analysis is capable of depicting molecular events closely linked to phenotype formation. Therefore, in this study, we used comparative proteomics analysis to differentiate Scots pine genotypes originating from geographically distant populations in Europe, which show distinct growth and cold adaptation phenotypes. Needles were collected from 3-month-old seedlings originating from populations in Spain, Lithuania and Finland. Under active growth-promoting conditions and upon acclimation treatment, 65 and 53 differentially expressed proteins were identified, respectively. Constitutive protein expression differences detected during active growth were associated with cell metabolism and stress response, and conveyed a population-specific adaptation to the distinct climatic conditions. Acclimation-induced protein expression patterns suggested the presence of a similar cold adaptation mechanism among the populations. Variation of adaptive capacity among the genotypes was potentially represented by a constitutive low level of expression of the Ser/Thr-protein phosphatase, the negative regulator of the adaptive response. Also, overall less pronounced acclimation-induced response in seedlings from the Spanish population was observed. Thus, our study demonstrates that comparative proteomic analysis of young conifer seedlings is capable of providing insights into adaptation processes at the cellular level, which could help to infer variability of adaptive capacity within the plant species.

1. Introduction

Sustainability of future forests is threatened by anthropogenic impacts, mainly due to pollution and climate change [1,2]. Therefore, special attention is being given to understanding of the molecular basis of adaptive traits in tree species (e.g., [3,4]). High genetic diversity, phenotypic plasticity and epigenetic mechanisms are the key factors of the adaptive capacity of forest tree populations (e.g., [5,6]). Trees have a low evolutionary rate in terms of nucleotide substitution and diversification rates as compared to herbaceous species [7]. On the other hand, tree populations are larger, less genetically structured and contain larger genetic diversity [7,8]. Most forest tree species show substantial variability in adaptive traits at levels much higher than is observed at neutral genetic markers and form a collection of highly differentiated populations with contrasting adaptation to local climates [9]. Intra-specific genetic variation of tree species is linked to the high adaptability under constantly changing environments [3] providing diversity and plasticity of genetic variants for evolution, where phenotypic plasticity is a property of a genotype to produce different phenotypes in response to particular environmental conditions that vary both temporally and geographically. Among the main evolutionary forces influencing the variability patterns, natural selection primarily acts to increase genetic differentiation among populations, while gene flow and phenotypic plasticity tend to reduce it [10]. In wide-spread wind pollinated conifers such as Scots pine, the effects of genetic drift and mutations are weak. Although natural variation in epigenetic marks and the relation to phenotypic traits in trees is still an under-explored area, epigenetic variation has been suggested to contribute to the phenotypic plasticity and adaptive potential of forest trees [11].
Gene expression is a dynamic process and depends on the specific tissue, developmental stage and environmental signals. Gene expression analysis at both transcriptome and proteome level could provide details about population variability, and the mechanisms involved in evolution and adaptation [12,13]. Protein expression is the final step in gene expression process, thus, proteomic analysis bridges the gap between the genotype and phenotype, providing a depiction of molecular events that are most directly linked to the phenotype at cellular or organism level. Quantitative variation of protein abundance is the result of a complex network of regulatory mechanisms that integrate genetic variation and interaction with the environment of individual plants [14,15]. A study on Arabidopsis ecotypes showed that the proteomics approach could constitute a powerful tool to mine the biodiversity between ecotypes of a single plant species [16]. Proteomics was used to study phenotypic and molecular response to water deficit of two Eucalyptus genotypes [17] and adaptive variation in Picea abies (L.) H.Karst. ecotypes [18,19].
Scots pine (Pinus sylvestris L.) is the most widespread coniferous tree species occupying an extensive area in Eurasia. The complex biogeographic history of Scots pine in Europe has been well studied [20,21,22]. At the neutral part of the genome, in open pollinated wide-spread conifers such as Scots pine, the phylogenic structure is weak and population differentiation is low (e.g., [23,24]). Based on the maternally inherited mitochondrial DNA markers, the European range of Scots pine is grouped into three major haplotypes: south western (refugia in Pyrenees), central (refugia in Balkan and Italy) and south eastern (refugium in the southern edge of the Ural mountains) [20,25,26,27]. Gene flow and isolation by distance were other major factors affecting population differentiation in the Scots pine that were defined using neutral DNA markers [28]. Numerous studies based on indoor and field trials revealed a significant clinal variation among the Scots pine populations over latitude gradients in adaptive traits such as phenology [29,30].
Timing of cold acclimation is of adaptive significance for northerly conifers such as the Scots pine [31,32,33]. The seedlings cease active growth and enter the endodormancy stage in response to elongating nights towards the autumn [34], followed by shoot hardiness development in response to chilling temperatures in the early autumn [35,36]. In trees, cold acclimation treatment results in complex structural, biochemical and genetic shifts from elongating meristems to completely dormant and frost hardy tissues [37]. In brief, it causes protoplasmic dehydration, concentration of solutes, reduction of osmotic potential, upregulation of abscisic acid, a change in carbohydrate metabolism along with changes in membrane permeability and fatty acid composition [38,39]. The most prominent changes are induced in the chloroplasts and manifest as structural reorganization of the organelar membranes [40,41]. At protein level, the increased amount of antioxidant enzymes [42,43,44] and accumulation of dehydrins mostly localized in chloroplasts and the mitochondria membrane system was reported during the acclimation [45]. Gene expression studies showed upregulation of enzymes involved with synthesis of sugars and sugar derivatives, dehydration stress and other stress tolerance enzymes, such as heat-shock proteins [46].
Owing to the photoperiod and temperature-triggered adaptation in the northern regions, conifers have population-specific requirement for critical night length to budset and chilling temperatures to induce variable levels of shoot frost hardiness in early autumn [36]. When transferred to a single location, the critical night length for bud set and shoot hardiness development vary markedly between geographically distant populations within a latitudinal gradient [29,32,47]. Therefore, after a cold acclimation treatment, markedly different gene expression profiles could be expected in seedlings from the geographically distant Scots pine populations.
To probe an adaptive capacity within the Scots pine species, our study assessed phenotypic traits and quantitative and qualitative protein expression patterns in the Scots pine seedlings originating from the geographically distant populations of northern and southern Europe. Protein expression patterns were assessed by two-dimensional protein electrophoresis. Firstly, differences among the three populations were assessed during the active growth period. Furthermore, an acclimation treatment of the seedlings was used to extend the analysis of protein expression profiles with specific emphasis on cold adaptation processes in the distinct Scots pine genotypes.

2. Materials and Methods

2.1. Seedling Production and Cold Acclimation Treatment

We used seeds of Scots pine collected from the three autochthonous populations representing the Mediterranean, temporal and boreal climatic zones in Europe: central Spain (Valsain Forest mountains, Sierra de Guadarrama province, 40°49′ lat., 3°57′ long., 1828 m.a.s.l. alt.), central Lithuania (Kazlu Ruda forest tact, 54°46′ lat., 23°32′ long., 79 m.a.s.l. alt.) and northern Finland (Suomussalmi forest, 65°07′ lat., 28°49′ long., 233 m.a.s.l. alt.), respectively (Figure 1). The seeds were collected in native Scots pine stands naturally-born over generations in large forest tracts, where the influence of non-autochthonous sources is negligible. These populations mainly represent a geographically large latitudinal gradient of genepools adapted to distinct photoperiodic and temperate requirements to initiate dormancy and cold hardiness processes that have a major effect on timing of cold acclimation of the seedlings. The seed lots were mixtures from more than 20 trees within a forest stand, which was considered as representative of a population.
Seeds were soaked overnight in sterile deionized water and surface sterilized for 30 min with 30% hydrogen peroxide. Floating seeds were discarded, and the remaining seeds were sown in sterile plastic containers (length 50 cm × width 17 cm × height 15 cm) on sterilized peat treated with chalk and soaked with sterile nutritive solution. Seeds were germinated and seedlings were grown in the climatic growth chamber LT-36VL (Percival) during a 3-month period. Constant temperature of 24 °C, 100% humidity and a long-day (16 h, 30 µmol m−2 s−1; fluorescent light 5500 K photoperiod) was used to promote active growth. Seedlings were watered at a two-day interval and fertilizer containing 44 g/L nitrogen, 11 g/L P2O5, 33 g/L K2O, 2.2–3.3 g/L MgO, was used every second week. Two containers for each population were used, and the containers were randomly relocated within the chamber once per week to reduce undesired external effects within in the chamber. After germination, seedlings were thinned out to maintain comparable plant density (approx. 100 seedlings per container). In 3 months after germination, samples for the active growth experimental group were collected, and the remaining seedlings were exposed to the 2-week cold acclimation treatment in the same growth chamber as follows: for one week, the plants were maintained at 10 °C 100% humidity and 8 h photoperiod, followed by one week of darkness with the same temperature and humidity conditions.

2.2. Seedling Morphology Measurements

Before the cold acclimation treatment, morphometric traits of 40 seedlings per population were measured to the nearest millimeter. To assess growth cessation during the acclimation treatment, the needle weight of pooled sample representing four repeats for the active growth and acclimated experimental groups of each population was measured. The pooled samples were prepared by combining equal amount of material from first crowns of needles from 25 seedlings.

2.3. Two-Dimensional Gel Electrophoresis and Protein Identification

Four biological repeats representing the active growth and acclimated experimental groups of each population were prepared as pooled needle samples as described in the previous section, and were frozen in liquid nitrogen and stored at −70 °C. Protein samples were prepared using phenol extraction and ammonium acetate precipitation, as described previously [48]. Briefly, samples were ground in liquid nitrogen and suspended in 500 μL of the extraction buffer (0.5 M Tris-HCl, pH −7.5, 0.1 M KCl, 0.05 M ethylenediaminetetraacetic acid, 0.7 M sucrose, 2% polyvinylpolypyrrolidone, 2% β-mercaptoethanol, 1 mM phenylmethanesulfonylfluoride). Equal volume of Tris-buffered phenol (pH −7.5) was added and incubated with shaking for 30 min at 4 °C. The tubes were centrifugated for 15 min at 15,000× g at 4 °C, the upper phase was recovered, and the phenol extraction was repeated. Proteins were precipitated by addition of equal volume of 3M ammonium acetate in methanol, incubation at −20 °C overnight and centrifugation as before. Protein pellets were washed twice by resuspending in cold methanol and once in cold acetone and centrifugation for 5 min as before. After aspiration of acetone, protein pellets were dried for few minutes and solubilised in two-dimensional electrophoresis lysis buffer (6 M urea, 2 M thiourea, 4% 3-cholamidopropyl dimethylammonio 1-propanesulfonate, 10 mM Tris-HCl, pH−8.5). Protein concentrations were measured using a Bradford assay [49]. Four biological repeats were prepared for each treatment. Internal standards were prepared from a pooled mixture of all protein extracts.
Protein separation and detection was performed using a differential gel electrophoresis procedure as described previously [50]. Briefly, sample aliquots of 50 mg were used for labeling with Cy3 and Cy5 fluorescent dyes, and the internal standard was labeled with Cy2 dye (Lumiprobe). For the preparative gel, 500 mg of unlabeled internal standard was mixed with 50 mg of the labeled internal standard. Proteins were applied to 24 cm immobilized pH gradient strips with a linear gradient of pH 4–7 and isoelectric focusing was performed with an Ettan IPGphor (GE Healthcare, Chicago, IL, USA). Further, the proteins were separated on 1-mm thick 12.5% polyacrylamide gels with Ettan DALTsix electrophoresis apparatus (GE Healthcare) and gels were scanned using a FLA 9000 fluorescence scanner (GE Healthcare). Relative protein quantifications were performed using DeCyder 2-D Differential Analysis Software, v.7.0 (GE Healthcare).
Preparative gel was fixed in 50% methanol and 10% acetic acid. Protein spots were excised manually and subjected to protein digestion with trypsin, according to a method described previously [51]. Protein digests were loaded and desalted on a 100 mm × 20 mm Acclaim PepMap C18 trap column and separated on a 75 mm × 150 mm Acclaim PepMap C18 column using an Ultimate 3000 rapid separation liquid chromatography system (Thermo-Scientific, Waltham, MA, USA), coupled to a Maxis G4 Q-TOF mass spectrometer detector with a Captive Spray nano-electrospray ionization source (Bruker Daltonics, Bremen, Germany). Peptide identification was performed using the MASCOT server (Matrix Science, Boston, MA, USA) against the Pinus taeda L. genome database [52]. The mass error tolerance for peptide matching was limited to 5 ppm. The threshold value for the identification of proteins was a Mascot score >50 and at least two peptides. Where more than one protein was identified from the same spot, the exponentially modified protein abundance index (emPAI) [53] was employed for the relative quantitation of the proteins. Spots were designated as protein mixture if the emPAI indexes were similar for several of the identified proteins.

2.4. Data Analysis

A t-test was used to assess to assess the differences in the seedling morphology traits among the populations (40 seedlings per population).
The Biological Variation Analysis module of the DeCyder software was used to match protein spots of the four biological repeats across different gels and to perform statistical analysis of differential protein expression. A two-way analysis of variance (ANOVA) test was used to assess statistically significant (p ≤ 0.01) population and cold acclimation treatment effects. The threshold value of at least 2-fold difference in protein abundance was used.
To find patterns in the expression profiles, a cluster analysis was performed using Euclidian distance matrix and complete linkage method at the Extended Data Analysis module of the DeCyder software.
Blast2GO software [54] was used for the annotation and gene ontology analysis of the identified protein sequences. The identified Pinus taeda sequences were searched with Basic Local Alignment Search Tool (BLAST) against Swissprot database and protein sequences and annotated with Gene Ontology (GO) terms.

3. Results

3.1. Seedling Morphology

After 3-month cultivation under growth-promoting conditions (16 h light, high humidity and constant temperature of 24 °C), the seedlings of the Finnish population had the shortest length of stem and cotyledons (Table 1). There was no significant difference in the seedling height between the Lithuanian and Spanish populations. Seedlings of the Spanish population had more cotyledons and produced the largest needle biomass, as estimated by fresh weight of the sampled needles. During the acclimation treatment, there was no significant change of needle fresh weight for the seedlings from the Spanish population, indicating growth cessation. Meanwhile, the needle mass increased approximately 1,5-fold for the two northern populations.

3.2. Differentially Expressed Proteins

The average number of detected protein spots was 1464 ± 145 per gel after alignment (Figure S1). Analysis of biological variance revealed statistically significant and >2-fold variation in the abundance of 76 proteoforms among the six experimental groups (non-aclimated and acclimated seedlings from the Spanish, Lithuanian and Finnish populations) (Figure 2 and Figure 3). The largest number of unique differentially expressed protein spots were a characteristic of the Finnish population (25 and 20 in the non-acclimated and acclimated experimental groups, respectively), followed by the Spanish population (15 and 19). The Lithuanian population had the lowest number of unique spots (2 and 1) (Figure 2a,b). The number of proteoforms differentially expressed between populations was similar for the non-acclimated and acclimated experimental group. The abundance of 53 proteoforms changed after cold acclimation treatment (Figure 2c), and approx. 72% of the detected changes were independent of the seedling origin.
Through liquid chromatography-tandem mass spectrometry fingerprinting of trypsin-digested peptides, 26 proteoforms were unequivocally identified and annotated (Table A1). Among the identified proteins, 14 were associated with metabolic process (0044237), 17 were associated with a related term of cellular metabolic process (0008152) and 15 were linked to response to stress (0006950). At a lower GO hierarchy level, the terms included small molecule metabolic process (0044281), cellular protein metabolic process (0044267), cellular component organization (0016043), regulation of biological process (0050789), oxidation-reduction process (0055114), defense response (0006952) and response to salt stress (0009651), as indicated in Figure 4. Three proteins annotated as thaumatin-like protein (PITA_000066768-RA), chloroplastic TIC 62 protein (PITA_000059572-RA) and NAD-binding Rossmann-fold superfamily protein (PITA_000059572-RA) (Table S1) had no GO terms assigned but they are known to be involved in stress response processes (Albrecht and Bowman, 2012; Dutta et al., 2009; Rahjam et al., 2007). Another 5 proteoforms were identified as proteins coded by the genes of unknown function, and 29 spots contained a mixture of 2–4 proteins (Table S1).

3.3. Protein Expression Patterns

Cluster analysis of variation of abundance of the 76 differentially expressed protein spots revealed protein expression patterns associated with pine seedling origin and/or their response to acclimation treatment (Figure 4). Cold acclimation treatment had little or no significant effect on the expression of the proteins assigned to the clusters 4–5, 7, 9 and 11–12 and these mostly reflected population-depended protein expression differences. The most prominent difference among the three populations was observed for the two proteoforms assigned to the cluster 12. The lactoylglutathione lyase had approximately 5- to 10-fold higher abundance in the seedlings from the Spanish population as compared to the other two populations. Four proteoforms included in the cluster 11 had higher abundance in the Spanish population (none of them could be identified). Furthermore, six proteoforms (including tubulin beta chain-like and tau class glutathione s-transferase (GST) were assigned to the cluster 9 and had higher abundance in the Spanish and Lithuanian populations.
In contrast, a group of 22 proteins (including calreticulin-3-like isoform x1, malate dehydrogenase (MDH), proteasome subunit beta type-5 (PSMB5), alcohol dehydrogenase, 40 s ribosomal protein sa-like, chaperone protein 1, Ser/Thr-protein phosphatase 2A (PP2A) catalytic subunit) had higher abundance either in the Finnish population or both northern populations and were assigned to the clusters 4, 5 and 7.
Expression pattern of the remaining clusters 1–3, 6, 8, 10 and 13 reflected population- and/or acclimation-dependent differences (Figure 4). Before the acclimation, the majority of the 18 proteoforms separated into the clusters 1, 3 and 13 (including CP29-like, thaumatine-like protein (TLP), cyanate hydratase and two proteoforms of pathogenesis-related 10 (PR-10) were significantly more abundant in the seedlings of the Finnish or both northern populations. These proteins were significantly upregulated in the acclimated seedlings from all the three populations, however, for cluster 3, the response was more prominent in the Finnish and/or in both northern populations. Cyanate hydratase (cluster 13) was upregulated approximately 10-fold in the acclimated seedlings of the Finnish population compared to the non-acclimated control, representing the strongest response to the acclimation treatment.
Among the 10 proteoforms included in cluster 2, only two significant differences (unidentified proteoforms 7 and 8) were detected among the populations before the acclimation treatment. All proteins of this cluster were upregulated by the acclimation treatment, and the response was more prominent in the seedlings for the northern populations, leading to increase of significant differences among the populations (protein disulfide-isomerase (PDI), nucleoside diphosphate kinase 1 (NDK1) and three unidentified proteoforms).
Furthermore, an abundance of 14 proteoforms included in the clusters 6, 8 and 10 was markedly decreased upon the acclimation treatment. The response to acclimation had similar effect on protein expression for all three populations, however, the protein abundance pattern was different. The northern populations had higher abundance of the guanosine-diphosphate (GDP)-mannose epimerase 2 and proteoform 49 among the three proteoforms of the cluster 6. Meanwhile the remaining proteins (clusters 8 and 10) including probable linoleate 9s-LOX 5, aconitase, TIC 62, oxygen-evolving enhancer protein (OEE), nicotinamide adenine dinucleotide (NAD)-binding RF superfamily protein, abscisic acid water deficit stress and ripening-like (ASR-like) and water deficit inducible lipoprotein 3 (LP3)-like, were more abundant in the samples from the Spanish population.
Notably, only a few proteoforms (including four identified proteins) showed a strictly population-specific regulation upon the acclimation treatment. Malate dehydrogenase was upregulated in both northern populations, although the acclimation-induced expression level in the Lithuanian population was still below the constitutive expression level in the Spanish population. Calreticulin-3-like isoform and alcohol dehydrogenase were upregulated in the Finnish and Lithuanian populations, respectively. Meanwhile, the catalytic subunit of Ser/Thr-PP2A was downregulated in the northern populations, and the protein abundance was reduced to the level that was comparable to the constitutive expression observed in the Spanish population.

4. Discussion

An adaptive capacity of species in response to specific environmental cues could be inferred from mechanistic information at molecular level and its integration with ecologically relevant phenotypes [55]. A previous study with conifer tree species by [19] has demonstrated that proteomic analysis of young seedlings imparts a variation of gene expression that could be of importance for adaptation to contrasting altitude conditions in the closely related Norway spruce ecotypes, and the adaptive capacity could be further explored by heat stress treatment of the seedlings [18]. Our study revealed a similar quantitative variation of protein expression profiles in the needles of young Scots pine seedlings from geographically distant populations, although it reflects a variation of different physiological traits of potentially adaptive significance that could be species specific. Conifer seedlings develop an acclimation capacity in response to a short photoperiod and low temperatures at an early stage of of development [56], therefore cold acclimation treatment of the seedlings have stimulated additional differences in protein expression profiles that reflect a variation in cold adaptation capacity of the pine genotypes.

4.1. Stress Response Proteins Upregulated in the Northern Populations under Active Growth Conditions

Under active growth promoting conditions, both northern populations and in particular the Finnish population had higher abundance of proteins involved in protein metabolism (synthesis or degradation), such as heat shock protein (HSP), calreticulin, structural constituents of 40S ribosomes and 26S proteasome (Figure 4), that are also known as constituents of the stress tolerance mechanism. The chaperone protein 1 is a 100 kDa HSP that is critical for thermotolerance [57] and plant development [58]. Studies with Arabidopsis mutants revealed that constitutive expression of Hsp100 leads to preformed adaptation that ensures effective protection under heat stress [59]. Calreticulin is another chaperone protein which folds newly synthesized proteins and is also involved in regulation of Ca2+ homeostasis in the endoplasmic reticulum [60]. Upregulation of calreticulin has been observed under salt and osmotic stress conditions in Arabidopsis [61], and a function of the signaling component involved in response to cold stress in rice was proposed for the protein [62]. The proteasome subunit beta type-5 is a component of the 26S proteasome that mediates ubiquitin-dependent proteolysis of proteins and plays an important role in the stress-induced degradation of misfolded proteins, as well as the regulation of stress response signaling [63]. The 40s ribosomal protein sa-like is a structural constituent of ribosomes that regulates protein synthesis. A function of the p40 protein homologue is crucial for the active tissue growth [64] and its transcriptional upregulation under salt stress has been described in Arabidopsis [65].
Furthermore, upregulation of the TLP and PR-10, has been unique to the Finnish population. TLP proteins belong to the PR5 protein family and are induced under biotic, osmotic or cold stress in a variety of plants [66,67,68]. It is also well established that the PR-10 proteins are involved in biotic and abiotic stress response [69]. In perennial plants, expression of the PR-10 is induced during winter dormancy in mulbery and it provides either cryoprotectection for freeze-labile enzymes or serves as a nitrogen-storage protein [70].
In addition, the seedlings from the Finnish population have had higher abundance of alcohol dehydrogenase, malate dehydrogenase, GDP-mannose-epimerase 2, and cyanate hydratase enzymes that are involved in the primary cell metabolism, but also play role in the stress tolerance. Upregulation of alcohol dehydrogenase is essential in anaerobic respiration under cold stress conditions [71,72]. Malate is important intermediate and regulator of metabolic processes. Plants have several malate dehydrogenases that catalyze conversion of malate and oxaloacetate coupled to reduction or oxidation of the NAD pool, and this function has been linked to an abiotic stress tolerance in alfalfa and apple [73,74]. Upregulation of the cyanate hydratase upon salt stress was observed in tomato [75] and Suaeda aegyptiaca (Hasselq.) Zohary plants [76]. An accumulation of cyanide compounds is associated with elevated levels of ethylene synthesis as a consequence of plant stress response, therefore cyanate hydratase is considered a key enzyme responsible for the cyanide detoxification [77,78]. GDP-mannose-epimerase converts GDP-D-mannose to GDP-L-galactose, a precursor of ascorbate and cell wall polysaccharides [79]. This branching metabolic pathway plays a dual role in neutralization of the reactive oxygen species [80] as well as plant development [81].
The seedlings of the northern populations showed shorter phenotype and accumulated less needle biomass as compared to the Spanish population under active growth conditions (Table 1). It could be presumed that the reduced growth and upregulation of the stress response proteins might be an indicator of onset of the cold acclimation processes triggered by photoperiod used for the active growth experiments (16 h light and 8 h dark), owing to the short night length critical for growth cessation for the northern populations [34]. However, the seedlings had a relatively high level of the catalytic subunit of PP2A enzyme under active growth conditions that was effectively downregulated by the acclimation treatment. The PP2A is involved in cold stress signaling as a negative modulator of adaptive responses [82], and transcript mRNA levels of the PP2A has been reduced in tomato and alfalfa plants upon cold stress treatment [83,84]. This would suggest that the seedlings of the northern populations maintained active growth under the conditions used in the experiment, and the observed growth and protein expression differences were rather a consequence of physiological differences among the genotypes. Higher abundance of the proteins involved in stress tolerance could be linked to a constitutive cold adaptation mechanism characteristic of the northern genotypes.

4.2. Stress Response Proteins Upregulated in the Lithuanian and/or Spanish Population under Active Growth Conditions

The clusters 9 through 12 shown in Figure 4 include 20 proteoforms that have been upregulated mostly in the Lithuanian and/or Spanish populations under active growth conditions. Lactoylglutathione lyase, also known as glyoxalase I, is highly expressed in the Spanish population (approx. 10 and 5-fold compared to the Lithuanian and Finnish populations, respectively). The enzyme detoxifies methylglyoxal, a cytotoxic compound that is produced through glycolysis intermediates and generates free radicals [85]. Levels of methylglyoxal have been shown to rise under salinity stress and are regulated by glutathione and glyoxalase I [86].
As it could be expected due to geographical origin, the Lithuanian population tends to follow a similar protein expression pattern as the Finnish population located further north. However, expression of the proteins assigned to the clusters 9 and 10 (including tubulin beta-chain like, tau class GST, chloroplastic OEE, ASR- and LP3-like) in the Lithuanian population is similar to the Spanish population. The tau class GST, is involved in the reactive oxygen species (ROS) detoxification reactions that conjugate reduced gluthatione with a variety of target compounds. Expression of GSTs in plants is responsive to both biotic and abiotic stresses [87] and overexpression of the enzyme leads to chilling and osmotic stress tolerance in tobacco seedlings [88]. Beta-tubulin is one of the two structural components of cell microtubules and it is intimately involved in regulation of cell morphogenesis. In addition, it has been shown that the ROS signaling-mediated rearrangement in tubulin cytoskeleton plays crucial part in cell adaptation to stress [89]. Further, the abscisic acid induced ASR-like protein plays an important role in drought and salinity stress [90], and overexpression of ASR genes results in an increased cold tolerance [91,92]. LP3-like protein is a water-deficit-induced protein, which is homologous to the ASR and accumulates in loblolly pine roots upon water-deficit stress [93]. OEE is involved in the assembly of the Photosystem II and has a protective effect on the manganese cluster [94]. Increased expression of OEE1 was shown under salt or water deficiency stress [95,96].
The observed variation of stress-related protein expression among the populations from the southern and northern boundaries of the continent indicates presence of the constitutive stress tolerance mechanisms that are pertinent to the climatic conditions of the geographically distant habitats. The Lithuanian population combines a mixture of the traits.

4.3. Effect of Acclimation Treatment on Protein Expression

Cold acclimation treatment revealed additional protein expression differences that could be separated into three major groups that include proteins either (i) upregulated or (ii) downregulated in seedlings of all the three populations, or (iii) regulated in a population-specific manner. These protein expression differences suggest the presence of distinct elements of the cold acclimation mechanism and/or variation of their regulation among the geographically distant pine populations.
The function of the first group of proteins was linked to stress response. HSC70-2 is a protein folding catalyst involved in protein stabilization under abiotic stress conditions [97]. PDI belongs to another well-known family of chaperones that mediates formation, isomerization and reduction/oxidation of disulfide bonds during protein synthesis [98]. Ribonucleoprotein chloroplastic-like is required for normal chloroplast development under cold stress conditions by stabilizing transcripts of mRNAs [99]. Nucleoside diphosphate kinase 1 is a cytosolic enzyme that is essential in homeostasis of cellular nucleoside triphosphate pools [100], and its role under cold stress conditions has been described [101]. The stress response related proteins PR-10, TLP and cyanate hydratase have been also highly upregulated upon acclimation treatment, however, their abundance in the Spanish and Lithuanian populations have not exceed the level observed in the Finnish population under active growth conditions.
Acclimation treatment downregulates chloroplastic Tic62, a component of Tic translocon that mediated translocation of nuclear-encoded precursor proteins across the inner envelope of chloroplasts, and chloroplastic OEE were downregulated upon the acclimation treatment. Further downregulation of several enzymes, including cytoplasmic aconitase involved in the Krebs cycle, probable linoleate 9S-lipoxygenase 5 of the oxylipin synthesis pathway and GDP-mannose-epimerase 2 enzyme involved in the ascorbate and cell wall polysaccharide synthesis might be the result of reduced metabolic and photosynthetic activity in the seedlings maintained under the acclimation conditions of low temperature and reduced light. However, seasonal changes in chloroplast ultrastructure [102] and depression of photosynthesis [103] is also a part of the winter adaptation process in conifers.
Abundance of ASR-like and homologous LP3-like protein have been also downregulated in the seedling needle samples. This may appear to contradict the previous studies demonstrating that the ASP and LP3 play essential role in stress response [90,91,93]. However, the study by Padmanabhan et al. [93] has shown that the transcript of the LP3 protein is preferentially induced in roots with a constitutive basal level of expression observed in stems and needles, which would explain low expression levels in our experimental set up.
Among the proteins regulated in a population-specific manner, PP2A, the negative modulator of adaptive response [82], has been downregulated, and MDH and calreticulin-3-like isoforms have been upregulated in the northern populations upon the acclimation treatment. In the Spanish population, the response to acclimation treatment has a distinctive reduced intensity compared to the northern populations.

5. Conclusions

Our study demonstrates that comparative proteomic analysis of young Scots pine seedlings is capable of reflecting on adaptive capacity among the geographically distant populations. Under active growth conditions, quantitative protein expression differences were related to biological processes of cell metabolism and stress response and could represent elements of the adaptive mechanism tailored to the specific growth environment. Furthermore, acclimation treatment revealed an additional stratum of protein expression differences that could provide insight into molecular basis of the cold adaptation capacity among the genotypes. As could be expected, the acclimation-induced protein expression patterns confirmed the presence of a similar cold acclimation mechanism among the pine populations, however the analysis also revealed a discrete variation of the adaptive traits. The most prominent differences were substantiated by the constitutive low levels of the PP2A enzyme expression and the overall less-pronounced response to the acclimation treatment of seedlings from the Spanish population. The identified proteins involved in the constitutive or acclimation-induced expression profiles presented interesting candidates for further studies on cold adaptation physiology in pine.

Supplementary Materials

The following are available online at https://www.mdpi.com/1999-4907/11/1/89/s1: Figure S1: Two-dimensional electrophoresis gel of the needle protein samples isolated from Scots pine seedlings; Table S1: Proteins differentially expressed in seedlings of the different populations of Scots pine under active growth conditions and upon acclimation treatment that were identified as unknown protein sequences of the Pinus taeda genome or contained several protein mixtures.

Author Contributions

Conceptualization, D.B. and D.D.; investigation, D.B., M.S., P.H., I.T.; writing—original draft preparation, D.B., D.D.; writing—review and editing, M.S., P.H., I.T.; funding acquisition, D.D. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Ministry of Science and Education of Lithuania, grant number VP1-3.1-ŠMM-08-K-01-025.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

Table
Table A1. Proteins differentially expressed in seedlings of the different populations of Scots pine under active growth conditions (white bars) and upon acclimation treatment (grey bars).
Table A1. Proteins differentially expressed in seedlings of the different populations of Scots pine under active growth conditions (white bars) and upon acclimation treatment (grey bars).
Spot No.AccessionMascot Score
/p. no.
/SC, %		Thr.
M.W.
/pI		Sequence DescriptionProtein Expression Profile
5PITA_000033197-RA a413
/7
/35.6		17.8
/6.2		pathogenesis-related 10 (PR-10) protein Forests 11 00089 i001
6PITA_000003453-RA a725
/8
/28.3		32.3
/4.8		29 kda ribonucleoprotein (CP29) chloroplastic-like Forests 11 00089 i002
9PITA_000000569-RA1258
/19
/32.8		64.4
/5.2		heat shock cognate 70 kda protein 2
(HSP 70-2)		 Forests 11 00089 i003
10PITA_000048817-RA1513
/23
/42.6		55.8
/4.8		protein disulfide-isomerase (PDI) Forests 11 00089 i004
15PITA_000069617-RA a231
/5
/35.8		16.3
/5.9		nucleoside diphosphate kinase 1
(NDK 1)		 Forests 11 00089 i005
20PITA_000066768-RA a348
/5
/33.6		25.1
/5.0		thaumatin-like partial Forests 11 00089 i006
27PITA_000060113-RA a479
/9
/25.9		49.8
/5.7		calreticulin-3-like isoform x1 Forests 11 00089 i007
29PITA_000042977-RA a828
/12
/35.9		37.7
/5.8		malate dehydrogenase (MDH) Forests 11 00089 i008
32PITA_000037611-RA a802
/11
/41.0		30.7
/6.2		proteasome subunit beta type-5 (PSMB5) Forests 11 00089 i009
36PITA_000004184-RA601
/8
/17.6		57.6
/7.2		alcohol dehydrogenase partial Forests 11 00089 i010
37PITA_000055561-RA a926
/15
/36.1		34.9
/5.0		40s ribosomal protein sa-like Forests 11 00089 i011
39PITA_000052060-RA1085
/18
/24.0		101.2
/5.7		chaperone protein 1 Forests 11 00089 i012
42PITA_000009823-RA a429
/8
/28.4		35.1
/4.9		Ser/Thr-protein phosphatase pp2a (PP2A) catalytic subunit Forests 11 00089 i013
48PITA_000019887-RA a477
/10
/28.9		44.9
/5.6		guanosine-diphosphate (GDP)-mannose -epimerase 2 Forests 11 00089 i014
52PITA_000011094-RA a673
/10
/15.1		105
/6.0		probable linoleate 9s-lipoxygenase (LOX) 5 Forests 11 00089 i015
53PITA_000020461-RA a912
/19
/20.0		111.7
/6.4		aconitase cytoplasmic Forests 11 00089 i016
55PITA_000016807-RA a619
/11
/40.8		29.9
/4.2		tubulin beta chain-like Forests 11 00089 i017
60PITA_000006003-RA a599
/12
/25.0		32.4
/6.7		tau class glutathione s-transferase (GST) Forests 11 00089 i018
61PITA_000059572-RA a1289
/18
/41.0		56.2
/5.4		protein Tic62 chloroplastic Forests 11 00089 i019
64PITA_000047327-RA a383
/7
/23.5		43.5
/6.0		oxygen-evolving enhancer (OEE) protein chloroplastic Forests 11 00089 i020
65PITA_000055142-RA a566
/9
/31.2		30.5
/8.8		nicotinamide adenine dinucleotide (NAD)-binding Rossmann-fold superfamily protein Forests 11 00089 i021
67PITA_000008754-RA a328.69
/4
/21.9		23.6
/5.1		abscisic acid water deficit stress and ripening (ASR) inducible-like partial Forests 11 00089 i022
68PITA_000022327-RA446
/5
/50.0		17.6
/5.6		water deficit inducible lipoprotein 3 (LP3)-like partial Forests 11 00089 i023
73PITA_000023770-RA a448
/9
/23.9		32.6
/5.0		lactoylglutathione lyase Forests 11 00089 i024
75PITA_000008205-RA a768
/14
/58.0		22.8
/8.8		PR-10 protein Forests 11 00089 i025
76PITA_000077110-RA a157
/3
/22.6		10.6
/5.1		cyanate hydratase Forests 11 00089 i026
a Sequences that were selected as most abundant sequences from protein mixture using emPAI index. Data presented in graphs as log mean and standard error of the mean of standardized abundance, the same letters indicate statistically significant (p < 0.01) differences between the three populations, the same number of stars indicate significant difference between non-acclimated and acclimated experimental groups within each population. Abbreviations: p. no.—peptide number; SC—sequence coverage; Theoretical MW/pI—theoretical molecular weight and pI values; SP, LT, FI—Scots pine populations from Spain, Lithuania and Finland.

References

  1. Paoletti, E.; Bytnerowicz, A.; Andersen, C.; Augustaitis, A.; Ferretti, M.; Grulke, N.; Gunthardt-Goerg, M.S.; Innes, J.; Johnson, D.; Karnosky, D.; et al. Impacts of air pollution and climate change on forest ecosystems—Emerging research needs. Sci. World J. 2007, 7 (Suppl. 1), 1–8. [Google Scholar] [CrossRef]
  2. Perie, C.; de Blois, S. Dominant forest tree species are potentially vulnerable to climate change over large portions of their range even at high latitudes. PeerJ 2016, 4, e2218. [Google Scholar] [CrossRef] [Green Version]
  3. Aspinwall, M.J.; Loik, M.E.; Resco de Dios, V.; Tjoelker, M.G.; Payton, P.R.; Tissue, D.T. Utilizing intraspecific variation in phenotypic plasticity to bolster agricultural and forest productivity under climate change. Plant Cell Environ. 2014, 38, 1752–1764. [Google Scholar] [CrossRef] [PubMed]
  4. Neale, D.B.; Ingvarsson, P.K. Population, quantitative and comparative genomics of adaptation in forest trees. Curr. Opin. Plant Biol. 2008, 11, 149–155. [Google Scholar] [CrossRef] [PubMed]
  5. Breda, N.; Huc, R.; Granier, A.; Dreyer, E. Temperate forest trees and stands under severe drought: A review of ecophysiological responses, adaptation processes and long-term consequences. Ann. For. Sci. 2006, 63, 625–644. [Google Scholar] [CrossRef] [Green Version]
  6. Kremer, A.; Ronce, O.; Robledo-Arnuncio, J.J.; Guillaume, F.; Bohrer, G.; Nathan, R.; Bridle, J.R.; Gomulkiewicz, R.; Klein, E.K.; Ritland, K.; et al. Long-distance gene flow and adaptation of forest trees to rapid climate change. Ecol. Lett. 2012, 15, 378–392. [Google Scholar] [CrossRef] [Green Version]
  7. Petit, R.J.; Hampe, A. Some evolutionary consequences of being a tree. Annu. Rev. Ecol. Evol. Syst. 2006, 37, 187–214. [Google Scholar] [CrossRef] [Green Version]
  8. Nybom, H. Comparison of different nuclear DNA markers for estimating intraspecific genetic diversity in plants. Mol. Ecol. 2004, 13, 1143–1155. [Google Scholar] [CrossRef] [PubMed]
  9. Savolainen, O.; Pyhäjärvi, T.; Knürr, T. Gene flow and local adaptation in trees. Annu. Rev. Ecol. Evol. Syst. 2007, 38, 595–619. [Google Scholar] [CrossRef]
  10. Eriksson, O. Evolutionary forces influencing variation among populations of Pinus sylvestris. Silva Fennica 1998, 32, 173–184. [Google Scholar] [CrossRef] [Green Version]
  11. Bräutigam, K.; Vining, K.J.; Lafon-Placette, C.; Fossdal, C.G.; Mirouze, M.; Marcos, J.G.; Fluch, S.; Fraga, M.F.; Guevara, M.Á.; Abarca, D.; et al. Epigenetic regulation of adaptive responses of forest tree species to the environment. Ecol. Evol. 2013, 3, 399–415. [Google Scholar] [CrossRef] [PubMed]
  12. Diz, A.P.; Martinez-Fernandez, M.; Rolan-Alvarez, E. Proteomics in evolutionary ecology: Linking the genotype with the phenotype. Mol. Ecol. 2012, 21, 1060–1080. [Google Scholar] [CrossRef] [PubMed]
  13. Karr, T.L. Application of proteomics to ecology and population biology. Heredity 2008, 100, 200–206. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  14. Rockman, M.V.; Kruglyak, L. Genetics of global gene expression. Nat. Rev. Genet. 2006, 7, 862–872. [Google Scholar] [CrossRef]
  15. Thiellement, H.; Zivy, M.; Plomion, C. Combining proteomic and genetic studies in plants. J. Chromatogr. B 2002, 782, 137–149. [Google Scholar] [CrossRef]
  16. Chevalier, F.; Martin, O.; Rofidal, V.; Devauchelle, A.D.; Barteau, S.; Sommerer, N.; Rossignol, M. Proteomic investigation of natural variation between Arabidopsis ecotypes. Proteomics 2004, 4, 1372–1381. [Google Scholar] [CrossRef]
  17. Bedon, F.; Villar, E.; Vincent, D.; Dupuy, J.W.; Lomenech, A.M.; Mabialangoma, A.; Chaumeil, P.; Barre, A.; Plomion, C.; Gion, J.M. Proteomic plasticity of two Eucalyptus genotypes under contrasted water regimes in the field. Plant Cell Environ. 2012, 35, 790–805. [Google Scholar] [CrossRef]
  18. Valcu, C.M.; Lalanne, C.; Plomion, C.; Schlink, K. Heat induced changes in protein expression profiles of Norway spruce (Picea abies) ecotypes from different elevations. Proteomics 2008, 8, 4287–4302. [Google Scholar] [CrossRef]
  19. Valcu, C.M.; Lalanne, C.; Muller-Starck, G.; Plomion, C.; Schlink, K. Protein polymorphism between 2 Picea abies populations revealed by 2-dimensional gel rlectrophoresis and tandem mass spectrometry. J. Hered. 2008, 99, 364–375. [Google Scholar] [CrossRef] [Green Version]
  20. Cheddadi, R.; Vendramin, G.G.; Litt, T.; Francois, L.; Kageyama, M.; Lorentz, S.; Laurent, J.M.; De Beaulieu, J.L.; Sadori, L.; Jost, A.; et al. Imprints of glacial refugia in the modern genetic diversity of Pinus sylvestris. Glob. Ecol. Biogeogr. 2006, 15, 271–282. [Google Scholar] [CrossRef]
  21. Sinclair, W.T.; Morman, J.D.; Ennos, R. The postglacial history of Scots pine (Pinus sylvestris L.) in Western Europe: Evidence from mitochondrial DNA variation. Mol. Ecol. 1999, 8, 83–88. [Google Scholar] [CrossRef]
  22. Soranzo, N.; Alia, R.; Provan, J.; Powell, W. Patterns of variation at a mitochondrial sequence-tagged-site locus provides new insights into the postglacial history of European Pinus sylvestris populations. Mol. Ecol. 2000, 9, 1205–1211. [Google Scholar] [CrossRef]
  23. Dvornyk, V.; Sirvio, A.; Mikkonen, M.; Savolainen, O. Low nucleotide diversity at the pal1 locus in the widely distributed Pinus sylvestris. Mol. Biol. Evol. 2002, 19, 179–188. [Google Scholar] [CrossRef] [Green Version]
  24. Garcia-Gil, M.R.; Mikkonen, M.; Savolainen, O. Nucleotide diversity at two phytochrome loci along a latitudinal cline in Pinus sylvestris. Mol. Ecol. 2003, 12, 1195–1206. [Google Scholar] [CrossRef]
  25. Naydenov, K.; Senneville, S.; Beaulieu, J.; Tremblay, F.; Bousquet, J. Glacial vicariance in Eurasia: Mitochondrial DNA evidence from Scots pine for a complex heritage involving genetically distinct refugia at mid-northern latitudes and in Asia Minor. BMC Evol. Biol. 2007, 7, 233. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  26. Pyhajarvi, T.; Salmela, M.J.; Savolainen, O. Colonization routes of Pinus sylvestris inferred from distribution of mitochondrial DNA variation. Tree Genet. Genomes 2008, 4, 247–254. [Google Scholar] [CrossRef]
  27. Willis, K.J.; van Andel, T.H. Trees or no trees? The environments of central and eastern Europe during the last glaciation. Quat. Sci. Rev. 2004, 23, 2369–2387. [Google Scholar] [CrossRef]
  28. Scalfi, M.; Piotti, A.; Rossi, M.; Piovani, P. Genetic variability of Italian southern Scots pine (Pinus sylvestris L.) populations: The rear edge of the range. Eur. J. For. Res. 2009, 128, 377. [Google Scholar] [CrossRef]
  29. Hannerz, M. Genetic and Seasonal Variation in Hardiness and Growth Rhythm in Boreal and Temperate Conifers—A Review and Annotated Bibliography; The Forestry Research Institute of Sweden Report; SkogForsk: Uppsala, Sweden, 1998; pp. 1–140. [Google Scholar]
  30. Wachowiak, W.; Perry, A.; Donnelly, K.; Cavers, S. Early phenology and growth trait variation in closely related European pine species. Ecol. Evol. 2018, 8, 655–666. [Google Scholar] [CrossRef] [Green Version]
  31. Dietrichson, J. A summary of studies on genetic variation in forest trees grown in Scandinavia with special reference to the adaptation problem. Norway Skogforsoksv Medd 1971, 29, 25–56. [Google Scholar]
  32. Eriksson, G.; Andersson, S.; Eiche, V.; Ifver, J.; Persson, A. Severity index and transfer effects on survival and volume production of Pinus sylvestris in northern Sweden. Stud. For. Suecica 1980, 156, 1–32. [Google Scholar]
  33. Strimbeck, G.R.; Schaberg, P.G.; Fossdal, C.G.; Schroder, W.P.; Kjellsen, T.D. Extreme low temperature tolerance in woody plants. Front. Plant Sci. 2015, 6, 884. [Google Scholar] [CrossRef] [PubMed]
  34. Dormling, I. Bud dormancy, frost hardiness and frost drought in seedlings of Pinus sylvestris and Picea abies. In Advances in Cold Hardiness; Li, P.H., Christersson, L., Eds.; CRC Press: Boca Raton, FL, USA, 1993; pp. 285–298. [Google Scholar]
  35. Dormling, I. Influence of light intensity and temperature on photoperiodic response of Norway spruce provenances. In Proceedings of the IUFRO Meeting of WP Norway Spruce Provenances (S 2.03.11) and Norway Spruce Breeding (S 2.02.11), Bucharest, Romania, 24 September–1 October 1979; pp. 398–407. [Google Scholar]
  36. Repo, T. Seasonal changes in frost hardiness in Picea abies and Pinus sylvestris in Finland. Can. J. For. Res. 1992, 22, 1949–1957. [Google Scholar] [CrossRef]
  37. Wisniewski, M.; Nassuth, A.; Arora, R. Cold hardiness in trees: A mini-review. Front. Plant Sci. 2018, 9, 1394. [Google Scholar] [CrossRef]
  38. Kontunen-Soppela, S.; Lankila, J.; Lähdesmäki, L.; Laine, K. Response of protein and carbohydrate metabolism of Scots pine seedlings to low temperature. J. Plant Physiol. 2002, 159, 175–180. [Google Scholar] [CrossRef]
  39. Pukacki, P.M.; Kamińska-Rożek, E. Differential effects of spring reacclimation and deacclimation on cell membranes of Norway spruce seedlings. Acta Soc. Bot. Pol. 2013, 82, 77–84. [Google Scholar] [CrossRef] [Green Version]
  40. Vogg, G.; Heim, R.; Hansen, J.; Schäfer, C.; Beck, E. Frost hardening and photosynthetic performance of Scots pine (Pinus sylvestris L.) needles. I. Seasonal changes in the photosynthetic apparatus and its function. Planta 1998, 204, 193–200. [Google Scholar] [CrossRef]
  41. Vogg, G.; Heim, R.; Gotschy, B.; Beck, E.; Hansen, J. Frost hardening and photosynthetic performance of Scots pine (Pinus sylvestris L.) needles. II. Seasonal changes in the fluidity of thylakoid membranes. Planta 1998, 204, 201–206. [Google Scholar] [CrossRef]
  42. Pukacki, P.M.; Kamińska-Rożek, E. Reactive species, antioxidants and cold tolerance during deacclimation of Picea abies populations. Acta Physiol. Plant. 2013, 35, 129–138. [Google Scholar] [CrossRef]
  43. Bae, J.J.; Choo, Y.S.; Ono, K.; Sumida, A.; Hara, T. Photoprotective mechanisms in cold-acclimated and nonacclimated needles of Picea glehnii. Photosynthetica 2010, 48, 110–116. [Google Scholar] [CrossRef]
  44. Zheng, Y.; Yang, Q.; Xu, M.; Chi, Y.; Shen, R.; Li, P.; Dai, H. Responses of Pinus massoniana and Pinus taeda to freezing in temperate forests in central China. Scand. J. For. Res. 2012, 27, 520–531. [Google Scholar] [CrossRef]
  45. Korotaeva, N.; Romanenko, A.; Suvorova, G.; Ivanova, M.V.; Lomovatskaya, L.; Borovskii, G.; Voinikov, V. Seasonal changes in the content of dehydrins in mesophyll cells of common pine needles. Photosynth. Res. 2015, 124, 159–169. [Google Scholar] [CrossRef] [PubMed]
  46. Wisniewski, M.; Basset, C.I.; Gusta, L.V. An overview of cold hardiness in woody plants: Seeing the forest through the trees. Hortscience 2003, 38, 952–959. [Google Scholar] [CrossRef]
  47. Danusevicius, D.; Gabrilavicius, R. Variation in juvenile growth rhythm among Picea abies provenances from the Baltic states and adjacent regions. Scand. J. For. Res. 2001, 16, 305–317. [Google Scholar] [CrossRef]
  48. Isaacson, T.; Damasceno, C.M.; Saravanan, R.S.; He, Y.; Catala, C.; Saladie, M.; Rose, J.K. Sample extraction techniques for enhanced proteomic analysis of plant tissues. Nat. Protoc. 2006, 1, 769–774. [Google Scholar] [CrossRef]
  49. Bradford, M.M. A rapid and sensitive method for the quantitation of microgram quantities of protein utilizing the principle of protein-dye binding. Anal. Biochem. 1976, 72, 248–254. [Google Scholar] [CrossRef]
  50. Tamosiune, I.; Staniene, G.; Haimi, P.; Stanys, V.; Rugienius, R.; Baniulis, D. Endophytic Bacillus and Pseudomonas spp. modulate apple shoot growth, cellular redox balance, and protein expression under in vitro conditions. Front. Plant Sci. 2018, 9, 889. [Google Scholar] [CrossRef] [Green Version]
  51. Shevchenko, A.; Tomas, H.; Havlis, J.; Olsen, J.V.; Mann, M. In-gel digestion for mass spectrometric characterization of proteins and proteomes. Nat. Protoc. 2006, 1, 2856–2860. [Google Scholar] [CrossRef]
  52. Neale, D.B.; Wegrzyn, J.L.; Stevens, K.A.; Zimin, A.V.; Puiu, D.; Crepeau, M.W.; Cardeno, C.; Koriabine, M.; Holtz-Morris, A.E.; Liechty, J.D.; et al. Decoding the massive genome of loblolly pine using haploid DNA and novel assembly strategies. Genome Biol. 2014, 15, R59. [Google Scholar] [CrossRef] [Green Version]
  53. Ishihama, Y.; Oda, Y.; Tabata, T.; Sato, T.; Nagasu, T.; Rappsilber, J.; Mann, M. Exponentially modified protein abundance index (emPAI) for estimation of absolute protein amount in proteomics by the number of sequenced peptides per protein. Mol. Cell Proteom. 2005, 4, 1265–1272. [Google Scholar] [CrossRef] [Green Version]
  54. Conesa, A.; Gotz, S.; Garcia-Gomez, J.M.; Terol, J.; Talon, M.; Robles, M. Blast2GO: A universal tool for annotation, visualization and analysis in functional genomics research. Bioinformatics 2005, 21, 3674–3676. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  55. Schirmer, K.; Fischer, B.B.; Madureira, D.J.; Pillai, S. Transcriptomics in ecotoxicology. Anal. Bioanal. Chem. 2010, 397, 917–923. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  56. Bigras, F.J.; Ryyppo, A.; Lindstrom, A.; Stattin, E. Cold acclimation and deacclimation of shoots and roots of conifer seedlings. In Conifer Cold Hardiness; Bigras, F.J., Colombo, S.J., Eds.; Springer: Dordrecht, The Netherlands, 2001; pp. 57–88. [Google Scholar]
  57. Agarwal, M.; Katiyar-Agarwal, S.; Grover, A. Plant Hsp100 proteins: Structure, function and regulation. Plant Sci. 2002, 163, 397–405. [Google Scholar] [CrossRef]
  58. Singla, S.L.; Pareek, A.; Grover, A. Yeast HSP 104 homologue rice HSP 110 is developmentally- and stress-regulated. Plant Sci. 1997, 125, 211–219. [Google Scholar] [CrossRef]
  59. Hong, S.W.; Vierling, E. Hsp101 is necessary for heat tolerance but dispensable for development and germination in the absence of stress. Plant J. 2001, 27, 25–35. [Google Scholar] [CrossRef]
  60. Thelin, L.; Mutwil, M.; Sommarin, M.; Persson, S. Diverging functions among calreticulin isoforms in higher plants. Plant Signal Behav. 2011, 6, 905–910. [Google Scholar] [CrossRef] [Green Version]
  61. Ndimba, B.K.; Chivasa, S.; Simon, W.J.; Slabas, A.R. Identification of Arabidopsis salt and osmotic stress responsive proteins using two-dimensional difference gel electrophoresis and mass spectrometry. Proteomics 2005, 5, 4185–4196. [Google Scholar] [CrossRef]
  62. Komatsu, S.; Yang, G.; Khan, M.; Onodera, H.; Toki, S.; Yamaguchi, M. Over-expression of calcium-dependent protein kinase 13 and calreticulin interacting protein 1 confers cold tolerance on rice plants. Mol. Genet. Genom. 2007, 277, 713–723. [Google Scholar] [CrossRef]
  63. Stone, S.L. The role of ubiquitin and the 26S proteasome in plant abiotic stress signaling. Front. Plant Sci. 2014, 5, 135. [Google Scholar] [CrossRef]
  64. Garcia-Hernandez, M.; Davies, E.; Baskin, T.I.; Staswick, P.E. Association of plant p40 protein with ribosomes is enhanced when polyribosomes form during periods of active tissue growth. Plant Physiol. 1996, 111, 559–568. [Google Scholar] [CrossRef] [Green Version]
  65. Omidbakhshfard, M.A.; Omranian, N.; Ahmadi, F.S.; Nikoloski, Z.; Mueller-Roeber, B. Effect of salt stress on genes encoding translation-associated proteins in Arabidopsis thaliana. Plant Signal Behav. 2012, 7, 1095–1102. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  66. Anzlovar, S.; Dermastia, M. The comparative analysis of osmotins and osmotin-like PR-5 proteins. Plant Biol. 2003, 5, 116–124. [Google Scholar] [CrossRef]
  67. Liu, J.J.; Sturrock, R.; Ekramoddoullah, A.K. The superfamily of thaumatin-like proteins: Its origin, evolution, and expression towards biological function. Plant Cell Rep. 2010, 29, 419–436. [Google Scholar] [CrossRef]
  68. Zhu, B.; Chen, T.H.; Li, P.H. Expression of an ABA-responsive osmotin-like gene during the induction of freezing tolerance in Solanum commersonii. Plant Mol. Biol. 1993, 21, 729–735. [Google Scholar] [CrossRef]
  69. Liu, J.J.; Ekramoddoullah, A.K.M. The family 10 of plant pathogenesis-related proteins: Their structure, regulation, and function in response to biotic and abiotic stresses. Physiol. Mol. Plant Pathol. 2006, 68, 3–13. [Google Scholar] [CrossRef]
  70. Ukaji, N.; Kuwabara, C.; Takezawa, D.; Arakawa, K.; Fujikawa, S. Accumulation of pathogenesis-related (PR) 10/Bet v 1 protein homologues in mulberry (Morus bombycis Koidz.) tree during winter. Plant Cell Environ. 2004, 27, 1112–1121. [Google Scholar] [CrossRef]
  71. Christie, P.J.; Hahn, M.; Walbot, V. Low-temperature accumulation of alcohol dehydrogenase-1 mRNA and protein activity in maize and rice seedlings. Plant Physiol. 1991, 95, 699–706. [Google Scholar] [CrossRef] [Green Version]
  72. Jarillo, J.A.; Leyva, A.; Salinas, J.; Martinez-Zapater, J.M. Low temperature induces the accumulation of alcohol dehydrogenase mRNA in Arabidopsis thaliana, a chilling-tolerant plant. Plant Physiol. 1993, 101, 833–837. [Google Scholar] [CrossRef] [Green Version]
  73. Tesfaye, M.; Temple, S.J.; Allan, D.L.; Vance, C.P.; Samac, D.A. Overexpression of malate dehydrogenase in transgenic alfalfa enhances organic acid synthesis and confers tolerance to aluminum. Plant Physiol. 2001, 127, 1836–1844. [Google Scholar] [CrossRef]
  74. Yao, Y.X.; Dong, Q.L.; Zhai, H.; You, C.X.; Hao, Y.J. The functions of an apple cytosolic malate dehydrogenase gene in growth and tolerance to cold and salt stresses. Plant Physiol. Biochem. 2011, 49, 257–264. [Google Scholar] [CrossRef] [PubMed]
  75. Nveawiah-Yoho, P.; Zhou, J.; Palmer, M.; Sauve, R.; Zhou, S. Identification of proteins for salt tolerance using a comparative proteomics analysis of tomato accessions with contrasting salt tolerance. J. Am. Soc. Hortic. Sci. 2013, 138, 382–394. [Google Scholar] [CrossRef] [Green Version]
  76. Askari, H.; Edqvist, J.; Hajheidari, M.; Kafi, M.; Salekdeh, G.H. Effects of salinity levels on proteome of Suaeda aegyptiaca leaves. Proteomics 2006, 6, 2542–2554. [Google Scholar] [CrossRef]
  77. Miller, J.M.; Conn, E.E. Metabolism of hydrogen cyanide by higher plants. Plant Physiol. 1980, 65, 1199–1202. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  78. Peiser, G.D.; Wang, T.T.; Hoffman, N.E.; Yang, S.F.; Liu, H.W.; Walsh, C.T. Formation of cyanide from carbon 1 of 1-aminocyclopropane-1-carboxylic acid during its conversion to ethylene. Proc. Natl. Acad. Sci. USA 1984, 81, 3059–3063. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  79. Conklin, P.L.; Norris, S.R.; Wheeler, G.L.; Williams, E.H.; Smirnoff, N.; Last, R.L. Genetic evidence for the role of GDP-mannose in plant ascorbic acid (vitamin C) biosynthesis. Proc. Natl. Acad. Sci. USA 1999, 96, 4198. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  80. Akram, N.A.; Shafiq, F.; Ashraf, M. Ascorbic acid—A potential oxidant scavenger and its role in plant development and abiotic stress tolerance. Front. Plant Sci. 2017, 8, 613. [Google Scholar] [CrossRef] [PubMed]
  81. Qi, T.; Liu, Z.; Fan, M.; Chen, Y.; Tian, H.; Wu, D.; Gao, H.; Ren, C.; Song, S.; Xie, D. GDP-D-mannose epimerase regulates male gametophyte development, plant growth and leaf senescence in Arabidopsis. Sci. Rep. 2017, 7, 10309. [Google Scholar] [CrossRef] [Green Version]
  82. Luo, J.; Shen, G.; Yan, J.; He, C.; Zhang, H. AtCHIP functions as an E3 ubiquitin ligase of protein phosphatase 2A subunits and alters plant response to abscisic acid treatment. Plant J. 2006, 46, 649–657. [Google Scholar] [CrossRef]
  83. Monroy, A.; Sangwan, V.; Dhindsa, R.S. Low temperature signal transduction during cold acclimation: Protein phosphatase 2A as an early target for cold-inactivation. Plant J. 1998, 13, 653–660. [Google Scholar] [CrossRef] [Green Version]
  84. Pais, S.M.; Gonzalez, M.A.; Tellez-Inon, M.T.; Capiati, D.A. Characterization of potato (Solanum tuberosum) and tomato (Solanum lycopersicum) protein phosphatases type 2A catalytic subunits and their involvement in stress responses. Planta 2009, 230, 13–25. [Google Scholar] [CrossRef]
  85. Kalapos, M.P. The tandem of free radicals and methylglyoxal. Chem. Biol. Interact. 2008, 171, 251–271. [Google Scholar] [CrossRef] [PubMed]
  86. Yadav, S.K.; Singla-Pareek, S.L.; Ray, M.; Reddy, M.K.; Sopory, S.K. Methylglyoxal levels in plants under salinity stress are dependent on glyoxalase I and glutathione. Biochem. Biophys. Res. Commun. 2005, 337, 61–67. [Google Scholar] [CrossRef] [PubMed]
  87. Droog, F. Plant glutathione S-transferase, a tale of theta and tau. J. Plant Growth Regul. 1997, 16, 95–107. [Google Scholar] [CrossRef]
  88. Roxas, V.P.; Smith, R.K., Jr.; Allen, E.R.; Allen, R.D. Overexpression of glutathione S-transferase/glutathione peroxidase enhances the growth of transgenic tobacco seedlings during stress. Nat. Biotechnol. 1997, 15, 988–991. [Google Scholar] [CrossRef]
  89. Livanos, P.; Galatis, B.; Apostolakos, P. The interplay between ROS and tubulin cytoskeleton in plants. Plant Signal Behav. 2014, 9, e28069. [Google Scholar] [CrossRef]
  90. Golan, I.; Dominguez, P.G.; Konrad, Z.; Shkolnik-Inbar, D.; Carrari, F.; Bar-Zvi, D. Tomato ABSCISIC ACID STRESS RIPENING (ASR) gene family revisited. PLoS ONE 2014, 9, e107117. [Google Scholar] [CrossRef] [Green Version]
  91. Hsu, Y.F.; Yu, S.C.; Yang, C.Y.; Wang, C.S. Lily ASR protein-conferred cold and freezing resistance in Arabidopsis. Plant Physiol. Biochem. 2011, 49, 937–945. [Google Scholar] [CrossRef]
  92. Kim, S.J.; Lee, S.C.; Hong, S.K.; An, K.; An, G.; Kim, S.R. Ectopic expression of a cold-responsive OsAsr1 cDNA gives enhanced cold tolerance in transgenic rice plants. Mol. Cells 2009, 27, 449–458. [Google Scholar] [CrossRef]
  93. Padmanabhan, V.; Dias, D.M.; Newton, R.J. Expression analysis of a gene family in loblolly pine (Pinus taeda L.) induced by water deficit stress. Plant Mol. Biol. 1997, 35, 801–807. [Google Scholar] [CrossRef]
  94. Seidler, A.; Rutherford, A.W.; Michel, H. On the role of the N-terminus of the extrinsic 33 kDa protein of Photosystem II. Plant Mol. Biol. 1996, 31, 183–188. [Google Scholar] [CrossRef]
  95. Riccardi, F.; Gazeau, P.; Jacquemot, M.P.; Vincent, D.; Zivy, M. Deciphering genetic variations of proteome responses to water deficit in maize leaves. Plant Physiol. Biochem. 2004, 42, 1003–1011. [Google Scholar] [CrossRef] [PubMed]
  96. Sugihara, K.; Hanagata, N.; Dubinsky, Z.; Baba, S.; Karube, I. Molecular characterization of cDNA encoding oxygen evolving enhancer protein 1 increased by salt treatment in the mangrove Bruguiera gymnorrhiza. Plant Cell Physiol. 2000, 41, 1279–1285. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  97. Wang, W.; Vinocur, B.; Shoseyov, O.; Altman, A. Role of plant heat-shock proteins and molecular chaperones in the abiotic stress response. Trends Plant Sci. 2004, 9, 244–252. [Google Scholar] [CrossRef]
  98. Wilkinson, B.; Gilbert, H.F. Protein disulfide isomerase. Biochim. Biophys. Acta 2004, 1699, 35–44. [Google Scholar] [CrossRef]
  99. Kupsch, C.; Ruwe, H.; Gusewski, S.; Tillich, M.; Small, I.; Schmitz-Linneweber, C. Arabidopsis chloroplast RNA binding proteins CP31A and CP29A associate with large transcript pools and confer cold stress tolerance by influencing multiple chloroplast RNA processing steps. Plant Cell 2012, 24, 4266–4280. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  100. Bernard, M.A.; Ray, N.B.; Olcott, M.C.; Hendricks, S.P.; Mathews, C.K. Metabolic functions of microbial nucleoside diphosphate kinases. J. Bioenerg. Biomembr. 2000, 32, 259–267. [Google Scholar] [CrossRef] [PubMed]
  101. Lee, D.G.; Ahsan, N.; Lee, S.H.; Kang, K.Y.; Lee, J.J.; Lee, B.H. An approach to identify cold-induced low-abundant proteins in rice leaf. C. R. Biol. 2007, 330, 215–225. [Google Scholar] [CrossRef]
  102. Krol, M.; Hurry, V.M.; Maxwell, D.P.; Malek, L.; Ivanov, A.G.; Huner, N.P.A. Low growth temperature inhibition of photosynthesis in cotyledons of jack pine seedlings (Pinus banksianna) is due to impaired chloroplast development. Can. J. Bot. 2002, 80, 1042–1051. [Google Scholar] [CrossRef]
  103. Oquist, G.; Huner, N.P. Photosynthesis of overwintering evergreen plants. Annu. Rev. Plant Biol. 2003, 54, 329–355. [Google Scholar] [CrossRef]
Forests 11 00089 g001 550
Figure 1. Location of the Finnish (1), Lithuanian (2) and Spanish (3) populations used in the study.
Figure 1. Location of the Finnish (1), Lithuanian (2) and Spanish (3) populations used in the study.
Forests 11 00089 g001
Forests 11 00089 g002 550
Figure 2. Quantitative distribution of differentially expressed protein spots identified by two-dimensional gel electrophoresis analysis of needle protein samples isolated from seedlings originating from the Spanish (SP), Finish (FI) and Lithuanian (LT) Scots pine populations. Number of differentially expressed protein spots between the populations under growth promoting conditions (a) and upon acclimation (b) is shown. Number of upregulated and downregulated protein spots upon acclimation treatment are indicated by numbers separated by slash (c). The portion of the cycles overlapping between two or all three populations indicate a number of protein spots that are differentially expressed between the two or among all three populations, respectively.
Figure 2. Quantitative distribution of differentially expressed protein spots identified by two-dimensional gel electrophoresis analysis of needle protein samples isolated from seedlings originating from the Spanish (SP), Finish (FI) and Lithuanian (LT) Scots pine populations. Number of differentially expressed protein spots between the populations under growth promoting conditions (a) and upon acclimation (b) is shown. Number of upregulated and downregulated protein spots upon acclimation treatment are indicated by numbers separated by slash (c). The portion of the cycles overlapping between two or all three populations indicate a number of protein spots that are differentially expressed between the two or among all three populations, respectively.
Forests 11 00089 g002
Forests 11 00089 g003 550
Figure 3. Quantitative differences in protein abundance between non-aclimated and acclimated Scots pine seedlings originating from the geographically separated populations. Representative images of ten identified spots are shown. Complete gel image is presented in Figure S1.
Figure 3. Quantitative differences in protein abundance between non-aclimated and acclimated Scots pine seedlings originating from the geographically separated populations. Representative images of ten identified spots are shown. Complete gel image is presented in Figure S1.
Forests 11 00089 g003
Forests 11 00089 g004 550
Figure 4. Cluster analysis of abundance of the 76 protein spots differentially expressed in needles of the Scots pine seedlings among the three populations and/or under cold acclimation treatment. The standardized abundance scale is show at the lower left corner. Thirteen major clusters are marked by numbers. Star indicates spots that were identified as mixtures of 2–4 proteins (listed in Table S1). GO terms associated with the identified unique protein sequences are shown on the right (the terms corresponding to the indicated numbers are listed in the Results section).
Figure 4. Cluster analysis of abundance of the 76 protein spots differentially expressed in needles of the Scots pine seedlings among the three populations and/or under cold acclimation treatment. The standardized abundance scale is show at the lower left corner. Thirteen major clusters are marked by numbers. Star indicates spots that were identified as mixtures of 2–4 proteins (listed in Table S1). GO terms associated with the identified unique protein sequences are shown on the right (the terms corresponding to the indicated numbers are listed in the Results section).
Forests 11 00089 g004
Table
Table 1. Quantitative traits of pine seedlings before and after the cold acclimation treatment.
Table 1. Quantitative traits of pine seedlings before and after the cold acclimation treatment.
TraitPopulation
SpanishLithuanianFinnish
Number of cotyledons (4 weeks after germination)7.0 ± 0.12 a6.0 ± 0.13 b6.0 ± 0.13 b
Seedling height before the acclimation treatment (12 weeks after seeding), mm58.6 ± 0.23 a63.6 ± 0.23 a31.9 ± 0.11 b
Length of hypocotyl before the acclimation treatment (12 weeks after seeding), mm30.2 ± 0.16 a30.8 ± 0.16 a17.0 ± 0.06 b
Needle weight before the acclimation treatment (12 weeks after seeding), g4.2 ± 0.51 a3.2 ± 0.26 b1.5 ± 0.06 c
Needle weight at the end of acclimation treatment (14 weeks after seeding), g4.5 ± 0.16 a4.7 ± 0.14 a2.2 ± 0.34 b
Seedling height and cotyledon number were measured on 40 seedlings per population. Needle weight represents fresh weight of four pooled samples of 25 seedlings per population. The data are presented as mean and standard error of the mean. Means followed by the same letter are not significantly different (p < 0.05 by Tukey least significant difference (LSD) test).

Share and Cite

MDPI and ACS Style

Baniulis, D.; Sirgėdienė, M.; Haimi, P.; Tamošiūnė, I.; Danusevičius, D. Constitutive and Cold Acclimation-Regulated Protein Expression Profiles of Scots Pine Seedlings Reveal Potential for Adaptive Capacity of Geographically Distant Populations. Forests 2020, 11, 89. https://doi.org/10.3390/f11010089

AMA Style

Baniulis D, Sirgėdienė M, Haimi P, Tamošiūnė I, Danusevičius D. Constitutive and Cold Acclimation-Regulated Protein Expression Profiles of Scots Pine Seedlings Reveal Potential for Adaptive Capacity of Geographically Distant Populations. Forests. 2020; 11(1):89. https://doi.org/10.3390/f11010089

Chicago/Turabian Style

Baniulis, Danas, Monika Sirgėdienė, Perttu Haimi, Inga Tamošiūnė, and Darius Danusevičius. 2020. "Constitutive and Cold Acclimation-Regulated Protein Expression Profiles of Scots Pine Seedlings Reveal Potential for Adaptive Capacity of Geographically Distant Populations" Forests 11, no. 1: 89. https://doi.org/10.3390/f11010089

APA Style

Baniulis, D., Sirgėdienė, M., Haimi, P., Tamošiūnė, I., & Danusevičius, D. (2020). Constitutive and Cold Acclimation-Regulated Protein Expression Profiles of Scots Pine Seedlings Reveal Potential for Adaptive Capacity of Geographically Distant Populations. Forests, 11(1), 89. https://doi.org/10.3390/f11010089

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop