# Mendelian randomization

In epidemiology, Mendelian randomization (commonly abbreviated to MR) is a method using measured variation in genes to interrogate the causal effect of an exposure on an outcome. Under key assumptions (see below), the design reduces both reverse causation and confounding, which often substantially impede or mislead the interpretation of results from epidemiological studies.

Gregor Mendel. The term Mendelian randomization was termed because the random assignment of genetic variants from parents to offspring is fundamental to the method.

The study design was first proposed in 1986 [1] and subsequently described by Gray and Wheatley[2] as a method for obtaining unbiased estimates of the effects of a putative causal variable without conducting a traditional randomized controlled trial (i.e. the "gold standard" in epidemiology for establishing causality). These authors also coined the term Mendelian randomization.

## Motivation

One of the predominant aims of epidemiology is to identify modifiable causes of health outcomes and disease especially those of public health concern. In order to ascertain whether modifying a particular trait (e.g. via an intervention, treatment or policy change) will convey a beneficial effect within a population, firm evidence that this trait causes the outcome of interest is required. However, many observational epidemiological study designs are limited in the ability to discern correlation from causation - specifically whether a particular trait causes an outcome of interest, is simply related to that outcome (but does not cause it) or is a consequence of the outcome itself. Only the former will be beneficial within a public health setting where the aim is to modify that trait to reduce the burden of disease. There are many epidemiological study designs that aim to understand relationships between traits within a population sample, each with shared and unique advantages and limitations in terms of providing causal evidence, with the "gold standard" being randomized controlled trials.

Well-known successful demonstrations of causal evidence consistent across multiple studies with different designs include the identified causal links between smoking and lung cancer, and between blood pressure and stroke. However, there have also been notable failures when exposures hypothesized to be a causal risk factor for a particular outcome were later shown by well conducted randomized controlled trials not to be causal. For instance, it was previously thought that hormone replacement therapy would prevent cardiovascular disease, but it is now known to have no such benefit and may even adversely affect health.[3] Another notable example is that of selenium and prostate cancer. Some observational studies found an association between higher circulating selenium levels (usually acquired through various foods and dietary supplements ) and lower risk of prostate cancer. However, the Selenium and Vitamin E Cancer Prevention Trial (SELECT) showed evidence that dietary selenium supplementation actually increased the risk of prostate and advanced prostate cancer and had an additional off-target effect on increasing type 2 diabetes risk.[4]

Such inconsistencies between observational epidemiological studies and randomized controlled trials are likely a function of social, behavioral, or physiological confounding factors in many observational epidemiological designs, which are particularly difficult to measure accurately and difficult to control for. Moreover, randomized controlled trials are usually expensive, time-consuming and laborious and many epidemiological findings cannot be ethically replicated in clinical trials.

## Definition

Mendelian randomization (MR) is fundamentally an Instrumental variables estimation method hailing from Econometrics. The method uses the properties of germline genetic variation (usually in the form of single nucleotide polymorphisms or SNPs) strongly associated with a putative exposure as a "proxy" or "instrument" for that exposure to test for and estimate a causal effect of the exposure on an outcome of interest from observational data. The genetic variation used will have either well-understood effects on exposure patterns (e.g. propensity to smoke heavily) or effects that mimic those produced by modifiable exposures (e.g., raised blood cholesterol[1]). Importantly, the genotype must only affect the disease status indirectly via its effect on the exposure of interest.[5]

Directed acyclic graph traditionally used to represent the Mendelian randomization framework and its core assumptions. ${\displaystyle Z}$ is the genetic variants, ${\displaystyle X}$ is the exposure, ${\displaystyle Y}$ is the outcome of interest, and ${\displaystyle U}$ are possible confounders.

As genotypes are assigned randomly when passed from parents to offspring during meiosis, then groups of individuals defined by genetic variation associated with an exposure at a population level should be largely unrelated to the confounding factors that typically plague observational epidemiology studies. Germline genetic variation (i.e. that which can be inherited) is also temporarily fixed at conception and not modified by the onset of any outcome or disease, precluding reverse causation. Additionally, given improvements in modern genotyping technologies, measurement error and systematic misclassification is often low with genetic data. In this regard Mendelian randomization can be thought of as analogous to "nature's randomized controlled trial".

Mendelian randomization requires three core instrumental variable assumptions.[6] Namely:

1. The genetic variant(s) being used as an instrument for the exposure is associated with the exposure. This is known as the "relevance" assumption.
2. There are no common causes (i.e. confounders) of the genetic variant(s) and the outcome of interest. This is known as the "independence" or "exchangeability" assumption.
3. There is no independent pathway between the genetic variant(s) and the outcome other than through the exposure. This is known as the "exclusion restriction" or "no horizontal pleiotropy" assumption.

To ensure that the first core assumption is validated, Mendelian randomization requires characterized associations between genetic variation and exposures of interest. These are usually sourced from genome-wide association studies though can also be candidate gene studies. The second assumption relies on there being no population substructure (e.g. geographical factors that induce an association between the genotype and outcome), mate choice that is not associated with genotype (i.e. random mating or panmixia) and no dynastic effects (i.e. where the expression of parental genotype in the parental phenotype directly affects the offspring phenotype).

## Statistical analysis

Mendelian randomization is usually applied through the use of instrumental variables estimation with genetic variants acting as instruments for the exposure of interest.[7] This can be implemented using data on the genetic variants, exposure and outcome of interest for a set of individuals in a single dataset or using summary data on the association between the genetic variants and the exposure and the association between the genetic variants and the outcome in separate datasets. The method has also been used in economic research studying the effects of obesity on earnings, and other labor market outcomes.[8]

When a single dataset is used the methods of estimation applied are those frequently used elsewhere in instrumental variable estimation, such as two-stage least squares.[9] If multiple genetic variants are associated with the exposure they can either be used individually as instruments or combined to create an allele score which is used as a single instrument.

Analysis using summary data often applies data from genome-wide association studies. In this case the association between genetic variants and the exposure is taken from the summary results produced by a genome-wide association study for the exposure. The association between the same genetic variants and the outcome is then taken from the summary results produced by a genome-wide association study for the outcome. These two sets of summary results are then used to obtain the MR estimate. Given the following notation:

${\displaystyle {\hat {\pi }}_{g}=}$ effect of genetic variant ${\displaystyle g}$ on the exposure ${\displaystyle (X)}$;

${\displaystyle {\hat {\Gamma }}_{g}=}$ estimated effect of genetic variant g on the outcome ${\displaystyle (Y)}$;

${\displaystyle {\hat {\sigma }}_{g}=}$ estimated standard error of this estimated effect;

${\displaystyle {\hat {\beta }}_{MR}=}$ MR estimate of the causal effect of the exposure ${\displaystyle X}$ on the outcome ${\displaystyle Y}$;

and considering the effect of a single genetic variant, the MR estimate can be obtained from the Wald ratio:

${\displaystyle {\hat {\beta }}_{MR}={\frac {{\hat {\Gamma }}_{g}}{{\hat {\pi }}_{g}}}}$

When multiple genetic variants are used, the individual ratios for each genetic variants are combined using inverse variance weighting where each individual ratio is weighted by the uncertainty in their estimation.[10] This gives the IVW estimate which can be calculated as:

${\displaystyle {\hat {\beta }}_{IVW}={\frac {\Sigma _{g=1}^{G}{\hat {\pi }}_{g}{\hat {\Gamma }}_{g}\sigma _{y,g}^{2}}{\Sigma _{g=1}^{G}{\hat {\pi }}_{g}^{2}\sigma _{y,g}^{2}}}}$

Alternatively, the same estimate can be obtained from a linear regression which used the genetic variant-outcome association as the outcome and the genetic variant-exposure association as the exposure. This linear regression is weighted by the uncertainty in the genetic-variant outcome association and does not include a constant.

${\displaystyle {\hat {\Gamma }}_{g}=\beta _{IVW}{\hat {\pi }}_{g}+u_{g}\qquad weighted\ by\ 1/{\hat {\sigma }}_{y,g}^{2}}$

These methods only provide reliable estimates of the causal effect of the exposure on the outcome under the core instrumental variable assumptions. Alternative methods are available that are robust to a violation of the third assumption, i.e. that provide reliable results under some types of horizontal pleiotropy.[11] Additionally some biases that arise from violations of the second IV assumption, such as dynastic effects, can be overcome through the use of data which includes siblings or parents and their offspring.[12]

## History

The Mendelian randomization method depends on two principles derived from the original work by Gregor Mendel on genetic inheritance. Its foundation come from Mendel’s laws namely 1) the law of segregation in which there is complete segregation of the two allelomorphs in equal number of germ-cells of a heterozygote and 2) separate pairs of allelomorphs segregate independently of one another and which were first published as such in 1906 by Robert Heath Lock. Another progenitor of Mendelian randomization is Sewall Wright who introduced path analysis, a form of causal diagram used for making causal inference from non-experimental data. The method relies on causal anchors, and the anchors in the majority of his examples were provided by Mendelian inheritance, as is the basis of MR.[13] Another component of the logic of MR is the instrumental gene, the concept of which was introduced by Thomas Hunt Morgan.[14] This is important as it removed the need to understand the physiology of the gene for making the inference about genetic processes.[citation needed]

Since that time the literature includes examples of research using molecular genetics to make inference about modifiable risk factors, which is the essence of MR. One example is the work of Gerry Lower and colleagues in 1979 who used the N-acetyltransferase phenotype as an anchor to draw inference about various exposures including smoking and amine dyes as risk factors for bladder cancer. [15] Another example is the work of Martijn Katan (then of Wageningen University & Research, Netherlands) in which he advocated a study design using Apolipoprotein E allele as an instrumental variable anchor to study the observed relationship between low blood cholesterol levels and increased risk of cancer.[1] In fact, the term “Mendelian randomization” was first used in print by Richard Gray and Keith Wheatley (both of Radcliffe Infirmary, Oxford, UK) in 1991 in a somewhat different context; in a method allowing instrumental variable estimation but in relation to an approach relying on Mendelian inheritance rather than genotype.[2] In their 2003 paper, Shah Ebrahim and George Davey Smith use the term again to describe the method of using germline genetic variants for understanding causality in an instrumental variable analysis, and it is this methodology that is now widely used and to which the meaning is ascribed.[16] The Mendelian randomization method is now widely adopted in causal epidemiology, and the number of MR studies reported in the scientific literature has grown every year since the 2003 paper. In 2021 STROBE-MR guidelines were published to assist readers and reviewers of Mendelian randomization studies to evaluate the validity and utility of published studies.[17]

## References

1. ^ a b c Katan MB (March 1986). "Apolipoprotein E isoforms, serum cholesterol, and cancer". Lancet. 1 (8479): 507–8. doi:10.1016/s0140-6736(86)92972-7. PMID 2869248. S2CID 38327985.
2. ^ a b Gray R, Wheatley K (1991). "How to avoid bias when comparing bone marrow transplantation with chemotherapy". Bone Marrow Transplantation. 7 Suppl 3: 9–12. PMID 1855097.
3. ^ Rossouw JE, Anderson GL, Prentice RL, LaCroix AZ, Kooperberg C, Stefanick ML, et al. (July 2002). "Risks and benefits of estrogen plus progestin in healthy postmenopausal women: principal results From the Women's Health Initiative randomized controlled trial". JAMA. 288 (3): 321–33. doi:10.1001/jama.288.3.321. PMID 12117397.
4. ^ Klein EA, Thompson IM, Tangen CM, Crowley JJ, Lucia MS, Goodman PJ, et al. (October 2011). "Vitamin E and the risk of prostate cancer: the Selenium and Vitamin E Cancer Prevention Trial (SELECT)". JAMA. 306 (14): 1549–56. doi:10.1001/jama.2011.1437. PMC 4169010. PMID 21990298.
5. ^ Holmes MV, Ala-Korpela M, Smith GD (October 2017). "Mendelian randomization in cardiometabolic disease: challenges in evaluating causality". Nature Reviews. Cardiology. 14 (10): 577–590. doi:10.1038/nrcardio.2017.78. PMC 5600813. PMID 28569269.
6. ^ Wade, Kaitlin (2021). "MR Dictionary". MR Dictionary.{{cite web}}: CS1 maint: url-status (link)
7. ^ Didelez, Vanessa; Sheehan, Nuala (August 2007). "Mendelian randomization as an instrumental variable approach to causal inference". Statistical Methods in Medical Research. 16 (4): 309–330. doi:10.1177/0962280206077743. ISSN 0962-2802. S2CID 6236517.
8. ^ Böckerman P, Cawley J, Viinikainen J, Lehtimäki T, Rovio S, Seppälä I, et al. (January 2019). "The effect of weight on labor market outcomes: An application of genetic instrumental variables". Health Economics. 28 (1): 65–77. doi:10.1002/hec.3828. PMC 6585973. PMID 30240095.
9. ^ Wooldridge, Jeffrey M. (2010). Econometric analysis of cross section and panel data (2nd ed.). Cambridge, Mass.: MIT Press. ISBN 978-0-262-23258-6. OCLC 627701062.
10. ^ Burgess, Stephen; Butterworth, Adam; Thompson, Simon G. (2013). "Mendelian Randomization Analysis With Multiple Genetic Variants Using Summarized Data". Genetic Epidemiology. 37 (7): 658–665. doi:10.1002/gepi.21758. ISSN 0741-0395. PMC 4377079. PMID 24114802.
11. ^ Hemani, Gibran; Bowden, Jack; Davey Smith, George (2018-08-01). "Evaluating the potential role of pleiotropy in Mendelian randomization studies". Human Molecular Genetics. 27 (R2): R195–R208. doi:10.1093/hmg/ddy163. ISSN 0964-6906. PMC 6061876. PMID 29771313.
12. ^ Brumpton, Ben; Sanderson, Eleanor; Heilbron, Karl; Hartwig, Fernando Pires; Harrison, Sean; Vie, Gunnhild Åberge; Cho, Yoonsu; Howe, Laura D.; Hughes, Amanda; Boomsma, Dorret I.; Havdahl, Alexandra (2020). "Avoiding dynastic, assortative mating, and population stratification biases in Mendelian randomization through within-family analyses". Nature Communications. 11 (1): 3519. Bibcode:2020NatCo..11.3519B. doi:10.1038/s41467-020-17117-4. ISSN 2041-1723. PMC 7360778. PMID 32665587.
13. ^ Wright, S. (1921). "Correlation and causation". J. Agricultural Research. 20: 557–585.
14. ^ Morgan, T. H. (1917). "The Theory of the Gene". The American Naturalist. 51 (609): 513–544. doi:10.1086/279629. ISSN 0003-0147. JSTOR 2456204. S2CID 84050307.
15. ^ Lower, G. M.; Nilsson, T.; Nelson, C. E.; Wolf, H.; Gamsky, T. E.; Bryan, G. T. (1979). "N-Acetyltransferase phenotype and risk in urinary bladder cancer: approaches in molecular epidemiology. Preliminary results in Sweden and Denmark". Environmental Health Perspectives. 29: 71–79. doi:10.1289/ehp.792971. ISSN 0091-6765. PMC 1637362. PMID 510245.
16. ^ Smith, George Davey; Ebrahim, Shah (2003). "'Mendelian randomization': can genetic epidemiology contribute to understanding environmental determinants of disease?". International Journal of Epidemiology. 32 (1): 1–22. doi:10.1093/ije/dyg070. ISSN 0300-5771. PMID 12689998.
17. ^ Skrivankova, Veronika W.; Richmond, Rebecca C.; Woolf, Benjamin A. R.; Davies, Neil M.; Swanson, Sonja A.; VanderWeele, Tyler J.; Timpson, Nicholas J.; Higgins, Julian P. T.; Dimou, Niki; Langenberg, Claudia; Loder, Elizabeth W. (2021-10-26). "Strengthening the reporting of observational studies in epidemiology using mendelian randomisation (STROBE-MR): explanation and elaboration". BMJ. 375: n2233. doi:10.1136/bmj.n2233. ISSN 1756-1833. PMC 8546498. PMID 34702754.