In order to investigate the applicability of routine 10s electrocardiogram (ECG) recordings for time-domain heart rate variability (HRV) calculation we explored to what extent these (ultra-)short ...recordings capture the "actual" HRV.
The standard deviation of normal-to-normal intervals (SDNN) and the root mean square of successive differences (RMSSD) were measured in 3,387 adults. SDNN and RMSSD were assessed from (ultra)short recordings of 10s(3x), 30s, and 120s and compared to 240s-300s (gold standard) measurements. Pearson's correlation coefficients (r), Bland-Altman 95% limits of agreement and Cohen's d statistics were used as agreement analysis techniques.
Agreement between the separate 10s recordings and the 240s-300s recording was already substantial (r = 0.758-0.764/Bias = 0.398-0.416/d = 0.855-0.894 for SDNN; r = 0.853-0.862/Bias = 0.079-0.096/d = 0.150-0.171 for RMSSD), and improved further when three 10s periods were averaged (r = 0.863/Bias = 0.406/d = 0.874 for SDNN; r = 0.941/Bias = 0.088/d = 0.167 for RMSSD). Agreement increased with recording length and reached near perfect agreement at 120s (r = 0.956/Bias = 0.064/d = 0.137 for SDNN; r = 0.986/Bias = 0.014/d = 0.027 for RMSSD). For all recording lengths and agreement measures, RMSSD outperformed SDNN.
Our results confirm that it is unnecessary to use recordings longer than 120s to obtain accurate measures of RMSSD and SDNN in the time domain. Even a single 10s (standard ECG) recording yields a valid RMSSD measurement, although an average over multiple 10s ECGs is preferable. For SDNN we would recommend either 30s or multiple 10s ECGs. Future research projects using time-domain HRV parameters, e.g. genetic epidemiological studies, could calculate HRV from (ultra-)short ECGs enabling such projects to be performed at a large scale.
Abstract
Summary
Summary statistics from a meta‐analysis of genome‐wide association studies (meta-GWAS) can be used for many follow-up analyses. One valuable application is the creation of polygenic ...scores. However, if polygenic scores are calculated in a validation cohort that was part of the meta-GWAS consortium, this cohort is not independent and analyses will therefore yield inflated results. The R package ‘MetaSubtract’ was developed to subtract the results of the validation cohort from meta‐GWAS summary statistics analytically. The statistical formulas for a meta‐analysis were inverted to compute corrected summary statistics of a meta‐GWAS leaving one (or more) cohort(s) out. These formulas have been implemented in MetaSubtract for different meta‐analyses methods (fixed effects inverse variance or square root sample size weighted z‐score) accounting for no, single or double genomic control correction. Results obtained by MetaSubtract correlate very well to those calculated using the traditional way, i.e. by performing a meta‐analysis leaving out the validation cohort. In conclusion, MetaSubtract allows researchers to compute meta‐GWAS summary statistics that are independent of the GWAS results of the validation cohort without requiring access to the cohort level GWAS results of the corresponding meta-GWAS consortium.
Availability and implementation
https://cran.r-project.org/web/packages/MetaSubtract.
Supplementary information
Supplementary data are available at Bioinformatics online.
For human complex traits, non-additive genetic variation has been invoked to explain “missing heritability,” but its discovery is often neglected in genome-wide association studies. Here we propose a ...method of using SNP data to partition and estimate the proportion of phenotypic variance attributed to additive and dominance genetic variation at all SNPs (hSNP2 and δSNP2) in unrelated individuals based on an orthogonal model where the estimate of hSNP2 is independent of that of δSNP2. With this method, we analyzed 79 quantitative traits in 6,715 unrelated European Americans. The estimate of δSNP2 averaged across all the 79 quantitative traits was 0.03, approximately a fifth of that for additive variation (average hSNP2 = 0.15). There were a few traits that showed substantial estimates of δSNP2, none of which were replicated in a larger sample of 11,965 individuals. We further performed genome-wide association analyses of the 79 quantitative traits and detected SNPs with genome-wide significant dominance effects only at the ABO locus for factor VIII and von Willebrand factor. All these results suggest that dominance variation at common SNPs explains only a small fraction of phenotypic variation for human complex traits and contributes little to the missing narrow-sense heritability problem.
Display omitted
•Unique genome-wide association study of multiple EDCs in 24-hour urine.•Four independent SNPs associated with EDCs through multi-trait GWAS analysis.•SNPs in CYP450 and SLC genes ...associated with MECPP and MEHHP.•CYP2C9 and SLC17A1 likely play causal roles in phthalate metabolism and excretion.
Ubiquitous exposure to environmental endocrine disrupting chemicals (EDCs) instigates a major public health problem, but much remains unknown on the inter-individual differences in metabolism and excretion of EDCs. To examine this we performed a two-stage genome-wide association study (GWAS) for 24-hour urinary excretions of four parabens, two bisphenols, and nine phthalate metabolites. Results showed five genome-wide significant (p-value < 5x10-8) and replicated single nucleotide polymorphisms (SNPs) representing four independent signals that associated with mono-(2-ethyl-5-carboxypentyl) phthalate (MECPP) and mono-(2-ethyl-5-hydroxyhexyl) phthalate (MEHHP). Three of the four signals were located on chromosome 10 in a locus harboring the cytochrome P450 (CYP) genes CYP2C9, CYP2C58P, and CYP2C19 (rs117529685, pMECPP = 5.38x10-25; rs117033379, pMECPP = 1.96x10-19; rs4918798, pMECPP = 4.01x10-71; rs7895726, pMEHHP = 1.37x10-15, r2 with rs4918798 = 0.93). The other signal was on chromosome 6 close to the solute carrier (SLC) genes SLC17A1, SLC17A3, SLC17A4, and SCGN (rs1359232, pMECPP = 7.6x10-16). These four SNPs explained a substantial part (8.3 % - 9.2 %) of the variance in MECPP in the replication cohort. Bioinformatics analyses supported a likely causal role of CYP2C9 and SLC17A1 in metabolism and excretion of MECPP and MEHHP. Our results provide biological insights into mechanisms of phthalate metabolism and excretion with a likely causal role for CYP2C9 and SLC17A1.
Obesity is a worldwide epidemic, with major health and economic costs. Here we estimate heritability for body mass index (BMI) in 172,000 sibling pairs and 150,832 unrelated individuals and explore ...the contribution of genotype-covariate interaction effects at common SNP loci. We find evidence for genotype-age interaction (likelihood ratio test (LRT) = 73.58, degrees of freedom (df) = 1, P = 4.83 × 10
), which contributed 8.1% (1.4% s.e.) to BMI variation. Across eight self-reported lifestyle factors, including diet and exercise, we find genotype-environment interaction only for smoking behavior (LRT = 19.70, P = 5.03 × 10
and LRT = 30.80, P = 1.42 × 10
), which contributed 4.0% (0.8% s.e.) to BMI variation. Bayesian association analysis suggests that BMI is highly polygenic, with 75% of the SNP heritability attributable to loci that each explain <0.01% of the phenotypic variance. Our findings imply that substantially larger sample sizes across ages and lifestyles are required to understand the full genetic architecture of BMI.
We examined whether the effect of maternal smoking during pregnancy on birthweight of the offspring was mediated by smoking-induced changes to DNA methylation in cord blood.
First, we used cord blood ...of 129 Dutch children exposed to maternal smoking vs 126 unexposed to maternal and paternal smoking (53% male) participating in the GECKO Drenthe birth cohort. DNA methylation was measured using the Illumina HumanMethylation450 Beadchip. We performed an epigenome-wide association study for the association between maternal smoking and methylation followed by a mediation analysis of the top signals false-discovery rate (FDR) < 0.05. We adjusted both analyses for maternal age, education, pre-pregnancy BMI, offspring's sex, gestational age and white blood cell composition. Secondly, in 175 exposed and 1248 unexposed newborns from two independent birth cohorts, we replicated and meta-analysed results of eight cytosine-phosphate-guanine (CpG) sites in the GFI1 gene, which showed the most robust mediation. Finally, we performed functional network and enrichment analysis.
We found 35 differentially methylated CpGs (FDR < 0.05) in newborns exposed vs unexposed to smoking, of which 23 survived Bonferroni correction (P < 1 × 10(-7)). These 23 CpGs mapped to eight genes: AHRR, GFI1, MYO1G, CYP1A1, NEUROG1, CNTNAP2, FRMD4A and LRP5. We observed partial confirmation as three of the eight CpGs in GFI1 replicated. These CpGs partly mediated the effect of maternal smoking on birthweight (Sobel P < 0.05) in meta-analysis of GECKO and the two replication cohorts. Differential methylation of these three GFI1 CpGs explained 12-19% of the 202 g lower birthweight in smoking mothers. Functional enrichment analysis pointed towards activation of cell-mediated immunity.
Maternal smoking during pregnancy was associated with cord blood methylation differences. We observed a potentially mediating role of methylation in the association between maternal smoking during pregnancy and birthweight of the offspring. Functional network analysis suggested a role in activating the immune system.
We propose a method (GREML-LDMS) to estimate heritability for human complex traits in unrelated individuals using whole-genome sequencing data. We demonstrate using simulations based on whole-genome ...sequencing data that ∼97% and ∼68% of variation at common and rare variants, respectively, can be captured by imputation. Using the GREML-LDMS method, we estimate from 44,126 unrelated individuals that all ∼17 million imputed variants explain 56% (standard error (s.e.) = 2.3%) of variance for height and 27% (s.e. = 2.5%) of variance for body mass index (BMI), and we find evidence that height- and BMI-associated variants have been under natural selection. Considering the imperfect tagging of imputation and potential overestimation of heritability from previous family-based studies, heritability is likely to be 60-70% for height and 30-40% for BMI. Therefore, the missing heritability is small for both traits. For further discovery of genes associated with complex traits, a study design with SNP arrays followed by imputation is more cost-effective than whole-genome sequencing at current prices.
Heavy metals are common in our environment, and all individuals are exposed to them to some extent. These toxic metals have several harmful effects on the body, including the kidney, which is a very ...sensitive organ. Indeed, heavy metal exposure has been linked to an increased risk of chronic kidney disease (CKD) and its progression, which may be explained by the well-established nephrotoxic effects of these metals. In this hypothesis and narrative literature review, we will shed light on the potential role that another highly common problem in patients with CKD, iron deficiency, may play in the damaging effects of heavy metal exposure in this patient group. Iron deficiency has previously been linked with an increased uptake of heavy metals in the intestine due to the upregulation of iron receptors that also take up other metals. Furthermore, recent research suggests a role of iron deficiency in the retention of heavy metals in the kidney. Therefore, we hypothesize that iron deficiency plays a crucial role in the damaging effects of heavy metal exposure in patients with CKD and that iron supplementation might be a strategy to combat these detrimental processes.
Iron deficiency (ID) is independently associated with an increased risk of death in renal transplant recipients (RTRs). ID promotes production and cleavage of intact fibroblast growth factor 23 ...(iFGF23) into C-terminal fibroblast growth factor 23 (cFGF23), elevated levels of which are also prospectively associated with adverse outcomes. We hypothesized that in RTRs, the relationship between ID and mortality is mediated by FGF23. We measured plasma iFGF23 and cFGF23 levels in 700 stable RTRs at a median of 5.4 years after transplant. RTRs with ID had median (interquartile range) cFGF23 concentrations higher than those of RTRs without ID (223 131-361 versus 124 88-180 RU/ml;
<0.001), whereas iFGF23 concentrations were similar between groups. In multivariable-adjusted Cox regression analyses, ID associated with increased mortality (81 events; hazard ratio, 1.95; 95% confidence interval, 1.22 to 3.10;
<0.01). However, this association lost significance after additional adjustment for cFGF23 levels (hazard ratio, 1.45; 95% confidence interval, 0.87 to 2.51;
=0.15). In further mediation analysis, cFGF23 explained 46% of the association between ID and mortality, whereas iFGF23 did not mediate this association. In conclusion, we found that cFGF23 levels are increased in iron-deficient RTRs and that the underlying biologic process driving production and cleavage of iFGF23, or alternatively the increased level of cFGF23 fragments, probably is an important mediator of the association between ID and mortality. Our results underline the strong relationship between iron and FGF23 physiology, and provide a potential mechanism explaining the relationship between ID and adverse outcome in RTRs.
Background The cause of heart failure with preserved ejection fraction (HFpEF) is poorly understood, and specific therapies are lacking. Previous studies suggested that inflammation plays a role in ...the development of HFpEF. Herein, we aimed to investigate in community-dwelling individuals whether a higher plasma interleukin 6 (IL-6) level is associated with an increased risk of developing new-onset heart failure (HF) over time, and specifically HFpEF. Methods and Results We performed a case-cohort study based on the PREVEND (Prevention of Renal and Vascular End-Stage Disease) study, a prospective general population-based cohort study. We included 961 participants, comprising 200 participants who developed HF and a random group of 761 controls. HF with reduced ejection fraction or HFpEF was defined on the basis of the left ventricular ejection fraction of ≤40% or >40%, respectively. In Cox proportional hazard regression analyses, IL-6 levels were statistically significantly associated with the development of HF (hazard ratio HR, 1.28; 95% CI, 1.02-1.61;
=0.03) after adjustment for key risk factors. Specifically, IL-6 levels were significantly associated with the development of HFpEF (HR, 1.59; 95% CI, 1.16-2.19;
=0.004), whereas the association with HF with reduced ejection fraction was nonsignificant (HR, 1.05; 95% CI, 0.75-1.47;
=0.77). In sensitivity analyses, defining HFpEF as left ventricular ejection fraction ≥50%, IL-6 levels were also significantly associated with the development of HFpEF (HR, 1.47; 95% CI, 1.04-2.06;
=0.03) after adjustment for key risk factors. Conclusions IL-6 is associated with new-onset HFpEF in community-dwelling individuals, independent of potential confounders. Our findings warrant further research to investigate whether IL-6 might be a novel treatment target to prevent HFpEF.