Bi-Gaussian analysis reveals distinct education-related alterations in spherical equivalent and axial length-results from the Gutenberg Health Study.

PURPOSE: The aim of this study is to investigate the distribution of spherical equivalent and axial length in the general population and to analyze the influence of education on spherical equivalent with a focus on ocular biometric parameters. METHODS: The Gutenberg Health Study is a population-based cohort study in Mainz, Germany. Participants underwent comprehensive ophthalmologic examinations as part of the 5-year follow-up examination in 2012-2017 including genotyping. The spherical equivalent and axial length distributions were modeled with gaussian mixture models. Regression analysis (on person-individual level) was performed to analyze associations between biometric parameters and educational factors. Mendelian randomization analysis explored the causal effect between spherical equivalent, axial length, and education. Additionally, effect mediation analysis examined the link between spherical equivalent and education. RESULTS: A total of 8532 study participants were included (median age: 57 years, 49% female). The distribution of spherical equivalent and axial length follows a bi-Gaussian function, partially explained by the length of education (i.e., < 11 years education vs. 11-20 years). Mendelian randomization indicated an effect of education on refractive error using a genetic risk score of education as an instrument variable (- 0.35 diopters per SD increase in the instrument, 95% CI, - 0.64-0.05, p = 0.02) and an effect of education on axial length (0.63 mm per SD increase in the instrument, 95% CI, 0.22-1.04, p = 0.003). Spherical equivalent, axial length and anterior chamber depth were associated with length of education in regression analyses. Mediation analysis revealed that the association between spherical equivalent and education is mainly driven (70%) by alteration in axial length. CONCLUSIONS: The distribution of axial length and spherical equivalent is represented by subgroups of the population (bi-Gaussian). This distribution can be partially explained by length of education. The impact of education on spherical equivalent is mainly driven by alteration in axial length.

Joint Linkage and Association Analysis Using GENEHUNTER-MODSCORE with an Application to Familial Pancreatic Cancer.

INTRODUCTION: Joint linkage and association (JLA) analysis combines two disease gene mapping strategies: linkage information contained in families and association information contained in populations. Such a JLA analysis can increase mapping power, especially when the evidence for both linkage and association is low to moderate. Similarly, an association analysis based on haplotypes instead of single markers can increase mapping power when the association pattern is complex. METHODS: In this paper, we present an extension to the GENEHUNTER-MODSCORE software package that enables a JLA analysis based on haplotypes and uses information from arbitrary pedigree types and unrelated individuals. Our new JLA method is an extension of the MOD score approach for linkage analysis, which allows the estimation of trait-model and linkage disequilibrium (LD) parameters, i.e., penetrance, disease-allele frequency, and haplotype frequencies. LD is modeled between alleles at a single diallelic disease locus and up to three diallelic test markers. Linkage information is contributed by additional multi-allelic flanking markers. We investigated the statistical properties of our JLA implementation using extensive simulations, and we compared our approach to another commonly used single-marker JLA test. To demonstrate the applicability of our new method in practice, we analyzed pedigree data from the German National Case Collection for Familial Pancreatic Cancer (FaPaCa). RESULTS: Based on the simulated data, we demonstrated the validity of our JLA-MOD score analysis implementation and identified scenarios in which haplotype-based tests outperformed the single-marker test. The estimated trait-model and LD parameters were in good accordance with the simulated values. Our method outperformed another commonly used JLA single-marker test when the LD pattern was complex. The exploratory analysis of the FaPaCa families led to the identification of a promising genetic region on chromosome 22q13.33, which can serve as a starting point for future mutation analysis and molecular research in pancreatic cancer. CONCLUSION: Our newly proposed JLA-MOD score method proves to be a valuable gene mapping and characterization tool, especially when either linkage or association information alone provide insufficient power to identify the disease-causing genetic variants.

GWAS of random glucose in 476,326 individuals provide insights into diabetes pathophysiology, complications and treatment stratification.

Conventional measurements of fasting and postprandial blood glucose levels investigated in genome-wide association studies (GWAS) cannot capture the effects of DNA variability on 'around the clock' glucoregulatory processes. Here we show that GWAS meta-analysis of glucose measurements under nonstandardized conditions (random glucose (RG)) in 476,326 individuals of diverse ancestries and without diabetes enables locus discovery and innovative pathophysiological observations. We discovered 120 RG loci represented by 150 distinct signals, including 13 with sex-dimorphic effects, two cross-ancestry and seven rare frequency signals. Of these, 44 loci are new for glycemic traits. Regulatory, glycosylation and metagenomic annotations highlight ileum and colon tissues, indicating an underappreciated role of the gastrointestinal tract in controlling blood glucose. Functional follow-up and molecular dynamics simulations of lower frequency coding variants in glucagon-like peptide-1 receptor (GLP1R), a type 2 diabetes treatment target, reveal that optimal selection of GLP-1R agonist therapy will benefit from tailored genetic stratification. We also provide evidence from Mendelian randomization that lung function is modulated by blood glucose and that pulmonary dysfunction is a diabetes complication. Our investigation yields new insights into the biology of glucose regulation, diabetes complications and pathways for treatment stratification.

Genetic insights into resting heart rate and its role in cardiovascular disease.

Resting heart rate is associated with cardiovascular diseases and mortality in observational and Mendelian randomization studies. The aims of this study are to extend the number of resting heart rate associated genetic variants and to obtain further insights in resting heart rate biology and its clinical consequences. A genome-wide meta-analysis of 100 studies in up to 835,465 individuals reveals 493 independent genetic variants in 352 loci, including 68 genetic variants outside previously identified resting heart rate associated loci. We prioritize 670 genes and in silico annotations point to their enrichment in cardiomyocytes and provide insights in their ECG signature. Two-sample Mendelian randomization analyses indicate that higher genetically predicted resting heart rate increases risk of dilated cardiomyopathy, but decreases risk of developing atrial fibrillation, ischemic stroke, and cardio-embolic stroke. We do not find evidence for a linear or non-linear genetic association between resting heart rate and all-cause mortality in contrast to our previous Mendelian randomization study. Systematic alteration of key differences between the current and previous Mendelian randomization study indicates that the most likely cause of the discrepancy between these studies arises from false positive findings in previous one-sample MR analyses caused by weak-instrument bias at lower P-value thresholds. The results extend our understanding of resting heart rate biology and give additional insights in its role in cardiovascular disease development.

Proteomics biomarker discovery for individualized prevention of familial pancreatic cancer using statistical learning.

BACKGROUND: The low five-year survival rate of pancreatic ductal adenocarcinoma (PDAC) and the low diagnostic rate of early-stage PDAC via imaging highlight the need to discover novel biomarkers and improve the current screening procedures for early diagnosis. Familial pancreatic cancer (FPC) describes the cases of PDAC that are present in two or more individuals within a circle of first-degree relatives. Using innovative high-throughput proteomics, we were able to quantify the protein profiles of individuals at risk from FPC families in different potential pre-cancer stages. However, the high-dimensional proteomics data structure challenges the use of traditional statistical analysis tools. Hence, we applied advanced statistical learning methods to enhance the analysis and improve the results' interpretability. METHODS: We applied model-based gradient boosting and adaptive lasso to deal with the small, unbalanced study design via simultaneous variable selection and model fitting. In addition, we used stability selection to identify a stable subset of selected biomarkers and, as a result, obtain even more interpretable results. In each step, we compared the performance of the different analytical pipelines and validated our approaches via simulation scenarios. RESULTS: In the simulation study, model-based gradient boosting showed a more accurate prediction performance in the small, unbalanced, and high-dimensional datasets than adaptive lasso and could identify more relevant variables. Furthermore, using model-based gradient boosting, we discovered a subset of promising serum biomarkers that may potentially improve the current screening procedure of FPC. CONCLUSION: Advanced statistical learning methods helped us overcome the shortcomings of an unbalanced study design in a valuable clinical dataset. The discovered serum biomarkers provide us with a clear direction for further investigations and more precise clinical hypotheses regarding the development of FPC and optimal strategies for its early detection.

Network reconstruction for trans acting genetic loci using multi-omics data and prior information.

BACKGROUND: Molecular measurements of the genome, the transcriptome, and the epigenome, often termed multi-omics data, provide an in-depth view on biological systems and their integration is crucial for gaining insights in complex regulatory processes. These data can be used to explain disease related genetic variants by linking them to intermediate molecular traits (quantitative trait loci, QTL). Molecular networks regulating cellular processes leave footprints in QTL results as so-called trans-QTL hotspots. Reconstructing these networks is a complex endeavor and use of biological prior information can improve network inference. However, previous efforts were limited in the types of priors used or have only been applied to model systems. In this study, we reconstruct the regulatory networks underlying trans-QTL hotspots using human cohort data and data-driven prior information. METHODS: We devised a new strategy to integrate QTL with human population scale multi-omics data. State-of-the art network inference methods including BDgraph and glasso were applied to these data. Comprehensive prior information to guide network inference was manually curated from large-scale biological databases. The inference approach was extensively benchmarked using simulated data and cross-cohort replication analyses. Best performing methods were subsequently applied to real-world human cohort data. RESULTS: Our benchmarks showed that prior-based strategies outperform methods without prior information in simulated data and show better replication across datasets. Application of our approach to human cohort data highlighted two novel regulatory networks related to schizophrenia and lean body mass for which we generated novel functional hypotheses. CONCLUSIONS: We demonstrate that existing biological knowledge can improve the integrative analysis of networks underlying trans associations and generate novel hypotheses about regulatory mechanisms.

Genetic analyses of the electrocardiographic QT interval and its components identify additional loci and pathways.

The QT interval is an electrocardiographic measure representing the sum of ventricular depolarization and repolarization, estimated by QRS duration and JT interval, respectively. QT interval abnormalities are associated with potentially fatal ventricular arrhythmia. Using genome-wide multi-ancestry analyses (>250,000 individuals) we identify 177, 156 and 121 independent loci for QT, JT and QRS, respectively, including a male-specific X-chromosome locus. Using gene-based rare-variant methods, we identify associations with Mendelian disease genes. Enrichments are observed in established pathways for QT and JT, and previously unreported genes indicated in insulin-receptor signalling and cardiac energy metabolism. In contrast for QRS, connective tissue components and processes for cell growth and extracellular matrix interactions are significantly enriched. We demonstrate polygenic risk score associations with atrial fibrillation, conduction disease and sudden cardiac death. Prioritization of druggable genes highlight potential therapeutic targets for arrhythmia. Together, these results substantially advance our understanding of the genetic architecture of ventricular depolarization and repolarization.

Common electrocardiogram measures are not associated with telomere length.

AIMS: Aging is accompanied by telomere shortening. Increased telomere shortening is considered a marker of premature aging. Cardiac aging results in the development of cardiac pathologies. Electrocardiogram (ECG) measures reflect cardiac excitation, conduction, and repolarization. ECG measures also prolong with aging and are associated with cardiac pathologies including atrial fibrillation. As premature prolongation of ECG measures is observed, we hypothesized that such prolongation may be associated with telomere length. METHODS AND RESULTS: We studied the large, community-based KORA F4 Study. Of 3,080 participants enrolled between 2006 and 2007 with detailed information on demographic, anthropometric, clinical, and ECG characteristics, 2,575 presented with available data on leukocyte telomere length. Telomere length was determined by real-time quantitative PCR and expressed relative to a single copy gene. We fitted multivariable adjusted linear regression models to associate the ECG measures RR-interval, PR-interval, QRS-duration, and heart rate corrected QTc with telomere length. In our cohort, the mean age was 54.9±12.9 years and 46.6% were men. Increased age was associated with shorter telomere length (p<0.01), and men had shorter telomere length than women (p<0.05). In unadjusted models, heart rate (p=0.023), PR-interval (p<0.01), and QTc-interval (p<0.01) were significantly associated with shorter telomere length. However, no significant associations remained after accounting for age, sex, and covariates. CONCLUSIONS: ECG measures are age-dependent, but not associated with shortened telomere length as a marker of biological aging. Further research is warranted to clarify if shortened telomeres are associated with clinical cardiac pathologies including atrial fibrillation.

ExomeChip-based rare variant association study in restless legs syndrome.

Restless legs syndrome (RLS) is a common sleep-related movement disorder in populations of European descent and disease risk is strongly influenced by genetic factors. Common variants have been assessed extensively in several genome-wide association studies, but the contribution of rarer genetic variation has not been investigated at this scale. We therefore genotyped a case-control set of 9246 individuals for mainly rare and low frequency exonic variants using the Illumina ExomeChip. However, standard single variant and gene-level association tests were negative. This does not preclude a role of rare variants in RLS, but is likely due to the small sample size and the limited selection of rare genetic variation captured on the array. Therefore, exome or whole genome sequencing should be performed rather than increasing the sample size of ExomeChip studies in order to identify rare risk variants for RLS.

Genome-wide meta-analysis of phytosterols reveals five novel loci and a detrimental effect on coronary atherosclerosis.

Phytosterol serum concentrations are under tight genetic control. The relationship between phytosterols and coronary artery disease (CAD) is controversially discussed. We perform a genome-wide meta-analysis of 32 phytosterol traits reflecting resorption, cholesterol synthesis and esterification in six studies with up to 9758 subjects and detect ten independent genome-wide significant SNPs at seven genomic loci. We confirm previously established associations at ABCG5/8 and ABO and demonstrate an extended locus heterogeneity at ABCG5/8 with different functional mechanisms. New loci comprise HMGCR, NPC1L1, PNLIPRP2, SCARB1 and APOE. Based on these results, we perform Mendelian Randomization analyses (MR) revealing a risk-increasing causal relationship of sitosterol serum concentrations and CAD, which is partly mediated by cholesterol. Here we report that phytosterols are polygenic traits. MR add evidence of both, direct and indirect causal effects of sitosterol on CAD.

Chronically elevated branched chain amino acid levels are pro-arrhythmic.

AIMS: Cardiac arrhythmias comprise a major health and economic burden and are associated with significant morbidity and mortality, including cardiac failure, stroke, and sudden cardiac death (SCD). Development of efficient preventive and therapeutic strategies is hampered by incomplete knowledge of disease mechanisms and pathways. Our aim is to identify novel mechanisms underlying cardiac arrhythmia and SCD using an unbiased approach. METHODS AND RESULTS: We employed a phenotype-driven N-ethyl-N-nitrosourea mutagenesis screen and identified a mouse line with a high incidence of sudden death at young age (6-9 weeks) in the absence of prior symptoms. Affected mice were found to be homozygous for the nonsense mutation Bcat2p.Q300*/p.Q300* in the Bcat2 gene encoding branched chain amino acid transaminase 2. At the age of 4-5 weeks, Bcat2p.Q300*/p.Q300* mice displayed drastic increase of plasma levels of branch chain amino acids (BCAAs-leucine, isoleucine, valine) due to the incomplete catabolism of BCAAs, in addition to inducible arrhythmias ex vivo as well as cardiac conduction and repolarization disturbances. In line with these findings, plasma BCAA levels were positively correlated to electrocardiogram indices of conduction and repolarization in the German community-based KORA F4 Study. Isolated cardiomyocytes from Bcat2p.Q300*/p.Q300* mice revealed action potential (AP) prolongation, pro-arrhythmic events (early and late afterdepolarizations, triggered APs), and dysregulated calcium homeostasis. Incubation of human pluripotent stem cell-derived cardiomyocytes with elevated concentration of BCAAs induced similar calcium dysregulation and pro-arrhythmic events which were prevented by rapamycin, demonstrating the crucial involvement of mTOR pathway activation. CONCLUSIONS: Our findings identify for the first time a causative link between elevated BCAAs and arrhythmia, which has implications for arrhythmogenesis in conditions associated with BCAA metabolism dysregulation such as diabetes, metabolic syndrome, and heart failure.

Cis-epistasis at the LPA locus and risk of cardiovascular diseases.

AIMS: Coronary artery disease (CAD) has a strong genetic predisposition. However, despite substantial discoveries made by genome-wide association studies (GWAS), a large proportion of heritability awaits identification. Non-additive genetic effects might be responsible for part of the unaccounted genetic variance. Here, we attempted a proof-of-concept study to identify non-additive genetic effects, namely epistatic interactions, associated with CAD. METHODS AND RESULTS: We tested for epistatic interactions in 10 CAD case-control studies and UK Biobank with focus on 8068 SNPs at 56 loci with known associations with CAD risk. We identified a SNP pair located in cis at the LPA locus, rs1800769 and rs9458001, to be jointly associated with risk for CAD [odds ratio (OR) = 1.37, P = 1.07 × 10-11], peripheral arterial disease (OR = 1.22, P = 2.32 × 10-4), aortic stenosis (OR = 1.47, P = 6.95 × 10-7), hepatic lipoprotein(a) (Lp(a)) transcript levels (beta = 0.39, P = 1.41 × 10-8), and Lp(a) serum levels (beta = 0.58, P = 8.7 × 10-32), while individual SNPs displayed no association. Further exploration of the LPA locus revealed a strong dependency of these associations on a rare variant, rs140570886, that was previously associated with Lp(a) levels. We confirmed increased CAD risk for heterozygous (relative OR = 1.46, P = 9.97 × 10-32) and individuals homozygous for the minor allele (relative OR = 1.77, P = 0.09) of rs140570886. Using forward model selection, we also show that epistatic interactions between rs140570886, rs9458001, and rs1800769 modulate the effects of the rs140570886 risk allele. CONCLUSIONS: These results demonstrate the feasibility of a large-scale knowledge-based epistasis scan and provide rare evidence of an epistatic interaction in a complex human disease. We were directed to a variant (rs140570886) influencing risk through additive genetic as well as epistatic effects. In summary, this study provides deeper insights into the genetic architecture of a locus important for cardiovascular diseases.

Common Genetic Variants Contribute to Risk of Transposition of the Great Arteries.

RATIONALE: Dextro-transposition of the great arteries (D-TGA) is a severe congenital heart defect which affects approximately 1 in 4,000 live births. While there are several reports of D-TGA patients with rare variants in individual genes, the majority of D-TGA cases remain genetically elusive. Familial recurrence patterns and the observation that most cases with D-TGA are sporadic suggest a polygenic inheritance for the disorder, yet this remains unexplored. OBJECTIVE: We sought to study the role of common single nucleotide polymorphisms (SNPs) in risk for D-TGA. METHODS AND RESULTS: We conducted a genome-wide association study in an international set of 1,237 patients with D-TGA and identified a genome-wide significant susceptibility locus on chromosome 3p14.3, which was subsequently replicated in an independent case-control set (rs56219800, meta-analysis P=8.6x10-10, OR=0.69 per C allele). SNP-based heritability analysis showed that 25% of variance in susceptibility to D-TGA may be explained by common variants. A genome-wide polygenic risk score derived from the discovery set was significantly associated to D-TGA in the replication set (P=4x10-5). The genome-wide significant locus (3p14.3) co-localizes with a putative regulatory element that interacts with the promoter of WNT5A, which encodes the Wnt Family Member 5A protein known for its role in cardiac development in mice. We show that this element drives reporter gene activity in the developing heart of mice and zebrafish and is bound by the developmental transcription factor TBX20. We further demonstrate that TBX20 attenuates Wnt5a expression levels in the developing mouse heart. CONCLUSIONS: This work provides support for a polygenic architecture in D-TGA and identifies a susceptibility locus on chromosome 3p14.3 near WNT5A. Genomic and functional data support a causal role of WNT5A at the locus.

Rare variant analysis in eczema identifies exonic variants in DUSP1, NOTCH4 and SLC9A4.

Previous genome-wide association studies revealed multiple common variants involved in eczema but the role of rare variants remains to be elucidated. Here, we investigate the role of rare variants in eczema susceptibility. We meta-analyze 21 study populations including 20,016 eczema cases and 380,433 controls. Rare variants are imputed with high accuracy using large population-based reference panels. We identify rare exonic variants in DUSP1, NOTCH4, and SLC9A4 to be associated with eczema. In DUSP1 and NOTCH4 missense variants are predicted to impact conserved functional domains. In addition, five novel common variants at SATB1-AS1/KCNH8, TRIB1/LINC00861, ZBTB1, TBX21/OSBPL7, and CSF2RB are discovered. While genes prioritized based on rare variants are significantly up-regulated in the skin, common variants point to immune cell function. Over 20% of the single nucleotide variant-based heritability is attributable to rare and low-frequency variants. The identified rare/low-frequency variants located in functional protein domains point to promising targets for novel therapeutic approaches to eczema.

Genetic Variation and Cardiovascular Risk Factors: A Cohort Study on Migrants from the Former Soviet Union and a Native German Population.

Resettlers are a large migrant group of more than 2 million people in Germany who migrated mainly from the former Soviet Union to Germany after 1989. We sought to compare the distribution of the major risk factors for cardiovascular disease (CVD) and to investigate the overall genetic differences in a study population which consisted of resettlers and native (autochthone) Germans. This was a joint analysis of two cohort studies which were performed in the region of Augsburg, Bavaria, Germany, with 3363 native Germans and 363 resettlers. Data from questionnaires and physical examinations were used to compare the risk factors for cardiovascular diseases between the resettlers and native Germans. A population-based genome-wide association analysis was performed in order to identify the genetic differences between the two groups. The distribution of the major risk factors for CVD differed between the two groups. The resettlers lead a less active lifestyle. While female resettlers smoked less than their German counterparts, the men showed similar smoking behavior. SNPs from three genes (BTNL2, DGKB, TGFBR3) indicated a difference in the two populations. In other studies, these genes have been shown to be associated with CVD, rheumatoid arthritis and osteoporosis, respectively.

A multi-ethnic epigenome-wide association study of leukocyte DNA methylation and blood lipids.

Here we examine the association between DNA methylation in circulating leukocytes and blood lipids in a multi-ethnic sample of 16,265 subjects. We identify 148, 35, and 4 novel associations among Europeans, African Americans, and Hispanics, respectively, and an additional 186 novel associations through a trans-ethnic meta-analysis. We observe a high concordance in the direction of effects across racial/ethnic groups, a high correlation of effect sizes between high-density lipoprotein and triglycerides, a modest overlap of associations with epigenome-wide association studies of other cardio-metabolic traits, and a largely non-overlap with lipid loci identified to date through genome-wide association studies. Thirty CpGs reached significance in at least 2 racial/ethnic groups including 7 that showed association with the expression of an annotated gene. CpGs annotated to CPT1A showed evidence of being influenced by triglycerides levels. DNA methylation levels of circulating leukocytes show robust and consistent association with blood lipid levels across multiple racial/ethnic groups.

Comparison of genetic risk prediction models to improve prediction of coronary heart disease in two large cohorts of the MONICA/KORA study.

It is still unclear how genetic information, provided as single-nucleotide polymorphisms (SNPs), can be most effectively integrated into risk prediction models for coronary heart disease (CHD) to add significant predictive value beyond clinical risk models. For the present study, a population-based case-cohort was used as a trainingset (451 incident cases, 1488 noncases) and an independent cohort as testset (160 incident cases, 2749 noncases). The following strategies to quantify genetic information were compared: A weighted genetic risk score including Metabochip SNPs associated with CHD in the literature (GRSMetabo ); selection of the most predictive SNPs among these literature-confirmed variants using priority-Lasso (PLMetabo ); validation of two comprehensive polygenic risk scores: GRSGola based on Metabochip data, and GRSKhera (available in the testset only) based on cross-validated genome-wide genotyping data. We used Cox regression to assess associations with incident CHD. C-index, category-free net reclassification index (cfNRI) and relative integrated discrimination improvement (IDIrel ) were used to quantify the predictive performance of genetic information beyond Framingham risk score variables. In contrast to GRSMetabo and PLMetabo , GRSGola significantly improved the prediction (delta C-index [95% confidence interval]: 0.0087 [0.0044, 0.0130]; IDIrel : 0.0509 [0.0131, 0.0894]; cfNRI improved only in cases: 0.1761 [0.0253, 0.3219]). GRSKhera yielded slightly worse prediction results than GRSGola .

Deoxyribonuclease 1 Q222R single nucleotide polymorphism and long-term mortality after acute myocardial infarction.

Upon activation, neutrophils release neutrophil extracellular traps (NETs), which contribute to circulating DNA burden and thrombosis, including ST-segment elevation myocardial infarction (STEMI). Deoxyribonuclease (DNase) 1 degrades circulating DNA and NETs. Lower DNase activity correlates with NET burden and infarct size. The DNase 1 Q222R single nucleotide polymorphism (SNP), impairing DNase 1 function, is linked with myocardial infarction. We assessed whether the Q222R SNP is connected to increased NET burden in STEMI and influences long-term outcomes. We enrolled 711 STEMI patients undergoing primary percutaneous coronary intervention (pPCI), and 1422 controls. Genotyping was performed for DNase 1 Q222R SNP. DNase activity, double-stranded (ds)DNA and citrullinated histone H3 were determined in culprit site and peripheral plasma during pPCI. The association of the Q222R variant on cardiovascular and all-cause mortality was assessed by multivariable Cox regression adjusted for cardiovascular risk factors. Homozygous Q222R DNase 1 variant was present in 64 (9.0%) STEMI patients, at the same frequency as in controls. Patients homozygous for Q222R displayed less DNase activity and increased circulating DNA burden. In overall patients, median survival was 60 months. Homozygous Q222R variant was independently associated with cardiovascular and all-cause mortality after STEMI. dsDNA/DNase ratio independently predicted cardiovascular and all-cause mortality. These findings highlight that the Q222R DNase 1 SNP is associated with increased NET burden and decreased compensatory DNase activity, and may serve as an independent risk factor for poor outcome after STEMI.