Analysis of cytoplasmic and maternal effects. II. Genetic models for triploid endosperms

Genetic models for quantitative traits of triploid endosperms are proposed for the analysis of direct gene effects, cytoplasmic effects, and maternal gene effects. The maternal effect is partitioned into maternal additive and dominance components. In the full genetic model, the direct effect is partitioned into direct additive and dominance components and high-order dominance component, which are the cumulative effects of three-allele interactions. If the high-order dominance effects are of no importance, a reduced genetic model can be used. Monte Carlo simulations were conducted in this study for demonstrating unbiasedness of estimated variance and covariance components from the MINQUE (0/1) procedure, which is a minimum norm quadratic unbiased estimation (MINQUE) method setting 0 for all the prior covariances and 1 for all the prior variances. Robustness of estimating variance and covariance components for the genetic models was tested by simulations. Both full and reduced genetic models are shown to be robust for estimating variance and covariance components under several situations of no specific effects. Efficiency of predicting random genetic effects for the genetic models by the MINQUE (0/1) procedure was compared with the best linear unbiased prediction (BLUP). A worked example is given to illustrate the use of the reduced genetic model for kernel growth characteristics in corn (Zea mays L.).     

Analysis of cytoplasmic and maternal effects I. A genetic model for diploid plant seeds and animals

A genetic model for modified diallel crosses is proposed for estimating variance and covariance components of cytoplasmic, maternal additive and dominance effects, as well as direct additive and dominance effects. Monte Carlo simulations were conducted to compare the efficiencies of minimum norm quadratic unbiased estimation (MINQUE) methods. For both balanced and unbalanced mating designs, MINQUE (0/1), which has 0 for all the prior covariances and 1 for all the prior variances, has similar efficiency to MINQUE(), which has parameter values for the prior values. Unbiased estimates of variance and covariance components and their sampling variances could be obtained with MINQUE(0/1) and jackknifing. A t-test following jackknifing is applicable to test hypotheses for zero variance and covariance components. The genetic model is robust for estimating variance and covariance components under several situations of no specific effects. A MINQUE(0/1) procedure is suggested for unbiased estimation of covariance components between two traits with equal design matrices. Methods of unbiased prediction for random genetic effects are discussed. A linear unbiased prediction (LUP) method is shown to be efficient for the genetic model. An example is given for a demonstration of estimating variance and covariance components and predicting genetic effects.     

Consequences of paternally inherited effects on the genetic evaluation of maternal effects

Background

Mixed models are commonly used for the estimation of variance components and genetic evaluation of livestock populations. Some evaluation models include two types of additive genetic effects, direct and maternal. Estimates of variance components obtained with models that account for maternal effects have been the subject of a long-standing controversy about strong negative estimates of the covariance between direct and maternal effects. Genomic imprinting is known to be in some cases statistically confounded with maternal effects. In this study, we analysed the consequences of ignoring paternally inherited effects on the partitioning of genetic variance.

Results

We showed that the existence of paternal parent-of-origin effects can bias the estimation of variance components when maternal effects are included in the evaluation model. Specifically, we demonstrated that adding a constraint on the genetic parameters of a maternal model resulted in correlations between relatives that were the same as those obtained with a model that fits only paternally inherited effects for most pairs of individuals, as in livestock pedigrees. The main consequence is an upward bias in the estimates of the direct and maternal additive genetic variances and a downward bias in the direct-maternal genetic covariance. This was confirmed by a simulation study that investigated five scenarios, with the trait affected by (1) only additive genetic effects, (2) only paternally inherited effects, (3) additive genetic and paternally inherited effects, (4) direct and maternal additive genetic effects and (5) direct and maternal additive genetic plus paternally inherited effects. For each scenario, the existence of a paternally inherited effect not accounted for by the estimation model resulted in a partitioning of the genetic variance according to the predicted pattern. In addition, a model comparison test confirmed that direct and maternal additive models and paternally inherited models provided an equivalent fit.

Conclusions

Ignoring paternally inherited effects in the maternal models for genetic evaluation can lead to a specific pattern of bias in variance component estimates, which may account for the unexpectedly strong negative direct-maternal genetic correlations that are typically reported in the literature.

Electronic supplementary material

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

Fixed and random effects in the variation of the finger ridge count: a study of fragile-X families.

A maximum-likelihood scoring technique for analysis of pedigree data allows for the concurrent estimation of random and of fixed effects in a quantitative trait. We included both types of effects in genetic models, to study the sources of variation in finger ridge count in 54 large families affected with the fragile-X disorder. The fixed effects were represented by fragile X and sex, and the random effects by environmental and genetic variance. We found a significant effect of fragile X in the mean of the finger ridge count on the thumb (finger 1) and index finger (finger 2), which had the lowest heritability and a negligible nonadditive component of genetic variance. This was in contrast with ridge counts on fingers 3 and 4, which showed little fragile-X effect, but high heritability and a significant nonadditive component. A contrast in genetic properties for ridge counts on fingers 1 and 2, compared with these counts on the remaining three fingers, may be relevant to increased selection pressures on functions of the thumb and of the index finger in evolution of modern man. We have also demonstrated an important effect of fragile X in increasing the additive variance in covariance, especially between male pairs. These findings suggest that the effect of the fragile-X genotype in finger ridge count is additive and superimposed on the normal hereditary variations in this trait.     

The use of linear mixed models to estimate variance components from data on twin pairs by maximum likelihood.

It is shown that maximum likelihood estimation of variance components from twin data can be parameterized in the framework of linear mixed models. Standard statistical packages can be used to analyze univariate or multivariate data for simple models such as the ACE and CE models. Furthermore, specialized variance component estimation software that can handle pedigree data and user-defined covariance structures can be used to analyze multivariate data for simple and complex models, including those where dominance and/or QTL effects are fitted. The linear mixed model framework is particularly useful for analyzing multiple traits in extended (twin) families with a large number of random effects.     

The effect of genetic drift on the variance/covariance components generated by multilocus additive x additive epistatic systems

The effect of population bottlenecks on the components of the genetic variance/covariance generated by n neutral independent additive x additive loci has been studied theoretically. In its simplest version, this situation can be modelled by specifying the allele frequencies and homozygous effects at each locus, and an additional factor measuring the strength of the n-th order epistatic interaction. The variance/covariance components in an infinitely large panmictic population (ancestral components) were compared with their expected values at equilibrium over replicates randomly derived from the base population, after t bottlenecks of size N (derived components). Formulae were obtained giving the derived components (and the between-line variance) as functions of the ancestral ones (alternatively, in terms of allele frequencies and effects) and the corresponding inbreeding coefficient F(t). The n-th order derived component of the genetic variance/covariance is continuously eroded by inbreeding, but the remaining components may increase initially until a critical F(t) value is attained, which is inversely related to the order of the pertinent component, and subsequently decline to zero. These changes can be assigned to the between-line variances/covariances of gene substitution and epistatic effects induced by drift. Numerical examples indicate that: (1) the derived additive variance/covariance component will generally exceed its ancestral value unless epistasis is weak; (2) the derived epistatic variance/covariance components will generally exceed their ancestral values unless allele frequencies are extreme; (3) for systems showing equal ancestral additive and total non-additive variance/covariance components, those including a smaller number of epistatic loci may generate a larger excess in additive variance/covariance after bottlenecks than others involving a larger number of loci, provided that F(t) is low. Our results indicate that it is unlikely that the rate of evolution may be significantly accelerated after population bottlenecks, in spite of occasional increments of the derived additive variance over its ancestral value.     

Bayesian adaptive Markov chain Monte Carlo estimation of genetic parameters

Accurate and fast estimation of genetic parameters that underlie quantitative traits using mixed linear models with additive and dominance effects is of great importance in both natural and breeding populations. Here, we propose a new fast adaptive Markov chain Monte Carlo (MCMC) sampling algorithm for the estimation of genetic parameters in the linear mixed model with several random effects. In the learning phase of our algorithm, we use the hybrid Gibbs sampler to learn the covariance structure of the variance components. In the second phase of the algorithm, we use this covariance structure to formulate an effective proposal distribution for a Metropolis-Hastings algorithm, which uses a likelihood function in which the random effects have been integrated out. Compared with the hybrid Gibbs sampler, the new algorithm had better mixing properties and was approximately twice as fast to run. Our new algorithm was able to detect different modes in the posterior distribution. In addition, the posterior mode estimates from the adaptive MCMC method were close to the REML (residual maximum likelihood) estimates. Moreover, our exponential prior for inverse variance components was vague and enabled the estimated mode of the posterior variance to be practically zero, which was in agreement with the support from the likelihood (in the case of no dominance). The method performance is illustrated using simulated data sets with replicates and field data in barley.     

Inbreeding and the genetic variance in floral traits of Mimulus guttatus

The additive genetic variance, V(A), is frequently used as a measure of evolutionary potential in natural plant populations. Many plants inbreed to some extent; a notable observation given that random mating is essential to the model that predicts evolutionary change from V(A). With inbreeding, V(A) is not the only relevant component of genetic variation. Several nonadditive components emerge from the combined effects of inbreeding and genetic dominance. An important empirical question is whether these components are quantitatively significant. We use maximum likelihood estimation to extract estimates for V(A) and the nonadditive 'inbreeding components' from an experimental study of the wildflower Mimulus guttatus. The inbreeding components contribute significantly to four of five floral traits, including several measures of flower size and stigma-anther separation. These results indicate that inbreeding will substantially alter the evolutionary response to natural selection on floral characters.     

Application of the covariance function approach with an iterative two-stage algorithm to the estimation of parameters of a random regression test day model for dairy production traits

The covariance function approach with an iterative two-stage algorithm of LIU et al. (2000) was applied to estimate parameters for the Polish Black-and-White dairy population based on a sample of 338 808 test day records for milk, fat, and protein yields. A multiple trait sire model was used to estimate covariances of lactation stages. A third-order Legendre polynomial was subsequently fitted to the estimated (co)variances to derive (co)variances of random regression coefficients for both additive genetic and permanent environment effects. Daily and 305-day heritability estimates obtained are consistent with several studies which used both fixed and random regression test day models. Genetic correlations between any two days in milk (DIM) of the same lactation as well as genetic correlations between the same DIM of two lactations were within a biologically acceptable range. It was shown that the applied estimation procedure can utilise very large data sets and give plausible estimates of (co)variance components.     

Model comparison for genetic evaluation of milk yield in crossbred Holsteins in the tropics

The objective of this study was to compare models for appropriate genetic parameter estimation for milk yield (305-day) in crossbred Holsteins in the tropics, where only records from crossbred cows were available. Eleven models with different effects of contemporary group (CG) at calving (herd-year-season or herd-year-month as fixed, and herd-year-month as random), age at calving (as linear or quadratic covariates, age-class, and age-class x lactation), and dominance were considered. On-farm records from small herds (n < 50) were included or excluded to validate the parameter estimates. Average Information Restricted Maximum Likelihood (AIREML) and Best Linear Unbiased Prediction (BLUP) were used to estimate variance components and breeding values. R-square (R2) and standard error of heritability (h2) were used to determine the appropriate model. The estimates of heritability from most models ranged from 0.18 to 0.22. CG formation of herd-year-month as a random effect slightly lowered the additive genetic variance but considerably decreased the permanent environmental variance. The model with age-class x lactation gave better R2 than other age adjustments. The models including records from smallholders gave similar estimates of heritability and a lower standard error than the models excluding them. The estimate of dominance variance as a proportion of total variance was close to zero. The low ratio of dominance to additive genetic variance suggested that the inclusion of dominance effects in the model was unjustified. In conclusion, the model including the effects of herd-year-month, age-class x lactation, as well as additive genetic, permanent environmental and residual effects, was the most appropriate for genetic evaluation in crossbred Holsteins, where records from smallholders could be included.     

Random regression analyses using B-spline functions to model growth of Nellore cattle

The objective of this study was to estimate (co)variance components using random regression on B-spline functions to weight records obtained from birth to adulthood. A total of 82 064 weight records of 8145 females obtained from the data bank of the Nellore Breeding Program (PMGRN/Nellore Brazil) which started in 1987, were used. The models included direct additive and maternal genetic effects and animal and maternal permanent environmental effects as random. Contemporary group and dam age at calving (linear and quadratic effect) were included as fixed effects, and orthogonal Legendre polynomials of age (cubic regression) were considered as random covariate. The random effects were modeled using B-spline functions considering linear, quadratic and cubic polynomials for each individual segment. Residual variances were grouped in five age classes. Direct additive genetic and animal permanent environmental effects were modeled using up to seven knots (six segments). A single segment with two knots at the end points of the curve was used for the estimation of maternal genetic and maternal permanent environmental effects. A total of 15 models were studied, with the number of parameters ranging from 17 to 81. The models that used B-splines were compared with multi-trait analyses with nine weight traits and to a random regression model that used orthogonal Legendre polynomials. A model fitting quadratic B-splines, with four knots or three segments for direct additive genetic effect and animal permanent environmental effect and two knots for maternal additive genetic effect and maternal permanent environmental effect, was the most appropriate and parsimonious model to describe the covariance structure of the data. Selection for higher weight, such as at young ages, should be performed taking into account an increase in mature cow weight. Particularly, this is important in most of Nellore beef cattle production systems, where the cow herd is maintained on range conditions. There is limited modification of the growth curve of Nellore cattle with respect to the aim of selecting them for rapid growth at young ages while maintaining constant adult weight.     

On flexible finite polygenic models for multiple-trait evaluation

Finite polygenic models (FPM) might be an alternative to the infinitesimal model (TIM) for the genetic evaluation of pedigreed multiple-generation populations for multiple quantitative traits. I present a general flexible Bayesian method that includes the number of genes in the FPM as an additional random variable. Markov-chain Monte Carlo techniques such as Gibbs sampling and the reversible jump sampler are used for implementation. Sampling of genotypes of all genes in the FPM is done via the use of segregation indicators. A broad range of FPM models, some combined with TIM, are empirically tested for the estimation of variance components and the number of genes in the FPM. Four simulation scenarios were studied, including genetic models with 5 or 50 additive independent diallelic genes affecting the traits, and random selection or selection on one of the traits was performed. The results in this study were based on ten replicates per simulation scenario. In the case of random selection, uniform priors on additive gene effects led to posterior mean estimates of genetic variance that were positively correlated with the number of genes fitted in the FPM. In the case of trait selection, assuming normal priors on gene effects also led to genetic variance estimates for the selected trait that were negatively correlated with the number of genes in the FPM. This negative correlation was not observed for the unselected trait. Treating the number of genes in the FPM as random revealed a positive correlation between prior and posterior mean estimates of this number, but the prior hardly affected the posterior estimates of genetic variance. Posterior inferences about the number of genes should be considered to be indicative where trait selection seems to improve the power of distinguishing between TIM and FPM. Based on the results of this study, I suggest not replacing TIM by the FPM, but combining TIM and FPM with the number of genes treated as random, to facilitate a highly flexible and thereby robust method for variance component estimation in pedigreed populations. Further study is required to explore the full potential of these models under different genetic model assumptions.     

Complex genetic effects in quantitative trait locus identification: a computationally tractable random model for use in F(2) populations

Methodology for mapping quantitative trait loci (QTL) has focused primarily on treating the QTL as a fixed effect. These methods differ from the usual models of genetic variation that treat genetic effects as random. Computationally expensive methods that allow QTL to be treated as random have been explicitly developed for additive genetic and dominance effects. By extending these methods with a variance component method (VCM), multiple QTL can be mapped. We focused on an F(2) crossbred population derived from inbred lines and estimated effects for each individual and their corresponding marker-derived genetic covariances. We present extensions to pairwise epistatic effects, which are computationally intensive because a great many individual effects must be estimated. But by replacing individual genetic effects with average genetic effects for each marker class, genetic covariances are approximated. This substantially reduces the computational burden by reducing the dimensions of covariance matrices of genetic effects, resulting in a remarkable gain in the speed of estimating the variance components and evaluating the residual log-likelihood. Preliminary results from simulations indicate competitiveness of the reduced model with multiple-interval mapping, regression interval mapping, and VCM with individual genetic effects in its estimated QTL positions and experimental power.     

MIXED MODEL APPROACHES FOR ESTIMATING GENETIC VARIANCES AND COVARIANCES

The limitations of methods for analysis of variance(ANOVA)in estimating genetic variances are discussed. Among the three methods(maximum likelihood ML, restricted maximum likelihood REML, and minimum norm quadratic unbiased estimation MINQUE)for mixed linear models, MINQUE method is presented with formulae for estimating variance components and covariances components and for predicting genetic effects. Several genetic models, which cannot be appropriately analyzed by ANOVA methods, are introduced in forms of mixed linear models. Genetic models with independent random effects can be analyzed by MINQUE(1)method whieh is a MINQUE method with all prior values setting 1. MINQUE(1)method can give unbiased estimation for variance components and covariance components, and linear unbiased prediction (LUP) for genetic effects. There are more complicate genetic models for plant seeds which involve correlated random effects. MINQUE(0/1)method, which is a MINQUE method with all prior covariances setting 0 and all prior variances setting 1, is suitable for estimating variance and covariance components in these models. Mixed model approaches have advantage over ANOVA methods for the capacity of analyzing unbalanced data and complicated models. Some problems about estimation and hypothesis test by MINQUE method are discussed.     

THE MIXED-MODEL ANALYSIS OF VARIANCE APPLIED TO QUANTITATIVE GENETICS: BIOLOGICAL MEANING OF THE PARAMETERS

The mixed-model factorial analysis of variance has been used in many recent studies in evolutionary quantitative genetics. Two competing formulations of the mixed-model ANOVA are commonly used, the “Scheffe” model and the “SAS” model; these models differ in both their assumptions and in the way in which variance components due to the main effect of random factors are defined. The biological meanings of the two variance component definitions have often been unappreciated, however. A full understanding of these meanings leads to the conclusion that the mixed-model ANOVA could have been used to much greater effect by many recent authors. The variance component due to the random main effect under the two-way SAS model is the covariance in true means associated with a level of the random factor (e.g., families) across levels of the fixed factor (e.g., environments). Therefore the SAS model has a natural application for estimating the genetic correlation between a character expressed in different environments and testing whether it differs from zero. The variance component due to the random main effect under the two-way Scheffe model is the variance in marginal means (i.e., means over levels of the fixed factor) among levels of the random factor. Therefore the Scheffe model has a natural application for estimating genetic variances and heritabilities in populations using a defined mixture of environments. Procedures and assumptions necessary for these applications of the models are discussed. While exact significance tests under the SAS model require balanced data and the assumptions that family effects are normally distributed with equal variances in the different environments, the model can be useful even when these conditions are not met (e.g., for providing an unbiased estimate of the across-environment genetic covariance). Contrary to statements in a recent paper, exact significance tests regarding the variance in marginal means as well as unbiased estimates can be readily obtained from unbalanced designs with no restrictive assumptions about the distributions or variance-covariance structure of family effects.     

A penalized spline approach to functional mixed effects model analysis

In this article, we propose penalized spline (P-spline)-based methods for functional mixed effects models with varying coefficients. We decompose longitudinal outcomes as a sum of several terms: a population mean function, covariates with time-varying coefficients, functional subject-specific random effects, and residual measurement error processes. Using P-splines, we propose nonparametric estimation of the population mean function, varying coefficient, random subject-specific curves, and the associated covariance function that represents between-subject variation and the variance function of the residual measurement errors which represents within-subject variation. Proposed methods offer flexible estimation of both the population- and subject-level curves. In addition, decomposing variability of the outcomes as a between- and within-subject source is useful in identifying the dominant variance component therefore optimally model a covariance function. We use a likelihood-based method to select multiple smoothing parameters. Furthermore, we study the asymptotics of the baseline P-spline estimator with longitudinal data. We conduct simulation studies to investigate performance of the proposed methods. The benefit of the between- and within-subject covariance decomposition is illustrated through an analysis of Berkeley growth data, where we identified clearly distinct patterns of the between- and within-subject covariance functions of children's heights. We also apply the proposed methods to estimate the effect of antihypertensive treatment from the Framingham Heart Study data.     

The use of linear mixed model theory for the genetic analysis of repeated measures from clonal tests of forest trees. I. A focus on spatially repeated data

The paper reviews the linear mixed models (LMM) methodology that is suitable for the statistical and genetic analyses of spatially repeated measures collected from clonal progeny tests. For example, we consider a poplar clonal trial where progenies of different families are propagated by cuttings, and only one ramet per clone is planted on each block. Modeling covariance structures following the LMM theory allows improving genetic parameter estimation based on clonal testing. Besides variance components, we also obtained an estimate of the covariance between residuals (within clonal effects in two different blocks). This covariance is due to planting more than one ramet from the same genotype in the same trial, which generates correlated residual effects from different blocks. Its estimation can significantly improve the comparison among clones within a progeny test or between tests in a clonal testing network. Results indicate that the covariance is also a component of the genetic variance estimator and plays a significant role in assessing the variance of specific (micro) environmental effects. A positive covariance implies that ramets show a similar performance in more than one block. Thus, a larger and more positive covariance implies a stronger genetic effect controlling the expression of the trait in the local environment and a smaller variance of specific environmental effects. On the contrary, a negative covariance implies that the performance of individual ramets is affected by strong microenvironmental effects, specific to one or more blocks, which can also directly increase the within-clone variability.     

Random regression models using different functions to model test-day milk yield of Brazilian Holstein cows

We analyzed 152,145 test-day records from 7317 first lactations of Holstein cows recorded from 1995 to 2003. Our objective was to model variations in test-day milk yield during the first lactation of Holstein cows by random regression model (RRM), using various functions in order to obtain adequate and parsimonious models for the estimation of genetic parameters. Test-day milk yields were grouped into weekly classes of days in milk, ranging from 1 to 44 weeks. The contemporary groups were defined as herd-test-day. The analyses were performed using a single-trait RRM, including the direct additive, permanent environmental and residual random effects. In addition, contemporary group and linear and quadratic effects of the age of cow at calving were included as fixed effects. The mean trend of milk yield was modeled with a fourth-order orthogonal Legendre polynomial. The additive genetic and permanent environmental covariance functions were estimated by random regression on two parametric functions, Ali and Schaeffer and Wilmink, and on B-spline functions of days in milk. The covariance components and the genetic parameters were estimated by the restricted maximum likelihood method. Results from RRM parametric and B-spline functions were compared to RRM on Legendre polynomials and with a multi-trait analysis, using the same data set. Heritability estimates presented similar trends during mid-lactation (13 to 31 weeks) and between week 37 and the end of lactation, for all RRM. Heritabilities obtained by multi-trait analysis were of a lower magnitude than those estimated by RRM. The RRMs with a higher number of parameters were more useful to describe the genetic variation of test-day milk yield throughout the lactation. RRM using B-spline and Legendre polynomials as base functions appears to be the most adequate to describe the covariance structure of the data.     

A Building Block Model for Quantitative Genetics

We introduce a quantitative genetic model for multiple alleles which permits the parameterization of the degree, D, of dominance of favorable or unfavorable alleles. We assume gene effects to be random from some distribution and independent of the D's. We then fit the usual least-squares population genetic model of additive and dominance effects in an infinite equilibrium population to determine the five genetic components--additive variance sigma 2 a, dominance variance sigma 2 d, variance of homozygous dominance effects d2, covariance of additive and homozygous dominance effects d1, and the square of the inbreeding depression h--required to treat finite populations and large populations that have been through a bottleneck or in which there is inbreeding. The effects of dominance can be summarized as functions of the average, D, and the variance, sigma 2 D. An important distinction arises between symmetrical and nonsymmetrical distributions of gene effects. With symmetrical distributions d1 = -d2/2 which is always negative, and the contribution of dominance to sigma 2 a is equal to d2/2. With nonsymmetrical distributions there is an additional contribution H to sigma 2 a and -H/2 to d1, the sign of H being determined by D and the skew of the distribution. Some numerical evaluations are presented for the normal and exponential distributions of gene effects, illustrating the effects of the number of alleles and of the variation in allelic frequencies. Random additive by additive (a*a) epistatic effects contribute to sigma 2 a and to the a*a variance, sigma 2/aa, the relative contributions depending on the number of alleles and the variation in allelic frequencies.(ABSTRACT TRUNCATED AT 250 WORDS)     

Diallel analysis for sex-linked and maternal effects

Genetic models including sex-linked and maternal effects as well as autosomal gene effects are described. Monte Carlo simulations were conducted to compare efficiencies of estimation by minimum norm quadratic unbiased estimation (MINQUE) and restricted maximum likelihood (REML) methods. MINQUE(1), which has 1 for all prior values, has a similar efficiency to MINQUE(), which requires prior estimates of parameter values. MINQUE(1) has the advantage over REML of unbiased estimation and convenient computation. An adjusted unbiased prediction (AUP) method is developed for predicting random genetic effects. AUP is desirable for its easy computation and unbiasedness of both mean and variance of predictors. The jackknife procedure is appropriate for estimating the sampling variances of estimated variances (or covariances) and of predicted genetic effects. A t-test based on jackknife variances is applicable for detecting significance of variation. Worked examples from mice and silkworm data are given in order to demonstrate variance and covariance estimation and genetic effect prediction.