OSBPL3 promotes colorectal cancer progression through modulation of TGF-β signaling and immune cell infiltration
Highlight box
Key findings
• Oxysterol-binding protein-like 3 (OSBPL3) is highly expressed in colorectal cancer (CRC) and is associated with increased proliferation, migration, and invasion of CRC cells.
• OSBPL3 expression correlates with TGF-β/SMAD signaling activation, epithelial-mesenchymal transition (EMT) markers, and variations in tumor immune cell infiltration.
What is known and what is new?
• OSBPL3 is implicated in cancer progression in several tumor types, but its role in CRC and interactions with the tumor microenvironment are not well understood.
• This study demonstrates that OSBPL3 promotes CRC progression through TGF-β-mediated EMT and influences the immune landscape, highlighting its potential as a therapeutic target.
What is the implication, and what should change now?
• OSBPL3 may serve as a clinically relevant biomarker for CRC aggressiveness and as a target for intervention strategies aiming to inhibit tumor progression and modulate the tumor immune microenvironment.
• Further validation in multi-center clinical cohorts and mechanistic studies are warranted to translate these findings into clinical practice.
Introduction
Colorectal cancer (CRC), as the third most common cancer, has gained increasing attention all over the world during the past decades (1). The development of CRC is closely related to unhealthy lifestyle habits such as a high-fat, low-fiber diet, lack of physical activity, obesity, smoking, and alcohol consumption. Additionally, about 10% to 15% of CRC patients have a genetic predisposition (2). The most common hereditary syndromes include hereditary nonpolyposis CRC (HNPCC), also known as Lynch syndrome (3), and familial adenomatous polyposis (FAP) (4,5). The 5-year survival rate for CRC varies significantly depending on the stage at diagnosis. In high-income countries where early detection is more common, 5-year survival rates range between 60% and 65% (6). In contrast, in low-income countries, where CRC is often diagnosed at a later stage, survival rates are much lower, typically below 40%. The incidence of CRC remains generally higher in men than in women, with men being more susceptible to lifestyle-related risk factors (such as alcohol consumption, smoking, and high-fat, low-fiber diets) (1,7).
Previous studies have shown that CRC patients with microsatellite instability-high (MSI-H) have better prognoses and respond well to immune checkpoint inhibitors and other immunotherapies (8,9). In recent years, an increasing number of CRC treatment strategies have been developed to offer personalized therapies specifically for MSI-H patient groups.
Additionally, KRAS, NRAS, and BRAF mutations are common molecular features in CRC (10,11). In recent years, testing for these mutations has increased, providing more guidance for targeted therapies. RAS mutations are often associated with resistance to anti-EGFR therapies, while BRAF mutations (especially BRAF V600E) typically indicate a poorer prognosis (12,13).
Oxysterol-binding protein-like 3 (OSBPL3) is a protein that belongs to the oxysterol-binding protein (OSBP) family, which is involved in regulating intracellular lipid transport, metabolism, and signaling (14). OSBPL3 plays an important role in lipid exchange between cell membranes and organelles, especially in the transport and metabolism of cholesterol and phospholipids (15). Emerging evidence suggests OSBPL3 exerts context-dependent functions in oncology. Studies indicate a tumor-suppressive role in specific malignancies, including colon cancer (16), lymphoma (17) and bladder cancer (18). Conversely, OSBPL3 overexpression is associated with invasiveness and disease progression in hepatic and gastric cancers (19,20). In patients with lung and gastric malignancies, elevated OSBPL3 expression correlates with metastasis, suggesting its potential utility as a prognostic marker (20,21). Despite these associations, the specific role of OSBPL3 in CRC remains relatively underexplored, and further research is needed to better understand its function in this context.
Given to these, we conduct a comprehensive study to explore the expression of OSBPL3 in patients with CRC through bioinformatic analysis and cellular experiments, detecting the role of OSBPL3 as well as related signaling pathway, especially the TGF-β signaling pathway, in CRC progression and development, providing potential therapeutic targets for clinical treatment. We present this article in accordance with the MDAR and ARRIVE reporting checklists (available at https://tcr.amegroups.com/article/view/10.21037/tcr-2025-aw-2464/rc).
Methods
Cell culture
CRC cell lines HT29 and SW480 and normal colon epithelial cell line NCM460 were purchased from the Cell Bank of the Chinese Academy of Sciences. All cells were cultured in Dulbecco’s modified Eagle’s medium (DMEM, Beijing Zhongshan Jinqiao, Beijing, China) containing 10% fetal bovine serum (FBS, Hangzhou Sijiqing, Hangzhou, China) and 1% penicillin-streptomycin solution (Shanghai Yuanye Biotechnology, Shanghai, China) in a constant temperature incubator at 37 °C and 5% CO2.
To study the function of OSBPL3, we used lentiviral vectors to construct stable knockdown (sh-OSBPL3-1, sh-OSBPL3-2) and overexpression of OSBPL3 (oe-OSBPL3) cell lines, and the negative controls were sh-NC and empty vector, respectively. The design and synthesis of lentiviral vectors were completed by Guangzhou Jisai Biotechnology, and cell infection was carried out according to the manufacturer’s instructions. Puromycin (Shanghai Kaichuang Biotechnology, Shanghai, China) was used to screen positive clones 72 hours after infection. The primers used in this study were designed based on sequences obtained from the NCBI database and synthesized commercially. The primer sequences for quantitative reverse transcription polymerase chain reaction (qRT-PCR) are provided in Table S1.
Immunohistochemistry (IHC) analysis
Human CRC tissue samples used for IHC were collected from patients at Jiangmen Central Hospital after written informed consent was obtained from all patients. The human study was approved by the Institutional Research Ethics Committee of Jiangmen Central Hospital (decision No. JXY2023110) and was conducted in accordance with the Declaration of Helsinki and its subsequent amendments.
IHC was performed according to our previous publication (18). Sections were baked at 60 °C for 1 hour, followed by xylene dewaxing and gradient ethanol hydration. Antigen retrieval was performed using citrate buffer (pH 6.0) heated in a microwave oven and treated with 3% hydrogen peroxide for 15 minutes at room temperature to block endogenous peroxidase activity.
Sections were incubated with OSBPL3 primary antibody (1:200, Thermo Fisher, Waltham, USA) overnight, and incubated with biotin-labeled secondary antibody (Abcam, Cambridge, UK) for 1 hour the next day. DAB reagent (Abcam) was used for color development, and images were observed and taken under a microscope at high and low magnification fields.
Western blot (WB) analysis
Total protein was extracted from cells using Racial and Identity Profiling Act (RIPA) lysis buffer (Solarbio, Beijing, China), and protein concentration was quantified using a bicinchoninic acid (BCA) kit (Solarbio). Equal amounts of protein samples were separated by sodium dodecyl sulfate-polyacrylamide gel electrophoresis (SDS-PAGE) and transferred to polyvinylidene difluoride (PVDF) membranes (Solarbio). The membranes were blocked in 5% skim milk powder for 1 hour at room temperature and then incubated overnight with OSBPL3 (1:1,000, Thermo Fisher Scientific, Inc.) and glyceraldehyde-3-phosphate dehydrogenase (GAPDH) (1:3,000, Thermo Fisher Scientific, Inc.) antibodies. After incubation with secondary antibodies for 1 hour, the membranes were developed using electrochemiluminescence (ECL) luminescent solution (Thermo Fisher Scientific, Inc.), and the grayscale values of the bands were analyzed using ImageJ software. The original uncropped WB images are provided in Figure S1.
Real-time quantitative PCR (RT-qPCR)
Total RNA was extracted from cells using Trizol reagent (Solarbio), and cDNA was synthesized using a reverse transcription kit (Solarbio). qPCR was detected on a Bio-Rad CFX96 system using SYBR Green Master Mix (Qingdao Patel, Qingdao, China). The relative mRNA expression level of OSBPL3 was calculated by the ΔΔCt method, and GAPDH was used as an internal reference gene.
Cell Counting Kit-8 (CCK-8) proliferation assay
Cells were seeded in 96-well plates (2,000 cells per well), and CCK-8 reagent (Solarbio) was used on the 1st, 2nd, 3rd, 4th, and 5th days after culture. The absorbance [optical density (OD) value] was measured at a wavelength of 450 nm after 2 hours of incubation. The experiment was repeated three times, with three replicates in each group.
Ethynyl deoxyuridine (EdU) incorporation assay
Using the EdU kit (Guangzhou Cell Probe, Guangzhou, China), cells were seeded in 24-well plates, and EdU solution was added for incubation for 2 hours, followed by fixation and fluorescent staining. Fluorescence microscopy (Nikon, Shanghai, China) was used to detect red fluorescently labeled EdUpositive cells to assess the DNA synthesis activity of cells.
Flow cytometry for apoptosis
Cells were labeled using Annexin V-FITC/PI double staining kit (Nanjing Keygen Biotechnology, Nanjing, China) and analyzed by flow cytometry (BD Biosciences, Franklin Lakes, USA). At least 10,000 events were detected for each group, and data were analyzed using FlowJo software.
Wound healing assay
Cells were seeded in 6-well plates. When the cell confluence reached more than 90%, a scratch was made on the cell monolayer with a sterile pipette tip. After 24 hours of culture, images of the scratch area were taken using a microscope (Olympus, Shanghai, China), and the wound closure rate was measured using ImageJ to evaluate cell migration ability.
Transwell invasion assay
The invasion assay used a 24-well Transwell chamber (8.0 µm pore size, Corning, Shanghai, China); 100 µL of Matrigel matrix gel (BD Biosciences) was added to the chamber and incubated for 30 minutes to solidify. The treated cells were seeded on the upper layer of the chamber and cultured in medium containing 10% FBS for 24 hours. The cells that passed through the membrane in the lower layer were fixed with methanol, stained with 0.1% crystal violet, photographed and counted under a microscope.
Gene set enrichment analysis (GSEA)
In order to explore the molecular mechanism of OSBPL3 in CRC, we downloaded the RNA-seq expression matrix and corresponding clinical data from The Cancer Genome Atlas (TCGA) CRC projects, including COAD and READ, using the TCGAbiolinks package in the R programming environment, focusing on the relationship between OSBPL3 and signaling pathways. According to the high and low expression groups of OSBPL3, the gene expression profile was analyzed to determine the relevant Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways, especially the association with the TGF-β signaling pathway. All bioinformatic analyses were completed in the R language environment, and the ClusterProfiler package was used for pathway enrichment analysis and visualization of the results.
Bioinformatic analysis
CRC and pan-cancer RNA sequencing data were obtained from the TCGA (https://www.cancer.gov/tcga) and Genotype-Tissue Expression (GTEx) database (https://gtexportal.org/), and the “limma” package was used for differentially expressed genes (DEGs) screening. The screening conditions were |log2 fold change (FC)| >1 and p.adj value <0.05 to determine the significantly DEGs between high and low expression of OSBPL3. The R package corrplot was used to perform inter-gene co-expression network analysis to further explore the potential regulatory mechanism of OSBPL3. Immune infiltration analysis was performed using the CIBERSORTxR package. Protein interaction networks were completed using the STRING website (https://string-db.org/).
Nude mouse xenograft model
All animal experiments complied with the 2020 AVMA Guidelines for the Euthanasia of Animals and were approved by the Animal Ethics Committee of Guangdong Medical University (approval No. GDY2302537), in compliance with the institutional guidelines for the care and use of animals. Female BALB/c nude mice (4–6 weeks; Beijing Weitong Lihua, Beijing, China) were housed in a specific pathogen-free facility at the Laboratory Animal Center of Guangdong Medical University in individually ventilated cages (3–5 mice per cage) with corncob bedding, sterilized chow and water ad libitum, a 12-h light/12-h dark cycle, ambient temperature 22±2 °C, and relative humidity 45–65%; environmental enrichment (nesting material, shelters, and gnawing blocks) was provided, and mice were socially housed unless welfare or scientific considerations required temporary single housing. Animals were acclimatized for at least 7 days and checked at least once daily by trained staff. Mice were randomized to groups, and tumor measurements were performed by investigators blinded to allocation. For inoculation, control or OSBPL3-knockdown HT29 or SW480 cells (2×106 per mouse) were injected subcutaneously under light isoflurane anesthesia (2–3% in oxygen via nose cone). Tumor length and width were measured at least 2–3 times per week with a caliper, and volume was calculated as V = 0.5 × length × width2. At week 4 or upon reaching predefined humane endpoints, mice were euthanized by isoflurane overdose (≥5% in oxygen in an induction chamber); exposure was maintained for at least 1 min after cessation of respiration, followed by cervical dislocation to ensure death. Death was confirmed by the absence of a heartbeat and corneal reflexes before tissue collection. Ether, chloroform, and chloral hydrate were not used at any time.
Statistical analysis
All experiments were repeated at least three times. Data are presented as mean ± standard deviation (SD). Statistical analyses were performed using GraphPad Prism 8.0 (GraphPad Software, San Diego, CA, USA). Comparisons between two groups were conducted using Student’s t-test, and comparisons among multiple groups were performed using one-way analysis of variance (ANOVA). A P value <0.05 was considered statistically significant.
Results
Expression and clinical significance of OSBPL3 in CRC
We first examined the expression of OSBPL3 in tumors. Pan-cancer analysis of unmatched samples based on the TCGA database revealed that OSBPL3 was differentially expressed across various cancers, with significant differences observed in COAD (Figure 1A). Analysis of matched samples further confirmed that OSBPL3 expression levels were markedly different in COAD and READ (Figure 1B), and survival analysis indicated that high OSBPL3 expression was significantly associated with poor prognosis in COAD patients (Figure 1C). Similar protein expression patterns were observed in COAD in the unmatched and matched analyses (Figure 1D,1E). Additionally, receiver operating characteristic (ROC) curve analysis demonstrated that when OSBPL3 was used as a diagnostic marker, the area under the curve (AUC) reached 0.979 (Figure 1F), highlighting its excellent diagnostic accuracy. This suggests that OSBPL3 could effectively distinguish CRC patients from healthy individuals and holds potential as a biomarker for early diagnosis and monitoring of CRC. To further evaluate the clinical relevance of OSBPL3, we conducted analyses based on clinicopathological parameters (Figure 1G-1N). The results showed that OSBPL3 expression varied significantly (P<0.050) across different pathological N stages (N0 vs. N1&N2), indicating a potential association with lymph node metastasis (Figure 1K). Additionally, OSBPL3 expression differences between pathological M stages (M0 vs. M1) approached statistical significance (P=0.06), suggesting a possible role in metastatic CRC (Figure 1L). The remaining parameters, including age, body mass index, gender, pathologic T stage, histologic type, and carcinoembryonic antigen (CEA) level, did not show significant differences (Figure 1G-1J,1M,1N). In addition, by dividing the patients into high and low expression groups according to the expression level of OSBPL3, subgroup analysis of clinical indicators showed that there were statistical differences in N stage (P=0.05) and pathologic 214 stage (P=0.04) among different expression groups (Table 1).
Table 1
| Characteristics | Low expression of OSBPL3 (n=239) | High expression of OSBPL3 (n=239) | P value |
|---|---|---|---|
| Age, years | 0.71 | ||
| ≤65 | 99 (20.7) | 95 (19.9) | |
| >65 | 140 (29.3) | 144 (30.1) | |
| Gender | 0.36 | ||
| Female | 118(24.7) | 108 (22.6) | |
| Male | 121 (25.3) | 131 (27.4) | |
| BMI, kg/m2 | 0.93 | ||
| ≤25 | 35 (13.7) | 52 (20.3) | |
| >25 | 69 (27) | 100 (39.1) | |
| CEA level, ng/mL | 0.53 | ||
| ≤5 | 99 (32.7) | 97 (32.0) | |
| >5 | 50 (16.5) | 57 (18.8) | |
| OS event | 0.22 | ||
| Alive | 193 (40.4) | 182 (38.1) | |
| Dead | 46 (9.6) | 57 (11.9) | |
| DSS event | 0.40 | ||
| No | 203 (43.9) | 195 (42.2) | |
| Yes | 29 (6.3) | 35 (7.6) | |
| Pathologic T stage | 0.15 | ||
| T1&T2 | 47 (9.9) | 47 (9.9) | |
| T3 | 168 (35.2) | 155 (32.5) | |
| T4 | 23 (4.8) | 37 (7.8) | |
| Pathologic N stage | 0.047 | ||
| N0 | 154 (32.2) | 130 (27.2) | |
| N1 | 51 (10.7) | 57 (11.9) | |
| N2 | 34 (7.1) | 52 (10.9) | |
| Pathologic M stage | 0.07 | ||
| M0 | 185 (44.6) | 164 (39.5) | |
| M1 | 27 (6.5) | 39 (9.4) | |
| Pathologic stage | 0.040 | ||
| Stage I | 39 (8.4) | 42 (9.0) | |
| Stage II | 109 (23.3) | 78 (16.7) | |
| Stage III | 61 (13.1) | 72 (15.4) | |
| Stage IV | 27 (5.8) | 39 (8.4) |
Data are presented as n (%). BMI, body mass index; CEA, carcinoembryonic antigen; DSS, disease-specific survival; M, metastasis; N, node; OS, overall survival; T, tumor.
Comprehensive analysis of DEGs in CRC and its clinical significance
To elucidate the comprehensive differences in gene expression in CRC, we analyzed data from the TCGA database and identified significantly upregulated and downregulated DEGs, including a total of 4,694 upregulated and 1,032 downregulated genes (logFC >1, p.adj <0.05) (Figure 2A). Given that tumor progression often involves complex regulatory networks, we performed a co-expression analysis of OSBPL3 with these DEGs (Figure 2B-2F). The results revealed a strong positive correlation between OSBPL3 and ZBTB20, and a strong negative correlation with SNHG25. Gene Ontology (GO) analysis indicated that these genes are mainly involved in biological processes (BPs) such as chemical stimulus detection, oxygen transport, and cellular oxide detoxification (BP, Figure 2G), and are enriched in cellular components (CCs) like hemoglobin complexes and DNA packaging complexes (CC, Figure 2H). Molecular function (MF) analysis demonstrated associations with oxygen binding, antioxidant activity, and related functions (Figure 2I). KEGG pathway analysis revealed significant enrichment of these genes in pathways such as taste transduction, spliceosome, and systemic lupus erythematosus (Figure 2J).
To further explore the expression patterns of these genes across different clinicopathological characteristics, we conducted a co-expression subgroup analysis. The results showed that OSBPL3 was closely associated with multiple DEGs in the T2&3, N0, and M0 stages (Figure 2K-2M). Additionally, OSBPL3 exhibited stronger associations with differential genes in the survival group compared to the non-survival group (Figure 2N). Clustering analysis of gene expression patterns with overall survival (OS) and tumor-node-metastasis (TNM) staging indicated that key genes associated with OSBPL3 were significantly correlated with patient survival status and TNM pathological stages (Figure 2O,2P). These findings suggest that these genes may play critical roles in CRC staging and progression.
Correlation analysis between OSBPL3 expression and immune cell infiltration in tumor microenvironment (TME)
The TME is widely recognized as an important contributor to tumor progression. To examine the relationship between OSBPL3 expression and immune infiltration, we performed an in silico pan-cancer analysis based on CIBERSORT-derived immune cell estimates. This analysis suggested that OSBPL3 expression is associated with the inferred infiltration levels of multiple immune cell subsets across cancers (Figure 3A). In CRC, Spearman correlation analysis identified statistically significant associations between OSBPL3 expression and several immune cell types, including central memory T cells (Tcm) (r=0.478, P<0.001), Th2 cells (r=0.140, P=0.002), and Treg cells (r=−0.196, P<0.001) (Figure 3B). Notably, while the correlations with Th2 and Treg cells were statistically significant, their effect sizes were weak, indicating limited biological magnitude and warranting cautious interpretation.
To further characterize these associations, CRC samples were stratified into high and low OSBPL3 expression groups. Differences in estimated infiltration were observed for several immune subsets, including Treg cells, Th17 cells, and Tcm cells (Figure 3C,3D), supporting that OSBPL3 expression is linked to variations in the immune infiltration landscape. To explore potential molecular features related to these patterns, we assessed the co-expression relationships between OSBPL3-associated DEGs and immune cell subsets. Overall, OSBPL3 expression tended to show negative associations with multiple inferred immune populations, consistent with immune-excluded or immunosuppressive TME features observed in some CRC contexts (Figure 3E,3F). We also observed positive associations with Th cells (predominantly Th2), macrophages, and Tcm cells. Importantly, these findings are correlative and computationally inferred; they do not establish functional immune regulation and should be interpreted as hypotheses that require experimental validation.
Interaction and correlation analysis of OSBPL3 with key proteins in CRC
STRING database analysis revealed that OSBPL3 interacts with multiple key proteins in CRC, including NFE2L2, VAPA, GAK, and HNRNPA2B1, suggesting that OSBPL3 may participate in various BPs through these protein interactions (Figure 4A). To further clarify the association between OSBPL3 and these related proteins, we performed Spearman correlation analysis on the expression levels of OSBPL3 and several related proteins. The results showed that OSBPL3 expression was significantly correlated with the expression of multiple proteins, such as VAPB (r=0.377, P<0.001, Figure 4B), VAPA (r=0.356, P<0.001, Figure 4C), RRAS (r=−0.090, P=0.054, Figure 4D), NFE2L3 (r=0.334, P<0.001, Figure 4E), NFE2L2 (r=0.481, P<0.001, Figure 4F), HNRNPA3 (r=0.417, P<0.001, Figure 4G), and HNRNPA2B1 (r=0.518, P<0.001, Figure 4H). These correlations suggest that OSBPL3 may influence the biological behavior of tumor cells by regulating the expression of these genes. However, the analysis also found that certain genes, such as GAK (r=0.046, P=0.33, Figure 4I) and PDE11A (r=−0.090, P=0.05, Figure 4J), did not show significant correlation with OSBPL3 expression, suggesting that the role of OSBPL3 in CRC may not be mediated through these specific molecular mechanisms.
In summary, the significant positive correlations between OSBPL3 and several key proteins in CRC indicate that OSBPL3 may have a crucial regulatory role in the biological behavior of tumor cells. However, its mechanisms may rely on specific molecular pathways, warranting further investigation.
OSBPL3 is highly expressed in CRC
To investigate the expression of OSBPL3 in CRC, we performed IHC and WB analyses. Using data from the Human Protein Atlas (HPA) database (https://www.proteinatlas.org/about/licence), the IHC results revealed low OSBPL3 expression in normal colon tissues but significantly elevated expression in CRC tissues. Under both high-magnification (10×10) and low-magnification (4×10) views, CRC samples exhibited markedly stronger positive staining for OSBPL3, confirming its high expression in these tissues (Figure 5A). Additionally, WB and PCR analyses of various CRC cell lines demonstrated that OSBPL3 protein levels (Figure 5B) and mRNA (Figure 5C) were significantly higher in CRC cell lines (SW620, SW480, HT29, HT115) compared to the normal colon cell line (NCM460).
OSBPL3 promotes proliferation, migration and invasion and inhibits apoptosis in CRC cells
To comprehensively assess the function of OSBPL3 in CRC cells, we conducted knockdown and overexpression experiments to examine its effects on cell proliferation, apoptosis, migration, and invasion (Figures 6,7). CCK-8 assays showed that OSBPL3 knockdown significantly inhibited the proliferation of HT29 and SW480 cells (Figure 6A, P<0.001), whereas OSBPL3 overexpression markedly enhanced their proliferation (Figure 7A, P<0.001). These findings were further supported by EdU incorporation assays, which demonstrated a significant reduction in the proportion of EdU-positive cells following OSBPL3 knockdown (Figure 6B, P<0.001) and a significant increase after OSBPL3 overexpression (Figure 7B, P<0.001). Flow cytometry analysis using Annexin V/PI staining showed that OSBPL3 knockdown was associated with an increased apoptotic fraction in HT29 and SW480 cells (Figure 6C), while OSBPL3 overexpression showed the opposite trend (Figure 7C). The bar graphs report the total apoptosis rate, defined as the sum of early apoptosis (Annexin V+/PI−) and late apoptosis (Annexin V+/PI+), with the corresponding quadrants and experimental conditions clearly labeled in the revised figures. Although the differences reached statistical significance, the absolute change was modest (≤~5%); therefore, the biological relevance of this effect should be interpreted cautiously and may warrant further validation. Wound healing assays demonstrated that OSBPL3 knockdown significantly impaired the migration ability of HT29 and SW480 cells (Figure 6D, P<0.001), whereas OSBPL3 overexpression significantly enhanced their migration (Figure 7D, P<0.001). Similarly, Transwell invasion assays showed that the invasive capacity of HT29 and SW480 cells was significantly diminished after OSBPL3 knockdown (Figure 6E, P<0.001), whereas OSBPL3 overexpression led to a marked increase in invasion (Figure 7E, P<0.001). These results suggest that OSBPL3 plays a crucial positive regulatory role in CRC cell proliferation, migration, and invasion, and may facilitate tumor cell survival by inhibiting apoptosis.
OSBPL3 is associated with CRC cell behavior and correlates with TGF-β/SMAD signaling and epithelial-mesenchymal transition (EMT)-related changes
To explore the relationship between OSBPL3 and CRC cell behavior, we performed in vitro knockdown and overexpression experiments. Stable cell lines with reduced or elevated OSBPL3 expression were generated and confirmed by PCR and WB analyses (Figure 8A,8B). GSEA suggested that OSBPL3 expression is enriched in gene signatures related to the TGF-β signaling pathway, indicating a potential association with TGF-β-related downstream signaling in CRC (Figure 8C). Consistent with this association, WB analysis showed that OSBPL3 knockdown was accompanied by decreased levels of TGF-β, TGF-βR1, and p-SMAD2, together with increased E-cadherin and reduced N-cadherin expression. In contrast, OSBPL3 overexpression was accompanied by increased expression of TGF-β pathway components and corresponding changes in EMT markers (Figure 8D). Taken together, these results support a close relationship between OSBPL3 expression and TGF-β/SMAD pathway activity, as well as EMT-associated molecular alterations, in CRC cells.
OSBPL3 knockdown inhibits tumor growth of CRC cells in vivo
To evaluate the effect of OSBPL3 on CRC tumor growth, we performed in vivo experiments in a nude mouse model. OSBPL3 knockdown CRC cells and untreated control CRC cells were injected into nude mice respectively, and the tumor volume was monitored.
The experimental results showed that knocking down OSBPL3 significantly inhibited the tumor growth of CRC cells compared with the control group (Figure 9A). Statistical analysis of tumor volume further showed that the tumor volume of the sh-OSBPL3 group was significantly smaller than that of the control group (P<0.001) (Figure 9B). These results indicate that OSBPL3 plays an important role in promoting the growth of CRC tumors. As summarized in Figure 10, our data suggest that OSBPL3 is linked to TGF-β/SMAD signaling activity and EMT-related molecular changes, which may contribute to CRC invasion and metastasis.
Discussion
OSBPL3, a member of the OSBP family, plays a crucial role in regulating intracellular lipid transport, metabolism, and signaling. Previous studies have indicated that OSBPL3 expression is closely associated with the invasiveness and progression of certain cancer types, though the specific mechanisms may vary across different cancers.
Evidence suggests that OSBPL3 exerts distinct, and sometimes opposing, effects across tumor types. For example, it has been reported to function as a tumor suppressor in lymphoma and bladder cancer, whereas in gastric and liver cancers it may facilitate tumor development. These discrepancies may reflect differences in tumor type and stage, molecular background, detection platforms, cohort composition, and microenvironmental influences, including context-specific shifts in TGF-β signaling. Notably, OSBPL3 is upregulated in gastric cancer and is associated with poor prognosis (20). In contrast, OSBPL3 is expressed at low levels in hepatocellular carcinoma and pancreatic cancer, and reduced expression has been linked to altered cholesterol and lipid metabolism that may influence hepatocellular carcinoma growth (19,22,23). These findings suggest that OSBPL3 exhibits diverse functions and regulatory mechanisms across different types of cancer. In this study, we performed pan-cancer analyses to explore OSBPL3 expression in various tumor cells and adjacent normal tissues. We found that OSBPL3 expression was significantly elevated in CRC patients, suggesting that its overexpression may be closely linked to the onset and progression of CRC. To further investigate the role of OSBPL3 in CRC, we conducted subgroup analyses using data from the TCGA database, which revealed that OSBPL3 expression was associated with the pathological N stage in CRC patients. Additionally, our analysis showed a correlation between OSBPL3 expression and tumor lymph node metastasis, as well as more advanced tumor stages.
Then we conducted the co-expression analysis in patients with CRC to explore relative DEGs. Database results showed that OSBPL3 was positively related to ZBTB20 significantly and negatively associated with SNHG25 on the opposite. ZBTB20 (zinc finger and BTB domain containing 20) is a transcription factor belonging to the zinc finger protein family, containing the BTB/POZ domain, regulating gluconeogenesis and insulin signaling pathways, maintaining normal blood glucose levels (24,25). The dysfunction of ZBTB20 is associated with metabolic diseases such as diabetes (26). Researchers have found that it can act as either an oncogene or a tumor suppressor, depending on the type of tumor and the environment (27). ZBTB20 modulates key cancer-related signaling pathways, such as the PI3K/AKT, Wnt/β-catenin, and TGF-β/SMAD pathways, thereby regulating the growth, differentiation, and metastasis of tumor cells (28). Both ZBTB20 and OSBPL3 are associated with the PI3K/AKT pathway, a critical oncogenic signaling pathway in various cancers. ZBTB20 activates this pathway to promote tumor cell proliferation and survival, while OSBPL3 may indirectly influence it by regulating cell migration and lipid metabolism. Consequently, these two proteins may overlap functionally within the PI3K/AKT pathway, collaboratively promoting tumor progression. In addition, SNHG25 (small nucleolar RNA host gene 25), a long non-coding RNA (lncRNA), is abnormally expressed in various tumors and may play a significant role in tumor development and progression (19,29). In hepatocellular carcinoma, elevated SNHG25 expression is associated with enhanced tumor proliferation, migration, and invasion (30). SNHG25 promotes cancer cell growth by acting as a “sponge” to sequester tumor-suppressive miRNAs, such as miR-497 (31,32). However, there is currently no direct research or clear evidence supporting a specific interaction between SNHG25 and OSBPL3, and further investigation is needed.
Immune infiltration is an important component of the TME and is closely linked to tumor progression (33). Based on CIBERSORT-derived estimates, Spearman correlation analysis suggested that OSBPL3 expression was positively associated with Tcm and Th2 cells, and negatively associated with Th17 cells and regulatory T cells (Tregs). Importantly, although the correlations with Th2 (r=0.140, P=0.002) and Treg cells (r=−0.196, P<0.001) reached statistical significance, their effect sizes were small, indicating that OSBPL3 is unlikely to be a dominant determinant of these immune subsets on its own. Instead, the biological relevance may lie in the directionality of these shifts and their consistency with broader immune-context features of CRC: even modest but reproducible changes in Th2- or Treg-associated signals can reflect a tendency toward immune polarization or immunoregulatory tone when integrated with other TME components and tumor-intrinsic programs. Tcm represent a long-lived T-cell subset that can contribute to anti-tumor immunity by sustaining immune memory and supporting rapid responses to tumor antigens (29,30). In CRC, the MSI-H subtype is often characterized by increased T-cell infiltration and has been associated with favorable responses to immunotherapy, and Tcm enrichment has been reported in this context (31). However, immune-suppressive features within the TME, including the presence of Tregs, may attenuate effective anti-tumor activity even when memory T-cell subsets are abundant. Th2 cells secrete cytokines such as IL-4, IL-5, and IL-13 and can shape immune polarization and macrophage function; in several tumor settings, Th2-skewed responses have been linked to immunosuppressive phenotypes and immune evasion (32). Th17 cells produce pro-inflammatory cytokines such as IL-17, and their roles in cancer are context dependent, with evidence supporting both anti-tumor and pro-tumor effects through inflammation and angiogenesis (34-36). For example, Liu et al. reported that Th17-related signaling may promote EMT and enhance migration and invasion in CRC through cytokine-mediated interactions involving Fn14 (37). Tregs are an immunosuppressive T-cell subset that can limit anti-tumor immunity via inhibitory cytokines including IL-10 and TGF-β (38,39), and elevated Treg infiltration has been associated with immune evasion and unfavorable outcomes in multiple cancers (40,41). Overall, these results support an association between OSBPL3 expression and variability in the inferred immune infiltration profile in CRC, with the strongest relationship observed for Tcm. The weaker Th2 and Treg associations should be viewed as subtle trends that may contribute to, rather than define, the immune milieu, and they are best interpreted alongside other immune populations and tumor-intrinsic pathways. Since these estimates were obtained through computational deconvolution and the analyses are correlational, they provide descriptive evidence of immune-context differences rather than direct proof of functional immune regulation.
Furthermore, we verified the role of OSBPL3 in CRC through cellular experiments and animal model. IHC staining showed the expression of OSBPL3 in tumor samples was significantly higher than that in normal tissues. In vitro knockdown of OSBPL3 in HT29 and SW480 cell lines resulted in a significant reduction in proliferation. Flow cytometry analysis revealed that both early and late apoptosis ratios in the experimental group were significantly higher compared to the control group. The wound healing assay demonstrated that the experimental group showed a marked reduction in cell healing ability after 24 hours, and the Transwell assay indicated a decrease in cell invasion ability in the OSBPL3 knockdown group compared to the control group. These findings suggest that in vitro knockdown of OSBPL3 in CRC cells can inhibit cell proliferation, invasion, and migration, while promoting tumor cell apoptosis. Further studies using a tumor-bearing animal model showed that after OSBPL3 knockdown, the subcutaneous tumor volume in nude mice was significantly reduced, confirming the results of the cell experiments. In vitro overexpression of OSBPL3 in CRC cells exhibited the opposite results, promoting cell proliferation, invasion, and migration, while reducing tumor cell apoptosis.
KEGG pathway analysis suggested that OSBPL3 is enriched in the TGF-β signaling pathway, and our WB results showed that changes in OSBPL3 expression were accompanied by altered levels of TGF-β pathway-related proteins and EMT markers. TGF-β signaling has well-recognized stage-dependent roles in cancer. In early tumorigenesis, TGF-β can restrain tumor development by inhibiting cell-cycle progression, promoting apoptosis, and limiting inflammatory responses. In contrast, at later stages, TGF-β signaling is frequently associated with tumor progression, invasion, and metastasis, partly through EMT induction and immunosuppressive effects (42,43). This functional shift reflects a transition of TGF-β signaling from predominantly tumor-suppressive to tumor-promoting activity in advanced disease (44). In addition, TGF-β can suppress anti-tumor immunity by inhibiting T-cell proliferation and activation and by promoting Treg differentiation (45). EMT facilitates tumor dissemination by enabling cancer cells to breach the basement membrane, invade surrounding tissues, and enter blood or lymphatic vessels, thereby promoting metastasis to distant organs. EMT is also linked to increased resistance to apoptosis and chemotherapy, complicating treatment (46), and may endow tumor cells with stem cell-like properties that enhance self-renewal and drug resistance (47). Mechanistically, TGF-β is a key inducer of EMT and can regulate EMT-related genes via both SMAD-dependent and SMAD-independent pathways (11,43,46). Other pathways, such as Wnt/β-catenin signaling, also contribute to EMT programs by promoting β-catenin nuclear translocation and inducing EMT-related gene expression (48). Collectively, these observations support the clinical relevance of EMT as a process associated with metastasis and therapy resistance; accordingly, EMT inhibition and EMT-related markers have potential value for therapeutic development and prognostic assessment.
In summary, OSBPL3, a pivotal regulator of lipid transport and a contributor to cancer progression, is increasingly recognized for its role in CRC. Elevated OSBPL3 expression is linked to heightened tumor aggressiveness, enhanced metastatic potential, and advanced disease stages. It positively correlates with the oncogene ZBTB20 and negatively with lncRNA SNHG25, thereby modulating the PI3K/AKT pathway and influencing tumor progression. OSBPL3 also interacts with immune cells, including Tcm, Th2, Th17, and Tregs, shaping the TME and regulating EMT. Functional assays reveal that OSBPL3 knockdown reduces CRC cell proliferation, migration, and invasion, while simultaneously inducing apoptosis in vitro and in vivo; conversely, its overexpression accelerates tumor growth. Furthermore, OSBPL3 is associated with the TGF-β signaling pathway and EMT—pathways integral to metastasis and contributors to drug resistance. Nonetheless, the study is constrained by limitations such as a small sample size, bioinformatic biases, and the necessity of further validating OSBPL3’s clinical implications and therapeutic potential.
There also remains several limitations about or study. First, because this was a single-center study with a relatively small sample size, the generalizability of our findings may be limited; validation in larger, multi-center cohorts is therefore needed to confirm the clinical relevance of OSBPL3 in CRC. Second, much of our analysis relied on public databases and bioinformatic pipelines, and cross-dataset heterogeneity—such as differences in platforms, batch effects, and cohort composition—may introduce systematic bias. Third, although we observed associations between OSBPL3 expression and TGF-β/SMAD signaling activity as well as EMT-related molecular changes, the evidence remains primarily correlative and does not establish a definitive mechanistic link. Finally, immune infiltration patterns were computationally inferred (e.g., using CIBERSORT-based deconvolution) and were not validated by experimental methods such as IHC, flow cytometry, cytokine profiling, or functional immune assays. In addition, the potential relationship between OSBPL3 and drug sensitivity, as well as its value for prognostic stratification, requires further confirmation in well-designed translational studies.
Conclusions
In summary, OSBPL3 is highly expressed in CRC and contributes to tumor cell proliferation, migration and invasion. Its expression is closely associated with activation of the TGF-β/SMAD signaling pathway, EMT-related molecular changes, and modulation of the tumor immune microenvironment. These findings highlight OSBPL3 as a potential biomarker for CRC aggressiveness and a promising therapeutic target. Future studies in larger, multi-center cohorts and functional investigations are warranted to further validate its clinical significance and therapeutic potential.
Acknowledgments
All authors would like to express their sincere thanks to the participants who provided colorectal tissue samples for this study and ensuring the successful completion of the study.
Footnote
Reporting Checklist: The authors have completed the MDAR and ARRIVE reporting checklists. Available at https://tcr.amegroups.com/article/view/10.21037/tcr-2025-aw-2464/rc
Data Sharing Statement: Available at https://tcr.amegroups.com/article/view/10.21037/tcr-2025-aw-2464/dss
Peer Review File: Available at https://tcr.amegroups.com/article/view/10.21037/tcr-2025-aw-2464/prf
Funding: This work was supported by grants from
Conflicts of Interest: All authors have completed the ICMJE uniform disclosure form (available at https://tcr.amegroups.com/article/view/10.21037/tcr-2025-aw-2464/coif). The authors have no conflicts of interest to declare.
Ethical Statement: The authors are accountable for all aspects of the work in ensuring that questions related to the accuracy or integrity of any part of the work are appropriately investigated and resolved. This study was conducted in accordance with the Declaration of Helsinki and its subsequent amendments. The study was approved by the Institutional Research Ethics Committee of Jiangmen Central Hospital (decision No. JXY2023110) and informed consent was obtained from all patients. Animal experiments were performed under a project license (No. GDY2302537) granted by the Animal Ethics Committee of Guangdong Medical University, in compliance with the institutional guidelines for the care and use of animals.
Open Access Statement: This is an Open Access article distributed in accordance with the Creative Commons Attribution-NonCommercial-NoDerivs 4.0 International License (CC BY-NC-ND 4.0), which permits the non-commercial replication and distribution of the article with the strict proviso that no changes or edits are made and the original work is properly cited (including links to both the formal publication through the relevant DOI and the license). See: https://creativecommons.org/licenses/by-nc-nd/4.0/.
References
- Dekker E, Tanis PJ, Vleugels JLA, et al. Colorectal cancer. Lancet 2019;394:1467-80. [Crossref] [PubMed]
- Rawla P, Sunkara T, Barsouk A. Epidemiology of colorectal cancer: incidence, mortality, survival, and risk factors. Prz Gastroenterol 2019;14:89-103. [Crossref] [PubMed]
- Win AK, Jenkins MA, Dowty JG, et al. Prevalence and Penetrance of Major Genes and Polygenes for Colorectal Cancer. Cancer Epidemiol Biomarkers Prev 2017;26:404-12. [Crossref] [PubMed]
- Thanikachalam K, Khan G. Colorectal Cancer and Nutrition. Nutrients 2019;11:164. [Crossref] [PubMed]
- Patel SG, Ahnen DJ. Familial colon cancer syndromes: an update of a rapidly evolving field. Curr Gastroenterol Rep 2012;14:428-38. [Crossref] [PubMed]
- Sung H, Ferlay J, Siegel RL, et al. Global Cancer Statistics 2020: GLOBOCAN Estimates of Incidence and Mortality Worldwide for 36 Cancers in 185 Countries. CA Cancer J Clin 2021;71:209-49. [Crossref] [PubMed]
- Graff RE, Möller S, Passarelli MN, et al. Familial Risk and Heritability of Colorectal Cancer in the Nordic Twin Study of Cancer. Clin Gastroenterol Hepatol 2017;15:1256-64. [Crossref] [PubMed]
- Albuquerque C, Bakker ER, van Veelen W, et al. Colorectal cancers choosing sides. Biochim Biophys Acta 2011;1816:219-31. [PubMed]
- Willis JA, Reyes-Uribe L, Chang K, et al. Immune Activation in Mismatch Repair-Deficient Carcinogenesis: More Than Just Mutational Rate. Clin Cancer Res 2020;26:11-7. [Crossref] [PubMed]
- Liao W, Overman MJ, Boutin AT, et al. KRAS-IRF2 Axis Drives Immune Suppression and Immune Therapy Resistance in Colorectal Cancer. Cancer Cell 2019;35:559-572.e7. [Crossref] [PubMed]
- Li J, Ma X, Chakravarti D, et al. Genetic and biological hallmarks of colorectal cancer. Genes Dev 2021;35:787-820. [Crossref] [PubMed]
- Vaughn CP, Zobell SD, Furtado LV, et al. Frequency of KRAS, BRAF, and NRAS mutations in colorectal cancer. Genes Chromosomes Cancer 2011;50:307-12. [Crossref] [PubMed]
- Kopetz S, Grothey A, Yaeger R, et al. Encorafenib, Binimetinib, and Cetuximab in BRAF V600E-Mutated Colorectal Cancer. N Engl J Med 2019;381:1632-43. [Crossref] [PubMed]
- Olkkonen VM. OSBP-Related Protein Family in Lipid Transport Over Membrane Contact Sites. Lipid Insights 2015;8:1-9. [Crossref] [PubMed]
- Aibara D, Sakaguchi A, Matsusue K. Oxysterol-binding protein-like 3 is a novel target gene of peroxisome proliferator-activated receptor γ in fatty liver disease. Mol Cell Endocrinol 2023;565:111887. [Crossref] [PubMed]
- Xu P, Richter J, Blatz A, et al. Downregulation of ORP3 Correlates with Reduced Survival of Colon Cancer Patients with Advanced Nodal Metastasis and of Female Patients with Grade 3 Colon Cancer. Int J Mol Sci 2020;21:5894. [Crossref] [PubMed]
- Njeru SN, Kraus J, Meena JK, et al. Aneuploidy-inducing gene knockdowns overlap with cancer mutations and identify Orp3 as a B-cell lymphoma suppressor. Oncogene 2020;39:1445-65. [Crossref] [PubMed]
- Wang X, Liu J, Azoitei A, et al. Loss of ORP3 induces aneuploidy and promotes bladder cancer cell invasion through deregulated microtubule and actin dynamics. Cell Mol Life Sci 2023;80:299. [Crossref] [PubMed]
- Tian K, Ying Y, Huang J, et al. The expression, immune infiltration, prognosis, and experimental validation of OSBPL family genes in liver cancer. BMC Cancer 2023;23:244. [Crossref] [PubMed]
- Hu Q, Masuda T, Koike K, et al. Oxysterol binding protein-like 3 (OSBPL3) is a novel driver gene that promotes tumor growth in part through R-Ras/Akt signaling in gastric cancer. Sci Rep 2021;11:19178. [Crossref] [PubMed]
- Zhang Q, Zhang J, Li J, et al. OSBPL3 promotes LUAD metastasis and glycolysis through transcriptional upregulation of PLAU via NFE2L2. Transl Oncol 2026;66:102693. [Crossref] [PubMed]
- Sun Q, Zhu X, Zou Q, et al. OSBPL3 modulates the immunosuppressive microenvironment and predicts therapeutic outcomes in pancreatic cancer. Biol Direct 2025;20:5. [Crossref] [PubMed]
- Su Y, Xue C, Gu X, et al. Integrated bioinformatics analysis reveals the function and prognostic value of OSBPL3 in hepatocellular carcinoma. Heliyon 2023;9:e17223. [Crossref] [PubMed]
- Zhang W, Mi J, Li N, et al. Identification and characterization of DPZF, a novel human BTB/POZ zinc finger protein sharing homology to BCL-6. Biochem Biophys Res Commun 2001;282:1067-73. [Crossref] [PubMed]
- Mitchelmore C, Kjaerulff KM, Pedersen HC, et al. Characterization of two novel nuclear BTB/POZ domain zinc finger isoforms. Association with differentiation of hippocampal neurons, cerebellar granule cells, and macroglia. J Biol Chem 2002;277:7598-609. [Crossref] [PubMed]
- Sutherland AP, Zhang H, Zhang Y, et al. Zinc finger protein Zbtb20 is essential for postnatal survival and glucose homeostasis. Mol Cell Biol 2009;29:2804-15. [Crossref] [PubMed]
- Xu X, Xie Q, Xie M, et al. LncRNA SNHG8 Serves as an Oncogene in Breast Cancer Through miR-634/ZBTB20 Axis. Cancer Manag Res 2021;13:3017-28. [Crossref] [PubMed]
- Zhang Y, Zhou X, Zhang M, et al. ZBTB20 promotes cell migration and invasion of gastric cancer by inhibiting IκBα to induce NF-κB activation. Artif Cells Nanomed Biotechnol 2019;47:3862-72. [Crossref] [PubMed]
- Terekhova M, Swain A, Bohacova P, et al. Single-cell atlas of healthy human blood unveils age-related loss of NKG2C(+)GZMB(-)CD8(+) memory T cells and accumulation of type 2 memory T cells. Immunity 2023;56:2836-2854.e9. [Crossref] [PubMed]
- Zhao X, Hu W, Park SR, et al. The transcriptional cofactor Tle3 reciprocally controls effector and central memory CD8(+) T cell fates. Nat Immunol 2024;25:294-306. [Crossref] [PubMed]
- Liu S, Wang P, Wang P, et al. Tissue-resident memory CD103+CD8+ T cells in colorectal cancer: its implication as a prognostic and predictive liver metastasis biomarker. Cancer Immunol Immunother 2024;73:176. [Crossref] [PubMed]
- Boonpiyathad T, Sözener ZC, Satitsuksanoa P, et al. Immunologic mechanisms in asthma. Semin Immunol 2019;46:101333. [Crossref] [PubMed]
- Yu Y. Multi-target combinatory strategy to overcome tumor immune escape. Front Med 2022;16:208-15. [Crossref] [PubMed]
- Zhang L, Yu X, Zheng L, et al. Lineage tracking reveals dynamic relationships of T cells in colorectal cancer. Nature 2018;564:268-72. [Crossref] [PubMed]
- León B, Ballesteros-Tato A. Modulating Th2 Cell Immunity for the Treatment of Asthma. Front Immunol 2021;12:637948. [Crossref] [PubMed]
- Schreiber S, Hammers CM, Kaasch AJ, et al. Metabolic Interdependency of Th2 Cell-Mediated Type 2 Immunity and the Tumor Microenvironment. Front Immunol 2021;12:632581. [Crossref] [PubMed]
- Liu X, Wang X, Yang Q, et al. Th17 Cells Secrete TWEAK to Trigger Epithelial-Mesenchymal Transition and Promote Colorectal Cancer Liver Metastasis. Cancer Res 2024;84:1352-71. [Crossref] [PubMed]
- Qin D, Zhang Y, Shu P, et al. Targeting tumor-infiltrating tregs for improved antitumor responses. Front Immunol 2024;15:1325946. [Crossref] [PubMed]
- Ohue Y, Nishikawa H, Regulatory T. Treg) cells in cancer: Can Treg cells be a new therapeutic target? Cancer Sci 2019;110:2080-9. [Crossref] [PubMed]
- Fischer K, Hoffmann P, Voelkl S, et al. Inhibitory effect of tumor cell-derived lactic acid on human T cells. Blood 2007;109:3812-9. [Crossref] [PubMed]
- Husain Z, Huang Y, Seth P, et al. Tumor-derived lactate modifies antitumor immune response: effect on myeloid-derived suppressor cells and NK cells. J Immunol 2013;191:1486-95. [Crossref] [PubMed]
- Su J, Morgani SM, David CJ, et al. TGF-β orchestrates fibrogenic and developmental EMTs via the RAS effector RREB1. Nature 2020;577:566-71. [Crossref] [PubMed]
- Chakravarthy A, Khan L, Bensler NP, et al. TGF-β-associated extracellular matrix genes link cancer-associated fibroblasts to immune evasion and immunotherapy failure. Nat Commun 2018;9:4692. [Crossref] [PubMed]
- Yang L, Moses HL. Transforming growth factor beta: tumor suppressor or promoter? Are host immune cells the answer? Cancer Res 2008;68:9107-11. [Crossref] [PubMed]
- Wang J, Zhao X, Wan YY. Intricacies of TGF-β signaling in Treg and Th17 cell biology. Cell Mol Immunol 2023;20:1002-22. [Crossref] [PubMed]
- Derynck R, Zhang YE. Smad-dependent and Smad-independent pathways in TGF-beta family signalling. Nature 2003;425:577-84. [Crossref] [PubMed]
- Fabregat I, Malfettone A, Soukupova J. New Insights into the Crossroads between EMT and Stemness in the Context of Cancer. J Clin Med 2016;5:37. [Crossref] [PubMed]
- Xue W, Yang L, Chen C, et al. Wnt/β-catenin-driven EMT regulation in human cancers. Cell Mol Life Sci 2024;81:79. [Crossref] [PubMed]

