- Open Access
Estimating causal effects of atherogenic lipid-related traits on COVID-19 susceptibility and severity using a two-sample Mendelian randomization approach
BMC Medical Genomics volume 14, Article number: 269 (2021)
As the number of COVID-19 deaths continues to rise worldwide, the identification of risk factors for the disease is an urgent issue, and it remains controversial whether atherogenic lipid-related traits including serum apolipoprotein B, low-density lipoprotein (LDL)-cholesterol, and triglyceride levels, are risk factors. The aim of this study was to estimate causal effects of lipid-related traits on COVID-19 risk in the European population using a two-sample Mendelian randomization (MR) approach.
We used summary statistics from a genome-wide association study (GWAS) that included 441,016 participants from the UK Biobank as the exposure dataset of lipid-related traits and from COVID-19 Host Genetics Initiative GWAS meta-analyses of European ancestry as the outcome dataset for COVID-19 susceptibility (32,494 cases and 1,316,207 controls), hospitalization (8316 cases and 1,549,095 controls), and severity (4792 cases and 1,054,664 controls). We performed two-sample MR analyses using the inverse variance weighted (IVW) method. As sensitivity analyses, the MR-Egger regression, weighted median, and weighted mode methods were conducted as were leave-one-out sensitivity analysis, the MR-PRESSO global test, PhenoScanner searches, and IVW multivariable MR analyses. A P value below 0.0055 with Bonferroni correction was considered statistically significant.
This MR study suggested that serum apolipoprotein B or LDL-cholesterol levels were not significantly associated with COVID-19 risk. On the other hand, we inferred that higher serum triglyceride levels were suggestively associated with higher risks of COVID-19 susceptibility (odds ratio [OR] per standard deviation increase in lifelong triglyceride levels, 1.065; 95% confidence interval [CI], 1.001–1.13; P = 0.045) and hospitalization (OR, 1.174; 95% CI, 1.04–1.33; P = 0.012), and were significantly associated with COVID-19 severity (OR, 1.274; 95% CI, 1.08–1.50; P = 0.004). Sensitivity and bidirectional MR analyses suggested that horizontal pleiotropy and reverse causation were unlikely.
Our MR study indicates a causal effect of higher serum triglyceride levels on a greater risk of COVID-19 severity in the European population using the latest and largest GWAS datasets to date. However, as the underlying mechanisms remain unclear and our study might be still biased due to possible horizontal pleiotropy, further studies are warranted to validate our findings and investigate underlying mechanisms.
The World Health Organization has reported that the number of deaths from coronavirus disease 2019 (COVID-19) caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) continues to rise worldwide (over 4.6 million as of September 2021) . Therefore, in addition to the establishment of effective therapies, the identification of risk factors for COVID-19 is an urgent issue. Observational studies have reported that the severity of COVID-19 depends on risk factors such as obesity, coronary artery disease (CAD), and diabetes [2, 3]. Although dyslipidemia is associated with these risk factors [4,5,6], different results from observational studies have been reported with regard to an association between dyslipidemia and COVID-19 risk [7,8,9,10,11]. Moreover, observational studies tend to suffer from bias due to possible confounders and reverse causation .
The Mendelian randomization (MR) method mimics a study with a randomized controlled design using single-nucleotide variants (SNVs) (also called single-nucleotide polymorphisms [SNPs]) as instrumental variables (IVs) and can estimate causal effects of risk factors on diseases of interest. According to Mendel’s law, genetic variants are randomly assigned at meiosis. Therefore, MR studies are less likely to suffer from possible confounders or reverse causation, which as stated above are limitations of observational studies . A recent genome-wide association study (GWAS) using UK Biobank (UKBB) data and MR analysis  reported that among atherogenic lipid-related traits (apolipoprotein B [Apo-B], low-density lipoprotein cholesterol [LDL-C], and triglyceride [TG]), Apo-B accounted for the causal effect on CAD risk, independently of LDL-C and TGs. Nevertheless, the estimated effects of atherogenic lipid-related traits on COVID-19 risk are inconsistent even among MR studies [14,15,16,17]. The aim of the present study was to estimate causal effects of serum Apo-B, LDL-C, and TG levels on risk of COVID-19 susceptibility, hospitalization, and severity in the European population using a two-sample MR approach.
We performed two-sample univariable MR analyses using summary-level GWAS datasets to estimate causal effects of circulating atherogenic lipid-related traits on COVID-19 using genetically predicted serum Apo-B, LDL-C, and TG levels as exposures and risk of COVID-19 susceptibility, hospitalization, and severity as outcomes. To examine reverse causation, we also performed bidirectional two-sample univariable MR analyses using genetically predicted risk of COVID-19 susceptibility, hospitalization, and severity as exposures and serum Apo-B, LDL-C, and TG levels as outcomes.
All analyses were conducted using the TwoSampleMR package (version 0.5.6) in R software (version 4.0.3) .
A P value below 0.0055 (0.005/3/3 by Bonferroni correction) was considered statistically significant and a P value between 0.0055 and 0.05 was considered suggestively significant in the MR analyses.
For the exposure dataset of genetically predicted serum Apo-B, LDL-C, and TG levels, summary statistics were available from a GWAS  that included up to 441,016 participants from UKBB. To our knowledge, this is to date the largest GWAS in sample size of atherogenic lipid-related traits. The data were standardized and normalized such that the mean was 0 and the standard deviation (SD) was 1. However, the mean and SD of TG were not provided in the GWAS; instead, the median TG was given as 1.50 (IQR = 1.11) mmol/L. The mean (SD) TG of 470,434 participants in UKBB was 1.75 (1.02) mmol/L . The GWAS datasets are publicly available from the IEU open GWAS database  as GWAS-IDs of “ieu-b-108” for serum Apo-B level, “ieu-b-110” for serum LDL-C level, and “ieu-b-111” for serum TG level.
For the outcome dataset of genetically predicted risk of COVID-19 susceptibility, hospitalization, and severity, summary statistics were available from the GWAS meta-analyses by COVID-19 Host Genetics Initiative (COVID-19-HGI)  (Round 5) of European ancestry, which excluded UKBB data. Therefore, in our two-sample MR study, possible study overlap between the exposure and outcome datasets was unlikely. The outcome dataset included 1,348,701 participants (32,494 laboratory confirmed cases of SARS-CoV-2 infection and 1,316,207 population controls) for COVID-19 susceptibility, 1,557,411 participants (8316 hospitalized COVID-19 patients and 1,549,095 population controls) for COVID-19 hospitalization, and 1,059,456 participants (4792 very severe respiratory confirmed COVID-19 cases and 1,054,664 controls) for COVID-19 severity. COVID-19-HGI defined very severe respiratory confirmed COVID-19 cases as patients hospitalized for laboratory-confirmed SARS-CoV-2 infection who died or were given respiratory support. The dataset (Round 5) was the latest version with the largest sample size involving European ancestry released on January 18, 2021 .
For the dataset of genetically predicted body mass index (BMI) trait, summary statistics are publicly available from a GWAS of European ancestry (a meta-analysis of GIANT [Genetic Investigation of ANthropometric Traits] consortium studies and UKBB with a total of 681,275 participants) , and from the IEU open GWAS database  as GWAS-ID of “ieu-b-40”.
Selection of instrumental variables
MR analyses use SNVs as IVs that must satisfy the following three assumptions : the IVs are associated with the exposure; the IVs affect the outcome only via exposure; and the IVs are not associated with confounders.
To estimate causal effects, we selected SNVs from the exposure GWAS dataset as IVs by clumping together all SNVs associated with the trait with P < 5.0 × 10–8 (a genome-wide significance level) and not in linkage disequilibrium with other SNVs (r2 < 0.001, and distance > 10,000 kb) using the clump_data function (population = "EUR"). We extracted the summary statistics for each SNV from both the exposure and outcome GWAS datasets and then harmonized them. We did not include proxy SNVs in the analysis [25, 26]. We excluded palindromic SNVs with an intermediate minor allele frequency (MAF) > 0.42 [24, 27]. To evaluate the strength of the IVs, the F-statistic for each SNV was calculated using the following formula: F-statistic = R2 × (N − 2)/(1 − R2), where R2 is the proportion of variance in phenotype explained by each SNV in exposure, and N is the sample size. We calculated R2 using the following formula: R2 = 2 × (Beta/SD)2 × (MAF) × (1 − MAF), where Beta is the per allele effect size of the association between each SNV and phenotype . IVs with an F-statistic below 10 (if any existed) were considered weak instruments .
Two-sample Mendelian randomization and sensitivity analyses
The Wald ratio, which estimates the causal effect of each IV exposure on outcome, was calculated as the ratio of Beta for the corresponding SNV in the outcome dataset divided by Beta for the same SNV in the exposure dataset . We conducted a meta-analysis of each Wald ratio using the inverse variance weighted (IVW) method as a main analysis to estimate the overall causal effect of genetically predicted values of the exposure on the outcome. For the IVW method, we used a multiplicative random-effects model when Cochran’s Q statistic (as described below) was significant (P < 0.05) ; otherwise, a fixed-effects model was used. Based on IVW results, we inferred the causal effect of the lifelong change in exposure on the outcome .
In addition, we conducted sensitivity analyses using the MR-Egger regression method, the weighted median method, the weighted mode method, and leave-one-out sensitivity analysis. The MR-Egger regression method can detect horizontal pleiotropy. The MR-Egger intercept is nonzero with statistical significance (P < 0.05) if possible horizontal pleiotropy of IVs exists . The weighted median method can generate a valid causal estimate if at least 50% of the instrument SNVs satisfy the IV assumptions . The weighted mode method forms clusters of individual SNVs and estimates the causal effect from the largest cluster . Leave-one-out sensitivity analysis removed each SNV from the IVW method and re-estimated the causal effect to assess the reliability of the analysis . We also measured heterogeneity among causal estimates across all SNVs in the IVW method by calculating Cochran’s Q statistic and the corresponding P value. Low heterogeneity (P > 0.05) provides more reliability for causal effects . Regarding MR to estimate the causal effect of TGs on COVID-19 severity, we added the following two analyses: we searched for SNVs associated with P < 5.0 × 10–8 with pleiotropic effects on BMI and any type of leukocyte or inflammatory marker using the web tool PhenoScanner (version 2) [35, 36] and then excluded them from the IVW method; we conducted the MR-PRESSO (Pleiotropy RESidual Sum and Outlier) global test with the run_mr_presso function to detect possible horizontal pleiotropy . Moreover, we conducted IVW multivariable MR analyses with the mv_ivw function to estimate the direct effects of genetically predicted TG levels on risk of COVID-19 independent of the effects of other exposure traits using genetically predicted Apo-B, LDL-C, or BMI traits as a covariate [13, 38,39,40,41].
We calculated statistical power in the MR analyses at a type-I error rate of 0.05 using the web tool mRnd [42, 43], as shown in Additional file 1: Table S1. For example, we achieved 80% power to detect an odds ratio (OR) of 1.066 (or 0.934) for the causal effect of genetically predicted Apo-B levels on COVID-19 susceptibility.
The characteristics of all SNVs included in our univariable MR analyses are shown in Additional file 1: Tables S2–S4. The F-statistic for every instrument was > 25, indicating no weak instrument bias. During harmonization, several SNVs were excluded because they were palindromic with intermediate MAFs (“palindromic ambiguous” was “TRUE” in Additional file 1: Tables S2–S4). The overall univariable MR results are shown in Table 1, Fig. 1, and Additional file 2: Figures S1–S3. None of the MR methods, including IVW, indicated any causal effect of genetically predicted Apo-B or LDL-C levels on risk of COVID-19, whereas some MR methods did indicate causal effects of genetically predicted TG levels (Table 1). By employing the IVW methods, we inferred that lifelong elevated TG levels had suggestive causal effects on a higher risk of COVID-19 susceptibility (OR per 1-SD increase in TGs, 1.065; 95% CI, 1.001–1.13; P = 0.045) and hospitalization (OR, 1.174; 95% CI, 1.04–1.33; P = 0.012). Moreover, IVW allowed us to infer that lifelong elevated TG levels had a significant causal effect on a higher risk of COVID-19 severity (OR, 1.274; 95% CI, 1.08–1.50; P = 0.004). The MR-Egger intercept indicated little evidence of horizontal pleiotropy (P = 0.75), and Cochran’s Q statistic for the IVW method indicated low heterogeneity (P = 0.43) and reliability of the causal effect. Leave-one-out sensitivity analysis (Additional file 2: Figure S3a) revealed the reliability of the IVW analysis, and the MR-PRESSO global test (P = 0.44) suggested a lack of possible horizontal pleiotropy. Furthermore, a funnel plot (Additional file 2: Figure S3b) depicted general symmetry, suggesting little evidence of heterogeneity or horizontal pleiotropy . The weighted median and weighted mode methods also showed OR scales and directions consistent with the IVW method; however, the effects were not significant, raising the possibility of horizontal pleiotropy and confounding. Therefore, we performed PhenoScanner searches to identify SNVs associated with possible pleiotropic effects on other risk factors for COVID-19 at P < 5.0 × 10–8. Among 208 SNVs used as IVs for serum TG levels, we identified 17 SNVs associated with BMI (rs10797119, rs10811662, rs10883026, rs11030107, rs13389219, rs2937124, rs2943645, rs3808477, rs3814883, rs394872, rs60856912, rs62135012, rs684773, rs7239575, rs7947951, rs921971, and rs998584), and 20 SNVs associated with any type of leukocyte or inflammatory marker (rs11187019, rs12185242, rs1292065, rs174566, rs2382825, rs28383314, rs326222, rs3775228, rs3814883, rs4665972, rs4761234, rs60856912, rs61905078, rs62427982, rs6432622, rs6999569, rs7239575, rs72603744, rs483082 [C-reactive protein (CRP)], and rs79287178 [TRAIL (tumor necrosis factor-related apoptosis inducing ligand) level]). We obtained results comparable to those of the original IVW method when we excluded the 17 SNVs or 20 SNVs (for the former: OR, 1.271; 95% CI, 1.07–1.51; P = 0.0059; number of SNVs, 191; for the latter: OR, 1.341; 95% CI, 1.10–1.64; P = 0.004; number of SNVs, 188).
Moreover, we carried out IVW multivariable MR analyses as sensitivity analyses to estimate the direct causal effect of genetically predicted TG levels on COVID-19 risk adjusted for each of the genetically predicted Apo-B, LDL-C, and BMI traits. The IVW multivariable MR results are shown in Table 2. The suggested causal effect of TGs on COVID-19 susceptibility was eliminated upon adjustment for each trait. We obtained comparable results regarding the causal effect of TGs on COVID-19 hospitalization and severity after adjustment for each trait; however, the latter had only suggestive significance.
Finally, we examined reverse causation by performing bidirectional two-sample univariable MR analyses using genetically predicted risks of COVID-19 as exposures and atherogenic lipid-related traits as outcomes. The results are shown in Additional file 3: Table S5. Overall, numbers of instrumental SNVs were low, and heterogeneities were high. We found that none of the COVID-19 risks had any significant causal effects on atherogenic lipid-related traits.
From univariable MR studies, we inferred suggestive causal effects of lifelong higher TG levels on higher risk of COVID-19 susceptibility (OR, 1.065; 95% CI, 1.001–1.13; P = 0.045) and hospitalization (OR, 1.174; 95% CI, 1.04–1.33; P = 0.012) and a significant causal effect of TGs on COVID-19 severity (OR, 1.274; 95% CI, 1.08–1.50; P = 0.004) but could not find any causal effect of Apo-B or LDL-C levels. The suggested effect of TGs on COVID-19 susceptibility was eliminated and the significant effect on COVID-19 severity was attenuated with suggestive significance after adjustment each for Apo-B, LDL-C, and BMI traits in IVW multivariable MR analyses. Nonetheless, our bidirectional MR study indicated that reverse causation of COVID-19 risk on atherogenic lipid-related traits was unlikely.
Observational studies have reported contradictory results. A study among 9005 UKBB participants (1508 patients testing positive for SARS-CoV-2 and 7497 controls) reported that Apo-B, LDL-C, and TGs were not significantly associated with SARS-CoV-2 infection . A systematic review and meta-analysis including 23 studies involving 10,122 COVID-19 patients showed that hospitalized patients with severe disease or non-survivor status had significantly lower serum LDL-C but not TG levels compared to patients with milder disease or survivor status; however, only a few studies of those with European ancestry were included . A retrospective single-center study with 654 patients in Spain showed that LDL-C < 69 mg/dl at admission was independently associated with a greater risk of 30-day mortality from COVID-19 (hazard ratio, 1.94; 95% CI, 1.14–3.31, P = 0.014) . However, a prospective single-center study with 48 COVID-19 patients in France did not detect a significant relationship between LDL-C and 28-day mortality . A retrospective single-center study with 600 COVID-19 patients in the United States reported hypertriglyceridemia as being associated with mortality (OR, 2.3; 95% CI, 1.4–3.7; P = 0.001) independent of obesity, high CRP, and high leukocyte count . Some observational studies of European-ancestry subjects were consistent with our MR results, though others were not. Such discrepancies may be due to the small sample size and retrospective and/or single-center designs of observational studies. Moreover, observational studies tend to suffer from bias by possible confounders and reverse causation . For example, serum lipid levels might be lower as a result of poor nutrition status due to COVID-19 severity . Associations between serum lipid levels and COVID-19 risks might be secondary to immune-inflammatory responses that could worsen COVID-19 outcomes [8, 9, 11], though the association of LDL-C with COVID-19 mortality was independent of inflammatory markers in the above Spanish study . MR studies can overcome such limitations of observational studies . Regardless, estimates of MR studies tend to be larger than those of observational studies, as the former estimate lifelong rather than short-term effects .
Some MR studies have investigated the causal effects of atherogenic lipid-related traits on COVID-19 risk. Leong A et al.  did not find any causal effect of LDL-C or TGs on COVID-19 susceptibility or severity that was significant after Bonferroni correction (LDL-C and susceptibility: OR, 1.09; 95% CI, 1.01–1.18; P = 0.04; LDL-C and severity: OR, 1.11; 95% CI, 0.98–1.26; P = 0.11; TG and susceptibility: OR, 1.05; 95% CI, 0.99–1.11; P = 0.08; TG and severity: OR, 1.05; 95% CI, 0.95–1.15; P = 0.36) using Global Lipids Genetic Consortium (GLGC) GWAS data for LDL-C and TG traits for 188,577 participants of European ancestry  and COVID-19-HGI GWAS meta-analyses data (Round 4). Aung N et al.  used two-sample MR with IVW methods to investigate causal effects of LDL-C and TGs on COVID-19 susceptibility (for LDL-C: OR, 1.20; 95% CI, 1.06–1.37, P = 0.006; for TG: OR, 1.10; 95% CI, 0.79–1.53; P = 0.58) using GLGC GWAS data  and COVID-19-HGI GWAS meta-analyses data (Round 2). Ponsford MJ et al.  reported no effect of LDL-C on COVID-19 with respiratory failure (OR, 0.85; 95% CI, 0.61–1.20) using GLGC GWAS data  and COVID-19 GWAS data involving 1610 cases with respiratory failure and 2205 controls in Italy and Spain . Recently, Zhang et al.  reported that Apo-B had a causal effect on COVID-19 susceptibility (OR, 1.18; 95% CI, 1.07–1.29; P = 0.001) using GWAS data of the Apo-B trait for 24,925 participants of European ancestry and UKBB data with 1221 COVID-19 cases and 4117 controls. However, they did not find any causal effect of Apo-B, LDL-C or TGs on COVID-19 severity (for Apo-B: OR, 0.95; 95% CI, 0.84–1.07; P = 0.36; for LDL-C: OR, 0.96; 95% CI, 0.85–1.09; P = 0.51; for TG: OR, 0.94; 95% CI, 0.78–1.15; P = 0.56) using LDL-C and TG traits GWAS for less than 300,000 participants of European ancestry and COVID-19 GWAS data obtained in Italy and Spain . We performed our MR study using the latest GWAS datasets with the largest sample sizes to date, which might explain the discrepancies between our and previous MR studies. In a preprint article, Sun Y et al.  performed an MR study using the UKBB dataset for the TG trait (given as GWAS-ID of “ukb-d-30870_irnt”), COVID-19-HGI GWAS meta-analyses data, and COVID-19 GWAS data obtained in Italy and Spain . Consistent with our MR study, they found that higher TGs had causal effects on the risk of severe COVID-19 (for 2972 very severe respiratory confirmed covid cases: beta, 0.219; P = 0.021; for 6492 hospitalized cases: beta, 0.224; P = 0.001; for 1610 cases with respiratory failure: beta, 0.446; P = 0.006).
The underlying mechanisms by which hypertriglyceridemia worsens COVID-19 outcomes clinically and pathologically remain unclear. Although an MR study showed that genetically predicted lower counts of basophils and myeloid white blood cells had causal effects on COVID-19 severity , an observational study found positive correlations of higher leukocyte counts and CRP with TGs in COVID-19 patients and suggested that hypertriglyceridemia might have a direct effect on COVID-19 severity due to an enhanced inflammatory response . In fact, studies suggested that hypertriglyceridemia could promote inflammation through leukocyte activation , macrophage accumulation in several organs  and increased sensitivity to cytokine stimulation of aortic endothelial cells . Moreover, an MR study indicated causal effects of cardiometabolic exposures, including BMI and TGs, on circulating proteins that might contribute to severe COVID-19 . For example, univariable MR indicated that both BMI and TGs had causal effects on reducing immunoglobulin G (IgG), a class of antibodies that help protect against infection. Nonetheless, multivariable MR indicated that BMI indirectly lowers IgG due to its influence on raising serum TG levels . Therefore, we infer that hypertriglyceridemia may worsen COVID-19 at least partly through the direct causal effect of TGs on inflammatory responses.
There are several major limitations to be noted in the present MR study. First, our MR analyses estimating a causal effect of TGs on COVID-19 severity using the weighted median and weighted mode methods had directionally consistent results but no statistical significance. Although so did other MR analyses [15, 37, 52], we must pay attention to possible horizontal pleiotropy and confounding. Therefore, we attempted to exclude possible pleiotropic effects by performing Phenoscanner searches. Previous MR studies have indicated that BMI is a risk factor for COVID-19 severity [14, 15, 37, 52], but CAD and diabetes are unlikely [14, 53]. Therefore, we excluded SNVs associated with BMI as well as inflammatory responses from the IVW method. Moreover, we conducted IVW multivariable MR analyses to eliminate the effect of Apo-B, LDL-C, or BMI traits. Regardless, we obtained results comparable to those of original IVW methods estimating causal effects of TGs on COVID-19 hospitalization and/or severity. Although the multivariable MR results had only suggestive significance, Bonferroni correction can be considered overly conservative, given the high correlation between lipid-related traits . Therefore, it is suggested that hypertriglyceridemia increases risk of COVID-19 hospitalization and severity to some extent independently of the effects of BMI, inflammatory responses, and other atherogenic lipid-related traits. Second, as described in the Methods section and Additional file 1: Table S1, our MR analyses estimating the causal effects of Apo-B and LDL-C on risk of COVID-19 might not possess sufficient statistical power to detect, if any existed, a weak association. Third, an observational study found a U-shaped association between LDL-C and COVID-19 severity . However, as two-sample MR analysis based on summary-level data assumes a linear relationship between exposure and outcome, we could not test for a nonlinear relationship between LDL-C and COVID-19 severity . Fourth, our MR analysis was based on populations of European ancestry, and the findings are unlikely to be generalized to other populations. Fifth, the GWAS of lipid-related traits was conducted, including blood samples during a non-fasting status that might affect serum lipid or lipoprotein levels [13, 55], which may have affected our causal estimate of TGs on the risks of COVID-19. Nonetheless, the GWAS and MR study adjusted for fasting time led to negligible changes in the effect estimates of Apo-B, LDL-C, and TGs on a higher risk of CAD .
Our two-sample MR approach indicated the causal effect of higher serum TG levels on a higher risk of COVID-19 severity in the European population using the most recent and largest GWAS datasets to date, suggesting that hypertriglyceridemia is a risk factor for COVID-19 severity. However, as the underlying mechanisms remain unclear and our MR study might be biased due to possible horizontal pleiotropy, further studies are warranted to validate our MR findings and investigate underlying mechanisms.
Availability of data and materials
The exposure GWAS dataset is publicly available from the IEU open GWAS database (https://gwas.mrcieu.ac.uk/datasets/) as GWAS-IDs of “ieu-b-108” for serum Apo-B level, “ieu-b-110” for serum LDL-C level, “ieu-b-111” for serum TG level, and “ieu-b-40” for BMI. The exposure GWAS dataset is publicly available from COVID-19-HGI GWAS meta-analyses (Round 5) (https://www.covid19hg.org/results/r5/). The functions used in the TwoSampleMR package are available from the website https://mrcieu.github.io/TwoSampleMR/reference/index.html.
Body mass index
Coronary artery disease
Coronavirus disease 2019
COVID-19 Host Genetics Initiative
Global Lipids Genetic Consortium
Genome-wide association study
Inverse variance weighted
Low-density lipoprotein cholesterol
Minor allele frequency
MR-Pleiotropy RESidual Sum and Outlier
Severe acute respiratory syndrome coronavirus 2
WHO Coronavirus Disease (COVID-19) Dashboard. https://covid19.who.int/ Accessed 13 Sept 2021.
Williamson EJ, Walker AJ, Bhaskaran K, Bacon S, Bates C, Morton CE, et al. Factors associated with COVID-19-related death using OpenSAFELY. Nature. 2020;584:430–6.
Gupta S, Hayek SS, Wang W, Chan L, Mathews KS, Melamed ML, et al. Factors associated with death in critically ill patients with coronavirus disease 2019 in the US. JAMA Intern Med. 2020;180:1436–47.
Juonala M, Magnussen CG, Berenson GS, Venn A, Burns TL, Sabin MA, et al. Childhood adiposity, adult adiposity, and cardiovascular risk factors. N Engl J Med. 2011;365:1876–85.
Ripatti P, Rämö JT, Mars NJ, Fu Y, Lin J, Söderlund S, et al. Polygenic hyperlipidemias and coronary artery disease risk. Circ Genom Precis Med. 2020;13:e002725.
Mooradian AD. Dyslipidemia in type 2 diabetes mellitus. Nat Clin Pract Endocrinol Metab. 2009;5:150–9.
Scalsky RJ, Chen Y-J, Desai K, O’Connell JR, Perry JA, Hong CC. Baseline cardiometabolic profiles and SARS-CoV-2 infection in the UK Biobank. PLoS ONE. 2021;16:e0248602.
Zinellu A, Paliogiannis P, Fois AG, Solidoro P, Carru C, Mangoni AA. Cholesterol and triglyceride concentrations, COVID-19 severity, and mortality: a systematic review and meta-analysis with meta-regression. Front Public Health. 2021;9:705916.
Aparisi Á, Iglesias-Echeverría C, Ybarra-Falcón C, Cusácovich I, Uribarri A, García-Gómez M, et al. Low-density lipoprotein cholesterol levels are associated with poor clinical outcomes in COVID-19. Nutr Metab Cardiovasc Dis. 2021;31:2619–27.
Tanaka S, De Tymowski C, Assadi M, Zappella N, Jean-Baptiste S, Robert T, et al. Lipoprotein concentrations over time in the intensive care unit COVID-19 patients: results from the ApoCOVID study. PLoS ONE. 2020;15:e0239573.
Dai W, Lund H, Chen Y, Zhang J, Osinski K, Jones SZ, et al. Hypertriglyceridemia during hospitalization independently associates with mortality in patients with COVID-19. J Clin Lipidol. 2021. https://doi.org/10.1016/j.jacl.2021.08.002.
Davey Smith G, Hemani G. Mendelian randomization: genetic anchors for causal inference in epidemiological studies. Hum Mol Genet. 2014;23:R89-98.
Richardson TG, Sanderson E, Palmer TM, Ala-Korpela M, Ference BA, Davey Smith G, et al. Evaluating the relationship between circulating lipoprotein lipids and apolipoproteins with risk of coronary heart disease: a multivariable Mendelian randomisation analysis. PLoS Med. 2020;17:e1003062.
Leong A, Cole JB, Brenner LN, Meigs JB, Florez JC, Mercader JM. Cardiometabolic risk factors for COVID-19 susceptibility and severity: a Mendelian randomization analysis. PLoS Med. 2021;18:e1003553.
Aung N, Khanji MY, Munroe PB, Petersen SE. causal inference for genetic obesity, cardiometabolic profile and COVID-19 susceptibility: a mendelian randomization study. Front Genet. 2020;11:586308.
Ponsford MJ, Gkatzionis A, Walker VM, Grant AJ, Wootton RE, Moore LSP, et al. Cardiometabolic traits, sepsis, and severe COVID-19: a mendelian randomization investigation. Circulation. 2020;142:1791–3.
Zhang K, Dong SS, Guo Y, Tang SH, Wu H, Yao S, et al. Causal associations between blood lipids and COVID-19 risk: a two-sample Mendelian randomization study. Arterioscler Thromb Vasc Biol. 2021. https://doi.org/10.1161/ATVBAHA.121.316324.
R Core Team. R: a language and environment for statistical computing. Vienna: R Foundation for Statistical Computing; 2020.
UK Biobank: Data-Field 30890. https://biobank.ctsu.ox.ac.uk/crystal/field.cgi?id=30870. Accessed 6 Sept 2021.
MRC IEU Open GWAS Project. https://gwas.mrcieu.ac.uk/. Accessed 1 Sept 2021.
COVID-19 Host Genetics Initiative. Mapping the human genetic architecture of COVID-19. Nature; 2021. https://doi.org/10.1038/s41586-021-03767-x.
The COVID-19 Host Genetics Initiative. https://www.covid19hg.org/results/r5/ Accessed 3 Feb 2021.
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.
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.
Xiuyun W, Qian WX, Weidong L, Lizhen L. Education and stroke: evidence from epidemiology and Mendelian randomization study. Sci Rep. 2020;10:21208.
Gill D, Karhunen V, Malik R, Dichgans M, Sofat N. Cardiometabolic traits mediating the effect of education on osteoarthritis risk: a Mendelian randomization study. Osteoarthritis Cartilage. 2021;29:365–71.
Liu HM, Hu Q, Zhang Q, Su GY, Xiao HM, Li BY, et al. Causal effects of genetically predicted cardiovascular risk factors on chronic kidney disease: a two-sample Mendelian randomization study. Front Genet. 2019;10:415.
Burgess S, Dudbridge F, Thompson SG. Combining information on multiple instrumental variables in Mendelian randomization: comparison of allele score and summarized data methods. Stat Med. 2016;35:1880–906.
Burgess S, Small D, Thompson S. A review of instrumental variable estimators for Mendelian randomization. Stat Methods Med Res. 2017;26:2333–55.
Zhao JV, Schooling CM. Effect of linoleic acid on ischemic heart disease and its risk factors: a Mendelian randomization study. BMC Med. 2019;17:61.
Burgess S, O’Donnell CJ, Gill D. Expressing results from a mendelian randomization analysis: separating results from inferences. JAMA Cardiol. 2021;6:7–8.
Walker VM, Davies NM, Hemani G, Zheng J, Haycock PC, Gaunt TR, et al. Using the MR-Base platform to investigate risk factors and drug targets for thousands of phenotypes. Wellcome Open Res. 2019;4:113.
Burgess S, Davey Smith G, Davies NM, Dudbridge F, Gill D, Glymour MM, et al. Guidelines for performing Mendelian randomization investigations. Wellcome Open Res. 2020;4:186.
Greco MFD, Minelli C, Sheehan NA, Thompson JR. Detecting pleiotropy in Mendelian randomisation studies with summary data and a continuous outcome. Stat Med. 2015;34:2926–40.
PhenoScanner V2. A database of human genotype-phenotype associations. http://www.phenoscanner.medschl.cam.ac.uk/. Accessed 6 Sept 2021.
Staley JR, Blackshaw J, Kamat MA, Ellis S, Surendran P, Sun BB, et al. PhenoScanner: a database of human genotype-phenotype associations. Bioinformatics. 2016;32:3207–9.
Freuer D, Linseisen J, Meisinger C. Impact of body composition on COVID-19 susceptibility and severity: a two-sample multivariable Mendelian randomization study. Metabolism. 2021;118:154732.
Yuan S, Michaëlsson K, Wan Z, Larsson SC. Associations of smoking and alcohol and coffee intake with fracture and bone mineral density: a Mendelian randomization study. Calcif Tissue Int. 2019;105:582–8.
Xu H, Jin C, Guan Q. Causal effects of overall and abdominal obesity on insulin resistance and the risk of type 2 diabetes mellitus: a two-sample Mendelian randomization study. Front Genet. 2020;2(11):603.
Adams CD, Boutwell BB. Can increasing years of schooling reduce type 2 diabetes (T2D)?: Evidence from a Mendelian randomization of T2D and 10 of its risk factors. Sci Rep. 2020;10:12908.
Lu Y, Guo Y, Lin H, Wang Z, Zheng L. Genetically determined tobacco and alcohol use and risk of atrial fibrillation. BMC Med Genomics. 2021;14:73.
mRnd: Power calculations for Mendelian Randomization. https://shiny.cnsgenomics.com/mRnd/ Accessed 18 Oct 2021.
Brion MJ, Shakhbazov K, Visscher PM. Calculating statistical power in Mendelian randomization studies. Int J Epidemiol. 2013;42:1497–501.
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–83.
Severe Covid-19 GWAS Group, Ellinghaus D, Degenhardt F, Bujanda L, Buti M, Albillos A, et al. Genomewide association study of severe Covid-19 with respiratory failure. N Engl J Med. 2020;383:1522–34.
Sun Y, Zhou J, Ye K. Prioritizing causal risk factors for severe COVID-19: an exhaustive Mendelian randomization study. 25 Jan 2021 (Version 1). https://doi.org/10.21203/rs.3.rs-149087/v1
Sun Y, Zhou J, Ye K. White blood cells and severe COVID-19: a Mendelian randomization study. J Pers Med. 2021;11:195.
Alipour A, van Oostrom AJ, Izraeljan A, Verseyden C, Collins JM, Frayn KN, et al. Leukocyte activation by triglyceride-rich lipoproteins. Arterioscler Thromb Vasc Biol. 2008;28:792–7.
Saja MF, Baudino L, Jackson WD, Cook HT, Malik TH, Fossati-Jimack L, et al. Triglyceride-rich lipoproteins modulate the distribution and extravasation of Ly6C/Gr1(low) monocytes. Cell Rep. 2015;12:1802–15.
Ting HJ, Stice JP, Schaff UY, Hui DY, Rutledge JC, Knowlton AA, et al. Triglyceride-rich lipoproteins prime aortic endothelium for an enhanced inflammatory response to tumor necrosis factor-alpha. Circ Res. 2007;100:381–90.
Richardson TG, Fang S, Mitchell RE, Holmes MV, Davey SG. Evaluating the effects of cardiometabolic exposures on circulating proteins which may contribute to severe SARS-CoV-2. EBioMedicine. 2021. https://doi.org/10.1016/j.ebiom.2021.103228.
Li S, Hua X. Modifiable lifestyle factors and severe COVID-19 risk: a Mendelian randomisation study. BMC Med Genomics. 2021;14:38.
Au Yeung SL, Zhao JV, Schooling CM. Evaluation of glycemic traits in susceptibility to COVID-19 risk: a Mendelian randomization study. BMC Med. 2021;19:72.
Gong J, Chen Y, Jie Y, Tan M, Jiang Z, Yuan L, et al. U-shaped relationship of low-density lipoprotein cholesterol with risk of severe COVID-19 from a multicenter pooled analysis. Front Cardiovasc Med. 2021;8:604736.
Yuan S, Tang B, Zheng J, Larsson SC. Circulating lipoprotein lipids, apolipoproteins and ischemic stroke. Ann Neurol. 2020;88:1229–36.
We highly appreciate the GWAS  conducted by Richardson TG et al. using UK Biobank data. We would like to thank the MRC IEU Open GWAS database and the COVID-19 Host Genetics Initiative for making GWAS datasets publicly available.
Ethics approval and consent to participate
As all data used in the present study were derived from publicly available summary-level GWAS datasets and no individual-level data were used, additional ethical approval and patient consent were not necessary.
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Table S1. Statistical power for estimating effects of genetically predicted lipid-lipoprotein traits on COVID-19 risk in MR analyses. Table S2. The characteristics of all SNVs included in MR analyses estimating the causal effect of genetically predicted serum apolipoprotein B levels (exposure) on risk of COVID-19 susceptibility (outcome1), hospitalization (outcome2), and severity (outcome3). Table S3. The characteristics of all SNVs included in MR analyses estimating the causal effect of genetically predicted serum LDL-cholesterol levels (exposure) on risk of COVID-19 susceptibility (outcome1), hospitalization (outcome2), and severity (outcome3). Table S4. The characteristics of all SNVs included in MR analyses estimating the causal effect of genetically predicted serum triglyceride levels (exposure) on risk of COVID-19 susceptibility (outcome1), hospitalization (outcome2), and severity (outcome3).
Figure S1. Scatter plots for estimating causal effects of genetically predicted serum Apo-B levels on risk of (a) susceptibility, (b) hospitalization, and (c) severity. Figure S2. Scatter plots for estimating causal effects of genetically predicted serum LDL-C levels on risk of (a) susceptibility, (b) hospitalization, and (c) severity. Figure S3. (a) Leave-one-out sensitivity analysis for estimating causal effect of genetically predicted serum TG levels on risk of COVID-19 severity. (b) Funnel plot for estimating causal effect of serum TG levels on risk of COVID-19 severity.
Table S5. Univariable MR results of the effect of risk of COVID-19 on atherogenic lipid-related traits.
About this article
Cite this article
Yoshikawa, M., Asaba, K. & Nakayama, T. Estimating causal effects of atherogenic lipid-related traits on COVID-19 susceptibility and severity using a two-sample Mendelian randomization approach. BMC Med Genomics 14, 269 (2021). https://doi.org/10.1186/s12920-021-01127-2
- Apolipoprotein B
- Mendelian randomization