Cited by 144Genetic Associations and EpidemiologyGenetic and phenotypic traits in livestockGenetic Mapping and Diversity in Plants and Animals

Mendelian imputation of parental genotypes improves estimates of direct genetic effects

Open access
Young, Alexander I., Nehzati, Seyed Moeen, Benonisdottir, Stefania, Okbay, Aysu, Jayashankar, Hariharan, Lee, Chanwook, Cesarini, David, Benjamin, Daniel J. +2 more·2022·Nature Genetics

Published version of a preprint

Abstract

Effects estimated by genome-wide association studies (GWASs) include effects of alleles in an individual on that individual (direct genetic effects), indirect genetic effects (for example, effects of alleles in parents on offspring through the environment) and bias from confounding. Within-family genetic variation is random, enabling unbiased estimation of direct genetic effects when parents are genotyped. However, parental genotypes are often missing. We introduce a method that imputes missing parental genotypes and estimates direct genetic effects. Our method, implemented in the software package snipar (single-nucleotide imputation of parents), gives more precise estimates of direct genetic effects than existing approaches. Using 39,614 individuals from the UK Biobank with at least one genotyped sibling/parent, we estimate the correlation between direct genetic effects and effects from standard GWASs for nine phenotypes, including educational attainment ( r = 0.739, standard error (s.e.) = 0.086) and cognitive ability ( r = 0.490, s.e. = 0.086). Our results demonstrate substantial confounding bias in standard GWASs for some phenotypes.

The claims and statistics of this paper haven’t been checked yet.

Run the check — free, no account needed. It reads the open-access text and checks the sources it cites.

A few free checks a day · papers already verified open instantly · results are public

Integrity signals

Registry

Crossref · 10.1038/s41588-022-01085-0

Cites

0 citation errors · 0 not checkable

Statistics

Claims the paper makes

The claims and statistics of this paper haven’t been checked yet.

Run the check — free, no account needed. It reads the open-access text and checks the sources it cites.

A few free checks a day · papers already verified open instantly · results are public

Claims other papers make about this one

Each is checked against this paper’s own text: does it say what the citing paper says it does?

  • Estimates of DGEs are less precise than estimates of population effects because they only use within-family genotype variation, which is half of the variation in a random-mating population.

    “Furthermore, even when parental genotypes are available, estimates of DGEs are less precise than estimates of population effects because they only use within-family genotype variation, which is half of the variation in a random-mating population3.”

    From Family-GWAS reveals effects of environment and mating on genetic associations (2024) · cites this paper as [3]

    The Main section says that random segregation within families enables unbiased estimates of direct genetic effects, and the Results section on genome-wide analyses says power to analyze individual-SNP direct effects was limited (page numbers unavailable). But the paper does not state the claimed comparison with population-effect precision or that within-family genotype variation is half of population variation. Its reported one-half effective sample size compares parent–offspring pairs with complete parental-genotype data, not direct-effect estimates with population-effect estimates. (LLM, conf=0.82)

    Experimental — accuracy not yet measured
  • Although population structure can introduce bias into imputed parental genotypes, the resulting bias in DGE estimates is negligible for samples of relatively homogeneous genetic ancestry typically used in GWAS.

    “Although population structure can introduce bias into imputed parental genotypes, the resulting bias in DGE estimates is negligible for the samples of relatively homogeneous genetic ancestry typically used in GWAS3,20.”

    From Family-GWAS reveals effects of environment and mating on genetic associations (2024) · cites this paper as [3]

    The Methods section, “Effect of population structure” (page not available), says population structure biases imputed parental genotypes and can bias direct-effect estimates. It says that for relatively homogeneous samples this bias is *likely* negligible at the individual-SNP level; the claim states negligibility without that uncertainty or scope limitation. The passage also does not establish that such samples are typically used in GWAS. (LLM, conf=0.85)

    Experimental — accuracy not yet measured
  • Because segregation of chromosomes during meiosis is independent of environment, estimates of DGEs are free from confounding due to gene-environment correlation, including from IGEs and population-stratification.

    “Because the segregation of chromosomes during meiosis is independent of environment, estimates of DGEs are free from confounding due to gene-environment correlation, including from IGEs and population-stratification3,9,15,16.”

    From Family-GWAS reveals effects of environment and mating on genetic associations (2024) · cites this paper as [3]

    The Main section says random segregation during meiosis enables parent–offspring trio analyses to obtain unbiased, stratification-free estimates of direct genetic effects. But it also warns that sibling-based estimates can be biased by indirect genetic effects from siblings; the unqualified claim that estimates of direct genetic effects are free of this confounding is too broad. (LLM, conf=0.84)

    Experimental — accuracy not yet measured
  • LDSC was designed for GWAS summary statistics on population effects derived from samples of unrelated individuals.

    “Since LDSC was designed for GWAS summary statistics on population effects derived from samples of unrelated individuals, we use the effective sample size for the required sample size input3,9.”

    From Family-GWAS reveals effects of environment and mating on genetic associations (2024) · cites this paper as [3]

    The Results section describes population effects as estimates from standard GWAS, and Methods > GWASs in unrelated individuals describes an unrelated GWAS sample. The paper does not state that LDSC was designed for summary statistics from such samples; Methods > IBD inference mentions LD Score Regression software only as a way to calculate LD scores. No page numbers are available in our copy. (LLM, conf=0.82)

    Experimental — accuracy not yet measured
  • FGWAS uses the natural experiment of random assignment of genotype within families to separate out the contribution of direct genetic effects (DGEs) from other factors contributing to genetic associations.

    “Family-based GWAS (FGWAS) uses the natural experiment of random assignment of genotype within families to separate out the contribution of direct genetic effects (DGEs) — causal effects of alleles in an individual on an individual — from other factors contributing to genetic associations.”

Self-check

Percentages, sums and reported differences re-computed from the paper’s own numbers.

The self-check needs the paper’s full text.

Cites

Flagged first, then by citations

Research works printed but not linked: 33

  1. [1]Mills MC Rahal C A scientometric review of genome-wide association studies Commun. Biol. 2019 2 1 11 10.1038/s42003-018-0261-x 30623105 PMC6323052 · DOI
  2. [2]Meuwissen TH Hayes BJ Goddard ME Prediction of total genetic value using genome-wide dense marker maps Genetics 2001 157 1819 1829 10.1093/genetics/157.4.1819 11290733 PMC1461589 · DOI · Paper page
  3. [4]Kong A The nature of nurture: effects of parental genotypes Science 2018 359 424 428 10.1126/science.aan6877 29371463 · DOI · Paper page
  4. [5]Price AL Principal components analysis corrects for stratification in genome-wide association studies Nat. Genet. 2006 38 904 909 10.1038/ng1847 16862161 · DOI · Paper page
  5. [7]Berg JJ Reduced signal for polygenic adaptation of height in UK Biobank Elife 2019 8 1 47 10.7554/eLife.39725 PMC6428572 30895923 · DOI · Paper page
  6. [8]Zaidi AA Mathieson I Demographic history mediates the effect of stratification on polygenic scores Elife 2020 9 e61548 10.7554/eLife.61548 33200985 PMC7758063 · DOI · Paper page
  7. [10]Brumpton B Avoiding dynastic, assortative mating, and population stratification biases in Mendelian randomization through within-family analyses Nat. Commun. 2020 11 1 13 10.1038/s41467-020-17117-4 32665587 PMC7360778 · DOI · Paper page
  8. [11]Sohail M Polygenic adaptation on height is overestimated due to uncorrected stratification in genome-wide association studies Elife 2019 8 8 11 10.7554/eLife.39702 PMC6428571 30895926 · DOI · Paper page
  9. [13]Visscher PM Assumption-free estimation of heritability from genome-wide identity-by-descent sharing between full siblings PLoS Genet. 2006 2 e41 10.1371/journal.pgen.0020041 16565746 PMC1413498 · DOI
  10. [14]Benyamin B Family-based genome-wide association studies Pharmacogenomics 2009 10 181 190 10.2217/14622416.10.2.181 19207019 · DOI
  11. [15]Weiner DJ Polygenic transmission disequilibrium confirms that common and rare variation act additively to create risk for autism spectrum disorders Nat. Genet. 2017 49 978 985 10.1038/ng.3863 28504703 PMC5552240 · DOI · Paper page
  12. [17]Mostafavi H Variable prediction accuracy of polygenic scores within an ancestry group Elife 2020 9 e48376 10.7554/eLife.48376 31999256 PMC7067566 · DOI · Paper page
  13. [18]Fulker DW Cherny SS Sham PC Hewitt JK Combined linkage and association sib-pair analysis for quantitative traits Am. J. Hum. Genet. 1999 64 259 267 10.1086/302193 9915965 PMC1377724 · DOI · Paper page
  14. [19]Hickey JM Kinghorn BP Tier B van der Werf JHJ Cleveland MA A phasing and imputation method for pedigreed populations that results in a single-stage genomic evaluation Genet. Sel. Evol. 2012 44 9 10.1186/1297-9686-44-9 22462519 PMC3378429 · DOI
  15. [20]Sargolzaei M Chesnais JP Schenkel FS A new approach for efficient genotype imputation using information from relatives BMC Genomics 2014 15 478 10.1186/1471-2164-15-478 24935670 PMC4076979 · DOI
  16. [21]Whalen A Gorjanc G Hickey JM AlphaFamImpute: high-accuracy imputation in full-sib families from genotype-by-sequencing data Bioinformatics 2020 36 4369 4371 10.1093/bioinformatics/btaa499 32467963 PMC7520044 · DOI
  17. [22]Bycroft C The UK Biobank resource with deep phenotyping and genomic data Nature 2018 562 203 209 10.1038/s41586-018-0579-z 30305743 PMC6786975 · DOI · Paper page
  18. [23]Manichaikul A Robust relationship inference in genome-wide association studies Bioinformatics 2010 26 2867 2873 10.1093/bioinformatics/btq559 20926424 PMC3025716 · DOI · Paper page
  19. [24]Cavalli-Sforza LL Feldman MW Cultural versus biological inheritance: phenotypic transmission from parents to children (a theory of the effect of parental phenotypes on children’s phenotypes) Am. J. Hum. Genet. 1973 25 618 637 4797967 PMC1762580
  20. [25]Rice J Cloninger CR Reich T Multifactorial inheritance with cultural transmission and assortative mating. I. Description and basic properties of the unitary models Am. J. Hum. Genet. 1978 30 618 643 747189 PMC1685878
  21. [27]Cook JP Mahajan A Morris AP Fine-scale population structure in the UK Biobank: implications for genome-wide association studies Hum. Mol. Genet. 2020 29 2803 2811 10.1093/hmg/ddaa157 32691046 · DOI
  22. [29]Kruglyak L Daly MJ Reeve-Daly MP Lander ES Parametric and nonparametric linkage analysis: a unified multipoint approach Am. J. Hum. Genet. 1996 58 1347 8651312 PMC1915045
  23. [31]Gudbjartsson DF A frameshift deletion in the sarcomere gene MYL4 causes early-onset familial atrial fibrillation. Eur. Heart J. 2016 38 27 34 10.1093/eurheartj/ehw379 27742809 · DOI
  24. [32]Rubin, D. B. Multiple Imputation for Nonresponse in Surveys (John Wiley & Sons, 2004).
  25. [33]Munafò MR Tilling K Taylor AE Evans DM Davey Smith G Collider scope: when selection bias can substantially influence observed associations Int. J. Epidemiol. 2017 47 226 235 10.1093/ije/dyx206 PMC5837306 29040562 · DOI · Paper page
  26. [35]Benonisdottir, S. & Kong, A. The genetics of participation: method and analysis. Preprint in bioRxiv10.1101/2022.02.11.480067 (2022).
  27. [38]Border, R. et al. Cross-trait assortative mating is widespread and inflates genetic correlation estimates. Preprint at bioRxiv10.1101/2022.03.21.485215 (2022). PMC9901291 36395242
  28. [39]Hwang L-D Estimating indirect parental genetic effects on offspring phenotypes using virtual parental genotypes derived from sibling and half sibling pairs. PLoS Genet. 2020 16 e1009154 10.1371/journal.pgen.1009154 33104719 PMC7646364 · DOI
  29. [41]Nagylaki T Assortative mating for a quantitative character J. Math. Biol. 1982 16 57 74 10.1007/BF00275161 7161578 · DOI
  30. [42]Balbona JV Kim Y Keller MC Estimation of parental effects using polygenic scores Behav. Genet. 2021 51 264 278 10.1007/s10519-020-10032-w 33387133 PMC8093180 · DOI
  31. [44]McCarthy S A reference panel of 64,976 haplotypes for genotype imputation Nat. Genet. 2016 48 1279 1283 10.1038/ng.3643 27548312 PMC5388176 · DOI · Paper page
  32. [47]Vilhjálmsson BJ Modeling linkage disequilibrium increases accuracy of polygenic risk scores Am. J. Hum. Genet. 2015 97 576 592 10.1016/j.ajhg.2015.09.001 26430803 PMC4596916 · DOI · Paper page
  33. [48]Loh P-R Efficient Bayesian mixed-model analysis increases association power in large cohorts Nat. Genet. 2015 47 284 290 10.1038/ng.3190 25642633 PMC4342297 · DOI · Paper page

Cited by

Sorted by citations

Authors
Young, A. I., Nehzati, S. M., Benonisdottir, S., Okbay, A., Jayashankar, H., Lee, C., Cesarini, D., Benjamin, D. J., Turley, P., & Kong, A.
Affiliations
  • University of Southern California (US)
  • Vrije Universiteit Amsterdam (NL)
  • National Bureau of Economic Research (US)
  • University of Oxford (GB)
Corresponding author
  • Young, Alexander I.
  • Kong, Augustine
Venue
Nature Genetics
Published
2022-06
Volume
54
Issue
6
Pages
897-905
ISSN
  • 1061-4036 (ISSN-L)
  • 1546-1718
OpenAlex
W4281917178
Cited by
144
References
0
Citations by year
2026: 32 · 2025: 30 · 2024: 29 · 2023: 19 · 2022: 19 · 2021: 12 · 2020: 3
FWCI
14.62
Type
journal article
Language
en
Open access
Locations
  • Nature Genetics (publishedVersion, cc-by)
  • bioRxiv (Cold Spring Harbor Laboratory) (acceptedVersion, cc-by-nc)
  • PubMed Central (publishedVersion, cc-by)
  • PubMed (publishedVersion)
  • VU Research Portal (publishedVersion, other-oa)
  • eScholarship (California Digital Library) (submittedVersion, cc-by)
  • Data Archiving and Networked Services (DANS) (submittedVersion, other-oa)
Funders
  • Li Ka Shing Foundation
  • Netherlands Organisation for Scientific Research | Nationaal Regieorgaan Praktijkgericht Onderzoek SIA (016.Veni.198.058)
  • Foundation for the National Institutes of Health (R24-AG065184, R01-AG042568, K99-AG062787-01)
Topics
  • Life Sciences › Biochemistry, Genetics and Molecular Biology › Genetics › Genetic Associations and Epidemiology
  • Life Sciences › Biochemistry, Genetics and Molecular Biology › Genetics › Genetic and phenotypic traits in livestock
  • Life Sciences › Biochemistry, Genetics and Molecular Biology › Genetics › Genetic Mapping and Diversity in Plants and Animals
Keywords
Mendelian randomization, direct genetic effects, indirect genetic effects, confounding, educational attainment, cognitive ability
MeSH terms
Genotype, Humans, Parents, Software, Polymorphism, Single Nucleotide, Genome-Wide Association Study

About these checks

Registry, retraction, Cites and statistics results are computed, not judged, with one exception: a printed reference that no index lists is read once by a model before it counts as a citation error. Claims are extracted by a model and kept only when their quote appears verbatim in the text; each cited source is then checked, with a model’s help, for whether it says what the claim relies on. Those verdicts can be wrong — read the claim’s quote and the cited source before relying on one.

Statistics and arithmetic are recomputed, not judged.

Confidence bands are shown only for fields calibrated on our annotated test set. None is calibrated yet, so no band is shown.

What counts as a claim ›How we measure confidence ›How we check this ›

Are you an author, or do you see a problem with this page? Request a review