Skip to main content

Causal associations between liver traits and Colorectal cancer: a Mendelian randomization study

Abstract

Objective

This study aimed to investigate the causal associations between several liver traits (liver iron content, percent liver fat, alanine transaminase levels, and liver volume) and colorectal cancer (CRC) risk using a Mendelian randomization (MR) approach to improve our understanding of the disease and its management.

Methods

Genetic variants were used as instrumental variables, extracted from genome-wide association studies (GWAS) datasets of liver traits and CRC. The Two-Sample MR package in R was used to conduct inverse variance weighted (IVW), MR Egger, Maximum likelihood, Weighted median, and Inverse variance weighted (multiplicative random effects) MR approaches to generate overall estimates of the effect. MR analysis was conducted with Benjamini-Hochberg method-corrected P values to account for multiple testing (P < 0.013). MR-PRESSO was used to identify and remove outlier genetic variants in Mendelian randomization (MR) analysis. The MR Steiger test was used to assess the validity of the assumption that exposure causes outcomes. Leave-one-out validation, pleiotropy, and heterogeneity testing were also conducted to ensure the reliability of the results. Multivariable MR was utilized for validation of our findings using the IVW method while also adjusting for potential confounding or pleiotropy bias.

Results

The MR analysis suggested a causal effect between liver volume and a reduced risk of CRC (OR 0.60; 95% CI, 0.44–0.82; P = 0.0010) but did not provide evidence for causal effects of liver iron content, percent liver fat, or liver alanine transaminase levels. The MR-PRESSO method did not identify any outliers, and the MR Steiger test confirmed that the causal direction of the analysis results was correct in the Mendelian randomization analysis. MR results were consistent with heterogeneity and pleiotropy analyses, and leave-one-out analysis demonstrated the overall values obtained were consistent with estimates obtained when all available SNPs were included in the analysis. Multivariable MR was utilized for validation of our findings using the IVW method while also adjusting for potential confounding or pleiotropy bias.

Conclusion

The study provides tentative evidence for a causal role of liver volume in CRC, while genetically predicted levels of liver iron content, percent liver fat, and liver alanine transaminase levels were not associated with CRC risk. The findings may inform the development of targeted therapeutic interventions for colorectal liver metastasis (CRLM) patients, and the study highlights the importance of MR as a powerful epidemiological tool for investigating causal associations between exposures and outcomes.

Peer Review reports

Introduction

Colorectal cancer(CRC) is the third most common type of cancer in the world and more than 50% of patients with CRC will develop colorectal liver metastases (CRLMs) [1, 2]. The liver is the most common site of colorectal metastasis, with 15–20% of patients being candidates for hepatectomy [3, 4]. Approximately one-third of patients will develop liver metastases within three years of diagnosis, which contributes to the overall poor prognosis and survival rates [5]. Surgical resection is currently the treatment of choice for colorectal liver metastasis (CRLM) and has been shown to be the only potentially curative therapy [6, 7]. The combination of advances in medical therapy, such as systemic chemotherapy (CTX), and the role of surgery in the treatment of metastatic disease has had a positive impact on the prognosis, with an increase in median survival and cure rates [8]. However, the biology of liver metastases is extremely complex and involves many interactions between tumor cells and the liver microenvironment, leading to significant challenges in the management of CRLM [9]. Consequently, there is an urgent need to explore novel therapeutic strategies and elucidate the underlying mechanisms linking liver traits to CRC development and progression.

The liver plays a crucial role in the natural course of CRC due to its anatomical and physiological importance. Metastatic cells in CRC acquire properties in addition to those required to become neoplastic because they are able to successfully dissociate, disseminate and colonize secondary sites: motility and invasion, the ability to modulate the secondary site or the local micro-environment, plasticity and the ability to colonize secondary tissues [10, 11]. In addition to the mucus and intestinal epithelial barrier, the intestine also has an intestinal vascular barrier (GVB), which acts as a gatekeeper to control the entry of molecules and microorganisms into the systemic circulation [12, 13]. Under hazardous conditions, harmful intestinal pathogens can cross the epithelial barrier, damaging the GVB and ultimately allowing the bacteria or their components to spread into the bloodstream, eventually reaching the liver [12, 13].

The liver and lung are the most common sites for CRC metastasis, with the liver’s predominant blood supply arising from the confluence of the gastrointestinal (GI) tract, which supplies blood vessels via the hepatic portal vein. For metastases arising from the colon and proximal parts of the rectum, the portal system directs blood flow directly to the liver [14]. In hematogenous dissemination, platelets and neutrophils protect circulating tumour cells (CTCs) from being eliminated by NK cells. Subsequent entry of CTCs into the liver microvasculature may trigger a proinflammatory cascade that induces the secretion of chemokines by Kupffer and stellate cells and upregulation of vascular adhesion receptors, resulting in the adhesion of CTCs to the liver microvasculature [15]. Furthermore, tumour cell interaction with the tumour microenvironment (TME) plays an important role in CRLM [16]. By remodeling the TME, cancer cells can induce the formation of pre-metastatic niches [17, 18]. For example, activation of hepatic stellate cells (α-HSC) is the most common biological process in secondary or primary liver cancer [19]. By remodeling and depositing extracellular matrix (ECM), α-HSCs can influence CRC cell growth and invasion [20].

Mendelian randomization (MR) is a powerful epidemiological approach that leverages genetic variants as instrumental variables to investigate causal associations between exposures and outcomes [21]. By capitalizing on the random allocation of genetic variants during meiosis, MR studies can provide robust evidence for causality, while minimizing the influence of potential confounding factors and reverse causation [22, 23]. Based on this, in the context of CRC and liver traits, MR presents a unique opportunity to untangle the intricate relationship between liver characteristics and colorectal cancer risk, which could, in turn, shed light on the biological mechanisms underpinning CRLM and inform the development of targeted therapeutic interventions, as depicted in Fig. 1.

Fig. 1
figure 1

Study design of Mendelian randomization between liver traits and CRC. The solid lines represent the association between the instrumental variables and exposure as well as the association between exposure and outcome. Dash lines with a cross means that the association meets two basic assumptions of Mendelian randomization: (i) the genetic variants are independent of confounders between exposure and outcomes, (ii) the genetic variants only influence the outcome via exposure

This MR study aims to investigate the causal associations between several liver traits (liver iron content, percent liver fat, alanine transaminase levels and liver volume) and CRC risk, with the ultimate goal of improving our understanding of the disease and its management. For this aim, we established an effective and reliable data processing framework, which made it possible to generate instrumental variables (IVs) for MR analysis. By selecting exposure-related genetic variations as IVs, MR eliminates confounding factors [24]. And then, using an inverse variance weighted MR approach, we calculated the pool result based on IVs. To further guarantee the reliability of our method, MR-Egger and Weighted median analysis and leave-one-out validation were conducted in this step too. MR-Egger is able to detect some violations of the standard assumptions of the instrumental variables and provide an estimate of the effect that is not subject to these violations. The weighted median analysis is used to combine data from multiple genetic variants into a single causal estimate [25]. Moreover, pleiotropy and heterogeneity testing were used to minimize false-positive conclusions and the risk of producing unreliable results [26]. What’s more, multivariable MR was utilized for validation of our findings using the IVW method, while also adjusting for potential confounding or pleiotropy bias. Finally, a full report was generated to provide an overview of the results of our study.

Materials and methods

Genetic instruments

Summary-level data from the GWAS dataset was fundamental for MR analysis. IVs were extracted from several liver traits (liver iron content, percent liver fat, alanine transaminase levels and liver volume) related GWAS datasets, and CRC related GWAS datasets were then used in the further analysis. To meet the MR assumptions and reduce the bias, the summarized GWAS data was processed. Subsequently, MR analysis involving an inverse variance weighted MR approach, leave-one-out validation, and MR-Egger analysis was used to comprehensively assess the causal effect of several liver traits on the risk of the development of CRC.

Data sources

We used data from the UK Biobank for the study of liver iron content, percent liver fat, alanine transaminase levels and liver volume [27]. The UK Biobank has approval from the North West Multicenter Research Ethics Committee (https://www.ukbiobank.ac.uk/ethics/), and these ethics regulations cover the work in this study. Summary statistics data for CRC were available from a genome-wide association meta-analysis of 177,028 European-descent individuals (3,022 cases and 174,006 noncases) from FinnGen [28]. A summary of the characterization and covariates adjusted for the datasets employed in this study can be found in the supplementary information (Table S1).

Our study analyzed the summarized data from genome-wide association studies (GWAS) to create instrumental variables (IVs) for MR analysis. In this context, IVs are genetic variants [29]. [30] The criteria we followed for selecting these instrumental variables are as follows: (1) Each individual SNP must show a significant correlation with several liver traits (liver iron content, percent liver fat, alanine transaminase levels, and liver volume) across the entire genomic region, with a significance threshold of P < 5.0 × 10− 8. (2) We assessed the linkage disequilibrium (LD) among SNPs through the OpenGWAS API, which contains LD reference panels for the five super-populations in the 1000 Genomes reference panel. The reference panel includes only bi-allelic SNPs with a minor allele frequency (MAF) > 0.010. The entire process utilizes the PLINK clumping method, where SNPs in linkage disequilibrium (LD) within a specific window are pruned. The SNP with the lowest p-value is kept. In this study, only those SNPs that produced the most significant results, among those with an R2 < 0.0010 within a clumping window of 10,000 kb, were kept. 7 SNPs were utilized as instrumental variables for liver iron content, 10 SNPs for percent liver fat, 228 SNPs for liver volume, and 10 SNPs for liver alanine transaminase levels. All instrumental variables in the supplementary information.

Two-sample mendelian randomization

Using genetic variants as instrumental variables, MR is a technique that is gaining popularity for determining the causal relationships between exposures and outcomes. In order to guarantee the accuracy of causal reasoning, MR is predicated on a number of essential premises. Assumption 1: There is a significant association between the genetic variation used as an instrumental variable and the relevant exposure. This premise is crucial because it guarantees that the genetic variation is a reliable substitute for the exposure, enabling the estimation of a causal impact. Estimates of the causal impact may be skewed if this assumption is broken. Assumption 2: The genetic variation is not linked to any confounding elements that might skew the relationship between exposure and result. The second fundamental tenet of MR is that there are no confounding variables linked with the genetic variant used as an instrumental variable that might skew the relationship between exposure and outcome. Assumption 3: The genetic variant has no other impact on the outcome besides how it affects the exposure [31]. The genetic variant used as an instrumental variable only influences the result through its impact on the exposure, and not through any other pathway, according to the third fundamental tenet of MR.

In this study, the Two-Sample MR package (version 0.5.6) in R 4.1.3 (R Foundation for Statistical Computing, Vienna, Austria) was used to conduct MR analysis, evaluating the causal effect of various liver traits (including liver iron content, percent liver fat, alanine transaminase levels, and liver volume) on CRC risk mediated by each individual instrumental SNP. Five methods, including inverse variance-weighted (IVW), MR-Egger, Maximum likelihood, Weighted median, and Inverse variance-weighted (multiplicative random effects), were employed to generate overall estimates of the effect. In order to test whether the associations were likely to be causal, we used inverse variance weighted (IVW) two-sample MR as our primary analysis to estimate the effect of the exposure factors on CRC [23]. The IVW approach utilized a meta-analysis framework to combine Wald estimates for each SNP to derive an overall estimate of the effect, assuming the absence of horizontal pleiotropy, with the IVW results being unbiased. However, if heterogeneity exists, random-effects IVW models are employed. The Maximum likelihood method, similar to IVW, assumes the absence of heterogeneity and horizontal pleiotropy. In the presence of these assumptions, the results will be unbiased, with smaller standard errors compared to IVW. To assess the presence of pleiotropy, the MR-Egger regression approach is utilized, based on the assumption of instrument strength independent of direct effect (InSIDE), allowing for the assessment of the intercept term. A zero-intercept term indicates the absence of horizontal pleiotropy, making the result of the MR-Egger regression consistent with IVW. The weighted median method permits accurate estimation of causal association even when up to 50% of instrumental variables are invalid. When the InSIDE hypothesis is violated, the weighted model estimate has greater power to detect a causal effect, less bias, and lower type I error rates than MR-Egger regression.

The strength of instrumental variables (IVs) is evaluated by computing the F-statistic using the formula F = R2×(N-1-K)/[(1 - R2)×K], where R2 represents the proportion of variance in the exposure explained by the genetic variants, N represents the sample size, and K represents the number of instruments. A corresponding F-statistic greater than 10 is considered indicative of the absence of significant weak instrumental bias [32]. Furthermore, MR-PRESSO corrects for pleiotropy and outlier effects, while the MR Steiger test evaluates the directionality of the causal effect [33]. Heterogeneity and pleiotropy tests assess the homogeneity and pleiotropic effects of instrumental variables, respectively [34]. Finally, leave-one-out analysis tests the robustness of the causal estimates obtained from the MR analysis by iteratively removing each instrumental variable to identify potential heterogeneous SNPs. Together, these methods help to improve the validity and reliability of causal inference in our MR analysis.

IVW multivariable mendelian randomization

IVW multivariable MR was utilized to assess the causal impact of different liver traits, encompassing liver iron content, percent liver fat, alanine transaminase (ALT) levels, and liver volume, on the risk of colorectal cancer. In this approach, instruments were chosen for each exposure variable. Subsequently, all exposures corresponding to those genetic variants were collectively regressed against the outcome variable, with weighting based on the inverse variance of the outcome. The threshold for instrument inclusion, denoted as pval_threshold, was set at 5.0 × 10− 8. Moreover, decisions were made about whether to estimate intercepts or use instrument-specific estimates for each exposure variable. This comprehensive analysis seeks to unravel the potential causal relationships between liver traits and CRC risk while considering the complexities of confounding and bias.

Results

Causal effect of liver traits on CRC

Summary statistics of instrumental SNPs as genetic instrumental variables (IVs) for several kinds of liver traits (liver iron content, percent liver fat, alanine transaminase levels and liver volume) were presented in Table S2. These SNPs were not associated with CRC and have no linkage disequilibrium (LD) associations. Each line of the table included 9items related to the SNP, such as the SNPs, effect allele, other allele, beta coefficients and standard error of the SNP on the risk of the corresponding liver traits and CRC.

As shown in Fig. 2, our analysis suggested a causal effect between liver volume and a reduced risk of colorectal cancer (the five MR methods demonstrated consistent effect directions and statistical significance with a a Benjamini-Hochberg method-corrected P value < 0.013 and did not provide evidence for causal effects of liver iron content, percent liver fat, or liver alanine transaminase levels. Large liver volumes may represent a decreased risk of CRC (Beta − 0.51; OR 0.60; 95% CI, 0.44–0.82; P = 0.0010). The risk of CRC lacked evidence of association with liver iron content, percent liver fat, and liver alanine transaminase levels. Liver iron content (Beta 0.16; OR 1.2; 95% CI, 0.97–1.4; P = 0.10), percent liver fat (Beta 0.020; OR 1.0; 95% CI, 0.82–1.2; P = 0.79), liver alanine transalanine transase levels (Beta − 0.060; OR 0.94; 95% CI, 0.32–2.8; P = 0.91) are provided in Supplemental Fig. 1.

Fig. 2
figure 2

Associations between genetically predicted liver and CRC. CI indicates confidence interval; OR, odds ratio. MR-Egger is able to detect some violations of the standard assumptions of the instrumental variables and provide an estimate of the effect that is not subject to these violations. The weighted median analysis is used to combine data from multiple genetic variants into a single causal estimate.( Liver iron content in mg/g, percent liver fat in %, liver volume in L, Liver enzyme levels (alanine transaminase) in IU/L)

For the causal analysis of liver volume and colorectal cancer risk, the MR-PRESSO method did not identify any outliers (Beta − 0.51; P = 0.011; Global Test P = 0.47). The MR Steiger test confirmed that the causal direction of the analysis results was correct (P = 6.0 × 10− 101). And no horizontal pleiotropy was detected according to MR-Egger intercepts, and no evident heterogeneity was identified; detailed information is provided in Table 1. In conclusion, our findings provide tentative evidence for a causal role of liver volume in CRC, although further research is needed to confirm these results.

Table 1 Heterogeneity and pleiotropy of individual single nucleotide polymorphisms for Mendelian randomization

The leave-one-out analysis was performed in order to assess the influence of individual SNPs on the estimates. The results demonstrate that the overall values obtained by the leave-one-out analysis were consistent with estimates obtained when all available SNPs were included in the analysis. Further information on the leave-one-out validation results can be found in Supplemental Fig. 2.

The findings derived from the application of IVW multivariable MR exhibited coherence with the outcomes obtained through the Two-sample MR. SNPs employed in IVW multivariable MR and their beta, standard error, and P value for exposure and outcome can be found in Additional File 1. Notably, the analysis revealed a significant and consistent trend for liver volume’s association with colorectal cancer risk. Specifically, liver volume demonstrated a noteworthy protective effect against colorectal cancer onset, as indicated by the estimated Beta coefficient of -0.45 (SE = 0.16, P = 0.0050). A graphical representation of the results is available in Supplemental Fig. 3, showcasing scatter plots generated from the IVW multivariable MR analysis. This alignment between the IVW multivariable MR and the Two-sample MR reinforces the observed connection between liver volume and colorectal cancer risk, further validating the significance of this relationship.

Discussion

This study provided genetic evidence for a causal relationship between liver volume and CRC risk and found a possible causal role for liver volume in reducing CRC risk. However, there was no evidence for a causal role of liver iron content, percent liver fat, or liver alanine transaminase levels in CRC risk.

The characteristics of the liver are closely related to the onset and progression of many diseases. In diabetic ketoacidosis (DKA) disease, it has been shown that alanine aminotransferase (ALT) has a greater correlation with the disease in patients aged 65 years and older and in obese patients [35]. Other studies have shown a strong causal relationship between higher levels of liver fat and the risk of type 2 diabetes, which is consistent with recent MR studies showing a causal relationship between non-alcoholic fatty liver disease or its markers (ALT and AST) and a higher risk of type 2 diabetes [30].

However, the iron content of the liver does influence the course of the disease. High iron levels are an independent inducer or cofactor of hepatocellular carcinoma (HCC). And there is a correlation between hepatic iron level and hepatocellular carcinoma among patients with end-stage liver disease, while the strongest correlation between hepatic iron level and hepatocellular carcinoma is found among patients with biliary cirrhosis and hepatitis C [36]. Patients with cancer have a history of chemotherapy, which puts them at increased risk of liver toxicity and pancytopenia, leading to elevated liver fat and iron. Moreover,liver volume is another factor that influences disease progression. Despite limitations in younger patients, liver volume remains a traditional imaging biomarker for polycystic liver [37]. In patients with diagnosis of hilar or distal bile duct cancer, liver volume is a useful tool to assess the efficacy of biliary stenting [38]. After gastrectomy in patients with gastric adenocarcinoma, the second and third segments of the liver show significant atrophy compared with the rest of the liver and the whole liver, and the volume reduction is poorly correlated with time [39].

The liver, the most common site of metastasis in CRC, has seen an increase in the number of patients with resectable colorectal liver metastases (CRLM) over the past two decades. However, the reasons for the increase in numbers are due to advances in intensive systemic chemotherapy regimens associated with targeted therapies or immunotherapy, on the one hand, and technical improvements that increase the size of future liver remnant (FLR) and facilitate parenchyma-sparing surgery, on the other [40].

Our results showed that liver volume has a protective effect against CRC. Liver volume, which reflects the number of liver parenchymal cells and is a quantitative measure of liver reserve function from a morphological point of view, is an indicator for assessing disease severity and predicting outcome [41]. The liver consists of approximately 80% hepatocytes and 20% non-hepatocytes [42], and there is a direct relationship between liver volume and cell content. As the most regenerative organ in the body, the liver’s powerful regenerative capacity comes at the expense of the liver’s reserve function [43]. When the liver is injured, there is massive necrosis, shrinkage, and structural collapse of the liver tissue, leading to a reduction in volume and dysfunction of the hepatic reserve [41].

Hepatocytes are the major parenchymal cells responsible for metabolic functions and most circulating plasma proteins such as albumin, transporters, protease inhibitors, coagulation factors and modulators of immune complexes and inflammation are expressed by hepatocytes [44]. They metabolise amino acids, metals and endogenous compounds such as haem and bilirubin, and control the homeostasis of molecules such as glucose/glycogen, triglycerides, cholesterol, bile acids and vitamins A and D [45]. Amino acids are not only building blocks of proteins, they are also intermediate metabolites that power a variety of biosynthetic pathways [46]. Glutamine is an essential amino acid that is used by cancer cells for biosynthetic, bioenergetic and antioxidant purposes [47]. For the bioenergetic and biosynthetic needs of cells, glutamate serves as an important carbon source. When proliferating cells use glutamine-derived carbon to enter the TCA cycle, most of the citrate produced is exported to the cytosol where it is converted to acetyl-CoA, a precursor for the biosynthesis of fatty acids and cholesterol. Oxaloacetate from the TCA cycle is also used to synthesise aspartate and asparagine [48].

In addition to glutamine metabolism, serine and glycine metabolism are also important mediators in cancer cell development [49]. Glycine is a precursor of serine and is involved in a wide range of metabolic processes in both humans and animals as a component of glutathione and as a substrate in the synthesis of purines and proteins [50]. The consumption of glycine and the expression of the mitochondrial glycine biosynthesis pathway are strongly correlated with the rate of proliferation of cancer cells in general and are a source of amino acids required for cancer growth [51, 52]. The primary route for the disposal of glycine is catalysed by the glycine cleavage system, which is located in the inner membrane of the mitochondria in the liver [53]. In contrast, the glycine decarboxylation system of cancer cells can break down glycine in the presence of the glycine decarboxylase complex to form ammonia, carbon dioxide, and methylenetetrahydrofolate [54]. Serine catabolism is initiated by the interconversion of serine to glycine, producing a large number of one-carbon units that activate one-carbon metabolism, forming an intracellular double-loop pathway linking the methionine and folate cycles via methionine synthase. The glycine cleavage system and choline catabolism are additional pathways to one-carbon metabolism, and the preferred carbon source driving one-carbon metabolism in cancer cells is the uptake and subsequent breakdown of glucose into various biomass precursors, including the synthesis of serine and glycine [55].In macromolecular synthesis, serine can load transfer RNAs (tRNAs) for protein synthesis, serve as a precursor for amino acids such as cysteine and glycine, provide head groups for sphingolipid and phospholipid synthesis, or donate 1 C units for nucleotide synthesis [56].Thus, the volume of the liver is determined by the hepatocytes, which have the role of metabolizing amino acids. Amino acid metabolism is an important mediator of cancer development and may have the same effect on colorectal cancer.

In summary, these findings provide insights into the biological mechanisms underlying the complex relationship between liver characteristics and CRC development, and may inform the development of targeted therapeutic interventions for patients with CRLM. Our study highlights the importance of MR as a powerful epidemiological tool for investigating causal associations between exposures and outcomes. However, there are some limitations to this study that should be considered. First, although we examined the effect of certain circulating liver characteristics on CRC, it was not possible to examine possible non-linear relationships with CRC. Another limitation is that the associations between liver traits and CRC were derived from a European cohort, making it difficult to extrapolate to other populations. An increase in the diversity of the populations studied is required for even modest differences between populations in the contribution of shared variation to complex traits [57]. Therefore, the results we present can only be extended to European populations. What’s more, in contrasting the covariate adjustments between the cohorts, it’s evident that the “Liver iron content” cohort incorporates specific imaging-related covariates, such as the imaging center, scan date, and scan time, which are absent in other cohorts, because it is a radiological feature. Additionally, while the “Liver alanine transaminase levels” cohort adjusts for a notable 40 genetic principal components, the “Colorectal cancer” cohort includes adjustments for just 10 genetic principal components and “Liver iron content” doesn’t consider PCs as a covariate. Such differences in covariate adjustments can have implications. For instance, the variance in the number of genetic principal components used between cohorts suggests that there might be different degrees of control for potential genetic confounding. These disparities underscore the importance of being cautious when interpreting results, especially in cross-cohort comparisons. Additionally, the demographic disparities within the dataset present a critical issue for consideration. It is evident that three of the datasets (Liver iron content, Percent liver fat, and Liver volume) predominantly encompass cohorts from the UK Biobank, consisting of individuals of European ancestry, meticulously ascertained based on self-identified “White British” ancestry. In contrast, the cohort for Liver alanine transaminase levels consists of a discovery sample drawn from individuals of European descent from the UK Biobank, while the replication sample encompasses individuals of European lineage from both the Netherlands and the United States. Furthermore, the cohorts for Colorectal cancer are exclusively composed of Finnish individuals. Notably, while the exposure groups exhibit uniformity in composition and origin, predominantly from the UK Biobank (excluding the replication sample), the outcome population originates from Finland. Even though both populations are of European descent, the potential biases introduced by their distinct backgrounds in Mendelian randomization cannot be dismissed. For future endeavors, it would be prudent to employ datasets with more congruent demographic backgrounds to validate these findings, ensuring their robustness.

Larger and more comprehensive MR studies are warranted to better understand the intrinsic relationship between liver traits and CRC. More research is needed to clarify the association between hepatic characteristics and the risk of CRC and to explore the underlying mechanisms.

Data Availability

The datasets generated and/or analyzed during the current study are available in the [GWAS summary data] repository, [Liver iron content: https://gwas.mrcieu.ac.uk/datasets/ebi-a-GCST90016674/, percent liver fat: https://gwas.mrcieu.ac.uk/datasets/ebi-a-GCST90016673/, alanine transaminase levels: https://gwas.mrcieu.ac.uk/datasets/ebi-a-GCST004940/, liver volume: https://gwas.mrcieu.ac.uk/datasets/ebi-a-GCST90016666/, Colorectal Cancer: https://gwas.mrcieu.ac.uk/datasets/finn-b-C3_COLORECTAL_EXALLC/]

References

  1. Hernandez-Alejandro R, Ruffolo LI, Sasaki K, Tomiyama K, Orloff MS, Pineda-Solis K, et al. Recipient and Donor outcomes after living-Donor Liver Transplant for Unresectable Colorectal Liver metastases. JAMA Surg. 2022;157:524.

    Article  PubMed  PubMed Central  Google Scholar 

  2. Li W, Lin S, He Y, Wang J, Pan Y. Deep learning survival model for Colorectal cancer patients (DeepCRC) with Asian clinical data compared with different theories. Arch Med Sci. 2023;19:264–9.

    Article  PubMed  PubMed Central  Google Scholar 

  3. Adam R, Vinet E. Regional treatment of Metastasis: Surgery of colorectal liver metastases. Ann Oncol off J Eur Soc Med Oncol. 2004;15(Suppl 4):iv103–106.

    Article  Google Scholar 

  4. Parkin DM, Bray F, Ferlay J, Pisani P. Global cancer statistics, 2002. CA Cancer J Clin. 2005;55:74–108.

    Article  PubMed  Google Scholar 

  5. Leporrier J, Maurel J, Chiche L, Bara S, Segol P, Launoy G. A population-based study of the incidence, management and prognosis of hepatic metastases from Colorectal cancer. Br J Surg. 2006;93:465–74.

    Article  CAS  PubMed  Google Scholar 

  6. Fong Y, Fortner J, Sun RL, Brennan MF, Blumgart LH. Clinical score for predicting recurrence after hepatic resection for metastatic Colorectal cancer: analysis of 1001 consecutive cases. Ann Surg. 1999;230:309–18. discussion 318–321.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  7. Ito K, Govindarajan A, Ito H, Fong Y. Surgical treatment of hepatic colorectal Metastasis: evolving role in the setting of improving systemic therapies and ablative treatments in the 21st century. Cancer J Sudbury Mass. 2010;16:103–10.

    Article  Google Scholar 

  8. Gallagher DJ, Kemeny N. Metastatic Colorectal cancer: from improved survival to potential cure. Oncology. 2010;78:237–48.

    Article  PubMed  Google Scholar 

  9. Van den Eynden GG, Majeed AW, Illemann M, Vermeulen PB, Bird NC, Høyer-Hansen G, et al. The multifaceted role of the microenvironment in liver Metastasis: biology and clinical implications. Cancer Res. 2013;73:2031–43.

    Article  PubMed  Google Scholar 

  10. Welch DR, Hurst DR. Defining the hallmarks of Metastasis. Cancer Res. 2019;79:3011–27.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  11. Lambert AW, Pattabiraman DR, Weinberg RA. Emerg Biol Principles Metastasis Cell. 2017;168:670–91.

    CAS  Google Scholar 

  12. Spadoni I, Pietrelli A, Pesole G, Rescigno M. Gene expression profile of endothelial cells during perturbation of the gut vascular barrier. Gut Microbes. 2016;7:540–8.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  13. Spadoni I, Zagato E, Bertocchi A, Paolinelli R, Hot E, Di Sabatino A, et al. A gut-vascular barrier controls the systemic dissemination of bacteria. Science. 2015;350:830–4.

    Article  CAS  PubMed  Google Scholar 

  14. Poturnajova M, Furielova T, Balintova S, Schmidtova S, Kucerova L, Matuskova M. Molecular features and gene expression signature of metastatic Colorectal cancer (review). Oncol Rep. 2021;45:10.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  15. Dou R, Liu K, Yang C, Zheng J, Shi D, Lin X, et al. EMT-cancer cells-derived exosomal miR-27b-3p promotes circulating tumour cells-mediated Metastasis by modulating vascular permeability in Colorectal cancer. Clin Transl Med. 2021;11:e595.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  16. Zhao S, Mi Y, Zheng B, Wei P, Gu Y, Zhang Z et al. Highly-metastatic Colorectal cancer cell released miR‐181a‐5p‐rich extracellular vesicles promote liver Metastasis by activating hepatic stellate cells and remodelling the tumour microenvironment. J Extracell Vesicles. 2022;11.

  17. Illemann M, Eefsen RHL, Bird NC, Majeed A, Osterlind K, Laerum OD, et al. Tissue inhibitor of matrix metalloproteinase-1 expression in Colorectal cancer liver metastases is associated with vascular structures. Mol Carcinog. 2016;55:193–208.

    Article  CAS  PubMed  Google Scholar 

  18. Schütte M, Risch T, Abdavi-Azar N, Boehnke K, Schumacher D, Keil M, et al. Molecular dissection of Colorectal cancer in pre-clinical models identifies biomarkers predicting sensitivity to EGFR inhibitors. Nat Commun. 2017;8:14262.

    Article  PubMed  PubMed Central  Google Scholar 

  19. Affo S, Yu L-X, Schwabe RF. The role of Cancer-Associated fibroblasts and fibrosis in Liver Cancer. Annu Rev Pathol. 2017;12:153–86.

    Article  CAS  PubMed  Google Scholar 

  20. Berman RS, Portera CA, Ellis LM. Biology of Liver metastases. In: Talamonti MS, Pappas SG, editors. Liver-Directed Therapy for primary and metastatic liver tumors. Boston, MA: Springer US; 2001. pp. 183–206.

    Chapter  Google Scholar 

  21. van der Graaf A, Claringbould A, Rimbert A, Westra H-J, Li Y, Wijmenga C, et al. Mendelian randomization while jointly modeling cis genetics identifies causal relationships between gene expression and lipids. Nat Commun. 2020;11:4930.

    Article  PubMed  PubMed Central  Google Scholar 

  22. Burgess S, Foley CN, Zuber V. Inferring Causal relationships between Risk factors and outcomes from Genome-Wide Association Study Data. Annu Rev Genomics Hum Genet. 2018;19:303–27.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  23. Lawlor DA, Harbord RM, Sterne JAC, Timpson N, Davey Smith G. Mendelian randomization: using genes as instruments for making causal inferences in epidemiology. Stat Med. 2008;27:1133–63.

    Article  PubMed  Google Scholar 

  24. Li Z, Xia L, Li X, He H. Assessing the causality between thyroid and breast Neoplasms: a bidirectional mendelian randomization study. Front Oncol. 2022;12:973161.

    Article  PubMed  PubMed Central  Google Scholar 

  25. Bowden J, Davey Smith G, Haycock PC, Burgess S. Consistent estimation in mendelian randomization with some Invalid instruments using a weighted median estimator. Genet Epidemiol. 2016;40:304–14.

    Article  PubMed  PubMed Central  Google Scholar 

  26. Cho Y, Haycock PC, Sanderson E, Gaunt TR, Zheng J, Morris AP, et al. Exploiting horizontal pleiotropy to search for causal pathways within a mendelian randomization framework. Nat Commun. 2020;11:1010.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  27. Collins R. What makes UK Biobank special? Lancet Lond Engl. 2012;379:1173–4.

    Article  Google Scholar 

  28. Fadista J, Skotte L, Karjalainen J, Abner E, Sørensen E, Ullum H, et al. Comprehensive genome-wide association study of different forms of hernia identifies more than 80 associated loci. Nat Commun. 2022;13:3200.

    Article  PubMed  PubMed Central  Google Scholar 

  29. Zhao H, Jin X. Causal associations between dietary antioxidant vitamin intake and Lung cancer: a mendelian randomization study. Front Nutr. 2022;9:965911.

    Article  PubMed  PubMed Central  Google Scholar 

  30. Martin S, Sorokin EP, Thomas EL, Sattar N, Cule M, Bell JD, et al. Estimating the effect of liver and pancreas volume and Fat Content on risk of Diabetes: a mendelian randomization study. Diabetes Care. 2022;45:460–8.

    Article  PubMed  Google Scholar 

  31. Peng H, Wang S, Wang M, Ye Y, Xue E, Chen X, et al. Nonalcoholic fatty Liver Disease and Cardiovascular Diseases: a mendelian randomization study. Metabolism. 2022;133:155220.

    Article  CAS  PubMed  Google Scholar 

  32. Bowden J, Holmes MV. Meta-analysis and mendelian randomization: a review. Res Synth Methods. 2019;10:486–96.

    Article  PubMed  PubMed Central  Google Scholar 

  33. Luo J, Xu Z, Noordam R, van Heemst D, Li-Gao R. Depression and Inflammatory Bowel Disease: a bidirectional two-sample mendelian randomization study. J Crohns Colitis. 2022;16:633–42.

    Article  PubMed  Google Scholar 

  34. Zhou D, Jiang Y, Zhong X, Cox NJ, Liu C, Gamazon ER. A unified framework for joint-tissue transcriptome-wide association and mendelian randomization analysis. Nat Genet. 2020;52:1239–46.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  35. Liu Q, Gong C, Geng Y, You J. Elevated alanine transaminase is nonlinearly associated with in-hospital death in ICU-admitted diabetic ketoacidosis patients. Diabetes Res Clin Pract. 2023;197:110555.

    Article  CAS  PubMed  Google Scholar 

  36. Ko C, Siddaiah N, Berger J, Gish R, Brandhagen D, Sterling RK, et al. Prevalence of hepatic iron overload and association with hepatocellular cancer in end-stage Liver Disease: results from the National Hemochromatosis Transplant Registry. Liver Int off J Int Assoc Study Liver. 2007;27:1394–401.

    Google Scholar 

  37. Bae KT, Tao C, Feldman R, Yu ASL, Torres VE, Perrone RD, et al. Volume progression and imaging classification of polycystic liver in early autosomal Dominant polycystic Kidney Disease. Clin J Am Soc Nephrol CJASN. 2022;17:374–84.

    Article  PubMed  Google Scholar 

  38. Lee CH, Kim SH, Kim IH, Kim SW, Lee ST, Kim DG, et al. Endoscopic stenting in bile duct cancer increases liver volume. Gastrointest Endosc. 2014;80:447–55.

    Article  PubMed  Google Scholar 

  39. Ozutemiz C, Obuz F, Taylan A, Atila K, Bora S, Ellidokuz H. Volume change of segments II and III of the liver after gastrectomy in patients with gastric cancer. Diagn Interv Radiol Ank Turk. 2016;22:109–15.

    Article  Google Scholar 

  40. Cusumano C, Deshayes E, Guiu B, De Meeus G, Carrère S, Bouillin A, et al. Prevention of post-hepatectomy Liver Failure after major resection of colorectal liver metastases: is hepato-biliary scintigraphy the optimal tool? Ann Transl Med. 2023;11:202.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  41. Hu C, Jiang N, Zheng J, Li C, Huang H, Li J, et al. Liver volume based prediction model for patients with Hepatitis B virus-related acute-on-chronic Liver Failure. J Hepato-Biliary-Pancreat Sci. 2022;29:1253–63.

    Article  Google Scholar 

  42. Godoy P, Hewitt NJ, Albrecht U, Andersen ME, Ansari N, Bhattacharya S, et al. Recent advances in 2D and 3D in vitro systems using primary hepatocytes, alternative hepatocyte sources and non-parenchymal liver cells and their use in investigating mechanisms of hepatotoxicity, cell signaling and ADME. Arch Toxicol. 2013;87:1315–530.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  43. Thorgersen EB, Barratt-Due A, Haugaa H, Harboe M, Pischke SE, Nilsson PH, et al. The role of complement in Liver Injury, Regeneration, and transplantation. Hepatol Baltim Md. 2019;70:725–36.

    Article  Google Scholar 

  44. Michalopoulos GK. Liver regeneration. J Cell Physiol. 2007;213:286–300.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  45. Tanaka M, Itoh T, Tanimizu N, Miyajima A. Liver stem/progenitor cells: their characteristics and regulatory mechanisms. J Biochem (Tokyo). 2011;149:231–9.

    Article  CAS  PubMed  Google Scholar 

  46. Tabe Y, Lorenzi PL, Konopleva M. Amino acid metabolism in hematologic malignancies and the era of targeted therapy. Blood. 2019;134:1014–23.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  47. Zhang J, Pavlova NN, Thompson CB. Cancer cell metabolism: the essential role of the nonessential amino acid, glutamine. EMBO J. 2017;36:1302–15.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  48. DeBerardinis RJ, Mancuso A, Daikhin E, Nissim I, Yudkoff M, Wehrli S, et al. Beyond aerobic glycolysis: transformed cells can engage in glutamine metabolism that exceeds the requirement for protein and nucleotide synthesis. Proc Natl Acad Sci U S A. 2007;104:19345–50.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  49. Bernhardt S, Bayerlová M, Vetter M, Wachter A, Mitra D, Hanf V, et al. Proteomic profiling of Breast cancer metabolism identifies SHMT2 and ASCT2 as prognostic factors. Breast Cancer Res BCR. 2017;19:112.

    Article  PubMed  Google Scholar 

  50. Dasarathy S, Kasumov T, Edmison JM, Gruca LL, Bennett C, Duenas C, et al. Glycine and urea kinetics in nonalcoholic steatohepatitis in human: effect of intralipid infusion. Am J Physiol - Gastrointest Liver Physiol. 2009;297:G567–75.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  51. Jain M, Nilsson R, Sharma S, Madhusudhan N, Kitami T, Souza AL, et al. Metabolite Profiling identifies a key role for Glycine in Rapid Cancer Cell Proliferation. Science. 2012;336:1040–4.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  52. AbouAitah K, Hassan HA, Swiderska-Sroda A, Gohar L, Shaker OG, Wojnarowicz J, et al. Targeted Nano-Drug delivery of colchicine against Colon Cancer cells by means of mesoporous silica nanoparticles. Cancers. 2020;12:144.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  53. Kikuchi G, Motokawa Y, Yoshida T, Hiraga K. Glycine cleavage system: reaction mechanism, physiological significance, and hyperglycinemia. Proc Jpn Acad Ser B Phys Biol Sci. 2008;84:246–63.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  54. Sun W, Liu R, Gao X, Lin Z, Tang H, Cui H, et al. Targeting serine-glycine-one-carbon metabolism as a vulnerability in cancers. Biomark Res. 2023;11:48.

    Article  PubMed  PubMed Central  Google Scholar 

  55. Geeraerts SL, Heylen E, De Keersmaecker K, Kampen KR. The ins and outs of serine and glycine metabolism in cancer. Nat Metab. 2021;3:131–41.

    Article  CAS  PubMed  Google Scholar 

  56. Li AM, Ye J. Reprogramming of serine, glycine and one-carbon metabolism in cancer. Biochim Biophys Acta Mol Basis Dis. 2020;1866:165841.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  57. Wojcik GL, Graff M, Nishimura KK, Tao R, Haessler J, Gignoux CR, et al. Genetic analyses of diverse populations improves discovery for complex traits. Nature. 2019;570:514–8.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

Download references

Acknowledgements

The authors express their gratitude to the participants and investigators of the FinnGen study. The authors thank all investigators for sharing this data.

Funding

This research was funded by the Guangdong Basic and Applied Basic Research Foundation(2022A1515110163)。.

Author information

Authors and Affiliations

Authors

Contributions

Conceptualization: [Ying Ni] and [Wenkai Wang]; Methodology: [Ying Ni], [Wenkai Wang] and [Yongming Liu]; software: [Ying Ni]and [Yongming Liu]; Formal analysis and investigation: [Ying Ni] and [Wenkai Wang]; Writing—original draft preparation: [Ying Ni] and [Wenkai Wang]; Writing—review and editing: [Ying Ni]; Funding acquisition: [Yun Jiang]; Resources: [Yun Jiang]; Supervision: [Yun Jiang]. All authors have read and agreed to the published version of the manuscript.

Corresponding author

Correspondence to Yun Jiang.

Ethics declarations

Consent for publication

Not Applicable.

Conflict of interest

No.

Ethical Statement

Not applicable.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Electronic supplementary material

Below is the link to the electronic supplementary material.

Supplementary Material 1

Supplementary Material 2

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Ni, Y., Wang, W., Liu, Y. et al. Causal associations between liver traits and Colorectal cancer: a Mendelian randomization study. BMC Med Genomics 16, 316 (2023). https://doi.org/10.1186/s12920-023-01755-w

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12920-023-01755-w

Keywords