Mendelian randomization analysis of a time-varying exposure for binary disease outcomes using functional data analysis methods

Approaches to estimate bidirectional causal effects using Mendelian randomization with application to body mass index and fasting glucose

Mendelian randomization (MR) is an epidemiological framework using genetic variants as instrumental variables (IVs) to examine the causal effect of exposures on outcomes. Statistical methods based on unidirectional MR (UMR) are widely used to estimate the causal effects of exposures on outcomes in observational studies. To estimate the bidirectional causal effects between two phenotypes, investigators have naively applied UMR methods separately in each direction. However, bidirectional causal effects between two phenotypes create a feedback loop that biases the estimation when UMR methods are naively applied. To overcome this limitation, we proposed two novel approaches to estimate bidirectional causal effects using MR: BiRatio and BiLIML, which are extensions of the standard ratio, and limited information maximum likelihood (LIML) methods, respectively. We compared the performance of the two proposed methods with the naive application of UMR methods through extensive simulations of several scenarios involving varying numbers of strong and weak IVs. Our simulation results showed that when multiple strong IVs are used, the proposed methods provided accurate bidirectional causal effect estimation in terms of median absolute bias and relative median absolute bias. Furthermore, compared to the BiRatio method, the BiLIML method provided a more accurate estimation of causal effects when weak IVs were used. Therefore, based on our simulations, we concluded that the BiLIML should be used for bidirectional causal effect estimation. We applied the proposed methods to investigate the potential bidirectional relationship between obesity and diabetes using the data from the Multi-Ethnic Study of Atherosclerosis cohort. We used body mass index (BMI) and fasting glucose (FG) as measures of obesity and type 2 diabetes, respectively. Our results from the BiLIML method revealed the bidirectional causal relationship between BMI and FG in across all racial populations. Specifically, in the White/Caucasian population, a 1 kg/m2 increase in BMI increased FG by 0.70 mg/dL (95% confidence interval [CI]: 0.3517-1.0489; p = 8.43×10-5), and 1 mg/dL increase in FG increased BMI by 0.10 kg/m2 (95% CI: 0.0441-0.1640; p = 6.79×10-4). Our study provides novel findings and quantifies the effect sizes of the bidirectional causal relationship between BMI and FG. However, further studies are needed to understand the biological and functional mechanisms underlying the bidirectional pathway.

Methodological approaches, challenges, and opportunities in the application of Mendelian randomisation to lifecourse epidemiology: A systematic literature review

Diseases diagnosed in adulthood may have antecedents throughout (including prenatal) life. Gaining a better understanding of how exposures at different stages in the lifecourse influence health outcomes is key to elucidating the potential benefits of disease prevention strategies. Mendelian randomisation (MR) is increasingly used to estimate causal effects of exposures across the lifecourse on later life outcomes. This systematic literature review explores MR methods used to perform lifecourse investigations and reviews previous work that has utilised MR to elucidate the effects of factors acting at different stages of the lifecourse. We conducted searches in PubMed, Embase, Medline and MedRXiv databases. Thirteen methodological studies were identified. Four studies focused on the impact of time-varying exposures in the interpretation of "standard" MR techniques, five presented methods for repeat measures of the same exposure, and four described methodological approaches to handling multigenerational exposures. A further 127 studies presented the results of an applied research question. Over half of these estimated effects in a single generation and were largely confined to the exploration of questions regarding body composition. The remaining mostly estimated maternal effects. There is a growing body of research focused on the development and application of MR methods to address lifecourse research questions. The underlying assumptions require careful consideration and the interpretation of results rely on select conditions. Whilst we do not advocate for a particular strategy, we encourage practitioners to make informed decisions on how to approach a research question in this field with a solid understanding of the limitations present and how these may be affected by the research question, modelling approach, instrument selection, and data availability.

Landmark mediation survival analysis using longitudinal surrogate

Clinical cancer trials are designed to collect radiographic measurements of each patient's baseline and residual tumor burden at regular intervals over the course of study. For solid tumors, the extent of reduction in tumor size following treatment is used as a measure of a drug's antitumor activity. Statistical estimation of treatment efficacy routinely reduce the longitudinal assessment of tumor burden to a binary outcome describing the presence versus absence of an objective tumor response as defined by RECIST criteria. The objective response rate (ORR) is the predominate method for evaluating an experimental therapy in a single-arm trial. Additionally, ORR is routinely compared against a control therapy in phase III randomized controlled trials. The longitudinal assessments of tumor burden are seldom integrated into a formal statistical model, nor integrated into mediation analysis to characterize the relationships among treatment, residual tumor burden, and survival. This article presents a frameworkfor landmark mediation survival analyses devised to incorporate longitudinal assessment of tumor burden. R ² effect-size measures are developed to quantify the survival treatment mediation effects using longitudinal predictors. Analyses are demonstrated with applications to two colorectal cancer trials. Survival prediction is compared in the presence versus absence of longitudinal analysis. Simulation studies elucidate settings wherein patterns of tumor burden dynamics require longitudinal analysis.

Interpretation of Mendelian randomization using a single measure of an exposure that varies over time

BACKGROUND

Mendelian randomization (MR) is a powerful tool through which the causal effects of modifiable exposures on outcomes can be estimated from observational data. Most exposures vary throughout the life course, but MR is commonly applied to one measurement of an exposure (e.g. weight measured once between ages 40 and 60 years). It has been argued that MR provides biased causal effect estimates when applied to one measure of an exposure that varies over time.

METHODS

We propose an approach that emphasizes the liability that causes the entire exposure trajectory. We demonstrate this approach using simulations and an applied example.

RESULTS

We show that rather than estimating the direct or total causal effect of changing the exposure value at a given time, MR estimates the causal effect of changing the underlying liability for the exposure, scaled to the effect of the liability on the exposure at that time. As such, results from MR conducted at different time points are expected to differ (unless the effect of the liability on exposure is constant over time), as we illustrate by estimating the effect of body mass index measured at different ages on systolic blood pressure.

CONCLUSION

Univariable MR results should not be interpreted as time-point-specific direct or total causal effects, but as the effect of changing the liability for the exposure. Estimates of how the effects of a genetic variant on an exposure vary over time, together with biological knowledge that provides evidence regarding likely effective exposure periods, are required to interpret time-point-specific causal effects.

Strengthening Causal Inference in Exposomics Research: Application of Genetic Data and Methods

Advances in technologies to measure a broad set of exposures have led to a range of exposome research efforts. Yet, these efforts have insufficiently integrated methods that incorporate genetic data to strengthen causal inference, despite evidence that many exposome-associated phenotypes are heritable. Objective: We demonstrate how integration of methods and study designs that incorporate genetic data can strengthen causal inference in exposomics research by helping address six challenges: reverse causation and unmeasured confounding, comprehensive examination of phenotypic effects, low efficiency, replication, multilevel data integration, and characterization of tissue-specific effects. Examples are drawn from studies of biomarkers and health behaviors, exposure domains where the causal inference methods we describe are most often applied. Discussion: Technological, computational, and statistical advances in genotyping, imputation, and analysis, combined with broad data sharing and cross-study collaborations, offer multiple opportunities to strengthen causal inference in exposomics research. Full application of these opportunities will require an expanded understanding of genetic variants that predict exposome phenotypes as well as an appreciation that the utility of genetic variants for causal inference will vary by exposure and may depend on large sample sizes. However, several of these challenges can be addressed through international scientific collaborations that prioritize data sharing. Ultimately, we anticipate that efforts to better integrate methods that incorporate genetic data will extend the reach of exposomics research by helping address the challenges of comprehensively measuring the exposome and its health effects across studies, the life course, and in varied contexts and diverse populations. https://doi.org/10.1289/EHP9098.

Functional principal component based landmark analysis for the effects of longitudinal cholesterol profiles on the risk of coronary heart disease

Patients' longitudinal biomarker changing patterns are crucial factors for their disease progression. In this research, we apply functional principal component analysis techniques to extract these changing patterns and use them as predictors in landmark models for dynamic prediction. The time-varying effects of risk factors along a sequence of landmark times are smoothed by a supermodel to borrow information from neighbor time intervals. This results in more stable estimation and more clear demonstration of the time-varying effects. Compared with the traditional landmark analysis, simulation studies show our proposed approach results in lower prediction error rates and higher area under receiver operating characteristic curve (AUC) values, which indicate better ability to discriminate between subjects with different risk levels. We apply our method to data from the Framingham Heart Study, using longitudinal total cholesterol (TC) levels to predict future coronary heart disease (CHD) risk profiles. Our approach not only obtains the overall trend of biomarker-related risk profiles, but also reveals different risk patterns that are not available from the traditional landmark analyses. Our results show that high cholesterol levels during young ages are more harmful than those in old ages. This demonstrates the importance of analyzing the age-dependent effects of TC on CHD risk.

Smoking and multiple sclerosis risk: a Mendelian randomization study

BACKGROUND

Striking changes in the demographic pattern of multiple sclerosis (MS) strongly indicate an influence of modifiable exposures, which lend themselves well to intervention. It is important to pinpoint which of the many environmental, lifestyle, and sociodemographic changes that have occurred over the past decades, such as higher smoking and obesity rates, are responsible. Mendelian randomization (MR) is an elegant tool to overcome limitations inherent to observational studies and leverage human genetics to inform prevention strategies in MS.

METHODS

We use genetic variants from the largest genome-wide association study for smoking phenotypes (initiation: N = 378, heaviness: N = 55, lifetime smoking: N = 126) and body mass index (BMI, N = 656) and apply these as instrumental variables in a two-sample MR analysis to the most recent meta-analysis for MS. We adjust for the genetic correlation between smoking and BMI in a multivariable MR.

RESULTS

In univariable and multivariable MR, smoking does not have an effect on MS risk nor explains part of the association between BMI and MS risk. In contrast, in both analyses each standard deviation increase in BMI, corresponding to roughly 5 kg/m² units, confers a 30% increase in MS risk.

CONCLUSION

Despite observational studies repeatedly reporting an association between smoking and increased risk for MS, MR analyses on smoking phenotypes and MS risk could not confirm a causal relationship. This is in contrast with BMI, where observational studies and MR agree on a causal contribution. The reasons for the discrepancy between observational studies and our MR study concerning smoking and MS require further investigation.

Birthweight, time-varying adiposity growth and early menarche in girls: A Mendelian randomisation and mediation analysis

OBJECTIVE

To explore the causal effect of time-varying z-BMI growth on early menarche using Mendelian randomisation (MR); to identify critical adiposity predictors of early menarche; to compare the effects of birthweight and time-varying z-BMI growth as mediators of the path from genes to early menarche using mediation analysis.

METHODS

We used data from the Taiwan Children Health Study with 21 obesity-related single-nucleotide polymorphisms (SNPs) to yield genetic (instrumental variable)IVs for adiposity. Children with available data on genotyping, birthweight, adiposity, and menarcheal age were included.

RESULTS

In MR analyses, results based on the time-varying z-BMI growth show more statistical power and capture more information of adiposity growth (p=0.01) than those based on single point z-BMI (p=0.02). Among adiposity measures, critical predictors of early menarche are fat free mass (RR=1.33, 95% CI 1.07-1.65) and waist/height ratio (RR=1.27, 95% CI 1.03-1.56). Other potential predictors of early menarche are sum of skinfold (RR=1.24, 95% CI 1.03-1.48) and total body fat (RR=1.20, 95% CI 1.05-1.38). In both one-mediation and multi-mediation analyses, time-varying z-BMI growth in the prepubertal years plays a crucial mediator in the pathway from the genes to early menarche.

CONCLUSIONS

This study discovered that greater prepubertal adiposity growth is a crucial mediator in the path from genes to early menarche. For girls with genes positively associated with obesity; and/or of lower birthweight, a strategy to prevent childhood adiposity should be implemented in order to avoid early menarche development.

Software Application Profile: PHESANT: a tool for performing automated phenome scans in UK Biobank

Motivation

Epidemiological cohorts typically contain a diverse set of phenotypes such that automation of phenome scans is non-trivial, because they require highly heterogeneous models. For this reason, phenome scans have to date tended to use a smaller homogeneous set of phenotypes that can be analysed in a consistent fashion. We present PHESANT (PHEnome Scan ANalysis Tool), a software package for performing comprehensive phenome scans in UK Biobank.

General features

PHESANT tests the association of a specified trait with all continuous, integer and categorical variables in UK Biobank, or a specified subset. PHESANT uses a novel rule-based algorithm to determine how to appropriately test each trait, then performs the analyses and produces plots and summary tables.

Implementation

The PHESANT phenome scan is implemented in R. PHESANT includes a novel Javascript D3.js visualization and accompanying Java code that converts the phenome scan results to the required JavaScript Object Notation (JSON) format.

Availability

PHESANT is available on GitHub at [https://github.com/MRCIEU/PHESANT]. Git tag v0.5 corresponds to the version presented here.