TY - JOUR
T1 - When drug treatments bias genetic studies
T2 - Mediation and interaction
AU - Schmidt, Amand F.
AU - Heerspink, Hiddo J. L.
AU - Denig, Petra
AU - Finan, Chris
AU - Groenwold, Rolf H. H.
PY - 2019/8/1
Y1 - 2019/8/1
N2 - Background Increasingly, genetic analyses are conducted using information from subjects with established disease, who often receive concomitant treatment. We determined when treatment may bias genetic associations with a quantitative trait. Methods Graph theory and simulated data were used to explore the impact of drug prescriptions on (longitudinal) genetic effect estimates. Analytic derivations of longitudinal genetic effects are presented, accounting for the following scenarios: 1) treatment allocated independently of a genetic variant, 2) treatment that mediates the genetic effect, 3) treatment that modifies the genetic effect. We additionally evaluate treatment modelling strategies on bias, the root mean squared error (RMSE), coverage, and rejection rate. Results We show that in the absence of treatment by gene effect modification or mediation, genetic effect estimates will be unbiased. In simulated data we found that conditional models accounting for treatment, confounding, and effect modification were generally unbiased with appropriate levels of confidence interval coverage. Ignoring the longitudinal nature of treatment prescription, however (e.g. because of incomplete records in longitudinal data), biased these conditional models to a similar degree (or worse) as simply ignoring treatment. Conclusion The mere presence of (drug) treatment affecting a GWAS phenotype is insufficient to bias genetic associations with quantitative traits. While treatment may bias associations through effect modification and mediation, this might not occur frequently enough to warrant general concern at the presence of treated subjects in GWAS. Should treatment by gene effect modification or mediation be present however, current GWAS approaches attempting to adjust for treatment insufficiently account for the multivariable and longitudinal nature of treatment trajectories and hence genetic estimates may still be biased.
AB - Background Increasingly, genetic analyses are conducted using information from subjects with established disease, who often receive concomitant treatment. We determined when treatment may bias genetic associations with a quantitative trait. Methods Graph theory and simulated data were used to explore the impact of drug prescriptions on (longitudinal) genetic effect estimates. Analytic derivations of longitudinal genetic effects are presented, accounting for the following scenarios: 1) treatment allocated independently of a genetic variant, 2) treatment that mediates the genetic effect, 3) treatment that modifies the genetic effect. We additionally evaluate treatment modelling strategies on bias, the root mean squared error (RMSE), coverage, and rejection rate. Results We show that in the absence of treatment by gene effect modification or mediation, genetic effect estimates will be unbiased. In simulated data we found that conditional models accounting for treatment, confounding, and effect modification were generally unbiased with appropriate levels of confidence interval coverage. Ignoring the longitudinal nature of treatment prescription, however (e.g. because of incomplete records in longitudinal data), biased these conditional models to a similar degree (or worse) as simply ignoring treatment. Conclusion The mere presence of (drug) treatment affecting a GWAS phenotype is insufficient to bias genetic associations with quantitative traits. While treatment may bias associations through effect modification and mediation, this might not occur frequently enough to warrant general concern at the presence of treated subjects in GWAS. Should treatment by gene effect modification or mediation be present however, current GWAS approaches attempting to adjust for treatment insufficiently account for the multivariable and longitudinal nature of treatment trajectories and hence genetic estimates may still be biased.
UR - https://www.scopus.com/inward/record.uri?partnerID=HzOxMe3b&scp=85071425770&origin=inward
UR - https://www.ncbi.nlm.nih.gov/pubmed/31461463
U2 - https://doi.org/10.1371/journal.pone.0221209
DO - https://doi.org/10.1371/journal.pone.0221209
M3 - Article
C2 - 31461463
SN - 1932-6203
VL - 14
JO - PLOS ONE
JF - PLOS ONE
IS - 8
M1 - e0221209
ER -