Next Article in Journal
Wastewater-Fertigated Short-Rotation Coppice, a Combined Scheme of Wastewater Treatment and Biomass Production: A State-of-the-Art Review
Next Article in Special Issue
Genetic Evaluation in Natural Populations of the Threatened Conifer Amentotaxus argotaenia (Hance) Pilg. (Taxaceae) Using Microsatellites
Previous Article in Journal
The Market Evolution of Medicinal and Aromatic Plants: A Global Supply Chain Analysis and an Application of the Delphi Method in the Mediterranean Area
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Heritable and Climatic Sources of Variation in Juvenile Tree Growth in an Austrian Common Garden Experiment of Central European Norway Spruce Populations

1
Gregor Mendel Institute, Austrian Academy of Sciences, 1030 Vienna, Austria
2
Department of Agrobiotechnology, Institute of Biotechnology in Plant Production, University of Natural Resources and Life Sciences, Vienna, 3430 Tulln an der Donau, Austria
3
Max Perutz Labs, University of Vienna, 1030 Vienna, Austria
*
Author to whom correspondence should be addressed.
Forests 2022, 13(5), 809; https://doi.org/10.3390/f13050809
Submission received: 27 April 2022 / Revised: 18 May 2022 / Accepted: 19 May 2022 / Published: 22 May 2022
(This article belongs to the Special Issue Genetic Variation and Phenotypic Plasticity in Trees)

Abstract

:
We leveraged publicly available data on juvenile tree height of 299 Central European Norway spruce populations grown in a common garden experiment across 24 diverse trial locations in Austria and weather data from the trial locations and population provenances to parse the heritable and climatic components of juvenile tree height variation. Principal component analysis of geospatial and weather variables demonstrated high interannual variation among trial environments, largely driven by differences in precipitation, and separation of population provenances based on altitude, temperature, and snowfall. Tree height was highly heritable and modeling the covariance between populations and trial environments based on climatic data led to more stable estimation of heritability and population × environment variance. Climatic similarity among population provenances was highly predictive of population × environment estimates for tree height.

1. Introduction

Under rapid anthropogenic climate change, ecosystem-dominant conifer species, such as Norway spruce (Picea abies), are facing environmental conditions to which they have not been adapted [1,2]. Although Norway spruce has broad genetic diversity and geographic adaptation, long generation times and the confounded effects of population structure and local adaptation make breeding populations for new environments challenging [3,4,5]. Common garden experiments can be used to evaluate a select number of diverse individuals or populations in different environments [6,7,8]. Assisted migration can then be applied, where well-adapted individuals are introduced into locations with compatible (current or forecasted) climates [8,9].
Precise estimation of genetic and environmental effects on tree growth phenotypes is critical for assisted migration and breeding program selection [8,9]. In addition to proper experimental design, germplasm choice, and phenotyping methods, statistical approaches can also increase the genetic signal and prediction accuracy of phenotypes [10,11,12]. For example, leveraging genetic marker, geospatial, and/or climatic data to model relationships between individuals, trial locations, and genotype × environment interactions is a common approach in agricultural and forest genetics and breeding [13,14,15,16,17].
In this study, we combined publicly available datasets on (a) juvenile tree height from a common garden experiment of Norway spruce, in which 299 Central European populations were evaluated across five years in 24 locations in Austria and (b) historical monthly precipitation and temperature measurements from the common garden trial locations and population provenances [6,18,19,20,21,22]. Previous studies have used this dataset to understand varying selection and intraspecific variation in climate response from an ecological perspective [6,18]. Here, we employed agricultural statistical methods to parse and to predict the heritable (genetic) and environmental components of phenotypic variance in tree height.

2. Materials and Methods

2.1. Common Garden Experiment

We used a publicly available dataset to investigate population and environmental effects on juvenile tree growth in Norway spruce (Picea abies) [6,18,23,24]. The plant material was derived from 299 Norway spruce populations collected across Central Europe during commercial seed harvests in 1971. Seeds from each population were germinated, and seedlings were then transplanted into a nursery in Austria in 1973. In 1978, a total of 65,534 five-year-old trees were transplanted into 24 trial locations across Austria in a balanced incomplete block design. On average, each population was grown in two locations, 27 populations were grown in each location, and 102 trees (replications) from each population were split into 3 randomized blocks within each location. The trials were conducted under natural conditions, receiving no management interventions.
The height (cm) of each tree was recorded at the ages of 7–10 and 15 years in 1980–1983 and 1988, respectively, resulting in a total of 300,310 observations. The previously published dataset also included the altitude, latitude, and longitude for all trial locations and for 278 population provenances.

2.2. Weather Data

We accessed monthly precipitation (total, liquid, solid) and temperature data from the publicly available Historical Instrumental Climatological Surface Time Series of the Greater Alpine Region (HISTALP) resource [19,20,21,22]. The HISTALP dataset includes monthly temperature and precipitation grids from 1780–2009 and 1801–2003, respectively, at 5 min × 5 min resolution ranging from 4–19° E and 43–46° N. We extracted data from the original Network Common Data Form (NetCDF) formatted HISTALP files using the “netcdf4” package in R [25,26].

2.3. Climatic Analysis

All altitude, latitude, longitude, and HISTALP weather variables used in the analyses described in this subsection were centered and scaled using the “scale” function in R [25]. We conducted principal component analysis (PCA) using the “FactoMineR” package in R [25,27]. To investigate spatiotemporal variation across the common garden experiment, we extracted monthly temperature and precipitation data from the grids closest to the trial locations in 1980–1983 and 1988 from the HISTALP dataset. We then conducted PCA using the altitude, latitude, longitude, and monthly weather data from each trial location in each year, henceforth referred to as trial environments (PCAEnv). To explore climatic differences among population provenances, we extracted monthly temperature and precipitation data from the grids closest to the population provenances from 1801 (earliest year with both temperature and precipitation data in HISTALP) to 1970 (year before seed harvests). Weather data from this time period should reflect the climates experienced by the parental populations of the trees evaluated in the common garden experiment. We then conducted a PCA using the altitude, latitude, longitude, and monthly weather data from each population provenance (PCAProv).
We used the altitude, latitude, longitude, and HISTALP weather data described above to model the relationships among population provenances and environments (trial location-years). We first calculated Euclidian distance matrices using said geospatial and weather variables for all pairs of population provenances ( D P r o v ) and trial environments ( D E n v ) with the “dist” function in R [25]. We then transformed the population provenance and trial environment distance matrices into population provenance ( S P r o v ) and trial environment ( S E n v ) similarity climatic matrices as S P r o v = max ( D P r o v ) D P r o v max ( D P r o v ) and S E n v = max ( D E n v ) D E n v max ( D E n v ) for further use in phenotypic modeling.

2.4. Phenotypic Analysis Using Climatic Data

Due to the large number of observations (N = 300,310) in the common garden experiment, it was not computationally feasible to specify complex variance–covariance matrices for the entire dataset in one model. As such, we used a two-stage approach common in multi-environment plant breeding and genetics studies [28,29], with mixed linear models fit with the “breedR” package in R [25,30]. In the first stage, we fit the following mixed model within each environment (location-year):
H e i g h t i j = μ + P o p i + B l o c k j + ε i j
where H e i g h t i j is the tree height response, μ is the overall mean, P o p i is the fixed effect of population i , B l o c k j is the random effect of block j , and ε i j is the random error. The random effects were independent and identically normally distributed as B l o c k j ~ N ( 0 , I σ b 2 ) and ε i j ~ N ( 0 , I σ ε 2 ) . Because height was measured on each tree within each location-year environment only once, a main effect for tree was not included in the model. We then extracted the fitted values ( y ^ i j ) for each population i within each block j from each within-environment Model (1) as:
y ^ i j = μ + p o p B L U E i + b l o c k B L U P j
where μ is the intercept, p o p B L U E i is the best linear unbiased estimate (BLUE) of population i , and b l o c k B L U P j is the best linear unbiased predictor (BLUP) of block j . We combined the fitted values across environments, resulting in a reduced and refined set of phenotypes for further analysis (N = 7314). Only fitted values for populations with climatic information were included in the second stage of the analysis (278/299 populations).
In the second stage, we fit a naïve across-year mixed model to understand the environmental and genetic (population-level) effects on juvenile tree growth as:
y ^ i j = μ + Y e a r + P o p i + E n v j + P o p E n v i j + ε i j
where y ^ i j is the fitted values from each within-environment Model (2), μ is the overall mean, Y e a r is the continuous fixed effect of year, P o p i is the random effect of population i , E n v j is the random effect of trial environment j , P o p E n v i j is the random effect of the interaction between population i and trial environment j , and ε i j is the random error. The random effects were assumed to be independent and identically normally distributed, where P o p i ~ N ( 0 , I σ p 2 ) , E n v j ~ N ( 0 , I σ e 2 ) , P o p E n v i j ~ N ( 0 , I σ p e 2 ) , and ε i j ~ N ( 0 , I σ ε 2 ) .
We also fit an across-year mixed model with form similar to Model (3) using historical weather data to model similarity among trial environments and populations provenances as:
y ^ i j = μ + Y e a r + P r o v i + E n v j + P r o v E n v i j + ε i j
where y ^ i j is the fitted values from each within-environment Model (2), μ is the overall mean, Y e a r is the continuous fixed effect of year, P r o v i is the random effect of population provenance i , E n v j is the random effect of trial environment j , P r o v E n v i j is the random effect of the interaction between population provenance i and trial environment j , and ε i j is the random error. The random effects were distributed as P r o v i ~ N ( 0 , S P r o v σ p 2 ) , E n v j ~ N ( 0 , S E n v σ e 2 ) , P r o v E n v i j ~ N ( 0 , I σ p e 2 ) , and ε i j ~ N ( 0 , I σ ε 2 ) , where S P r o v and S E n v are the population provenance and trial environment climatic similarity matrices, respectively.
We fit naïve within-year mixed models with the form in order to understand the environmental and genetic effects on juvenile tree height within developmental stage (tree age) as:
y ^ i j = μ + P o p i + E n v j + P o p E n v i j + ε i j
where y ^ i j is the fitted values from each within-environment Model (2) from a given year (1980, 1981, 1982, or 1983), μ is the overall mean, P o p i is the random effect of population i , E n v j is the random effect of trial environment j , P o p E n v i j is the random effect of the interaction between population i and trial environment j , and ε i j is the random error. The random effects were assumed to be independent and identically normally distributed, where P o p i ~ N ( 0 , I σ p 2 ) , E n v j ~ N ( 0 , I σ e 2 ) , P o p E n v i j ~ N ( 0 , I σ p e 2 ) , and ε i j ~ N ( 0 , I σ ε 2 ) .
We also fit a within-year mixed model with a form similar to Model (5) using environmental and population provenance weather data as:
y ^ i j = μ + P r o v i + E n v j + P r o v E n v i j + ε i j
where y ^ i j is the fitted values from each within-environment Model (2) from a given year (1980, 1981, 1982, or 1983), μ is the overall mean, P r o v i is the random effect of population provenance i , E n v j is the random effect of trial environment j , P r o v E n v i j is the random effect of the interaction between population provenance i and trial environment j , and ε i j is the random error. The random effects were distributed as P r o v i ~ N ( 0 , S P r o v σ p 2 ) , E n v j ~ N ( 0 , S E n v σ e 2 ) , P r o v E n v i j ~ N ( 0 , I σ p e 2 ) , and ε i j ~ N ( 0 , I σ ε 2 ) , where S P r o v and S E n v are the population provenance and trial environment climatic similarity matrices, respectively, for the provenances and environments observed within each year.
In Models (4) and (6), we used the term “population provenance” rather than “population” for two reasons: (i) to distinguish Models (4) and (6) from Models (3) and (5), and (ii) the population variance was modeled as the similarity matrix estimated from climatic data from the provenances (seed source locations) of the populations. Because of insufficient population × environment replication in 1988, we did not fit within-year models for 1988.
For each second stage Model (3)–(6), we estimated broad-sense heritability ( H 2 ) as:
H 2 = σ P 2 σ p 2 + ( σ p e 2 / n e ) + ( σ ε 2 / n e n r )
where n e and n r are the number of trial environments ( n e = 10 across years; n e = 2 within years) and replications per population/provenance per trial environment ( n r = 3), respectively, and σ p 2 , σ p e 2 , and σ ε 2 are the variance components for population provenance, population/provenance × environment, and error, respectively [31,32].

2.5. Predictions Using Climatic Covariance

Genomic prediction, in which phenotypes are predicted based on genetic similarity (kinship), is a widely used approach in plant and animal breeding for quantitative traits (controlled by many genes with small individual effects) [33]. In addition, genotype × environment interactions can be modeled with both genetic and weather-based covariance structures to improve genomic prediction accuracy [34,35,36,37]. We extracted the population, environment, and population × environment BLUPs from the naïve across-year phenotypic Model (3) and then fit the following mixed models:
p o p B L U P i = μ + P r o v i + ε i
e n v B L U P j = μ + E n v j + ε i
p o p E n v B L U P i = μ + P r o v i + ε i
p o p E n v B L U P j = μ + E n v j + ε i
p o p E n v B L U P i j = μ + P r o v i + E n v j + ε i
( p o p B L U P + p o p E n v B L U P ) i = μ + P r o v i + ε i
( p o p B L U P + p o p E n v B L U P ) j = μ + E n v j + ε i
( p o p B L U P + p o p E n v B L U P ) i j = μ + P r o v i + E n v j + ε i
where p o p B L U P i is the population BLUPs from Model (3) as the response for Model (8); e n v B L U P j is the environment BLUPs from Model (3) as the response for Model (9); p o p E n v B L U P i , p o p E n v B L U P j , and p o p E n v B L U P i j are the population × environment BLUPs from Model (3) as the response for Models (10)–(12), respectively; ( p o p B L U P + p o p E n v B L U P ) i , ( p o p B L U P + p o p E n v B L U P ) j , and ( p o p B L U P + p o p E n v B L U P ) i j are the sum of the population and population × environment BLUPs from Model (3) as the response for Models (13)–(15), respectively; μ is the overall mean for Models (8)–(15); P r o v i is the random effect of population provenance i in Models (8), (10), (12), (13), and (15); E n v j is the random effect of trial environment j in Models (9), (11), (12), (14), and (15); and ε i is the random error for Models (8)–(15). The random effects were distributed as P r o v i ~ N ( 0 , S P r o v σ p 2 ) , E n v j ~ N ( 0 , S E n v σ e 2 ) , and ε i ~ N ( 0 , I σ ε 2 ) . Each model was five-fold cross-validated with 20 replications. Predictive ability was estimated within each fold of each replication as the Pearson’s correlation between the observed and predicted values of the test set.

3. Results

3.1. Phenotypic and Climatic Variation among Trial Environments and Population Provenances

The trials were conducted across 24 locations in Austria, which varied with respect to altitude, latitude, longitude, and climate (Figure 1). In a PCA using geospatial and weather data from trial environments (location-years), the first and second principal components (PCs) explained 28% and 21% of the variance, respectively (Figure 1). The first PC demonstrated a separation between environments based on altitude and summer/winter temperature (Figure 1). On the second PC, environments were separated by year, spring/winter temperature, and summer precipitation (Figure 1).
Seeds from each population were sourced from 278 provenances across Central Europe with varying climate and geography (Figure 2). We conducted PCA for population provenances using geospatial and weather data. The first and second PCs explained 44% and 23% of the variance, respectively (Figure 2). The eigenvectors of the PCA predictor variables demonstrated that variation in PC 1 was driven by altitude, temperature, and snowfall (solid precipitation), while variation in PC 2 was driven by rainfall (liquid precipitation). Hierarchical clustering of the PCA revealed six clusters, which were geographically distributed as Carinthia, Austria, and Slovenia; Styria, Austria; Lower Austria, Germany, Czech Republic, Slovakia, and Poland; and three alpine regions (Figure 2).
Tree height, after controlling for population and spatial effects within trials, was normally distributed within years from 1980 to 1983 and positively skewed within 1988 and across years (1980–1988) (Figure 3). The coefficient of variation (CV) increased linearly with tree age, with the lowest CV in 1980 (seven-year-old trees, CV = 0.15) and the highest CV in 1988 (15-year-old trees, CV = 0.49) (Figure 3). The inter-population variance for tree height also varied by age, with higher variance components estimated with increasing age (1980 population variance = 14 cm2; 1983 population variance = 30 cm2).

3.2. Climatic Modeling Increases Heritable Signal for Juvenile Tree Growth

Within and across years, we found moderate to high heritability for tree height. In the “naïve” across-year phenotypic Model (3), we found that population and population × environment explained a smaller proportion of the total variance in tree height (9% each) than the environment, which accounted for 75% of the total variance (Table 1). Although population variance was low (9%), the highly replicated experimental design allowed for a high estimation of broad-sense heritability for tree height (H2 = 0.89) (Table 1). Year had an effect estimate of 13.46 ± 0.73 cm.
We also fit a “climatic” across-year phenotypic Model (4) in which the variance–covariance of the population provenance and environment random effects were modeled as the climatic similarity matrices between population provenances and trial environments, respectively. The climatic across-year model increased the proportion of the total variance in tree height explained by environment (92%) compared to the naïve across-year model (75%) (Table 1). Although population and population × environment variance were smaller (population = 5%, population × environment = 2%) than in the naïve across-year model (population = 9%; population × environment = 9%), the broad-sense heritability increased (H2naive = 0.89; H2climatic = 0.96) (Table 1). Year had no effect on tree height in the across-year climatic model.
In the naïve within-year analysis (Model (5)), the environmental variance increased while the broad-sense heritability and population variance tended to decrease over time from 1980 (seven-year-old trees, H2 = 0.81, genetic variance = 33%, environmental variance = 46%) to 1983 (10-year-old trees, H2 = 0.49, genetic variance = 4%, environmental variance = 83%) (Table 1). In contrast, specifying the variance-covariance of the population and environment terms with climatic similarity matrices (“climatic” Model (6)) led to more stable estimates of broad-sense heritability (H21980→1983 = 0.94→0.77) and population × environment variance (population × environment variance1980→1983 = 4→2%) over time (Table 1). Within-year climatic modeling (6) demonstrated an increase in environmental variance and a decrease in population variance over time (genetic variance1980→1983 = 36→4%, environmental variance1980→1983 = 58→92%) similar to that of the naïve models (Table 1).

3.3. Climate Explains Environmental Variance in Tree Height

Because the common garden experiment consisted of even-aged stands planted in the same year across trial locations, tree age (physiological growth stage) and trial environment (location-year) may have been confounded. As such, we modeled the environment BLUPs for tree height with climatic similarity between trial environments (Model (9)). The environmental similarity matrix, estimated from climatic and geospatial data across trial location-year environments, explained 71% of the variance in environmental BLUP estimates for tree height and had a predictive ability of r = 0.48 (Table 2). Complementarily, climatic phenotypic modeling including climatic variance–covariance matrices within (Model (6)) and across (Model (4)) years increased the variance explained by the environment term (Table 2) and removed the effect of year on tree height when compared the naïve phenotypic modeling (Models (3) and (5)). Combined, these results suggest that the variation in tree height estimates among trial location-year environments was largely driven by climatic variation rather than by the confounding effects of tree growth stage.

3.4. Provenance Climate Is Highly Predictive of Inter-Population and Population × Environment Variation in Tree Height

The population provenance similarity matrix, estimated from historical climatic and geospatial data across population provenances, demonstrated strong population structure (Figure 1) and explained 45% of the variance in populations BLUPs for tree height (Model (8)), with a predictive ability of r = 0.44 (Table 2). Although population × environment BLUPs were poorly modeled (r < 0; explained variance ≤ 1%) by provenance and/or environment climatic similarity (Models (10)–(12)), provenance climatic similarity explained nearly all of the variance (>91%) in and was highly predictive (r = 0.68) of the sum of population BLUPs plus population × environment BLUPs (Models (13) and (15)) (Table 2). In addition, climatic phenotypic modeling (Models (4) and (6)) demonstrated more stable estimation of population × environment effects when compared to the naïve phenotypic Models (3) and (5). These results indicate a strong signal of local adaptation, where populations from provenances with similar geography/climate have similar tree growth, regardless of the trial environment in which they are observed.

4. Discussion

We used agricultural statistical methods to dissect the heritable and environmental components of juvenile tree growth in Central European Norway Spruce. The highly replicated design of the common garden experiment allowed for a precise estimation of the genetic (population and population × environment) and environmental (trial location-year) effects. Population explained a relatively small proportion of the total variance in tree height compared to environment, similar to previous findings in this dataset [18,38,39] and consistent with low estimates of population differentiation between Norway spruce populations [38,40]. However, broad-sense heritability was moderate to high and previous studies on tree height in clonal and progeny trials of Norway spruce [37,39] have reported similar heritability estimates. Modeling the variance–covariance of the population provenance and environment terms using climate data led to higher estimates of heritability and environmental variance and more stable estimation of heritability and population × environment for tree height across developmental time. The higher heritability estimates with climatic modeling were likely the result of increased environmental signal and subsequently decreased error variance. Tree height becomes less predictive of total biomass with age, and growth habits (branching pattern, etc.) may be a better indicator of tree growth after the transition from the juvenile to adult stages [41], which may partially explain the reduction in heritability and genetic variance for tree height over time in this study. Drought conditions in Central Europe increased in severity from the beginning to the end of the trials (1983–1988), with a large historical drought event beginning in 1987 [42], supporting the PCA separation between the 1988 and 1981–1983 trial years reported here. In addition, the increasing environmental variance for tree height with developmental age may be partially explained by increasing drought severity across the trial period.
Although genetic marker data were not available in this dataset, we modeled inter-population relationships using historical climatic data from the population provenances. Plant and animal breeders frequently use climate and/or genetic marker data to model relationships among individuals and/or environments in order to improve the prediction accuracy of phenotypes in new (genetically related) material and even in new environments [13,14,15,34,35,36,37,43,44]. Previous studies have demonstrated that local adaptation is prevalent in tree species, including Norway spruce, and that ancestral environment is a strong predictor of growth phenotypes in Norway spruce [8,45,46,47,48]. We found that climatic relationships among population provenances (a) were highly predictive of population + population × environment tree height variation and (b) demonstrated population structure nested within geography, indicating that these populations have strong local adaptation [8]. Our findings complement previous results from ecological modeling of the dataset, where within-population variation was driven by population provenance temperature [18]. The modeling used in this study results in refined tree growth estimates for population and population × environment, which could be used to identify potential sources of material adapted to new environments.

Author Contributions

Conceptualization, L.M.; data curation, L.M.; formal analysis, L.M.; funding acquisition, K.S.; methodology, L.M.; project administration, L.M.; software, L.M.; validation, L.M.; visualization, L.M.; writing—original draft, L.M.; writing—review and editing, L.M. and K.S. All authors have read and agreed to the published version of the manuscript.

Funding

This work was partially supported by funding from the Austrian Academy of Sciences.

Data Availability Statement

This study used publicly available datasets. Data from the common garden experiment can be accessed at https://doi.org/10.5061/dryad.877ts. Climate data can be accessed at http://www.zamg.ac.at/histalp/dataset/grid/five_min.php (accessed on 26 April 2022).

Acknowledgments

Open Access Funding by the University of Vienna.

Conflicts of Interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

References

  1. Lindner, M.; Maroschek, M.; Netherer, S.; Kremer, A.; Barbati, A.; Garcia-Gonzalo, J.; Seidl, R.; Delzon, S.; Corona, P.; Kolström, M.; et al. Climate Change Impacts, Adaptive Capacity, and Vulnerability of European Forest Ecosystems. For. Ecol. Manag. 2010, 259, 698–709. [Google Scholar] [CrossRef]
  2. Seidl, R.; Thom, D.; Kautz, M.; Martin-Benito, D.; Peltoniemi, M.; Vacchiano, G.; Wild, J.; Ascoli, D.; Petr, M.; Honkaniemi, J.; et al. Forest Disturbances under Climate Change Europe PMC Funders Group. Nat Clim Chang. 2017, 7, 395–402. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  3. Leslie, A.B.; Beaulieu, J.M.; Rai, H.S.; Crane, P.R.; Donoghue, M.J.; Mathews, S. Hemisphere-Scale Differences in Conifer Evolutionary Dynamics. Proc. Natl. Acad. Sci. USA 2012, 109, 16217–16221. [Google Scholar] [CrossRef] [Green Version]
  4. Wang, X.; Bernhardsson, C.; Ingvarsson, P.K. Demography and Natural Selection Have Shaped Genetic Variation in the Widely Distributed Conifer Norway Spruce (Picea abies). Genome Biol. Evol. 2020, 12, 3803–3817. [Google Scholar] [CrossRef] [PubMed]
  5. Yeaman, S.; Hodgins, K.A.; Lotterhos, K.E.; Suren, H.; Nadeau, S.; Degner, J.C.; Nurkowksi, K.A.; Smets, P.; Wang, T.; Gray, L.K.; et al. Convergent Local Adaptation to Climate in Distantly Related Conifers. Science 2016, 353, 1431–1433. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  6. Kapeller, S.; Lexer, M.J.; Geburek, T.; Hiebl, J.; Schueler, S. Intraspecific Variation in Climate Response of Norway Spruce in the Eastern Alpine Range: Selecting Appropriate Provenances for Future Climate. For. Ecol. Manag. 2012, 271, 46–57. [Google Scholar] [CrossRef]
  7. Oleksyn, J.; Modrzyński, J.; Tjoelker, M.G.; Zytkowiak, R.; Reich, P.B.; Karolewski, P. Growth and Physiology of Picea abies Populations from Elevational Transects: Common Garden Evidence for Altitudinal Ecotypes and Cold Adaptation. Funct. Ecol. 1998, 12, 573–590. [Google Scholar] [CrossRef]
  8. Aitken, S.N.; Bemmels, J.B. Time to Get Moving: Assisted Gene Flow of Forest Trees. Evol. Appl. 2016, 9, 271–290. [Google Scholar] [CrossRef]
  9. Williams, M.I.; Dumroese, R.K. Preparing for Climate Change: Forestry and Assisted Migration. J. For. 2013, 111, 287–297. [Google Scholar] [CrossRef]
  10. Isidro y Sánchez, J.; Akdemir, D. Training Set Optimization for Sparse Phenotyping in Genomic Selection: A Conceptual Overview. Front. Plant Sci. 2021, 12, 715910. [Google Scholar] [CrossRef]
  11. Heslot, N.; Jannink, J.-L.; Sorrells, M.E. Perspectives for Genomic Selection Applications and Research in Plants. Crop Sci. 2015, 55, 1–12. [Google Scholar] [CrossRef]
  12. Crossa, J.; Pérez-Rodríguez, P.; Cuevas, J.; Montesinos-López, O.; Jarquín, D.; de los Campos, G.; Burgueño, J.; González-Camacho, J.M.; Pérez-Elizalde, S.; Beyene, Y.; et al. Genomic Selection in Plant Breeding: Methods, Models, and Perspectives. Trends Plant Sci. 2017, 22, 961–975. [Google Scholar] [CrossRef]
  13. Heslot, N.; Akdemir, D.; Sorrells, M.E.; Jannink, J.L. Integrating Environmental Covariates and Crop Modeling into the Genomic Selection Framework to Predict Genotype by Environment Interactions. Theor. Appl. Genet. 2014, 127, 463–480. [Google Scholar] [CrossRef] [PubMed]
  14. Bustos-Korts, D.; Boer, M.P.; Malosetti, M.; Chapman, S.; Chenu, K.; Zheng, B.; van Eeuwijk, F.A. Combining Crop Growth Modeling and Statistical Genetic Modeling to Evaluate Phenotyping Strategies. Front. Plant Sci. 2019, 10, 1491. [Google Scholar] [CrossRef] [PubMed]
  15. Crossa, J.; De Los Campos, G.; Maccaferri, M.; Tuberosa, R.; Burgueño, J.; Pérez-Rodríguez, P. Extending the Marker × Environment Interaction Model for Genomic-Enabled Prediction and Genome-Wide Association Analysis in Durum Wheat. Crop Sci. 2016, 56, 2193–2209. [Google Scholar] [CrossRef] [Green Version]
  16. Avanzi, C.; Piermattei, A.; Piotti, A.; Büntgen, U.; Heer, K.; Opgenoorth, L.; Spanu, I.; Urbinati, C.; Vendramin, G.G.; Leonardi, S. Disentangling the Effects of Spatial Proximity and Genetic Similarity on Individual Growth Performances in Norway Spruce Natural Populations. Sci. Total Environ. 2019, 650, 493–504. [Google Scholar] [CrossRef]
  17. Rodríguez-Álvarez, M.X.; Boer, M.P.; van Eeuwijk, F.A.; Eilers, P.H.C. Correcting for Spatial Heterogeneity in Plant Breeding Experiments with P-Splines. Spat. Stat. 2018, 23, 52–71. [Google Scholar] [CrossRef] [Green Version]
  18. Kapeller, S.; Dieckmann, U.; Schueler, S. Varying Selection Differential throughout the Climatic Range of Norway Spruce in Central Europe. Evol. Appl. 2017, 10, 25–38. [Google Scholar] [CrossRef] [Green Version]
  19. Chimani, B.; Matulla, C.; Böhm, R.; Hofstätter, M. A New High Resolution Absolute Temperature Grid for the Greater Alpine Region Back to 1780. Int. J. Climatol. 2013, 33, 2129–2141. [Google Scholar] [CrossRef]
  20. Chimani, B.; Böhm, R.; Matulla, C.; Ganekind, M. Development of a Longterm Dataset of Solid/Liquid Precipitation. Adv. Sci. Res. 2011, 6, 39–43. [Google Scholar] [CrossRef]
  21. Hiebl, J.; Auer, I.; Böhm, R.; Schöner, W.; Maugeri, M.; Lentini, G.; Spinoni, J.; Brunetti, M.; Nanni, T.; Tadić, M.P.; et al. A High-Resolution 1961-1990 Monthly Temperature Climatology for the Greater Alpine Region. Meteorol. Z. 2009, 18, 507–530. [Google Scholar] [CrossRef]
  22. Efthymiadis, D.; Jones, P.D.; Briffa, K.R.; Auer, I.; Böhm, R.; Schöner, W.; Frei, C.; Schmidli, J. Construction of a 10-Min-Gridded Precipitation Data Set for the Greater Alpine Region for 1800-2003. J. Geophys. Res. Atmos. 2006, 111, 6120. [Google Scholar] [CrossRef] [Green Version]
  23. Nather, J.; Holzer, K. Über Die Bedeutung Und Die Anlage von Kontrollflächen Zur Prüfung von Anerkanntem Fichtenpflanzgut. In Informationsdienst Forstliche Bundesversuchsanstalt Wien 181; Forstliche Bundesversuchsanstalt: Vienna, Austria, 1979. [Google Scholar]
  24. Schulze, U. Fichtenherkunftsprüfung 1978. In Internationaler Verband Forstlicher Versuchsanstalten IUFRO, Arbeitsgruppe S 2.02-11 Fichtenherkünfte; Forstliche Bundesversuchsanstalt: Vienna, Austria, 1985; pp. 1–12. [Google Scholar]
  25. R Core Team. R: A Language and Environment for Statistical Computing 2020; R Foundation for Statistical Computing: Vienna, Austria, 2013; Software; Available online: https://www.R-project.org/ (accessed on 20 May 2022).
  26. Pierce, D. ncdf4: Interface to Unidata NetCDF (Version 4 or Earlier) Format Data Files 2021. Software. Available online: http://cirrus.ucsd.edu/~pierce/ncdf/ (accessed on 20 May 2022).
  27. Lê, S.; Josse, J.; Husson, F. FactoMineR: An R Package for Multivariate Analysis. J. Stat. Softw. 2008, 25, 1–18. [Google Scholar] [CrossRef] [Green Version]
  28. Piepho, H.P.; Möhring, J.; Schulz-Streeck, T.; Ogutu, J.O. A Stage-Wise Approach for the Analysis of Multi-Environment Trials. Biometrical J. 2012, 54, 844–860. [Google Scholar] [CrossRef] [PubMed]
  29. Piepho, H.P.; Tempelman, R.J.; Williams, E.R. Guest Editors’ Introduction to the Special Issue on “Recent Advances in Design and Analysis of Experiments and Observational Studies in Agriculture”. J. Agric. Biol. Environ. Stat. 2020, 25, 453–456. [Google Scholar] [CrossRef]
  30. Muñoz, F.; Sanchez, L. BreedR: Statistical Methods for Forest Genetic Resources Analysts 2020. Software. Available online: http://famuvie.github.io/breedR/ (accessed on 20 May 2022).
  31. Piepho, H.P.; Moehring, J. Computing Heritability and Selection Response from Unbalanced Plant Breeding Trials. Genetics 2007, 177, 1881–1888. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  32. Holland, J.B.; Nyquist, W.E.; Cervantes-Martinez, C.T. Estimating and Interpreting Heritability for Plant Breeding. Plant Breed. Rev. 2003, 22, 9–112. [Google Scholar]
  33. Hickey, J.M.; Chiurugwi, T.; Mackay, I.; Powell, W. Genomic Prediction Unifies Animal and Plant Breeding Programs to Form Platforms for Biological Discovery. Nat. Genet. 2017, 49, 1297–1303. [Google Scholar] [CrossRef]
  34. Gillberg, J.; Marttinen, P.; Mamitsuka, H.; Kaski, S. Modelling G3E with Historical Weather Information Improves Genomic Prediction in New Environments. Bioinformatics 2019, 35, 4045–4052. [Google Scholar] [CrossRef] [Green Version]
  35. Jarquín, D.; Crossa, J.; Lacaze, X.; Du Cheyron, P.; Daucourt, J.; Lorgeou, J.; Piraux, F.; Guerreiro, L.; Pérez, P.; Calus, M.; et al. A Reaction Norm Model for Genomic Selection Using High-Dimensional Genomic and Environmental Data. Theor. Appl. Genet. 2014, 127, 595–607. [Google Scholar] [CrossRef] [Green Version]
  36. Costa-Neto, G.; Fritsche-Neto, R.; Crossa, J. Nonlinear Kernels, Dominance, and Envirotyping Data Increase the Accuracy of Genome-Based Prediction in Multi-Environment Trials. Heredity 2021, 126, 92–106. [Google Scholar] [CrossRef] [PubMed]
  37. Chen, Z.Q.; Karlsson, B.; Wu, H.X. Patterns of Additive Genotype-by-Environment Interaction in Tree Height of Norway Spruce in Southern and Central Sweden. Tree Genet. Genomes 2017, 13, 25. [Google Scholar] [CrossRef] [Green Version]
  38. Korecký, J.; Čepl, J.; Stejskal, J.; Faltinová, Z.; Dvořák, J.; Lstibůrek, M.; El-Kassaby, Y.A. Genetic Diversity of Norway Spruce Ecotypes Assessed by GBS-Derived SNPs. Sci. Rep. 2021, 11, 23119. [Google Scholar] [CrossRef]
  39. Hannrup, B.; Cahalan, C.; Chantre, G.; Grabner, M.; Karlsson, B.; Le Bayon, I.; Jones, G.L.; Müller, U.; Pereira, H.; Rodrigues, J.C.; et al. Genetic Parameters of Growth and Wood Quality Traits in Picea abies. Scand. J. For. Res. 2004, 19, 14–29. [Google Scholar] [CrossRef]
  40. Androsiuk, P.; Shimono, A.; Westin, J.; Lindgren, D.; Fries, A.; Wang, X.R. Genetic Status of Norway Spruce (Picea abies) Breeding Populations for Northern Sweden. Silvae Genet. 2013, 62, 127–136. [Google Scholar] [CrossRef] [Green Version]
  41. Zianis, D.; Muukkonen, P.; Mäkipää, R.; Mencuccini, M. Biomass and Stem Volume Equations for Tree Species in Europe; Tammer-Paino Oy: Tampere, Finland, 2005; Volume 4, ISBN 9514019830. [Google Scholar]
  42. Spinoni, J.; Naumann, G.; Vogt, J.V.; Barbosa, P. The Biggest Drought Events in Europe from 1950 to 2012. J. Hydrol. Reg. Stud. 2015, 3, 509–524. [Google Scholar] [CrossRef]
  43. Washburn, J.D.; Burch, M.B.; Franco, J.A.V. Predictive Breeding for Maize: Making Use of Molecular Phenotypes, Machine Learning, and Physiological Crop Models. Crop Sci. 2020, 60, 622–638. [Google Scholar] [CrossRef]
  44. Messina, C.D.; Technow, F.; Tang, T.; Totir, R.; Gho, C.; Cooper, M. Leveraging Biological Insight and Environmental Variation to Improve Phenotypic Prediction: Integrating Crop Growth Models (CGM) with Whole Genome Prediction (WGP). Eur. J. Agron. 2018, 100, 151–162. [Google Scholar] [CrossRef]
  45. Berg, J.J.; Coop, G. A Population Genetic Signal of Polygenic Adaptation. PLoS Genet. 2014, 10, e1004412. [Google Scholar] [CrossRef] [Green Version]
  46. Alberto, F.J.; Aitken, S.N.; Alía, R.; González-Martínez, S.C.; Hänninen, H.; Kremer, A.; Lefèvre, F.; Lenormand, T.; Yeaman, S.; Whetten, R.; et al. Potential for Evolutionary Responses to Climate Change - Evidence from Tree Populations. Glob. Chang. Biol. 2013, 19, 1645–1661. [Google Scholar] [CrossRef] [Green Version]
  47. Milesi, P.; Berlin, M.; Chen, J.; Orsucci, M.; Li, L.; Jansson, G.; Karlsson, B.; Lascoux, M. Assessing the Potential for Assisted Gene Flow Using Past Introduction of Norway Spruce in Southern Sweden: Local Adaptation and Genetic Basis of Quantitative Traits in Trees. Evol. Appl. 2019, 12, 1946–1959. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  48. Chen, Z.Q.; Zan, Y.; Milesi, P.; Zhou, L.; Chen, J.; Li, L.; Cui, B.B.; Niu, S.; Westin, J.; Karlsson, B.; et al. Leveraging Breeding Programs and Genomic Data in Norway Spruce (Picea abies L. Karst) for GWAS Analysis. Genome Biol. 2021, 22, 179. [Google Scholar] [CrossRef] [PubMed]
Figure 1. (a) Map of trial locations across Austria, labeled based on a combination of altitude, longitude, and latitude. (b) Plot of the first and second principal components (PCs) from principal component analysis (PCA) of trial environments (location-year) based on geospatial and climatic data, with trial environments labeled based on a combination of altitude, longitude, and latitude. (c) Plot of the first and second PCs from PCA of trial environments based on geospatial and climatic data, with trial environments labeled by year, and PCA eigenvectors for the geospatial and weather predictors plotted as arrows.
Figure 1. (a) Map of trial locations across Austria, labeled based on a combination of altitude, longitude, and latitude. (b) Plot of the first and second principal components (PCs) from principal component analysis (PCA) of trial environments (location-year) based on geospatial and climatic data, with trial environments labeled based on a combination of altitude, longitude, and latitude. (c) Plot of the first and second PCs from PCA of trial environments based on geospatial and climatic data, with trial environments labeled by year, and PCA eigenvectors for the geospatial and weather predictors plotted as arrows.
Forests 13 00809 g001
Figure 2. (a) Plot of the first and second principal components (PCs) from principal component analysis (PCA) of population provenances based on geospatial and climatic data. (b) Map of population provenances. Points in (a,b) were labeled based on the five clusters identified by hierarchical clustering. (c) Population provenance climatic similarity matrix, with off-diagonal cells representing the similarity between pairs of populations. Pairwise similarities are colored on a heatmap. The matrix was organized by hierarchical clustering.
Figure 2. (a) Plot of the first and second principal components (PCs) from principal component analysis (PCA) of population provenances based on geospatial and climatic data. (b) Map of population provenances. Points in (a,b) were labeled based on the five clusters identified by hierarchical clustering. (c) Population provenance climatic similarity matrix, with off-diagonal cells representing the similarity between pairs of populations. Pairwise similarities are colored on a heatmap. The matrix was organized by hierarchical clustering.
Forests 13 00809 g002
Figure 3. Histograms of tree height estimates within and across years (1980–1983, 1988), with coefficients of variation (CV) displayed above each histogram.
Figure 3. Histograms of tree height estimates within and across years (1980–1983, 1988), with coefficients of variation (CV) displayed above each histogram.
Forests 13 00809 g003
Table 1. Broad-sense heritability (H2) and number of levels and estimated variances for population/provenance, environment, and residual random effects from tree height models across (1980–1988) and within (1980–1983) years.
Table 1. Broad-sense heritability (H2) and number of levels and estimated variances for population/provenance, environment, and residual random effects from tree height models across (1980–1988) and within (1980–1983) years.
YearModelTermEstimated Variance
± Standard Error
Proportion of Total VarianceLevelsH2
AcrossNaivePopulation35.67 ± 3.670.0942780.89
(1980– Environment284.38 ± 39.620.747107
1988) Pop × Env32.55 ± 1.420.0862702
Error27.91 ± 0.630.073
ClimaticProvenance289.49 ± 44.300.0502780.96
Environment5325.40 ± 717.780.920107
Prov × Env115.60 ± 4.150.0202702
Error59.77 ± 1.290.010
1980NaivePopulation14.06 ± 1.630.3332780.81
Environment19.33 ± 5.850.45724
Pop × Env5.78 ± 0.590.137612
Error3.09 ± 0.140.073
ClimaticProvenance60.57 ± 9.810.3622780.94
Environment96.66 ± 29.470.57824
Prov × Env6.84 ± 0.600.041612
Error3.08 ± 0.140.018
1981NaivePopulation14.91 ± 1.990.1692780.72
Environment56.16 ± 16.800.63824
Pop × Env8.50 ± 1.010.097612
Error8.46 ± 0.390.096
ClimaticProvenance60.17 ± 11.200.1602780.91
Environment297.05 ± 89.450.79224
Prov × Env9.57 ± 0.960.026612
Error8.43 ± 0.390.022
1982NaivePopulation18.42 ± 3.230.0732780.57
Environment191.41 ± 56.970.76224
Pop × Env20.31 ± 2.420.081612
Error20.97 ± 0.960.083
ClimaticProvenance73.21 ± 16.980.0692780.84
Environment938.58 ± 281.110.89024
Prov × Env21.55 ± 2.210.020612
Error20.93 ± 0.970.020
1983NaivePopulation29.73 ± 6.240.0422780.49
Environment582.74 ± 173.030.82524
Pop × Env47.06 ± 5.490.067612
Error46.68 ± 2.170.066
ClimaticProvenance109.13 ± 30.340.0412780.77
Environment2427.60 ± 724.660.92224
Prov × Env50.90 ± 5.070.019612
Error46.56 ± 2.160.018
Table 2. Predictive ability and proportion of total variance in population, environment, and/or population × environment best linear unbiased predictors (BLUPs) for tree height explained by population provenance and/or trial environment climatic similarity matrices.
Table 2. Predictive ability and proportion of total variance in population, environment, and/or population × environment best linear unbiased predictors (BLUPs) for tree height explained by population provenance and/or trial environment climatic similarity matrices.
Response
(BLUPs)
ModelProvenance
(prop. var.)
Environment
(prop. var.)
Predictive Ability
Population(8)0.450 ± 0.068 0.441 ± 0.083
Environment(9) 0.706 ± 0.1150.476 ± 0.157
Population × Environment(10)7 × 10−7 ± 5 × 10−7 −0.068 ± 0.061 b
(11)1 × 10−6 ± 6 × 10−7 −0.169 ± 0.037 a
(12)0.014 ± 0.0495 × 10−17 ± 3 × 10−17−0.090 ± 0.073 b
Population + (13)0.911 ± 0.091 0.678 ± 0.064 a
Population × Environment(14) 0.024 ± 0.0900.063 ± 0.084 b
(15)0.921 ± 0.0096 × 10−18 ± 6 × 10−170.676 ± 0.028 a
All models were five-fold cross-validated with 20 replications. Predictive ability was defined as the Pearson’s correlation between the observed and predicted values of the test set. Groups that are not connected by the same letter have significantly different predictive ability for population × environment BLUPs or population + population × environment BLUPs (Tukey’s HSD test, p < 0.05).
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Morales, L.; Swarts, K. Heritable and Climatic Sources of Variation in Juvenile Tree Growth in an Austrian Common Garden Experiment of Central European Norway Spruce Populations. Forests 2022, 13, 809. https://doi.org/10.3390/f13050809

AMA Style

Morales L, Swarts K. Heritable and Climatic Sources of Variation in Juvenile Tree Growth in an Austrian Common Garden Experiment of Central European Norway Spruce Populations. Forests. 2022; 13(5):809. https://doi.org/10.3390/f13050809

Chicago/Turabian Style

Morales, Laura, and Kelly Swarts. 2022. "Heritable and Climatic Sources of Variation in Juvenile Tree Growth in an Austrian Common Garden Experiment of Central European Norway Spruce Populations" Forests 13, no. 5: 809. https://doi.org/10.3390/f13050809

APA Style

Morales, L., & Swarts, K. (2022). Heritable and Climatic Sources of Variation in Juvenile Tree Growth in an Austrian Common Garden Experiment of Central European Norway Spruce Populations. Forests, 13(5), 809. https://doi.org/10.3390/f13050809

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