Low lung function has been associated with increased body mass index (BMI). The aim of this study was to investigate whether the effect of BMI on lung function is mediated by DNA methylation.
We used ...individual data from 285,495 participants in four population-based cohorts: the European Community Respiratory Health Survey, the Northern Finland Birth Cohort 1966, the Swiss Study on Air Pollution and Lung Disease in Adults, and the UK Biobank. We carried out Mendelian randomisation (MR) analyses in two steps using a two-sample approach with SNPs as instrumental variables (IVs) in each step. In step 1 MR, we estimated the causal effect of BMI on peripheral blood DNA methylation (measured at genome-wide level) using 95 BMI-associated SNPs as IVs. In step 2 MR, we estimated the causal effect of DNA methylation on FEV
, FVC, and FEV
/FVC using two SNPs acting as methQTLs occurring close (in cis) to CpGs identified in the first step. These analyses were conducted after exclusion of weak IVs (F statistic < 10) and MR estimates were derived using the Wald ratio, with standard error from the delta method. Individuals whose data were used in step 1 were not included in step 2.
In step 1, we found that BMI might have a small causal effect on DNA methylation levels (less than 1% change in methylation per 1 kg/m2 increase in BMI) at two CpGs (cg09046979 and cg12580248). In step 2, we found no evidence of a causal effect of DNA methylation at cg09046979 on lung function. We could not estimate the causal effect of DNA methylation at cg12580248 on lung function as we could not find publicly available data on the association of this CpG with SNPs.
To our knowledge, this is the first paper to report the use of a two-step MR approach to assess the role of DNA methylation in mediating the effect of a non-genetic factor on lung function. Our findings do not support a mediating effect of DNA methylation in the association of lung function with BMI.
BackgroundClassifying individuals at high chronic obstructive pulmonary disease (COPD)-risk creates opportunities for early COPD detection and active intervention.ObjectiveTo develop and validate a ...statistical model to predict 10-year probabilities of COPD defined by post-bronchodilator airflow obstruction (post-BD-AO; forced expiratory volume in 1 s/forced vital capacity<5th percentile).SettingGeneral Caucasian populations from Australia and Europe, 10 and 27 centres, respectively.ParticipantsFor the development cohort, questionnaire data on respiratory symptoms, smoking, asthma, occupation and participant sex were from the Tasmanian Longitudinal Health Study (TAHS) participants at age 41–45 years (n=5729) who did not have self-reported COPD/emphysema at baseline but had post-BD spirometry and smoking status at age 51–55 years (n=2407). The validation cohort comprised participants from the European Community Respiratory Health Survey (ECRHS) II and III (n=5970), restricted to those of age 40–49 and 50–59 with complete questionnaire and spirometry/smoking data, respectively (n=1407).Statistical methodRisk-prediction models were developed using randomForest then externally validated.ResultsArea under the receiver operating characteristic curve (AUCROC) of the final model was 80.8% (95% CI 80.0% to 81.6%), sensitivity 80.3% (77.7% to 82.9%), specificity 69.1% (68.7% to 69.5%), positive predictive value (PPV) 11.1% (10.3% to 11.9%) and negative predictive value (NPV) 98.7% (98.5% to 98.9%). The external validation was fair (AUCROC 75.6%), with the PPV increasing to 17.9% and NPV still 97.5% for adults aged 40–49 years with ≥1 respiratory symptom. To illustrate the model output using hypothetical case scenarios, a 43-year-old female unskilled worker who smoked 20 cigarettes/day for 30 years had a 27% predicted probability for post-BD-AO at age 53 if she continued to smoke. The predicted risk was 42% if she had coexistent active asthma, but only 4.5% if she had quit after age 43.ConclusionThis novel and validated risk-prediction model could identify adults aged in their 40s at high 10-year COPD-risk in the general population with potential to facilitate active monitoring/intervention in predicted ‘COPD cases’ at a much earlier age.
: MR-Egger regression has recently been proposed as a method for Mendelian randomization (MR) analyses incorporating summary data estimates of causal effect from multiple individual variants, which ...is robust to invalid instruments. It can be used to test for directional pleiotropy and provides an estimate of the causal effect adjusted for its presence. MR-Egger regression provides a useful additional sensitivity analysis to the standard inverse variance weighted (IVW) approach that assumes all variants are valid instruments. Both methods use weights that consider the single nucleotide polymorphism (SNP)-exposure associations to be known, rather than estimated. We call this the `NO Measurement Error' (NOME) assumption. Causal effect estimates from the IVW approach exhibit weak instrument bias whenever the genetic variants utilized violate the NOME assumption, which can be reliably measured using the F-statistic. The effect of NOME violation on MR-Egger regression has yet to be studied.
An adaptation of the I2 statistic from the field of meta-analysis is proposed to quantify the strength of NOME violation for MR-Egger. It lies between 0 and 1, and indicates the expected relative bias (or dilution) of the MR-Egger causal estimate in the two-sample MR context. We call it IGX2 . The method of simulation extrapolation is also explored to counteract the dilution. Their joint utility is evaluated using simulated data and applied to a real MR example.
In simulated two-sample MR analyses we show that, when a causal effect exists, the MR-Egger estimate of causal effect is biased towards the null when NOME is violated, and the stronger the violation (as indicated by lower values of IGX2 ), the stronger the dilution. When additionally all genetic variants are valid instruments, the type I error rate of the MR-Egger test for pleiotropy is inflated and the causal effect underestimated. Simulation extrapolation is shown to substantially mitigate these adverse effects. We demonstrate our proposed approach for a two-sample summary data MR analysis to estimate the causal effect of low-density lipoprotein on heart disease risk. A high value of IGX2 close to 1 indicates that dilution does not materially affect the standard MR-Egger analyses for these data.
: Care must be taken to assess the NOME assumption via the IGX2 statistic before implementing standard MR-Egger regression in the two-sample summary data context. If IGX2 is sufficiently low (less than 90%), inferences from the method should be interpreted with caution and adjustment methods considered.
Abstract
Background
Two-sample summary-data Mendelian randomization (MR) incorporating multiple genetic variants within a meta-analysis framework is a popular technique for assessing causality in ...epidemiology. If all genetic variants satisfy the instrumental variable (IV) and necessary modelling assumptions, then their individual ratio estimates of causal effect should be homogeneous. Observed heterogeneity signals that one or more of these assumptions could have been violated.
Methods
Causal estimation and heterogeneity assessment in MR require an approximation for the variance, or equivalently the inverse-variance weight, of each ratio estimate. We show that the most popular ‘first-order’ weights can lead to an inflation in the chances of detecting heterogeneity when in fact it is not present. Conversely, ostensibly more accurate ‘second-order’ weights can dramatically increase the chances of failing to detect heterogeneity when it is truly present. We derive modified weights to mitigate both of these adverse effects.
Results
Using Monte Carlo simulations, we show that the modified weights outperform first- and second-order weights in terms of heterogeneity quantification. Modified weights are also shown to remove the phenomenon of regression dilution bias in MR estimates obtained from weak instruments, unlike those obtained using first- and second-order weights. However, with small numbers of weak instruments, this comes at the cost of a reduction in estimate precision and power to detect a causal effect compared with first-order weighting. Moreover, first-order weights always furnish unbiased estimates and preserve the type I error rate under the causal null. We illustrate the utility of the new method using data from a recent two-sample summary-data MR analysis to assess the causal role of systolic blood pressure on coronary heart disease risk.
Conclusions
We propose the use of modified weights within two-sample summary-data MR studies for accurately quantifying heterogeneity and detecting outliers in the presence of weak instruments. Modified weights also have an important role to play in terms of causal estimation (in tandem with first-order weights) but further research is required to understand their strengths and weaknesses in specific settings.
Abstract
Background
With genome-wide association data for many exposures and outcomes now available from large biobanks, one-sample Mendelian randomization (MR) is increasingly used to investigate ...causal relationships. Many robust MR methods are available to address pleiotropy, but these assume independence between the gene-exposure and gene-outcome association estimates. Unlike in two-sample MR, in one-sample MR the two estimates are obtained from the same individuals, and the assumption of independence does not hold in the presence of confounding.
Methods
With simulations mimicking a typical study in UK Biobank, we assessed the performance, in terms of bias and precision of the MR estimate, of the fixed-effect and (multiplicative) random-effects meta-analysis method, weighted median estimator, weighted mode estimator and MR-Egger regression, used in both one-sample and two-sample data. We considered scenarios differing by the: presence/absence of a true causal effect; amount of confounding; and presence and type of pleiotropy (none, balanced or directional).
Results
Even in the presence of substantial correlation due to confounding, all two-sample methods used in one-sample MR performed similarly to when used in two-sample MR, except for MR-Egger which resulted in bias reflecting direction and magnitude of the confounding. Such bias was much reduced in the presence of very high variability in instrument strength across variants (IGX2 of 97%).
Conclusions
Two-sample MR methods can be safely used for one-sample MR performed within large biobanks, expect for MR-Egger. MR-Egger is not recommended for one-sample MR unless the correlation between the gene-exposure and gene-outcome estimates due to confounding can be kept low, or the variability in instrument strength is very high.
Summary
Background
There is uncertainty about the clinical usefulness of currently available asthma predictive tools. Validation of predictive tools in different populations and clinical settings is ...an essential requirement for the assessment of their predictive performance, reproducibility and generalizability. We aimed to critically appraise asthma predictive tools which have been validated in external studies.
Methods
We searched MEDLINE and EMBASE (1946‐2017) for all available childhood asthma prediction models and focused on externally validated predictive tools alongside the studies in which they were originally developed. We excluded non‐English and non‐original studies. PROSPERO registration number is CRD42016035727.
Results
From 946 screened papers, eight were included in the review. Statistical approaches for creation of prediction tools included chi‐square tests, logistic regression models and the least absolute shrinkage and selection operator. Predictive models were developed and validated in general and high‐risk populations. Only three prediction tools were externally validated: the Asthma Predictive Index, the PIAMA and the Leicester asthma prediction tool. A variety of predictors has been tested, but no studies examined the same combination. There was heterogeneity in definition of the primary outcome among development and validation studies, and no objective measurements were used for asthma diagnosis. The performance of tools varied at different ages of outcome assessment. We observed a discrepancy between the development and validation studies in the tools’ predictive performance in terms of sensitivity and positive predictive values.
Conclusions
Validated asthma predictive tools, reviewed in this paper, provided poor predictive accuracy with performance variation in sensitivity and positive predictive value.
Abbreviations: COPD, chronic obstructive pulmonary disease; NICE, National Institute for Health and Care Excellence; QALY, Quality-Adjusted Life Year; VOI, Value of Information Provenance: ...Commissioned; not externally peer reviewed At a time when the scale of investments has raised justifiable concerns about the ability of ongoing research to fulfill expectations 1, the long-run sustainability of research programs will depend on demonstration of value for money. ...if decision-making is about implementation of a cancer screening program, relevant parameters would include sensitivity and specificity of the screening test, cancer prevalence, and health benefits (relative to mortality, morbidity, and quality of life) and costs associated with the clinical pathway with and without the proposed program.
This paper provides guidelines for performing Mendelian randomization investigations. It is aimed at practitioners seeking to undertake analyses and write up their findings, and at journal editors ...and reviewers seeking to assess Mendelian randomization manuscripts. The guidelines are divided into nine sections: motivation and scope, data sources, choice of genetic variants, variant harmonization, primary analysis, supplementary and sensitivity analyses (one section on robust methods and one on other approaches), data presentation, and interpretation. These guidelines will be updated based on feedback from the community and advances in the field. Updates will be made periodically as needed, and at least every 18 months.