Search PubMed⌕ Search

Biomedical subjects

S C Heath

Publications and source records attributed to S C Heath.

23 records · Page 2Linked to original sources

Markov chain Monte Carlo segregation and linkage analysis for oligogenic models.

A new method for segregation and linkage analysis, with pedigree data, is described. Reversible jump Markov chain Monte Carlo methods are used to implement a sampling scheme in which the Markov chain can jump between parameter subspaces corresponding to models with different numbers of quantitative-trait loci (QTL's). Joint estimation of QTL number, position, and effects is possible, avoiding the problems that can arise from misspecification of the number of QTL's in a linkage analysis. The method is illustrated by use of a data set simulated for the 9th Genetic Analysis Workshop; this data set had several oligogenic traits, generated by use of a 1,497-member pedigree. The mixing characteristics of the method appear to be good, and the method correctly recovers the simulated model from the test data set. The approach appears to have great potential both for robust linkage analysis and for the answering of more general questions regarding the genetic control of complex traits.

Computer Simulation↗

Markov chain Monte Carlo methods for radiation hybrid mapping.

The ordering of genetic loci is central to genetic mapping at all levels. Markov chain Monte Carlo (MCMC) techniques can provide estimates of the posterior density of orders while accounting naturally for missing data, data errors, and unknown parameters. MCMC sampling schemes have been proposed for mapping problems such as linkage mapping and radiation hybrid mapping. The sampling schemes tend, however, to suffer from poor mixing caused by strong correlations between the model parameters. The method described here investigates the effect of using a modified sampling scheme, simulated tempering, on the mixing characteristics of the Markov chain. The method is illustrated by the analysis of haploid radiation hybrid mapping data; the principles are, however, applicable to a range of mapping problems. The results demonstrate that simulated tempering greatly improves the performance of the MCMC sampling scheme. For the radiation hybrid problem, the approach is probably not suitable for simultaneously ordering very large number of loci (> 100); it could, however, be useful for fine scale mapping of subsections of chromosomes.

Algorithms↗

Imperfect genes, Fisherian mutation and the evolution of sex.

In this paper we present a mathematical model of mutation and selection that allows for the coexistence of multiple alleles at a locus with very small selective differences between alleles. The model also allows for the determination of fitness by multiple loci. Models of this sort are biologically plausible. However, some previous attempts to construct similar models have assumed that all mutations produce a decrease in fitness, and this has led to a tendency for the average fitness of population members to decline when population numbers are finite. In our model we incorporate some of the ideas of R. A. FISHER, so that both deleterious and beneficial mutations are possible. As a result, average fitness tends to approach a stationary distribution. We have used computer simulation methods to apply the Fisherian mutation model to the problem of the evolution of sex and recombination. The results suggest that sex and recombination can provide very large benefits in terms of average fitness. The results also suggest that obligately sexual species will win ecological competitions with species that produce a substantial fraction of their offspring asexually, so long as the number of sites under selection within the genomes of the competing species is not too small and the population sizes are not too large. Our model focuses on fertility selection in an hermaphroditic plant. However, the results are likely to generalize to a wide variety of other situations as well.

Biological Evolution↗

Rates of change of genetic parameters of body weight in selected mouse lines.

A method based on the animal model is described which allows the estimation of continuous changes in variance components over time using restricted maximum likelihood (REML). The method was applied to the analysis of a selection experiment in which a foundation population formed from a cross between two inbred strains of mice (C57BL/6J and DBA/2J) was divergently selected for 6 week body weight over 20 generations. The analysis suggested that there was an increase in phenotypic variance of about 50% in the low selected lines over the course of the experiment which was attributed to increases in the environmental and additive variance components. Variance changes in the High selected lines were generally smaller than in the Low lines, although there was an estimated 20% increase in the environmental variance. Simple models to explain these effects involving dominance, linkage and epistasis were explored. Testing which of these was responsible for the variance changes noted in this experiment (if any) is difficult, although the epistasis and dominance models require less stringent conditions than the linkage model, and the dominance model is supported by evidence of heterosis in the F1.

Analysis of Variance↗

Generating consistent genotypic configurations for multi-allelic loci and large complex pedigrees.

Markov chain Monte Carlo techniques have been applied to many different genetic analysis problems. Gibbs sampling in particular has been widely used because of its simplicity and because it can be extended to arbitrarily complex pedigrees and genetic models (albeit with modifications for multi-allelic loci). Gibbs sampling requires an initial genotypic configuration, consistent with observed data, the generation of which is not trivial with large complex pedigrees and multi-allelic loci. A method to generate feasible genotype configurations in these circumstances using a combination of peeling and genotype elimination is described. The method is illustrated using two complex multi-generation pedigrees, one real and one simulated, each partially typed for one highly polymorphic marker locus.

Alleles↗