Unraveling the link between sterol ester and colorectal cancer: a two-sample Mendelian randomization study

doi:10.21203/rs.3.rs-4369169/v1

Download PDF

Research Article

Unraveling the link between sterol ester and colorectal cancer: a two-sample Mendelian randomization study

https://doi.org/10.21203/rs.3.rs-4369169/v1

This work is licensed under a CC BY 4.0 License

You are reading this latest preprint version

Background

Several studies reported the sterol ester (SE), one subclass of subtype of cholesterol esters (CE), is associated with the incidence of Colorectal cancer (CRC). Nevertheless, the causal relationship of SE on CRC remains unknown.

Methods

A two-sample Mendelian randomization study was performed with the summary statistics of sterol ester (27:1/14:0) from the largest available genome-wide association study meta-analysis(n = 377277) conducted by FinnGen consortium. The summary data were obtained from UK Biobank repository (377673 cases and 372016 controls). And we used relative filter (p < 5 x 10^− 6 and LD r² < 0.01) of instrumental variables to explore the causal effect and complete the sensitive analysis with the threshold p < 5 x 10^− 8 and LD r² < 0.01. Inverse variance weighted, MR-Egger, weighted median, Simple Mode and weighted model, were used to examine the causal association between SE (27:1/14:0) and CRC. Cochran’s Q statistics were used to quantify the heterogeneity of instrumental variables.

Results

The IVW results showed that SE (27:1/14:0) (OR = 1.004; 95% CI 1.002, 1.005; p < 0.001) have genetic causal relationship with CRC. The results of Weighted median, Weighted mode and Simple mode are all consistent with IVW models. Though, the result from the MR-Egger method (OR = 1.005; 95% CI 1.004, 1.009; p = 0.052) didn’t demonstrate a significant result. There was no heterogeneity, horizontal pleiotropy or outliers, and results were normally distributed. The MR analysis results were not driven by a single SNP. And results from two filter threshold is consistent.

Conclusion

Altogether, genetically predicted sterol ester (27:1/14:0) plays a causal association role in the incidence of CRC. This finding will provide a new screening and diagnosis indicator of CRC in the future.

sterol ester

causal

genetic

colorectal cancer

mendelian randomization

Colorectal cancer (CRC) is a highly prevalent malignant tumor of the digestive tract, ranking third in incidence and becoming the leading cause of cancer-related death in men and the second leading cause in women recently[1]. In 2020, the annual number of new colorectal cancer cases reached 1.9 million, and presently, the incidence rate is even higher, significantly contributing to the global cancer burden[2]. Typically, significant delays occur from the onset of symptoms to the diagnosis of early-onset colorectal cancer cause the asymptomatic at the beginning[3]. Therefore, it is crucial to prevent the occurrence of colorectal cancer through targeted health interventions aimed at addressing its etiology.

Colorectal cancer (CRC) is a multifactorial disease influenced by various risk factors, such as excess ethanol intake in alcoholic drinks[4], smoking[5], overweight[6], bacterial species[7] in gut including Fusobacterium nucleatum, enterotoxigenic Bacteroides fragilis, and pks + E. coli etc. Furthermore, certain lipid metabolism pathways are linked to an elevated risk of CRC[8]. And The disordered levels of lipids in colorectal cancer tissue and patient serum may be correlated with initiation of colorectal cancer[9]. Additionally, lipids such as phosphatidylcholine, phosphatidylethanolamine, and sphingomyelin may serve as diagnostic markers for various stages of colorectal cancer[10]. Nevertheless, the association between lipids and CRC warrants further investigation and validation in a larger cohort. Sterol esters(SE)[11], is a subtype of cholesteryl ester (CE), are a group of lipids implicated in the onset and progression of CRC. Liu et al. examined the expression levels of enzymes associated with cholesteryl ester metabolism and found a strong correlation between the expression of enzymes involved in cholesteryl ester synthesis (ACAT1 and ACAT2) and the occurrence of CRC[12]. However, Sadek et al[13]. discovered that free-form plant sterol esters have been shown to inhibit colon cancer via suppressing inflammation and inducing apoptosis. Given the contentious association between sterol esters and CRC, additional investigations are warranted to elucidate this relationship.

Mendelian randomization (MR) serves as an analytical method for making causal inferences within the realm of epidemiological etiology[14]. To unveil the causal relationship between the CRC and sterol ester (27:1/14:0), we conducted a two-sample Mendelian randomization utilizing methods including "MR-Egger," "Weighted median," "Inverse variance weighted," "Simple mode," and "Weighted mode." Genome-wide association studies (GWAS) have pinpointed numerous variants linked to diseases and traits located within noncoding regions of the genome.[15] By incorporating instrumental variables as genetic predictors, the causal relationship between genes and diseases remains unaffected by common confounding factors such as environmental influences, socioeconomic variables, and individual behaviors. Thus, this study is to investigate the causal relationship between the CRC and sterol ester (27:1/14:0) and provide some clinical helpful insights for the diagnosis of CRC patients.

2.1 Exposures: Ottensmann et al. analysis

We utilized data from Ottensmann et al.'s genome-wide association analysis of plasma lipidome, which identified 495 genetic associations with 179 lipids[16]. Ottensmann et al reported 11318730 sterol ester (27:1/14:0) associated genome-wide SNPs from 377277 FinnGen participants. Initially, we applied a threshold (p < 5 x 10^− 8) to identify independent SNPs. Subsequently, we clustered the SNPs within a 5000-kb region based on linkage disequilibrium (LD) with an r² > 0.01. Due to the limited number of SNPs, we set a relatively relaxed threshold (p < 5 x 10^− 6) and clumped LD r² > 0.01 in the 5000-kb region to obtain enough IVs to do sensitivity analysis.

2.2 Outcomes

Summary statistics of CRC were obtained from a GWAS comprising 377673 cases and 372016 controls of European population[17]. Our exposure and outcome population samples did not overlap.

We use the significant independent sterol ester (27:1/14:0) associated SNPs to assess the causal effect between CRC and sterol ester (27:1/14:0). All of the SNPs are available in the Summary data of CRC. The SNPs were used as IVs to perform the MR analysis (Supplementary Table 1). Ethical approval was obtained for each of the original GWAS, and detailed information can be found in the respective publications.

2.3 MR analysis

To investigate the causal relationship via the MR analysis, the IVs should satisfied the three assumptions[18]: (1) the selected IVs should be directly associated with sterol ester (27:1/14:0); (2) the selected IVs should not be associated with confounders;(3) the selected IVs must exert no effects on the CRC other than through the sterol ester (27:1/14:0) (Instrument Strength Independent of Direct Effect (INSIDE) assumption).

We utilized the classic Inverse Variance Weighted model (IVW) to perform the primary MR analysis, known for providing a stable and accurate assessment of causal relationships in the absence of directional pleiotropy[19, 20]. Furthermore, MR-Egger regression, Weighted Median (WM), Simple Mode, and Weighted Mode were employed as complementary methods to evaluate the causal effect. Under the INSIDE assumption, the MR-Egger method can estimate the horizontal average pleiotropic effect by conducting a weighted linear regression that relies on instrument strength independent of direct effects[18, 21]. The WM method plays a crucial role when the majority of genetic variants are invalid. It allows for obtaining a robust overall causal estimate, ensuring reliability even in the presence of invalid variants[22]. The Simple Mode and Weighted Mode are both mode-based methods, capable of estimating the causal effect of individual SNPs to form clusters. To elaborate, the Simple Mode selects the causal estimation from the largest cluster of SNPs, while the Weighted Mode assigns weights to each SNP[23].

The TwoSampleMR package (version 0.5.10) in R (version 4.2.3) was utilized to conduct the five MR methods.

2.4 Pleiotropy, heterogeneity, and sensitivity evaluation

The results from all five methods are required same directions, with p-values less than 0.05 indicating significance. We used the MR-Egger method to return intercept values for testing horizontal pleiotropy. Cochran’s Q statistic from IVW was taken as the assessment of heterogeneity. We also performed Leave-one-out and single SNP analyses to determine if a single SNP was driving the causal estimates. Additionally, another MR analysis was conducted with a stricter threshold of p < 5 x 10^− 8) and LD with r² < 0.01 for the complement of sensitivity analysis (supplementary materials).

After filtering of the SNPs from the Finngen at a genome-wide significance threshold of p < 5 x 10^− 6 and LD r² < 0.01, 14 SNPs associated with sterol ester (27:1/14:0) were identified (supplementary table 1). We conducted random-effects IVW models, along with two-sample MR analyses, utilizing these SNPs as instruments. The results indicate a genetically predicted causal effect of sterol ester on the risk of CRC (OR = 1.004; 95% CI 1.002, 1.005; p < 0.001) (Fig. 1A-1C). Though, the result from the MR-Egger method (OR = 1.005; 95% CI 1.004, 1.009; p = 0.052) didn’t demonstrate a significant result (supplementary table 2).

The results of Weighted median, Weighted mode, and Simple mode are all consistent with IVW models (Table). There was no heterogeneity observed among the selected SNPs, as demonstrated by both the MR-Egger (Cochran’s Q = 9.951; p = 0.620) and IVW methods (Cochran’s Q = 10.089; P = 0.687). Additionally, MR-Egger analysis was conducted to assess horizontal pleiotropy among the selected SNPs, revealing no evidence of pleiotropy (β intercept = -8.63E-05; SE = 0.000233; p = 0.717) influencing the results (Figure). The leave-one-out sensitivity analysis revealed that the causal association between sterol ester (27:1/14:0) and CRC is not driven by a single SNP (Fig. 1D). The associations of each variant with sterol ester and the risk of CRC are depicted in Figure S3.

In this study, a two-sample MR analysis was conducted using European GWAS data to investigate the relationship between CRC and sterol ester (27:1/14:0). Our findings indicate that sterol ester (27:1/14:0) exhibits a causal effect on CRC risk. To validate the conclusion, an MR analysis was also conducted using the 2 instrumental variables (IVs) meeting the threshold criteria (p < 5 x 10^− 8 and LD r² < 0.01). We employed the IVW method and obtained a similar result (OR = 1.007; 95% CI 1.002, 1.11; p = 0.007), reinforcing the same conclusion (supplementary materials).

Sterol ester (27:1/14:0), as a subtype of cholesterol esters, has been demonstrated to be associated with CRC[24, 25]. In a study conducted by Bershteĭn et al[24]. on CRC tissue, they observed a significant increase in cholesterol esters among CRC patients. Additionally, Munir et al[25]. compared isogenic primary and metastatic colon cancer cell lines and observed that, compared to SW480, CE content is significantly decreased in SW620 cells. This suggests that CE may play an important role in CRC progression. Enzymes involved in CE metabolism are strongly correlated with the pathogenesis, progression, and heterogeneity of CRC[12]. However, the precise causal relationship between sterol ester (27:1/14:0) and CRC, which could be beneficial for CRC diagnosis, remains unclear.

There is some experimental evidence for the association between sterol ester and CRC. In a previous study, adenomatous polyposis coli mice were fed with plant sterol ester, revealing that a high intake of plant sterol ester accelerates intestinal tumorigenesis in females. Gene expression analysis from the mucosa of the small intestine illustrated up-regulated genes associated with cell cycle control and cholesterol biosynthesis in female mice fed with plant sterol ester[26]. Sterol ester (27:1/14:0) may influence the cell cycle of epithelial cells, potentially leading to the development of CRC.

This study explores the cause relationship between sterol ester (27:1/14:0) and CRC without via MR methods. All the methods indicate a significant result except MR-Egger method while MR-Egger analysis revealed no evidence of pleiotropy (β intercept = -8.63E-05; SE = 0.000233; p = 0.717). Thus, the non-significant from MR-Egger may due to the limited number of samples. In addition, all the SNPs were filter at the threshold of p < 5 x 10^− 6 and LD r² < 0.01, indicating a lower likelihood of weak instrument bias. And the MR result from the threshold of p < 5 x 10^− 8 and LD r² < 0.01 is consist with the previous one which is further decrease potential the bias. What’s more, we used several approaches to assess the sensitive and pleiotropy, indicating the similar result.

Nonetheless, this study has some limitations that need to be acknowledged. Firstly, all SNP data are from European cohorts, and further validation of our findings is needed in other ethnicities. Additionally, we employed a relatively relaxed threshold (p < 5 x 10^− 6 and LD r² < 0.01), as only 2 SNPs were identified under the more stringent criteria of p < 5 x 10^− 8 and LD r² < 0.01. Furthermore, the experimental validation is warrant for the causal relationship between sterol ester (27:1/14:0).

Altogether, the association between genetically predicted sterol ester (27:1/14:0) and CRC has been established. This finding will contribute to the screening and diagnosis of CRC in the future.

SE, sterol ester; CE, cholesterol esters; CRC, Colorectal cancer, LD, linkage disequilibrium, IVW, Variance Weighted model; WM, Weighted Median; IVs, instrumental variables.

Ethics approval and consent to participate

Not applicable

Consent for publication

Not applicable

Availability of data and materials

The datasets analyzed during the current study are available in the FinnGen repository (https://ftp.ebi.ac.uk/pub/databases/gwas/summary_statistics/GCST90277001-GCST90278000/GCST90277238/) and the UK Biobank repository, https://gwas.mrcieu.ac.uk/datasets/ieu-b-4965/.

Competing interests

The authors hereby declare that they have no conflicts of interest to disclose.

Funding

None

Authors' contributions

CYL: Methodology, data curation, formal analysis, visualization, original draft writing. JFX, BLY, JQZ: revised the manuscript, and XX: design of the work. All authors reviewed the manuscript.

Acknowledgements

We would like to express our gratitude to the participants and investigators involved in the FinnGen study, as well as the collaborators of the UK Biobank.

Siegel RL, Giaquinto AN, Jemal A. Cancer Stat 2024 CA: cancer J Clin. 2024;74(1):12–49.
Sung H, et al. Global Cancer Statistics 2020: GLOBOCAN Estimates of Incidence and Mortality Worldwide for 36 Cancers in 185 Countries. Cancer J Clin. 2021;71(3):209–49.
Patel SG, et al. The rising tide of early-onset colorectal cancer: a comprehensive review of epidemiology, clinical features, biology, risk factors, prevention, and early detection. The lancet. Gastroenterol Hepatol. 2022;7(3):262–74.
Vieira AR, et al. Foods and beverages and colorectal cancer risk: a systematic review and meta-analysis of cohort studies, an update of the evidence of the WCRF-AICR Continuous Update Project. Annals oncology: official J Eur Soc Med Oncol. 2017;28(8):1788–802.
Botteri E, et al. Smoking and colorectal cancer: a meta-analysis. JAMA. 2008;300(23):2765–78.
Kyrgiou M et al. Adiposity and cancer at major anatomical sites: umbrella review of the literature. BMJ (Clinical research ed.), 2017. 356: p. j477.
Dougherty MW, Jobin C. Intestinal bacteria and colorectal cancer: etiology and treatment. Gut Microbes. 2023;15(1):2185028.
O'Keefe SJD. Diet, microorganisms and their metabolites, and colon cancer. Nat Rev Gastroenterol Hepatol. 2016;13(12):691–706.
Zhang X, et al. Lipid levels in serum and cancerous tissues of colorectal cancer patients. World J Gastroenterol. 2014;20(26):8646–52.
Elmallah MIY, et al. Lipidomic profiling of exosomes from colorectal cancer cells and patients reveals potential biomarkers. Mol Oncol. 2022;16(14):2710–8.
Lampe MA, et al. Human stratum corneum lipids: characterization and regional variations. J Lipid Res. 1983;24(2):120–30.
Liu Z, et al. Correlation of cholesteryl ester metabolism to pathogenesis, progression and disparities in colorectal Cancer. Lipids Health Dis. 2022;21(1):22.
Sadek NF, et al. Plant Sterol Esters in Extruded Food Model Inhibits Colon Carcinogenesis by Suppressing Inflammation and Stimulating Apoptosis. J Med Food. 2017;20(7):659–66.
Smith GD, Ebrahim S. 'Mendelian randomization': can genetic epidemiology contribute to understanding environmental determinants of disease? 2003. p. 1–22.
Farh KK, et al. Genetic and epigenetic fine mapping of causal autoimmune disease variants. Nature. 2015;518(7539):337–43.
Ottensmann L, et al. Genome-wide association analysis of plasma lipidome identifies 495 genetic associations. Nat Commun. 2023;14(1):6934.
Tassi A, Mavromatis I, Piechocki RJR. A dataset of full-stack ITS-G5 DSRC communications over licensed and unlicensed bands using a large-scale urban testbed. Data brief. 2019;25:104368.
Bowden J, Davey Smith G, Burgess S. Mendelian randomization with invalid instruments: effect estimation and bias detection through Egger regression. Int J Epidemiol. 2015;44(2):512–25.
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(11):1880–906.
Burgess S, et al. Using published data in Mendelian randomization: a blueprint for efficient identification of causal risk factors. Eur J Epidemiol. 2015;30(7):543–52.
Burgess S, Thompson SG. Interpreting findings from Mendelian randomization using the MR-Egger method. Eur J Epidemiol. 2017;32(5):377–89.
Bowden J, et al. Consistent Estimation in Mendelian Randomization with Some Invalid Instruments Using a Weighted Median Estimator. Genet Epidemiol. 2016;40(4):304–14.
Hartwig FP, Davey Smith G, Bowden J. Robust inference in summary data Mendelian randomization via the zero modal pleiotropy assumption. Int J Epidemiol. 2017;46(6):1985–98.
Bershteĭn LM, et al. [Content of cyclic nucleotides, cholesterol and phospholipids in tumors of the large intestine]. Vopr Onkol. 1985;31(2):50–4.
Munir R, et al. Abundance, fatty acid composition and saturation index of neutral lipids in colorectal cancer cell lines. Acta Biochim Pol. 2021;68(1):115–8.
Marttinen M, et al. Plant sterol feeding induces tumor formation and alters sterol metabolism in the intestine of Apc(Min) mice. Nutr Cancer. 2014;66(2):259–69.

No competing interests reported.

Download PDF

Editorial decision: Revision requested
19 Aug, 2024
Reviews received at journal
18 Aug, 2024
Reviewers agreed at journal
08 Aug, 2024
Reviews received at journal
25 Jul, 2024
Reviewers agreed at journal
21 Jul, 2024
Reviewers invited by journal
13 May, 2024
Editor invited by journal
13 May, 2024
Editor assigned by journal
10 May, 2024
Submission checks completed at journal
10 May, 2024
First submitted to journal
04 May, 2024

You are reading this latest preprint version

Unraveling the link between sterol ester and colorectal cancer: a two-sample Mendelian randomization study

Status:

Version 1

Abstract

Background

Methods

Results

Conclusion

Figures

1. INTRODUCTION

2. METHODS

2.1 Exposures: Ottensmann et al. analysis

2.2 Outcomes

2.3 MR analysis

2.4 Pleiotropy, heterogeneity, and sensitivity evaluation

3. Results

4. Disscussion

Abbreviations

Declarations

References

Additional Declarations

Supplementary Files

Status:

Version 1