Mendelian randomization rules out the causal relationship between serum lipids and cholecystitis
BMC Medical Genomics volume 14, Article number: 224 (2021)
The relationship between serum lipids and cholecystitis is still under investigation. To examine the causal effect of serum lipids on cholecystitis using the Mendelian randomization method.
We conducted univariable Mendelian randomization (MR) analyses using summary statistics from two independent genome-wide association studies (GWAS) on serum lipids (n = 132,908) and cholecystitis (n = 361,194). Mainly, the inverse-variance weighted (IVW) method was utilized to combine each SNP’s causal estimation, and the MR-Egger was adopted as a complementary method, together with the weighted median. Cochrane’s Q value was employed to appraise heterogeneity. The MR-Egger intercept and MR-PRESSO were used to detect the horizontal pleiotropy.
Our univariable results displayed a minor protective effect of serum low-density lipoprotein (LDL) cholesterol (OR [95% CI] = 0.9984483 [0.9984499, 0.9984468]; p = 0.008) on cholecystitis. No significant causal effect of total cholesterol (TC) (OR [95% CI] = 0.9994228 [0.9994222, 0.9994233]; p = 0.296), triglycerides (OR [95% CI] = 0.9990893 [0.9990882, 0.9990903]; p = 0.238) and high-density lipoprotein (HDL) cholesterol (OR [95% CI] = 0.9997020 [0.9997017, 0.9997023]; p = 0.565) was found on cholecystitis.
These findings suggest that LDL cholesterolhas a slight protective effect on cholecystitis, which can be easily affected by confounding factors. TC, triglycerides and HDL cholesterol don’t have causal effect on cholecystitis. The protective effect of serum lipids on cholecystitis, though possible, remain less certain.
Cholecystitis is one of the most common diseases of the gallbladder. According to a comprehensive survey, 262,411 hospitalizations occurred in 2000 for cholecystitis . Fortunately, the mortality of cholecystitis is only at 0.6% . Cholecystitis could be roughly classified as acute cholecystitis acalculous cholecystitis and chronic cholecystitis. Cholecystectomy is the gold standard treatment of all kinds of cholecystitis. Although it is not a big surgery, it puts a financial and physical burden on patients. Therefore, the prevention of cholecystitis is of vital importance. Many risk factors were reported by the previous study such as obesity, old age, female gender, long-term total parenteral nutrition, diabetes, etc. . However, whether serum lipids are risk factors for cholecystitis is still under investigation. Furthermore, the current results are controversial. For example, in univariable analyses conducted by Mohr et al. , they found a negative association between high-density lipoprotein (HDL) cholesterol, low-density lipoprotein (LDL) cholesterol, and cholecystitis while a case–control study conducted by Wang et al.  reported that HDL cholesterol was positively associated with cholecystitis. Thus, we hope to find out the causal relationship between blood lipids and cholecystitis. However, it is difficult to conduct a randomized control study (RCT) to investigate the causal relationship. Therefore, we carried out a Mendelian randomization analysis to overcome the limitations of the previous study.
Mendelian randomization (MR) is a method using genetic variants to determine whether there exists a causal relationship between the exposure (usually a risk factor) and the outcome [5, 6]. The random assortment of genetic variants follows the law of Mendelian. As a result, the allocation of exposure for each individual is achieved randomly, which is the fundament of MR study. Numerous MR studies have been published in recent years where some of them even overturned widely accepted opinions. For example, the MR analysis of Holmes et al.  ruled out the causal effect of HDL cholesterol on coronary heart disease. MR analysis is a good way to overcome difficulties of conventional research such as confounders, loss of follow-up, time-consuming. The effect of the instrumental variables (IV) lasts for the whole life even before the birth of the individual. It is exactly suitable for the investigation of the lifetime exposure’s effect on the target disease. Compared to conventional observational studies, MR analysis is more convincing and reliable.
Here, we used single nucleotide polymorphisms (SNP) associated with TC, HDL cholesterol, LDL cholesterol, and triglyceride as instrument variables to conduct MR analysis. We extracted SNPs from the genome-wide association studies (GWAS) from the GLGC consortium to investigate the effect of blood lipids on cholecystitis . The levels of TC, HDL cholesterol, LDL cholesterol, and triglycerides were treated as the exposure separately in the univariable MR analysis to investigate its causal effect on cholecystitis. Thus, we hope to disentangle the complex relationship between blood lipids and cholecystitis.
Data involved in this study are publicly available. We extracted the exposure data from the Global Lipids Genetics Consortium (GLGC, http://lipidgenetics.org/). This GWAS study contained summary statistics from 45 studies. Among them, 37 studies consisted of individuals of European ancestry (n = 114,230). Another 9 studies consisted primarily of individuals with non-European ancestry: two studies of South Asian descent (n = 1516, n = 3385); two studies of East Asian descent (n = 1771, n = 7044); five studies of recent African ancestry, from Uganda (n = 1687), from the Caribbean (n = 426), and the United States (n = 1614, n = 397, n = 838). The exposure data was given in per SD unit. To identify outlier studies, average standard errors for association statistics from each study were plotted against the study sample size. Allele frequencies were scrutinized to guarantee all analyses using the same strand assignment of alleles. Reported statistics and allelic effects were consistent with published findings for known loci. Genomic control values for study-specific analyses were < 1.20. Variants whose minor alleles were observed < 7 times were excluded. Phenotype distribution, proportion of phenotype variance explained by instruments were provided in Table 1.
The summary statistics associated with cholecystitis were extracted from the UK biobank (www.nealelab.is/uk-biobank). They filtered from 487,409 individuals down to 337,199 individuals. Individuals included in the GWAS analysis of the UK biobank were restricted to European descent. Individuals who met any of the following criteria were excluded: (1) closely related individuals, (2) individuals with sex chromosome aneuploidies, (3) individuals who had withdrawn consent from the UK biobank study. Over 92 million imputed autosomal SNPs were available for the GWAS analysis. SNPs with minor allele frequency > 0.1% and Hardy–Weinberg equilibrium (HWE) p value > 1 × 10−10 were included. Finally, 10.8 million SNPs remained for analysis.
IV selection (F statistics)
All the SNPs were tested for association with the trait of interest at a genome-wide significant level (p < 5 × 10−8). SNPs missing in summary statistics of outcome were removed. For the sake of reducing the missing SNPs we selected proxy SNPs (European populations) in high LD (r2 > 0.8). We excluded the weak instrument variable. The weak instrument variable was defined as the F statistic less than 10. The F statistic is usually quoted as a measure of the strength of an instrument variable. SNPs that failed to pass the leave-one-out analysis were eliminated. All the SNPs were scrutinized for linkage disequilibrium (LD) according to the criterion as r2 < 0.01, kilobase (kb) > 10,000, clumping threshold < 5 × 10−8. Palindromic SNPs were excluded from this research. MR-PRESSO global test, outlier test, and distortion test were conducted to identify and remove SNPs with horizontal pleiotropy. Compared to other tests of pleiotropy, all three tests of MR-PRESSO have higher power to detect horizonal pleiotropy except for the perfectly overlapping samples . After series of strict filtration the rest of SNPs were regarded as qualified IVs.
Mendelian randomization relies on three assumptions (Fig. 1): (1) the instrument variable (genetic variant) is associated with the risk factor; (2) the instrument variable is not associated with confounders; (3) the instrument variable influences the outcome only by the risk factor. A genetic variant satisfying these assumptions is known as an instrumental variable. SNPs involved in this study could easily fit in assumption 1. We extracted SNPs related to TC, HDL cholesterol, LDL cholesterol and triglycerides as instrumental variables while the SNPs related to cholecystitis were treated as the outcome. With a SNP used as an IV and a biniary outcome, assuming all associations are log-linear, the causal effect of the exposure on the outcome can be estimated as the ratio of the change in the outcome per additional variant allele divided by the change in the risk factor per additional variant allele . Analyses were performed using R, version 4.0.3 (http: //www.r-project.org). R package ‘TwoSampleMR’ and ‘MRPRESSO’ were employed in this study. The R code is provided in this Github repository (https://github.com/YangHQ638/MR-research).
Many specific methods of MR analysis are available. We adopted the following three methods: inverse-variance weighted (IVW), MR-Egger, and weighted median. IVW is the main method adopted in both MR analysis and multivariable MR analysis. IVW estimation is the inverse variance weighted combination of ratio estimates  where the random effect model was adopted if there was heterogeneity. Wald ratio is used to obtain a ratio from each SNP. MR-Egger can test and estimate the causal effect. It can also test the directional pleiotropy, making it preferred in the case of pleiotropy. The estimation of MR-Egger is close to that of IVW when the intercept term of MR Egger regression is close to 0 . Compared to IVW (fixed-effect model), the weighted median method could correctly estimate the causality when some of the instrumental variables are invalid . Thus, it is complementary to MR Egger. We have adopted all methods for all exposures. The statistical power of this study was tested by an online tool (https://cnsgenomics.shinyapps.io/mRnd/). To address the issue of multiple testing, results were considered statistically significant at the 0.05 level after a false discovery rate approach of the significance level. We conducted a Cochran Q test to examine the heterogeneity of these analyses. We also tested the pleiotropy by collecting the intercept values from MR Egger regression. The leave-one-out sensitivity analysis was performed to test the robustness of this study.
We made a scatter plot for the result of each exposure. Different lines of scatter plots represent the causal estimation of different methods. Results of sensitivity analysis were visualized by a leave-one-out plot. The pooled result of all SNPs was shown in the forest plot. Biased was tested by funnel plot .
The study identified 111 SNPs associated with serum TC level, 123 SNPs associated with HDL cholesterol level, 78 SNPs associated with LDL cholesterol level, and 71 SNPs associated with triglycerides level. All the SNPs were associated with the trait of interest at a genome-wide significant level (p < 5 * 10−8). Since the summary statistics were missing in the outcome, 2 SNPs related to TC and 2 SNPs related to LDL cholesterol were removed. None of the SNPs was weak instrumental variable (F statistics < 10). 16 SNPs associated with LDL cholesterol were omitted by leave-one-out analysis. A total of 5 palindromic SNPs were excluded. Due to pleiotropy, we removed 3 SNPs related to TC and 2 SNPs related to LDL in the MR-PRESSO test. Excluded SNPs, numbers of identified SNPs, were showed in Table 2. Specific identified SNPs of each phenotype were provided in the Additional materials [see Additional file 1].
The causal estimations of MR analysis were provided in Table 3. The minor causal effect of LDL cholesterol on cholecystitis was demonstrated by IVW (multiplicative random effect model) (OR [95% CI] = 0.9984483 [0.9984499, 0.9984468]; p = 0.008). There was an evidence of heterogeneity (Q = 116.489, p value = 0.004). No evidence was found about pleiotropy (egger intercept = − 2.1 × 10−5, p value = 0.659). Although it was not significant at the level of p = 0.05, a minor protective effect was observed in TC (OR [95% CI] = 0.9994228 [0.9994222, 0.9994233]; p = 0.296), triglycerides (OR [95% CI] = 0.9990893 [0.9990882, 0.9990903]; p = 0.238) and HDL cholesterol (OR [95% CI] = 0.9997020 [0.9997017, 0.9997023]; p = 0.565). Heterogeneity was detected in HDL cholesterol (Q = 188.737, p value = 0.000) and triglycerides (Q = 96.010, p value = 0.021) while it was not detected in TC (Q = 135.009, p value = 0.053). The MR-Egger intercept suggested there was no evidence of pleiotropy in TC (egger intercept = 4.06 × 10−5, p value = 0.335), HDL cholesterol (egger intercept = 4.46 × 10−5, p value = 0.319), and triglycerides (egger intercept = 6.55 × 10−5, p value = 0.168). The statistical power of these four exposures was relatively low (TC: 0.08, HDL cholesterol: 0.10, LDL cholesterol: 0.08, triglycerides: 0.08). A Scatter plot of MR results was shown in Fig. 2. Forest plots of each SNP were provided in the Additional materials [see Additional file 2].
The funnel plots showed the extent of heterogeneity among the individual Wald ratio estimates. The forest plots displayed the effect estimates of each SNP. The leave-one-out plots showed each SNP was robust in MR analysis. The funnel plots, forest plots and leave-one-out plots were provided in the Additional materials [see Additional file 2].
We conducted univariable MR analyses to investigate the causal effect of serum lipids on cholecystitis. The current study found that LDL cholesterol was protective factors against cholecystitis but such protective effect may be subtle because odds ratios were close to 1. The results of HDL cholesterolTC and triglycerides were not statistically significant but all of them showed a protective tendency. The causal association of serum lipids with cholecystitis remains possible, but less certain. Statistical power was relatively low, but nevertheless, the protective effect of LDL cholesterol could still be demonstrated in current study. LDL cholesterol and triglycerides showed a similar trend as a protective factor of cholecystitis, which was partially consistent with some of the previous findings. In univariate analyses conducted by Mohr et al. , they found an inverse relationship between HDL cholesterol, LDL cholesterol, and cholecystitis but they also found LDL cholesterol and triglycerides were positively related while HDL cholesterol was not related to cholecystitis in multiple regression. A case–control study conducted by Wang et al.  reported that LDL cholesterol was reversely related while HDL cholesterol was positively related. Wang et al. also found a low level of TC increased the risk of gallstone disease . Gallstones can directly induce cholecystitis. Several reports have shown that cholesterol played an important role in the process of gallstone formation [3, 14, 15]. Jiang et al. also reported that patients with gallstone disease had lower serum TC . However, Another case–control study failed to find an association between HDL cholesterol and gallbladder disease . The reason why there are different opinions is that the protective effect of specific serum lipid is really small (but they do exist which is demonstrated in this research). Even some confounding factors are enough to affect the result.
A possible explanation for the protective effect might be the differences in cholesterol transport and secretion. Clinical studies have demonstrated an inverse relationship between serum HDL cholesterol level and bile cholesterol saturation . Lithogenic bile with supersaturated cholesterol is the initial event of the pathogenesis of gallstones and chronic cholecystitis [19, 20]. Supersaturated concentrations of biliary cholesterol cause abnormalities in the mucosa and muscle layers of the gallbladder . They are the results of increase transport of cholesterol across the gallbladder mucosa due to the high concentration in gallbladder , which is the result of increasing cholesterol transport from serum to gallbladder. This phenomenon may cause the decrease of serum cholesterol and increase of supersaturated cholesterol biliary simultaneously. Analysis of inbred mice also showed a similar result . Some biochemical analyses have found many molecules related to the transport of cholesterol. SR-BI (scavenger receptor class B type I) is an important determinant of plasma-to-bile transport of HDL cholesterol. Mice with overexpression of SR-BI enhanced biliary cholesterol levels while a decrease of serum cholesterol level was observed in SR-BI deficiency mice [22,23,24,25]. SR-BI might upregulate in cholecystitis patients which needs further research to verify. Juvonen et al. reported TaqBI CEPT gene polymorphism, which is correlated with lower serum TC, was associated with cholesterol gallstone disease . It has been demonstrated that gallstone patients secret more cholesterol than normal people [19, 26]. Gallstones could directly lead to cholecystitis. The absolute decrease in serum TC might be a crucial factor in determining cholesterol saturation in bile . However, a long-term high cholesterol diet might cause an increase in serum TC and indirectly increased the secretion of cholesterol . Nevertheless, the way of instrumental variable affect exposure was not the same as the way of diet. So, if we could exclude people with high cholesterol diet, then the result would be clearer. This might be the reason that odds ratio of serum lipids was near to 1. Along with instrumental variables, confounders and environmental factors may co-influence the final result. It is important to conduct a stratified analysis for future research.
Another explanation for our result might be that some of these lipids accelerate gallbladder emptying. Hopman et al.  have demonstrated that long-chain triglycerides increased the concentration of plasma cholecystokinin (CCK), which could enhance gallbladder motility while medium-chain triglycerides were not. This might explain that the protective effect of triglycerides was very small. Our results have demonstrated that triglyceride was a minor protective factor. There might be some unknown interactions between different types of cholesterol. Further analyses need to find out the effect of different types of cholesterol and the effect of triglycerides on the gallbladder.
There were still a few drawbacks to our study. Firstly, canalization might have an impact on our result. Due to the whole life-span exposure to low levels of plasma lipids, the body might emerge a negative feedback mechanism to compensate. Secondly, because the instrument variable only could slightly affect the exposure, our result was based on extrapolation. Thirdly, our study was based on a mixed population. Because of different races, cultures, diets, and climate, a study based on a specific population might more persuasive. Moreover, recent epidemiological studies reported women are at a higher risk than men , while others believed not . Due to limitations of the Mendelian randomization study, we didn’t investigate the differences between genders. Further studies should be based on a stratified population to overcome this difficulty. Finally, limited by the data source, the statistical power of this research was quite low, hence typeIIerror might occurs in this research. However, it didn’t affect the conclusion that LDL cholesterol was a protective factor against cholecystitis.
Furthermore, previous studies have investigated or stated hypotheses a few underlying mechanisms of lipids affecting cholecystitis. Measuring serum lipids is much easier than measuring cholesterol in the gallbladder in clinical practice. Future studies should not only focus on the relationship between cholesterol saturation and cholecystitis but also investigate the relationship between serum lipids and cholecystitis.
In conclusion, the minor protective effect of LDL cholesterol has been demonstrated in this MR analysis. The protective effect of HDL cholesterol, triglycerides and TC on cholecystitis were also found but not statistically significant in this study. However, the protective effect demonstrated in this study is negligible. The protective effect of serum lipids on cholecystitis, though possible, remain less certain.
Availability of data and materials
The data used to support the findings of this study are included within the article.
Single nucleotide polymorphisms
Genome-wide association studies
Shaffer EA. Gallstone disease: epidemiology of gallbladder stone disease. Best Pract Res Clin Gastroenterol. 2006;20(6):981–96.
Elwood DR. Cholecystitis. Surg Clin North Am. 2008;88(6):1241–52, viii.
Mohr GC, Kritz-Silverstein D, Barrett-Connor E. Plasma lipids and gallbladder disease. Am J Epidemiol. 1991;134(1):78–85.
Wang J, Shen S, Wang B, Ni X, Liu H, Ni X, Yu R, Suo T, Liu H. Serum lipid levels are the risk factors of gallbladder stones: a population-based study in China. Lipids Health Dis. 2020;19(1):50.
Emdin CA, Khera AV, Kathiresan S. Mendelian randomization. JAMA. 2017;318(19):1925–6.
Smith GD, Ebrahim S. “Mendelian randomization”: can genetic epidemiology contribute to understanding environmental determinants of disease? Int J Epidemiol. 2003;32(1):1–22.
Holmes MV, Asselbergs FW, Palmer TM, Drenos F, Lanktree MB, Nelson CP, Dale CE, Padmanabhan S, Finan C, Swerdlow DI, et al. Mendelian randomization of blood lipids for coronary heart disease. Eur Heart J. 2015;36(9):539–50.
Willer CJ, Schmidt EM, Sengupta S, Peloso GM, Gustafsson S, Kanoni S, Ganna A, Chen J, Buchkovich ML, Mora S, et al. Discovery and refinement of loci associated with lipid levels. Nat Genet. 2013;45(11):1274–83.
Verbanck M, Chen CY, 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(5):693–8.
Thomas DC, Conti DV. Commentary: the concept of “Mendelian randomization.” Int J Epidemiol. 2004;33(1):21–5.
Burgess S, Butterworth A, Thompson SG. Mendelian randomization analysis with multiple genetic variants using summarized data. Genet Epidemiol. 2013;37(7):658–65.
Burgess S, Thompson SG. Interpreting findings from Mendelian randomization using the MR-Egger method. Eur J Epidemiol. 2017;32(5):377–89.
Egger M, Davey Smith G, Schneider M, Minder C. Bias in meta-analysis detected by a simple, graphical test. BMJ. 1997;315(7109):629–34.
Juvonen T, Savolainen MJ, Kairaluoma MI, Lajunen LH, Humphries SE, Kesäniemi YA. Polymorphisms at the apoB, apoA-I, and cholesteryl ester transfer protein gene loci in patients with gallbladder disease. J Lipid Res. 1995;36(4):804–12.
Fuchs M, Lammert F, Wang DQ, Paigen B, Carey MC, Cohen DE. Sterol carrier protein 2 participates in hypersecretion of biliary cholesterol during gallstone formation in genetically gallstone-susceptible mice. Biochem J. 1998;336(Pt 1):33–7.
Jiang ZY, Han TQ, Suo GJ, Feng DX, Chen S, Cai XX, Jiang ZH, Shang J, Zhang Y, Jiang Y, et al. Polymorphisms at cholesterol 7alpha-hydroxylase, apolipoproteins B and E and low density lipoprotein receptor genes in patients with gallbladder stone disease. World J Gastroenterol. 2004;10(10):1508–12.
Scragg RK, Calvert GD, Oliver JR. Plasma lipids and insulin in gall stone disease: a case-control study. Br Med J (Clin Res Ed). 1984;289(6444):521–5.
Thornton JR, Heaton KW, Macfarlane DG. A relation between high-density-lipoprotein cholesterol and bile cholesterol saturation. Br Med J (Clin Res Ed). 1981;283(6303):1352–4.
van Erpecum KJ, Wang DQ, Moschetta A, Ferri D, Svelto M, Portincasa P, Hendrickx JJ, Schipper M, Calamita G. Gallbladder histopathology during murine gallstone formation: relation to motility and concentrating function. J Lipid Res. 2006;47(1):32–41.
Jacyna MR, Ross PE, Bakar MA, Hopwood D, Bouchier IA. Characteristics of cholesterol absorption by human gall bladder: relevance to cholesterolosis. J Clin Pathol. 1987;40(5):524–9.
Khanuja B, Cheah YC, Hunt M, Nishina PM, Wang DQ, Chen HW, Billheimer JT, Carey MC, Paigen B. Lith1, a major gene affecting cholesterol gallstone formation among inbred strains of mice. Proc Natl Acad Sci U S A. 1995;92(17):7729–33.
Mardones P, Quiñones V, Amigo L, Moreno M, Miquel JF, Schwarz M, Miettinen HE, Trigatti B, Krieger M, VanPatten S, et al. Hepatic cholesterol and bile acid metabolism and intestinal cholesterol absorption in scavenger receptor class B type I-deficient mice. J Lipid Res. 2001;42(2):170–80.
Trigatti B, Rayburn H, Viñals M, Braun A, Miettinen H, Penman M, Hertz M, Schrenzel M, Amigo L, Rigotti A, et al. Influence of the high density lipoprotein receptor SR-BI on reproductive and cardiovascular pathophysiology. Proc Natl Acad Sci U S A. 1999;96(16):9322–7.
Van Eck M, Twisk J, Hoekstra M, Van Rij BT, Van der Lans CA, Bos IS, Kruijt JK, Kuipers F, Van Berkel TJ. Differential effects of scavenger receptor BI deficiency on lipid metabolism in cells of the arterial wall and in the liver. J Biol Chem. 2003;278(26):23699–705.
Wang DQ, Carey MC. Susceptibility to murine cholesterol gallstone formation is not affected by partial disruption of the HDL receptor SR-BI. Biochim Biophys Acta. 2002;1583(2):141–50.
Salen G, Nicolau G, Shefer S, Mosbach EH. Hepatic cholesterol metabolism in patients with gallstones. Gastroenterology. 1975;69(3):676–84.
von Bergmann K, Leiss O. Effect of short-term treatment with bezafibrate and fenofibrate on biliary lipid metabolism in patients with hyperlipoproteinaemia. Eur J Clin Investig. 1984;14(2):150–4.
Kakimoto T, Kanemoto H, Fukushima K, Ohno K, Tsujimoto H. Effect of a high-fat-high-cholesterol diet on gallbladder bile acid composition and gallbladder motility in dogs. Am J Vet Res. 2017;78(12):1406–13.
Hopman WP, Jansen JB, Rosenbusch G, Lamers CB. Effect of equimolar amounts of long-chain triglycerides and medium-chain triglycerides on plasma cholecystokinin and gallbladder contraction. Am J Clin Nutr. 1984;39(3):356–9.
Behar J, Mawe GM, Carey MC. Roles of cholesterol and bile salts in the pathogenesis of gallbladder hypomotility and inflammation: cholecystitis is not caused by cystic duct obstruction. Neurogastroenterol Motil. 2013;25(4):283–90.
Support in writing manuscript and data analysis was provided by Bai Ji and Qiong Yu. We thank C. J. Willer et al. who kindly provided the GWAS data for our work.
This research did not receive any specific grant from funding agencies in the public, commercial, or not-for-profit sectors.
Ethics approval and consent to participate
Consent for publication
All the authors have consented for publication.
The authors declare that they have no competing interest.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Yang, H., Chen, L., Liu, K. et al. Mendelian randomization rules out the causal relationship between serum lipids and cholecystitis. BMC Med Genomics 14, 224 (2021). https://doi.org/10.1186/s12920-021-01082-y