KMT2C is a potential biomarker of prognosis and chemotherapy sensitivity in breast cancer

Pan-cancer analysis and experimental verification of cytochrome B561 as a prognostic and therapeutic biomarker in breast cancer

OBJECTIVE

This study investigates Cytochrome B561 (CYB561) expression in Pan-Cancer, its relationship with immune invasion, and its prognostic value in Breast Cancer (BRCA) patients.

METHODS

Data from The Cancer Genome Atlas (TCGA) were analyzed. CYB561 expression in normal and tumor tissues was examined, with correlations to immune invasion, mutation, and immune checkpoints. Wilcoxon rank-sum test assessed expression differences. Gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses were conducted. Logistic regression, Kaplan-Meier, and Cox regression analyses evaluated clinicopathological features and survival outcomes. A Cox multivariate analysis-based Nomogram predicted CYB561's prognostic impact. CYB561 knockout in breast cancer cells assessed functional effects. Single-cell RNA sequencing identified prognostic biomarkers.

RESULTS

CYB561 was highly expressed in most tumors. BRCA showed the highest correlation with ESTIMATE scores and significant negative correlation with immune checkpoints. High CYB561 expression correlated with specific clinicopathological features and survival outcomes. The nomogram predicted BRCA prognosis. CYB561 knockout inhibited breast cancer cell proliferation. Seven predictive agents for CYB561 inhibition were identified.

CONCLUSIONS

CYB561 exhibits aberrant expression in tumors, particularly in BRCA, and serves as a predictive marker for immune-related therapies and a prognostic indicator in BRCA.

m6A related metabolic genes in breast cancer and their relationship with prognosis

Breast cancer (BC) is the most prevalent malignancy among women, with incidence rates rising annually. N6-methyladenosine (m6A) modification has been recognized as a key regulator in the onset and progression of BC. Nevertheless, the role of m6A-associated metabolic genes (mMGs) in BC regulation remains insufficiently understood. In this study, we first analyzed and clustered single-cell transcriptomic (scRNA-seq) data from the peripheral blood of BC patients. Differentially expressed genes (DEGs) across various cell populations were intersected with mMGs to identify differentially expressed mMGs (DEmMGs). The AUCell algorithm was employed to score DEmMGs across cell populations, followed by subgroup clustering of high-scoring cell types. Additionally, DEGs from BC and control transcriptomic (RNA-seq) data in The Cancer Genome Atlas (TCGA) were intersected with DEmMGs. BC subtypes were identified based on the expression levels of overlapping genes, and differences in survival rates and immune microenvironment characteristics were examined across subtypes. A BC risk model was constructed using Lasso-Cox regression, and variations in prognosis, tumor mutational burden, immune cell infiltration, and drug sensitivity were explored. Finally, real-time quantitative PCR (qRT-PCR) and Western blot (WB) were used to validate the identified prognostic genes. NDUFAB1, VDAC1, TYMP, UGDH, ATP6AP1, and ALDH2 showed consistent and significant expression differences between the BC and control groups. This study's identification of key prognostic genes and the development of a risk model offer potential new targets for therapeutic intervention and clinical management of BC.

Clinical, pathological and gene expression profiling of estrogen receptor discordance in breast cancer

BACKGROUND

Breast cancer (BC) is the world's largest tumor species in which hormone receptor-positive patients have relatively good prognosis. However, majority of patients will develop late resistance, one of the important factors is due to the loss of the original estrogen receptor (ER) expression.

METHODS

We conducted this study in 115 patients with BC who experienced second biopsy at Jiangsu Province Hospital (JSPH) and divided patients into two subgroups ER + to - and ER + to + . First, clinicopathological characteristics between two groups were evaluated. Second, we explored candidate genes related to BC ER intratumor heterogeneity by applying next-generation sequencing (NGS) in 42 patients. Multi-omics integrative analysis of tumor transcriptomic, cancer-related pathway, diagnostic and prognostic value and immune profile were conducted. Besides, preliminary assay were also used to evaluate the correlation between KMT2C and ERα (ESR1) expression. The CCK-8, 5-Ethynyl-2'-deoxyuridine (EdU) assays, Transwell assays and the wound scratch tests were applied to explore the cellular interactions between KMT2C and BC.

RESULTS

We find the histological type (p = 0.008) and disease-free survival (DFS) (p = 0.004) were significantly different in two subgroups. In Cox survival analysis, metastasis (Hazard ratio (HR) > 1, p = 0.007) and neo-adjuvant (HR < 1, p < 0.001) are independent prognostic factors of DFS. Besides, by analyzing NGS results, we found four genes KMT2C, FGFR19, FGF1 and FGF4 were highly mutated genes in ER + to - subgroup. Furthermore, the gene KMT2C displayed significant diagnostic value and prognostic value in BC and pan-cancer. In addition, a positive correlation between KMT2C expression and immune infiltrating levels of T cell CD4 + , macrophage and neutrophil was found. In the end, Western blot and RT-qPCR assay were used and found KMT2C and ERα (ESR1) expressions are strongly positive correlated in mRNA and protein level. Inhibition of KMT2C significantly reduced proliferation, invasion, and migration of MCF7 cells.

CONCLUSION

People in two cohorts from JSPH presented different clinical characteristics and prognosis. The gene KMT2C may affect the progression of BC by regulating the molecular, epigenetic activity and immune infiltration. It may also serve as a novel prognostic biomarker for BC patients who underwent ER status converted from positive to negative.

A long-term prognosis study of human USP8-mutated ACTH-secreting pituitary neuroendocrine tumours

OBJECTIVE

Somatic variants in the ubiquitin-specific protease 8 (USP8) gene are the most common genetic cause of Cushing disease. We aimed to explore the relationship between clinical outcomes and USP8 status in a single centre.

DESIGN, PATIENTS AND MEASUREMENTS

We investigated the USP8 status in 48 patients with pituitary corticotroph tumours. A median of 62 months of follow-up was conducted after surgery from November 2013 to January 2015. The clinical, biochemical and imaging features were collected and analysed.

RESULTS

Seven USP8 variants (p.Ser718Pro, p.Ser719del, p.Pro720Arg, p.Pro720Gln, p.Ser718del, p.Ser718Phe, p.Lys713Arg) were identified in 24 patients (50%). USP8 variants showed a female predominance (100% vs. 75% in wild type [WT], p = .022). Patients with p.Ser719del showed an older age at surgery compared to patients with the p.Pro720Arg variant (47- vs. 24-year-olds, p = .033). Patients with p.Pro720Arg showed a higher rate of macroadenoma compared to patients harbouring the p.Ser718Pro variant (60% vs. 0%, p = .037). No significant differences were observed in serum and urinary cortisol and adrenocorticotropin hormone (ACTH) levels. Immediate surgical remission (79% vs. 75%) and long-term hormone remission (79% vs. 67%) were not significantly different between the two groups. The recurrence rate was 21% (4/19) in patients harbouring USP8 variants and 13% (2/16) in WT patients. Recurrence-free survival presented a tendency to be shorter in USP8-mutated individuals (76.7 vs. 109.2 months, p = .068).

CONCLUSIONS

Somatic USP8 variants accounted for 50% of the genetic causes in this cohort with a significant female frequency. A long-term follow-up revealed a tendency toward shorter recurrence-free survival in USP8-mutant patients.

KMT2C and KMT2D aberrations in breast cancer

KMT2C and KMT2D are histone lysine methyltransferases responsible for the monomethylation of histone 3 lysine 4 (H3K4) residues at gene enhancer sites. KMT2C/D are the most frequently mutated histone methyltransferases (HMTs) in breast cancer, occurring at frequencies of 10-20% collectively. Frequent damaging and truncating somatic mutations indicate a tumour-suppressive role of KMT2C/D in breast oncogenesis. Recent studies using cell lines and mouse models to replicate KMT2C/D loss show that these genes contribute to oestrogen receptor (ER)-driven transcription in ER+ breast cancers through the priming of gene enhancer regions. This review provides an overview of the functions of KMT2C/D and outlines the recent clinical and experimental evidence of the roles of KMT2C and KMT2D in breast cancer development.

Germline variants in early and late-onset Brazilian prostate cancer patients

BACKGROUND

The median age for Prostate Cancer (PCa) diagnosis is 66 years, but 10% are diagnosed before 55 years. Studies on early-onset PCa remain both limited and controversial. This investigation sought to identify and characterize germline variants within Brazilian PCa patients classified as either early or later onset disease.

METHODS

Peripheral blood DNA from 71 PCa patients: 18 younger (≤ 55 years) and 53 older (≥ 60 years) was used for Targeted DNA sequencing of 20 genes linked to DNA damage response, transcriptional regulation, cell cycle, and epigenetic control. Subsequent genetic variant identification was performed and variant functional impacts were analyzed with in silico prediction.

RESULTS

A higher frequency of variants in the BRCA2 and KMT2C genes across both age groups. KMT2C has been linked to the epigenetic dysregulation observed during disease progression in PCa. We present the first instance of KMT2C mutation within the blood of Brazilian PCa patients. Furthermore, out of the recognized variants within the KMT2C gene, 7 were designated as deleterious. Thirteen deleterious variants were exclusively detected in the younger group, while the older group exhibited 37 variants. Within these findings, 4 novel variants emerged, including 1 designated as pathogenic.

CONCLUSIONS

Our findings contribute to a deeper understanding of the genetic factors associated with PCa susceptibility in different age groups, especially among the Brazilian population. This is the first investigation to explore germline variants specifically in younger Brazilian PCa patients, with high relevance given the genetic diversity of the population in Brazil. Additionally, our work presents evidence of functionally deleterious germline variants within the KMT2C gene among Brazilian PCa patients. The identification of novel and functionally significant variants in the KMT2C gene emphasizes its potential role in PCa development and warrants further investigation.

Pathogenic alterations in PIK3CA and KMT2C are frequent and independent prognostic factors in anal squamous cell carcinoma treated with salvage abdominoperineal resection

The management of anal squamous cell carcinoma (ASCC) has yet to experience the transformative impact of precision medicine. Conducting genomic analyses may uncover novel prognostic biomarkers and offer potential directions for the development of targeted therapies. To that end, we assessed the prognostic and theragnostic implications of pathogenic variants identified in 571 cancer-related genes from surgical samples collected from a homogeneous, multicentric French cohort of 158 ASCC patients who underwent abdominoperineal resection treatment. Alterations in PI3K/AKT/mTOR, chromatin remodeling, and Notch pathways were frequent in HPV-positive tumors, while HPV-negative tumors often harbored variants in cell cycle regulation and genome integrity maintenance genes (e.g., frequent TP53 and TERT promoter mutations). In patients with HPV-positive tumors, KMT2C and PIK3CA exon 9/20 pathogenic variants were associated with worse overall survival in multivariate analysis (Hazard ratio (HR)KMT2C  = 2.54, 95%CI = [1.25,5.17], P value = .010; HRPIK3CA  = 2.43, 95%CI = [1.3,4.56], P value = .006). Alterations with theragnostic value in another cancer type was detected in 43% of patients. These results suggest that PIK3CA and KMT2C pathogenic variants are independent prognostic factors in patients with ASCC with HPV-positive tumors treated by abdominoperineal resection. And, importantly, the high prevalence of alterations bearing potential theragnostic value strongly supports the use of genomic profiling to allow patient enrollment in precision medicine clinical trials.

Somatic mutations in a multigene panel and impact on prognosis based on TP53 status in Chinese HER2-positive patients undergoing neoadjuvant therapy: A single-institution retrospective cohort

BACKGROUND

Gene mutations play a crucial role in the occurrence and development of tumors, particularly in breast cancer (BC). Neoadjuvant therapy (NAT) has shown greater clinical benefit in HER2-positive breast cancer. However, further clinical investigation is needed to fully understand the correlation between genetic mutations and NAT efficacy and the long-term prognosis in HER2-positive BC.

METHODS

This was a retrospective cohort study of 222 patients receiving NAT between 2017 and 2021 in the Department of Breast Surgery of Fudan University Shanghai Cancer Center. Tumor samples from these patients were subjected to Next Generation Sequencing (NGS) to analyze mutations in 513 cancer-related genes. This study aimed to investigate the association between these genetic mutations and postoperative pathological complete response (pCR), as well as their impact on disease-free survival (DFS).

RESULTS

In total, 48.65% patients reached pCR, ER-negative status (p < 0.001), PR-negative status (p < 0.001), Ki67 ≥ 20 (p = 0.011), and dual-targeted therapy (p < 0.001) were all associated with enhanced pCR rates. The frequency of somatic alterations in TP53 (60%), PIK3CA (15%), and ERBB2 (11%) was highest. In the HER2+/HR- cohort, patients who achieved pCR had a significant benefit in prognosis (HR = 3.049, p = 0.0498). KMT2C (p = 0.036) and TP53 (p = 0.037) mutations were significantly increased in patients with DFS events. Moreover, TP53 mutations had prognostic significance in HER2-positive BC patients with HR-negative (HR = 3.712, p = 0.027) and pCR (HR = 6.253, p = 0.027) status and who received herceptin-only targeted therapy (HR = 4.145, p = 0.011).

CONCLUSIONS

The genetic mutation profiles of Chinese HER2+ patients who received NAT were discrepant with respect to HR status or DFS events. TP53 mutations have significant prognostic value in patients with NAT for HER2-positive BC and patients benefit differently depending on HR status, the neoadjuvant regimen and response, which highlights the significance of genetic factors in treatment customization based on individual genetic and clinical characteristics.

The Expression of ZNF268 and Its Role in The Cisplatin-based Chemoresistance of Breast Cancer

Objective

Breast cancer is one of the most prevalent cancers in females worldwide and is one of the leading causes of cancer death and disability in women. Multiple therapies have been applied to breast cancer treatment; however, the long-term survival rate remains low. Although cisplatin has been widely utilized for cancer therapy, chemoresistance still influences the outcome.

Methods

After collecting the breast cancer cell line MDA-MB-231 treated with or without cisplatin and sample information from The Cancer Genome Atlas Program (TCGA), we screened out their common parameters and influences on the prognoses of patients' potential targets. Surgical excisional tissue sections of patients with breast cancer who were admitted and treated in the Department of Breast and Thyroid Surgery, Liuzhou People's Hospital from 2017 to 2020 was collected and follow up. After a series of assays combined with clinical information, we tested the reliability of the target.

Results

We found that a high expression level of ZNF268 in breast cancer cell lines significantly enhances the sensitivity to cisplatin, contrary to the effects of low expression. Furthermore, a significantly worse prognosis was observed in patients with a high expression of ZNF268 after cisplatin chemotherapy.

Conclusion

The expression level of ZNF268 in breast cancer patients after cisplatin chemotherapy may become a potential target to predict the chemoresistance of patients to cisplatin. This study provides a novel idea for improving breast cancer treatment and survival rates.

Deciphering Brain Metastasis Stem Cell Properties From Colorectal Cancer Highlights Specific Stemness Signature and Shared Molecular Features

BACKGROUND & AIMS

Brain metastases (BMs) from colorectal cancer (CRC) are associated with significant morbidity and mortality, with chemoresistance and short overall survival. Migrating cancer stem cells with the ability to initiate BM have been described in breast and lung cancers. In this study, we describe the identification and characterization of cancer stem cells in BM from CRC.

METHODS

Four brain metastasis stem cell lines from patients with colorectal cancer (BM-SC-CRC1 to BM-SC-CRC4) were obtained by mechanical dissociation of patient's tumors and selection of cancer stem cells by appropriate culture conditions. BM-SC-CRCs were characterized in vitro by clonogenic and limiting-dilution assays, as well as immunofluorescence and Western blot analyses. In ovo, a chicken chorioallantoic membrane (CAM) model and in vivo, xenograft experiments using BALB/c-nude mice were realized. Finally, a whole exome and RNA sequencing analyses were performed.

RESULTS

BM-SC-CRC formed metaspheres and contained tumor-initiating cells with self-renewal properties. They expressed stem cell surface markers (CD44v6, CD44, and EpCAM) in serum-free medium and CRC markers (CK19, CK20 and CDX-2) in fetal bovine serum-enriched medium. The CAM model demonstrated their invasive and migratory capabilities. Moreover, mice intracranial xenotransplantation of BM-SC-CRCs adequately recapitulated the original patient BM phenotype. Finally, transcriptomic and genomic approaches showed a significant enrichment of invasiveness and specific stemness signatures and highlighted KMT2C as a potential candidate gene to potentially identify high-risk CRC patients.

CONCLUSIONS

This original study represents the first step in CRC BM initiation and progression comprehension, and further investigation could open the way to new therapeutics avenues to improve patient prognosis.

Molecular profiling of a real-world breast cancer cohort with genetically inferred ancestries reveals actionable tumor biology differences between European ancestry and African ancestry patient populations

BACKGROUND

Endocrine-resistant HR+/HER2- breast cancer (BC) and triple-negative BC (TNBC) are of interest for molecularly informed treatment due to their aggressive natures and limited treatment profiles. Patients of African Ancestry (AA) experience higher rates of TNBC and mortality than European Ancestry (EA) patients, despite lower overall BC incidence. Here, we compare the molecular landscapes of AA and EA patients with HR+/HER2- BC and TNBC in a real-world cohort to promote equity in precision oncology by illuminating the heterogeneity of potentially druggable genomic and transcriptomic pathways.

METHODS

De-identified records from patients with TNBC or HR+/HER2- BC in the Tempus Database were randomly selected (N = 5000), with most having stage IV disease. Mutations, gene expression, and transcriptional signatures were evaluated from next-generation sequencing data. Genetic ancestry was estimated from DNA-seq. Differences in mutational prevalence, gene expression, and transcriptional signatures between AA and EA were compared. EA patients were used as the reference population for log fold-changes (logFC) in expression.

RESULTS

After applying inclusion criteria, 3433 samples were evaluated (n = 623 AA and n = 2810 EA). Observed patterns of dysregulated pathways demonstrated significant heterogeneity among the two groups. Notably, PIK3CA mutations were significantly lower in AA HR+/HER2- tumors (AA = 34% vs. EA = 42%, P < 0.05) and the overall cohort (AA = 28% vs. EA = 37%, P = 2.08e-05). Conversely, KMT2C mutation was significantly more frequent in AA than EA TNBC (23% vs. 12%, P < 0.05) and HR+/HER2- (24% vs. 15%, P = 3e-03) tumors. Across all subtypes and stages, over 8000 genes were differentially expressed between the two ancestral groups including RPL10 (logFC = 2.26, P = 1.70e-162), HSPA1A (logFC = - 2.73, P = 2.43e-49), ATRX (logFC = - 1.93, P = 5.89e-83), and NUTM2F (logFC = 2.28, P = 3.22e-196). Ten differentially expressed gene sets were identified among stage IV HR+/HER2- tumors, of which four were considered relevant to BC treatment and were significantly enriched in EA: ERBB2_UP.V1_UP (P = 3.95e-06), LTE2_UP.V1_UP (P = 2.90e-05), HALLMARK_FATTY_ACID_METABOLISM (P = 0.0073), and HALLMARK_ANDROGEN_RESPONSE (P = 0.0074).

CONCLUSIONS

We observed significant differences in mutational spectra, gene expression, and relevant transcriptional signatures between patients with genetically determined African and European ancestries, particularly within the HR+/HER2- BC and TNBC subtypes. These findings could guide future development of treatment strategies by providing opportunities for biomarker-informed research and, ultimately, clinical decisions for precision oncology care in diverse populations.

Activated Mast Cells Combined with NRF2 Predict Prognosis for Esophageal Cancer

Background

Esophageal cancer (EC) had the sixth-highest mortality rate of all cancers due to its poor prognosis. Immune cells and mutation genes influenced the prognosis of EC, but their combined effect on predicting EC prognosis was unknown. In this study, we comprehensively analyzed the immune cell infiltration (ICI) and mutation genes and their combined effects for predicting prognosis in EC.

Methods

The CIBERSORT and ESTIMATE algorithms were used to analyse the ICI scape based on the TCGA and GEO databases. EC tissues and pathologic sections from Huai'an, China, were used to verify the key immune cells and mutation genes and their interactions.

Results

Stromal/immune score patterns and ICI/gene had no statistical significance in overall survival (OS) (p > 0.05). The combination of ICI and tumor mutation burden (TMB) showed that the high TMB and high ICI score group had the shortest OS (p = 0.004). We recognized that the key mutation gene NRF2 was significantly different in the high/low ICI score subgroups (p = 0.002) and positivity with mast cells (MCs) (p < 0.05). Through experimental validation, we found that the MCs and activated mast cells (AC-MCs) were more infiltration in stage II/III (p = 0.032; p = 0.013) of EC patients and that NRF2 expression was upregulated in EC (p = 0.045). AC-MCs combined with NRF2 had a poor prognosis, according to survival analysis (p = 0.056) and interactive analysis (p = 0.032).

Conclusions

We presume that NRF2 combined with AC-MCs could be a marker to predict prognosis and could influence immunotherapy through regulating PD-L1 in the EC.

Analysis of matched primary and recurrent BRCA1/2 mutation-associated tumors identifies recurrence-specific drivers

Recurrence is a major cause of death among BRCA1/2 mutation carriers with breast (BrCa) and ovarian cancers (OvCa). Herein we perform multi-omic sequencing on 67 paired primary and recurrent BrCa and OvCa from 27 BRCA1/2 mutation carriers to identify potential recurrence-specific drivers. PARP1 amplifications are identified in recurrences (False Discovery Rate q = 0.05), and PARP1 is significantly overexpressed across primary BrCa and recurrent BrCa and OvCa, independent of amplification status. RNA sequencing analysis finds two BRCA2 isoforms, BRCA2-201/Long and BRCA2-001/Short, respectively predicted to be sensitive and insensitive to nonsense-mediated decay. BRCA2-001/Short is expressed more frequently in recurrences and associated with reduced overall survival in breast cancer (87 vs. 121 months; Hazard Ratio = 2.5 [1.18-5.5]). Loss of heterozygosity (LOH) status is discordant in 25% of patient's primary and recurrent tumors, with switching between both LOH and lack of LOH found. Our study reveals multiple potential drivers of recurrent disease in BRCA1/2 mutation-associated cancer, improving our understanding of tumor evolution and suggesting potential biomarkers.

KMT2C mutation in a Chinese man with primary multidrug-resistant metastatic adenocarcinoma of rete testis: a case report

Background

Adenocarcinoma of the rete testis (AORT) is an extremely rare malignant tumor with poor prognosis and limited responsiveness to traditional chemotherapy. Few previous studies have focused on the molecular mechanisms underlying therapy resistance in AORT and further scrutiny is required to enable searches for targeted drugs to guide treatment selection.

Case presentation

The current case concerns a 55-year-old man with AORT who presented with isolated bone metastasis at initial diagnosis and experienced rapid disease progression after multi-line platinum-based combination chemotherapy. Next-generation sequencing revealed a novel somatic lysine methyltransferase 2C (KMT2C) c.5605 T > C mutation in exon 36 with an abundance of 49.27%. The patient received antiangiogenic drug treatment for 2 months but this was discontinued due to unacceptable anorexia and nausea. He survived for 12 months after diagnosis.

Conclusion

A potential correlation between AORT primary multi-drug resistance and KMT2C mutations is implied. Further studies are needed to determine the efficacy of PARP1/2 inhibitors for tumors with KMT2C mutations.

Integration of metabolites from meta-analysis with transcriptome reveals enhanced SPHK1 in PDAC with a background of pancreatitis

BACKGROUND

Pathophysiology of transformation of inflammatory lesions in chronic pancreatitis (CP) to pancreatic ductal adenocarcinoma (PDAC) is not clear.

METHODS

We conducted a systematic review, meta-analysis of circulating metabolites, integrated this data with transcriptome analysis of human pancreatic tissues and validated using immunohistochemistry. Our aim was to establish biomarker signatures for early malignant transformation in patients with underlying CP and identify therapeutic targets.

RESULTS

Analysis of 19 studies revealed AUC of 0.86 (95% CI 0.81-0.91, P < 0.0001) for all the altered metabolites (n = 88). Among them, lipids showed higher differentiating efficacy between PDAC and CP; P-value (< 0.0001). Pathway enrichment analysis identified sphingomyelin metabolism (impact value-0.29, FDR of 0.45) and TCA cycle (impact value-0.18, FDR of 0.06) to be prominent pathways in differentiating PDAC from CP. Mapping circulating metabolites to corresponding genes revealed 517 altered genes. Integration of these genes with transcriptome data of CP and PDAC with a background of CP (PDAC-CP) identified three upregulated genes; PIGC, PPIB, PKM and three downregulated genes; AZGP1, EGLN1, GNMT. Comparison of CP to PDAC-CP and PDAC-CP to PDAC identified upregulation of SPHK1, a known oncogene.

CONCLUSIONS

Our analysis suggests plausible role for SPHK1 in development of pancreatic adenocarcinoma in long standing CP patients. SPHK1 could be further explored as diagnostic and potential therapeutic target.