#PAGE_PARAMS# #ADS_HEAD_SCRIPTS# #MICRODATA#

Genetics of height and risk of atrial fibrillation: A Mendelian randomization study


Authors: Michael G. Levin aff001;  Renae Judy aff004;  Dipender Gill aff005;  Marijana Vujkovic aff002;  Shefali S. Verma aff010;  Yuki Bradford aff010aff012;  Marylyn D. Ritchie aff010;  Matthew C. Hyman aff001;  Saman Nazarian aff001;  Daniel J. Rader aff002;  Benjamin F. Voight aff010;  Scott M. Damrauer aff003
Authors place of work: Division of Cardiovascular Medicine, University of Pennsylvania Perelman School of Medicine, Philadelphia, Pennsylvania, United States of America aff001;  Department of Medicine, University of Pennsylvania Perelman School of Medicine, Philadelphia, Pennsylvania, United States of America aff002;  Corporal Michael J. Crescenz VA Medical Center, Philadelphia, Pennsylvania, United States of America aff003;  Department of Surgery, University of Pennsylvania Perelman School of Medicine, Philadelphia, Pennsylvania, United States of America aff004;  Department of Epidemiology and Biostatistics, School of Public Health, Imperial College London, London, United Kingdom aff005;  Centre for Pharmacology & Therapeutics, Department of Medicine, Imperial College London, London, United Kingdom aff006;  Novo Nordisk Research Centre Oxford, Oxford, United Kingdom aff007;  Clinical Pharmacology and Therapeutics Section, Institute of Medical and Biomedical Education and Institute for Infection and Immunity, St George’s, University of London, London, United Kingdom aff008;  Clinical Pharmacology Group, Pharmacy and Medicines Directorate, St George’s University Hospitals NHS Foundation Trust, London, United Kingdom aff009;  Department of Genetics, University of Pennsylvania Perelman School of Medicine, Philadelphia, Pennsylvania, United States of America aff010;  Institute for Biomedical Informatics, University of Pennsylvania Perelman School of Medicine, Philadelphia, Pennsylvania, United States of America aff011;  Tarrytown, New York, United States of America aff012;  Institute for Translational Medicine and Therapeutics, University of Pennsylvania Perelman School of Medicine, Philadelphia, Pennsylvania, United States of America aff013;  Department of Systems Pharmacology and Translational Therapeutics, University of Pennsylvania Perelman School of Medicine, Philadelphia, Pennsylvania, United States of America aff014
Published in the journal: Genetics of height and risk of atrial fibrillation: A Mendelian randomization study. PLoS Med 17(10): e32767. doi:10.1371/journal.pmed.1003288
Category: Research Article
doi: https://doi.org/10.1371/journal.pmed.1003288

Summary

Background

Observational studies have identified height as a strong risk factor for atrial fibrillation, but this finding may be limited by residual confounding. We aimed to examine genetic variation in height within the Mendelian randomization (MR) framework to determine whether height has a causal effect on risk of atrial fibrillation.

Methods and findings

In summary-level analyses, MR was performed using summary statistics from genome-wide association studies of height (GIANT/UK Biobank; 693,529 individuals) and atrial fibrillation (AFGen; 65,446 cases and 522,744 controls), finding that each 1-SD increase in genetically predicted height increased the odds of atrial fibrillation (odds ratio [OR] 1.34; 95% CI 1.29 to 1.40; p = 5 × 10−42). This result remained consistent in sensitivity analyses with MR methods that make different assumptions about the presence of pleiotropy, and when accounting for the effects of traditional cardiovascular risk factors on atrial fibrillation. Individual-level phenome-wide association studies of height and a height genetic risk score were performed among 6,567 European-ancestry participants of the Penn Medicine Biobank (median age at enrollment 63 years, interquartile range 55–72; 38% female; recruitment 2008–2015), confirming prior observational associations between height and atrial fibrillation. Individual-level MR confirmed that each 1-SD increase in height increased the odds of atrial fibrillation, including adjustment for clinical and echocardiographic confounders (OR 1.89; 95% CI 1.50 to 2.40; p = 0.007). The main limitations of this study include potential bias from pleiotropic effects of genetic variants, and lack of generalizability of individual-level findings to non-European populations.

Conclusions

In this study, we observed evidence that height is likely a positive causal risk factor for atrial fibrillation. Further study is needed to determine whether risk prediction tools including height or anthropometric risk factors can be used to improve screening and primary prevention of atrial fibrillation, and whether biological pathways involved in height may offer new targets for treatment of atrial fibrillation.

Keywords:

Genetics – Genome-wide association studies – coronary heart disease – Genetics of disease – atrial fibrillation – Instrumental variable analysis – diabetes mellitus – Cardiovascular disease risk

Introduction

Atrial fibrillation is a common cardiac arrhythmia, with a population prevalence of 0.5%, affecting more than 33 million individuals worldwide [1]. A number of risk factors are associated with atrial fibrillation, including chronic diseases like chronic kidney disease, heart failure, thyroid disease, obesity, obstructive sleep apnea, sleep apnea, and valvular heart disease, as well as cardiac surgery, smoking, and anthropometric factors [25]. Recent studies have identified both common and rare genetic variants at more than 100 independent loci contributing to the incidence of atrial fibrillation, and heritability is estimated at 20% [68]. Even with treatment, affected individuals are at risk of cardioembolic stroke, heart failure, and death [2].

Height has been identified as a risk factor for a number of cardiometabolic diseases, including coronary artery disease, atrial fibrillation, and venous thromboembolism [9,10]. The relationship between height and atrial fibrillation in particular has been identified in large observational studies, with greater height strongly associated with increased risk of atrial fibrillation [4,1118]. These studies are limited in assessing the causality and potential mediators or confounders of this association by their observational design, and randomized controlled trials of anthropometric traits are not feasible. Although a number of factors influence height including childhood illness/development, diet/nutrition, and socioeconomic status, height has a strong genetic component. Large genome-wide association studies (GWASs) have provided heritability estimates of 60%–70%, and have identified more than 700 independent loci that contribute to height [1922].

In the current study, we utilize human genetic data within the Mendelian randomization (MR) framework to evaluate a potential causal association between height and atrial fibrillation. MR analysis exploits the random assortment of genetic variants during meiosis as an instrumental variable to estimate the causal relationship between a trait and an outcome of interest. Here, we use summary data from large, multiethnic GWASs of height (693,529 individuals) and atrial fibrillation (65,446 cases and 522,744 controls) to estimate the effect of genetically predicted height on risk of atrial fibrillation [7,19]. We then use participant-level data from the Penn Medicine Biobank within the observational phenome-wide association study (PheWAS) framework to identify other clinical associations with height, and within the single-sample MR framework to further assess the impact of height on atrial fibrillation after adjustment for clinical risk factors.

Methods

This study is reported per the Strengthening the Reporting of Observational Studies in Epidemiology (STROBE) guideline (S1 STROBE Checklist). The study did not have a pre-registered or published analysis plan. All analyses were planned prior to study initiation.

Two-sample MR analysis

Summary-level data for height was obtained from a 2018 meta-analysis of GWASs of height [19]. This analysis was a fixed-effects meta-analysis combining results from a GWAS of height performed among 456,426 participants from the UK Biobank (adjusted for age, sex, recruitment center, genotyping batch, and 10 genetic principal components) and results from a 2014 GWAS published by the GIANT (Genetic Investigation of ANthropometric Traits) Consortium, which included 253,288 participants from 79 studies (adjusted for age, sex, and study-specific covariates). Genetic variants associated with height at genome-wide significance (p < 5 × 10−8) were then LD-pruned (distance threshold = 10,000 kb, r2 = 0.001) using the clump_data command in the TwoSampleMR package in R to identify an independent set of variants to serve as a genetic instrument for height [23]. The independent variants associated with height at genome-wide significance (p < 5 × 10−8) were then harmonized with variants from the 2018 Roselli et al. [7] atrial fibrillation GWAS from the AFGen (Atrial Fibrillation Genetics) Consortium, using the default settings of the harmonize_data command in the TwoSampleMR package in R to ensure that effect estimates were aligned to the same allele. This study included 65,446 atrial fibrillation cases and 522,744 controls from more than 50 studies (84.2% European, 12.5% Japanese, 2% African American, and 1.3% Brazilian and Hispanic), including participants from UK Biobank, Biobank Japan, other international biobanks, and international cardiovascular cohort studies (adjusted for age, sex, and study-specific covariates) [7]. In total, 707 independent SNPs associated with height were available in the atrial fibrillation outcome GWAS. These 707 SNPs accounted for 11.2% of the variability in height (S1 Table).

Inverse-variance-weighted 2-sample MR was used as the primary analysis, with weighted median, MR-Egger, and MR-PRESSO performed as sensitivity analyses to account for potential violations of the instrumental variable assumptions [24]. Further sensitivity analysis was performed using a genetic instrument for height constructed excluding any SNPs nominally (p < 0.05) associated with coronary artery disease, high-density lipoprotein (HDL), low-density lipoprotein (LDL), total cholesterol, triglycerides, fasting glucose, fasting insulin, diabetes, BMI, waist-to-hip ratio, and systolic blood pressure, all identified using the MR-Base database [2528].

For each variant included in the genetic instruments, variance (R2) was calculated using the formula R2 = 2 × MAF × (1−MAF) × beta2 (where MAF represents the effect allele frequency and beta represents the effect estimate of the genetic variant in the exposure GWAS). F-statistics were then calculated for each variant using the formula F=R2×(N-2)1-R2 (where R2 represents the variance in exposure explained by the genetic variant, and N represents the number of individuals in the exposure GWAS) to assess the strength of the selected instruments [29]. Bias due to sample overlap between the exposure and outcome GWASs was estimated and found to be negligible (S1 Methods) [30]. The estimates for the effect of height on atrial fibrillation are reported per 1–standard deviation (SD) increase in height. A scaled estimate per 10-cm increase in height has also been calculated for the main analysis based on the February 2020 UK Biobank population standard deviation of height (9.28 cm).

To account for the possibility that the restrictive genetic instrument may introduce collider or ascertainment bias by conditioning on associated traits, multivariable MR was performed using the TwoSampleMR package in R. This method allows the direct effects of multiple traits on an outcome to be determined jointly. The effect of height on atrial fibrillation was estimated in analyses that individually adjusted for potential confounders, including coronary artery disease, HDL, LDL, total cholesterol, triglycerides, fasting glucose, fasting insulin, diabetes, BMI, waist-to-hip ratio, and systolic blood pressure, using genetic variants obtained from MR-Base as above. Potential confounders that had significant (p < 0.05) direct effects on atrial fibrillation in the individual multivariable models were then combined in a final model to jointly estimate their direct effects on atrial fibrillation.

Bidirectional MR was performed to assess for the possibility of reverse causality. A genetic instrument for atrial fibrillation containing 73 independent, genome-wide significant variants was constructed using the same methods as for the height instrument. Inverse-variance-weighted, weighted median, and MR-Egger analyses were performed.

Height genetic risk score

A standardized, weighted genetic risk score (GRS) for height was calculated for each individual in the Penn Medicine Biobank using imputed dosage information in an additive model, weighted by the effect size of each independent (distance threshold = 10,000 kb, r2 = 0.001), genome-wide significant (p < 5 × 10−8) variant from the 2018 Yengo et al. GWAS of height [19], using the clump_data command in the TwoSampleMR package in R. Scores were standardized, centered to mean = 0, and scaled.

PheWAS

PheWASs were performed to identify clinical diagnoses associated both measured height and the height GRS, using the default setting of the PheWAS package in R [31]. Briefly, international classification disease codes obtained from the electronic health record were mapped to “Phecodes,” and individuals were assigned case/control status or excluded using default mapping parameters. Association between height/height GRS and phenotypes was tested using logistic regression, adjusted for age, age2, sex, and 10 genetic principal components. The 10 genetic principal components included in the PheWAS analysis were selected a priori, without investigation of the variation of the data. For the more detailed individual-level MR analyses, 6 genetic principal components were identified qualitatively using a scree plot as explaining the majority of variation in the data (S4 Fig). The Bonferroni method was used to account for multiple testing (p = 0.05/1,816 phenotypes).

Individual-level MR

To estimate the population-averaged effect of increased height on risk of atrial fibrillation, individual-level MR was performed using the 2-stage method [32]. This analysis included 6,567 individuals with available height, genotype, and clinical covariate information. The instrumental variable was a GRS for height, computed from independent, genome-wide significant variants, weighted by effect on height in the 2018 GIANT/UK Biobank GWAS meta-analysis. Of the 707 genetic variants included in the main instrumental variable for height in the 2-sample MR analysis, 695 were available for use in the individual-level analysis (S3 Table). In the first stage of the 2-stage process, a linear regression was fitted with height as the dependent variable, and the GRS for height as the independent variable, among the subset of individuals without atrial fibrillation, adjusted for age, sex, and 6 genetic principal components. The 6 genetic principal components were selected qualitatively using a scree plot as explaining the majority of variation in the data (S4 Fig). In the second stage, a logistic regression model with robust standard errors was fit for both atrial fibrillation cases and controls, incorporating the fitted values of height from the first stage, adjusted for age, sex, and the 6 genetic principal components, with atrial fibrillation as the outcome. In additional models, each stage was adjusted for clinical diagnoses of hypertension, coronary artery disease, heart failure, hyperlipidemia, diabetes, chronic kidney disease, stroke, valvular heart disease, sleep apnea, thyroid disease, cardiac surgery, smoking, and echocardiographic left atrial size to account for potential cofounders/mediators of the height–atrial fibrillation relationship. These risk factors were selected for being common cardiometabolic risk factors or having prior associations with atrial fibrillation [2]. Sensitivity analysis was performed including only individuals who had available echocardiographic data.

MR assumptions

MR relies on the naturally random assortment of genetic variants at conception to provide unbiased estimates of the effect of an exposure (in this case height) on an outcome (atrial fibrillation). To provide unbiased estimates, MR depends on 3 key assumptions [33]. First, genetic instruments must be associated with the exposure of interest (relevance assumption). Weak instruments that are not strongly associated with the exposure of interest may result in biased estimates. To estimate the strength of the height genetic instrument, the F-statistic was calculated (S1 Table). The second assumption is that no confounders of the SNP–outcome association are present. This assumption is not readily testable because some confounders may not be measured, although at both the 2-sample and individual-level the MR analyses were adjusted for common cardiovascular risk factors. Further, because genetic variants are randomly allocated at conception, confounders of the SNP–outcome association should be distributed across the population. The third assumption is that the effect of the genetic instrument on atrial fibrillation is entirely through the effect on height (exclusion restriction). We performed a PheWAS using the GRS for height to assess for potentially pleiotropic effects of height. Further, we applied multiple MR methods including weighted median, MR-Egger, and MR-PRESSO, which are more robust to the presence of horizontal pleiotropy, a potential violation of the exclusion restriction assumption [24]. Finally, MR assumes that genetically mediated effects are similar to exogenous effects. In this case, height is highly heritable, although the effect of changes in stature related to environmental exposures (e.g., fetal growth restriction, malnutrition) may differ from the genetic effects assessed here.

Penn Medicine Biobank

The Penn Medicine Biobank is a longitudinal genomics and precision medicine study in which participants consent to linkage of genomic information and biospecimens to the electronic health record. More than 60,000 individuals are currently enrolled. This study included 6,567 individuals of European ancestry (genetically determined) who underwent genotyping and had available electronic health record data. For individual-level analyses in the Penn Medicine Biobank, continuous demographic variables were summarized by mean and standard deviation, with categorical variables summarized by count and group percent. The Wilcoxon rank sum test and Fisher’s exact test, respectively, were used to assess for significant differences between groups.

Phenotype ascertainment

For individual-level analyses in Penn Medicine Biobank, phenotypes were determined by querying the electronic health record. International Classification of Diseases (ICD)–9/10 and Current Procedural Terminology (CPT) codes, in addition to laboratory measurements and vital signs, were used to identify height, weight, BMI, smoking status, diagnoses of heart failure, hypertension, diabetes mellitus, chronic kidney disease, sleep apnea, stroke, thyroid disease, valvular heart disease, and cardiac surgery. We identified 6,567 participants with available high-quality genotype and electronic health record phenotype data, with transthoracic echocardiogram data available for a subset of 2,842 individuals. Atrial fibrillation was defined using ICD-9/10 codes 427.31, I48.0, I48.1, I48.2, and I48.91. Atrial fibrillation ascertainment in the 2018 Roselli et al. GWAS meta-analysis was study-specific [7]. Data were extracted as of January 2017 (S1 Methods).

Ethical approval

The investigators from the Penn Medicine Biobank (MGL, RJ, MV, SSV, YB, MDR, DJR, BFV, SMD) have received approval from the University of Pennsylvania Institutional Review Board. No additional approval was required for analyses of publicly available summary statistics.

Statistical analysis

All statistical analyses were performed using R version 3.5.1 [34]. For statistical analyses, p < 0.05 was considered statistically significant.

Results

Population-level MR

To characterize the relationship between increasing height and risk of atrial fibrillation, we constructed a genetic instrument for height using 707 independent SNPs associated with height at a genome-wide level of significance (p < 5 × 10−8), which accounted for 11.2% of the variability in height (S1 Table). The mean F-statistic was 110 (range 22–948), suggesting the risk of weak instrument bias was low [30]. We performed 2-sample MR using summary statistics from a GWAS of atrial fibrillation including 65,446 atrial fibrillation cases and 522,744 controls. Inverse-variance-weighted modeling identified a significant association between increasing height and atrial fibrillation (odds ratio [OR] 1.34 per 1-SD increase in height; 95% CI 1.29 to 1.40; p = 5 × 10−42) (Fig 1). This corresponds to an OR of 1.37 (95% CI 1.31 to 1.44) per 10-cm increase in height among UK Biobank participants. The intercept from Egger regression was −0.001 (p = 0.13), thus not providing evidence for significant pleiotropic bias. Similar estimates were obtained in sensitivity analyses from weighted median, MR-Egger, and MR-PRESSO models.

Two-sample Mendelian randomization (MR).
Fig. 1. Two-sample Mendelian randomization (MR).
Two-sample MR was performed using a genetic instrument containing 707 independent SNPs associated with height. (A) Each point represents the SNP effects on height and atrial fibrillation. Colored lines represent inverse-variance-weighted (red), weighted median (green), and MR-Egger (blue) estimates of the association between a 1-SD increase in height and risk of atrial fibrillation. (B) Odds ratios (ORs), 95% confidence intervals (CIs), and p-values for MR estimates.

An additional genetic instrument for height was constructed, excluding variants nominally associated (p < 0.05) with potentially pleiotropic risk factors for atrial fibrillation. The selection of SNPs using a nominal p-value association of 0.05 was an arbitrarily chosen threshold to more liberally account for potentially pleiotropic effects of height-associated SNPs. The restrictive instrument consisted of 224 independent SNPs, in sum explaining 2.8% of the variance in height (S2 Table). The mean F-statistic was 88 (range 22–867). MR results using this restrictive genetic instrument were similar (S1 Fig).

We next performed multivariable MR. Height remained significantly associated with atrial fibrillation after adjustment for the effect of genetic variants separately on each of coronary artery disease, HDL, LDL, total cholesterol, triglycerides, fasting glucose, fasting insulin, diabetes, BMI, waist-to-hip ratio, and systolic blood pressure (S4 Table). This analysis identified significant associations of BMI, systolic blood pressure, total cholesterol, and coronary artery disease with atrial fibrillation after adjustment for height. The effect of height on risk of atrial fibrillation was similar in models accounting for these risk factors individually, and in a combined model jointly considering genetic variants associated with height, BMI, systolic blood pressure, total cholesterol, and coronary artery disease (Fig 2 and S5 Table).

Multivariable Mendelian randomization analysis.
Fig. 2. Multivariable Mendelian randomization analysis.
Multivariable Mendelian randomization was performed to jointly consider the effect of genetic variants for cardiometabolic traits and height on atrial fibrillation. The effect of a 1-SD increase in height on risk of atrial fibrillation estimated by each model is displayed. ORs per 1-SD increase in height, 95% confidence intervals (CIs), and p-values for Mendelian randomization estimates are displayed.

To assess for the possibility of reverse causation, where genetic variants primarily associated with increased risk of atrial fibrillation may represent invalid instruments for height, a bidirectional MR analysis was performed. Using a genetic instrument including 73 independent, genome-wide significant variants associated with atrial fibrillation, there was no evidence for a causal association with height (S2 Fig).

Demographics

Individual-level analysis focused on 6,567 European-ancestry individuals from the Penn Medicine Biobank with available genotype data linked to the electronic health record, recruited between 2008 and 2015. Individuals with atrial fibrillation were older (mean 66 versus 60 years; p < 0.001) than individuals without atrial fibrillation. Individuals with atrial fibrillation were taller (mean 174 versus 170 cm; p < 0.001), and significantly more likely to have been diagnosed with chronic kidney disease, coronary artery disease, heart failure, hyperlipidemia, hypertension, prior cardiac surgery, sleep apnea, smoking, stroke, or valvular heart disease, compared to individuals without atrial fibrillation (Table 1 and S6S8). Among the subset of individuals with echocardiograms, individuals with atrial fibrillation had larger left atrial diameter (mean 4.39 versus 3.92 cm; p < 0.001).

Tab. 1. Demographics.
Demographics.

PheWASs of height and height GRS

To determine the association between genetically predicted height and clinical diagnoses, we first constructed a GRS for height using weights derived from the 2018 GIANT/UK Biobank height GWAS meta-analysis. PheWASs were performed to identify clinical diagnoses associated with both measured height and the height GRS across the Penn Medicine Biobank. Each 1-SD increase in measured height was associated with increased risk of atrial fibrillation (OR 1.55; 95% CI 1.41 to 1.71; p = 3.6 × 10−18) and decreased risk of coronary atherosclerosis (OR 0.66; p = 6.9 × 10−18). Each 1-SD increase in height GRS was associated with increased risk of atrial fibrillation and flutter (OR 1.16; 95% CI 1.09 to 1.23; p = 7 × 10−6) (Fig 3).

Phenome-wide associations of clinical diagnoses with height and height genetic risk score (GRS).
Fig. 3. Phenome-wide associations of clinical diagnoses with height and height genetic risk score (GRS).
Phenome-wide association studies were performed in Penn Medicine Biobank participants to identify clinical phenotypes associated with increased (A) height and (B) height GRS. The horizontal red line denotes the Bonferroni-adjusted level of significance (0.05/1,816 phenotypes), the blue line denotes the nominal level of significance (0.05), and triangles denote the direction of association between increasing height and risk of the phenotype (pointing upward = increased risk; pointing downward = decreased risk).

Individual-level MR

Individual-level MR was performed in Penn Medicine Biobank participants to further assess the association between height and atrial fibrillation. Using the height GRS as an instrumental variable, height was significantly associated with atrial fibrillation (OR 1.75 per 1-SD increase in height; 95% CI 1.53 to 2.0; p = 6 × 10−5) after adjusting for age, sex, and 6 genetic principal components (Fig 4). This corresponds to an OR of 1.66 (95% CI 1.46 to 1.86) per 10-cm increase in height among Penn Medicine Biobank participants. Results were similar in sex-stratified analysis. Height remained associated with atrial fibrillation after adjustment for weight, hypertension, coronary artery disease, heart failure, hyperlipidemia, diabetes, chronic kidney disease, sleep apnea, stroke, thyroid disease, smoking, cardiac surgery, and valvular heart disease. After further adjustment for left atrial size, height remained significantly associated with atrial fibrillation (OR 1.83; 95% CI 1.44 to 2.31; p = 0.01). When individual-level MR was restricted to the subset of individuals with complete data, the effect estimates were similar, with modest attenuation with sequential adjustment for clinical risk factors and left atrial size (S3 Fig).

Individual-level instrumental variable analysis.
Fig. 4. Individual-level instrumental variable analysis.
Individual-level instrumental variable analysis was performed in Penn Medicine Biobank participants, using a weighted genetic risk score for height as an instrumental variable for measured height. The base model was adjusted for age, sex, and 6 genetic principal components. Model 1 was additionally adjusted for weight, hypertension, coronary artery disease, heart failure, hyperlipidemia, diabetes, chronic kidney disease, sleep apnea, stroke, thyroid disease, smoking, cardiac surgery, and valvular heart disease. Model 2 was additionally adjusted for left atrial size as measured on transthoracic echocardiogram. Odds ratios (ORs) are reported per 1-SD increase in height.

Discussion

In this study we used both population- and individual-level genetic information to test the association between height and atrial fibrillation. At the population level, there was a strong causal association between genetic determinants of height and risk of atrial fibrillation. This finding was robust to multiple sensitivity analyses of the MR methods and the genetic instrument for height. Observational analysis at the individual level identified a strong association between height and increased risk of atrial fibrillation, and decreased risk of coronary artery disease. MR analysis at both the population and individual levels suggested that height remains a causal risk factor for atrial fibrillation even after adjustment for other traditional risk factors.

Our findings are consistent with prior observational analyses that have identified height as a risk factor for atrial fibrillation [4,1118]. These studies, including a large Swedish national cohort study of 1.5 million military conscripts recruited over the course of 28 years, have consistently identified a strong association between height and atrial fibrillation [4]. Observational designs have limited these studies, due to the possibility of residual confounding. The Helsinki Birth Cohort Study partially addressed this limitation by considering the effect of maternal height on risk of atrial fibrillation in offspring, but was limited by a small, homogenous sample [35]. We further build on those prior findings using the MR framework, considering both summary-level and individual-level genetic data, examining genetic variants associated with height and atrial fibrillation to provide a causal estimate for the effect of increasing height on risk of atrial fibrillation. Our population-level MR analysis reinforces a prior MR estimate for association between increasing height and atrial fibrillation (OR 1.33; 95% CI 1.26 to 1.40) using updated GWAS summary statistics for height and atrial fibrillation [9]. By including a large multiethnic study of atrial fibrillation, we increase the generalizability of the prior MR findings, which were limited only to white participants of the UK Biobank.

Several mechanisms have been proposed to explain the relationship between height and atrial fibrillation. Increased left atrial size is correlated with height, and has been identified as an independent predictor of atrial fibrillation [17,18,36,37]. Consistent with findings from the Cardiovascular Health Study, however, we found that the effect of height on risk of atrial fibrillation was not attenuated after adjustment for left atrial diameter [14]. It is possible that other, more nuanced markers of left atrial structure and function that have been associated with severity of atrial fibrillation, such as left atrial volume index, emptying fraction, expansion index, and contractile function, may also be affected by height and may better explain the association between height and atrial fibrillation [38]. Similarly, 2-dimensional echocardiography significantly underestimates left atrial size compared to MRI assessment [39]. Investigation of these factors in the 2-sample setting is limited by the lack of genetic studies of cardiac structure/function by echocardiography, MRI, or cardiac CT, and study of these parameters in large cohorts with genetics and imaging data is warranted. Bioimpedance and dual-energy X-ray absorptiometry measures of body composition have also been associated with atrial fibrillation [40]. While the current study focused primarily on common cardiometabolic risk factors used in clinical practice, the possibility remains that more advanced anthropometric screening beyond height, weight, body mass index, and waist-to-hip ratio may explain some of the effect of increased height on atrial fibrillation.

MR has previously identified an association between shorter stature and increased risk of coronary artery disease, mediated in part by increases in LDL cholesterol and triglyceride levels [10]. Our PheWAS of measured height similarly identified an association between height and decreased risk of coronary artery disease, and our PheWAS of a height GRS identified a nominally (p < 0.05) protective association between height and coronary artery disease. In the current MR analysis, we were unable to detect significant associations between these lipid traits and risk of atrial fibrillation when considered in multivariable models alongside height. Similarly, in multivariable models considering height alongside common cardiometabolic risk factors for atrial fibrillation, effect estimates were all similar, suggesting these factors may not substantially mediate the effect of height on risk of atrial fibrillation.

Although we detected no evidence of horizontal pleiotropy and our results remained robust to extensive sensitivity analyses, height is a highly polygenic trait, and the possibility remains that pleiotropic pathways mediate the association between height and atrial fibrillation. A recent GWAS found genes at atrial-fibrillation-associated loci to be enriched in pathways important for tissue formation [41]. Coupled with the finding that genes at height-associated loci are enriched for cardiovascular and endocrine tissue types, these results raise the possibility of a more complex shared genetic architecture affecting both height and atrial fibrillation [22]. Thus, we cannot exclude the possibility that genetic variants broadly associated with growth and development may simultaneously affect height and establish structural cardiovascular changes that may predispose to atrial fibrillation.

The findings of this study have several clinical implications. Anthropometric characteristics have been included in clinical models that predict incident atrial fibrillation, including BMI in a model derived from the Framingham Heart Study, and height in the CHARGE-AF (Cohorts for Heart and Aging Research in Genomic Epidemiology–Atrial Fibrillation) model from the CHARGE-AF Consortium [4244]. While height is not a readily modifiable risk factor, the recognition that taller individuals have increased risk of atrial fibrillation may prompt more aggressive management of modifiable cardiovascular risk factors like overweight/obesity, prediabetes/diabetes, hypertension, and alcohol/tobacco use. As both height and atrial fibrillation are heritable, it is possible that height might act as a surrogate for family history of atrial fibrillation when this information is not readily available. Further study may clarify the risks of incident atrial fibrillation attributable to height and family history. Similarly, further mechanistic investigation of the height–atrial fibrillation link may identify novel modifiable pathways for intervention. Future study is needed to determine whether risk prediction tools including height or other anthropometric factors can be used to improve screening and primary prevention of atrial fibrillation.

Our study has several limitations. First, observational analyses from the Penn Medicine Biobank may not be generalizable, as this population represents a cohort within a single academic health system. Second, despite our use of multiethnic GWASs of height and atrial fibrillation, the underlying studies focused primarily on individuals of European ancestry. Genetic studies in broader populations are warranted to further improve the generalizability of these findings. Third, our population-level MR analyses relied on publicly available summary statistics that contain some overlapping samples/cohorts. Two-sample MR tends to bias the causal effect estimates to the null, but sample overlap may make the estimate susceptible to weak instrument bias. In this study, however, despite overlapping cohorts, simulation suggests that the large sample sizes of the height and atrial fibrillation GWASs, and large F-statistics, make the risk of weak instrument bias low [30].

In conclusion, we find that increased height is associated with increased risk of atrial fibrillation, and this relationship is likely to be causal. These results raise the possibility of investigating height/growth-related pathways as a means for gaining novel mechanistic insights into atrial fibrillation, as well as the possibility of incorporating height into larger targeted screening strategies for atrial fibrillation.

Supporting information

S1 Fig [a]
Two-sample MR sensitivity analysis.

S2 Fig [a]
Two-sample MR directionality sensitivity analysis.

S3 Fig [ors]
Individual MR sensitivity analysis.

S4 Fig [tif]
Genetic principal component scree plot for Penn Medicine Biobank.

S1 Methods [docx]
Supplemental methods.

S1 STROBE Checklist [docx]
STROBE checklist.

S1 Table [xlsx]
Full genetic instrument for height used in 2-sample MR analysis including independent (distance threshold = 10,000 kb, = 0.001), genome-wide significant ( < 5 × 10) variants.

S2 Table [xlsx]
Restrictive genetic instrument for height used in 2-sample MR analysis including independent (distance threshold = 10,000 kb, = 0.001), genome-wide significant ( < 5 × 10) variants, excluding those SNPs nominally ( < 0.05) associated with coronary artery disease, HDL, LDL, total cholesterol, triglycerides, fasting glucose, fasting insulin, diabetes, BMI, waist-to-hip ratio, and systolic blood pressure in the MR-Base database.

S3 Table [xlsx]
Full genetic instrument for height used in individual-level MR analysis including independent (distance threshold = 10,000 kb, = 0.001), genome-wide significant ( < 5 × 10) variants.

S4 Table [xlsx]
Results of bivariate multivariable MR analysis considering the effect of height and each of coronary artery disease, HDL, LDL, total cholesterol, triglycerides, fasting glucose, fasting insulin, diabetes, BMI, waist-to-hip ratio, and systolic blood pressure on atrial fibrillation.

S5 Table [xlsx]
Results of combined multivariable MR analysis considering the effects of height, body mass index, total cholesterol, systolic blood pressure, and coronary heart disease.

S6 Table [xlsx]
Summary demographics for Penn Medicine Biobank participants overall.

S7 Table [xlsx]
Summary demographics for Penn Medicine Biobank participants stratified by quartile of height GRS.

S8 Table [xlsx]
Summary demographics for Penn Medicine Biobank participants limited to subset of individuals with complete phenotype and echocardiographic data.


Zdroje

1. Chugh SS, Roth GA, Gillum RF, Mensah GA. Global burden of atrial fibrillation in developed and developing nations. Global Heart. 2014;9:113–9. doi: 10.1016/j.gheart.2014.01.004 25432121

2. January CT, Wann LS, Alpert JS, Calkins H, Cigarroa JE, Cleveland JC, et al. 2014 AHA/ACC/HRS guideline for the management of patients with atrial fibrillation. J Am Coll Cardiol. 2014;64:e1–76. doi: 10.1016/j.jacc.2014.03.022 24685669

3. Tikkanen E, Gustafsson S, Ingelsson E. Associations of fitness, physical activity, strength, and genetic risk with cardiovascular disease: longitudinal analyses in the UK Biobank study. Circulation. 2018;137:2583–91. doi: 10.1161/CIRCULATIONAHA.117.032432 29632216

4. Crump C, Sundquist J, Winkleby MA, Sundquist K. Height, weight, and aerobic fitness level in relation to the risk of atrial fibrillation. Am J Epidemiol. 2018;187:417–26. doi: 10.1093/aje/kwx255 28641376

5. Fenger-Grøn M, Overvad K, Tjønneland A, Frost L. Lean body mass is the predominant anthropometric risk factor for atrial fibrillation. J Am Coll Cardiol. 2017;69:2488–97. doi: 10.1016/j.jacc.2017.03.558 28521886

6. Choi SH, Weng L-C, Roselli C, Lin H, Haggerty CM, Shoemaker MB, et al. Association between titin loss-of-function variants and early-onset atrial fibrillation. JAMA. 2018;320:2354–64. doi: 10.1001/jama.2018.18179 30535219

7. Roselli C, Chaffin MD, Weng L-C, Aeschbacher S, Ahlberg G, Albert CM, et al. Multi-ethnic genome-wide association study for atrial fibrillation. Nat Genet. 2018;50:1225–33. doi: 10.1038/s41588-018-0133-9 29892015

8. 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. 2018;50:1234–9. doi: 10.1038/s41588-018-0171-3 30061737

9. Lai FY, Nath M, Hamby SE, Thompson JR, Nelson CP, Samani NJ. Adult height and risk of 50 diseases: a combined epidemiological and genetic analysis. BMC Med. 2018;16:187. doi: 10.1186/s12916-018-1175-7 30355295

10. Nelson CP, Hamby SE, Saleheen D, Hopewell JC, Zeng L, Assimes TL, et al. Genetically determined height and coronary artery disease. N Engl J Med. 2015;372:1608–18. doi: 10.1056/NEJMoa1404881 25853659

11. Marott JL, Skielboe AK, Dixen U, Friberg JB, Schnohr P, Jensen GB. Increasing population height and risk of incident atrial fibrillation: the Copenhagen City Heart Study. Eur Heart J. 2018;39:4012–9. doi: 10.1093/eurheartj/ehy367 29961878

12. Siscovick DS, Ix JH, Gottdiener JS, Kizer JR, Mukamal KJ, Yee LM, et al. Measures of body size and composition and risk of incident atrial fibrillation in older people: the Cardiovascular Health Study. Am J Epidemiol. 2016;183:998–1007. doi: 10.1093/aje/kwv278 27188936

13. Schmidt M, Bøtker HE, Pedersen L, Sørensen HT. Adult height and risk of ischemic heart disease, atrial fibrillation, stroke, venous thromboembolism, and premature death: a population based 36-year follow-up study. Eur J Epidemiol. 2014;29:111–8. doi: 10.1007/s10654-013-9867-y 24337942

14. Chang JD, Zimetbaum PJ, Rosenberg MA, Patton KK, Sotoodehnia N, Karas MG, et al. The impact of height on the risk of atrial fibrillation: the Cardiovascular Health Study. Eur Heart J. 2012;33:2709–17. doi: 10.1093/eurheartj/ehs301 22977225

15. Suzuki S, Yamashita T, Ohtsuka T, Sagara K, Uejima T, Oikawa Y, et al. Body size and atrial fibrillation in Japanese outpatients. Circ J. 2010;74:66–70. doi: 10.1253/circj.cj-09-0431 19797821

16. Rosengren A, Hauptman PJ, Lappas G, Olsson L, Wilhelmsen L, Swedberg K. Big men and atrial fibrillation: effects of body size and weight gain on risk of atrial fibrillation in men. Eur Heart J. 2009;30:1113–20. doi: 10.1093/eurheartj/ehp076 19304990

17. Mont L, Tamborero D, Elosua R, Molina I, Coll-Vinent B, Sitges M, et al. Physical activity, height, and left atrial size are independent risk factors for lone atrial fibrillation in middle-aged healthy individuals. Europace. 2008;10:15–20. doi: 10.1093/europace/eum263 18178694

18. Hanna IR, Heeke B, Bush H, Brosius L, King-Hageman D, Beshai JF, et al. The relationship between stature and the prevalence of atrial fibrillation in patients with left ventricular dysfunction. J Am Coll Cardiol. 2006;47:1683–8. doi: 10.1016/j.jacc.2005.11.068 16631009

19. Yengo L, Sidorenko J, Kemper KE, Zheng Z, Wood AR, Weedon MN, et al. Meta-analysis of genome-wide association studies for height and body mass index in ∼700000 individuals of European ancestry. Hum Mol Genet. 2018;27:3641–9. doi: 10.1093/hmg/ddy271 30124842

20. Yang J, Bakshi A, Zhu Z, Hemani G, Vinkhuyzen AAE, Lee SH, et al. Genetic variance estimation with imputed variants finds negligible missing heritability for human height and body mass index. Nat Genet. 2015;47:1114–20. doi: 10.1038/ng.3390 26323059

21. Yang J, Zeng J, Goddard ME, Wray NR, Visscher PM. Concepts, estimation and interpretation of SNP-based heritability. Nat Genet. 2017;49:1304–10. doi: 10.1038/ng.3941 28854176

22. Wood AR, Esko T, Yang J, Vedantam S, Pers TH, Gustafsson S, et al. Defining the role of common variation in the genomic and biological architecture of adult human height. Nat Genet. 2014;46:1173–86. doi: 10.1038/ng.3097 25282103

23. Hemani G, Zheng J, Elsworth B, Wade KH, Haberland V, Baird D, et al. The MR-Base platform supports systematic causal inference across the human phenome. Elife. 2018;7:e34408. doi: 10.7554/eLife.34408 29846171

24. Verbanck M, Chen C-Y, Neale B, Do R. Detection of widespread horizontal pleiotropy in causal relationships inferred from Mendelian randomization between complex traits and diseases. Nat Genet. 2018;50:693–8. doi: 10.1038/s41588-018-0099-7 29686387

25. Willer CJ, Schmidt EM, Sengupta S, Peloso GM, Gustafsson S, Kanoni S, et al. Discovery and refinement of loci associated with lipid levels. Nat Genet. 2013;45:1274–85. doi: 10.1038/ng.2797 24097068

26. Shungin D, Winkler T, Croteau-Chonka DC, Ferreira T, Locke AE, Mägi R, et al. New genetic loci link adipose and insulin biology to body fat distribution. Nature. 2015;518:187–96. doi: 10.1038/nature14132 25673412

27. DIAbetes Genetics Replication And Meta-analysis (DIAGRAM) Consortium, Asian Genetic Epidemiology Network Type 2 Diabetes (AGEN-T2D) Consortium, South Asian Type 2 Diabetes (SAT2D) Consortium, Mexican American Type 2 Diabetes (MAT2D) Consortium, Type 2 Diabetes Genetic Exploration by Nex-generation sequencing in muylti-Ethnic Samples (T2D-GENES) Consortium, Mahajan A, et al. Genome-wide trans-ancestry meta-analysis provides insight into the genetic architecture of type 2 diabetes susceptibility. Nat Genet. 2014;46:234–44. doi: 10.1038/ng.2897 24509480

28. Dupuis J, Langenberg C, Prokopenko I, Saxena R, Soranzo N, Jackson AU, et al. New genetic loci implicated in fasting glucose homeostasis and their impact on type 2 diabetes risk. Nat Genet. 2010;42:105–16. doi: 10.1038/ng.520 20081858

29. Palmer TM, Lawlor DA, Harbord RM, Sheehan NA, Tobias JH, Timpson NJ, et al. Using multiple genetic variants as instrumental variables for modifiable risk factors. Stat Methods Med Res. 2012;21:223–42. doi: 10.1177/0962280210394459 21216802

30. Burgess S, Davies NM, Thompson SG. Bias due to participant overlap in two-sample Mendelian randomization. Genet Epidemiol. 2016;40:597–608. doi: 10.1002/gepi.21998 27625185

31. Carroll RJ, Bastarache L, Denny JC. R PheWAS: data analysis and plotting tools for phenome-wide association studies in the R environment. Bioinformatics. 2014;30:2375–6. doi: 10.1093/bioinformatics/btu197 24733291

32. Burgess S, Small DS, Thompson SG. A review of instrumental variable estimators for Mendelian randomization. Stat Methods Med Res. 2017;26:2333–55. doi: 10.1177/0962280215597579 26282889

33. Davies NM, Holmes M V, Davey Smith G. Reading Mendelian randomisation studies: a guide, glossary, and checklist for clinicians. BMJ. 2018;362:k601. doi: 10.1136/bmj.k601 30002074

34. R Development Core Team. R: a language and environment for statistical computing. Version 3.5.1. Vienna: R Foundation for Statistical Computing; 2017.

35. Johnson LSB, Salonen M, Kajantie E, Conen D, Healey JS, Osmond C, et al. Early life risk factors for incident atrial fibrillation in the Helsinki Birth Cohort Study. J Am Heart Assoc. 2017;6:e006036. doi: 10.1161/JAHA.117.006036 28649086

36. Pritchett AM, Jacobsen SJ, Mahoney DW, Rodeheffer RJ, Bailey KR, Redfield MM. Left atrial volume as an index of left atrial size: a population-based study. J Am Coll Cardiol. 2003;41:1036–43. doi: 10.1016/s0735-1097(02)02981-9 12651054

37. Leung DY, Chi C, Allman C, Boyd A, Ng AC, Kadappu KK, et al. Prognostic implications of left atrial volume index in patients in sinus rhythm. Am J Cardiol. 2010;105:1635–9. doi: 10.1016/j.amjcard.2010.01.027 20494675

38. Gupta DK, Shah AM, Giugliano RP, Ruff CT, Antman EM, Grip LT, et al. Left atrial structure and function in atrial fibrillation: ENGAGE AF-TIMI 48. Eur Heart J. 2014;35:1457–65. doi: 10.1093/eurheartj/eht500 24302269

39. Rodevand O, Bjornerheim R, Ljosland M, Maehle J, Smith HJ, Ihlen H. Left atrial volumes assessed by three- and two-dimensional echocardiography compared to MRI estimates. Int J Card Imaging. 1999;15:397–410. doi: 10.1023/a:1006276513186 10595406

40. Tikkanen E, Gustafsson S, Knowles JW, Perez M, Burgess S, Ingelsson E. Body composition and atrial fibrillation: a Mendelian randomization study. Eur Heart J. 2019;40:1277–82. doi: 10.1093/eurheartj/ehz003 30721963

41. Nielsen JB, Fritsche LG, Zhou W, Teslovich TM, Holmen OL, Gustafsson S, et al. Genome-wide study of atrial fibrillation identifies seven risk loci and highlights biological pathways and regulatory elements involved in cardiac development. Am J Hum Genet. 2018;102:103–15. doi: 10.1016/j.ajhg.2017.12.003 29290336

42. Alonso A, Krijthe BP, Aspelund T, Stepas KA, Pencina MJ, Moser CB, et al. Simple risk model predicts incidence of atrial fibrillation in a racially and geographically diverse population: the CHARGE-AF consortium. J Am Heart Assoc. 2013;2:e000102. doi: 10.1161/JAHA.112.000102 23537808

43. Alonso A, Roetker NS, Soliman EZ, Chen LY, Greenland P, Heckbert SR. Prediction of atrial fibrillation in a racially diverse cohort: the Multi-Ethnic Study of Atherosclerosis (MESA). J Am Heart Assoc. 2016;5:e003077. doi: 10.1161/JAHA.115.003077 26908413

44. Schnabel RB, Sullivan LM, Levy D, Pencina MJ, Massaro JM, D’Agostino RB, et al. Development of a risk score for atrial fibrillation (Framingham Heart Study): a community-based cohort study. Lancet. 2009;373:739–45. doi: 10.1016/S0140-6736(09)60443-8 19249635


Článek vyšel v časopise

PLOS Medicine


2020 Číslo 10
Nejčtenější tento týden
Nejčtenější v tomto čísle
Kurzy Podcasty Doporučená témata Časopisy
Přihlášení
Zapomenuté heslo

Zadejte e-mailovou adresu, se kterou jste vytvářel(a) účet, budou Vám na ni zaslány informace k nastavení nového hesla.

Přihlášení

Nemáte účet?  Registrujte se

#ADS_BOTTOM_SCRIPTS#