Abstract
Under gynogenetic reproduction, offspring receive genes only from their dams and completely homozygous offspring are produced within one generation. When gynogenetic reproduction is applied to fully inbred individuals, homozygous clone lines are produced. A mixed model method was developed for breeding value and variance component estimation in gynogenetic families, which requires the inverse of the numerator relationship matrix. A general method for creating the inverse for a population with unusual relationships between animals is presented, which reduces to simple rules as is illustrated for gynogenetic populations. The presence of clones in gynogenetic populations causes singularity of the numerator relationship matrix. However, clones can be regarded as repeated observations of the same genotype, which can be accommodated by modifying the incidence matrix, and by considering only unique genotypes in the estimation procedure. Optimum gynogenetic sib family sizes for estimating heritabilities and estimates of their accuracy were derived and compared to those for conventional full-sib designs. This was done by means of a deterministic derivation and by stochastic simulation using Gibbs sampling. Optimum family sizes were smallest for gynogenetic families. Only for low heritabilities, there was a small advantage in accuracy under the gynogenetic design.