Calculating the genetic contribution via heritable variation
We estimated components of variation across all the data from both
gardens to calculate broad-sense heritability and Qst.
The spatial autocorrelation corrected data were used to parse variation
in our nested hierarchical structure of population, genotype, and
environment via the following equations:
\begin{equation}
{\left(1\right)Y}_{\text{Gigp}}=\ {{\ ss}_{G}+\alpha_{G}+\ \alpha_{\text{Gp}}+\ \alpha}_{\text{Gpg}}+\varepsilon_{\text{Gipg}}\nonumber \\
\end{equation}Our goal is to quantify within-garden genetic variation, we therefore
parsed variation in branches grown in two different gardens (G) using
the fixed effect \({\ ss}_{G}\). The parameters p for population
(i.e. provenance of genotype), g for genotype, and i for theith individual tree sampled are all random
effects. All branch data were modeled as gamma distributions using
Bayesian mixed regression models via the package rstan v.2.18.2
in (T. S. D. Team, 2014). The random effects outcomes (\(\alpha^{\prime}s\)) and
fixed effect (\({\ ss}_{G})\ \)of equation (1) were estimated as the
mean of 6,000 random draws from the posterior distribution (Table 1 &
S3).
The resultant fixed and random effect estimates
(\({{\ ss}_{G}+\alpha_{G}+\ \alpha_{\text{Gp}}+\ \alpha}_{\text{Ggp}})\)from equation 1 were then used as our genetic estimate for each genotype
and are displayed in all graphical analyses (often referred to as a Best
Linear Unbiased Prediction or BLUP). The variation parameters estimated
from equation 1 were used to calculate broad-sense heritability,
H2, and Qst of each of the traits.
H2 was calculated for all traits in stems, branches,
and roots, using the random effects variances from equation (1) as:
(2)\((\sigma_{\text{Genotype}}^{2}+\sigma_{\text{Population}}^{2})\ /\ (\sigma_{\text{Population}}^{2}+\sigma_{\text{Genotype}}^{2}+\ \sigma_{\text{Microenvironment}}^{2})\)
Population variance was included with Genotype as it is also
representative of genetic differences between individuals. Genotype
variance was taken as the variance among replicates and
microenvironmental variance was taken as the residual variance of the
model. Qst was calculated via the formula (Spitze, 1993;
Whitlock & Gilbert, 2012):
(3)\(\sigma_{\text{Population}}^{2}\ /\ (2\sigma_{\text{Genotype}}^{2}+\ \sigma_{\text{Population}}^{2})\)
Fst was taken from a previous publication using the same
genotypes and calculated in 1-kb windows as
(πT-πS)/πT; where
πT is SNP diversity across all individuals and
πS is weighted within-population SNP diversity (Evans et
al., 2014).