 Proceedings
 Open Access
 Published:
Mendelian randomization in family data
BMC Proceedings volume 3, Article number: S45 (2009)
Abstract
The phrase "mendelian randomization" has become associated with the use of genetic polymorphisms to uncover causal relationships between phenotypic variables. The statistical methods useful in mendelian randomization are known as instrumental variable techniques. We present an approach to instrumental variable estimation that is useful in family data and is robust to the use of weak instruments. We illustrate our method to measure the causal influence of lowdensity lipoprotein on highdensity lipoprotein, body mass index, triglycerides, and systolic blood pressure. We use the Framingham Heart Study data as distributed to participants in the Genetics Analysis Workshop 16.
Background
In epidemiological studies, establishing and measuring causal relationships is of primary importance. Unfortunately, randomization, the most important tool for unraveling causal relationships, is not generally available. Despite advances in study design and statistical adjustment, the possibility of confounding and reverse causation continues to be problematic. In recent years it has been suggested that nature itself has already performed a set of randomized experiments by assigning genes according to Mendel's laws [1]. These genes affect the function or expression levels of specific gene products, that in a cascade of cause and effect, eventually lead to human disease. By utilizing the statistical concept of an instrumental variable (IV), it may be possible to use genetics to solve some of the problems that have plagued epidemiology for decades. The goal of this approach, known as "mendelian randomization," is not to detect genetic factors of disease, but rather to use genetic factors of disease to uncover the causal relationships between phenotypes.
Figure 1a shows a graphical depiction of a situation in which IV methods might be useful. Suppose that we wished to assess the relationship between X and Y. Because U is unmeasured, there is no way of estimating the strength of this relationship using ordinary epidemiological methodology such as linear regression. In simplistic terms, the logic behind MR is that G can only affect Y by affecting X. Therefore, an association between G and Y is best explained by a causal relationship between X and Y. In IV methods, we look only at how that part of X which is influenced by G affects Y. In order to do this we need to make at least the following three assumptions [2]:

1.
G is not independent of X.

2.
G is independent of U.

3.
G is independent of Y given U and X.
Only Assumption 1 is testable. Assumptions 2 and 3 must be accepted or rejected using subjectspecific knowledge. It is well known that there are many potential violations of these conditions when using MR. Such problems include: linkage disequilibrium, pleiotropy, population stratification, and canalization [2]. However, by careful selection of polymorphisms, many of these concerns may be minimized. A well done MR study has assumptions that cannot be tested, but these assumptions are easier to believe than those involved in a direct assessment of the relationship between X and Y.
IV techniques are well known tools in econometrics, and there is an extensive body of literature discussing the properties of various estimators. Popular methods of estimation include twostage leastsquares regression and limited information maximumlikelihood estimation [3]. These methods generally impose the assumption of linearity of the effects. Traditionally, a Wald type confidence interval of the form "estimate ± error" is used. However, a growing body of literature has shown that if the instrument is weak, this may result in confidence intervals of incorrect size [4]. A weak instrument is one which violates or comes close to violating Assumption 1 above. One approach to solving this problem is to invert tests that are robust to weak instruments [5]. We have adopted this approach below. The motivation for our approach is similar to that of Gestimation.
With few exceptions to date, discussions of MR have focused on population data. However, there is a large amount of existing family data, and there are well known advantages to familybased studies. For example, MR may provide the ability to check for mendelian errors, and it may provide protection against population stratification [6]. In this paper we suggest an approach to MR which is broadly applicable to family data.
We apply this approach to measure the causal relationship between lowdensity lipoprotein levels (LDL) and the variables highdensity lipoprotein (HDL), triglycerides (TG), body mass index (BMI), and systolic blood pressure (SBP) in the Framingham Heart Study data as distributed to participants in the Genetic Analysis Workshop 16 (GAW16) as Problem 2. For our IV we use singlenucleotide polymorphisms (SNPs) from the 50 k data set, which are in or near the LDLR (OMIM 606945) and APOB genes (OMIM 107730), because they have direct influences on LDL. We make some additional comments in justification of our approach in the discussion.
Methods
Data description
We utilized the 50 k genotype data for all three cohorts of the Framingham Heart Study data. Data were obtained and used in compliance with the data use agreement and the Case Western Reserve University Institutional Review Board. We found only one SNP (rs2738457) in the LDLR gene in this data set. To select SNPs from APOB we compared the available SNPs to those analyzed in an independent sample by Benn et al. [7]. There were five SNPs shared in common between the two studies. We chose the two SNPs from Benn et al. [7] that had a pvalue less than 0.001, appeared to be acting in an additive manner, and were shared in common with the Framingham Heart Study 50 k SNP data. These SNPs were rs1042031 and rs679899.
LDL was calculated using the Friedewald equation. In order to deal with possible heterogeneity of effect and make use of the multiple visits, each of the phenotypic observations was stratified by age at examination. The following strata were used: 029, 3044, 4560, and over 60 years of age. An approximate year of birth for each individual was calculated as the mean difference between the approximate exam date and the age at exam. All phenotype variables besides age and approximate year of birth were logtransformed for analysis in ASSOC (S.A.G.E. v5.4.1). In all analysis, the mean centered year of birth, age, and sex were used as covariates. SBP was adjusted by adding 10 to those on treatment [8].
Statistical method
Suppose there are K pedigrees and n_{ k }individuals in each pedigree. Suppose also that x_{ k }and y_{ k }are n_{ k }× 1 vectors representing the variables X and Y in Figure 1. Also, G_{ k }is a n_{ k }× t matrix representing coded genotypes, and A_{ k }and B_{ k }represent matrices of covariates for x_{ k }and y_{ k }respectively. Let and let y and G be defined analogously. We assume the following two linear structural equations:
Here, e_{ xk }and e_{ yk }are (possibly) correlated error terms. A variancecomponent model is assumed to describe the covariance structure . The variancecomponents model is described in some more detail below.
Note that β_{ xy }is the parameter that represents the causal effect size. Suppose that we wish to test the hypothesis H_{0} : β_{ xy }= β_{ H }. The basic approach we suggest here is to form a new trait
Combining Eq. (3) with Eqs. (2) and (1), we obtain
for suitably defined , and e_{ rk }. Under the model assumptions, r(β_{ H }) is independent of G given the covariates if and only if β_{ xy }= β_{ H }. Hence, we may test the hypothesis H_{0} : β_{ xy }= β_{ H }using some form of a familybased genetic association test. The test may be recomputed along an entire grid of values. Those values that are rejected at some alpha level lie outside the confidence interval. If the association between G and x is not statistically significant, then in most cases the confidence interval will cover the entire real line. For finite samples it is also possible that the confidence interval is the null set.
We have chosen to use the method implemented in the program ASSOC of S.A.G.E. (v5.4.1) to test for genetic association. This familybased test of association between markers and a continuous phenotype developed by George and Elston [9] allows for familial correlations by simultaneously estimating residual and multifactorial (polygenic, familial and marital) variance components. It is assumed that the data can be transformed via the GeorgeElston transformation into a multivariate distribution. ASSOC maximizes the likelihood conditional on the genotype values given the assumptions above. We supplied guesses as to the beginning and end points of the confidence interval for β_{ xy }. We then ran ASSOC iteratively over a grid of values for β_{ H }between the end points. Values for which the maximization algorithm clearly did not properly converge were not used. We used cubic spline interpolation as implemented in the R function "splinefun" to find the boundaries on the confidence interval and the point estimate. As a point estimate we used the maximum pvalue.
We also used ASSOC to assess the causal relationships directly. We have labelled this method as "regression" below, although it is not strictly a linear regression because the familial correlation is accounted for using the same variancecomponent model described above. The R^{2} values reported below are calculated as where and . Here N_{1} and N_{2} are the total number of nonmissing observations on x and G, respectively.
Results
As mentioned in the Introduction, we chose to look at the effect of LDL on a number of other variables (see Figure 1b). The first step in this analysis is to look at the effect of the chosen SNPs on LDL levels. There is good evidence for genetic association for all four age groups (Table 1). However, the chosen SNPs explain only a small proportion of the overall genetic variance.
There is little that can be said based on the IV analysis alone. In all cases the 95% confidence intervals are quite large. Only 1 out of 16 intervals exclude 0. It is interesting to note that in a greaterthanexpected number of cases, the signs of the IV estimates and the regression estimates are the same (Table 2). The observed negative correlation between LDL and HDL is to be expected because HDL is thought to help rid the body of LDL. However, if this is the case, the causal arrow runs from HDL to LDL. The regression estimates relating HDL to LDL and TG to LDL appear to change with age. However, the regression estimates should not be overanalyzed because of potential confounders and reverse causation.
Discussion
It may be difficult to define exactly what is meant by "causal effect" in a system that is evolving dynamically. Onetime interventions in such a system will have different effect sizes depending on the time since intervention. Eventually the system will settle back into an equilibrium state. We take the causal meaning of our results to reflect more upon the equilibrium state of the system than the instantaneous response to intervention. For example, our estimate of the causal effect of log [LDL] on log [TG] is 1.770 in the 029 age group. Essentially we interpret this to mean that if we could raise an individual's log [LDL] by one log [mg/dL], eventually the system would settle into a new state with log [TG] raised by about 1.770 log [mg/dL].
Every MR study should be accompanied by a biological argument that certain assumptions are met. Our argument is that polymorphisms on or near the APOB and LDLR genes will probably only affect the other phenotypes through their effect on LDL. This is because the APOB protein is one of the principal components of LDL. Similarly, LDLR is known to bind to LDL, allowing for endocytosis. Of course, the above claim could be disputed. The purpose of this study is mainly to demonstrate a statistical approach, and we do not necessarily claim that the results have a strong biological justification for Assumptions 23.
One concern that we believe is particularly relevant is whether the chosen markers are in linkage disequilibrium with polymorphisms that modify the function or the plasma levels of the LDL protein. For example, the causal effect size of LDL on HDL could differ depending on certain polymorphism in the APOB gene. In this case, Assumption 3 would be violated. While this would invalidate our estimates of effect size, it may not invalidate the qualitative conclusion. That is, if the confidence interval does not include 0, we still have evidence that LDL is causally related to HDL.
With regard to the selection of SNPs as instruments, we have several comments. First, we caution that in this type of study, the polymorphisms used should be justified a priori. Searching through a large number of possible instruments for statistical significance then performing IV analysis within the same data is not an acceptable procedure. Second, the selected polymorphisms should have biological justification as an IV. That is, it should meet Assumptions 13 mentioned in the introduction. Third, if there are no polymorphisms that explain a reasonably high percentage of the variance in the cause, then IVs are unlikely to be useful. Our IV estimates have very large confidence intervals. We believe this is ultimately a result of instruments that are weak predictors of LDL levels.
Conclusion
In this paper we suggest a method for performing MR in family data. While we have chosen to use the program ASSOC to implement this method, nearly any statistical test for association between a quantitative trait and a set of polymorphisms could be used. We used ASSOC because it allows for familial correlations while using information all individuals in the family. If concerns about population stratification are raised, a more robust test such as the transmissiondisequilibrium test could be used. Our approach is robust to the problem of weak instruments in the sense that it should maintain correct coverage rates. However, it will still have low power. Perhaps the biggest limitation to MR is the need for polymorphisms that explain a large percentage of the variation in a given trait.
Abbreviations
 BMI:

Body mass index
 GAW16:

Genetic Analysis Workshop 16
 HDL:

Highdensity lipoprotein
 IV:

Instrumental variable
 LDL:

Lowdensity lipoprotein
 MR:

Mendelian randomization
 SBP:

Systolic blood pressure
 SNP(s):

Singlenucleotide polymorphism
 TG:

Triglycerides.
References
 1.
Davey Smith G, Ebrahim S: 'Mendelian randomization': can genetic epidemiology contribute to understanding environmental determinants of disease?. Int J Epidemiol. 2003, 32: 122. 10.1093/ije/dyg070.
 2.
Didelez V, Sheehan N: Mendelian randomization as an instrumental variable approach to causal inference. Stat Methods Med Res. 2007, 16: 309330. 10.1177/0962280206077743.
 3.
Bowden RJ, Turkington DA: Instrumental Variables. 1984, New York, Cambridge University Press
 4.
Staiger D, Stock JH: Instrumental variables regression with weak instruments. Econometrica. 1997, 65: 557586. 10.2307/2171753.
 5.
Andrews DWK, Stock JH: Inference with Weak Instruments. 2005, NBER Technical Working Papers 0313, National Bureau of Economic Research, Inc
 6.
Visscher PM, Andrew T, Nyholt DR: Genomewide association studies of quantitative traits with related individuals: little (power) lost but much to be gained. Eur J Hum Genet. 2008, 16: 387390. 10.1038/sj.ejhg.5201990.
 7.
Benn M, Stene MCA, Nordestgaard BG, Jensen GB, Steffensen R, TybjaergHansen A: Common and rare alleles in apolipoprotein B contribute to plasma levels of lowdensity lipoprotein cholesterol in the general population. J Clin Endocrinol Metab. 2008, 93: 10381045. 10.1210/jc.20071365.
 8.
Cui JS, Hopper JL, Harrap SB: Antihypertensive treatments obscure familial contributions to blood pressure variation. Hypertension. 2003, 41: 207210. 10.1161/01.HYP.0000044938.94050.E3.
 9.
George VT, Elston RC: Testing the association between polymorphic markers and quantitative traits in pedigrees. Genet Epidemiol. 1987, 4: 193201. 10.1002/gepi.1370040304.
Acknowledgements
The Genetic Analysis Workshops are supported by NIH grant R01 GM031575 from the National Institute of General Medical Sciences. This work is supported by the National Center for Research Resources (NCRR) Human Genetic Analysis Resource (RR03655) (CGM), the NCRR Multidisciplinary Clinical Research Career Development Programs Grant (KL2RR024990) (CMS), and National Heart Lung and Blood Institute grant (HL0756) (NJM). Some of the results were obtained using the package S.A.G.E., which is supported by a grant (RR03655) from the NCRR.
This article has been published as part of BMC Proceedings Volume 3 Supplement 7, 2009: Genetic Analysis Workshop 16. The full contents of the supplement are available online at http://www.biomedcentral.com/17536561/3?issue=S7.
Author information
Affiliations
Corresponding author
Additional information
Competing interests
The authors declare that they have no competing interests.
Authors' contributions
NJM conceived of the statistical method, analyzed the data, and drafted the manuscript. CGM helped select the traits and advised on the epidemiological aspects of the study. CMS contributed to the statistical development and writing of the paper.
Rights and permissions
This article is published under license to BioMed Central Ltd. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
About this article
Cite this article
Morris, N.J., GrayMcGuire, C. & Stein, C.M. Mendelian randomization in family data. BMC Proc 3, S45 (2009). https://doi.org/10.1186/175365613S7S45
Published:
Keywords
 Instrumental Variable
 Population Stratification
 Mendelian Randomization
 Family Data
 Genetic Analysis Workshop