Search

Abstract Background A prediction model can be a useful tool to quantify the risk of a patient developing dementia in the next years and take risk-factor-targeted intervention. Numerous dementia prediction models have been developed, but few have been externally validated, likely limiting their clinical uptake. In our previous work, we had limited success in externally validating some of these existing models due to inadequate reporting. As a result, we are compelled to develop and externally validate novel models to predict dementia in the general population across a network of observational databases. We assess regularization methods to obtain parsimonious models that are of lower complexity and easier to implement. Methods Logistic regression models were developed across a network of five observational databases with electronic health records (EHRs) and claims data to predict 5-year dementia risk in persons aged 55–84. The regularization methods L1 and Broken Adaptive Ridge (BAR) as well as three candidate predictor sets to optimize prediction performance were assessed. The predictor sets include a baseline set using only age and sex, a full set including all available candidate predictors, and a phenotype set which includes a limited number of clinically relevant predictors. Results BAR can be used for variable selection, outperforming L1 when a parsimonious model is desired. Adding candidate predictors for disease diagnosis and drug exposure generally improves the performance of baseline models using only age and sex. While a model trained on German EHR data saw an increase in AUROC from 0.74 to 0.83 with additional predictors, a model trained on US EHR data showed only minimal improvement from 0.79 to 0.81 AUROC. Nevertheless, the latter model developed using BAR regularization on the clinically relevant predictor set was ultimately chosen as best performing model as it demonstrated more consistent external validation performance and improved calibration. Conclusions We developed and externally validated patient-level models to predict dementia. Our results show that although dementia prediction is highly driven by demographic age, adding predictors based on condition diagnoses and drug exposures further improves prediction performance. BAR regularization outperforms L1 regularization to yield the most parsimonious yet still well-performing prediction model for dementia.

Abstract Background There is currently no consensus on the impact of class imbalance methods on the performance of clinical prediction models. We aimed to empirically investigate the impact of random oversampling and random undersampling, two commonly used class imbalance methods, on the internal and external validation performance of prediction models developed using observational health data. Methods We developed and externally validated prediction models for various outcomes of interest within a target population of people with pharmaceutically treated depression across four large observational health databases. We used three different classifiers (lasso logistic regression, random forest, XGBoost) and varied the target imbalance ratio. We evaluated the impact on model performance in terms of discrimination and calibration. Discrimination was assessed using the area under the receiver operating characteristic curve (AUROC) and calibration was assessed using calibration plots. Results We developed and externally validated a total of 1,566 prediction models. On internal and external validation, random oversampling and random undersampling generally did not result in higher AUROCs. Moreover, we found overestimated risks, although this miscalibration could largely be corrected by recalibrating the models towards the imbalance ratios in the original dataset. Conclusions Overall, we found that random oversampling or random undersampling generally does not improve the internal and external validation performance of prediction models developed in large observational health databases. Based on our findings, we do not recommend applying random oversampling or random undersampling when developing prediction models in large observational health databases.

Abstract Background Neuropsychiatric symptoms (NPS) are prevalent in the early clinical stages of Alzheimer’s disease (AD) according to proxy-based instruments. Little is known about which NPS clinicians report and whether their judgment aligns with proxy-based instruments. We used natural language processing (NLP) to classify NPS in electronic health records (EHRs) to estimate the reporting of NPS in symptomatic AD at the memory clinic according to clinicians. Next, we compared NPS as reported in EHRs and NPS reported by caregivers on the Neuropsychiatric Inventory (NPI). Methods Two academic memory clinic cohorts were used: the Amsterdam UMC (n = 3001) and the Erasmus MC (n = 646). Patients included in these cohorts had MCI, AD dementia, or mixed AD/VaD dementia. Ten trained clinicians annotated 13 types of NPS in a randomly selected training set of n = 500 EHRs from the Amsterdam UMC cohort and in a test set of n = 250 EHRs from the Erasmus MC cohort. For each NPS, a generalized linear classifier was trained and internally and externally validated. Prevalence estimates of NPS were adjusted for the imperfect sensitivity and specificity of each classifier. Intra-individual comparison of the NPS classified in EHRs and NPS reported on the NPI were conducted in a subsample (59%). Results Internal validation performance of the classifiers was excellent (AUC range: 0.81–0.91), but external validation performance decreased (AUC range: 0.51–0.93). NPS were prevalent in EHRs from the Amsterdam UMC, especially apathy (adjusted prevalence = 69.4%), anxiety (adjusted prevalence = 53.7%), aberrant motor behavior (adjusted prevalence = 47.5%), irritability (adjusted prevalence = 42.6%), and depression (adjusted prevalence = 38.5%). The ranking of NPS was similar for EHRs from the Erasmus MC, although not all classifiers obtained valid prevalence estimates due to low specificity. In both cohorts, there was minimal agreement between NPS classified in the EHRs and NPS reported on the NPI (all kappa coefficients

The spatial and frontal QRS-T angles are electrocardiographic (ECG) predictors for arrhythmia. This work used genetic analyses to identify associated loci and pathways, and explore their relationships with other ECG traits and cardiovascular disease.

Introduction: Monoclonal antibodies (mAbs) targeting immunoglobulin E (IgE) [omalizumab], type 2 (T2) cytokine interleukin (IL) 5 [mepolizumab, reslizumab], IL-4 Receptor (R) α [dupilumab], and IL-5R [benralizumab]), improve quality of life in patients with T2-driven inflammatory diseases. However, there is a concern for an increased risk of helminth infections. The aim was to explore safety signals of parasitic infections for omalizumab, mepolizumab, reslizumab, dupilumab, and benralizumab.Methods: Spontaneous reports were used from the Food and Drug Administration’s Adverse Event Reporting System (FAERS) database from 2004 to 2021. Parasitic infections were defined as any type of parasitic infection term obtained from the Standardised Medical Dictionary for Regulatory Activities® (MedDRA®). Safety signal strength was assessed by the Reporting Odds Ratio (ROR).Results: 15,502,908 reports were eligible for analysis. Amongst 175,888 reports for omalizumab, mepolizumab, reslizumab, dupilumab, and benralizumab, there were 79 reports on parasitic infections. Median age was 55 years (interquartile range 24–63 years) and 59.5% were female. Indications were known in 26 (32.9%) reports; 14 (53.8%) biologicals were reportedly prescribed for asthma, 8 (30.7%) for various types of dermatitis, and 2 (7.6%) for urticaria. A safety signal was observed for each biological, except for reslizumab (due to lack of power), with the strongest signal attributed to benralizumab (ROR = 15.7, 95% Confidence Interval: 8.4–29.3).Conclusion: Parasitic infections were disproportionately reported for mAbs targeting IgE, T2 cytokines, or T2 cytokine receptors. While the number of adverse event reports on parasitic infections in the database was relatively low, resulting safety signals were disproportionate and warrant further investigation.

Abstract Background Many dementia prediction models have been developed, but only few have been externally validated, which hinders clinical uptake and may pose a risk if models are applied to actual patients regardless. Externally validating an existing prediction model is a difficult task, where we mostly rely on the completeness of model reporting in a published article. In this study, we aim to externally validate existing dementia prediction models. To that end, we define model reporting criteria, review published studies, and externally validate three well reported models using routinely collected health data from administrative claims and electronic health records. Methods We identified dementia prediction models that were developed between 2011 and 2020 and assessed if they could be externally validated given a set of model criteria. In addition, we externally validated three of these models (Walters’ Dementia Risk Score, Mehta’s RxDx-Dementia Risk Index, and Nori’s ADRD dementia prediction model) on a network of six observational health databases from the United States, United Kingdom, Germany and the Netherlands, including the original development databases of the models. Results We reviewed 59 dementia prediction models. All models reported the prediction method, development database, and target and outcome definitions. Less frequently reported by these 59 prediction models were predictor definitions (52 models) including the time window in which a predictor is assessed (21 models), predictor coefficients (20 models), and the time-at-risk (42 models). The validation of the model by Walters (development c-statistic: 0.84) showed moderate transportability (0.67–0.76 c-statistic). The Mehta model (development c-statistic: 0.81) transported well to some of the external databases (0.69–0.79 c-statistic). The Nori model (development AUROC: 0.69) transported well (0.62–0.68 AUROC) but performed modestly overall. Recalibration showed improvements for the Walters and Nori models, while recalibration could not be assessed for the Mehta model due to unreported baseline hazard. Conclusion We observed that reporting is mostly insufficient to fully externally validate published dementia prediction models, and therefore, it is uncertain how well these models would work in other clinical settings. We emphasize the importance of following established guidelines for reporting clinical prediction models. We recommend that reporting should be more explicit and have external validation in mind if the model is meant to be applied in different settings.

The QT interval is a heritable electrocardiographic measure associated with arrhythmia risk when prolonged. Here, the authors used a series of genetic analyses to identify genetic loci, pathways, therapeutic targets, and relationships with cardiovascular disease.

Abstract Background HRV has mostly shown associations with systolic dysfunction and more recently, with diastolic dysfunction in Heart failure (HF) patients. But the role of sympathetic nervous system in changes of left ventricular (LV) systolic and diastolic function and new-onset HF has not been extensively studied. Methods Among 3157 men and 4405 women free of HF and atrial fibrillation retrospectively included from the population-based Rotterdam Study, we used linear mixed models to examine associations of RR-interval differences and standard deviation of RR-intervals corrected for heart rate (RMSSDc and SDNNc) with longitudinal changes of LV ejection fraction (LVEF), E/A ratio, left atrial (LA) diameter, E/e’ ratio. Afterwards, using cox regressions, we examined their association with new-onset HF. Results Mean (SD) age was 65 (9.95) in men and 65.7 (10.2) in women. Every unit increase in log RMSSDc was accompanied by 0.75% (95%CI:-1.11%;-0.39%) and 0.31% (− 0.60%;-0.01%) lower LVEF among men and women each year, respectively. Higher log RMSSDc was linked to 0.03 (− 0.04;-0.01) and 0.02 (− 0.03;-0.003) lower E/A and also − 1.76 (− 2.77;− 0.75) and − 1.18 (− 1.99;-0.38) lower LVM index in both sexes and 0.72 mm (95% CI: − 1.20;-0.25) smaller LA diameters in women. The associations with LVEF in women diminished after excluding HF cases during the first 3 years of follow-up. During a median follow-up of 8.7 years, hazard ratios (95%CI) for incident HF were 1.34 (1.08;1.65) for log RMSSDc in men and 1.15 (0.93;1.42) in women. SDNNc showed similar associations. Conclusions Indices of HRV were associated with worse systolic function in men but mainly with improvement in LA size in women. Higher HRV was associated with higher risk of new-onset HF in men. Our findings highlight potential sex differences in autonomic function underlying cardiac dysfunction and heart failure in the general population.

Abstract Background We investigated whether we could use influenza data to develop prediction models for COVID-19 to increase the speed at which prediction models can reliably be developed and validated early in a pandemic. We developed COVID-19 Estimated Risk (COVER) scores that quantify a patient’s risk of hospital admission with pneumonia (COVER-H), hospitalization with pneumonia requiring intensive services or death (COVER-I), or fatality (COVER-F) in the 30-days following COVID-19 diagnosis using historical data from patients with influenza or flu-like symptoms and tested this in COVID-19 patients. Methods We analyzed a federated network of electronic medical records and administrative claims data from 14 data sources and 6 countries containing data collected on or before 4/27/2020. We used a 2-step process to develop 3 scores using historical data from patients with influenza or flu-like symptoms any time prior to 2020. The first step was to create a data-driven model using LASSO regularized logistic regression, the covariates of which were used to develop aggregate covariates for the second step where the COVER scores were developed using a smaller set of features. These 3 COVER scores were then externally validated on patients with 1) influenza or flu-like symptoms and 2) confirmed or suspected COVID-19 diagnosis across 5 databases from South Korea, Spain, and the United States. Outcomes included i) hospitalization with pneumonia, ii) hospitalization with pneumonia requiring intensive services or death, and iii) death in the 30 days after index date. Results Overall, 44,507 COVID-19 patients were included for model validation. We identified 7 predictors (history of cancer, chronic obstructive pulmonary disease, diabetes, heart disease, hypertension, hyperlipidemia, kidney disease) which combined with age and sex discriminated which patients would experience any of our three outcomes. The models achieved good performance in influenza and COVID-19 cohorts. For COVID-19 the AUC ranges were, COVER-H: 0.69–0.81, COVER-I: 0.73–0.91, and COVER-F: 0.72–0.90. Calibration varied across the validations with some of the COVID-19 validations being less well calibrated than the influenza validations. Conclusions This research demonstrated the utility of using a proxy disease to develop a prediction model. The 3 COVER models with 9-predictors that were developed using influenza data perform well for COVID-19 patients for predicting hospitalization, intensive services, and fatality. The scores showed good discriminatory performance which transferred well to the COVID-19 population. There was some miscalibration in the COVID-19 validations, which is potentially due to the difference in symptom severity between the two diseases. A possible solution for this is to recalibrate the models in each location before use.

Abstract Background Knowledge graphs can represent the contents of biomedical literature and databases as subject-predicate-object triples, thereby enabling comprehensive analyses that identify e.g. relationships between diseases. Some diseases are often diagnosed in patients in specific temporal sequences, which are referred to as disease trajectories. Here, we determine whether a sequence of two diseases forms a trajectory by leveraging the predicate information from paths between (disease) proteins in a knowledge graph. Furthermore, we determine the added value of directional information of predicates for this task. To do so, we create four feature sets, based on two methods for representing indirect paths, and both with and without directional information of predicates (i.e., which protein is considered subject and which object). The added value of the directional information of predicates is quantified by comparing the classification performance of the feature sets that include or exclude it. Results Our method achieved a maximum area under the ROC curve of 89.8% and 74.5% when evaluated with two different reference sets. Use of directional information of predicates significantly improved performance by 6.5 and 2.0 percentage points respectively. Conclusions Our work demonstrates that predicates between proteins can be used to identify disease trajectories. Using the directional information of predicates significantly improved performance over not using this information.

On the electrocardiogram, the PR interval reflects conduction from the atria to ventricles and also serves as risk indicator of cardiovascular morbidity and mortality. Here, the authors perform genome-wide meta-analyses for PR interval in multiple ancestries and identify 141 previously unreported genetic loci.

Genome-wide association studies (GWAS) have identified many single nucleotide polymorphisms (SNPs) that play important roles in the genetic heritability of traits and diseases. With most of these SNPs located on the non-coding part of the genome, it is currently assumed that these SNPs influence the expression of nearby genes on the genome. However, identifying which genes are targeted by these disease-associated SNPs remains challenging. In the past, protein knowledge graphs have often been used to identify genes that are associated with disease, also referred to as “disease genes”. Here, we explore whether protein knowledge graphs can be used to identify genes that are targeted by disease-associated non-coding SNPs by testing and comparing the performance of six existing methods for a protein knowledge graph, four of which were developed for disease gene identification. We compare our performance against two baselines: (1) an existing state-of-the-art method that is based on guilt-by-association, and (2) the leading assumption that SNPs target the nearest gene on the genome. We test these methods with four reference sets, three of which were obtained by different means. Furthermore, we combine methods to investigate whether their combination improves performance. We find that protein knowledge graphs that include predicate information perform comparable to the current state of the art, achieving an area under the receiver operating characteristic curve (AUC) of 79.6% on average across all four reference sets. Protein knowledge graphs that lack predicate information perform comparable to our other baseline (genetic distance) which achieved an AUC of 75.7% across all four reference sets. Combining multiple methods improved performance to 84.9% AUC. We conclude that methods for a protein knowledge graph can be used to identify which genes are targeted by disease-associated non-coding SNPs.

Background: Fabry disease (FD) is an X-linked, lysosomal storage disorder leading to severe cardiomyopathy in a significant proportion of patients. To identify ECG markers that reflect early cardiac involvement and disease progression, we conducted a long term retrospective study in a large cohort of FD patients. Methods: A total of 1995 ECGs from 133 patients with classical FD (64% females, 80% treated with enzyme replacement therapy), spanning 20 years of follow-up, were compared to ECGs from 3893 apparently healthy individuals. Generalized linear mixed models were used to evaluate the effect of age, FD and sex on: P-wave duration, PR-interval, QRS-duration, QTc, Cornell index, spatial QRS-T angle and frontal QRS-axis. Regression slopes and absolute values for each parameter were compared between FD patients and control subjects. Results: At a younger age (p < 0.05). For the other ECG parameters, the rate of change, more than the absolute value, was greater in FD patients compared to controls (p < 0.05). From the fifth decade (men) or sixth (women) onwards, absolute values for P-wave duration, QRS-duration, QTc and spatial QRS-T angle were longer and higher in FD patients compared to control subjects. Conclusions: ECG abnormalities indicative of FD are age and sex dependent. Tracking the rate of change in ECG parameters could be a good way to detect disease progression, guiding treatment initiation. Moreover, monitoring ECG changes in FD can be used to evaluate the effectiveness of treatment.

Background Early prehospital recognition of critical conditions such as ST‐segment–elevation myocardial infarction (STEMI) has prognostic relevance. Current international electrocardiographic STEMI thresholds are predominantly based on individuals of Western European descent. However, because of ethnic electrocardiographic variability both in health and disease, there is a need to reevaluate diagnostic ST‐segment elevation thresholds for different populations. We hypothesized that fulfillment of ST‐segment elevation thresholds of STEMI criteria (STE‐ECGs) in apparently healthy individuals is ethnicity dependent. Methods and Results HELIUS (Healthy Life in an Urban Setting) is a multiethnic cohort study including 10 783 apparently healthy subjects of 6 different ethnicities (African Surinamese, Dutch, Ghanaian, Moroccan, South Asian Surinamese, and Turkish). Prevalence of STE‐ECGs across ethnicities, sexes, and age groups was assessed with respect to the 2 international STEMI thresholds: sex and age specific versus sex specific. Mean prevalence of STE‐ECGs was 2.8% to 3.4% (age/sex‐specific and sex‐specific thresholds, respectively), although with large ethnicity‐dependent variability. Prevalences in Western European Dutch were 2.3% to 3.0%, but excessively higher in young (45% STE‐ECGs. Conclusions The accuracy of diagnostic tests partly relies on background prevalence in healthy individuals. In apparently healthy subjects, there is a highly variable ethnicity‐dependent prevalence of ECGs with ST‐segment elevations exceeding STEMI thresholds. This has potential consequences for STEMI evaluations in individuals who are not of Western European descent, putatively resulting in adverse outcomes with both over‐ and underdiagnosis of STEMI.

Abstract Background Biomedical knowledge graphs have become important tools to computationally analyse the comprehensive body of biomedical knowledge. They represent knowledge as subject-predicate-object triples, in which the predicate indicates the relationship between subject and object. A triple can also contain provenance information, which consists of references to the sources of the triple (e.g. scientific publications or database entries). Knowledge graphs have been used to classify drug-disease pairs for drug efficacy screening, but existing computational methods have often ignored predicate and provenance information. Using this information, we aimed to develop a supervised machine learning classifier and determine the added value of predicate and provenance information for drug efficacy screening. To ensure the biological plausibility of our method we performed our research on the protein level, where drugs are represented by their drug target proteins, and diseases by their disease proteins. Results Using random forests with repeated 10-fold cross-validation, our method achieved an area under the ROC curve (AUC) of 78.1% and 74.3% for two reference sets. We benchmarked against a state-of-the-art knowledge-graph technique that does not use predicate and provenance information, obtaining AUCs of 65.6% and 64.6%, respectively. Classifiers that only used predicate information performed superior to classifiers that only used provenance information, but using both performed best. Conclusion We conclude that both predicate and provenance information provide added value for drug efficacy screening.

Abnormal PR interval duration is associated with risk for atrial fibrillation and heart block. Here, van Setten et al. identify 44 PR interval loci in a genome-wide association study of over 92,000 individuals and find genetic overlap with QRS duration, heart rate and atrial fibrillation.

Abstract Background Genome-wide association studies conducted on QRS duration, an electrocardiographic measurement associated with heart failure and sudden cardiac death, have led to novel biological insights into cardiac function. However, the variants identified fall predominantly in non-coding regions and their underlying mechanisms remain unclear. Results Here, we identify putative functional coding variation associated with changes in the QRS interval duration by combining Illumina HumanExome BeadChip genotype data from 77,898 participants of European ancestry and 7695 of African descent in our discovery cohort, followed by replication in 111,874 individuals of European ancestry from the UK Biobank and deCODE cohorts. We identify ten novel loci, seven within coding regions, including ADAMTS6, significantly associated with QRS duration in gene-based analyses. ADAMTS6 encodes a secreted metalloprotease of currently unknown function. In vitro validation analysis shows that the QRS-associated variants lead to impaired ADAMTS6 secretion and loss-of function analysis in mice demonstrates a previously unappreciated role for ADAMTS6 in connexin 43 gap junction expression, which is essential for myocardial conduction. Conclusions Our approach identifies novel coding and non-coding variants underlying ventricular depolarization and provides a possible mechanism for the ADAMTS6-associated conduction changes.

Abstract Background Electrocardiographic measures of left ventricular hypertrophy (LVH) are used as predictors of cardiovascular risk. We combined linkage and association analyses to discover novel rare genetic variants involved in three such measures and two principal components derived from them. Methods The study was conducted among participants from the Erasmus Rucphen Family Study (ERF), a Dutch family-based sample from the southwestern Netherlands. Variance components linkage analyses were performed using Merlin. Regions of interest (LOD > 1.9) were fine-mapped using microarray and exome sequence data. Results We observed one significant LOD score for the second principal component on chromosome 15 (LOD score = 3.01) and 12 suggestive LOD scores. Several loci contained variants identified in GWAS for these traits; however, these did not explain the linkage peaks, nor did other common variants. Exome sequence data identified two associated variants after multiple testing corrections were applied. Conclusions We did not find common SNPs explaining these linkage signals. Exome sequencing uncovered a relatively rare variant in MAPK3K11 on chromosome 11 (MAF = 0.01) that helped account for the suggestive linkage peak observed for the first principal component. Conditional analysis revealed a drop in LOD from 2.01 to 0.88 for MAP3K11, suggesting that this variant may partially explain the linkage signal at this chromosomal location. MAP3K11 is related to the JNK pathway and is a pro-apoptotic kinase that plays an important role in the induction of cardiomyocyte apoptosis in various pathologies, including LVH.

Heart rate variability (HRV) describes the individual variation in cardiac cycle duration and is a measure of vagal control of heart rate. Here, the authors identify seventeen single-nucleotide polymorphisms associated with HRV, lending new insight into the vagal regulation of heart rhythm.

Aims: QT variability is a promising electrocardiographic marker. It has been studied as a screening tool for coronary artery disease and left ventricular hypertrophy, and increased QT variability is a known risk factor for sudden cardiac death. Considering that comprehensive normal values for QT variability were lacking, we set out to establish these in standard 10-s electrocardiograms (ECGs) covering both sexes and all ages.Methods: Ten-second, 12-lead ECGs were provided by five Dutch population studies (Pediatric Normal ECG Study, Leiden University Einthoven Science Project, Prevention of Renal and Vascular End-stage Disease Study, Utrecht Health Project, Rotterdam Study). ECGs were recorded digitally and processed by well-validated analysis software. We selected cardiologically healthy participants, 46% being women. Ages ranged from 11 days to 91 years. After quality control, 13,828 ECGs were available. We assessed three markers: standard deviation of QT intervals (SDqt), short-term QT variability (STVqt), and QT variability index (QTVI).Results: For SDqt and STVqt, the median and the lower limit of normal remained stable with age. The upper limit of normal declined until around age 45, and increased strongly in the elderly, notably so in women. This implies that a subset of the population, small enough not to have appreciable effect on the median, shows a high degree of QT variability with a possible risk of arrhythmias or worse, especially in women. Otherwise, sex differences were negligible in all three measurements. For QTVI, median, and normal limits decreased until age 20, and steadily went up afterwards except for the lower limit of normal, which flattens off after age 65.Conclusion: We report the first set of normal values for QT variability based on 10-s ECGs, for all ages and both sexes.

eTRANSAFE is a research project funded within the Innovative Medicines Initiative (IMI), which aims at developing integrated databases and computational tools (the eTRANSAFE ToxHub) that support the translational safety assessment of new drugs by using legacy data provided by the pharmaceutical companies that participate in the project. The project objectives include the development of databases containing preclinical and clinical data, computational systems for translational analysis including tools for data query, analysis and visualization, as well as computational models to explain and predict drug safety events.

Purpose: Heart-rate variability (HRV) measured on standard 10-s electrocardiograms (ECGs) has been associated with increased risk of cardiac and all-cause mortality, but age- and sex-dependent normal values have not been established. Since heart rate strongly affects HRV, its effect should be taken into account. We determined a comprehensive set of normal values of heart-rate corrected HRV derived from 10-s ECGs for both children and adults, covering both sexes.Methods: Five population studies in the Netherlands (Pediatric Normal ECG Study, Leiden University Einthoven Science Project, Prevention of Renal and Vascular End-stage Disease Study, Utrecht Health Project, Rotterdam Study) provided 10-s, 12-lead ECGs. ECGs were stored digitally and analyzed by well-validated analysis software. We included cardiologically healthy participants, 42% being men. Their ages ranged from 11 days to 91 years. After quality control, 13,943 ECGs were available. Heart-rate correction formulas were derived using an exponential model. Two time-domain HRV markers were analyzed: the corrected standard deviation of the normal-to-normal RR intervals (SDNNc) and corrected root mean square of successive RR-interval differences (RMSSDc).Results: There was a considerable age effect. For both SDNNc and RMSSDc, the median and the lower limit of normal decreased steadily from birth until old age. The upper limit of normal decreased until the age of 60, but increased markedly after that age. Differences of the median were minimal between men and women.Conclusion: We report the first comprehensive set of normal values for heart-rate corrected 10-s HRV, which can be of value in clinical practice and in further research.

Nature Communications 8: Article number: 15805 (2017); Published: 14 June 2017; Updated: 2 August 2017 In Supplementary Fig. 10 of this Article, images for panels a and b were inadvertently omitted. The correct version of Supplementary Fig. 10 is provided as Supplementary Information associated withthis Erratum.