Using human genetics to understand the causes and consequences of circulating cardiac troponin I in the general population
==========================================================================================================================

* Marta R Moksnes
* Helge Røsjø
* Anne Richmond
* Magnus N Lyngbakken
* Sarah E Graham
* Ailin Falkmo Hansen
* Brooke N Wolford
* Sarah A Gagliano Taliun
* Jonathon LeFaive
* Humaira Rasheed
* Laurent F Thomas
* Wei Zhou
* Nay Aung
* Ida Surakka
* Nicholas J Douville
* Archie Campbell
* David J Porteous
* Steffen E. Petersen
* Patricia B. Munroe
* Paul Welsh
* Naveed Sattar
* George Davey Smith
* Lars G Fritsche
* Jonas B Nielsen
* Bjørn Olav Åsvold
* Kristian Hveem
* Caroline Hayward
* Cristen J Willer
* Ben M Brumpton
* Torbjørn Omland

## Abstract

Circulating cardiac troponin proteins are associated with structural heart disease and predict incident cardiovascular disease in the general population. However, the genetic contribution to cardiac troponin I (cTnI) concentrations and its causal effect on cardiovascular phenotypes is unclear. We combine data from two large population-based studies, the Trøndelag Health Study and the Generation Scotland Scottish Family Health Study and perform a genome-wide association study of high-sensitivity cTnI concentrations with 48 115 individuals. We further used two-sample Mendelian randomization to investigate the causal effects of circulating cTnI on acute myocardial infarction and heart failure.

We identified 12 genetic loci (8 novel) associated with cTnI concentrations. Associated protein-altering variants highlighted putative functional genes: *CAND2, HABP2, ANO5, APOH, FHOD3, TNFAIP2, KLKB1* and *LMAN1*. Phenome-wide association tests in 1283 phecodes and 274 continuous traits in UK Biobank showed associations between a polygenic risk score for cTnI and cardiac arrhythmias, aspartate aminotransferase 1 and anthropometric measures. Excluding individuals with a known history of comorbidities did not materially change associations with cTnI. Using two-sample Mendelian randomization we confirmed the non-causal role of cTnI in acute myocardial infarction (5 948 cases, 355 246 controls). We found some indications for a causal role of cTnI in heart failure (47 309 cases and 930 014 controls), but this was not supported by secondary analyses using left ventricular mass as outcome (18 257 individuals).

Our findings clarify the biology underlying the heritable contribution to circulating cTnI and support cTnI as a non-causal biomarker for acute myocardial infarction and heart failure development in the general population. Using genetically informed methods for causal inference of cTnI helps inform the role and value of measuring cTnI in the general population.

## Introduction

Cardiac specific troponins I (cTnI) and T (cTnT) are structural proteins involved in cardiac muscle contraction. During acute myocardial ischemia, cardiac troponin proteins are released into the blood stream and are used to diagnose acute myocardial infarction (AMI)(1). In AMI, cardiac troponin release is a result of direct ischemic cardiomyocyte injury. In contrast, the pathobiology of cardiac troponin release outside of AMI is not known, although low-level cardiac troponin elevations are associated with future cardiovascular events in the general population, leading some clinicians to suggest a need for wider use of cardiac troponins in risk stratification in general(2). Hence, a better understanding of low-level cardiac troponin release is needed, and the study of large population-based cohorts with genetic data could bring new insights.

The genetic determinants of cardiac troponin in the general population have not been studied in detail, and the downstream biology, including possible causal effects on cardiovascular disease, is unclear. The cTnI and cTnT proteins are encoded by the *TNNI3* and *TNNT2* genes(3), respectively, and mutations in these genes are recognized as disease-causing mutations in patients with familial cardiomyopathies(4,5). Patients with cardiomyopathies demonstrate increased concentrations of circulating cardiac troponins(6), and these conditions are strongly associated with heart failure (HF) development(7). We recently demonstrated a close association between cTnI and left ventricular hypertrophy(8), which represents an intermediate stage prior to symptomatic HF(9). Moreover, as experimental studies have found immune responses targeting cTnI to aggravate acute cardiac damage and to induce myocardial dysfunctions and HF(10,11), cTnI may have a direct or indirect causal effect on HF development in humans. Circulating cardiac troponin concentrations are also influenced by protein degradation and clearance(12). A previous genome-wide association study (GWAS) in the general population from the Generation Scotland Scottish Family Health Study (GS:SFHS) reported five genetic loci associated with cTnI concentrations and four other loci associated with cTnT concentrations(13), however without replication in an independent cohort. More knowledge about the genetic determinants of circulating cardiac troponins could reveal new pathobiology important for HF development. In contrast, as AMI type 1 is a result of an acute epicardial coronary artery occlusion and cTnI is only expressed in cardiac myocytes, no causal role for cTnI is expected for incident AMI. Due to limited knowledge of the genetic contribution to circulating cardiac troponins, no previous studies have used genetically informed methods of causal inference, namely Mendelian randomization (MR), to assess the causal effect of cardiac troponin concentrations on AMI or HF(14).

High-sensitivity cTnI has been cross-sectionally and repeatedly measured in three surveys of the Trøndelag Health Study(15) (HUNT) from 1995 to 2019. Additionally, recent efforts to sequence 53 831 diverse genomes in the NHLBI TOPMed Program(16) have enabled the imputation of more than 26 million variants in HUNT and GS:SFHS combined, including less common loss-of-function variants that have not been previously investigated(13).

In order to validate previous findings and discover novel genetic variants associated with circulating cTnI in the general population, we conducted a genome-wide meta-analysis of cTnI concentration in 48 115 individuals from HUNT and GS:SFHS. Additionally, we assessed changes in cTnI with time in 5 178 HUNT participants with repeated measures. Finally, we performed MR to investigate the causal effects of circulating cTnI on AMI and HF.

## Methods

### Cohort Descriptions

The recruitment and design of the HUNT and GS:SFHS studies have been reported elsewhere(15,17). Participants (aged ≥ 20 years in HUNT and 35–65 years in GS:SFHS) have been genotyped using Illumina HumanCoreExome array 12 v.1.0 or v.1.1 in HUNT and Illumina HumanOmniExpressExome-8 v1_A or v1-2 in GS:SFHS. Sample and variant quality control (QC) was performed using standard practices and has been reported in detail for both studies elsewhere(18,19). All variants were imputed from the TOPMed reference panel(16) (freeze 5) using Minimac4 v1.0 (see Web Resources). The reference panel is composed of 53 831 multi-ethnic samples and 308 107 085 single nucleotide variants (SNVs) and indel variants at high depth (mean read depth 38.2X). The individuals included in the analyses were of European ancestry.

### Cardiac Troponin I

High-sensitivity cTnI concentrations (Architect STAT, Abbott Laboratories, Abbott Diagnostics) were measured from serum samples in both studies using the manufacturers’ protocol and quality assurance scheme.

In HUNT, non-fasting serum samples were drawn in 1995-1997 (HUNT2), 2006-2008 (HUNT3) and 2017-2019 (HUNT4). The HUNT2 and HUNT3 samples were centrifuged at room temperature and stored at -80°C after collection. The HUNT2 samples were thawed and refrozen once in 2008 and subsequently stored at -20°C. High-sensitivity cTnI concentrations were analyzed in 2014 from HUNT2, in 2015 from HUNT3 and was analyzed consecutively in 2017-2019 from fresh serum samples in HUNT4. We combined cTnI measurements from all studies into one variable for cTnI concentration: If an individual had measurements from multiple HUNT surveys, we used the earliest measurement in the analyses, assuming participants to be less exposed to potential confounding factors at a younger age. After prioritization, 29 839 individuals of European ancestry had measurements of cTnI in HUNT, of which 23 990 were above or equal to the lower limit of detection (1.2 ng/L) reported for the assay(20). Measurements below 1.2 ng/L are imprecise on an individual level but were included since they collectively could add information to the GWAS analysis.

We defined change in cTnI concentration per time as the difference in cTnI concentration from HUNT2 to HUNT3 divided by the follow-up time. The median follow-up time was 10 years and 7 months. In total, 5 178 participants had data of change in cTnI over time.

In GS:SFHS, fasting blood samples were collected between 2006 and 2010. Serum samples were centrifuged, and serum aliquots stored at –80°C until future biochemical analyses. A total of 18 276 individuals had measurements of cTnI in GS:SFHS.

Measurements below the reported limit of detection (1.2 ng/L)(20) were reported as half of the limit of detection (0.6 ng/L) since the original measurements were unavailable.

### Association Analyses

We used a linear mixed model regression under an additive genetic model for each variant as implemented in BOLT-LMM (v.2.3.4)(21) to control for relatedness between the samples in all analyses.

We included age, sex, genotyping batch and the first 10 genetic principal components of ancestry (PCs) as covariates in both HUNT and GS:SFHS. Additionally, we included as a covariate in HUNT the survey from which the serum sample was collected for cTnI measurement. We applied rank-based inverse normal transformation for the cTnI concentration after adjusting for age and sex using linear regression (additionally adjusting for survey in HUNT). Variants with a minor allele count (MAC) < 10 or imputation R2 < 0.3 were excluded, leaving 21 233 433 variants in the HUNT analysis and 21 677 569 variants in the GS:SFHS analysis. We performed genomic control correction after each analysis.

### Meta-analysis

We performed a fixed-effect inverse-variance weighted meta-analysis of 26 711 992 variants imputed in up to 48 115 individuals using METAL(22). We considered genetic loci reaching a p-value < 5×10−8 in the meta-analysis for follow-up analysis.

### Sensitivity Analyses

To investigate if the identified genetic loci could be explained by clinical endpoints associated with chronically elevated cardiac troponin levels, we performed a sensitivity analysis of cTnI concentration in HUNT, excluding 5 253 individuals with a current or known history of comorbidities (atrial fibrillation [ICD9 427.3, ICD10 I48], diabetes mellitus [self-reported], HF[ICD9 428, ICD10 I50], cardiomyopathy [ICD9 425, ICD10 I42] or impaired kidney function [estimated glomerular Filtration rate [eGFR] < 60 ml/min/1.73m2] in the HUNT survey where the serum sample was collected for cTnI measurement).

### Heritability Estimation

We estimated the narrow-sense (additive) SNV heritability of cTnI concentration in HUNT using genome-wide complex trait analysis (GCTA) (23,24). We created a genetic relationship matrix (GRM) based on all genotyped autosomal variants in the 29 839 HUNT participants with troponin measurements and used the GRM with GCTA-GREML (genomic-relatedness-based restricted maximum-likelihood) to estimate the variance in cTnI concentration explained by the variants. Age, sex, genotyping batch and HUNT survey were used as covariates in the analysis, and the cardiac troponin measurement was transformed to normality with rank based inverse normal transformation after regression on age, sex and HUNT survey prior to analysis.

### Variant Annotation and Fine-mapping

We annotated genetic variants using ANNOVAR (v.2019Oct24)(25). For fine-mapping of the cTnI loci we applied SuSiE (v0.7.1)(26), which uses iterative Bayesian stepwise selection to estimate the number of causal variants in each locus and create credible sets with a 95% cumulative posterior probability for containing a causal variant. The LD matrix used with SuSiE was calculated from the imputed variants of 5 000 unrelated individuals in HUNT, and variants only imputed in GS:SFHS were therefore excluded from the fine-mapping. We visualized the fine-mapped regions with LocusZoom(27).

### Tissue, Gene and Gene Set Enrichment Tests

We performed gene set enrichment, gene prioritization and tissue/cell type enrichment tests on the cTnI loci using DEPICT (v1.0)(28). To ensure that DEPICT could identify variants in all loci we filtered the meta-analysis results to the 3 906 791 variants overlapping with the DEPICT reference files. We used plink (v1.9)(29) to identify two sets of independent variants by clumping the filtered results into 500 kb regions of correlated variants (r2 < 0.05) meeting two different p-value cut-offs (p-value < 1×10−5 and p-value < 5×10−8). The two sets of independent variants were used as input for DEPICT and analyzed separately. All cTnI loci were included in the analyses except variants in the HLA region.

### Expression Quantitative Trait Locus Tests

To identify cTnI index variants associated with the expression level of a gene in specific tissue types, we downloaded single tissue cis-eQTL data from the Genotype-Tissue Expression (GTEx) portal, data set v8 (see Web Resources). We identified all significant associations (p-value < 1×10−9 [5×10−8/49 tissues]) between cTnI index variants and gene expression levels in any tissue type. To further assess if any cis-eQTL signals and GWAS meta-analysis signals were consistent with a shared causal variant for cTnI concentration and gene expression, we used Bayesian colocalization analysis as implemented in the R package coloc(30). For each tissue type, we analyzed all genes whose expression was associated (p-value < 5×10−8) with at least one cTnI associated variant (p-value < 5×10−8), using effect sizes and standard errors for each variant-trait association as input. The coloc software estimated the variance in each trait (cTnI concentration or gene expression level) from the sample sizes and minor allele frequencies. We set the prior probability of a genetic variant being associated with only cTnI concentrations, only gene expression or both traits to be 10−4, 10−4 and 10−6 respectively. We considered posterior probabilities (PP4) above 75% to give strong support for a common causal variant for the two traits in a given tissue, posterior probability between 50% and 75% to be inconclusive or to give weak support for a common causal variant, and posterior probabilities < 50% to indicate no support for a common causal variant between the two traits.

### Polygenic Risk Score for cTnI

We used R (v3.6.2)(31) to construct a polygenic risk score (PRS) for cTnI, by summing the product of the effect size and the estimated allele count for the 12 index variants in genome-wide significant loci (p-values < 5×10−8). To test the association between the PRS and cTnI in an independent sample, we performed the GWAS on cTnI concentrations in 26 415 HUNT participants, excluding 3 424 individuals who were unrelated (no 3rd degree relatives or closer) both to each other and to the individuals in the GWAS. We used KING (v2.1.5)(32) to estimate the relatedness. The smaller GWAS in HUNT (N = 26 415) was meta-analyzed with the full GWAS from GS:SFHS to obtain effect estimates independent of the excluded individuals. We constructed the PRS from the effect size estimates from this meta-analysis and used a linear regression model to test the association between the PRS and serum concentrations of cTnI (inverse rank transformed to normality) in the independent sample of 3 424 unrelated HUNT participants.

### Phenome-wide Association Tests (PheWAS)

We tested the association of the PRS (PRS-PheWAS), as well as the individual index variants (single variant PheWAS), with 1 688 phecodes and 83 continuous biomarkers and traits in participants of white British ancestry in UK Biobank. We calculated a PRS for each participant as described in the above paragraph, using TOPMed imputed estimated allele counts and effect sizes from the meta-analysis. We used a logistic or linear regression model respectively to assess the association of the single variant estimated allele counts or inverse normalized PRS and each phecode or continuous biomarker. For the PRS-pheWAS we included as covariates sex and birth year for binary traits, and sex and age at measurement for continuous traits. For the single variant pheWAS we used GWAS summary statistics generated with SAIGE (v.29.4.2)(33), with sex and the first four PCs as covariates in addition to age at initial assessment for quantitative traits and birth year for binary traits. To correct for multiple testing of phecodes and continuous variables we used a Bonferroni corrected significance cut-off of 2.8×10−5.

### Mendelian Randomization of cTnI on AMI and HF

To assess causal associations of circulating cTnI on AMI or HF, we performed inverse variance weighted (IVW) two-sample MR analyses using R (v3.6.2)(31) with the TwoSampleMR(34) and MRInstruments(35) libraries. We used the effect sizes of 11 index variants from the meta-analysis as genetic instrument for cTnI, while we found the respective associations and effect sizes of the instrument variants for AMI (ICD10 I21) and HF (ICD10 I50) from independent genome-wide summary level data (5 948 AMI cases vs 355 246 controls from UK Biobank(36) available from the Neale laboratory (see Web Resources) and 47 309 HF cases vs 930 014 controls from the Heart Failure Molecular Epidemiology for Therapeutic Targets consortium [HERMES] consortium(37) available from the Cardiovascular Disease Knowledge Portal (see Web Resources). The rare cTnI index variant rs151313792 was excluded from both analyses because it was not present in either outcome data set, and we were unable to identify a proxy variant. In the analysis of HF, the variant rs8039472 was used as a proxy for rs8024538 (r2HUNT = 0.98) because rs8024538 was not present in the HF data. We performed sensitivity analyses using Steiger filtered IVW, weighted mode, weighted median, penalized weighted median and MR Egger. To evaluate the reverse causal association of HF on cTnI concentration, we also performed an IVW two-sample MR analysis using the effect sizes of 11 index variants for HF, and the respective associations and effect sizes of the same variants from our GWAS meta-analysis of cTnI. We further performed sensitivity analyses using left ventricular mass, which is associated with HF development. We found the associations and effect sizes of these variants for left ventricular mass from an independent GWAS in 18 257 participants in UK Biobank(38).

### Ethics

Participation in HUNT is based on informed consent, and the study has been approved by the Norwegian Data Protection Authority and the Regional Committee for Medical and Health Research Ethics in Central Norway (REK Reference number: 2015/2294).

The GS:SFHS obtained written informed consent from all participants and received ethical approval from the National Health Service Tayside Committee on Medical Research Ethics (REC Reference number: 05/S1401/89).

This study was covered by the ethics approval for UK Biobank studies (application 24460) from the NHS National Research Ethics Service on 17th June 2011 (Ref 11/NW/0382) and extended on 10th May 2016 (Ref 16/NW/0274).

## Results

### Discovery of Genetic Loci Associated with cTnI

After a discovery genome-wide association meta-analysis of HUNT (n = 29 839) and GS:SFHS (n = 18 276), we identified 12 genetic loci associated (p-value < 5×10−8) with cTnI concentrations (Tables 1, S1, S2, Figures S1, S2). All 12 loci were well imputed in both HUNT and GS:SFHS and had consistent directions and effect sizes across both studies. Eight of the 12 loci were previously unreported, while four were previously reported, but not yet replicated(13). Additionally, we discovered 11 novel, rare variants (MAF < 0.001, p-value < 5×10−8) that were only imputed in HUNT and could not therefore be replicated (Table S1). Excluding individuals with a known history of comorbidities did not materially change the results (Table S4). No genetic loci were genome-wide significant for changes in cTnI concentrations with time, however three loci were suggestive (p-value < 10−6) (Tables S3, S5).

View this table:
[Table 1.](http://medrxiv.org/content/early/2020/12/04/2020.09.04.20187401/T1)

Table 1. 
Index variants in genome-wide significant cTnI loci

We found three protein-altering index variants in cTnI loci: rs7080536 (MAF = 0.04, effect size = -0.10, p-value = 1×10−12) in *HABP2*, rs1801690 (MAF = 0.06, effect size = 0.10, p-value = 1×10−14) in *APOH*, and rs7481951 (MAF = 0.41, effect size = 0.05, p-value = 2×10−19) in *ANO5*. Furthermore, we identified a rare protein-altering index variant (rs151313792, MAF = 0.004, effect size = 0.26, p-value = 1×10−8) in *FHOD3*. Two additional index variants were in perfect linkage disequilibrium (LD) (r2=1.00) with protein-altering variants in *KLKB1* and *LMAN1*, as described in the section on fine-mapping below.

### Fine-mapping for Causal Variants in the cTnI Loci

We identified 14 credible sets (Table S9) of genetic variants in the cTnI loci using SuSiE. For each credible set there was a 95% cumulative posterior probability that it contained at least one causal variant. The sizes of the credible sets ranged from 1 variant to 64 variants (in the HLA region and one region on chromosome 10), providing varying evidence for the importance of specific variants. Nine of the cTnI index variants (rs7650482, rs12331618, rs7080536, rs7481951, rs4283165, rs8024538, rs1801690, rs151313792 and rs12604424) were found in a credible set, and one of these credible sets only contained the index variant itself (rs7080536), suggesting that this was the most likely causal variant in the locus. The two loci with the lowest p-values (chromosomes 10 and 11) gave rise to two credible sets each, indicating multiple independent signals. However, the credible sets were in close proximity or overlapping and therefore the seemingly independent signals might have resulted from uncaptured LD structures.

Annotation of the variants in each credible set (Table S10) revealed four additional protein-altering variants in moderate to high LD with non-exonic index variants (Table 1) : rs11718898 in *CAND2*, correlation r2 in HUNT (r2HUNT) = 0.91 with index variant rs7650482 (Figure 1a); rs3733402 in *KLKB1*, r2HUNT = 1.00 with index variant rs12331618 (Figure 1b); rs1132339 in *TNFAIP2*, r2HUNT = 0.65 with index variant rs4283165 (Figure 1c); rs2298711 in *LMAN1*, r2HUNT = 1.00 with index variant rs12604424 (Figure 1d). The rs2298711 protein-altering variant was reported as index variant in the previous GWAS study in GS:SFHS(13).

![Figure 1:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/12/04/2020.09.04.20187401/F1.medium.gif)

[Figure 1:](http://medrxiv.org/content/early/2020/12/04/2020.09.04.20187401/F1)

Figure 1: 
LocusZoom plots of cTnI loci with protein-altering variants in moderate or high LD with the index variant: Chromosome positions (GRCh38) are given on the x-axes, and the meta-analysis -log10(p-value) for the variants are given on the y-axis. T The genome-wide significance p-value threshold (5 x 10−8) is indicated by the dotted line. The correlation r2 for each variant is indicated by colors relative to the index variants (LD Ref Var) in each panel: Panel a: The CAND2 locus on chromosome 3 with index variant rs7650482. A protein-altering variant in CAND2, rs11718898 (pVal77Ala), is in high LD with the index variant. Panel b: The KLKB1 locus on chromosome 4 with index variant rs12331618. A protein-altering variant in KLKB1, rs3733402 (p.Ser143Asn, p.Ser105Asn), is in perfect LD with the index variant. Panel c: The TNFAIP2 locus on chromosome 14 with index variant rs4283165. A protein-altering variant in TNFAIP2, rs1132339 (p.Gln282Glu), is in moderate LD with the index variant. Panel d: The LMAN1 locus on chromosome 18 with index variant rs12604424. A protein-altering variant in LMAN1, rs2298711 (p.Met410Leu), is in perfect LD with the index variant.

Furthermore, we found a rare exonic frameshift substitution (rs137854521 in *ANO5*, r2HUNT = 0.003, D’ = 1.00 with the more common index variant rs7481951) in the second credible set in the same region as rs7481951, and a variant in another credible set was located in the *F12* 5’ untranslated region (UTR).

### Prioritization of Genes, Pathways and Tissues

To further explore the underlying biology of the cTnI associated loci, we used DEPICT to identify candidate causal genes at each locus. For suggestive (p-value < 1×10−5) cTnI loci (excluding the HLA region), we identified gene sets enriched for genes in cTnI loci, prioritized genes in each locus that shared predicted function with the genes in the remaining cTnI loci more often than expected by chance, and identified tissues and cell types in which cTnI locus genes were highly expressed (Tables S6-S8). After removing results with false discovery rate above 5%, we found an enriched expression of suggestive cTnI genes in heart and heart ventricle tissues, and in tooth development. Among 14 461 reconstituted gene sets, one gene set (“abnormal placenta morphology”) was enriched with genes in significant cTnI loci. DEPICT prioritized the genes *KLHL21, THAP3, CAND2, AOPEP, CTNNA3, SYNPO2L, VCL, ANXA11, MYOF, NRAP, CASP7, ANO5, ALPK3* and *ZNF592* based on shared predicted functions with genes from the other associated loci. Most of these (*CAND2, VCL, ANO5, ALPK3, NRAP, CASP7, SYNPO2L, ZNF592*) were in significant (p-value < 1×10−8) cTnI loci, but not always near the index variant.

### Heritability of cTnI Concentration

We estimated the proportion of phenotypic variance explained by all genotyped autosomal variants using GCTA and found a narrow-sense SNV heritability of cTnI concentration in HUNT (variance explained, Vg/Vp ± 1SE) to be 0.15 ± 0.01, p-value 2.63×10−32. The HUNT cohort has a high degree of relatedness and is therefore more genetically homogeneous than cohorts with fewer relatives. Including more rare variants from the TOPMed panel would further increase the homogeneity of the dataset. This might, in part, explain the lower heritability observed in HUNT than in previous studies(13,39).

### Variants Associated with Tissue-Specific cis-eQTLs

Five of the cTnI index variants were significantly associated with the expression of one or more genes in at least one tissue type (Table S11). Two of these variant-gene associations were found in the cardiovascular system: rs7481951 was associated with the expression of *ANO5* in aortic and tibial arteries and rs12331618 was associated with the expression of *ADK* in tibial arteries.

To discover shared causal variants for cTnI concentration and expression levels of nearby genes in different tissues, we used the cis-eQTL data from GTEx v8 and performed a Bayesian colocalization analysis of the gene expression signals and cTnI concentration signals. Four cTnI loci were significantly colocalized (posterior probability of common causal variant > 75%) with the cis-eQTL signals for a gene in at least one tissue (Table S12): The expression of coagulation factor *F11* in tibial arteries (rs12331618), the expression of coagulation factor *F12* in liver (rs2731672), the expression of *BMS1P4* (pseudogene) in whole blood (rs2270552), the expression of *ADK* in the sigmoid colon (rs2270552), as well as the expression of lncRNA gene *Lnc*-*WDR73*-10 in the left heart ventricle, mammary breast tissue and brain frontal cortex (rs8024538). The association signal for *ADK* was also near the posterior probability cut-off of 75% in aortic tissue (rs2270552) (Table S12).

### PheWAS of cTnI Loci and cTnI Polygenic Risk Score

Four of the index variants associated with increased cTnI concentrations were significantly associated with diseases or continuous traits in UK Biobank (Figures 2 and 3, Tables S13, S14) after correction for multiple testing: The variant rs7650482 was associated with cardiac arrhythmias, in particular atrial fibrillation and flutter, as well as waist hip ratio, standing height and basal metabolic rate. In addition, the variant rs12331618 was associated with pulmonary heart disease and phlebitis/thrombophlebitis, the variant rs2731672 was associated with standing height and basal metabolic rate, and the variant rs9391746 was mainly associated with autoimmune inflammatory conditions.

![Figure 2:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/12/04/2020.09.04.20187401/F2.medium.gif)

[Figure 2:](http://medrxiv.org/content/early/2020/12/04/2020.09.04.20187401/F2)

Figure 2: 
PRS-PheWAS: Associations of PRS (increasing cTnI concentrations) with 1 688 phecodes in UK Biobank.

![Figure 3:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/12/04/2020.09.04.20187401/F3.medium.gif)

[Figure 3:](http://medrxiv.org/content/early/2020/12/04/2020.09.04.20187401/F3)

Figure 3: 
PRS-PheWAS: Associations of PRS (increasing cTnI concentrations) with 83 continuous traits and biomarkers in UK Biobank. ECG: Electrocardiogram, MRI: Magnetic resonance imaging.

The polygenic risk score (PRS) based on genome-wide significant index variants (n = 12, variance explained=1.2%) was associated with seven phecodes and continuous traits in UK Biobank after correction for multiple testing (Figures 2, 3, Tables S15, S16): cardiac arrhythmias in general, and the sub-phenotype atrial fibrillation and flutter in particular, as well as standing height, aspartate aminotransferase 1, whole body fat free mass, whole body water mass and basal metabolic rate.

### Mendelian Randomization

To investigate the causal effect of circulating cTnI on both AMI and HF, we performed two-sample MR using 11 of 12 index variants from our meta-analysis as an instrument for cTnI (F-statistic = 44, r2 = 0.012 in the independent sample [N=34 24]), and summary-level data on AMI (5 948 AMI cases, 355 246 controls from UK Biobank) and HF (47 309 cases and 930 014 controls from the HERMES consortium).

There was no evidence of a causal effect of circulating cTnI on AMI (Inverse Variance Weighted (IVW) odds ratio (OR) = 1.00, 95% Confidence Interval (CI) [0.99, 1.00]).

Bidirectional MR however suggested evidence for a potential causal effect of cTnI on HF (IVW OR = 1.20, 95% CI: [1.00, 1.45]) in addition to the well-known effect of HF on cTnI concentration (IVW effect size = 0.14, 95% CI: [0.00, 0.28]) (Table S17). Sensitivity analyses showed similar results between all methods (weighted mode, weighted median, penalized weighted median, MR-Egger), although wider confidence intervals due to low power, and Steiger filtering did not change the results (Table S17, Figures 4-5). The MR-Egger intercept was close to zero for both outcomes, giving no strong evidence of direct effects of the index variants on AMI or HF that are not mediated by cTnI (horizontal pleiotropy). However, there was some evidence for heterogeneity in the analysis of HF (heterogeneity statistic Q = 27, p-value = 0.003), but not for the analysis of AMI (Q = 9, p-value = 0.50). The heterogeneity in the HF analysis could possibly reflect horizontal pleiotropy, which would violate the instrumental variable analysis assumptions, or indicate multiple underlying mechanisms. Because HF development is commonly preceded by increase in left ventricular mass, we further investigated the association of cTnI on HF. We found a weak association of circulating cTnI on left ventricular mass (IVW effect size = 0.07, 95% CI: -0.11, 0.27), but the effect was attenuated and changed direction in sensitivity analyses (Table S17).

![Figure 4:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/12/04/2020.09.04.20187401/F4.medium.gif)

[Figure 4:](http://medrxiv.org/content/early/2020/12/04/2020.09.04.20187401/F4)

Figure 4: 
Mendelian randomization analysis of the causal effect of cTnI on AMI: The analysis was based on summary statistics from the meta-analysis of HUNT and GS:SFHS (n = 48 115) for cTnI and from a GWAS of AMI in UK Biobank from the Neale laboratory (cases/controls= 5 9948/355 246). Effect allele-cTnI associations given on x-axis and effect allele-AMI associations given on y-axis. All associations are given with 1 standard deviation error bars.

![Figure 5:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/12/04/2020.09.04.20187401/F5.medium.gif)

[Figure 5:](http://medrxiv.org/content/early/2020/12/04/2020.09.04.20187401/F5)

Figure 5: 
Mendelian randomization analysis of the causal effect of cTnI on HF: The analysis was based on summary statistics from the meta-analysis of HUNT and GS:SFHS (n = 48 115) for cTnI and from a meta-analysis of HF in the HERMES consortium (cases/controls = 47 309/930 014). Effect allele-cTnI associations given on x-axis and effect allele-HF associations given on y-axis. All associations are given with 1 standard deviation error bars. The outlier variant with a large effect size for both traits is rs9391746 from the HLA region.

## Discussion

Our GWAS meta-analysis of 48 115 individuals in HUNT and GS:SFHS revealed a genetic contribution to circulating protein concentrations of cTnI with associations to several biological processes. We identified 12 loci consistently associated with cTnI concentrations across both studies, confirming four of the five previously unreplicated loci(13). Furthermore, we suggest several novel loci for change in cTnI with time, confirm cTnI as a non-causal biomarker for AMI and observe indications that cTnI may have a causal effect on HF development.

Cardiac troponin I and T are important biomarkers for cardiovascular disease, and especially used to diagnose AMI in patients with symptoms of acute myocardial ischemia(1,40). Moreover, cTnI and cTnT have been found to provide strong and independent information regarding future risk for cardiac morbidity and mortality among subjects of the general population(41). Of note, while cardiac troponins are excellent markers of AMI in the acute setting, high-sensitivity troponins seem more related to incident heart failure events than to future AMI in the general population(42,43). Hence, data from several independent cohorts support a model of cardiac troponins as either directly or indirectly associated with heart failure development. This model is supported by high-sensitivity troponin concentrations being correlated with indices of structural heart disease, and most prominently left ventricular mass(42,44–46).

Using genetic data from two large cohorts from the general population (HUNT and GS:SFHS), we now provide new and additional information related to a role for cardiac troponins in cardiovascular disease. Four protein-altering index variants have previously been associated with cardiovascular health or with muscle development and injury: Firstly, a protein-altering variant in the fibrinolysis related gene *HABP2* (rs7080536) that reduces the activity of the translated protein(47), coagulation factor VII activating protease, and has been associated with carotid stenosis(48) and venous thromboembolism(49). Secondly, a protein-altering variant in *APOH* (rs1801690) that has been associated with the concentrations of creatine kinase(50), which is an unspecific biomarker for acute muscle injury, and with activated partial thromboplastin time(50), a measure of coagulation. Third, a previously reported protein-altering variant in *ANO5* (rs7481951) that has been suggested to be causal for adult-onset muscular dystrophy(51). Finally, we identified a rare protein-altering variant (rs151313792) in *FHOD3*, a gene that regulates actin assembly in cardiac sarcomeres(52) and has been associated with dilated cardiomyopathy(53) and left ventricular ejection fraction(54).

The remaining index variants were in close proximity to genes previously associated with the cardiovascular system, either through disease associations or links to cardiovascular development, structure and function: *ALPK3, FHOD3, VCL* and *CAND2* have been implicated in heart development and/or function(50,55–57), and *TNFAIP2* is central in blood vessel formation(58). Both *ALPK3* and *VCL* are associated with cardiomyopathy(4,53,55), and *CAND2* is associated with atrial fibrillation(59) and also creatine kinase concentrations(50). Furthermore, the *KLKB1*/*CYP4V2* locus has been associated with atherosclerosis(60), and the *KLKB1* and *F12* genes are both involved in several proteolytic reaction cascades in the cardiovascular system(61). These cascades are associated with coagulation/fibrinolysis, blood pressure regulation and inflammation(61). *LMAN1* has also been associated with coagulation(62), and both *TNFAIP2* and *HLA*-*DQB1* are associated with vascular inflammation(58,63). Both coagulation and atherosclerosis are pathobiological components crucial to the progress of acute exacerbations of coronary artery disease, and blood pressure regulation will impact left ventricular mass, and this is therefore comparable with previous data on cardiac troponins and left ventricular structure and function. The *F12*/*GRK6* locus has additionally been associated with chronic kidney disease(64). Several of the protein-altering and regulatory region variants that we identified had previously also been linked to some of the above-mentioned processes(47–51,60,65,66). Excluding individuals with impaired kidney function or a history of diabetes, atrial fibrillation and flutter, cardiomyopathy or HF, did not attenuate the effect sizes of any of the 12 index variants, supporting that the genetic associations with cTnI concentrations were not driven by these clinical endpoints.

Fine-mapping of the cTnI associated loci pointed to additional protein-altering and regulatory region variants, and further indicated two separate causal signals in each of the two most significant cTnI loci. A rare *ANO5* frameshift substitution (rs137854521) was possibly driving the second independent signal in the most significant locus as indicated by the fine-mapping analysis. A variant in another credible set was located in the *F12* 5’ untranslated region (UTR) and is associated with coagulation factor XII deficiency (rs1801020, r2HUNT = 0.99 with rs2731672). The protein-altering variant in *HABP2* was alone in its credible set, providing additional support for its contribution to circulating cTnI concentrations in the general population. On the other hand, a large credible set from the barely significant HLA region contained only intergenic variants, none of which were the index variant, weakening the evidence for a role of this specific variant in troponin biology. This is not surprising, given that variants in this region are known to be difficult to impute and interpret(67).

With a small number of significant cTnI loci, the DEPICT analysis was limited by the variant database. We filtered on variants in the DEPICT database so that all cTnI loci were represented, but with the cost that the results were not necessarily based on the most significant variants within the loci. Despite this limitation, we found that the expression of genes in cTnI loci with p-value < 10−5 was most significantly enriched in the heart (particularly in the ventricles), which is as expected given that cTnI is only expressed in cardiac myocytes, which is also the basis for the excellent performance of cardiac troponins as cardiovascular biomarkers(68). Additionally, DEPICT prioritized several genes in the cTnI loci that were associated with heart development, structure and function, and with cardiovascular diseases: *NRAP* is involved in sarcomere assembly during cardiomyocyte development(69) and is associated with familial cardiomyopathies(70), while *CASP7* is involved in activating apoptosis(71) and is associated with cardiac electrophysiology(71). Both are located in the *HABP2* locus. Likewise, *SYNPO2L*, which is involved in sarcomeric signal transduction(72) and cardiac development(37), is located in the second most significant locus on chromosome 10, but further away from the index variant than *VCL*. The associations with sarcomere structure and function are plausible based on the biology as cardiac troponins are essential components of the sarcomere.

While we identified protein-altering variants in more than half of the loci, several statistically significant variants were located outside of the coding regions. Non-coding variants might however have regulatory functions, and we therefore investigated their associations with variation in gene expression across 49 tissues. Four cTnI loci with non-coding index variants were highly colocalized with the expression signals for a nearby gene in at least one tissue, including the expression of coagulation factors in arteries and liver, and a protein degradation activator in blood. The colocalization with expression signals for adenosine kinase was also quite strong in the aorta. This protein has been implicated in cardiomyocyte microtubule dynamics and cardiac adaption to high blood pressure and is found to attenuate cardiac hypertrophy(73). These findings corroborate the role of these loci in cTnI biology, possibly in relation to the response to cardiac injury and myocardial stress, and to the breakdown of circulating cTnI.

The single-variant phenome-wide association tests in the UK Biobank indicated some level of pleiotropy of cTnI variants: Seven of 12 index variants were associated with other traits than cTnI. These traits included cardiovascular diseases, waist-hip ratio, standing height, aspartate aminotransferase 1 concentrations and basal metabolic rate. Waist-hip ratio has been linked to future risk of cardiovascular outcomes(74), standing height correlates with blood pressure, cardiac output and vascular resistance(75), and the basal metabolic rate could potentially affect cTnI degradation. Aspartate aminotransferase 1 is an enzyme released from the heart in patients with AMI(76), but is also related to amino acid metabolism in the liver and is used as an unspecific biomarker for liver cell damage(77). This enzyme could thus both be related to cardiovascular disease and possibly to degradation and clearance of circulating cTnI in the liver. However, there were no associations between cTnI variants and the more liver-specific enzyme alanine aminotransferase 1, which could indicate that the association is related to the cardiovascular system, rather than to liver damage. Hence, the mechanism whereby aspartate aminotransferase 1 influences circulating cTnI concentrations will require additional studies. The identified HLA region variant in our study was mainly associated with autoimmune and inflammatory diseases, and this is relevant for cardiovascular diseases(78), but more studies are needed, especially since interpretation of specific HLA region variants is difficult.

The PRS combining the 12 cTnI index variants explained 1.2% of variation in cTnI concentrations in an independent sample from HUNT. In UK Biobank it was associated with some of the same traits as the single variants, which is not surprising: standing height, basal metabolic rate, aspartate aminotransferase 1 concentrations and cardiac arrhythmias, including atrial fibrillation and flutter.

Furthermore, the PRS was associated with whole body fat-free mass and whole-body water mass which was not seen for any single variant. By performing two-sample MR, we confirmed a non-causal role of cTnI in AMI. In contrast, we found indications of a potential causal effect of cTnI in HF development, although the estimate was imprecise. This was surprising given the current opinion of cTnI as a passive biomarker of cardiac injury(68), but in line with some animal studies suggesting that circulating cTnI can trigger immune responses causing HF(10,11). Follow-up analyses of left ventricular mass, which is known to increase during HF development, however, did not support this finding.

In conclusion, our study highlights diverse genetic determinants of cTnI concentration in the general population and suggests cTnI as a non-causal biomarker for AMI and HF development in the general population. Furthermore, it suggests several novel loci for change in cTnI over time, although future studies are required to validate these findings and to explore possible mechanisms.

## Supporting information

Supplemental Tables S1, S5, S6, S9, S10, S12 [[supplements/187401_file03.xlsx]](pending:yes)

Supplemental Figures S1, S2. Supplemental Tables S2-S4, S7, S8, S11, S13-S17 [[supplements/187401_file04.pdf]](pending:yes)

## Data Availability

Data generated or analyzed during this study are available from the corresponding authors upon reasonable request.

## Supplemental Data

Supplemental Data include two figures and 17 tables.

## Declaration of interests

TO reports research support and honoraria from Abbott Diagnostics and Roche Diagnostics and honoraria from Siemens Healthineers, manufacturers of high sensitivity cardiac troponin assays. The other authors declare no competing interests.

## Web Resources

Minimac4, [https://genome.sph.umich.edu/wiki/Minimac4](https://genome.sph.umich.edu/wiki/Minimac4)

The GTEx Portal, [https://www.gtexportal.org](https://www.gtexportal.org)

The Neale Laboratory,[http://www.nealelab.is/ukbiobank](http://www.nealelab.is/ukbiobank)

The Cardiovascular Disease Knowledge Portal, [http://www.broadcvdi.org](http://www.broadcvdi.org)

## Data Availability

Data generated or analyzed during this study are available from the corresponding authors upon reasonable request.

## Author contributions

M.R.M analyzed the data and wrote the first draft of the manuscript. H.Røsjø contributed to the writing of the manuscript with knowledge on cardiac biology and cardiac troponin measurements. A.R, M.N.L, S.E.G, A.F.H, B.N.W, S.A.G.T, J.LF, H.Rasheed, L.F.T, W.Z, N.A., L.G.F, J.B.N and C.H contributed to analyses. B.M.B. and T.O. conceived and designed the study. All authors interpreted results and revised the paper.

## Acknowledgments

The Trøndelag Health Study (The HUNT Study) is a collaboration between HUNT Research Center (Faculty of Medicine and Health Sciences, NTNU, Norwegian University of Science and Technology), Trøndelag County Council, Central Norway Regional Health Authority, and the Norwegian Institute of Public Health. The genotyping in HUNT was financed by the National Institutes of Health; University of Michigan; the Research Council of Norway; the Liaison Committee for Education, Research and Innovation in Central Norway; and the Joint Research Committee between St Olav’s hospital and the Faculty of Medicine and Health Sciences, NTNU. The K.G. Jebsen Center for Genetic Epidemiology is funded by Stiftelsen Kristian Gerhard Jebsen; Faculty of Medicine and Health Sciences, NTNU; The Liaison Committee for education, research and innovation in Central Norway; and the Joint Research Committee between St. Olavs Hospital and the Faculty of Medicine and Health Sciences, NTNU. The cardiac troponin I analyses were sponsored by Abbott Diagnostics.

## Footnotes

*   Data on left ventricular mass added to sensitivity analyses; results and discussion revised according to additional analyses; Additional authors added

*   Received September 4, 2020.
*   Revision received December 4, 2020.
*   Accepted December 4, 2020.


*   © 2020, Posted by Cold Spring Harbor Laboratory

The copyright holder for this pre-print is the author. All rights reserved. The material may not be redistributed, re-used or adapted without the author's permission.

## References

1.  1.Thygesen K, Alpert JS, Jaffe AS, Chaitman BR, Bax JJ, Morrow DA, et al. Fourth Universal Definition of Myocardial Infarction (2018). Vol. 138, Circulation. 2018. 618–651 p.
    
    

2.  2.Farmakis D, Mueller C, Apple FS. High-sensitivity cardiac troponin assays for cardiovascular risk stratification in the general population. Eur Heart J. 2020;1–8.
    
    

3.  3.Sheng JJ, Jin JP. TNNI1, TNNI2 and TNNI3: Evolution, regulation, and protein structure-function relationships. Gene [Internet]. 2016;576(1):385–94. Available from: [http://dx.doi.org/10.1016/j.gene.2015.10.052](http://dx.doi.org/10.1016/j.gene.2015.10.052)
    
    

4.  4.Ingles J, Goldstein J, Thaxton C, Caleshu C, Corty EW, Crowley SB, et al. Evaluating the Clinical Validity of Hypertrophic Cardiomyopathy Genes. Circ Genomic Precis Med. 2019;12(2):57–64.
    
    

5.  5.Richard P, Charron P, Carrier L, Ledeuil C, Cheav T, Pichereau C, et al. Hypertrophic cardiomyopathy: Distribution of disease genes, spectrum of mutations, and implications for a molecular diagnosis strategy. Circulation. 2003;107(17):2227–32.
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6MTQ6ImNpcmN1bGF0aW9uYWhhIjtzOjU6InJlc2lkIjtzOjExOiIxMDcvMTcvMjIyNyI7czo0OiJhdG9tIjtzOjUwOiIvbWVkcnhpdi9lYXJseS8yMDIwLzEyLzA0LzIwMjAuMDkuMDQuMjAxODc0MDEuYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) 

6.  6.McCarthy CP, Raber I, Chapman AR, Sandoval Y, Apple FS, Mills NL, et al. Myocardial Injury in the Era of High-Sensitivity Cardiac Troponin Assays. JAMA Cardiol. 2019;4(10):1034–42.
    
    

7.  7.Seferović PM, Polovina M, Bauersachs J, Arad M,  Gal T Ben, Lund LH, et al. Heart failure in cardiomyopathies: a position paper from the Heart Failure Association of the European Society of Cardiology. Eur J Heart Fail. 2019;21(5):553–76.
    
    

8.  8.Lyngbakken MN, Røsjø H, Holmen OL, Nygård S, Dalen H, Hveem K, et al. Gender, High-Sensitivity Troponin I, and the Risk of Cardiovascular Events (from the Nord-Trøndelag Health Study). Am J Cardiol. 2016;118(6):816–21.
    
    

9.  9.Lazzeroni D, Rimoldi O, Camici PG. From left ventricular hypertrophy to dysfunction and failure. Circ J. 2016;80(3):555–64.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1253/circj.CJ-16-0062&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=26853555&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F12%2F04%2F2020.09.04.20187401.atom) 

10. 10.Göser S, Andrassy M, Buss SJ, Leuschner F, Volz CH, Öttl R, et al. Cardiac troponin I but not cardiac troponin T induces severe autoimmune inflammation in the myocardium. Circulation. 2006;114(16):1693–702.
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6MTQ6ImNpcmN1bGF0aW9uYWhhIjtzOjU6InJlc2lkIjtzOjExOiIxMTQvMTYvMTY5MyI7czo0OiJhdG9tIjtzOjUwOiIvbWVkcnhpdi9lYXJseS8yMDIwLzEyLzA0LzIwMjAuMDkuMDQuMjAxODc0MDEuYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) 

11. 11.Okazaki T, Tanaka Y, Nishio R, Mitsuiye T, Mizoguchi A, Wang J, et al. Autoantibodies against cardiac troponin I are responsible for dilated cardiomyopathy in PD-1-deficient mice. Nat Med. 2003;9(12):1477–83.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/nm955&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=14595408&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F12%2F04%2F2020.09.04.20187401.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000186845400030&link_type=ISI) 

12. 12.Hammarsten O, Mair J, Möckel M, Lindahl B, Jaffe AS. Possible mechanisms behind cardiac troponin elevations. Biomarkers [Internet]. 2018;23(8):725–34. Available from: [https://doi.org/10.1080/1354750X.2018.1490969](https://doi.org/10.1080/1354750X.2018.1490969)
    
    

13. 13.Welsh P, Preiss D, Hayward C, Shah ASV, McAllister D, Briggs A, et al. Cardiac Troponin T and Troponin I in the General Population. Circulation. 2019;139(24):2754–64.
    
    

14. 14.Smith GD, Ebrahim S. “Mendelian randomization”: Can genetic epidemiology contribute to understanding environmental determinants of disease? Int J Epidemiol. 2003;32(1):1–22.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/ije/dyg070&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=12689998&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F12%2F04%2F2020.09.04.20187401.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000182341300001&link_type=ISI) 

15. 15.Krokstad S, Langhammer A, Hveem K, Holmen TL, Midthjell K, Stene TR, et al. Cohort profile: The HUNT study, Norway. Int J Epidemiol. 2013;42(4):968–77.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/ije/dys095&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=22879362&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F12%2F04%2F2020.09.04.20187401.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000325167800009&link_type=ISI) 

16. 16.Taliun D, Harris DN, Kessler MD, Carlson J, Szpiech ZA, Torres R, et al. Sequencing of 53,831 diverse genomes from the NHLBI TOPMed Program. Biorxiv. 2019. p. 1–46.
    
    

17. 17.Smith BH, Campbell A, Linksted P, Fitzpatrick B, Jackson C, Kerr SM, et al. Cohort profile: Generation scotland: Scottish family health study (GS: SFHS). The study, its participants and their potential for genetic research on health and illness. Int J Epidemiol. 2013;42(3):689–700.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/ije/dys084&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=22786799&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F12%2F04%2F2020.09.04.20187401.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000322955900012&link_type=ISI) 

18. 18.Ferreira MA, Vonk JM, Baurecht H, Marenholz I, Tian C, Hoffman JD, et al. Shared genetic origin of asthma, hay fever and eczema elucidates allergic disease biology. Nat Genet. 2017;49(12):1752–7.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/ng.3985&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=29083406&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F12%2F04%2F2020.09.04.20187401.atom) 

19. 19.Nagy R, Boutin TS, Marten J, Huffman JE, Kerr SM, Campbell A, et al. Exploration of haplotype research consortium imputation for genome-wide association studies in 20,032 Generation Scotland participants. Genome Med. 2017;9(1):1.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1186/s13073-017-0403-7&link_type=DOI) 

20. 20.Apple FS, Collinson PO. Analytical characteristics of high-sensitivity cardiac troponin assays. Clin Chem. 2012;58(1):54–61.
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6ODoiY2xpbmNoZW0iO3M6NToicmVzaWQiO3M6NzoiNTgvMS81NCI7czo0OiJhdG9tIjtzOjUwOiIvbWVkcnhpdi9lYXJseS8yMDIwLzEyLzA0LzIwMjAuMDkuMDQuMjAxODc0MDEuYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) 

21. 21.Loh PR, Tucker G, Bulik-Sullivan BK, Vilhjálmsson BJ, Finucane HK, Salem RM, et al. Efficient Bayesian mixed-model analysis increases association power in large cohorts. Nat Genet. 2015;47(3):284–90.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/ng.3190&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=25642633&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F12%2F04%2F2020.09.04.20187401.atom) 

22. 22.Willer CJ, Li Y, Abecasis GR. METAL: Fast and efficient meta-analysis of genomewide association scans. Bioinformatics. 2010;26(17):2190–1.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/bioinformatics/btq340&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=20616382&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F12%2F04%2F2020.09.04.20187401.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000281738900017&link_type=ISI) 

23. 23.Yang J, Lee SH, Goddard ME, Visscher PM. GCTA: A tool for genome-wide complex trait analysis. Am J Hum Genet [Internet]. 2011;88(1):76–82. Available from: [http://dx.doi.org/10.1016/j.ajhg.2010.11.011](http://dx.doi.org/10.1016/j.ajhg.2010.11.011)
    
    

24. 24.Yang J, Benyamin B, McEvoy BP, Gordon S, Henders AK, Nyholt DR., et al. Common SNPs explain a large proportion of heritability for human height. Nat Genet. 2010;42(7):565–9.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/ng.608&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=20562875&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F12%2F04%2F2020.09.04.20187401.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000279242400007&link_type=ISI) 

25. 25.Wang K, Li M, Hakonarson H. ANNOVAR: Functional annotation of genetic variants from high-throughput sequencing data. Nucleic Acids Res. 2010;38(16):1– 7.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/nar/gkp829&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=19843612&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F12%2F04%2F2020.09.04.20187401.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000273229100007&link_type=ISI) 

26. 26.Wang G, Sarkar AK, Carbonetto P, Stephens M. A simple new approach to variable selection in regression, with application to genetic fine-mapping [Internet]. bioRxiv. 2018. Available from: [http://biorxiv.org/content/early/2018/12/19/501114.abstract](http://biorxiv.org/content/early/2018/12/19/501114.abstract)
    
    

27. 27.Pruim RJ, Welch RP, Sanna S, Teslovich TM, Chines PS, Gliedt TP, et al. LocusZoom: Regional visualization of genome-wide association scan results. Bioinformatics. 2011;27(13):2336–7.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/bioinformatics/ btq419&link_type=DOI) 

28. 28.Pers TH, Karjalainen JM, Chan Y, Westra HJ, Wood AR, Yang J, et al. Biological interpretation of genome-wide association studies using predicted gene functions. Nat Commun. 2015;6:5890.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/ncomms6890&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=25597830&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F12%2F04%2F2020.09.04.20187401.atom) 

29. 29.Chang CC, Chow CC, Tellier LCAM, Vattikuti S, Purcell SM, Lee JJ. Second- generation PLINK: Rising to the challenge of larger and richer datasets. Gigascience. 2015;4(1):1–16.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1186/2047-217X-4-1&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=25838885&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F12%2F04%2F2020.09.04.20187401.atom) 

30. 30.Giambartolomei C, Vukcevic D, Schadt EE, Franke L, Hingorani AD, Wallace C, et al. Bayesian Test for Colocalisation between Pairs of Genetic Association Studies Using Summary Statistics. PLoS Genet. 2014;10(5).
    
    

31. 31.R Core Team. R: A Language and Environment for Statistical Computing [Internet]. Vienna, Austria; 2019. Available from: [https://www.r-project.org/](https://www.r-project.org/)
    
    

32. 32.Manichaikul A, Mychaleckyj JC, Rich SS, Daly K, Sale M, Chen WM. Robust relationship inference in genome-wide association studies. Bioinformatics. 2010;26(22):2867–73.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/bioinformatics/btq559&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=20926424&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F12%2F04%2F2020.09.04.20187401.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000283919800010&link_type=ISI) 

33. 33.Zhou W, Nielsen JB, Fritsche LG, Dey R, Gabrielsen ME, Wolford BN, et al. Efficiently controlling for case-control imbalance and sample relatedness in large- scale genetic association studies. Nat Genet [Internet]. 2018;50(9):1335–41. Available from: [http://dx.doi.org/10.1038/s41588-018-0184-y](http://dx.doi.org/10.1038/s41588-018-0184-y)
    
    

34. 34.Hemani G, Haycock P, Zheng J, Gaunt T, Elsworth B. TwoSampleMR: Two Sample MR functions and interface to MR Base database. 2020.
    
    

35. 35.Hemani G. MRInstruments: Data sources for genetic instruments to be used in MR [Internet]. 2020. Available from: [https://github.com/MRCIEU/MRInstruments](https://github.com/MRCIEU/MRInstruments)
    
    

36. 36.Sudlow C, Gallacher J, Allen N, Beral V, Burton P, Danesh J, et al. UK Biobank: An Open Access Resource for Identifying the Causes of a Wide Range of Complex Diseases of Middle and Old Age. PLoS Med. 2015;12(3):1–10.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1371/journal.pmed.1001878&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=25710373&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F12%2F04%2F2020.09.04.20187401.atom) 

37. 37.Shah S, Henry A, Roselli C, Lin H, Sveinbjörnsson G, Fatemifar G, et al. Genome- wide association and Mendelian randomisation analysis provide insights into the pathogenesis of heart failure. Nat Commun. 2020;11(1):1–12.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41467-019-13889-6&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=31911652&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F12%2F04%2F2020.09.04.20187401.atom) 

38. 38.Aung N, Vargas JD, Yang C, Cabrera CP, Warren HR, Fung K, et al. Genome-wide analysis of left ventricular image-derived phenotypes identifies fourteen loci associated with cardiac morphogenesis and heart failure development. Circulation. 2019;140(16):1318–30.
    
    

39. 39.Bossard M, Thériault S, Aeschbacher S, Schoen T, Kunz S, von Rotz M, et al. Factors independently associated with cardiac troponin I levels in young and healthy adults from the general population. Clin Res Cardiol. 2017;106(2):96–104.
    
    

40. 40.Lyngbakken MN, Myhre PL, Røsjø H, Omland T. Novel biomarkers of cardiovascular disease: Applications in clinical practice. Crit Rev Clin Lab Sci [Internet]. 2019;56(1):33–60. Available from: [https://doi.org/10.1080/10408363.2018.1525335](https://doi.org/10.1080/10408363.2018.1525335)
    
    

41. 41.Willeit P, Welsh P, Evans JDW, Tschiderer L, Boachie C, Jukema JW, et al. High- Sensitivity Cardiac Troponin Concentration and Risk of First-Ever Cardiovascular Outcomes in 154,052 Participants. J Am Coll Cardiol. 2017;70(5):558–68.
    
    [FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6MzoiUERGIjtzOjExOiJqb3VybmFsQ29kZSI7czo0OiJhY2NqIjtzOjU6InJlc2lkIjtzOjg6IjcwLzUvNTU4IjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjAvMTIvMDQvMjAyMC4wOS4wNC4yMDE4NzQwMS5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 

42. 42.de Lemos JA, Drazner MH, Omland T, Ayers CR, Khera A, Rohatgi A, et al. Association of Troponin T Detected With a Highly Sensitive Assay and Cardiac Structure and Mortality Risk in the General Population. JAMA. 2010;304(22).
    
    

43. 43.Omland T, De Lemos JA, Sabatine MS, Christophi CA, Rice MM, Jablonski KA, et al. A sensitive cardiac troponin T assay in stable coronary artery disease. N Engl J Med. 2009;361(26):2538–47.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1056/NEJMoa0805299&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=19940289&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F12%2F04%2F2020.09.04.20187401.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000273003700009&link_type=ISI) 

44. 44.Omland T, Pfeffer MA, Solomon SD, De Lemos JA, Røsjø H, Benth JŠ, et al. Prognostic value of cardiac troponin I measured with a highly sensitive assay in patients with stable coronary artery disease. J Am Coll Cardiol. 2013;61(12):1240– 9.
    
    [FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6MzoiUERGIjtzOjExOiJqb3VybmFsQ29kZSI7czo0OiJhY2NqIjtzOjU6InJlc2lkIjtzOjEwOiI2MS8xMi8xMjQwIjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjAvMTIvMDQvMjAyMC4wOS4wNC4yMDE4NzQwMS5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 

45. 45.Chin CWL, Shah ASV, McAllister DA, Cowell SJ, Alam S, Langrish JP, et al. High-sensitivity troponin I concentrations are a marker of an advanced hypertrophic response and adverse outcomes in patients with aortic stenosis. Eur Heart J. 2014;35(34):2312–21.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/eurheartj/ehu189&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=24829362&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F12%2F04%2F2020.09.04.20187401.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000342233800014&link_type=ISI) 

46. 46.Røsjø H, Hijazi Z, Omland T, Westerbergh J, Lyngbakken MN, Alexander JH, et al. Cardiac troponin is associated with cardiac outcomes in men and women with atrial fibrillation, insights from the ARISTOTLE trial. J Intern Med. 2020;
    
    

47. 47.Olsson M, Stanne TM, Pedersen A, Lorentzen E, Kara E, Martinez-Palacian A, et al. Genome-wide analysis of genetic determinants of circulating factor VII-activating protease (FSAP) activity. J Thromb Haemost. 2018;16(10):2024–34.
    
    

48. 48.Willeit J, Kiechl S, Weimer T, Mair A, Santer P, Wiedermann CJ, et al. Marburg I polymorphism of factor VII-activating protease: A prominent risk predictor of carotid stenosis. Circulation. 2003;107(5):667–70.
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6MTQ6ImNpcmN1bGF0aW9uYWhhIjtzOjU6InJlc2lkIjtzOjk6IjEwNy81LzY2NyI7czo0OiJhdG9tIjtzOjUwOiIvbWVkcnhpdi9lYXJseS8yMDIwLzEyLzA0LzIwMjAuMDkuMDQuMjAxODc0MDEuYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) 

49. 49.Ahmad-Nejad P, Dempfle CE, Weiss C, Bugert P, Borggrefe M, Neumaier M. The G534E-polymorphism of the gene encoding the Factor VII-activating protease is a risk factor for venous thrombosis and recurrent events. Thromb Res. 2012;130(3):441–4.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.thromres.2012.02.009&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=22421107&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F12%2F04%2F2020.09.04.20187401.atom) 

50. 50.Kanai M, Akiyama M, Takahashi A, Matoba N, Momozawa Y, Ikeda M, et al. Genetic analysis of quantitative traits in the Japanese population links cell types to complex human diseases. Nat Genet [Internet]. 2018;50(3):390–400. Available from: [http://dx.doi.org/10.1038/s41588-018-0047-6](http://dx.doi.org/10.1038/s41588-018-0047-6)
    
    

51. 51.Hicks D, Sarkozy A, Muelas N, Köehler K, Huebner A, Chinnery PF, et al. Europe PMC Funders Group Europe PMC Funders Author Manuscripts A founder mutation in Anoctamin 5 is a major cause of limb girdle muscular dystrophy. 2014;134(0 1):171–82.
    
    

52. 52.Matsuyama S, Kage Y, Fujimoto N, Ushijima T, Tsuruda T, Kitamura K, et al. Interaction between cardiac myosin-binding protein C and formin Fhod3. Proc Natl Acad Sci U S A. 2018;115(19):E4386–95.
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NDoicG5hcyI7czo1OiJyZXNpZCI7czoxMjoiMTE1LzE5L0U0Mzg2IjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjAvMTIvMDQvMjAyMC4wOS4wNC4yMDE4NzQwMS5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 

53. 53.Esslinger U, Garnier S, Korniat A, Proust C, Kararigas G, Müller-Nurasyid M, et al. Exome-wide association study reveals novel susceptibility genes to sporadic dilated cardiomyopathy. PLoS One. 2017;12(3).
    
    

54. 54.Pirruccello JP, Bick A, Wang M, Chaffin M, Friedman S, Yao J, et al. Analysis of cardiac magnetic resonance imaging in 36,000 individuals yields genetic insights into dilated cardiomyopathy. Nat Commun [Internet]. 2020;11(1):2254. Available from: [http://www.ncbi.nlm.nih.gov/pubmed/32382064](http://www.ncbi.nlm.nih.gov/pubmed/32382064)
    
    

55. 55.Phelan DG, Anderson DJ, Howden SE, Wong RCB, Hickey PF, Pope K, et al. ALPK3-deficient cardiomyocytes generated from patient-derived induced pluripotent stem cells and mutant human embryonic stem cells display abnormal calcium handling and establish that ALPK3 deficiency underlies familial cardiomyopathy. Eur Heart J. 2016;37(33):2586–90.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/eurheartj/ehw160&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=27106955&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F12%2F04%2F2020.09.04.20187401.atom) 

56. 56.Norland K, Sveinbjornsson G, Thorolfsdottir RB, Davidsson OB, Tragante V, Rajamani S, et al. Sequence variants with large effects on cardiac electrophysiology and disease. Nat Commun. 2019;10(1):1–10.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41467-019-09078-0&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=30602773&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F12%2F04%2F2020.09.04.20187401.atom) 

57. 57.Shiraishi S, Zhou C, Aoki T, Sato N, Chiba T, Tanaka K, et al. TBP-interacting Protein 120B (TIP120B)/Cullin-associated and Neddylation-dissociated 2 (CAND2) inhibits SCF-dependent ubiquitination of myogenin and accelerates myogenic differentiation. J Biol Chem. 2007;282(12):9017–28.
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6MzoiamJjIjtzOjU6InJlc2lkIjtzOjExOiIyODIvMTIvOTAxNyI7czo0OiJhdG9tIjtzOjUwOiIvbWVkcnhpdi9lYXJseS8yMDIwLzEyLzA0LzIwMjAuMDkuMDQuMjAxODc0MDEuYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) 

58. 58.Jia L, Shi Y, Wen Y, Li W, Feng J, Chen C. The roles of TNFAIP2 in cancers and infectious diseases. J Cell Mol Med. 2018;22(11):5188–95.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1111/jcmm.13822&link_type=DOI) 

59. 59.Nielsen JB, Thorolfsdottir RB, Fritsche LG, Zhou W, Skov MW, Graham SE, et al. Biobank-driven genomic discovery yields new insight into atrial fibrillation biology. Nat Genet [Internet]. 2018;50(9):1234–9. Available from: [http://dx.doi.org/10.1038/s41588-018-0171-3](http://dx.doi.org/10.1038/s41588-018-0171-3)
    
    

60. 60.Inouye M, Ripatti S, Kettunen J, Lyytikäinen LP, Oksala N, Laurila PP, et al. Novel Loci for Metabolic Networks and Multi-Tissue Expression Studies Reveal Genes for Atherosclerosis. PLoS Genet. 2012;8(8):e1002907.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1371/journal.pgen.1002907&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=22916037&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F12%2F04%2F2020.09.04.20187401.atom) 

61. 61.Björkqvist J, Jämsä A,Renné T. Plasma kallikrein: The bradykinin-producing enzyme. Thromb Haemost. 2013;110(3):399–407.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1160/TH13-03-0258&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=23846131&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F12%2F04%2F2020.09.04.20187401.atom) 

62. 62.Nichols WM, Seligsohn U, Zivelin A, Terry VH, Hertel CE, Wheatley MA, et al. Mutations in the ER-Golgi intermediate compartment protein ERGIC-53 cause combined deficiency of coagulation factors V and VIII. Cell. 1998;93(1):61–70.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/S0092-8674(00)81146-0&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=9546392&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F12%2F04%2F2020.09.04.20187401.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000072910800010&link_type=ISI) 

63. 63.Xiong Y, Wang L, Mo P, Huang G, Li A, Chai R, et al. Association between HLA- DQB1 alleles and susceptibility to coronary artery disease in Southern Han Chinese. Hum Immunol. 2017;78(9):540–6.
    
    

64. 64.Köttgen A, Pattaro C, Böger CA, Fuchsberger C, Olden M, Glazer NL, et al. Multiple New Loci Associated with Kidney Function and ChronicKidney Disease: The CKDGen consortium. Nat Genet. 2010;42(5):376–84.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/ng.568&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=20383146&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F12%2F04%2F2020.09.04.20187401.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000277179500006&link_type=ISI) 

65. 65.Houlihan LM, Davies G, Tenesa A, Harris SE, Luciano M, Gow AJ, et al. Common Variants of Large Effect in F12, KNG1, and HRG Are Associated with Activated Partial Thromboplastin Time. Am J Hum Genet [Internet]. 2010;86(4):626–31. Available from: [http://dx.doi.org/10.1016/j.ajhg.2010.02.016](http://dx.doi.org/10.1016/j.ajhg.2010.02.016)
    
    

66. 66.Verweij N, Mahmud H, Leach IM, De Boer RA, Brouwers FP, Yu H, et al. Genome-wide association study on plasma levels of midregional- proadrenomedullin and C-terminal-pro-endothelin-1. Hypertension. 2013;61(3):602–8.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1161/HYPERTENSIONAHA.111.203117&link_type=DOI) 

67. 67.Williams TM. Human leukocyte antigen gene polymorphism and the histocompatibility laboratory. J Mol Diagnostics. 2001;3(3):98–104.
    
    

68. 68.Omland T. New features of troponin testing in different clinical settings. J Intern Med. 2010 Sep;268(3):207–17.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1111/j.1365-2796.2010.02253.x&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=20626552&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F12%2F04%2F2020.09.04.20187401.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000280643200001&link_type=ISI) 

69. 69.Manisastry SM, Zaal KJM, Horowits R. Myofibril assembly visualized by imaging N-RAP, alpha-actinin, and actin in living cardiomyocytes. Exp Cell Res [Internet]. 2009;315(12):2126–39. Available from: [http://dx.doi.org/10.1016/j.yexcr.2009.02.006](http://dx.doi.org/10.1016/j.yexcr.2009.02.006)
    
    

70. 70.Gehmlich K, Geier C, Osterziel KJ, Van Der Ven PFM, Fürst DO. Decreased interactions of mutant muscle LIM protein (MLP) with N-RAP and α-actinin and their implication for hypertrophic cardiomyopathy. Cell Tissue Res. 2004;317(2):129–36.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1007/s00441-004-0873-y&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=15205937&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F12%2F04%2F2020.09.04.20187401.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000222865200003&link_type=ISI) 

71. 71.Brentnall M, Rodriguez-Menocal L, De Guevara RL, Cepero E, Boise LH. Caspase- 9, caspase-3 and caspase-7 have distinct roles during intrinsic apoptosis. BMC Cell Biol. 2013;14(1).
    
    

72. 72.van Eldik W, Passier R. Signalling in sarcomeres in development and disease. Netherlands Hear J. 2013;21(–8):367–71.
    
    

73. 73.Fassett J, Xu X, Kwak D, Zhu G, Fassett EK, Zhang P, et al. Adenosine kinase attenuates cardiomyocyte microtubule stabilization and protects against pressure overload-induced hypertrophy and LV dysfunction. J Mol Cell Cardiol [Internet]. 2019;130(November 2018):49–58. Available from: [https://doi.org/10.1016/j.yjmcc.2019.03.015](https://doi.org/10.1016/j.yjmcc.2019.03.015)
    
    

74. 74.Huxley R, Mendis S, Zheleznyakov E, Reddy S, Chan J. Body mass index, waist circumference and waist:hip ratio as predictors of cardiovascular riska review of the literature. Eur J Clin Nutr. 2010;64(1):16–22.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/ejcn.2009.68&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=19654593&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F12%2F04%2F2020.09.04.20187401.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000273563700004&link_type=ISI) 

75. 75.Evans JM, Wang S, Greb C, Kostas V, Knapp CF, Zhang Q, et al. Body size predicts cardiac and vascular resistance effects on men’s and women’s blood pressure. Front Physiol. 2017;8(AUG):1–12.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.3389/fphys.2017.00788&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=28154536&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F12%2F04%2F2020.09.04.20187401.atom) 

76. 76.Panteghini M. Aspartate aminotransferase isoenzymes. Clin Biochem. 1990;23(4):311–9.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/0009-9120(90)80062-N&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=2225456&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F12%2F04%2F2020.09.04.20187401.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=A1990DX25400009&link_type=ISI) 

77. 77.Pratt DS, Kaplan MM. Evaluation of Abnormal Liver-Enzyme Results In Asumptomatic Patients. N Engl J Med. 2000;342(17):1266–71.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1056/NEJM200004273421707&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=10781624&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F12%2F04%2F2020.09.04.20187401.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000086660600007&link_type=ISI) 

78. 78.Black PH, Garbutt LD. Stress, inflammation and cardiovascular disease. J Psychosom Res. 2002;52(1):1–23.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/S0022-3999(01)00302-6&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=11801260&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F12%2F04%2F2020.09.04.20187401.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000173548800001&link_type=ISI)