Sessão Temática 2 Análise Bayesiana

1. Sessão Temática 2Análise Bayesiana Utilizando a abordagem Bayesiana no mapeamento de QTL´s Roseli Aparecida Leandro ESALQ/USP 11o SEAGRO / 50ª RBRAS Londrina, Paraná 04 a 08 de Julho de 2005

2. Colaboradores • Prof. Dr. Cláudio Lopes Souza Jr. • Prof. Dr. Antônio Augusto Franco Garcia • (Departamento de Genética ESALQ/USP) • Elisabeth Regina de Toledo • (PPG Estatística e Experimentação Agronômica, ESALQ/USP)

3. Qualitative trait Mendelian gene

4. Quantitative trait

5. Bayesian mapping of QTL Geneticists are often interested in locating regions in the chromosome contributing to phenotypic variation of a quantitative trait. Quantitative Trait Loci

6. Location • Effects : Additive, dominance QTL Genetics Markers

7. Escala dos Valores Genotípicos em que: d/a é o grau de dominância

8. Chromosomal regions of known location • Do not have a physiological causal association to the trait under study Genetic Markers

9. Genetics Markers

10. By studying the joint pattern of inheritance of the markers and trait Inferences can be made about the number, location and effects of the QTL affecting trait.

11. Experimental Design • Offspring data: Divergent inbred lines Backcross ( code 0=aa, 1=Aa ) (Recessive) F2 (code –1=aa, 0=Aa, 1=AA) • Reason: maximize linkage desiquilibrium

12. F2 Design

13. Data set

14. QTL phenotype model • One QTL

15. Multiple QTL phenotype Model

16. Our aim is to make joint inference about the number of QTL, their positions (loci) and the sizes of their effects. • Assume that a linkage map has been developed for the genome.

17. Genetic Map

18. 0 < r = fração de recombinação < 0.5

19. Classic approach • Interval mapping (Lander & Botstein,1989) • Least squares method (Haley & Knott,1992) • Composite interval mapping (Jansen, 1993; Jansen and Stam, 1994; Zeng 1993, 1994)

20. Bayesian approach • Satagopan et al. (1996) • Satagopan & Yandell (1998) • Sillanpää & Arjas (1998)

21. The joint posterior distribution of all unknowns (s, , Q, ) is proportional to

22. In practice, we observe the phenotypic trait . and the marker genotypes but NOT the QTL genotypes . • For convenience consider only one linkage group with ordered markers {1,2,...,m}. Assume that genotypes: • The markers are assumed to be at known distances

23. The conditional distribution * assuming the loci segregate independently ** under Haldane assumption of independent recombination

24. The marginal likelihood of the parameters s,  and  for the ith individual may be obtained from the joint distribution of traits and QTL genotypes. by summing over the set of all possible QTL genotypes for the ith individual,

25. Therefore, • When the data Y are n independent observations, the marginal likelihood for the trait data is the product over individuals, a familiar misture model likelihood,

26. The joint likelihood is a mixture of densities, and hence, is difficult to evaluate when there are multiple QTL. • The joint posterior distribution of all unknowns (s, , Q, ) is proportional to

27. A Bayesian approach combined with reversible jump MCMC is well suited for QTL studies

28. Random-sweep Metropolis-Hastings algorithm for general state spaces (Richardson and Green, 1997) • Suppose current state of the chain indexed by s.

29. The chain can (1) move to a “birth” step (number of loci s  s+1 ) (2) move to a “death” step (number of loci s  s-1 ) (3) continue with “current” number (s) of loci

30. S-1 S S+1

31. Simulation • Simulated F2 intercross • n=250 • 1 cromossome • 2 QTL

35. Referências • Satagopan, J. M.; Yandell, B. S. (1998) • Bayesian model determination for quantitative trait