A recursive procedure to compute the gametic relationship matrix and its inverse for marked QTL clusters

Recently, a variety of mixed linear models have been proposed for marker-assisted prediction of the effects of quantitative trait loci (QTLs) in outbred populations of animals. One of them addresses the effects of a cluster of linked QTLs, or those of a particular chromosomal segment, marked by DNA marker(s) and requires that the inverse of the corresponding gametic relationship matrix whose elements are the conditional expected values of the identity-by-descent (IBD) proportions between gametes for individuals be evaluated. Here, for a model of this type, utilizing the property of the IBD set and using the information on the joint gametogenesis processes at the flanking marker loci, we present a recursive method to systematically calculate the elements of the gametic relationship matrix and its inverse. A numerical example is given to illustrate the proposed computing procedure.