Modeling additive and non-additive effects in a hybrid population using genome-wide genotyping: prediction accuracy implications

Hybrids are broadly used in plant breeding and accurate estimation of variance components is crucial for optimizing genetic gain. Genome-wide information may be used to explore models designed to assess the extent of additive and non-additive variance and test their prediction accuracy for the genomic selection. Ten linear mixed models, involving pedigree- and marker-based relationship matrices among parents, were developed to estimate additive (A), dominance (D) and epistatic (AA, AD and DD) effects. Five complementary models, involving the gametic phase to estimate marker-based relationships among hybrid progenies, were developed to assess the same effects. The models were compared using tree height and 3303 single-nucleotide polymorphism markers from 1130 cloned individuals obtained via controlled crosses of 13 Eucalyptus urophylla females with 9 Eucalyptus grandis males. Akaike information criterion (AIC), variance ratios, asymptotic correlation matrices of estimates, goodness-of-fit, prediction accuracy and mean square error (MSE) were used for the comparisons. The variance components and variance ratios differed according to the model. Models with a parent marker-based relationship matrix performed better than those that were pedigree-based, that is, an absence of singularities, lower AIC, higher goodness-of-fit and accuracy and smaller MSE. However, AD and DD variances were estimated with high s.es. Using the same criteria, progeny gametic phase-based models performed better in fitting the observations and predicting genetic values. However, DD variance could not be separated from the dominance variance and null estimates were obtained for AA and AD effects. This study highlighted the advantages of progeny models using genome-wide information.     

On the Additive and Dominant Variance and Covariance of Individuals Within the Genomic Selection Scope

Genomic evaluation models can fit additive and dominant SNP effects. Under quantitative genetics theory, additive or “breeding” values of individuals are generated by substitution effects, which involve both “biological” additive and dominant effects of the markers. Dominance deviations include only a portion of the biological dominant effects of the markers. Additive variance includes variation due to the additive and dominant effects of the markers. We describe a matrix of dominant genomic relationships across individuals, D, which is similar to the G matrix used in genomic best linear unbiased prediction. This matrix can be used in a mixed-model context for genomic evaluations or to estimate dominant and additive variances in the population. From the “genotypic” value of individuals, an alternative parameterization defines additive and dominance as the parts attributable to the additive and dominant effect of the markers. This approach underestimates the additive genetic variance and overestimates the dominance variance. Transforming the variances from one model into the other is trivial if the distribution of allelic frequencies is known. We illustrate these results with mouse data (four traits, 1884 mice, and 10,946 markers) and simulated data (2100 individuals and 10,000 markers). Variance components were estimated correctly in the model, considering breeding values and dominance deviations. For the model considering genotypic values, the inclusion of dominant effects biased the estimate of additive variance. Genomic models were more accurate for the estimation of variance components than their pedigree-based counterparts.     

Estimating Additive and Non-Additive Genetic Variances and Predicting Genetic Merits Using Genome-Wide Dense Single Nucleotide Polymorphism Markers

Non-additive genetic variation is usually ignored when genome-wide markers are used to study the genetic architecture and genomic prediction of complex traits in human, wild life, model organisms or farm animals. However, non-additive genetic effects may have an important contribution to total genetic variation of complex traits. This study presented a genomic BLUP model including additive and non-additive genetic effects, in which additive and non-additive genetic relation matrices were constructed from information of genome-wide dense single nucleotide polymorphism (SNP) markers. In addition, this study for the first time proposed a method to construct dominance relationship matrix using SNP markers and demonstrated it in detail. The proposed model was implemented to investigate the amounts of additive genetic, dominance and epistatic variations, and assessed the accuracy and unbiasedness of genomic predictions for daily gain in pigs. In the analysis of daily gain, four linear models were used: 1) a simple additive genetic model (MA), 2) a model including both additive and additive by additive epistatic genetic effects (MAE), 3) a model including both additive and dominance genetic effects (MAD), and 4) a full model including all three genetic components (MAED). Estimates of narrow-sense heritability were 0.397, 0.373, 0.379 and 0.357 for models MA, MAE, MAD and MAED, respectively. Estimated dominance variance and additive by additive epistatic variance accounted for 5.6% and 9.5% of the total phenotypic variance, respectively. Based on model MAED, the estimate of broad-sense heritability was 0.506. Reliabilities of genomic predicted breeding values for the animals without performance records were 28.5%, 28.8%, 29.2% and 29.5% for models MA, MAE, MAD and MAED, respectively. In addition, models including non-additive genetic effects improved unbiasedness of genomic predictions.     

WILL POPULATION BOTTLENECKS AND MULTILOCUS EPISTASIS INCREASE ADDITIVE GENETIC VARIANCE?

We apply new analytical methods to understand the consequences of population bottlenecks for expected additive genetic variance. We analyze essentially all models for multilocus epistasis that have been numerically simulated to demonstrate increased additive variance. We conclude that for biologically plausible models, large increases in expected additive variance--attributable to epistasis rather than dominance--are unlikely. Naciri-Graven and Goudet (2003) found that as the number of epistatically interacting loci increases, additive variance tends to be inflated more after a bottleneck. We argue that this result reflects biologically unrealistic aspects of their models. Specifically, as the number of loci increases, higher-order epistatic interactions become increasingly important in these models, with an increasing fraction of the genetic variance becoming nonadditive, contrary to empirical observations. As shown by Barton and Turelli (2004), without dominance, conversion of nonadditive to additive variance depends only on the variance components and not on the number of loci per se. Numerical results indicating that more inbreeding is needed to produce maximal release of additive variance with more loci follow directly from our analytical results, which show that high levels of inbreeding (F > 0.5) are needed for significant conversion of higher-order components. We discuss alternative approaches to modeling multilocus epistasis and understanding its consequences.     

Genomic prediction of hybrid crops allows disentangling dominance and epistasis

We revisited, in a genomic context, the theory of hybrid genetic evaluation models of hybrid crosses of pure lines, as the current practice is largely based on infinitesimal model assumptions. Expressions for covariances between hybrids due to additive substitution effects and dominance and epistatic deviations were analytically derived. Using dense markers in a GBLUP analysis, it is possible to split specific combining ability into dominance and across-groups epistatic deviations, and to split general combining ability (GCA) into within-line additive effects and within-line additive by additive (and higher order) epistatic deviations. We analyzed a publicly available maize data set of Dent × Flint hybrids using our new model (called GCA-model) up to additive by additive epistasis. To model higher order interactions within GCAs, we also fitted “residual genetic” line effects. Our new GCA-model was compared with another genomic model which assumes a uniquely defined effect of genes across origins. Most variation in hybrids is accounted by GCA. Variances due to dominance and epistasis have similar magnitudes. Models based on defining effects either differently or identically across heterotic groups resulted in similar predictive abilities for hybrids. The currently used model inflates the estimated additive genetic variance. This is not important for hybrid predictions but has consequences for the breeding scheme—e.g. overestimation of the genetic gain within heterotic group. Therefore, we recommend using GCA-model, which is appropriate for genomic prediction and variance component estimation in hybrid crops using genomic data, and whose results can be practically interpreted and used for breeding purposes.     

Estimates of genetic variance in an F2 maize population

Maize (Zea mays L.) breeders have used several genetic-statistical models to study the inheritance of quantitative traits. These models provide information on the importance of additive, dominance, and epistatic genetic variance for a quantitative trait. Estimates of genetic variances are useful in understanding heterosis and determining the response to selection. The objectives of this study were to estimate additive and dominance genetic variances and the average level of dominance for an F2 population derived from the B73 x Mo17 hybrid and use weighted least squares to determine the importance of digenic epistatic variances relative to additive and dominance variances. Genetic variances were estimated using Design III and weighted least squares analyses. Both analyses determined that dominance variance was more important than additive variance for grain yield. For other traits, additive genetic variance was more important than dominance variance. The average level of dominance suggests either overdominant gene effects were present for grain yield or pseudo-overdominance because of linkage disequilibrium in the F2 population. Epistatic variances generally were not significantly different from zero and therefore were relatively less important than additive and dominance variances. For several traits estimates of additive by additive epistatic variance decreased estimates of additive genetic variance, but generally the decrease in additive genetic variance was not significant.     

Mixed Model Methods for Genomic Prediction and Variance Component Estimation of Additive and Dominance Effects Using SNP Markers

We established a genomic model of quantitative trait with genomic additive and dominance relationships that parallels the traditional quantitative genetics model, which partitions a genotypic value as breeding value plus dominance deviation and calculates additive and dominance relationships using pedigree information. Based on this genomic model, two sets of computationally complementary but mathematically identical mixed model methods were developed for genomic best linear unbiased prediction (GBLUP) and genomic restricted maximum likelihood estimation (GREML) of additive and dominance effects using SNP markers. These two sets are referred to as the CE and QM sets, where the CE set was designed for large numbers of markers and the QM set was designed for large numbers of individuals. GBLUP and associated accuracy formulations for individuals in training and validation data sets were derived for breeding values, dominance deviations and genotypic values. Simulation study showed that GREML and GBLUP generally were able to capture small additive and dominance effects that each accounted for 0.00005–0.0003 of the phenotypic variance and GREML was able to differentiate true additive and dominance heritability levels. GBLUP of the total genetic value as the summation of additive and dominance effects had higher prediction accuracy than either additive or dominance GBLUP, causal variants had the highest accuracy of GREML and GBLUP, and predicted accuracies were in agreement with observed accuracies. Genomic additive and dominance relationship matrices using SNP markers were consistent with theoretical expectations. The GREML and GBLUP methods can be an effective tool for assessing the type and magnitude of genetic effects affecting a phenotype and for predicting the total genetic value at the whole genome level.     

Genetics of tassel branch number in maize and its implications for a selection program for small tassel size

Summary Tassel branch numbers of six crosses of maize (Zea mays L.) were analyzed to determine inheritance of this trait. Generation mean analyses were used to estimate genetic effects, and additive and nonadditive components of variance were calculated and evaluated for bias due to linkage. Both narrow-sense and broad-sense heritabilities were estimated. Additive genetic variance estimates were significant in five of the six crosses, whereas estimates of variance due to nonadditive components were significant in only three crosses. Additionally, estimates of additive variance components usually were larger than corresponding nonadditive components. There was no evidence for linkage bias in these estimates. Estimates of additive genetic effects were significant in four of six crosses, but significant dominance, additive × additive and additive × dominance effects also were detected. Additive, dominance, and epistatic gene action, therefore, all influenced the inheritance of tassel branch number, but additive gene action was most important. Both narrow-sense and broadsense heritability estimates were larger than those reported for other physiological traits of maize and corroborated conclusions concerning the importance of additive gene action inferred from analyses of genetic effects and variances. We concluded that selection for smalltasseled inbreds could be accomplished most easily through a mass-selection and/or pedigree-selection system. Production of a small-tasseled hybrid would require crossing of two small-tasseled inbreds. We proposed two genetic models to explain unexpected results obtained for two crosses. One model involved five interacting loci and the other employed two loci displaying only additive and additive × additive gene action.Journal Paper No. J-9231 of the Iowa Agriculture and Home Economics Experiment Station, Ames, Iowa 50011. Project No. 2152     

Robust modeling of additive and nonadditive variation with intuitive inclusion of expert knowledge

We propose a novel Bayesian approach that robustifies genomic modeling by leveraging expert knowledge (EK) through prior distributions. The central component is the hierarchical decomposition of phenotypic variation into additive and nonadditive genetic variation, which leads to an intuitive model parameterization that can be visualized as a tree. The edges of the tree represent ratios of variances, for example broad-sense heritability, which are quantities for which EK is natural to exist. Penalized complexity priors are defined for all edges of the tree in a bottom-up procedure that respects the model structure and incorporates EK through all levels. We investigate models with different sources of variation and compare the performance of different priors implementing varying amounts of EK in the context of plant breeding. A simulation study shows that the proposed priors implementing EK improve the robustness of genomic modeling and the selection of the genetically best individuals in a breeding program. We observe this improvement in both variety selection on genetic values and parent selection on additive values; the variety selection benefited the most. In a real case study, EK increases phenotype prediction accuracy for cases in which the standard maximum likelihood approach did not find optimal estimates for the variance components. Finally, we discuss the importance of EK priors for genomic modeling and breeding, and point to future research areas of easy-to-use and parsimonious priors in genomic modeling.     

Genomic analysis of dominance effects on milk production and conformation traits in Fleckvieh cattle

Background

Estimates of dominance variance in dairy cattle based on pedigree data vary considerably across traits and amount to up to 50% of the total genetic variance for conformation traits and up to 43% for milk production traits. Using bovine SNP (single nucleotide polymorphism) genotypes, dominance variance can be estimated both at the marker level and at the animal level using genomic dominance effect relationship matrices. Yield deviations of high-density genotyped Fleckvieh cows were used to assess cross-validation accuracy of genomic predictions with additive and dominance models. The potential use of dominance variance in planned matings was also investigated.

Results

Variance components of nine milk production and conformation traits were estimated with additive and dominance models using yield deviations of 1996 Fleckvieh cows and ranged from 3.3% to 50.5% of the total genetic variance. REML and Gibbs sampling estimates showed good concordance. Although standard errors of estimates of dominance variance were rather large, estimates of dominance variance for milk, fat and protein yields, somatic cell score and milkability were significantly different from 0. Cross-validation accuracy of predicted breeding values was higher with genomic models than with the pedigree model. Inclusion of dominance effects did not increase the accuracy of the predicted breeding and total genetic values. Additive and dominance SNP effects for milk yield and protein yield were estimated with a BLUP (best linear unbiased prediction) model and used to calculate expectations of breeding values and total genetic values for putative offspring. Selection on total genetic value instead of breeding value would result in a larger expected total genetic superiority in progeny, i.e. 14.8% for milk yield and 27.8% for protein yield and reduce the expected additive genetic gain only by 4.5% for milk yield and 2.6% for protein yield.

Conclusions

Estimated dominance variance was substantial for most of the analyzed traits. Due to small dominance effect relationships between cows, predictions of individual dominance deviations were very inaccurate and including dominance in the model did not improve prediction accuracy in the cross-validation study. Exploitation of dominance variance in assortative matings was promising and did not appear to severely compromise additive genetic gain.     

The correlation of substitution effects across populations and generations in the presence of nonadditive functional gene action

Allele substitution effects at quantitative trait loci (QTL) are part of the basis of quantitative genetics theory and applications such as association analysis and genomic prediction. In the presence of nonadditive functional gene action, substitution effects are not constant across populations. We develop an original approach to model the difference in substitution effects across populations as a first order Taylor series expansion from a “focal” population. This expansion involves the difference in allele frequencies and second-order statistical effects (additive by additive and dominance). The change in allele frequencies is a function of relationships (or genetic distances) across populations. As a result, it is possible to estimate the correlation of substitution effects across two populations using three elements: magnitudes of additive, dominance, and additive by additive variances; relationships (Nei’s minimum distances or Fst indexes); and assumed heterozygosities. Similarly, the theory applies as well to distinct generations in a population, in which case the distance across generations is a function of increase of inbreeding. Simulation results confirmed our derivations. Slight biases were observed, depending on the nonadditive mechanism and the reference allele. Our derivations are useful to understand and forecast the possibility of prediction across populations and the similarity of GWAS effects.     

THE EVOLUTIONARY GENETICS OF MALE LIFE-HISTORY CHARACTERS IN DROSOPHILA MELANOGASTER

Alternative models of the maintenance of genetic variability, theories of life-history evolution, and theories of sexual selection and mate choice can be tested by measuring additive and nonadditive genetic variances of components of fitness. A quantitative genetic breeding design was used to produce estimates of genetic variances for male life-history traits in Drosophila melanogaster. Additive genetic covariances and correlations between traits were also estimated. Flies from a large, outbred, laboratory population were assayed for age-specific competitive mating ability, age-specific survivorship, body mass, and fertility. Variance-component analysis then allowed the decomposition of phenotypic variation into components associated with additive genetic, nonadditive genetic, and environmental variability. A comparison of dominance and additive components of genetic variation provides little support for an important role for balancing selection in maintaining genetic variance in this suite of traits. The results provide support for the mutation-accumulation theory, but not the antagonistic-pleiotropy theory of senescence. No evidence is found for the positive genetic correlations between mating success and offspring quality or quantity that are predicted by “good genes” models of sexual selection. Additive genetic coefficients of variation for life-history characters are larger than those for body weight. Finally, this set of male life-history characters exhibits a very low correspondence between estimates of genetic and phenotypic correlations.     

Single-step genomic prediction of Eucalyptus dunnii using different identity-by-descent and identity-by-state relationship matrices

Genomic selection based on the single-step genomic best linear unbiased prediction (ssGBLUP) approach is becoming an important tool in forest tree breeding. The quality of the variance components and the predictive ability of the estimated breeding values (GEBV) depends on how well marker-based genomic relationships describe the actual genetic relationships at unobserved causal loci. We investigated the performance of GEBV obtained when fitting models with genomic covariance matrices based on two identity-by-descent (IBD) and two identity-by-state (IBS) relationship measures. Multiple-trait multiple-site ssGBLUP models were fitted to diameter and stem straightness in five open-pollinated progeny trials of Eucalyptus dunnii, genotyped using the EUChip60K. We also fitted the conventional ABLUP model with a pedigree-based covariance matrix. Estimated relationships from the IBD estimators displayed consistently lower standard deviations than those from the IBS approaches. Although ssGBLUP based in IBS estimators resulted in higher trait-site heritabilities, the gain in accuracy of the relationships using IBD estimators has resulted in higher predictive ability and lower bias of GEBV, especially for low-heritability trait-site. ssGBLUP based on IBS and IBD approaches performed considerably better than the traditional ABLUP. In summary, our results advocate the use of the ssGBLUP approach jointly with the IBD relationship matrix in open-pollinated forest tree evaluation.Subject terms: Plant breeding, Genomics     

Genetic Models for the Analysis of Data from the Families of Identical Twins

Genetic models are described which exploit the unique relationships that exist within the families of identical twins to obtain weighted least squares estimates of additive, dominance and epistatic components of genetic variance as well as estimates of the contributions of X-linked genes, maternal effects and three sources of environmental variation. Since all of the relationships required to achieve a resolution of these variance components are contained within each family unit, the model would appear to be superior to previous approaches to the analysis of quantitative traits in man.     

Something Old and Something New: Wedding Recombinant Inbred Lines with Traditional Line Cross Analysis Increases Power to Describe Gene Interactions

In this paper we present a novel approach to quantifying genetic architecture that combines recombinant inbred lines (RIL) with line cross analysis (LCA). LCA is a method of quantifying directional genetic effects (i.e. summed effects of all loci) that differentiate two parental lines. Directional genetic effects are thought to be critical components of genetic architecture for the long term response to selection and as a cause of inbreeding depression. LCA typically begins with two inbred parental lines that are crossed to produce several generations such as F1, F2, and backcrosses to each parent. When a RIL population (founded from the same P1 and P2 as was used to found the line cross population) is added to the LCA, the sampling variance of several nonadditive genetic effect estimates is greatly reduced. Specifically, estimates of directional dominance, additive x additive, and dominance x dominance epistatic effects are reduced by 92%, 94%, and 56% respectively. The RIL population can be simultaneously used for QTL identification, thus uncovering the effects of specific loci or genomic regions as elements of genetic architecture. LCA and QTL mapping with RIL provide two qualitatively different measures of genetic architecture with the potential to overcome weaknesses of each approach alone. This approach provides cross-validation of the estimates of additive and additive x additive effects, much smaller confidence intervals on dominance, additive x additive and dominance x dominance estimates, qualitatively different measures of genetic architecture, and the potential when used together to balance the weaknesses of LCA or RIL QTL analyses when used alone.     

Population bottlenecks increase additive genetic variance but do not break a selection limit in rain forest Drosophila

According to neutral quantitative genetic theory, population bottlenecks are expected to decrease standing levels of additive genetic variance of quantitative traits. However, some empirical and theoretical results suggest that, if nonadditive genetic effects influence the trait, bottlenecks may actually increase additive genetic variance. This has been an important issue in conservation genetics where it has been suggested that small population size might actually experience an increase rather than a decrease in the rate of adaptation. Here we test if bottlenecks can break a selection limit for desiccation resistance in the rain forest-restricted fly Drosophila bunnanda. After one generation of single-pair mating, additive genetic variance for desiccation resistance increased to a significant level, on average higher than for the control lines. Line crosses revealed that both dominance and epistatic effects were responsible for the divergence in desiccation resistance between the original control and a bottlenecked line exhibiting increased additive genetic variance for desiccation resistance. However, when bottlenecked lines were selected for increased desiccation resistance, there was only a small shift in resistance, much less than predicted by the released additive genetic variance. The small selection response in the bottlenecked lines was no greater than that observed in the control lines. Thus bottlenecks might produce a statistically detectable change in additive genetic variance but this change has no impact on the response to selection.     

How does selfing affect the genetic variance of quantitative traits? An updated meta‐analysis on empirical results in angiosperm species

Most theoretical works predict that selfing should reduce the level of additive genetic variance available for quantitative traits within natural populations. Despite a growing number of quantitative genetic studies undertaken during the last two decades, this prediction is still not well supported empirically. To resolve this issue and confirm or reject theoretical predictions, we reviewed quantitative trait heritability estimates from natural plant populations with different rates of self‐fertilization and carried out a meta‐analysis. In accordance with models of polygenic traits under stabilizing selection, we found that the fraction of additive genetic variance is negatively correlated with the selfing rate. Although the mating system explains a moderate fraction of the variance, the mean reduction of narrow‐sense heritability values between strictly allogamous and predominantly selfing populations is strong, around 60%. Because some nonadditive components of genetic variance become selectable under inbreeding, we determine whether self‐fertilization affects the relative contribution of these components to genetic variance by comparing narrow‐sense heritability estimates from outcrossing populations with broad‐sense heritability estimated in autogamous populations. Results suggest that these nonadditive components of variance may restore some genetic variance in predominantly selfing populations; it remains, however, uncertain how these nonadditive components will contribute to adaptation.     

Multienvironment genomic variance decomposition analysis of open-pollinated Interior spruce (<Emphasis Type="Italic">Picea glauca</Emphasis> x <Emphasis Type="Italic">engelmannii</Emphasis>)

The advantages of open-pollinated (OP) family testing over controlled crossing (i.e., structured pedigree) are the potential to screen and rank a large number of parents and offspring with minimal cost and efforts; however, the method produces inflated genetic parameters as the actual sibling relatedness within OP families rarely meets the half-sib relatedness assumption. Here, we demonstrate the unsurpassed utility of OP testing after shifting the analytical mode from pedigree- (ABLUP) to genomic-based (GBLUP) relationship using phenotypic tree height (HT) and wood density (WD) and genotypic (30k SNPs) data for 1126 38-year-old Interior spruce (Picea glauca (Moench) Voss x P. engelmannii Parry ex Engelm.) trees, representing 25 OP families, growing on three sites in Interior British Columbia, Canada. The use of the genomic realized relationship permitted genetic variance decomposition to additive, dominance, and epistatic genetic variances, and their interactions with the environment, producing more accurate narrow-sense heritability and breeding value estimates as compared to the pedigree-based counterpart. The impact of retaining (random folding) vs. removing (family folding) genetic similarity between the training and validation populations on the predictive accuracy of genomic selection was illustrated and highlighted the former caveats and latter advantages. Moreover, GBLUP models allowed breeding value prediction for individuals from families that were not included in the developed models, which was not possible with the ABLUP. Response to selection differences between the ABLUP and GBLUP models indicated the presence of systematic genetic gain overestimation of 35 and 63% for HT and WD, respectively, mainly caused by the inflated estimates of additive genetic variance and individuals’ breeding values given by the ABLUP models. Extending the OP genomic-based models from single to multisite made the analysis applicable to existing OP testing programs.     

Detecting epistatic genetic variance with a clonally replicated design: models for lowvs high-order nonallelic interaction

A quantitative genetic model, that uses known family structure with clonal replicates to separate genetic variance into its additive, dominance and epistatic components, is available in the current literature. Making use of offspring testing, this model is based on the theory that components of variance from the linear model of an experimental design may be expressed in terms of expected covariances among relatives. However, if interactions between a pair of quantitative trait loci (QTLs) explain a large proportion of the total epistasis, it will seriously overestimate the additive and dominance variances but underestimate the epistatic variance. In the present paper, a new model is developed to manipulate this problem by combining parental and offspring material into the same test. Under the condition described above, the new model can provide an accurate estimate for additive x additive variances. Also, its accuracy in estimating dominance and total epistatic variances is much greater than the accuracy of the previous model. However, if there is obvious evidence showing the major contribution of high-order interactions, especially among 4QTLs, to the total epistasis, the previous model is more appropriate to partition the genetic variance for a quantitative trait. The re-analysis of an example from a factorial mating design in poplar shows large differences in estimating variance components between the new and previous models when two different assumptions (lowvs high-order epistatic interactions) are used. The new model will be an alternative to estimating the mode of quantitative inheritance for species, especially for longlived, predominantly outcrossing forest trees, that can be clonally replicated.