Real-world evidence (RWE) is increasingly considered in regulatory decision making. When, and to which extent, RWE is considered relevant by regulators likely depends on many factors. This review... Show moreReal-world evidence (RWE) is increasingly considered in regulatory decision making. When, and to which extent, RWE is considered relevant by regulators likely depends on many factors. This review aimed to identify factors that make RWE necessary or desirable to inform regulatory decision making. A scoping review was conducted using literature databases (PubMed, Embase, Emcare, Web of Science, and Cochrane Library) and websites of regulatory agencies, health technology assessment agencies, research institutes, and professional organizations involved with RWE. Articles were included if: (1) they discussed factors or contexts that impact whether RWE could be necessary or desirable in regulatory decision making; (2) focused on pharmacological or biological interventions in humans; and (3) considered decision making in Europe or North America, or without a focus on a specific region. We included 118 articles in the scoping review. Two major themes and six subthemes were identified. The first theme concerns questions addressable with RWE, with subthemes epidemiology and benefit-risk assessment. The second theme concerns contextual factors, with subthemes feasibility, ethical considerations, limitations of available evidence, and disease and treatment-specific aspects. Collectively, these themes encompassed 43 factors influencing the need for RWE in regulatory decisions. Although single factors may not make RWE fully necessary, their cumulative influence could make RWE essential and pivotal in regulatory decision making. This overview contributes to ongoing discussions emphasizing the nuanced interplay of factors influencing the necessity or desirability of RWE to inform regulatory decision making. Show less
Nemeth, B.; Smeets, M.; Pedersen, A.B.; Kristiansen, E.B.; Nelissen, R.; Whyte, M.; ... ; Arya, R. 2024
Background: The risk of venous thromboembolism (VTE) following total hip arthroplasty (THA) and total knee arthroplasty (TKA) is 1.0% to 1.5%, despite uniform thromboprophylaxis. Objectives: To... Show moreBackground: The risk of venous thromboembolism (VTE) following total hip arthroplasty (THA) and total knee arthroplasty (TKA) is 1.0% to 1.5%, despite uniform thromboprophylaxis. Objectives: To develop and validate a prediction model for 90-day VTE risk. Methods: A multinational cohort study was performed. For model development, records were used from the Oxford Royal College of General Practitioners Research and Surveillance Centre linked to Hospital Episode Statistics and Office of National Statistics UK routine data. For external validation, data were used from the Danish Hip and Knee Arthroplasty Registry, the National Patient Registry, and the National Prescription Registry. Binary multivariable logistic regression techniques were used for development. Results: In the UK data set, 64 032 THA/TKA procedures were performed and 1.4% developed VTE. The prediction model consisted of age, body mass index, sex, cystitis within 1 year before surgery, history of phlebitis, history of VTE, presence of varicose veins, presence of asthma, history of transient ischemic attack, history of myocardial infarction, presence of hypertension and THA or TKA. The area under the curve of the model was 0.65 (95% CI, 0.63-0.67). Furthermore, 36 169 procedures were performed in the Danish cohort, of whom 1.0% developed VTE. Here, the area under the curve was 0.64 (95% CI, 0.61-0.67). The calibration slope was 0.92 in the validation study and 1.00 in the development study. Conclusion: This clinical prediction model for 90-day VTE risk following THA and TKA performed well in both development and validation data. This model can be used to estimate an individual’s risk for VTE following THA/TKA. Show less
Tweel, M.M. van den; Munckhof, E.H.A. van den; Zanden, M. van der; Molijn, A.C.; Lith, J.M.M. van; Cessie, S. le; Boers, K.E. 2023
PurposeThis study investigates the role of bacterial vaginosis (BV) on pregnancy rates during various fertility treatments. BV is known to influence several obstetric outcomes, such as preterm... Show morePurposeThis study investigates the role of bacterial vaginosis (BV) on pregnancy rates during various fertility treatments. BV is known to influence several obstetric outcomes, such as preterm delivery and endometritis. Only few studies investigated the effect of BV in subfertile women, and studies found a negative effect on fecundity especially in the in vitro fertilisation population.MethodsObservational prospective study, 76 couples attending a fertility clinic in the Netherlands between July 2019 and June 2022, undergoing a total of 133 attempts of intra uterine insemination, in vitro fertilization or intra cytoplasmatic sperm injection. Vaginal samples taken at oocyte retrieval or insemination were analysed on qPCR BV and 16S rRNA gene microbiota analysis of V1-V2 region. Logistic regression with a Generalized Estimated Equations analysis was used to account for multiple observations per couples.ResultsA total of 26% of the 133 samples tested positive for BV. No significant differences were observed in ongoing pregnancy or live birth rates based on BV status (OR 0.50 (0.16-1.59), aOR 0.32 (0.09-1.23)) or microbiome community state type. There was a tendency of more miscarriages based on positive BV status (OR 4.22 (1.10-16.21), aOR 4.28 (0.65-28.11)) or community state type group III and IV. On baseline qPCR positive participants had significantly higher body mass index and smoked more often. Odds ratios were adjusted for smoking status, body mass index, and socioeconomic status.ConclusionBacterial vaginosis does not significantly impact ongoing pregnancy rates but could affect miscarriage rates. Show less
Heckert, S.L.; Maassen, J.M.; Cessie, S. le; Goekoop-Ruiterman, Y.P.M.; Güler-Yüksel, M.; Lems, W.; ... ; Allaart, C.F. 2023
Objectives To study long-term (up to 20-year) mortality of two treat-to-target trial cohorts in undifferentiated arthritis (UA) and early rheumatoid arthritis (RA).Methods The BeSt ... Show moreObjectives To study long-term (up to 20-year) mortality of two treat-to-target trial cohorts in undifferentiated arthritis (UA) and early rheumatoid arthritis (RA).Methods The BeSt (BehandelStrategieen) study (n=508, early RA) was performed between 2000 and 2012. For 10 years, patients were treated-to-target disease activity score (DAS)<= 2.4.The Induction therapy with Methotrexate and Prednisone in Rheumatoid Or Very Early arthritic Disease (IMPROVED) study (n=610, early RA/UA) was performed between 2007 and 2015. For 5 years, patients were treated-to-target DAS<1.6.Vital status of BeSt/IMPROVED participants was assessed up to and including 31 December 2021. Standardised mortality ratios (SMRs) were calculated. Stratified analyses for anticitrullinated protein antibody (ACPA) and smoking status were performed. Death causes and the potential effect of disease activity during the trial period on late mortality were assessed.Results Excess mortality was found in both BeSt (SMR 1.32, 95% CI 1.14 to 1.53) and IMPROVED (SMR 1.33, 95% CI 1.10 to 1.63) and became manifest after 10 years. Excess mortality was statistically significant in ACPA+ patients who smoked (BeSt: SMR 2.80, 95% CI 2.16 to 3.64; IMPROVED: 2.14, 95% CI 1.33 to 3.45). Mean survival time was 10 (95% CI 5 to 16) months shorter than expected in BeSt and 13 (95% CI 11 to 16) months in IMPROVED. The HR for mortality was 1.34 (95% CI 0.96 to 1.86; BeSt)/1.13 (95% CI 0.67 to 1.91; IMPROVED) per 1 point increase in mean DAS during the trial. The main cause of death was malignancy.Conclusions After long-term treatment-to-target, excess mortality occurred in patients with RA after>10 years since treatment start, with smoking as an important risk factor. Show less
Camilleri, E.; Smit, D.L.; Rein, N. van; Cessie, S. le; Hon, O. de; Heijer, M. den; ... ; Ronde, W. de 2023
Background: Anabolic androgenic steroids (AAS) are thought to increase venous thromboembolism (VTE) risk. Objectives: We investigated whether AAS influence coagulation parameters associated with... Show moreBackground: Anabolic androgenic steroids (AAS) are thought to increase venous thromboembolism (VTE) risk. Objectives: We investigated whether AAS influence coagulation parameters associated with VTE by assessing their changes during and after AAS use. Methods: The HAARLEM study enrolled 100 male amateur athletes voluntarily starting an AAS cycle between 2015 and 2018. We measured procoagulant and anticoagulant protein levels, D-dimer levels, endogenous thrombin potential (ETP), and clot lysis time (CLT) at baseline and during 2 years of follow-up. Changes in coagulation during AAS cycle, 3 months after its discontinuation, and 1 year after its inclusion compared with baseline were estimated using linear mixed models. The associations between AAS dose and duration of use with these outcomes were studied through adjusted multivariable linear regression. Results: Participants used AAS for a median of 13 weeks (IQR: 10-23) with a median weekly dose of 901 mg (IQR: 634-1345 mg). Mean levels of multiple coagulation factors (F) increased during use compared with baseline, whereas FVIII and von Willebrand factor levels remained unchanged. Protein S and D-dimer showed the biggest increase (22% [95% CI: 15-29] and 1.3-fold [95% CI: 1.2-1.5], respectively). CLT was 8 minutes longer (95% CI: 5-10) and ETP was 165 nM*min (95% CI: -205 to -124) lower during the AAS cycle. A high weekly AAS dose and short cycle duration were associated with changes in protein S and ETP during use. All parameters returned to baseline values 3 months after discontinuation and remained similar after. Conclusion: During AAS use, procoagulant and anticoagulant protein levels increased in a reversible manner. The overall balance did not suggest a clear procoagulant state. Show less
Rationale: Supplemental oxygen is widely administered to ICU patients, but appropriate oxygenation targets remain unclear.Objectives: This study aimed to determine whether a low-oxygenation... Show moreRationale: Supplemental oxygen is widely administered to ICU patients, but appropriate oxygenation targets remain unclear.Objectives: This study aimed to determine whether a low-oxygenation strategy would lower 28-day mortality compared with a high-oxygenation strategy.Methods: This randomized multicenter trial included mechanically ventilated ICU patients with an expected ventilation duration of at least 24 hours. Patients were randomized 1:1 to a low-oxygenation (Pa-O2, 55-80mmHg; or oxygen saturation as measured by pulse oximetry, 91-94%) or high-oxygenation (Pa-O2, 110-150mmHg; or oxygen saturation as measured by pulse oximetry, 96-100%) target until ICU discharge or 28 days after randomization, whichever came first. The primary outcome was 28-day mortality. The study was stopped prematurely because of the COVID-19 pandemic when 664 of the planned 1,512 patients were included.Measurements and Main Results: Between November 2018 and November 2021, a total of 664 patients were included in the trial: 335 in the low-oxygenation group and 329 in the high-oxygenation group. The median achieved Pa-O2 was 75mmHg (interquartile range, 70-84) and 115mmHg (interquartile range, 100-129) in the low- and high-oxygenation groups, respectively. At Day 28, 129 (38.5%) and 114 (34.7%) patients had died in the low- and high-oxygenation groups, respectively (risk ratio, 1.11; 95% confidence interval, 0.9-1.4; P = 0.30). At least one serious adverse event was reported in 12 (3.6%) and 17 (5.2%) patients in the low- and high-oxygenation groups, respectively.Conclusions: Among mechanically ventilated ICU patients with an expected mechanical ventilation duration of at least 24 hours, using a low-oxygenation strategy did not result in a reduction of 28-day mortality compared with a high-oxygenation strategy. Show less
Van't Oever, R.M.; Zwiers, C.; Haas, M. de; Cessie, S. le; Lopriore, E.; Oepkes, D.; Verweij, E.J.T. 2023
ObjectivePregnant women who received at least one intrauterine transfusion (IUT) for haemolytic disease of the fetus and newborn (HDFN) in the preceding pregnancy are presumed to have a high... Show moreObjectivePregnant women who received at least one intrauterine transfusion (IUT) for haemolytic disease of the fetus and newborn (HDFN) in the preceding pregnancy are presumed to have a high likelihood of requiring IUTs again, often starting at an earlier gestational age. Our aim was to quantify these risks in a large national cohort.DesignRetrospective cohort study of a nationwide Dutch database.SettingThe Netherlands.PopulationAll women treated in The Netherlands with IUTs for Rhesus D (RhD)- or Kell-mediated HDFN between 1999 and 2017 and their follow-up pregnancies were included. Pregnancies with an antigen-negative fetus were excluded.MethodsElectronic patient files were searched for the number and gestational age of each IUT, and analysed using descriptive statistics and linear regression.Main outcome measuresPercentage of women requiring one or more IUTs again in the subsequent pregnancy, and gestational age at first IUT in both pregnancies.ResultsOf the 321 women in our study population, 21% (69) had a subsequent ongoing pregnancy at risk. IUTs were administered in 86% (59/69) of cases. In subsequent pregnancies, the median gestational age at first IUT was 3 weeks earlier (interquartile range -6.8 to 0.4) than in the preceding pregnancy.ConclusionsOur study shows that pregnant women with a history of IUTs in the previous pregnancy are highly likely to require IUTs again, and on average 3 weeks earlier. Clinicians need to be aware of these risks and ensure timely referral, and close surveillance from early pregnancy onwards. Additionally, for women with a history of IUT and their caregivers, this information is essential to enable adequate preconception counselling. Show less
Martens, L.G.; Hamersveld, D. van; Cessie, S. le; Dijk, K.W. van; Heemst, D. van; Noordam, R. 2023
Objectives: Low socioeconomic status (SES) is associated with cardiovascular risk factors and increased coronary artery disease (CAD) risk. We tested whether SES is an effect modifier of the... Show moreObjectives: Low socioeconomic status (SES) is associated with cardiovascular risk factors and increased coronary artery disease (CAD) risk. We tested whether SES is an effect modifier of the association between classical cardiovascular risk factors and CAD using SES-stratified Mendelian Randomization in European-ancestry participants from UK Biobank.Study Design and Setting: We calculated weighted genetic risk scores (GRS) for the risk factors body mass index (BMI), systolic blood pressure, low-density lipoprotein cholesterol, and triglycerides. Participants were stratified by Townsend deprivation index score. Lo-gistic regression models were used to investigate associations between GRSs and CAD occurrence. Additionally, stratification based on GRS-adjusted Townsend deprivation index residuals was conducted to correct for possible collider-stratification bias.Results: In a total sample size of N 5 446,485, with 52,946 cases, the risk for CAD per standard deviation increase in genetically influenced BMI was highest in the group with the lowest 25% SES (odds ratio: 1.126, 95% confidence interval: 1.106-1.145; odds ratio: 1.081, 95% confidence interval: 1.059-1.103 in high SES), remaining similar after controlling for possible collider-stratification bias. The effects of genetically influenced systolic blood pressure, low-density lipoprotein cholesterol, and triglyceride on CAD were similar between SES groups.Conclusion: CAD risk attributable to increased BMI is not homogenous and could be modified by SES. This emphasizes the need of tailor-made approaches for BMI-associated CAD risk reduction. & COPY; 2023 The Author(s). Published by Elsevier Inc. This is an open access article under the CC BY license (http://creativecommons.org/licenses/by/4.0/). Show less
Rationale & Objective: Hypokalemia may accelerate kidney function decline. Both hypo- and hyperkalemia can cause sudden cardiac death. However, little is known about the relationship between... Show moreRationale & Objective: Hypokalemia may accelerate kidney function decline. Both hypo- and hyperkalemia can cause sudden cardiac death. However, little is known about the relationship between serum potassium and death or the occurrence of kidney failure requiring replacement therapy (KRT). We investigated this relationship in older people with chronic kidney disease (CKD) stage 4-5.Study Design: Prospective observational cohort study.Setting & Participants: We followed 1,714 patients (>= 65 years old) from the European Quality (EQUAL) study for 8 years from their first estimated glomerular filtration rate (eGFR) < 20 mL/min/1.73 m(2) measurement. Exposure: Serum potassium was measured every 3 to 6 months and categorized as <= 3.5, >3.5-<= 4.0, >4.0-<= 4.5, >4.5-<= 5.0 (reference), >5.0-<= 5.5, >5.5-<= 6.0, and >6.0 mmol/L.Outcome: The combined outcome death before KRT or start of KRT.Analytical Approach: The association between categorical and continuous time-varying potassium and death or KRT start was examined using Cox proportional hazards and restricted cubic spline analyses, adjusted for age, sex, diabetes, cardiovascular disease, renin-angiotensin-aldosterone system (RAAS) inhibition, eGFR, and subjective global assessment (SGA).Results: At baseline, 66% of participants were men, 42% had diabetes, 47% cardiovascular disease, and 54% used RAAS inhibitors. Their mean age was 76 +/- 7 (SD) years, mean eGFR was 17 +/- 5 (SD) mL/min/1.73 m(2), and mean SGA was 6.0 +/- 1.0 (SD). Over 8 years, 414 (24%) died before starting KRT, and 595 (35%) started KRT. Adjusted hazard ratios for death or KRT according to the potassium categories were 1.6 (95% CI, 1.1-2.3), 1.4 (95% CI, 1.1-1.7), 1.1 (95% CI, 1.0-1.4), 1 (reference), 1.1 (95% CI, 0.9-1.4), 1.8 (95% CI, 1.4-2.3), and 2.2 (95% CI, 1.5-3.3). Hazard ratios were lowest at a potassium of about 4.9 mmol/L.Limitations: Shorter intervals between potassium measurements would have allowed for more precise estimations.Conclusions: We observed a U-shaped relationship between serum potassium and death or KRT start among patients with incident CKD 4- 5, with a nadir risk at a potassium level of 4.9 mmol/L. These findings underscore the potential importance of preventing both high and low potassium in patients with CKD 4-5. Show less
Background We explored the potential of emerging and conventional urinary kidney injury biomarkers in recipients of living donor (LD) or donation after circulatory death (DCD) kidney... Show moreBackground We explored the potential of emerging and conventional urinary kidney injury biomarkers in recipients of living donor (LD) or donation after circulatory death (DCD) kidney transplantation, patients with chronic kidney disease (CKD), and individuals from the general population. Methods Urine samples from kidney allograft recipients with mild (LD; n = 199) or severe (DCD; n = 71) ischemia-reperfusion injury (IRI) were analyzed for neutrophil gelatinase-associated lipocalin (NGAL), insulin-like growth factor-binding protein 7 (IGFBP7), tissue inhibitor of metalloproteinases 2 (TIMP2), kidney injury molecule-1 (KIM-1), chemokine C-X-C motif (CXCL9), solute carrier family 22 member 2 (SLC22A2), nephrin, and uromodulin (UMOD) by quantitative multiplex LC-MS/MS analysis. The fold-change in biomarker levels was determined in mild and severe IRI and in patients with CKD stage 1-2 (n = 127) or stage & GE;3 (n = 132) in comparison to the general population (n = 1438). Relationships between the biomarkers and total protein, & beta;2-microglobulin (B2M), creatinine, and osmolality were assessed. Results NGAL, IGFBP7, TIMP2, KIM-1, CXCL9, and UMOD were quantifiable, whereas nephrin and SLC22A2 were below the limit of detection. Kidney injury biomarkers were increased up to 6.2-fold in allograft recipients with mild IRI and 8.3-fold in recipients with severe IRI, compared to the reference population, with the strongest response observed for NGAL and B2M. In CKD stage 1-2, B2M, NGAL, IGFBP7, TIMP2, KIM-1, UMOD, and CXCL9 were not altered, but in individuals with CKD stage & GE;3, B2M, NGAL, and KIM-1 were increased up to 1.3-fold. IGFBP7, TIMP2, NGAL, and CXCL9 were strongly correlated (all r & GE; 0.8); correlations with B2M and TP were smaller (all r & LE; 0.6). Conclusions IRI, but not stable CKD, was associated with increased urinary levels of kidney injury biomarkers determined by LC-MS/MS. Absolute and multiplexed protein quantitation by LC-MS/MS is an effective strategy for biomarker panel evaluation for translation toward the clinical laboratory. Show less
Verhagen, I.E.; Lentsch, S.D.; Arend, B.W.H. van der; Cessie, S. le; MaassenVanDenBrink, A.; Terwindt, G.M. 2023
Background and purpose Anti-calcitonin gene-related peptide (CGRP) (receptor) antibodies effectively reduce overall migraine attack frequency, but whether there are differences in effect between... Show moreBackground and purpose Anti-calcitonin gene-related peptide (CGRP) (receptor) antibodies effectively reduce overall migraine attack frequency, but whether there are differences in effect between perimenstrual and nonperimenstrual migraine days has not been investigated.Methods We performed a single-arm study among women with migraine. Participants were followed with electronic E-diaries during one (pretreatment) baseline month and 6 months of treatment with either erenumab or fremanezumab. Differences in treatment effect on perimenstrual and nonperimenstrual migraine days were assessed using a mixed effects logistic regression model, with migraine day as dependent variable; treatment, menstrual window, and an interaction term (treatment x menstrual window) as fixed effects; and patient as a random effect.Results There was no interaction between the menstrual window and treatment effect, indicating that the reduction in migraine days under treatment was similar during the menstrual window and the remainder of the menstrual cycle (odds ratio for treatment = 0.44, 95% confidence interval = 0.38-0.51).Conclusions Our findings support prophylactic use of anti-CGRP (receptor) antibodies for women with menstrual migraine, as this leads to consistent reductions in number of migraine days during the entire menstrual cycle. Show less
In clinical settings, the absolute risk reduction due to treatment that can be expected in a particular patient is of key interest. However, logistic regression, the default regression model for... Show moreIn clinical settings, the absolute risk reduction due to treatment that can be expected in a particular patient is of key interest. However, logistic regression, the default regression model for trials with a binary outcome, produces estimates of the effect of treatment measured as a difference in log odds. We explored options to estimate treatment effects directly as a difference in risk, specifically in the network meta-analysis setting. We propose a novel Bayesian (meta-)regression model for binary outcomes on the additive risk scale. The model allows treatment effects, covariate effects, interactions and variance parameters to be estimated directly on the linear scale of clinical interest. We compared effect estimates from this model to (1) a previously proposed additive risk model by Warn, Thompson and Spiegelhalter ( "WTS model ") and (2) backtransforming the predictions from a logistic model to the natural scale after regression. The models were compared in a network meta-analysis of 20 hepatitis C trials, as well as in the analysis of simulated single trial settings. The resulting estimates diverged, in particular for small sample sizes or true risks close to 0% or 100%. Researchers should be aware that modelling untransformed risk can yield very different results from default logistic models. The treatment effect in participants with such extreme predicted risks weighed more heavily on the overall treatment effect estimate from our proposed model compared to the WTS model. In our network meta-analysis, this sensitivity of our proposed model was needed to detect all information in the data. Show less
Verhagen, I.E.; Arend, B.W.H. van der; Casteren, D.S. van; Cessie, S. le; MaassenVanDenBrink, A.; Terwindt, G.M. 2023
Objective: In this prospective cohort study, characteristics of perimenstrual and non-perimenstrual migraine attacks in women were compared with migraine attacks in men.Background: Women report... Show moreObjective: In this prospective cohort study, characteristics of perimenstrual and non-perimenstrual migraine attacks in women were compared with migraine attacks in men.Background: Women report longer migraine attacks and more accompanying symptoms than men in cross-sectional questionnaire studies, but this has not been confirmed in longitudinal studies. Supposed differences could result from different characteristics specific to perimenstrual migraine attacks, or of attacks in women in general.Methods: This cohort study was performed among patients with migraine who were treated at the Leiden Headache Clinic. We assessed differences in migraine attack characteristics between men and women who were prospectively followed by a previously validated electronic headache diary. The primary outcome was "attack" duration. Differences between perimenstrual (Days -2 to +3 of the menstrual cycle) and non-perimenstrual attacks in women versus attacks in men were corrected for age, chronic migraine, and medication overuse headache.Results: A total of 1347 women and 284 men were included, reflecting the preponderance of women in migraine prevalence. Crude median (first and third quartile [Q1-Q3]) attack duration in men was 32.1 [17.7-53.6] h, compared to 36.7 [21.9-62.4] h for non-perimenstrual migraine attacks and 44.4 [17.9-79.0] h for perimenstrual migraine attacks in women. After correction for confounding, perimenstrual migraine attacks were 1.62 (95% confidence interval [CI] 1.47-1.79; p < 0.001) and non-perimenstrual 1.15 (95% CI 1.05-1.25; p = 0.003) times longer compared to migraine attacks in men. The mean relapse percentage in men was 9.2%, compared to 12.6% for non-perimenstrual migraine attacks, and 15.7% for perimenstrual migraine attacks. Relapse risk was greater for perimenstrual (odds ratio [OR] 2.39, 95% CI 1.93-2.95; p < 0.001), but not for non-perimenstrual (OR 1.18, 95% CI 0.97-1.45; p = 0.060) attacks. Migraine attacks in women were more often accompanied by photophobia, phonophobia, and nausea, but less often aura.Conclusion: Compared to attacks in men, both perimenstrual and non-perimenstrual migraine attacks are of longer duration and are more often accompanied by associated symptoms. A sex-specific approach to migraine treatment and research is needed. Show less
Ill-defined research questions could be particularly problematic in an epidemiological setting where measurements fluctuate over time due to intercurrent events, such as medication use. When a... Show moreIll-defined research questions could be particularly problematic in an epidemiological setting where measurements fluctuate over time due to intercurrent events, such as medication use. When a research question fails to specify how medication use should be handled methodologically, arbitrary decisions may be made during the analysis phase, which likely leads to a mismatch between the intended question and the performed analysis. The mismatch can result in vastly different or meaningless interpretations of estimated effects. Thus, a research question such as "what is the effect of X on Y? " requires further elaboration, and it should consider whether and how medication use has affected the measurements of interest. In our study, we will discuss how well-defined questions can be formulated when medication use is involved in observational studies. We will distinguish between a situation where an exposure is affected by medication use and where the outcome of interest is affected by medication use. For each setting, we will give examples of different research questions that could be asked depending on how medication use is considered in the estimand and discuss methodological considerations under each question. Show less
Lengton, R.; Willik, E.M. van der; Rooij, E.N.M. de; Meuleman, Y.; Cessie, S. le; Michels, W.M.; ... ; NECOSAD Study Grp 2022
Background Chronic kidney disease-associated pruritus (CKD-aP) is common in dialysis patients, and is associated with lower quality of life and increased risk of death. We investigated the... Show moreBackground Chronic kidney disease-associated pruritus (CKD-aP) is common in dialysis patients, and is associated with lower quality of life and increased risk of death. We investigated the association between residual estimated glomerular filtration rate (eGFR), dialysis adequacy or serum phosphate level and CKD-aP in incident dialysis patients. Methods A total of 1256 incident hemodialysis (HD) and 670 peritoneal dialysis (PD) patients (>18 years) from the Netherlands Cooperative Study on the Adequacy of Dialysis (NECOSAD) study were included (1997-2007) and followed until death, transplantation or a maximum of 10 years. CKD-aP was measured using a single item of the Kidney Disease Quality of Life Instrument-36. The associations were studied by logistic and linear regression analyses, adjusted for potential baseline confounders. Results At baseline mean (standard deviation) age was 60 (16) years, 62% were men and median (interquartile range) residual eGFR was 3.4 (1.7; 5.3) mL/min/1.73 m(2). The prevalence of CKD-aP (similar to 70%) was similar in HD and PD. It was observed that 12 months after starting dialysis (after multivariable adjustment) each 1 mL/min/1.73 m(2) higher residual eGFR, one unit higher total weekly Kt/V, or 1 mmol/L lower serum phosphate level was associated with lower burden of CKD-aP in HD and PD patients of -0.05 (95% CI -0.09; -0.02) and -0.09 (95% CI -0.13; -0.05), -0.15 (95% CI -0.26; -0.05) and -0.35 (95% CI -0.54; -0.16), and of -0.34 (95%CI: -0.51; -0.17) and -0.45 (95%CI: -0.71; -0.19), respectively. We found no association between dialysis Kt/V and CKD-aP. Conclusions Higher residual eGFR and lower serum phosphate level, but not the dialysis dose, were related with lower burden of CKD-aP in dialysis patients. Show less
Luo, J.; Noordam, R.; Jukema, J.W.; Dijk, K.W. van; Hagg, S.; Grassmann, F.; ... ; Heemst, D. van 2022
Aim: Mitochondrial DNA dysfunction has been implicated in the pathogenesis of cardiovascular diseases. We aimed to investigate the associations between leukocyte mitochondrial DNA (mtDNA) abundance... Show moreAim: Mitochondrial DNA dysfunction has been implicated in the pathogenesis of cardiovascular diseases. We aimed to investigate the associations between leukocyte mitochondrial DNA (mtDNA) abundance, as a proxy of mitochondrial function, and coronary artery disease (CAD) and heart failure (HF) in a cohort study and approximate the causal nature of these relationships using Mendelian randomization (MR) in genetic studies. Methods and results: Multivariable-adjusted Cox regression analyses were conducted in 273 619 unrelated participants of European ancestry from the UK Biobank (UKB). For genetic studies, we first performed MR analyses with individual-level data from the UKB using a weighted genetic risk score (GRS); two-sample MR analyses were subsequently performed using summary-level data from the publicly available three consortia/biobank for CAD and two for HF. MR analyses were performed per database separately and results were subsequently meta-analysed using fixed-effects models. During a median follow-up of 11.8 years, restricted cubic spline Cox regression analyses showed associations between lower mtDNA abundance and higher risk of CAD and HF. Hazard ratios for participants in the lowest quintile of mtDNA abundance compared with those in the highest quintile were 1.08 (95% confidence interval: 1.03, 1.14) and 1.15 (1.05, 1.24) for CAD and HF. Genetically, no evidence was observed for a possible non-linear causal effect using individual-level weighted genetic risk scores calculated in the UKB on the study outcomes; the pooled odds ratios (95% confidence interval) from two-sample MR of genetically predicted per one-SD decrease in mtDNA abundance were 1.09 (1.03, 1.16) for CAD and 0.99 (0.92, 1.08) for HF, respectively. Conclusion: Our findings support a possible causal role of lower leukocyte mtDNA abundance in higher CAD risk, but not in HF. Show less
Luo, J.; Noordam, R.; Jukema, J.W.; Dijk, K.W. van; Hägg, S.; Grassmann, F.; ... ; Heemst, D. van 2022
AimMitochondrial DNA dysfunction has been implicated in the pathogenesis of cardiovascular diseases. We aimed to investigate the associations between leukocyte mitochondrial DNA (mtDNA) abundance,... Show moreAimMitochondrial DNA dysfunction has been implicated in the pathogenesis of cardiovascular diseases. We aimed to investigate the associations between leukocyte mitochondrial DNA (mtDNA) abundance, as a proxy of mitochondrial function, and coronary artery disease (CAD) and heart failure (HF) in a cohort study and approximate the causal nature of these relationships using Mendelian randomization (MR) in genetic studies.Methods and resultsMultivariable-adjusted Cox regression analyses were conducted in 273 619 unrelated participants of European ancestry from the UK Biobank (UKB). For genetic studies, we first performed MR analyses with individual-level data from the UKB using a weighted genetic risk score (GRS); two-sample MR analyses were subsequently performed using summary-level data from the publicly available three consortia/biobank for CAD and two for HF. MR analyses were performed per database separately and results were subsequently meta-analysed using fixed-effects models. During a median follow-up of 11.8 years, restricted cubic spline Cox regression analyses showed associations between lower mtDNA abundance and higher risk of CAD and HF. Hazard ratios for participants in the lowest quintile of mtDNA abundance compared with those in the highest quintile were 1.08 (95% confidence interval: 1.03, 1.14) and 1.15 (1.05, 1.24) for CAD and HF. Genetically, no evidence was observed for a possible non-linear causal effect using individual-level weighted genetic risk scores calculated in the UKB on the study outcomes; the pooled odds ratios (95% confidence interval) from two-sample MR of genetically predicted per one-SD decrease in mtDNA abundance were 1.09 (1.03, 1.16) for CAD and 0.99 (0.92, 1.08) for HF, respectively.ConclusionOur findings support a possible causal role of lower leukocyte mtDNA abundance in higher CAD risk, but not in HF. Show less
Maurits, M.P.; Wouters, F.; Niemantsverdriet, E.; Huizinga, T.W.J.; Akker, E.B. van den; Cessie, S. le; ... ; Knevel, R. 2022
Objective. To investigate whether established genetic predictors for rheumatoid arthritis (RA) differentiate healthy controls, patients with clinically suspect arthralgia (CSA), and RA patients... Show moreObjective. To investigate whether established genetic predictors for rheumatoid arthritis (RA) differentiate healthy controls, patients with clinically suspect arthralgia (CSA), and RA patients.Methods. Using analyses of variance, chi-square tests, and mean risk difference analyses, we investigated the association of an RA polygenic risk score (PRS) and HLA shared epitope (HLA-SE) with all participant groups, both unstratified and stratified for anti-citrullinated protein antibody (ACPA) status. We used 3 separate data sets sampled from the same Dutch population (1,015 healthy controls, 479 CSA patients, and 1,146 early classified RA patients). CSA patients were assessed for conversion to inflammatory arthritis over a period of 2 years, after which they were classified as either CSA converters (n = 84) or CSA nonconverters (n = 395).Results. The PRS was increased in RA patients (mean +/- SD PRS 1.31 +/- 0.96) compared to the complete CSA group (1.07 +/- 0.94) and compared to CSA converters (1.12 +/- 0.94). In ACPA- strata, PRS distributions differed strongly when comparing the complete CSA group (mean +/- SD PRS 1.05 +/- 0.94) and CSA converters (0.97 +/- 0.87) to RA patients (1.20 +/- 0.94), while in the ACPA+ strata, the complete CSA group (1.25 +/- 0.99) differed clearly from healthy controls (1.05 +/- 0.94) and RA patients (1.41 +/- 0.96). HLA-SE was more prevalent in the RA group (prevalence 0.64) than the complete CSA group (0.45), with small differences between RA patients and CSA converters (0.64 versus 0.60) and larger differences between CSA converters and CSA nonconverters (0.60 versus 0.42). HLA-SE prevalence differed more strongly within the ACPA+ strata as follows: healthy controls (prevalence 0.43), CSA nonconverters (0.48), complete CSA group (0.59), CSA converters (0.66), and RA patients (0.79).Conclusion. We observed that genetic predisposition increased across pre-RA participant groups. The RA PRS differed in early classified RA and inflammatory pre-disease stages, regardless of ACPA stratification. HLA-SE prevalence differed between arthritis patients, particularly ACPA+ patients, and healthy controls. Genetics seem to fulfill different etiologic roles. Show less
Maurits, M.P.; Wouters, F.; Niemantsverdriet, E.; Huizinga, T.W.J.; Akker, E.B. van den; Cessie, S. le; ... ; Knevel, R. 2022
ObjectiveTo investigate whether established genetic predictors for rheumatoid arthritis (RA) differentiate healthy controls, patients with clinically suspect arthralgia (CSA), and RA patients... Show moreObjectiveTo investigate whether established genetic predictors for rheumatoid arthritis (RA) differentiate healthy controls, patients with clinically suspect arthralgia (CSA), and RA patients.MethodsUsing analyses of variance, chi-square tests, and mean risk difference analyses, we investigated the association of an RA polygenic risk score (PRS) and HLA shared epitope (HLA-SE) with all participant groups, both unstratified and stratified for anti–citrullinated protein antibody (ACPA) status. We used 3 separate data sets sampled from the same Dutch population (1,015 healthy controls, 479 CSA patients, and 1,146 early classified RA patients). CSA patients were assessed for conversion to inflammatory arthritis over a period of 2 years, after which they were classified as either CSA converters (n = 84) or CSA nonconverters (n = 395).ResultsThe PRS was increased in RA patients (mean ± SD PRS 1.31 ± 0.96) compared to the complete CSA group (1.07 ± 0.94) and compared to CSA converters (1.12 ± 0.94). In ACPA– strata, PRS distributions differed strongly when comparing the complete CSA group (mean ± SD PRS 1.05 ± 0.94) and CSA converters (0.97 ± 0.87) to RA patients (1.20 ± 0.94), while in the ACPA+ strata, the complete CSA group (1.25 ± 0.99) differed clearly from healthy controls (1.05 ± 0.94) and RA patients (1.41 ± 0.96). HLA-SE was more prevalent in the RA group (prevalence 0.64) than the complete CSA group (0.45), with small differences between RA patients and CSA converters (0.64 versus 0.60) and larger differences between CSA converters and CSA nonconverters (0.60 versus 0.42). HLA-SE prevalence differed more strongly within the ACPA+ strata as follows: healthy controls (prevalence 0.43), CSA nonconverters (0.48), complete CSA group (0.59), CSA converters (0.66), and RA patients (0.79).ConclusionWe observed that genetic predisposition increased across pre-RA participant groups. The RA PRS differed in early classified RA and inflammatory pre-disease stages, regardless of ACPA stratification. HLA-SE prevalence differed between arthritis patients, particularly ACPA+ patients, and healthy controls. Genetics seem to fulfill different etiologic roles. Show less
Alshehri, T.; Mook-Kanamori, D.O.; Mutsert, R. de; Penninx, B.W.J.H.; Rosendaal, F.R.; Cessie, S. le; Milaneschi, Y. 2022
Background: Adiposity has been shown to be linked with atypical energy-related symptoms (AES) of depression. We used genomics to separate the effect of adiposity from that of metabolic... Show moreBackground: Adiposity has been shown to be linked with atypical energy-related symptoms (AES) of depression. We used genomics to separate the effect of adiposity from that of metabolic dysregulations to examine whether the link between obesity and AES is dependent on the presence of metabolic dysregulations. Method: Data were from NEO (n = 5734 individuals) and NESDA (n = 2238 individuals) cohorts, in which the Inventory of Depressive Symptomatology (IDS-SR30) was assessed. AES profile was based on four symptoms: increased appetite, increased weight, low energy level, and leaden paralysis. We estimated associations between AES and two genetic risk scores (GRS) indexing increasing total body fat with (metabolically unhealthy adiposity, GRS-MUA) and without (metabolically healthy adiposity, GRS-MHA) metabolic dysregulations. Results: We validated that both GRS-MUA and GRS-MHA were associated with higher total body fat in NEO study, but divergently associated with biomarkers of metabolic health (e.g., fasting glucose and HDL-cholesterol) in both cohorts. In the pooled results, per standard deviation, GRS-MUA was specifically associated with a higher AES score (beta = 0.03, 95%CI: 0.01; 0.05), while there was no association between GRS-MHA and AES (beta =-0.01, 95%CI:-0.03; 0.01). Conclusion: These results suggest that the established link between adiposity and AES profile emerges in the presence of metabolic dysregulations, which may represent the connecting substrate between the two conditions. Show less