- 1Department of Oncology, First Teaching Hospital of Tianjin University of Traditional Chinese Medicine, Tianjin, China
- 2National Clinical Research Center for Chinese Medicine Acupuncture and Moxibustion, Tianjin, China
Background and Objective: Enhancing therapy choices for varying stages of esophageal cancer and improving patient survival depend on timely and precise diagnosis. Blood metabolites may play a role in either causing or preventing esophageal cancer, but further research is needed to determine whether blood metabolites constitute a genetic risk factor for the disease. In order to tackle these problems, we evaluated the causal association between esophageal cancer and 486 blood metabolites that functioned as genetic proxies using a two-sample Mendelian randomization (MR) study.
Methods: We utilized two-sample MR analyses to evaluate the causal links between blood metabolites and esophageal cancer. For the exposure, we used a genome-wide association study (GWAS) of 486 metabolites, and a GWAS study on esophageal cancer from Sakaue et al. was used for preliminary analyses. Causal analyses employed randomized inverse variance weighted (IVW) as the main method, supplemented by MR-Egger and weighted median (WM) analyses. Sensitivity analyses included the MR-Egger intercept test, Cochran Q test, MR-PRESSO, and leave-one-out analysis. Additionally, independent esophageal cancer GWAS data were utilized for replication and meta-analysis. FDR correction was applied to discern features with causal relationships. For conclusive metabolite identification, we conducted Steiger tests, linkage disequilibrium score regression, and colocalization analyses. Moreover, we utilized the program MetaboAnalyst 5.0 to analyze metabolic pathways.
Results: This study found an important association between esophageal cancer and three metabolites: 1-linoleoylglycerophosphoethanolamine* [odds ratio (OR) = 3.21, 95% confidence interval (CI): 1.42–7.26, p < 0.01], pyroglutamine* (OR = 1.92, 95% CI: 1.17–3.17, p < 0.01), and laurate (12:0) (OR = 3.06, 95% CI: 1.38–6.78, p < 0.01).
Conclusion: This study establishes a causal link between three defined blood metabolites and esophageal cancer, offering fresh insights into its pathogenesis.
1 Introduction
Esophageal cancer seventh in global cancer incidence and is the sixth leading cause of cancer-related deaths (Global Burden of Disease, 2019 Cancer Collaboration et al., 2022). Esophageal cancer develops insidiously, with early symptoms often unnoticed by patients. Diagnosis commonly occurs at an advanced stage, resulting in a poor prognosis (Global Burden of Disease, 2019 Cancer Collaboration et al., 2022). Globally, the incidence of esophageal cancer exhibits distinct geopolitical patterns. Typically, esophageal adenocarcinoma is more prevalent in Western countries, whereas esophageal squamous cell carcinoma (ESCC) is notably more common in Asia’s esophageal cancer belt, as well as in certain areas of East Africa and South America (Shaheen and Ransohoff, 2002). In 2020, an estimated 6,041,000 new cases and 5,440,762 deaths due to esophageal cancer were reported globally (Sung et al., 2021). The heterogeneity and complexity of esophageal cancer lead to most patients being diagnosed at an advanced stage, precluding the possibility of surgery and optimal treatment. Consequently, patients with esophageal cancer have a worse prognosis, increased mortality, and a worse quality of life, evidenced by 5-year survival rates under 20% (Zhang et al., 2023), thereby intensifying the global disease burden.
The incorporation of metabolomics with systems biology has recently offered a new approach in researching disease processes. It remains unclear whether metabolic dysregulation is a cause or a result of esophageal cancer progression. However, substantial evidence indicates that it plays a critical role in both the development and progression of the disease (Wang et al., 2013), beyond changes in glucose metabolism, dysregulation in amino acid, lipid, and protein metabolism has been observed in vitro and in vivo, as exemplified by the well-known Warburg effect (Gu et al., 2016; Kim et al., 2010; Matsunaga et al., 2018). Although Huang et al. (2020) described 45 blood metabolites associated with esophageal cancer, no thorough investigation has systematically examined the causal connection between esophageal cancer and blood metabolites. The possible causal relationship between alterations in blood metabolite levels and the risk of esophageal cancer requires more investigation. This will aid in understanding the causes of esophageal cancer and enhance its precise prevention and control.
The Mendelian randomization (MR) study is a well-established genetic epidemiological methodology that employs genetic variants as proxies for target exposures to assess the causal relationship between exposures and disease outcomes (Zuccolo and Holmes, 2017). This approach has distinct advantages over traditional observational studies. First, by utilizing random independent classification of DNA during meiosis of alleles, MR analysis avoids traditional confounding. Second, reverse causal effects bias does not occur due to the immutability of the genomes of the propagating germline (Burgess et al., 2015). Third, in most cases, genetic variation is usually accurately measured and reported and is less susceptible to bias and measurement error. Thus, it is particularly useful in assessing risk factors for long-term effects (Hu et al., 2019). In this context, we performed MR analysis with genome-wide association study (GWAS) pooled data to comprehensively explore the causal effects of 486 blood metabolites on esophageal cancer. Additionally, colocalization and metabolic pathway analyses were conducted to explore the gene and protein-level biological processes underlying esophageal cancer. The objective of this study is to uncover the metabolic etiology of esophageal cancer and offer insights into its biological mechanisms.
2 Methods
2.1 Study design
Three presumptions must be met by a valid MR study: (1) IVs are highly correlated with the relevant exposures; (2) IVs are unaffected by confounding variables; (3) IVs have no correlation to outcomes and only have an impact on outcomes through exposures (Boef et al., 2015). All of the MR analyses in this work were performed using R software (version 4.3.2) Two Sample MR and MRPRESSO packages. An overview of the study is shown in Figure 1.
Figure 1. Overview of the study design. Assumption 1, genetic instruments are strongly associated with the exposures of interest; Assumption 2, genetic instruments are independent of confounding factors; Assumption 3, genetic instruments are not associated with outcome and affect outcome only via exposures. LD, linkage disequilibrium; SNPs, single nucleotide polymorphisms; IVW, inverse variance weighted; WM, weighted median; LOO analysis, leave-one-out analysis; MR-PRESSO, MR-Pleiotropy RESidual sum and outlier.
2.2 GWAS summary statistics
All analyses of MR in this study were based on summary statistics from the GWAS study of blood metabolites and Esophageal cancer. The summary data for Esophageal cancer is obtained from Sakaue et al. (2021)’s trans-ethnic GWAS meta-analysis, which includes 2,386 cases (998 cases of European origin and 1,388 cases of East Asian origin) and 476,696 controls (475,308 controls of European origin and 159,201 controls of East Asian origin). From the Metabolomics GWAS Server (https://metabolomics.helmholtz-muenchen.de/gwas/), genetic information on blood metabolites was collected. This is the largest analysis on blood metabolite genetic loci to date, with approximately 2.1 million single nucleotide polymorphisms (SNPs) for 486 metabolites associated with genetic variations in humans. Shin et al. (2014)’s high-throughput metabolic profiling and genome-wide association searches were used to accomplish this. Supplementary Table S1 lists the 486 metabolites, with those labeled ‘X-’ indicating unknown chemical properties. The study encompassed 7,824 European participants, comprising 1,768 from Germany’s KORA F4 study and 6,056 from the UK Twin Study. 107 of the 486 metabolites had their chemical characteristics so poorly characterized that they were categorized as unknown. The Kyoto Encyclopedia of Genes and Genomes (KEGG) database also shows that 309 metabolites were chemically confirmed and classified into eight metabolic groups: lipid, nucleotide, energy, peptide, xenobiotic metabolism, carbohydrate, cofactors and vitamins, and amino acid (Kanehisa et al., 2012). As this study was based on publicly available data, no additional ethical approval or consent was required.
2.3 IVs selection
MR Research has three core instrumental variable assumptions. (1) Selection of genetic variants as instrumental variables (IVs) associated with exposure factors blood metabolites. To find SNPs associated with exposure factors and to ensure the veracity and accuracy of conclusions about the causal link between blood metabolites and Esophageal cancer risk, the following steps were used to select the best SNPs. After consulting the pertinent literature, the threshold was first modified (p < 1 × 10−5) due to the small number of SNPs related with metabolites in order to achieve comprehensive and trustworthy results. Then, we eliminated linkage disequilibrium (LD, R2 < 0.001, kb = 10,000) to clumped SNPs (Yang et al., 2020; Choi et al., 2019). The F-statistic was used to assess the correlation between SNPs and exposure statistical strength, the F-statistic is given by: F =
2.4 Statistical analysis and sensitivity analysis
The primary analysis utilized the random-effects inverse variance weighted (IVW) method to establish a significant causal link between serum metabolites and esophageal cancer (p < 0.05). When IVs meet all three main hypotheses, the IVW method yields the most accurate and unbiased estimates, making it the most effective method for MR (Pierce and Burgess, 2013). However, horizontal pleiotropy arises when IVs impact traits beyond the exposure pathway, directly affecting the outcomes and potentially leading to imprecise causal estimates (Thomas and Conti, 2004). MR-Egger regression provides consistent estimates in the presence of genetically pleiotropic IVs (Bowden et al., 2015). The WM method’s strength lies in its ability to yield consistent causal estimates even with over 50% null IVs (Bowden et al., 2016). In repeated comparisons, we used the Benjamini–Hochberg technique to control the false discovery rate (FDR) (Benjamini and Hochberg, 1995). The Benjamini-Hochberg-corrected p values (FDR p-values) were used as the threshold of significance. Associations that passed the FDR controlled significance levels (FDR p-values < 0.05) were considered to be strong evidence of associations, whereas results with p less than 0.05 but that failed to pass FDR correction were regarded as suggestive associations.
Sensitivity analysis is essential for detecting horizontal pleiotropy and heterogeneity in the data. Cochran’s Q statistic was employed to assess heterogeneity (Cohen et al., 2015). The MR-Pleiotropy RESidual Sum and Outlier (MR-PRESSO) method was used to identify outliers, which were then excluded for re-analysis if detected (Verbanck et al., 2018). Horizontal pleiotropy was assessed using MR-Egger, a statistically significant difference in intercept terms indicated notable horizontal pleiotropy in the study (Burgess and Thompson, 2017). Robustness of the MR findings was evaluated using used leave-one-out (LOO) analysis. By eliminating each SNP individually, after which MR analysis is conducted to see whether a single SNP significantly affected the outcome.
In conclusion, multiple criteria were used in our thorough screening for blood metabolites that may be causally linked to esophageal cancer: (1) Significant p-values (IVW-derived p < 0.05, FDR p < 0.05); (2) uniformity in magnitude and direction among the 3 MR methods; (3) No horizontal pleiotropy or heterogeneity in the MR results; (4) Lack of significant confounding in MR estimates by individual SNPs.
2.5 Replication and meta-analysis
To thoroughly evaluate the robustness of the identified candidate metabolites, we utilized data from the FinnGen consortium (https://www.finngen.fi/fi) for additional validation of the exposure results and their relation to esophageal cancer. FinnGen’s release R10 includes data on 169 esophageal cancer cases and 314,193 other cases. The causal association of blood metabolites with esophageal cancer was confirmed through meta-analysis of 2 MR studies. Using R software (4.3.2), the random effects IVW model was used to implement meta-analysis.
2.6 Confounding analysis
While we evaluated the horizontal pleiotropy of MR results through sensitivity analyses to identify SNPs violating MR assumptions, a small residual amount of confounding SNPs may remain. IVs were examined using the Phenoscanner V2 website (http://www.phenoscanner.medschl.cam.ac.uk/) (Staley et al., 2016). Each SNP was evaluated for its association with known esophageal cancer risk factors, including smoking (Doll et al., 1994; Freedman et al., 2007), alcohol consumption (Yokoyama et al., 2002; Prabhu et al., 2014), and obesity (Veugelers et al., 2006; Lagergren, 2011). If any SNPs were associated with these confounders (p < 1 × 10−5), we reran the MR analysis after their removal to confirm the results’ reliability.
2.7 Evaluation of directionality
The Steiger test was used to confirm if the observed causality was influenced by reverse causality. This test examined whether the SNPs included in the study explained more of the esophageal cancer than the metabolites identified (Hemani et al., 2017). When SNP combinations were shown to have no genetic risk for esophageal cancer when compared to metabolites, the results showed that causal inference was not biased (Steiger p > 0.05).
2.8 Colocalization analysis
After MR analyses identified causality of blood metabolites on esophageal cancer risk, we used co-localization analyses as a sensitivity analysis to s to evaluate whether blood metabolites and esophageal cancer risk share the same causal genetic variant at the GHRL locus (Giambartolomei et al., 2014). Such an analysis can indicate whether the phenotypes are influenced by different causal genetic variants that are in LD, indicative of horizontal pleiotropy and violation of the exclusion restriction assumption (Wallace, 2013). PP.H4 > 75% was considered supportive of colocalization of the two phenotypes.
2.9 Metabolic pathway analysis
Using MetaboAnalyst 5.0 (https://www.metaboanalyst.ca/), we performed a metabolic pathway analysis (Jewison et al., 2014) to elucidate the biological mechanisms through which blood metabolites causally affect esophageal cancer, thereby exploring its potential pathogenesis.
3 Results
3.1 Preliminary analysis
Following rigorous quality control of the IVs, our MR study identified 486 blood metabolites. Preliminary research depicted in Figure 2 points to a relationship between blood metabolites and the risk of esophageal cancer. The refined IVs comprised 276 SNPs. IVW analysis identified 14 metabolites potentially causally related to esophageal cancer, consisting of 9 metabolites with known chemical properties and 5 with unknown properties, as shown in Figure 3. Figure 3 illustrates that the nine identified metabolites are chemically classified into amino acids, lipids, and xenobiotics. After FDR, no metabolites were found that still had a significant effect. Palmitate (16:0) (odds ratio [OR] = = 0.19, 95% confidence interval [CI]: 0.07–0.57, p = 0.003, adjust p = 0.013), X-06226 (OR = 0.29, 95%CI: 0.11–0.79, p = 0.015, adjust p = 0.03); phenol sulfate (OR = 0.52, 95%CI: 0.29–0.95, p = 0.033, adjust p = 0.038); N-acetylglycine (OR = 0.52, 95%CI: 0.34–0.81, p = 0.004, adjust p = 0.014); homostachydrine* (OR = 0.38, 95%CI: 0.17–0.86, p = 0.020, adjust p = 0.035); X-12740 (OR = 0.73, 95%C: 0.55–0.9, p = 0.028, adjust p = 0.039) was negatively associated with the risk of esophageal cancer, whereas laurate (12:0) (OR = 4.03, 95%CI: 1.65–9.81, p = 0.002, adjust p = 0.03); X-05426 (OR = 2.11, 95%CI: 1.07–4.16, p = 0.031, adjust p = 0.039); X-11261 (OR = 1.99, 95%CI: 1.00–3.96 p = 0.049, adjust p = 0.049); 1-linoleoylglycerophosphoethanolamine*9 (OR = 4.36, 95%CI: 1.68–11.29, p = 0.002, adjust p = 0.017); pyroglutamine* (OR = 2.05, 95%CI: 1.18–3.56, p = 0.011, adjust p = 0.031); isobutyrylcarnitine (OR = 2.68, 95%CI: 1.13–6.36, p = 0.026, adjust p = 0.04); X-1349(OR = 5.01, 95%CI: 1.39–18.09, p = 0.014, adjust p = 0.033); 1-arachidonoylglycerophosphoinositol* (OR = 4.36, 95%CI: 1.68–11.29, p = 0.043, adjust p = 0.046) were positively correlated with the risk of esophageal cancer. For 1-arachidonoylglycerophosphoinositol* (OR = 4.36, 95%CI: 1.68–11.29, p = 0.043, β = 0.919), the MR-Egger method showed different results (OR = 0.743, 95%CI: 0.108–5.129, p = 0.767, β = −0.297). Despite no horizontal pleiotropy in MR-Egger regression (p_intercept = 0.184), potential outliers may cause inconsistencies between MR-Egger and IVW results. Thus, further investigations are needed to validate its association with esophageal cancer. In summary, IVW-derived estimates showed significance (p < 0.05), and the direction and magnitude of the IVW, MR-Egger, and WM estimates aligned consistently, as depicted in Figure 4. The Cochran Q test indicated no heterogeneity among the selected SNPs, as evidenced by Q_pval values for IVW and MR-Egger all exceeding 0.05. Homostachydrine* was excluded by the MR-PRESSO method and the rest were not detected with outliers (p > 0.05). MR-PRESSO results did not indicate heterogeneous SNPs post-outlier removal (Supplementary Table S3). Cochran Q test and MR-Egger cutoff test results (both p > 0.05) strongly suggest the lack of pleiotropy and heterogeneity (Table 1). The results of the LOO analysis indicated that no single SNP introduced bias into the MR estimates (Supplementary Figure S1). Twelve blood metabolites were identified as subjects for more investigation.
Figure 2. Preliminary MR estimates of the association between blood metabolites and esophageal cancer risk. The shades of color depict the magnitude of the p-value.
Figure 3. Forest plot based on inverse variance weighted (IVW) analysis was created to show the relationship between blood metabolites and esophageal cancer. The graph shows the correlation between a number of blood metabolites and the risk of esophageal cancer. Each horizontal line represents the OR and 95% CI for how each blood metabolite affects the risk of esophageal cancer. IVW, inverse variance weighted; CI, confidence interval; SNPs, single nucleotide polymorphisms; OR, odds ratio.
Figure 4. Scatterplot shows estimates that are strongly linked (IVW determined p < 0.05) and directionally consistent. SNP, single nucleotide polymorphisms.
Table 1. Supplementary and sensitivity analysis examine the relationship between blood metabolites and esophageal cancer.
3.2 Replication and meta-analysis
To bolster the credibility of our findings, we replicated the MR analysis using a separate GWAS dataset for esophageal cancer. The other GWAS dataset revealed a consistent trend for the candidate metabolites, albeit the results were not statistically significant, most likely due to major sample size differences. The meta-analysis identified 6 blood metabolites (3 known and 3 unknown) with potential effects on esophageal cancer, as shown in Figure 5. Higher levels of X-06226 (OR = 0.35, 95%CI: 0.14–0.89, p = 0.03) were associated with reduced esophageal cancer risk. Conversely, elevated levels of 1-linoleoylglycerophosphoethanolamine* (OR = 3.21, 95%CI: 1.42–7.26, p < 0.01), pyroglutamine* (OR = 1.92, 95%CI: 1.17–3.17, p < 0.01), laurate (12:0) (OR = 3.06, 95%CI: 1.38–6.78, p < 0.01), X-13496 (OR = 4.36, 95%CI: 1.43–13.25, p < 0.01), and X-11261 (OR = 1.84, 95%CI: 1.03–3.26, p = 0.04) appeared to increase the risk. Although palmitate (16:0), isobutyrylcarnitine, N-acetylglycine, phenol sulfate, X-05426, and X-12740 consistently aligned in both MR analyses, they were excluded due to non-significant estimates in the meta-analysis. Although the six metabolites mentioned above were not significant in the replication analysis, their consistent directionality suggests more than mere coincidence. This pattern in MR estimation might be due to differences in sample size. Therefore, by expanding the sample size and statistical robustness through meta-analysis, we identified metabolites causally linked to esophageal cancer.
Figure 5. Meta-analysis of significantly associated (IVW derived p < 0.05) between blood metabolites and esophageal cancer. OR, odds ratio; 95% CI, 95%confidence interval. (A) X-06226, (B) 1-linoleoylglycerophosphoethanolamine*, (C) pyroglutamine*, (D) laurate (12:0), (E) X-13496, (F) X-11261.
3.3 Confounding analysis
While sensitivity analysis excluded SNPs violating our estimates, to meet assumption 2 (IVs are independent of confounders), we checked in Phenoscanner if all SNPs linked to the six metabolites were unaffected by risk factors for esophageal cancer, including alcohol consumption, smoking, and obesity. Our analysis revealed no confounders for X-06226, pyroglutamine*, and X-1349. However, among other IVs related to three metabolites, five SNPs were associated with confounders (Supplementary Table S4). After excluding these SNPs and reapplying MR analysis, we determined that laurate (12:0) (OR = 4.02, 95%CI: 1.62–9.96, p = 0.003) and 1-linoleoylglycerophosphoethanolamine* (OR = 4.22, 95%CI: 1.73–10.34, p = 0.002) maintained a stable association with esophageal cancer. In contrast, X-11261 (OR = 1.67, 95%CI: 0.55–2.48, p = 0.690) showed an unstable association.
3.4 Evaluation of directionality
Furthermore, the Steger test was conducted to investigate the possibility of an inverse causal relationship between the three chemically identified metabolites and esophageal cancer. The Steger test results did not indicate a reverse causal effect between these metabolites and esophageal cancer (p < 0.05) (Supplementary Table S5).
3.5 Colocalization analysis
Colocalization analysis results revealed a low probability of shared genetic variation between the metabolites and esophageal cancer (laurate (12:0): PPH4 = 22%; pyroglutamine*: PPH4 = 3%; 1-linoleoylglycerophosphoethanolamine*: PPH4 = 21%) (Supplementary Table S6). This implies that the MR analysis is not bias by horizontal pleiotropy.
3.6 Metabolic pathway analysis
Utilizing three identified metabolites, we pinpointed three metabolic pathways potentially implicated in esophageal cancer pathogenesis (Supplementary Table S7). The pathways of Mitochondrial Beta-Oxidation of Medium Chain Saturated Fatty Acids, Beta Oxidation of Very Long Chain Fatty Acids, and Fatty Acid Biosynthesis are suggested as underlying biological mechanisms in the development of esophageal cancer. Significantly, laurate (12:0) is a component in all these metabolic pathways. This suggests a crucial role for laurate (12:0) and its associated metabolic pathways in the pathogenesis of esophageal cancer.
4 Discussion
In the current research, we merged two large-scale GWAS data to examine the causal effects of 486 blood metabolites on esophageal cancer by genetic proxy through a rigorous MR design. We confirmed causal associations between 3 known blood metabolites and esophageal cancer risk. Specifically high levels of 1-linoleoylglycerophosphoethanolamine*, pyroglutamine*, and laurate (12:0) were genetically predisposed to increase susceptibility to esophageal cancer. Three metabolic pathways that might have a role in the biological mechanisms behind esophageal cancer have been discovered. To the best of what we know, this is the first MR research to use colocalization analysis and metabolic pathways in conjunction with the largest available blood metabolite GWAS data to investigate the causal relationship with esophageal cancer. These results underline the critical role that blood metabolites play in the pathophysiology of esophageal cancer and offer insightful information for potential studies on the detection and avoidance of the disease.
The high morbidity and mortality rates of esophageal cancer have placed a tremendous burden on people globally in recent years, making early detection and earlier treatment of the illness an urgent goal. Esophageal cancer is associated with a variety of dietary habits. A meta-analysis showed that meat, red meat, saturated fat, salt intake, and the temperature of the food were positively associated with the incidence of esophageal cancer (Castro et al., 2018). An important part of preventing esophageal cancer is understanding the dietary exposure factors linked to the disease. In addition, animal, or saturated fats, including butter, have been positively associated with a high incidence of esophageal cancer (Tullio et al., 2020). However, the evidence to date is not sufficient to demonstrate exactly how they promote cancer development, and many studies are still needed to prove it. The perceived value of metabolites in esophageal tissue is being investigated with more attention due to the advent of metabolomics technology. Notably, because blood metabolites concurrently record endogenous and external activities, they visually offer a picture of biological operations (Rappaport et al., 2014). The contribution of metabolites to early screening and prevention of esophageal cancer is restricted because of the uncertain causal link between the two, despite prior research showing that they are advantageous to therapy and have a part in the molecular processes of esophageal cancer. In order to clarify the causative association between blood metabolites and esophageal cancer and to investigate the metabolic pathways involved, we thus carried out a critical MR study. This study served as a guide for future research aimed at fully understanding the pathophysiology of esophageal cancer.
The identification of biomarkers is a significant therapeutic contribution of this work. Our study supports the positive association between laurate (12:0), pyroglutamine*, 1-linoleoylglycerophosphoethanolamine* and the risk of esophageal cancer by combining genetics and metabolomics from a causality perspective. And laurate (12:0) is implicated in three significantly enhanced metabolic pathways: Mitochondrial Beta-Oxidation of Medium Chain Saturated Fatty Acids, Beta Oxidation of Very Long Chain Fatty Acids, and Fatty Acid Biosynthesis. Laurate (12:0), also known as dodecanoic acid, is one of the saturated fatty acids and a medium chain fatty acid (MCFA), which is the body’s energy supply fuel (Schönfeld and Wojtczak, 2016). Saturated fatty acids are a class of fatty acids (FA) that do not contain unsaturated double bonds and are one of the basic components that make up lipids (Schönfeld and Wojtczak, 2016). Saturated fatty acids are risk factors for several cancers. Laurate (12:0), one of the saturated fatty acids, can be ingested through the diet or produced through the digestion of medium and long chain fatty acids at the small intestine. These medium-chain fatty acids are absorbed through the intestinal wall and enter the blood circulatory system to reach all parts of the body to perform their physiological roles. Related studies have shown that lauric acid can stimulate mammary cell proliferation in mice by activating the GPR84 and PI3K/Akt signaling pathways (Meng et al., 2017). In in vitro studies, lauric acid can activate macrophages by promoting signaling at Toll-like receptors (TLR2 heterodimer and TLR4 homodimer) (Huang et al., 2012; Lee et al., 2004), and it can also play a pro-inflammatory role by activating macrophages via the MCFA receptor GPR84 (Suzuki et al., 2013).
Increased FA biosynthesis is one of the hallmarks of abnormal metabolism in tumor cells. Currie et al. (2013) and Cao (2019) have shown the significance of FAs production for cancer cell proliferation and survival. The increase in fatty acid biosynthesis may be a response to the high metabolic demands of cancer cells or an adaptation to the reduced availability of serum-derived lipids in the tumor microenvironment. Systemic mobilisation of lipids from adipose tissue during cancer cachexia promotes tumor growth (Argilés et al., 2014). Thus, the significant increase in laurate (12:0) in esophageal cancer patients may be an adaptation to the high energy expenditure of cancer patients. FA biosynthesis is considered a potential therapeutic target for cancer. A study showed that saturated free fatty acid palmitate induced EMT in hepatocellular carcinoma (HCC) cells through activation of the Wnt/catenin and TGF-AKT pathways (Nath et al., 2015). Furthermore, in studies related to breast cancer, free fatty acids (FFA) are circulating plasma factors that are associated with increased proliferation and invasiveness of estrogen receptor alpha (ERα)-positive (ER +) breast cancer cells. FFA activates ERα and mTOR pathways in breast cancer cells and alters metabolism (Madak-Erdogan et al., 2019). Therefore, we speculate that the development of esophageal cancer affects fatty acid synthesis in the body, which may contribute to cancer progression through increased uptake of various fatty acids such as lauric acid.
Though its function in esophageal cancer is yet unknown, 1-linoleoyglycerophosphoethanolamine* is substantially expressed in colorectal cancer tissues (Mika et al., 2020) and has been demonstrated in a MR experiment to be protective against colorectal cancer (Yun et al., 2023). Furthermore, 1-linoleoylglycerophosphoethanolamine is a significant component of the phosphatidy-lethanolamine (PE) (van der Veen et al., 2017). PE, a significant constituent of phospholipids in cell membranes, is linked to anxiety and plays a crucial role in preserving the stability of cell structure. According to Reichel’s study (Reichel et al., 2015), people with alcohol dependence frequently felt anxious when they abstained from drinking and had higher plasma PE concentrations when they did not. Furthermore, Yang et al. (2020)’s investigation into the connection between metabolites and certain mental conditions discovered that 1-linoleoylglycerophosphatidylethanolamine was linked to a higher risk of major depressive disorder; this finding was subsequently corroborated by another research (Xiao et al., 2022). This seems to suggest that 1-linoleoylglycerophosphatidylethanolamine may well be able to affect depression, which might therefore act as a mediator in the development of esophageal cancer. So far, there are no reports on Pyroglutamine and esophageal cancer, and there is only one study suggesting that Pyroglutamine may be associated with prostate carcinogenesis, but the mechanism is not clear, and more research is necessary to determine the connection between pyroglutamine and cancer.
Our study has several strengths. This is the first and most important point. This is the most thorough and organized investigation of the causative association between blood metabolites and esophageal cancer, and the first MR research to evaluate the relationship between blood metabolites and esophageal cancer. With this kind of design, the limits of confounders that are frequently present in conventional observational studies are mitigated, and there may be more proof of a causal association between exposure and result. Secondly, our results are compelling. In terms of direction and sensitivity analysis, all 3 MR estimations exhibit good consistency. Third, we conducted replication studies and meta-analyses using extra GWAS data and incorporated GWAS datasets from several populations in our study. As a result, our conclusions are now more trustworthy and thorough. Finally, our study provides new insights into the molecular pathways involved in esophageal carcinogenesis by combining genomics and metabolomics.
The present investigation is not without limits. First, our MR analysis set a little relaxed threshold due to the modest number of SNPs related with metabolites. Nonetheless, F-statistics more than 10 were found for every metabolite-associated SNP, indicating strong IVs efficacy; consistent causal direction support from the Steiger test results provides confidence to our relaxed threshold choice. Secondly, even if the database utilized for this research contained two sizable populations with East Asian and European ancestry, it is still not representative of the entire world’s population. Furthermore, although MR analyses provide valuable insights into etiology, it is important to note that our findings should be validated by rigorous randomized controlled trials and basic research before clinical application.
5 Conclusion
In conclusion, this MR study showed that 3 known blood metabolites are causally associated with esophageal cancer. Three metabolic pathways that may be associated with the development of esophageal cancer were also identified. The identification of these serum metabolites offers important new information on the early detection, avoidance, and therapy of esophageal cancer in addition to the planning of upcoming clinical trials, even if more validation of the experimental data is still required. The pathophysiology and etiology of esophageal cancer may also be explored in reference form using the combined genomic and metabolomic MR study.
Data availability statement
The original contributions presented in the study are included in the article/Supplementary Material, further inquiries can be directed to the corresponding authors.
Author contributions
CJ: Methodology, Writing–original draft, Writing–review and editing. DY: Supervision, Writing–original draft, Writing–review and editing. MM: Resources, Writing–original draft, Writing–review and editing. QX: Data curation, Writing–original draft, Writing–review and editing. YO: Conceptualization, Writing–original draft, Writing–review and editing. FK: Funding acquisition, Software, Writing–original draft, Writing–review and editing. YJ: Funding acquisition, Writing–original draft, Writing–review and editing.
Funding
The author(s) declare that financial support was received for the research, authorship, and/or publication of this article. This study was supported by the National Natural Science Foundation of China Grant Program (Grant/Award Number: 82104553).
Acknowledgments
The authors express their gratitude to all researchers for sharing their GWAS data, as their priceless contributions have advanced this field of study.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmolb.2024.1391419/full#supplementary-material
References
Argilés, J. M., Busquets, S., Stemmler, B., and López-Soriano, F. J. (2014). Cancer cachexia: understanding the molecular basis. Nat. Rev. Cancer 14 (11), 754–762. doi:10.1038/nrc3829
Benjamini, Y., and Hochberg, Y. (1995). Controlling the false discovery rate: a practical and powerful approach to multiple testing. J. R. Stat. Soc. Ser. B Stat. Methodol. 57 (1), 289–300. doi:10.1111/j.2517-6161.1995.tb02031.x
Boef, A. G., Dekkers, O. M., and le Cessie, S. (2015). Mendelian randomization studies: a review of the approaches used and the quality of reporting. Int. J. Epidemiol. 44 (2), 496–511. doi:10.1093/ije/dyv071
Bowden, J., Davey Smith, G., and Burgess, S. (2015). Mendelian randomization with invalid instruments: effect estimation and bias detection through Egger regression. Int. J. Epidemiol. 44, 512–525. doi:10.1093/ije/dyv080
Bowden, J., Davey Smith, G., Haycock, P. C., and Burgess, S. (2016). Consistent estimation in mendelian randomization with some invalid instruments using a weighted median estimator. Genet. Epidemiol. 40, 304–314. doi:10.1002/gepi.21965
Burgess, S., Butterworth, A., and Thompson, S. G. (2013). Mendelian randomization analysis with multiple genetic variants using summarized data. Genet. Epidemiol. 37 (7), 658–665. doi:10.1002/gepi.21758
Burgess, S., Daniel, R. M., Butterworth, A. S., and Thompson, S. G.EPIC-InterAct Consortium (2015). Network Mendelian randomization: using genetic variants as instrumental variables to investigate mediation in causal pathways. Int. J. Epidemiol. 44 (2), 484–495. doi:10.1093/ije/dyu176
Burgess, S., and Thompson, S. G. (2017). Interpreting findings from Mendelian randomization using the MR-egger method. Eur. J. Epidemiol. 32 (5), 377–389. doi:10.1007/s10654-017-0255-x
Cao, Y. (2019). Adipocyte and lipid metabolism in cancer drug resistance. J. Clin. Invest 129 (8), 3006–3017. doi:10.1172/JCI127201
Castro, C., Peleteiro, B., and Lunet, N. (2018). Modifiable factors and esophageal cancer: a systematic review of published meta-analyses. J. Gastroenterol. 53 (1), 37–51. doi:10.1007/s00535-017-1375-5
Choi, K. W., Chen, C. Y., Stein, M. B., Klimentidis, Y. C., Wang, M. J., Koenen, K. C., et al. (2019). Assessment of bidirectional relationships between physical activity and depression among adults: a 2-sample Mendelian randomization study. JAMA Psychiat 76 (4), 399–408. doi:10.1001/jamapsychiatry.2018.4175
Cohen, J. F., Chalumeau, M., Cohen, R., Korevaar, D. A., Khoshnood, B., and Bossuyt, P. M. (2015). Cochran's Q test was useful to assess heterogeneity in likelihood ratios in studies of diagnostic accuracy. J. Clin. Epidemiol. 68 (3), 299–306. doi:10.1016/j.jclinepi.2014.09.005
Currie, E., Schulze, A., Zechner, R., Walther, T. C., and Farese, R. V. (2013). Cellular fatty acid metabolism and cancer. Cell Metab. 18 (2), 153–161. doi:10.1016/j.cmet.2013.05.017
Doll, R., Peto, R., Wheatley, K., Gray, R., and Sutherland, I. (1994). Mortality in relation to smoking: 40 years' observations on male British doctors. BMJ 309, 901–911. doi:10.1136/bmj.309.6959.901
Freedman, N. D., Abnet, C. C., Leitzmann, M. F., Mouw, T., Subar, A. F., Hollenbeck, A. R., et al. (2007). A prospective study of tobacco, alcohol, and the risk of esophageal and gastric cancer subtypes. Am. J. Epidemiol. 165 (12), 1424–1433. doi:10.1093/aje/kwm051
Giambartolomei, C., Vukcevic, D., Schadt, E. E., Franke, L., Hingorani, A. D., Wallace, C., et al. (2014). Bayesian test for colocalisation between pairs of genetic association studies using summary statistics. PLoS Genet. 10 (5), e1004383. doi:10.1371/journal.pgen.1004383
Global Burden of Disease 2019 Cancer Collaboration Kocarnik, J. M., Compton, K., Dean, F. E., Fu, W., Gaw, B. L., et al. (2022). Cancer incidence, mortality, years of life lost, years lived with disability, and disability-adjusted life years for 29 cancer groups from 2010 to 2019: a systematic analysis for the global burden of disease study 2019. JAMA Oncol. 8 (3), 420–444. doi:10.1001/jamaoncol.2021.6987
Gu, J., Hu, X., Shao, W., Ji, T., Yang, W., Zhuo, H., et al. (2016). Metabolomic analysis reveals altered metabolic pathways in a rat model of gastric carcinogenesis. Oncotarget 7, 60053–60073. doi:10.18632/oncotarget.11049
Hemani, G., Tilling, K., and Davey, S. G. (2017). Orienting the causal relationship between imprecisely measured traits using GWAS summary data. PLoS Genet. 13 (11), e1007081. doi:10.1371/journal.pgen.1007081
Hu, Q., Hao, P., Liu, Q., Dong, M., Gong, Y., Zhang, C., et al. (2019). Mendelian randomization studies on atherosclerotic cardiovascular disease: evidence and limitations. Sci. China Life Sci. 62 (6), 758–770. doi:10.1007/s11427-019-9537-4
Huang, S., Guo, Y., Li, Z., Zhang, Y., Zhou, T., You, W., et al. (2020). A systematic review of metabolomic profiling of gastric cancer and esophageal cancer. Cancer Biol. Med. 17 (1), 181–198. doi:10.20892/j.issn.2095-3941.2019.0348
Huang, S., Rutkowsky, J. M., Snodgrass, R. G., Ono-Moore, K. D., Schneider, D. A., Newman, J. W., et al. (2012). Saturated fatty acids activate TLR-mediated proinflammatory signaling pathways. J. Lipid Res. 53 (9), 2002–2013. doi:10.1194/jlr.D029546
Jewison, T., Su, Y., Disfany, F. M., Liang, Y., Knox, C., Maciejewski, A., et al. (2014). SMPDB 2.0: big improvements to the small molecule pathway database. Nucleic Acids Res. 42 (D1), D478–D484. doi:10.1093/nar/gkt1067
Kanehisa, M., Goto, S., Sato, Y., Furumichi, M., and Tanabe, M. (2012). KEGG for integration and interpretation of large-scale molecular data sets. Nucleic Acids Res. 40 (D1), D109–D114. doi:10.1093/nar/gkr988
Kim, K. B., Yang, J. Y., Kwack, S. J., Park, K. L., Kim, H. S., Ryu, D. H., et al. (2010). Toxicometabolomics of urinary biomarkers for human gastric cancer in a mouse model. J. Toxicol. Environ. Health A 73 (21-22), 1420–1430. doi:10.1080/15287394.2010.511545
Lagergren, J. (2011). Influence of obesity on the risk of esophageal disorders. Nat. Rev. Gastroenterol. Hepatol. 8 (6), 340–347. doi:10.1038/nrgastro.2011.73
Lee, J. Y., Zhao, L., Youn, H. S., Weatherill, A. R., Tapping, R., Feng, L., et al. (2004). Saturated fatty acid activates but polyunsaturated fatty acid inhibits Toll-like receptor 2 dimerized with Toll-like receptor 6 or 1. J. Biol. Chem. 279 (17), 16971–16979. doi:10.1074/jbc.M312990200
Madak-Erdogan, Z., Band, S., Zhao, Y. C., Smith, B. P., Kulkoyluoglu-Cotul, E., Zuo, Q., et al. (2019). Free fatty acids rewire cancer metabolism in obesity-associated breast cancer via estrogen receptor and mTOR signaling. Cancer Res. 79 (10), 2494–2510. doi:10.1158/0008-5472.CAN-18-2849
Matsunaga, S., Nishiumi, S., Tagawa, R., and Yoshida, M. (2018). Alterations in metabolic pathways in gastric epithelial cells infected with Helicobacter pylori. Microb. Pathog. 124, 122–129. doi:10.1016/j.micpath.2018.08.033
Meng, Y., Zhang, J., Zhang, F., Ai, W., Zhu, X., Shu, G., et al. (2017). Lauric acid stimulates mammary gland development of pubertal mice through activation of GPR84 and PI3K/akt signaling pathway. J. Agric. Food Chem. 65 (1), 95–103. doi:10.1021/acs.jafc.6b04878
Mika, A., Pakiet, A., Czumaj, A., Kaczynski, Z., Liakh, I., Kobiela, J., et al. (2020). Decreased triacylglycerol content and elevated contents of cell membrane lipids in colorectal cancer tissue: a lipidomic study. J. Clin. Med. 9 (4), 1095. doi:10.3390/jcm9041095
Nath, A., Li, I., Roberts, L. R., and Chan, C. (2015). Elevated free fatty acid uptake via CD36 promotes epithelial-mesenchymal transition in hepatocellular carcinoma. Sci. Rep. 5, 14752. doi:10.1038/srep14752
Palmer, T. M., Lawlor, D. A., Harbord, R. M., Sheehan, N. A., Tobias, J. H., Timpson, N. J., et al. (2012). Using multiple genetic variants as instrumental variables for modifiable risk factors. Stat. Meth Med. Res. 21 (3), 223–242. doi:10.1177/0962280210394459
Pierce, B. L., and Burgess, S. (2013). Efficient design for Mendelian randomization studies: subsample and 2-sample instrumental variable estimators. Am. J. Epidemiol. 178 (7), 1177–1184. doi:10.1093/aje/kwt084
Prabhu, A., Obi, K. O., and Rubenstein, J. H. (2014). The synergistic effects of alcohol and tobacco consumption on the risk of esophageal squamous cell carcinoma: a meta-analysis. Am. J. Gastroenterol. 109 (6), 822–827. doi:10.1038/ajg.2014.71
Rappaport, S. M., Barupal, D. K., Wishart, D., Vineis, P., and Scalbert, A. (2014). The blood exposome and its role in discovering causes of disease. Environ. Health Perspect. 122 (8), 769–774. doi:10.1289/ehp.1308015
Reichel, M., Hönig, S., Liebisch, G., Lüth, A., Kleuser, B., Gulbins, E., et al. (2015). Alterations of plasma glycerophospholipid and sphingolipid species in male alcohol-dependent patients. Biochim. Biophys. Acta 1851 (11), 1501–1510. doi:10.1016/j.bbalip.2015.08.005
Sakaue, S., Kanai, M., Tanigawa, Y., Karjalainen, J., Kurki, M., Koshiba, S., et al. (2021). A cross-population atlas of genetic associations for 220 human phenotypes. Nat. Genet. 53 (10), 1415–1424. doi:10.1038/s41588-021-00931-x
Schönfeld, P., and Wojtczak, L. (2016). Short- and medium-chain fatty acids in energy metabolism: the cellular perspective. J. Lipid Res. 57 (6), 943–954. doi:10.1194/jlr.R067629
Shaheen, N., and Ransohoff, D. F. (2002). Gastroesophageal reflux, barrett esophagus, and esophageal cancer: scientific review. JAMA 287 (15), 1972–1981. doi:10.1001/jama.287.15.1972
Shin, S. Y., Fauman, E. B., Petersen, A. K., Krumsiek, J., Santos, R., Huang, J., et al. (2014). An atlas of genetic influences on human blood metabolites. Nat. Genet. 46 (6), 543–550. doi:10.1038/ng.2982
Staley, J. R., Blackshaw, J., Kamat, M. A., Ellis, S., Surendran, P., Sun, B. B., et al. (2016). PhenoScanner: a database of human genotype-phenotype associations. Bioinformatics 32, 3207–3209. doi:10.1093/bioinformatics/btw373
Sung, H., Ferlay, J., Siegel, R. L., Laversanne, M., Soerjomataram, I., Jemal, A., et al. (2021). Global Cancer Statistics 2020: GLOBOCAN estimates of incidence and mortality worldwide for 36 cancers in 185 countries. CA Cancer J. Clin. 71 (3), 209–249. doi:10.3322/caac.21660
Suzuki, M., Takaishi, S., Nagasaki, M., Onozawa, Y., Iino, I., Maeda, H., et al. (2013). Medium-chain fatty acid-sensing receptor, GPR84, is a proinflammatory receptor. J. Biol. Chem. 288 (15), 10684–10691. doi:10.1074/jbc.M112.420042
Thomas, D. C., and Conti, D. V. (2004). Commentary: the concept of mendelian randomization. Int. J. Epidemiol. 33, 21–25. doi:10.1093/ije/dyh048
Tullio, V., Gasperi, V., Catani, M. V., and Savini, I. (2020). The impact of whole grain intake on gastrointestinal tumors: a focus on colorectal, gastric, and esophageal cancers. Nutrients 13 (1), 81. doi:10.3390/nu13010081
van der Veen, J. N., Kennelly, J. P., Wan, S., Vance, J. E., Vance, D. E., and Jacobs, R. L. (2017). The critical role of phosphatidylcholine and phosphatidylethanolamine metabolism in health and disease. Biochim. Biophys. Acta Biomembr. 1859 (9 Pt B), 1558–1572. doi:10.1016/j.bbamem.2017.04.006
Verbanck, M., Chen, C. Y., Neale, B., and Do, R. (2018). Detection of widespread horizontal pleiotropy in causal relationships inferred from Mendelian randomization between complex traits and diseases. Nat. Genet. 50 (5), 693–698. doi:10.1038/s41588-018-0099-7
Veugelers, P. J., Porter, G. A., Guernsey, D. L., and Casson, A. G. (2006). Obesity and lifestyle risk factors for gastroesophageal reflux disease, Barrett esophagus and esophageal adenocarcinoma. Dis. Esophagus 19 (5), 321–328. doi:10.1111/j.1442-2050.2006.00602.x
Wallace, C. (2013). Statistical testing of shared genetic control for potentially related traits. Genet. Epidemiol. 37, 802–813. doi:10.1002/gepi.21765
Wang, L., Chen, J., Chen, L., Deng, P., Bu, Q., Xiang, P., et al. (2013). 1H-NMR based metabonomic profiling of human esophageal cancer tissue. Mol. Cancer 12, 25. doi:10.1186/1476-4598-12-25
Xiao, G., He, Q., Liu, L., Zhang, T., Zhou, M., Li, X., et al. (2022). Causality of genetically determined metabolites on anxiety disorders: a two-sample Mendelian randomization study. J. Transl. Med. 20 (1), 475. doi:10.1186/s12967-022-03691-2
Yang, J., Yan, B., Zhao, B., Fan, Y., He, X., Yang, L., et al. (2020). Assessing the causal effects of human serum metabolites on 5 major psychiatric disorders. Schizophr. Bull. 46 (4), 804–813. doi:10.1093/schbul/sbz138
Yokoyama, A., Kato, H., Yokoyama, T., Tsujinaka, T., Muto, M., Omori, T., et al. (2002). Genetic polymorphisms of alcohol and aldehyde dehydrogenases and glutathione S-transferase M1 and drinking, smoking, and diet in Japanese men with esophageal squamous cell carcinoma. Carcinogenesis 23 (11), 1851–1859. doi:10.1093/carcin/23.11.1851
Yun, Z., Guo, Z., Li, X., Shen, Y., Nan, M., Dong, Q., et al. (2023). Genetically predicted 486 blood metabolites in relation to risk of colorectal cancer: a Mendelian randomization study. Cancer Med. 12 (12), 13784–13799. doi:10.1002/cam4.6022
Zhang, Y., Kong, L., Guan, J., Li, J., Wang, Y., Wang, Y., et al. (2023). 16S rDNA sequence analysis of the characteristics of gut Flora in patients with esophageal squamous cell carcinoma. Chin. Gen. Pract. 26, 2496–2502. doi:10.12114/j.issn.1007-9572.2022.0832-1
Keywords: Mendelian randomization, blood metabolites, colocalization analysis, SNPs, esophageal cancer
Citation: Jia C, Yi D, Ma M, Xu Q, Ou Y, Kong F and Jia Y (2024) Genetically predicted 486 blood metabolites in relation to risk of esophageal cancer: a Mendelian randomization study. Front. Mol. Biosci. 11:1391419. doi: 10.3389/fmolb.2024.1391419
Received: 25 February 2024; Accepted: 16 September 2024;
Published: 02 October 2024.
Edited by:
Jacopo Troisi, University of Salerno, ItalyReviewed by:
Irina Nakashidze, Shota Rustaveli State University, GeorgiaMarcos Edgar Herkenhoff, University of São Paulo, Brazil
Copyright © 2024 Jia, Yi, Ma, Xu, Ou, Kong and Jia. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Fanming Kong, kongfanming08@163.com; Yingjie Jia, jiayingjie1616@sina.com
†Present addresses: Caiyan Jia, Tianjin Cancer Institute of Traditional Chinese Medicine, Tianjin, China; Dan Yi, Tianjin Cancer Institute of Traditional Chinese Medicine, Tianjin, China; Mingze Ma, Tianjin Cancer Institute of Traditional Chinese Medicine, Tianjin, China; Qian Xu, Tianjin Cancer Institute of Traditional Chinese Medicine, Tianjin, China; Yan Ou, Tianjin Cancer Institute of Traditional Chinese Medicine, Tianjin, China; Fanming Kong, Tianjin Cancer Institute of Traditional Chinese Medicine, Tianjin, China; Yingjie Jia, Tianjin Cancer Institute of Traditional Chinese Medicine, Tianjin, China
‡These authors have contributed equally to this work