Abstract
Background or Purpose The interplay between the sex-specific differences in tumor metabolome and colorectal cancer (CRC) prognosis has never been studied and represents an opportunity to improve patient outcomes. This study aims to examine the link between tumor metabolome and prognosis by sex for CRC patients.
Methods Using untargeted metabolomics analysis, abundances of 91 metabolites were obtained from primary tumor tissues from 197 patients (N=95 females, N=102 males) after surgical colectomy for stage I-III CRC. Cox Proportional Hazards (PH) regression models were applied to estimate the associations between tumor metabolome and 5-year overall survival (OS) and 5-year recurrence-free survival (RFS), and their interactions with sex.
Results Eleven metabolites had significant sex differences in their associations with 5-year OS, and five metabolites for 5-year RFS (Pinteraction < .05). The metabolites asparagine and serine had sex interactions for both OS and RFS. Furthermore, sex-specific differences were found in the associations between prognosis and metabolic pathways. Notably, in the asparagine synthetase (ASNS)-catalyzed asparagine synthesis pathway, asparagine was associated with substantially poorer OS (hazard ratio [HR] = 6.39, 95% confidence interval [CI] = 1.78-22.91, P = .004) and RFS (HR = 4.36, 95% CI = 1.39-13.68, P = .01) for female patients only (Pinteraction, OS = .02, Pinteraction, RFS = .003). Similar prognostic disadvantages in females were seen in lysophospholipid and polyamine synthesis.
Discussion or Conclusions Unique metabolite profiles indicated increased asparagine synthesis was associated with poorer prognosis for females only, providing insights into precision medicine for CRC treatment stratified by sex.
Introduction
According to the U.S. Centers for Disease Control and Detection (CDC), colorectal cancer (CRC) was the third leading cause of cancer-related deaths in the United States in 20211. Prognosis is partly dependent on stage at presentation according to TNM classification system based on the American Joint Committee on Cancer (AJCC)2 and Union for International Cancer Control (UICC)3. Additional clinicopathologic variables can also impact survival such as primary tumor location4,5, body mass index (BMI)6, and diabetes conditions7.
The application of biomarkers in a clinical context has potential to improve prognosis based on gene mutations8,9, microRNAs10, and DNA methylation11,12. During the last decade, cancer studies have benefited from the implementation of metabolomics13. CRC metabolomic studies have shown the correlation of certain metabolite profiles to cancer survival and recurrence using biological samples of urine and tissue14. Some plasma metabolites are associated with recurrence in stage II and III patients15, showing the potential value of metabolomics in finding prognostic biomarkers of CRC.
In addition, sex-related differences in CRC prognosis have become increasingly important in cancer research16,17. Female patients have a higher prevalence of right-sided colon cancer (RCC), which is associated with poorer overall survival (OS)18,19. However, male CRC patients have a survival disadvantage in most subgroups by anatomic location and age20. At present, female and male patients do not receive different sex-specific therapies.
Previously, our untargeted metabolomics analyses in a large series of CRC tumor tissues revealed sex-specific metabolic subphenotypes that could influence tumor progression and response to therapeutics21. Important metabolic pathways in these subphenotypes included asparagine synthesis, energy production (glycolysis and pentose phosphate pathway (PPP)), lysophospholipid synthesis, methionine metabolism, and polyamine synthesis21. Particularly, female patients with RCC showed increased asparagine levels21, which can be produced intracellularly by an ATP- and glutamine-dependent reaction catalyzed by asparagine synthetase (ASNS)22. By examining The Cancer Genome Atlas (TCGA) data, we discovered that high ASNS expression was associated with poorer OS for female CRC patients21. Additionally, increased ASNS expression has been linked to aggressive tumors and poorer prognosis in other cancer types23,24. These findings suggest that sex plays a vital role in CRC prognosis together with the covariates of metabolism and anatomic tumor location. A clear understanding of these sex-specific metabolic differences represents a golden opportunity to prolong the survival of CRC patients.
In this study, we examine whether there is a link between CRC tumor metabolites and prognosis stratified by sex of the patient. This is the first study to reveal that sex-related differences exist in primary tumor metabolite profiles that are associated with OS and recurrence-free survival (RFS) in patients who have undergone curative intent colectomy. We also identify metabolic pathways which are associated with poor survival only in female CRC.
Methods
Sample Collection and Metabolite Measurements
Fresh frozen tumor tissue samples were acquired after curative intent colectomy for patients (N=197) with RCCs (cecum, ascending colon and hepatic flexure) or LCCs (splenic flexure, descending colon, sigmoid, rectosigmoid) with stage I-III CRC (1991-2000). No patients received neoadjuvant chemotherapy. Finally, abundances of 91 metabolites were obtained as published previously21. Metabolites were extracted and analyzed by hydrophilic interaction liquid chromatography mass spectrometry (HILIC-MS) and reverse phase liquid chromatography mass spectrometry (RPLC-MS)-based metabolomics as previously described21. More information is presented in S1 Text. Yale University IRB determined that this study was not considered Human Subjects Research and did not require IRB review (IRB/HSC# 1612018746).
Statistical Analysis
We included age, sex, anatomic primary tumor location, chemotherapy history, and clinical stage as covariates. Missing values of metabolites abundances were replaced using median imputation. Multivariable Cox proportional hazard (Cox PH) regression models were constructed to evaluate the associations between prognosis with individual metabolite abundance per standard deviation (SD) increase on a log2-scale, adjusted for covariates in either all patients or stratified by sex. Assessment of associations between sex-specific metabolic pathways and prognosis was conducted by multivariate Cox PH models that included important metabolites involved in these pathways. Sex-specific difference was examined in multivariate Cox PH models by including an interaction between sex and the metabolite. 5-year OS, 5-year RFS were calculated. Due to the absence of death events in female patients at clinical stage I, we recoded the clinical stages I and II as “early stage”, and III as “late stage”. Patients with adjuvant chemotherapy were coded as having chemotherapy history. Survival analyses were conducted using package “survival” in R (version 4.0.4). A 2-sided P value less than .05 was considered statistically significant.
Results
Population Characteristics
Baseline characteristics of 197 patients (N=102 males, N=95 females), including 5-year OS and 5-year RFS, are provided in Table 1. Median follow-up time since date of surgery was 74.8 months (0.1-169.2). Older age, administration of chemotherapy, and advanced clinical stage were inversely related to OS (Table 1). Each subgroup by anatomic tumor location and clinical stage, demographic characteristics including race/ethnicity are displayed in S1 Table. Prognosis among different subgroups is shown in S2 Table.
Sex-specific differences in the associations between individual metabolite abundances and CRC prognosis
Eighteen metabolites were significantly associated with OS in males or females or in both sexes (S3 Table), adjusted for anatomic location, clinical stage, chemotherapy history, and age, similarly 25 metabolites had significant correlations with RFS (S4 Table). Only carnitine and hypoxanthine remained significantly associated with RFS for males after false discovery rate (FDR) adjustment.
Multivariate Cox PH models on single metabolites and their associations with OS or RFS revealed dichomotous findings by sex. Fig 1 presents the metabolites (adenosine, asparagine, citrulline, glycerol 3-phosphate (Gro3P), lysophosphatidylcholine (LysoPC)(16:0), ornithine, succinate, threonine, UDP-D-glucose, uracil, and xanthosine) that had sex differences in their associations with CRC prognosis (Pinteraction < .05). Among these 11 metabolites, succinate was significantly associated with improved OS for female patients (HROS = 0.35, 95% CI=0.12-0.99; P = .047), while it was associated with poorer OS for males (HROS = 1.91, 95% CI=1.23-2.96, P = .004) (Fig 1A). Five metabolites (argininosuccinic acid, asparagine, creatinine, hypoxanthine, and serine) showed significance for RFS and sex-specific differences (Fig 1B). Of these metabolites, asparagine exhibited sex differences for both OS and RFS (Pinteraction,OS = .03, Pinteraction,RFS = .007), and similarly for serine (Pinteraction,OS = .04, Pinteraction,RFS = .01) (S3 and S4 Table). Asparagine was significantly associated with better CRC prognosis in males (HROS = 0.72, 95% CI=0.54-0.96, P = .03; HRRFS = 0.74, 95% CI=0.56-0.97, P = .03), while its associations with worse prognosis in females trended towards significance (POS = .13, PRFS = .09, S3 and S4 Table). Similarly, serine was significantly associated with better prognosis in males (HROS = 0.55, 95% CI=0.37-0.81, P = .002. HRRFS = 0.57, 95% CI=0.40-0.83, P = .003), while it has no significant associations with prognosis in females (S3 and S4 Table).
(A) 5-year overall survival (OS) and (B) 5-year recurrence-free survival (RFS). LysoPC: Lysophosphatidylcholine. All metabolites were observed with significant sex interactions (Pinteraction < .05). Hazard ratio (HR) for each metabolite (log2-transformed abundance) was calculated by multivariate Cox PH regression adjusted for anatomic location, chemotherapy history, clinical stage, and age (continuous). Error bars represent 95% confidence intervals (CIs). 95% CIs marked with asterisks indicate significant associations between the metabolite and the corresponding prognosis (P values < .05): A metabolite with HR < 1 was associated with better prognosis; a metabolite with HR > 1 was associated with worse prognosis. The x-axes are log-scaled.
Sex-specific Differences in CRC Prognosis Associated with Asparagine Synthesis
We examined sex-specific associations between the ASNS-catalyzed asparagine synthesis pathway and CRC prognosis (S5 and S6 Table). Significant interactions with sex were detected in associations between asparagine abundance and both OS and RFS (Pinteraction,OS = .02, Pinteraction,RFS = .003) as summarized as Model 1, where asparagine abundance was treated as a continuous variable (Table 2). In Model 1, asparagine was significantly associated with poor OS (HROS = 6.39, 95% CI=1.78-22.91, P = .004) and poor RFS (HRRFS = 4.36, 95% CI=1.39-13.68, P = .01) only in females, while asparagine was only significantly associated with better OS only in males (HROS = 0.57, 95% CI=0.36-0.91, P = .02) after adjustment for clinical stages, chemotherapy history, age, anatomic location, and other crucial metabolites involved in ASNS expression (Table 2).
An additional multivariate Cox PH regression model was built; Model 2, where asparagine abundance was dichotomized to low levels and high levels. There was a significant interaction between sex and dichotomized asparagine levels for RFS and a similar interaction that was approaching significance for OS (Table 2) (Pinteraction,OS = .052, Pinteraction,RFS = .03).
Fig 2 displays the Cox adjusted survival curves for dichotomized asparagine levels based on Model 2, illustrating that female CRC patients with a high level of ASNS-catalyzed asparagine had significantly poorer OS (Fig 2A) and poorer RFS (Fig 2C) compared with female patients with a low level. Strikingly, levels of asparagine were not associated with OS (Fig 2B) or RFS (Fig 2D) for male CRC patients.
(A) OS of females, (B) OS of males, (C) RFS of females, and (D) RFS of males. Results were based on Model 2 (Table 2) using multivariate Cox PH regression, adjusted for anatomic location, chemotherapy history, clinical stages, age (continuous), and other crucial metabolites involved in ASNS-catalyzed asparagine synthesis (aspartate, glutamate, glutamine, and AMP were treated as continuous variables and were log2 transformed). P values for the associations between dichotomized asparagine levels and the prognosis in Model 2.
Sex-Specific Differences in CRC Prognosis Associated with Glycolysis and the Pentose Phosphate Pathway
Considering glycolysis and the PPP metabolism, the association between Gro3P and OS for females trended towards significance (P = .08), and there was no association in males (P = .76) (S8 Table). Gro3P was significantly associated with a lower risk of 5-year recurrence only for males (HRRFS = 0.27, 95% CI=0.08-0.88, P = .03), though no sex differences were detected (Pinteraction= .25) (S7 and S8 Table).
Sex-Specific Differences in CRC Prognosis Associated with Lysophospholipid Synthesis
We analyzed the relationships between CRC prognosis and important metabolites involved in lysophospholipid synthesis. The interactions between sex and LysoPC (16:0), LysoPC (18:1), lysophosphatidylethanolamine (LysoPE) (18:0), and LysoPE (22:5) in the association with OS were significant (S9 Table). In the multivariate analysis on OS for lysophospholipids synthesis by sex, LysoPC (16:0) was associated with poor OS only for females (HROS = 2.07, 95% CI=1.06-4.07, P = .03), while it was not significantly associated with OS for males (P = .32) (S10 Table). In the models of RFS for lysophospholipid synthesis by sex, LysoPE (22:5) presented a robust trend toward a significant sex difference in its association with RFS (Pinteraction = .052, S9 Table). LysoPE (22:5) was associated with increased risks of 5-year recurrence of CRC for females (HRRFS = 5.38, 95% CI=1.02-28.42, P = .048), while it was associated with reduced risk of recurrence for males (HRRFS = 0.10, 95% CI=0.02-0.62, P = .01) (S10 Table).
Sex-specific Differences in CRC Prognosis Associated with Methionine Metabolism
Analysis of methionine metabolism for all patients indicated that serine favored better OS (HROS = 0.56, 95% CI=0.34-0.92, P = .02) and RFS (HRRFS = 0.49, 95% CI=0.29-0.83, P = .008) (S11 Table). After further stratification by sex, serine was associated with better OS in males (HRRFS = 0.37, 95% CI=0.19-0.72, P = .004) and improved RFS in males as well (HRRFS = 0.49, 95% CI=0.24-0.98, P = .04), while its associations with CRC prognosis in females were not significant (S12 Table). The interaction between serine and sex in terms of OS and RFS were both significant (Pinteraction,OS = .04, Pinteraction,RFS = .045, S11 Table). Methionine had a significant sex interaction (S11 Table), but it was not significantly associated with OS or RFS in either sex (S12 Table). The additional metabolites in this pathway did not have significant sex interactions (S11 Table).
Sex-specific Differences in CRC Prognosis Associated with Polyamine Synthesis
Lastly, we examined metabolites involved in polyamine synthesis. Arginine had a marked trend towards a significant sex interaction with OS (Pinteraction= 0.07), and a significant sex interaction for RFS (Pinteraction = .04) (S13 Table). Further stratified by sex, arginine was associated with worse OS only for females (HROS = 5.63, 95% CI=1.33-23.78, P =.02), while its association with OS in males was not significant (P = .87) (S14 Table). For RFS, arginine was associated with higher risks of 5-year CRC recurrence only for females (HRRFS = 3.22, 95% CI=1.14-9.13, P = .03) while it had no significant association with RFS for males (P = .24) (S14 Table). Citrulline suggested significant sex interactions with both OS and RFS (Pinteraction,OS= .002, Pinteraction,RFS = .02) (S13 Table): it had a significant association with poorer OS for females (HROS = 2.57, 95% CI=1.04-6.33, P = .04) but not for males (P = .18), while it was not significantly associated with RFS in either sex (S14 Table). Other critical metabolites in polyamine synthesis did not have significant sex interactions (S13 Table).
Discussion
To our knowledge, this is the first study to show that sex-specific differences in the CRC tumor metabolome are linked to prognosis. It is also the first application of a novel approach to examine metabolic pathways and their associations to CRC prognosis. By applying state of the art mass spectrometry to fresh frozen specimens acquired operatively in a large colorectal biorepository, we show that females have a distinct metabolite profile, underscored by asparagine metabolism, that correlates with prognosis after surgical colectomy. This suggests that females and males have unique colorectal tumor metabolism which contributes to sex-specific clinical outcomes.
Our study revealed that 11 tumor metabolites were associated with CRC 5-year OS and five metabolites for 5-year RFS, and these associations were sex-specific. Asparagine and serine were associated with both 5-year OS and 5-year RFS in a sex-specific manner. Further examination of asparagine within the de novo asparagine synthesis pathway (catalyzed by ASNS), showed consistently poor OS and RFS with high asparagine abundance for females. This supports previous findings revealing poor female survival and high ASNS expression from the TCGA21. In addition, we previously hypothesized that asparagine is linked to tumor aggressiveness in females with RCCs as it correlates with increased threonine and serine uptake, which indicates aggressive tumor phenotypes21. ASNS-catalyzed asparagine production has also been shown to be crucial for in vitro cancer cell proliferation by amino acid homeostasis25, and the availability of asparagine in vitro and in vivo propagate the metastatic progression of breast cancer26.
In CRC, SRY-box 12 (SOX12) expression promotes cell proliferation and metastasis, and facilitates ASNS expression27. In addition, expression of the mutant KRAS gene correlates with a marked decrease in aspartate levels and increased asparagine levels due to upregulated ASNS expression, which indicated that ASNS might be a novel therapeutic target for KRAS-mutant CRC28,29. Moreover, SLC25A22 serves as an essential metabolic regulator for CRC progression30 by promoting the synthesis of aspartate-derived amino acids (asparagine) in KRAS-mutant CRC cells31. Furthermore, increased dietary asparagine in animals promotes metastatic progression in breast cancer, and L-asparaginase treatment or dietary asparagine restriction inhibits metastasis32. Hence, both internally produced asparagine and external exposure to asparagine through diet may influence CRC tumor progression. However, a review study pointed out that some studies are contradictory regarding whether ASNS expression favors or inhibits CRC progression33, while our findings indicate that such inconsistency may attribute to sex differences.
As mentioned, serine associated with both 5-year OS and RFS in a sex-specific manner. Our data showed that serine favored longer OS and RFS in males but independent of prognosis in females. Serine metabolism has been reported to contribute to CRC metabolism and growth34 and enzymes involved in serine synthesis such as phosphoglycerate dehydrogenase (PHGDH)35 are found in higher levels in colonic tumor tissue than paired normal tissue36. These studies did not consider sex differences which might lend to a potential sex-specific phenomenon.
Our results indicated that sex significantly modified the associations between lysophospholipid synthesis and CRC prognosis. LysoPC (16:0) was associated with worse OS for females in both individual metabolite analysis and the model for lysophospholipid synthesis. Our previous study identified that multiple LysoPCs and LysoPEs were upregulated in female patients with RCC (stage I) only, suggesting that the higher level of lysophospholipids in women with RCC would promote fatty acid supply that is essential for cancer cell growth at early stages21. It has also been discussed that sex-specific estrogen regulation may inhibit CRC cell survival by suppressed triglyceride biosynthesis, a vital lipid marker related to CRC progression17.
Lastly, sex interacts with the associations between polyamine metabolism and CRC prognosis. Arginine and citrulline was associated with worse RFS and OS respectively for females but not for males. A previous study investigated the L-arginine/nitric oxide (NO) pathway in CRC by clinical stage, location of a primary tumor, and sex, and they observed elevated arginine levels in LCCs compared with RCCs, and a drop in arginine levels in the early postoperative period in females exclusively37. Arginine and citrulline are important part of the urea cycle38, where ornithine is produced and metabolized into polyamines by ornithine decarboxylase (ODC)39. ODC is a target for inhibition by difluoromethylornithine (DFMO) (an anti-polyamine therapy for colorectal neoplasia)40, therefore, some of the polyamine precursors could be important for polyamine treatment responses, where sex can play a role. There are NCI funded cooperative group CRC chemoprevention trials targeting at polyamine synthesis (PACES/SWOG)41. Additionally, a review study showed that methionine dietary restriction can suppress tumor formation, potentially through reducing the availability of decarboxylated SAM which is a cofactor for spermidine (a polyamine) production42. Thus, the key metabolic pathways may have great applicability in sex-specific treatment or preventative interventions for CRC.
Some of the study limitations include the number of events in the cohorts compared to sample size, and the multiple categories such as three clinical stages, two sex groups, two anatomic locations, which might weaken statistical power. In addition, some potential confounders that may bias the results were not available in the clinical records for consideration in our analyses. Therefore, future studies from an independent cohort will be valuable.
In conclusion, we apply a novel methodology using tumor tissue to identify sex-specific metabolic pathways associated with CRC prognosis in a surgical cohort. We show that tumor metabolite levels have different associations with patient prognosis by sex, in particular ASNS-catalyzed asparagine production that are linked to poor prognosis only in females. Our study shows that sex-specific tumor metabolism is an important factor for understanding differences in patient prognosis and potentially differences in treatment responses. Accordingly, we recommend that the design of precision medicine approaches consider both tumor metabolism and sex to obtain more robust findings in clinical trials and improve patient outcomes.
Data Availability
The untargeted metabolomics data underlying this article is available at www.ebi.ac.uk/metabolights MTBLS1122, MTBLS1124, MTBLS1129, MTBLS1130. All other data and code will be shared on reasonable request to the corresponding author.
Acknowledgments
This work was funded by NIH 1R21CA223686-01 (CJ, SK) and American Cancer Society Research Scholar Grant 134273-RSG-20-065-01-TBE (CJ). This publication was also made possible by CTSA Grant Number UL1 TR001863 from the National Center for Advancing Translational Science (NCATS), components of the National Institutes of Health (NIH), and NIH roadmap for Medical Research. Its contents are solely the responsibility of the authors and do not necessarily represent the official view of NIH. This work was also supported by the Lampman Research Fund in Yale Surgical Oncology.
Footnotes
Disclosures The authors declare no conflicts of interest.
Reference List
