Bayesian inference on quantitative genetic (animal) model in BUGS & R

There is a paper by Lars Damgaard in Journal of Animal Science showing how to fit quantitative genetic (animal) model in WinBUGS and this of course extends also to other BUGS versions. This is really great. It is based on decomposition of additive relationship matrix A=TDT'. I will try to provide functions to ease the preparation of needed data for BUGS in my R package GeneticsPed. This will ease the exploration and first stage fitting of models on small or medium datasets. As noted by Damgaard, BUGS will not serve for genetic evaluation. I wonder if one can gain in speed with using Umacs or JAGS (site seems to be down at the moment).