- Academic Editor
Background: To investigate the gene regulation of tumor cells in the process of different organ metastasis on a xenograft mouse model and screen the genes involved in the organ-target metastasis of tumor cells. Methods: A multi-organ metastasis model was constructed with a human ovarian clear cell carcinoma cell line (ES-2) based on a severe immunodeficiency mouse strain (NCG). Differentially expressed tumor proteins among multi-organ metastases were successfully characterized by microliter liquid chromatography-high-resolution mass spectrometry, sequence-specific data analysis and multivariate statistical data analysis. Liver metastases were selected as typical for subsequent bioinformatic analysis. Selected liver metastasis-specific genes in ES-2 cells were validated by sequence-specific quantitation including high resolution-multiple reaction monitoring quantification at protein level and quantitative real-time polymerase chain reaction at mRNA level. Results: From the mass spectrometry data, a total of 4503 human proteins were identified using the sequence-specific data analysis strategy. Of them, 158 proteins were selected as specifically regulated genes in liver metastases for subsequent bioinformatics studies. Based on Ingenuity Pathway Analysis (IPA) pathway analysis and sequence-specific quantitation, Ferritin light chain (FTL), lactate dehydrogenase A (LDHA) and long-chain-fatty-acid–CoA ligase 1 (ACSL1) were finally validated as specifically upregulated proteins in liver metastases. Conclusions: Our work provides a new approach to analyze gene regulation in tumor metastasis in xenograft mouse model. In presence of a large number of mouse protein interference, we validated the up-regulation of human ACSL1, FTL and LDHA in ES-2 liver metastases, which reflects the adaptive regulation of tumor cells to the liver microenvironment through metabolic reprogramming.
Metastasis is a process in which tumor cells disperse from the primary site. The tumor cell disseminates to distant areas through blood vessels, lymph vessels or direct invasion [1, 2]. In other tissues or organs, secondary tumors form with a histological type similar to that of the primary tumor . Early-stage tumor metastasis-related studies mostly focus on a series of molecular events occurring in the early and middle stages of metastasis from tumor growth in situ to the transport survival stage of the vasculature . However, the mechanism of extravasation to select a secondary organ and colonize distant sites during the late steps is still poorly understood . The earliest understanding of the tendency of organs with tumor metastasis can be traced back to the “seed-soil” doctrine proposed by surgeon Stephen Paget in 1889 . It was called the three landmark theories in tumor metastasis research along with Ewing’s “theory of anatomical mechanisms”  and the “metastatic waterfall theory” proposed by Bross and Blumenson . In recent years, it has been gradually recognized that the interaction between the tumor microenvironment and tumor cells through bidirectional signaling regulation is an important cause and key rate-limiting step for the occurrence of molecular events in the later stage of tumor metastasis and the final formation of metastatic lesions in specific organs .
There is evidence that single or clusters of cells can spread from the primary tumor and that after extravasation to the vasculature, successful metastatic cells must be selected and survive the new tissue microenvironment conducive to survival . Thus, metastatic colonization is not only a product of primary tumor dissemination but also through a complex interaction between the disseminated tumor cells and the tissue microenvironment of the whole organism . In this process, tumor cells exist in a continuous phenotypic state, triggering their intrinsic programs through their interaction with the extracellular microenvironment, conferring their metastatic properties through epigenetic regulation, and overcoming various obstacles to eventually form metastatic lesions. It has been reported that metastatic lesions in different tumor types are often detected in different specific organs (liver, lung, bone and brain) . However, little is known about the tropism of different organs in cancer subtypes. Understanding the complex interplay between the complex molecular-level behavior of tumor cells and the tumor microenvironment leading to metastasis requires the integration of massive molecular signature data collected from in vitro and in vivo experimental models, processing high-throughput data using different computational methods, and integrating information to characterize the metastatic process of tumor cells. However, to date, our understanding of this dynamic complex process remains incomplete.
In recent years, quantitative proteomics has become a powerful tool to study the changes in the protein profiles in cancer metastasis . Mass spectrometry-based proteomics technology is a large-scale, high-throughput proteomics technique that enables the accurate analysis of complex protein mixtures . This technique enables the scale-up search of multiple sets of protein or protein populations of interest and illustrates the complex dynamics of proteins in cells at the overall level through both qualitative and quantitative analysis of proteins . Mass spectrometry, by combining proteomics and bioinformatics technologies, allows for a quantitative and precise systematic analysis of thousands of proteins, making it possible to analyze and characterize the changes in the protein spectrum involved in the metastatic colonization process of tumor cells .
In this research, to explore the phenomenon of gene regulation of tumor cells during metastatic colonization in different organs, a novel research strategy was established based on a sequence-specific proteomics analysis and a mouse model of multiple organ metastasis of human tumor cells. Herein, the ES-2 cell line, an ovarian clear cell carcinoma (OCCC) cell line, was applied to severely immunodeficient mice (NOD/ShiltJGpt-Prkdcem26Cd52 Il2rgem26Cd22/Gpt, NCG strain) to form multiple organ metastases, and an analytical method based on microliter liquid chromatography tandem high resolution mass spectrometry and a sequence-specific data analysis strategy was used to study the differences in protein expression among various organ metastases (Fig. 1) . Through this strategy, we screened out a series of differential proteins for the organ tropism of ES-2 cells, and specific regulatory proteins in liver metastasis were presented as typical in this study. Our results indicated that this sequence-specific proteomics analysis strategy can provide a powerful tool to explore the gene regulation of xenografts or hosts from proteomes by avoiding quantitative interference of homologous proteins. To the best of our knowledge, this is the first study of multiple organ metastasis from one tumor cell line by sequence-specific proteomics analysis. It is of great significance for us to comprehensively understand the molecular mechanism of organ-specific metastasis of tumors.
Experiment design and workflow.
The ES-2 human ovarian clear cell carcinoma cell line and HEK293T cell
line were purchased from the American Culture Collection (ATCC) and cultured in
Dulbecco’s modified Eagle’s medium (DMEM, Gibco, Grand Island, NY, USA) containing 10% Fetal Bovine Serum
(FBS, HyClone, Logan, UT, USA) and 1% penicillin/streptomycin (P/S) (Gibco, Grand Island, NY, USA) at 37 °C with
NCG (NOD/ShiltJGpt-Prkdcem26Cd52Il2rgem26Cd22/Gpt) strain mice, the severe
immunodeficient strain mice lacking mature T cells, B cells, and NK cells used in
the experiment (female, SPF grade,
The mice were killed and left open by cervical dislocation between 18 and 20 days postinoculation. The distribution of metastatic lesions in each organ of 12 mice was observed and then validated by fluorescence imaging. The static dissection of each organ was loaded into a frozen storage tube and frozen in liquid nitrogen.
ES-2 cells were washed twice with cold Phosphate buffered saline (PBS) and lysed in 1000
Sequential window acquisition of all theoretical mass spectra (SWATH-MS)
analysis was performed on a quadrupole time-of-flight spectrometer
(TripleTOF5600, AB/Sciex, Framingham, MA, USA) coupled to a
micro-HPLC (NanoLC 415, Eksigent, Framingham, MA, USA). The LC
flow phase consisted of mobile phase A (2% acetonitrile and 0.1% FA in
ultrapure water) and mobile phase B (98% acetonitrile and 0.1% FA in ultrapure
water) (all purchased from Thermo Fisher Scientific, Waltham, MA, USA). Samples
were first captured using a ChromXP c18 precolumn (5
For HR-MRM experiments, LC‒MS/MS was performed on a quadrupole time-of-flight
spectrometer (TripleTOF5600, AB/Sciex, Framingham, MA, USA) coupled to a micro-HPLC
(NanoLC 415, Eksigent, Framingham, MA, USA). The LC flow phase consisted of
mobile phase A (2% acetonitrile and 0.1% FA in ultrapure water) and mobile
phase B (0.1% FA and 2% ultrapure water in acetonitrile) (all purchased from
Thermo Fisher Scientific, Waltham, MA, USA). The sample was first loaded onto a
ChromXP c18 precolumn (5
Mass spectrometry acquisition data were processed using
Protein-Pilot software (AB/Sciex, v.5.0.1, Framingham, MA,
USA). The nonredundant human UniProtKB/Swiss-Prot protein database (https://www.uniprot.org/uniprotkb?facets=reviewed%3Atrue&query=proteome%3AUP000005640), BSA peptide sequence and mouse UniProtKB/Swiss-Prot protein database (https://www.uniprot.org/uniprotkb?facets=reviewed%3Atrue&query=proteome%3AUP000000589) were combined to form a heterozygous library (https://www.uniprot.org/uniprotkb/P02769/entry) for IDA
data alignment to identify proteins. Search parameters were set as follows:
sample type: identification; Cys alkylation: iodoacetamide; digestion: trypsin;
instrument: TripleTOF 5600; special factors: none; species: Homo sapiens; search
effort: thorough ID; results quality: 0.05 (detected protein threshold) performed
Paragon searches. The false discovery rate (FDR) for proteins and peptides was
set to 0.01 (minimum length set to seven amino acids) for analysis. Data from all
IDA method identification runs performed using ProteinPilot software were merged
into one batch and used for SWATH assay library construction. After database
retrieval, the mouse polypeptides were removed by using the self-compiled R
language software package (R 3.6.0 programming language, R Foundation for Statistical Computing, Vienna, Austria) to establish a human protein
quantification database. Protein quantification information was extracted using
PeakView software 2.2 and the swath plugin, where the BSA peptides were used to
calibrate the retention time calibration. To obtain better quantitative results
using the SWATH acquisition method, the following parameters were set: peptide
filter: 15 peptides per protein; six transitions per peptide; 95% peptide
confidence; 1% FDR threshold; exclusion of modified peptides; XIC option: 8 min
XIC extraction window; 50 ppm XIC width. Missing values in the results are
automatically populated by the PeakView software (AB/Sciex, version 2.2, Framingham, MA, USA) built-in SWAP plugin.
Subsequent statistical analysis of the data was performed using R statistical
software, version 3.6.0. The median method was chosen to normalize the data, and
the p value of protein expression was calculated with the Kruskal‒Wallis test. A
Benjamini‒Hochberg adjustment p value
Bioinformatics analysis was performed using the Gene Ontology (GO) database (http://geneontology.org/) and Kyoto Encyclopedia of Genes and Genome (KEGG) (http://www.genome.jp/kegg/pathway.html) based on the DEPs. In addition, Ingenuity Pathway Analysis (IPA) (http://www.ingenuity.com/products/ipa) was used to determine the canonical pathways among DEPs. Fisher’s exact test was used to calculate the p values.
Considering the inevitable mixing of human tumor cells with murine organ tissue during tumor sampling, we designed primers for differential bases in human and mouse mRNAs, which is used to solve the problem that the human metastasis samples cannot be completely separated from the tissues of NCG mice due to surgical operation during the RT‒PCR experiments. Primer synthesis was commissioned by Beijing Qingke Biological Company. The PCR primer sequences are shown in Supplementary Table 1.
Total RNA from ES-2 cells and metastatic foci of various organs was extracted
using TRIzol reagent (Gibco, Life Technologies, Carlsbad, CA, USA) and
reverse-transcribed into cDNA using a RevertAid™ First Strand cDNA
Synthesis Kit with a DNase I Reverse Transcription Kit (ThermoFisher, USA). PCR
amplification was performed using the TransStart Taq DNA Polymerase PCR Kit
(TransGen Biotech, Inc., Beijing, China). Fluorescence quantitative PCR experiments were performed using a
CFX96 dual-channel real-time PCR instrument (Bio-Rad, Hercules, CA, USA) with 2
Quantitative data are represented as the mean
To facilitate observation and verification of metastases, ES-2 cells were labeled using lentiviral transfection with green fluorescent protein (GFP), and distinctly fluorescent green was observed under fluorescence microscopy. At 18 to 22 weeks after cell injection, mice were sacrificed, and multiorgan metastases of ES-2 cells, such as in the liver, lung, ovary and peritoneum, were clearly observed in NCG mice. In addition, fluorescence imaging proved that the number of metastases formed by the ES-2 cell lines labeled with GFP varied significantly in different organs (Fig. 2A–C).
Tissue harvesting and fluorescence imaging of ES-2 cell line. (A) Fluorescent imaging of the metastasis formed by ES-2 cells. (B) In vivo optical imaging of liver, lung, and ovary metastases formed by ES-2 cells. (C) In vivo optical imaging of connective tissue metastases formed by ES-2 cells. The white arrowhead indicates metastases. The connective tissues in different organs are collectively referred to as Other.
From the data acquired by the SWATH strategy, a total of 6023 proteins were
identified from 32,602 different unique peptides by matching to the heterozygous
library containing human and mouse protein sequences and bovine serum albumin
peptide sequences. After excluding 17,809 mouse source peptides, the number of
remaining human peptides was 14,793, from which a total of 4503 human proteins
were identified. Ultimately, a human protein quantitative database was
established. p values were calculated using the Kruskal‒Wallis test, which were
corrected using the Benjamini‒Hochberg strategy with an FDR
To understand the gene regulation and protein profile expression changes from
circulating ES-2 cells in blood vessels to different organs in colonization and
metastasis, we first performed a hierarchical clustering analysis (HCA) on
protein data from ES-2 cell lines and the metastases of different organs. From
the heatmap (Fig. 3A), the tumors metastasized to organs showed significant
differences in protein expression profiles compared with the ES-2 cell lines.
Moreover, various profiles were also clearly observed among the metastases from
different organs. Subsequently, a principal component analysis (PCA) of the
protein expression data of different metastasis groups was carried out, and the
results showed obvious distinctions in the distribution of these metastatic
samples, each with significant organ-dependent separation. Interestingly, the
protein spectrum of liver metastases showed the most pronounced deviation
compared with metastases from other organs (Fig. 3B). Next, we performed
statistical analysis of differential proteins between the ES-2 cell line and each
organ metastatic tumor group separately. According to the screening criteria of
Multivariate statistical analysis of the total protein expressed
by different organs metastasis and ES-2 cell lines. (A) Hierarchical cluster
analysis of the 2868 DEPs. Heatmap showed that clustering of differential
proteins between ES-2 cells and different organ metastases. (B) Principal
Component Analysis of the remanent DEPs from metastases in different organs after
removing the proteins expressed by the cells. PCA analysis showed that there was
a significant outlier trend of proteins expressed by liver metastases. The
p-value was calculated using the Kruskal-Wallis test. A Benjamini-Hochberg
Since in all metastases, liver metastases showed the most significant differences in protein expression compared to original ES-2 cells and showed obvious deviation from other organ metastases in PCA scoring, we focused on the 158 proteins specifically expressed in liver metastases (Listed in Supplementary Table 2) as representatives in subsequent bioinformatical studies. We first performed GO analysis and KEGG pathway analysis to understand the protein expression changes observed in liver metastases. GO is used to describe gene functions and is divided into three categories: biological processes, cellular components, and molecular functions. GO enrichment analysis showed that 158 proteins specifically expressed in the liver were mainly located in the cytosol and nucleus and were mainly involved in DNA transcription, and molecular function focused on ATP binding (Fig. 4A). Further KEGG pathway analysis of the above 158 proteins showed that these proteins were mainly enriched in metabolic pathways (Fig. 4B).
Summary of the 158 proteins specifically expressed in liver metastases. (A) GO annotation of the 158 proteins specifically expressed in liver metastases. Only the top 10 words are shown according to the order of the p-adjusted. BP, biological process; CC, cellular component; MF, molecular function. (B) KEGG pathway enrichment analysis of 158 proteins specifically expressed by liver metastases. (C) Diagram of classical overlapping pathway analysis of 158 protein genes specifically expressed in liver metastases (Classical overlapping pathways share one or more genes by interconnection, where the deeper the red, the more significant the contribution).
Next, we introduced these protein genes into the IPA for classical overlapping canonical pathway analysis and screened genes related to the liver and ovarian cancer background as well as tumor colonization and metastasis. Pathway analysis showed that the genes specifically expressed by liver metastases were mainly concentrated in the oxidative phosphorylation and glutathione biosynthetic pathways (4C). Based on IPA pathway analysis results and SWATH quantitative data, we screened 17 proteins that were upregulated in liver metastases with functional importance and sufficient MS response intensity for further verification by targeted quantification using HR-MRM acquisition. Six proteins (Listed in Table 1), including ferritin light chain (FTL), lactate dehydrogenase A (LDHA), eukaryotic translation initiation factor 2 subunit 1 (EIF2S1), carbamoyl-phosphate synthase (CPS1), long-chain-fatty-acid–CoA ligase 1 (ACSL1) and glutamate–cysteine ligase catalytic subunit (GCLC), which could be effectively quantified by HR-MRM acquisition and showed exact upregulation in liver metastases, were screened out (Fig. 5A), and this result further narrowed the range of target proteins.
|Entry name||Protein names||Gene names|
|FRIL_HUMAN||Ferritin light chain||FTL|
|LDHA_HUMAN||Lactate Dehydrogenase A||LDHA, PIG19|
|IF2A_HUMAN||Eukaryotic translation initiation factor 2 subunit 1||EIF2S1|
|ACSL1_HUMAN||long-chain-fatty-acid–CoA ligase 1||ACSL1, FACL1, FACL2, LACS, LACS1, LACS2|
|GSH1_HUMAN||Glutamate–cysteine ligase catalytic subunit||GCLC, GLCL, GLCLC|
RT-qPCR results. (A) Exact upregulation of six proteins
quantified by HR-MRM acquisition. (B) PCR preliminary experiment of specific
primer for target protein gene (ACSL1, LDHA and FTL), negative control for normal
liver tissue of NCG mice. (C) Real-time PCR analysis of genes upregulated in
liver metastases (ACSL1, LDHA, and FTL) expression in various organ metastases of
NCG mice (n = 3). The positive control (POS) was the ES-2 cell line, and the
negative control (CON) was the liver tissue of normal NCG mice, GAPDH was used as
a reference gene. *p
To distinguish between human and mRNAs, we designed specific primers aimed at differential bases for distinguishing human- and murine-derived cDNA and performed PCR preexperiments. PCR experiments demonstrated that this series of primers are targeted to human-derived cDNA (Fig. 5B). To examine the reliability of the proteomic data in the study, six proteins that were completely upregulated in liver metastases resulting from the HR-MRM validation were selected for RT‒PCR validation. The results showed that the three differential protein gene expression levels derived from the six candidate target proteins were consistent with the proteomic data, including ACSL1, LDHA, and FTL, which were significantly higher in liver metastases than in lung, ovary, and other metastases (Fig. 5C).
Despite our deeper understanding of the process of tumor metastasis in recent years, distant metastasis remains the leading cause of cancer death . Our understanding of the complex interaction of tumor metastasis leading to a series of molecular events between tumor cells and the microenvironment is still insufficient. This process persists and is in dynamic changes, requiring the integration of large molecular data from in vivo and in vitro experimental models and the information therein to characterize the metastatic process of tumor cells. Here, we employed NCG mice, an immunodeficient mouse strain, to construct a human-derived tumor cell line xenograft model. NCG mice are a severely immunodeficient strain with defects in mature T cells, B cells, and NK cells, which makes NCG mice one of the best human-derived tumor cell line xenograft model animals to date . In preliminary experiments, we found that ES-2 cell line can present very stable liver, lung, ovary, and subcutaneous multisite metastases in NCG mice, providing an excellent source of model samples for this study. At the same time, to solve the difficulty of separating human tumors from mouse tissues, we innovatively removed mouse peptides at the level of omics data processing and successfully constructed a protein quantitative database containing only human specificity. In the choice of methodology, we used a combination of untargeted and targeted protein quantification techniques to identify and verify the target proteins. This strategy integrates the advantages of targeted and nontargeted protein quantification techniques to make the final proteomics quantification results more realistic . As a quantitative proteome technology based on a data-independent acquisition mode, SWATH has high throughput, high sensitivity and no bias and is more suitable for studies requiring the identification of large-scale unknown proteins . To further validate the differentially expressed proteins obtained based on SWATH data screening, we chose HR-MRM quantification based on high-resolution mass spectrometry to validate the target protein quantified by SWATH acquisition mode [20, 21]. It is worth mentioning that to detect the expression of human proteins specifically, quantitative strategies based on sequence specificity were used throughout this study. In SWATH analysis, we first formed a heterozygous library by combining human and mouse protein databases, and only unique peptides matched in this heterozygous library were chosen for protein quantification; thus, the coexisting peptides in human and mouse protein databases were excluded to minimize the interference between two species. This strategy made up for the disadvantage that the organ metastases obtained from the experimental operation cannot completely eliminate the mouse tissue and accurately distinguish the human and mouse proteins during the protein identification process, through which the subsequent data analysis could reflect the true regulation of protein expression in human tumor cells.
In this ES-2 metastasis model, compared with ES-2 cells cultured in vitro, the multi-organ metastatic tissues showed plenty of similar gene regulation, which was reflected by the 1020 proteins in the intersection of multiple organ metastases in Venn diagram. We analyzed 227 proteins that were consistently upregulated in various organ metastases and found a lot of proteins are involved in the cell adhesion with ATP binding and cadherin participation. For example, the upregulation of Integrin beta-4 (ITB4) and its downstream signaling molecules mitogen-activated protein kinase kinase 2 (MAPK2) confirmed that the Ras/Raf/MEK/ERK pathway associated with the tumor cell adhesion process is activated during the multiorgan metastasis of ovarian clear cell carcinoma . It has been reported that the activation of the Ras/Raf/MEK/ERK pathway mediated by the upregulation of ITB4 is a switch for cancer cells to initiate a series of sequential metastatic steps, including extracellular matrix remodeling and invasion, the active intracellular protein transport and the hydrolysis process [23, 24]. In addition, insulin-like growth factor II mRNA-binding protein 3 (IGF2BP3), was also found to be upregulated in all organ metastases . IGF2BP3 is known to be synthesized de novo in cancer, and it acts as oncogenes to stimulate the growth and migration of tumor cells in lung cancer, gastrointestinal cancer and ovarian cancer .
On the other hand, our result also reflected the obvious differential gene regulations among different organ metastases of ES-2 cells, especially, the protein expression profile of liver metastases showed a significant trend of outliers compared with that of other organ metastases. Therefore, liver metastases were chosen as typical research objects for differential protein screening in this study. Through multivariate statistical and bioinformatic analysis as well as sequence-specific validation, including HR-MRM quantification at the protein level and qRT‒PCR at the mRNA level, LDHA (lactate dehydrogenase A), ACSL1 (long-chain fatty acid CoA ligase 1) and FTL (ferritin light chain) were chosen as liver metastasis-specific genes of ES-2 cells. Notably, maybe due to the heterogeneity of tumor tissue colonizing in different mice, larger standard deviations were observed in some groups in HR-MRM assay (six biological replicates for each target protein), therefore, the expressions of target proteins were doubly verified by semi-quantitative PCR and real-time qPCR at mRNA level. These results clearly exhibited the truth of the regulation on target protein expression during liver metastasis of ES-2 cells. LDHA is known to be the A chain of LDH lactate dehydrogenase. During glycolysis, LDH (L-lactate dehydrogenase) is the key rate-limiting enzyme in the metabolic process of pyruvate to lactate [27, 28]. LDH plays a crucial role in active glycolysis in cancer cells; moreover, the high expression of LDHA as a marker of tissue damage in OCCC has been demonstrated . Furthermore, ACSL1 activates long-chain fatty acids by catalyzing the formation of long-chain fatty acyl-CoA and initiates fatty acid metabolism , and its high expression is closely related to the dysregulation of lipid metabolism and active fatty acid uptake , which is consistent with our finding that genes specifically expressed in liver metastases are concentrated in metabolic pathways such as oxidative phosphorylation and glutathione biosynthesis (Fig. 6). A major function of ferritin is the storage of iron in a soluble and nontoxic state. Previous studies have shown that the particular microenvironment of persistent oxidative stress caused by high concentrations of free iron during the progression of its disease course has been shown to be important in the malignant degeneration of ovarian epithelial cells. The resulting metabolic reprogramming and the formation of a stress microenvironment are prominent features of OCCC . It is well known that cellular metabolic reprogramming is an important marker of tumorigenesis . By ingesting large amounts of nutrients from the microenvironment, tumor cells use aerobic glycolysis (Warburg effect) to provide more macromolecular intermediates while supporting their own proliferation and invasion through amino acids and lipid anabolites [15, 34]. Interestingly, the changes in metabolite concentrations resulting from the metabolic reprogramming of tumor cells also stimulate the tumor microenvironment to secrete specific cytokines together with nutrients to maintain the malignant phenotype of cancer cells and to drive the invasion and metastasis of tumor cells through extracellular matrix remodeling . It is well known that the organ tropism of metastatic tumors largely depends on the similarity of driver mutations in the metastatic site to the primary tumor , and studies have been reported predicting the properties of bone metastasis by the gene expression profile of primary breast cancer . At the same time, similar expression trends have been reported of these three selected proteins in other cancer models, which may predict a homologous regulatory role in other multiple cancer models . Secondary distal metastasis sites create a microenvironment similar to the primary tumor to support the metastatic colonization of the tumor cells through their interaction with the tumor cells. Actually, the liver is the most important organ for glucose and lipid metabolism and iron storage. On the one hand, its active metabolic microenvironment and colonization in the tumor cell complex interaction and metabolic coupling drive the metabolic conversion of tumor cells, making the key metabolic enzymes of anabolic pathway appear highly expressed and synthesizing more nutrients to maintain the abnormal proliferation of cancer cells and invasion and metastasis . On the other hand, increased expression of ACSL1, FTL and LDHA in liver metastases reflects the adaptation of tumor cells to the liver microenvironment through metabolic reprogramming. The above results indicate that the liver has a microenvironment similar to that of OCCC due to metabolic reprogramming and oxidative stress, which might be one of the key reasons why OCCC cells tend to colonize the liver to form distal metastases.
Schematic representation of the metabolic features of the cancer
cells in liver metastases. The figure illustrates the pathways of glycolysis,
OXPHOS, and lipid metabolism in the tumor cells. Key enzymes involved in
different metabolic pathways are marked in different colors: blue-glycolysis;
orange-lipid metabolism. TCA, tricarboxylic acid cycle; OXPHOS, oxidative
phosphorylation; ATP, Adenosine triphosphate; ROS, Reactive oxygen species; PGC1,
In summary, in this research, on a xenograft mouse model of a human-derived tumor cell line, the protein expression profile study of multiorgan metastases from ES-2 cells was realized by a sequence-specific proteomic analysis strategy. Our approach provides a new idea for the study of the mechanism of tumor target organ metastasis and colonization, by which the confusion of tumor cell proteins and target organ proteins can be effectively avoided and realize the analysis of the true protein expression profile of tumor cells.
The datasets used and/or analyzed during the current study are available from the corresponding author on reasonable request.
NY and MG designed this study. MG, JC, SW, HY and WZ processed, analyzed and interpreted data. NY, XL, YL and TS contributed to the experiments. NY and MG wrote the draft. JC and MG review and edited the draft. All authors read and approved the final version of the manuscript.
Study protocol obeyed the Declaration of Helsinki for all animal experimental investigations and was approved and filed by the Ethics Committee of the West China Hospital, Sichuan University (approval number: 20220222033).
This work was supported by the Science and Technology Department of Sichuan
Province (No. 2021YFH0061) and the 1
The authors declare no conflict of interest.
Publisher’s Note: IMR Press stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.