Quantitative Genetics and Twin Studies

Similar documents
Midterm#1 comments. Overview- chapter 6. Recombination. Recombination 1 st sense

Course Information. Computational Genetics Lecture 1. Course Prerequisites. Course Goals

STAT 536: Genetic Drift

Student Activity 3: Single Factor ANOVA

Technical Appendix: Childhood Family Structure and Schooling Outcomes: Evidence for Germany

Tests for the Ratio of Two Poisson Rates

Non-Linear & Logistic Regression

Lect 11: Inbreeding, evolution at multiple loci. Exam 1, Monday, 2 October. Consequences of Inbreeding: Genotype frequencies

Discrete Mathematics and Probability Theory Spring 2013 Anant Sahai Lecture 17

Population bottleneck : dramatic reduction of population size followed by rapid expansion,

AS a powerful design in the case of quantitative traits the ESP method remains powerful when used with caution.

Haplotype Frequencies and Linkage Disequilibrium. Biostatistics 666

Continuous Random Variable X:

Discrete Mathematics and Probability Theory Summer 2014 James Cook Note 17

Prediction of Rates of Inbreeding in Populations Selected on Best Linear Unbiased Prediction of Breeding Value

Mendelian Genetics. References for lectures 1 and 2. Overview of Genes I. Overview of Genes I. Variation is widespread in nature

Lecture INF4350 October 12008

Lecture 3 Gaussian Probability Distribution

CAPITAL ASSET PRICING MODEL (CAPM)

For the percentage of full time students at RCC the symbols would be:

Lecture 1: Introduction to integration theory and bounded variation

Data Assimilation. Alan O Neill Data Assimilation Research Centre University of Reading

Maximum Likelihood Estimation for Allele Frequencies. Biostatistics 666

A Brief Review on Akkar, Sandikkaya and Bommer (ASB13) GMPE

Tutorial 4. b a. h(f) = a b a ln 1. b a dx = ln(b a) nats = log(b a) bits. = ln λ + 1 nats. = log e λ bits. = ln 1 2 ln λ + 1. nats. = ln 2e. bits.

Lecture 2: January 27

Chapter 9: Inferences based on Two samples: Confidence intervals and tests of hypotheses

4.1. Probability Density Functions

1B40 Practical Skills

HETEROSIS, a term to describe the superiority of

Is there an easy way to find examples of such triples? Why yes! Just look at an ordinary multiplication table to find them!

Jackson 2.26 Homework Problem Solution Dr. Christopher S. Baird University of Massachusetts Lowell

The steps of the hypothesis test

THE detection of quantitative trait loci (QTL) in

A-Level Mathematics Transition Task (compulsory for all maths students and all further maths student)

Genetic Programming. Outline. Evolutionary Strategies. Evolutionary strategies Genetic programming Summary

LECTURE # How does one test whether a population is in the HW equilibrium? (i) try the following example: Genotype Observed AA 50 Aa 0 aa 50

Fault Modeling. EE5375 ADD II Prof. MacDonald

1. Weak acids. For a weak acid HA, there is less than 100% dissociation to ions. The B-L equilibrium is:

Section 11.5 Estimation of difference of two proportions

The area under the graph of f and above the x-axis between a and b is denoted by. f(x) dx. π O

Math 31S. Rumbos Fall Solutions to Assignment #16

Factors affecting the phonation threshold pressure and frequency

W. We shall do so one by one, starting with I 1, and we shall do it greedily, trying

1.3 Regular Expressions

Riemann is the Mann! (But Lebesgue may besgue to differ.)

1.4 Nonregular Languages

Lecture 21: Order statistics

Continuous probability distributions

Numerical integration

Affected Sibling Pairs. Biostatistics 666

Operations with Polynomials

Intro to Nuclear and Particle Physics (5110)

Comparison Procedures

Robust Predictions in Games with Incomplete Information

Model building & assumptions. Matt Keller, Sarah Medland, Hermine Maes TC21 March 2008

Partitioning Genetic Variance

I1 = I2 I1 = I2 + I3 I1 + I2 = I3 + I4 I 3

BS 50 Genetics and Genomics Week of Oct 3 Additional Practice Problems for Section. A/a ; B/B ; d/d X A/a ; b/b ; D/d

Solution for Assignment 1 : Intro to Probability and Statistics, PAC learning

SUMMER KNOWHOW STUDY AND LEARNING CENTRE

Continuous Random Variables

Downloaded From:

SOLUTIONS FOR ADMISSIONS TEST IN MATHEMATICS, COMPUTER SCIENCE AND JOINT SCHOOLS WEDNESDAY 5 NOVEMBER 2014

5.1 How do we Measure Distance Traveled given Velocity? Student Notes

1 Probability Density Functions

THERMAL EXPANSION COEFFICIENT OF WATER FOR VOLUMETRIC CALIBRATION

Reference. Vector Analysis Chapter 2

Reading from Young & Freedman: For this topic, read the introduction to chapter 24 and sections 24.1 to 24.5.

Reversals of Signal-Posterior Monotonicity for Any Bounded Prior

CS667 Lecture 6: Monte Carlo Integration 02/10/05

Decision Science Letters

Properties of Integrals, Indefinite Integrals. Goals: Definition of the Definite Integral Integral Calculations using Antiderivatives

ECO 317 Economics of Uncertainty Fall Term 2007 Notes for lectures 4. Stochastic Dominance

Legend: S spotted Genotypes: P1 SS & ss F1 Ss ss plain F2 (with ratio) 1SS :2 WSs: 1ss. Legend W white White bull 1 Ww red cows ww ww red

p-adic Egyptian Fractions

13: Diffusion in 2 Energy Groups

8 Laplace s Method and Local Limit Theorems

Modeling IBD for Pairs of Relatives. Biostatistics 666 Lecture 17

1 PYTHAGORAS THEOREM 1. Given a right angled triangle, the square of the hypotenuse is equal to the sum of the squares of the other two sides.

( ) as a fraction. Determine location of the highest

5 Probability densities

Testing categorized bivariate normality with two-stage. polychoric correlation estimates

AB Calculus Review Sheet

Fig. 1. Open-Loop and Closed-Loop Systems with Plant Variations

Math Calculus with Analytic Geometry II

14.4. Lengths of curves and surfaces of revolution. Introduction. Prerequisites. Learning Outcomes

Section 6: Area, Volume, and Average Value

( ) where f ( x ) is a. AB Calculus Exam Review Sheet. A. Precalculus Type problems. Find the zeros of f ( x).

A signalling model of school grades: centralized versus decentralized examinations

8.6 The Hyperbola. and F 2. is a constant. P F 2. P =k The two fixed points, F 1. , are called the foci of the hyperbola. The line segments F 1

Breeding Values and Inbreeding. Breeding Values and Inbreeding

Math 4310 Solutions to homework 1 Due 9/1/16

Robust Predictions in Games with Incomplete Information

Today s Outline. Inheritance. Darwin & Mendel near miss. Gregor Johann Mendel. Terms Punnett-square. Modern Synthesis

For convenience, we rewrite m2 s m2 = m m m ; where m is repeted m times. Since xyz = m m m nd jxyj»m, we hve tht the string y is substring of the fir

List all of the possible rational roots of each equation. Then find all solutions (both real and imaginary) of the equation. 1.

Tries and suffixes trees

Diffraction. Diffraction and Polarization. Diffraction. Diffraction. Diffraction I. Angular Spread: θ ~ λ/a

Chapter 2 Fundamental Concepts

Transcription:

Count Count Count Count Quntittive Genetics nd Twin Studies n Introduction! co de Geus -Dept. Biologicl Psychology -Netherlnds Twin Register msterdm, the Netherlnds 600 N = 6602 M = 48,27 SD = 25,0 75 N = 794 M = 84,76 SD=,32 400 50 200 25 0 Individul differences 40 30 20 25 50 75 00 PQ neuroticism N =423 M = 9,9 SD= 7,02 0 60 80 00 20 mbultory hert rte (morning) 60 40 N = 698 M = 8,22 SD= 8,25 0 20-9 0 9 8 27 36 Distolic BP rectivity (Tone voidnce RT tsk) 0 0,00 20,00 30,00 40,00 Cortisol t wkening

Vrince Components Phenotype = µ + G + + G* Vrince (P) = vr (µ +G + + G*) = vr(g) + vr() + vr(g*) + 2Cov(G,)+2Cov(G,G*)+2Cov(,G*) ll G covs set to zero: = vr (G) + vr () + vr (G*) vr(g*) set to zero: Vr(P) = vr(g) + vr() llelic vrition llele of Fther Mother 2 4 (C) 3 C3 (C) 9 C9 Gene = unit of heredity Locus = loction (of gene) on chromosome llele = different vrints of the gene t the sme locus Homozygous genotype = pternl nd mternl lleles re identicl Heterozygous genotype = pternl nd mternl lleles re different

Dichotomous trits: Huntington s disese x H h h h Genotypes: ¼Hh ¼Hh ¼hh ¼hh ffected Unffected MIM4300 (993) 4p6.3 dysfunctionl Huntingtin protein Dominnt disese Dichotomous trits: PKU x P p P p Genotypes: ¼PP ¼Pp ¼ pp ¼ pp Unffected, 2 Crriers ffected MIM26600 (984) 2q24. dysfunctionl Phenyllnine dehydroxylse Recessive disese

Quntittive trits llele of Fther Mother b B 2 C3 4 C9 Regultory zone b..b incresing lleles for better trnscription: more protein synthesis xon 2..4 incresing lleles for more efficient isoform of the protein xon 5 C..C9 incresing lleles for dysfunctionl protein - d = -4 67 75 bb Bb m (79) 9 = -2 = +2 BB xmple: stture gene with bi-llelic polymorphism in the regultory zone of growth protein yields three types of persons with three different heights: 67, 75, 9 cm. breeding (thought!) experiment P homozygote seletion (BB fthers * bb mothers) => F, ll Bb F intercross (Bb fthers * Bb mothers) => F2, BB, Bb, bb Genotype (i) BB Bb bb Frequency (f) ¼ ½ ¼ Genotypic effect (x) d - f * x ¼ ½ d -¼ Men genetic effect (µ g ) = f i *x i = 2*¼ + -4*½ + -2*¼ = -2 Genetic vrince = f i *(x-µ g ) 2 vr(p) = vr() Vr + (G) vr (D) + vr() =¼*(2--2) 2 + ½*(-4 - -2) 2 + ¼*(-2 - -2) 2 =76 The contribution of the stture locus to the genetic vrince in stture

Polygenetic trits For multiple genes, dditive or dominnt, with vrying llelic frequencies, genetic effects on the sme trit my be summed nd cn be shown - under the centrl limit theorem - to yield normlly distributed men genetic effect with vrince Vr() + Vr(D). = B = C = smll effect 2 = B2 = C2 = lrge effect Heritbility How do we estblish the reltive contribution of genetic influences to trit vrince (.k.. Heritbility or H 2 )?? Fmily studies to test the strength of the correltion between genetic reltedness nd trit resemblnce!

Fmilil resemblnce in disese trits Popultion Hlf Sibling Grnd child Nephew/Niece Uncle/unt Full Sibling Children Prents 0% genetic resemblnce 25% genetic resemblnce 50% genetic resemblnce 0 2 4 6 8 0 2 4 Schizophreni Prevlence Fmilil resemblnce in quntittive trits Tmbs et l. (992) collected resting blood pressure in 74994 Norwegins from three genertions

Shred (fmily) environment Fmily members shre prt of their environment (SS, Neighborhood, Diet, Prentl rering style, School, etc.). This common environment cn lso contribute to fmily resemblnce!! Vr(P) = vr () + vr (D) + vr () = vr () + vr (D) + vr (C) + vr () C = common environment = unique environment Monozygotic MZ twins re 00% geneticlly identicl

Dizygotic DZ twins (nd singleton siblings) shre on verge 50% of their genetic vrition. Clssicl Twin design The clssicl twin design for disese trits dtes from the 920's, when dignosis of monozygotic (MZ) nd dizygotic (DZ) twins becme relible. It is ssumed tht twins in the sme household re exposed to shred environmentl risk fctors to the sme extent whether they re MZ or DZ. Under the equl environment ssumption, differences in intrpir resemblnce by zygosity reflect genetic rther thn common environmentl fctors.

Intrpir twin correltion for stture MZ mles 20 200 90 80 70 rmz= 0,95 Rndom sme-ged mle pirs 20 stture twin 2 60 50 40 50 60 70 80 90 200 20 200 90 stture twin DZ mles 80 70 200 90 rdz = 0,60 60 80 stture 50 40 70 40 50 60 70 80 90 200 stture dolescent Dutch twins, 5-9 yr stture twin 2 60 50 50 60 stture twin 70 80 90 200 Twin correltions: rules of thumb Vr () / totl Vr (P) = 2 = 2 (rmz - rdz) Vr () / totl Vr (P) = e 2 = rmz Vr (C) / totl Vr (P) = c 2 = 2rDZ - rmz Vr (D) / totl Vr (P) = d 2 = 2rMZ - 4rDZ compliction: Vrince = + D + C + ; MZ covrince: + D + C; DZ covrince: ½ + ¼ D + C 3 equtions, 4 unknowns D&C cnnot be estimted simultneously, either C or D hs to be ssumed to be 0. Usully, D is modeled only if MZ >> 2*DZ.

Heritbility of stture Intrpir correltions: Monozygotic: rmz= 0,95 Dizygotic: rdz = 0,60 = 2 * (.0.95 0.60) = 0.70 C = 2 * 0.60 0.95 = 0.25 = 0.95 = 0.05 D = 2 * 0.95 4 * 0.60 = -0.50! % individul vrince explined by Heritbility: 70 % nvironmentl influences: 30 % (shred 25%, unique 5%) Univrite pth model for twin studies r MZ = r DZ/sib = r MZ =, r DZ/sib = 0.5 C C e c c e Stture Twin Stture Twin 2 Cov (MZ) = (**) + (c**c) = 2 +c 2 Cov (DZ) = (*0.5*) + (c**c) = 0.5 2 + c 2 Vr (stture) = (**) + (c**c) + (e**e) = 2 + c 2 + e 2

Model fitting Model C 0.5/ C lterntive models compred by likelihoodrtio tests. e c P c e P2 Model 2 0.5/ e P ge Model 3 P2 e Model 4: No sex differences Genetic nlysis of stture: Wht is the best model? Model Chi df p. ge, Genes, C, (sex dif) 23.9 2.30 2. ge, Genes, (sex diff) 26.6 23.27 3. ge, Genes, (no sex diff) 78.2 25.00 Best model: (no C) with sex differences: Mles Femles Genes 75% 89% Unique nvironment 5% 0% ge 20% % (NB 5-9!)

s: r MZ =, r DZ/sib = 0.5 C s: r MZ = r DZ/sib = Multivrite extension C C 2 2 C C 2 2 c c 22 c 2 2 c 2 22 2 22 c c 22 Stture Twin e Weight Twin e 22 Stture Twin 2 e e 2 e 2 Weight Twin 2 e 22 2 2 Genetics of trcking: mplifiction or emergence of new genes I = 2 = n = II = 2 = 2 n = 0.5 III 22 > 0 nn > 0 2 n 2 22 n2 nn n.

sy to dd observed genotypes on mrkers or cndidte genes r MZ = r DZ/sib = r MZ =, r DZ/sib = 0.5 e C c r MZ =, r DZ/sib = π^ C c e Q Q q q P twin P twin2 Q = Quntittive Trit Locus, mesured genotype π =IBD/2, mesure for sib resemblnce t the QTL Not limited to twin reserch! r sib = 0.5 e r sib = π^ e Q Q q q P sib P sib2