epistemology twins

Twin Studies

The Basic Models

Suppose we're interested in how important genes or familial environment are in determining people's heights. For simplicity, suppose (1) each chromosome has a "height score" (2) person's height is a weighted sum of (a) all their chromosomes' "height scores", (b) their familial environment "height score", and (c) random noise. These are the assumptions that underlie the ACE model ACE model, which is, as far as I can tell, the most commonly used model in analyzing twin study data.

Computing the estimates of the amount of variance in height due to genes, common environment, and unshared environment is straightforward in this model if you have the correlation of monozygotic (identical) twins' heights and the correlation between dizygotic (fraternal) twins' heights.

To compute the genetic contribution, just take the difference between the two correlations and double it. To get the common environment contribution, just take the dizygotic twin correlation and subtract half the genetic contribution. The unshared environment is just whatever variance is left over, which includes things like measurement error and luck.

The other common model is the ADE model ADE model. The downside of the ADE model is that it requires the assumption that there is zero common environment. The upside is that it relaxes the assumption that all genetic effects are linear by allowing for dominance effects.

Because ACE ignores genetic dominance effects, its biased against genetic explanations and in favor of common environment ones. Because ADE ignores common environment, it's biased the other way. In my experience, ACE is far more common in the literature. There is a model that includes both common environment and genetic dominance effects: the ACDE model, but its typically not used because it is difficult/impossible to estimate all four parameters using a sample of only monozygotic and dizygotic twins.

Measurement Error

Suppose you're trying to estimate the heritability of height. Each person has some platonic height, but your dataset, sadly, only has measured height. Measured height can differ from platonic height due to a number of factors: rounding, human error, height's tendency to shrink throughout the day, etc.

Presumably, your genes affect your platonic height but not the sources of error that end up mediating measured height. That is, the only way genes affect your measured height is through your platonic height.

When we are computing the monozygotic and dizygotic correlations, we can adjust for this source of measurement error. Here's how we do so:

Let $r$ be the correlation between true height and measured height. Note, since true height is impossible to observed directly, we can't observe this correlation directly either.

To get around this issue, we can measure someone's height twice, the correlation between the two measured heights should be $r^2$. Hence, the square root of the "test-retest" correlation should be $r$.

If we want to get the true $r_{MZ}$, it turns out we can just divide it by $r^2$. Ditto for $r_{DZ}$. Hence, the heritability estimate becomes

$$ 2 \left( \frac{r_{MZ}}{r^2} - \frac{r_{DZ}}{r^2} \right) $$

Or just

$$ \frac{1}{r^2} 2 \left( r_{MZ} - r_{DZ} \right) $$

The hawk-eyed reader will note that this just means our adjusted estimate of heritability is the naive estimate of heritability divided by $r^2$.

Lest you think I'm a crank, here's people with credentials talking about the same issue Frisell, and here is a study where this kind of thing played a major issue Riemann.

Assortative Mating

Assortative mating will cause naive twin studies to over-estimate the effects of shared environment and underestimate the effects of additive genes. Sadly, the vast majority of twin studies ignore this issue.

How strong is this effect?

Suppose there are no common environment effect. If we know (a) the true heritability and (b) the observed phenotype correlation between the mom and dad, we can compute the expected conclusion of a naive twin study that assumes no assortative mating. We can do so using some simple code.

Running this code with different parameters yields this table:

h2Father-Mother Correlation
0.000.100.200.300.400.500.600.700.800.901.00
0%0.0020.0020.0020.0020.0020.0020.0020.0020.0020.0020.002
10%0.1020.1020.1010.1010.1010.1000.1000.0990.0990.0980.098
20%0.2030.2010.1990.1970.1950.1930.1920.1900.1880.1860.184
30%0.3030.2980.2940.2900.2860.2820.2780.2740.2710.2670.263
40%0.4020.3950.3870.3800.3730.3670.3600.3540.3480.3420.335
50%0.5020.4900.4790.4680.4570.4470.4370.4280.4190.4110.402
60%0.6020.5850.5680.5530.5380.5240.5110.4980.4860.4750.463
70%0.7020.6780.6560.6360.6160.5980.5810.5640.5490.5340.520
80%0.8010.7710.7420.7160.6910.6680.6470.6270.6080.5900.572
90%0.9010.8620.8270.7940.7640.7360.7100.6860.6630.6420.621
100%0.9990.9520.9090.8690.8330.8000.7690.7410.7150.6900.667

Note: Even though I used a sample of two million hypothetical twins, there is still some sampling error, which is why (for instance) the first row is close to zero but doesn't equal it precisely. This is likewise why the first column is close to the actual heritability but not exactly the same.

We can use this table to adjust twin studies for assortative mating. For instance, if a twin study finds a heritability of 72% but there is a correlation of 0.3 between fathers and mothers for that trait, the true heritability is probably a smidge above 80%.

Age

In general, the heritability of psychiatric traits increase with age Bergen and the common environment of traits decrease with age Meta-analysis of the heritability of human traits based on fifty years of twin studies . For this reason, estimates derived from studies of, say, young children should generally not be used to make conclusions about adults (and vice versa).

One place this shows up is educational attainment, where the years of education obtained has a large common environment component - probably because it is the result largely of events and decisions that occur in childhood, when parents still have large influence over their children. It notably contrasts with standardized high school test scores, which generally have approximately zero common environment component, and with IQ, which is probably the trait most famous for having a heritability that shifts with age.

As an aside, I am forced to wonder whether this has anything to do with the fadeout of educational interventions over time.

Other Methodological Shortcomings

In addition to the issues with using ACE or ADE, twin studies typically suffer from other methodological shortcomings

  • Twin studies often use questionnaires, often not verifying the reported information including whether the pair of people in question are even twins . Misclassifying people a fraternal twin as identical (or vice-versa) biases estimates against genes and towards common environment. The unreliability of self-reported data is rando and is similar to the measurement error discussed above and will generally bias estimates of genes and shared environment downwards - the exception is if the unreliability is non-random (e.g. twins are likely to misreport in similar ways).
  • Twin studies frequently draw on (non-random) convenience samples.
  • When twin studies use siblings rather than (or in addition to) dizygotic twins, this causes birth-order effects to be interpreted as genetic effects, biasing heritability estimates upwards and common environment estimates downwards.
  • If monozygotic twins are treated more similarly than dizygotic twins (or siblings if they're used), this biases common and unshared environment downwards while biasing genetic estimates upwards.
  • Maybe twins make choices to be more (or less) similar. This would bias estimates of unshared environment upwards. Alternatively, maybe this force is stronger in identical twins than fraternal twins, which would bias estimates of heritability upwards. This hypothesis is difficult to evaluate, but here are some related studies Fortuna Mark.
  • To the literature's credit, publication bias seems generally minor, with one major exception: the choice of ACE or ADE. Generally, studies start with an ACE model and, if this results in a negative estimate for common environment, change to an ADE model. This use of choice causes a positive bias for dominance genetic effects, a positive bias for common environment effects, and a negative bias for additive genetic effects - especially in small sample sizes Meta-analysis of the heritability of human traits based on fifty years of twin studies (see also the supplement).
  • TODO: gene-environment interactions Kaprio Avinun
  • TODO Grotzinger

Despite the widespread belief that smaller populations should have less additive genetic variance, the empirical evidence is scant Wood.

todo Røysamb

Appendix: Assortative Mating II

If a study is assuming zero common environment, then it is common to estimate heritability as just $h^2 = r_{MZ}$. In this case, assortative mating will bias heritability estimates upwards, and the corrections needed to account for it are different:

h2Father-Mother Correlation
0.000.100.200.300.400.500.600.700.800.901.00
0%0.0010.0010.0010.0010.0010.0010.0010.0010.0010.0010.001
10%0.1000.1040.1090.1130.1170.1220.1260.1300.1340.1390.143
20%0.1990.2070.2150.2230.2300.2370.2450.2520.2590.2660.273
30%0.2990.3090.3200.3290.3390.3480.3570.3660.3750.3830.391
40%0.3990.4110.4220.4330.4440.4540.4640.4730.4820.4910.500
50%0.4990.5110.5230.5340.5450.5550.5650.5740.5830.5920.600
60%0.5990.6110.6220.6320.6420.6520.6610.6690.6770.6850.692
70%0.6990.7100.7190.7280.7360.7440.7520.7590.7650.7720.778
80%0.8000.8070.8140.8210.8270.8330.8380.8440.8480.8530.857
90%0.9000.9040.9080.9120.9150.9180.9210.9240.9260.9290.931
100%1.0001.0001.0001.0001.0001.0001.0001.0001.0001.0001.000
Elks, C. E., Den Hoed, M., Zhao, J. H., Sharp, S. J., Wareham, N. J., Loos, R. J., & Ong, K. K. (2012). Variability in the heritability of body mass index: a systematic review and meta-regression. Frontiers in endocrinology, 3, 29. https://doi.org/10.3389/fendo.2012.00029 Wikipedia contributors. (2021, March 30). ACE model. In Wikipedia, The Free Encyclopedia. Retrieved 23:22, October 10, 2021, from https://en.wikipedia.org/w/index.php?title=ACE_model&oldid=1015154064 Wikipedia contributors. (2018, December 26). ADE model. In Wikipedia, The Free Encyclopedia. Retrieved 23:23, October 10, 2021, from https://en.wikipedia.org/w/index.php?title=ADE_model&oldid=875467243 Ozaki, K., Toyoda, H., Iwama, N., Kubo, S., & Ando, J. (2011). Using non-normal SEM to resolve the ACDE model in the classical twin design. Behavior Genetics, 41(2), 329-339. https://doi.org/10.1007/s10519-010-9386-5 Silver, I. A., Nedelec, J. L., Segal, N. L., & Lonergan, H. (2021). Heteropaternal Siblings Misclassified as Dizygotic Twins: A Potential Biasing Factor for Heritability Estimates?. Behavior Genetics, 51(2), 137-143. Fortuna, K., Goldner, I., & Knafo, A. (2010). Twin relationships: A comparison across monozygotic twins, dizygotic twins, and nontwin siblings in early childhood. Family Science, 1(3-4), 205-211. https://doi.org/10.1080/19424620.2010.569367 Mark, K. M., Pike, A., Latham, R. M., & Oliver, B. R. (2017). Using twins to better understand sibling relationships. Behavior genetics, 47(2), 202-214. https://doi.org/10.1007/s10519-016-9825-z Polderman, T. J., Benyamin, B., De Leeuw, C. A., Sullivan, P. F., Van Bochoven, A., Visscher, P. M., & Posthuma, D. (2015). Meta-analysis of the heritability of human traits based on fifty years of twin studies. Nature genetics, 47(7), 702-709. https://doi.org/10.1038/ng.3285 Polderman, T. J., Benyamin, B., De Leeuw, C. A., Sullivan, P. F., Van Bochoven, A., Visscher, P. M., & Posthuma, D. (2015). Supplementary Figures 1–12, Supplemen tary Note and Supplementary Tables 1–19, 22–24 and 26–31. Nature genetics, 47(7), 702-709. https://static-content.springer.com/esm/art%3A10.1038%2Fng.3285/MediaObjects/41588_2015_BFng3285_MOESM4_ESM.pdf Bergen, S. E., Gardner, C. O., & Kendler, K. S. (2007). Age-related changes in heritability of behavioral phenotypes over adolescence and young adulthood: a meta-analysis. Twin Research and Human Genetics, 10(3), 423-433. https://doi.org/10.1375/twin.10.3.423 Wood, J. L., Yates, M. C., & Fraser, D. J. (2016). Are heritability and selection related to population size in nature? Meta‐analysis and conservation implications. Evolutionary applications, 9(5), 640-657. https://doi.org/10.1111/eva.12375 Røysamb, E., & Tambs, K. (2016). The beauty, logic and limitations of twin studies. https://fhi.brage.unit.no/fhi-xmlui/bitstream/handle/11250/2431737/Roysamb_2016_The.pdf?sequence=2 Kaprio, J. (2012). Twins and the mystery of missing heritability: the contribution of gene–environment interactions. Journal of internal medicine, 272(5), 440-448. https://doi.org/10.1111/j.1365-2796.2012.02587.x Grotzinger, A. D., de la Fuente, J., Nivard, M. G., & Tucker-Drob, E. M. (2021). Pervasive Downward Bias in Estimates of Liability Scale Heritability in GWAS Meta-Analysis: A Simple Solution. medRxiv. https://doi.org/10.1101/2021.09.22.21263909 Avinun, R., & Knafo, A. (2014). Parenting as a reaction evoked by children’s genotype: A meta-analysis of children-as-twins studies. Personality and Social Psychology Review, 18(1), 87-102. https://doi.org/10.1177/1088868313498308 Frisell, T., Öberg, S., Kuja-Halkola, R., & Sjölander, A. (2012). Sibling comparison designs: Bios from non-shored confounders and measurement error. Epidemiology, 713-720. https://doi.org/10.1097/EDE.0b013e31825fa230 Riemann, R., & Kandler, C. (2010). Construct validation using multitrait‐multimethod‐twin data: The case of a general factor of personality. European Journal of Personality, 24(3), 258-277. https://doi.org/10.1002/per.760