Functional Roles of SPINK1 in Cancers

High expression of serine protease inhibitor kazal type 1 predicts poor prognosis and promotes the progression and invasion of oral tongue squamous cell carcinoma

OBJECTIVE

This study aimed to investigate the expression of serine protease inhibitor kazal type 1 (SPINK1) and its carcinogenic effect in oral tongue squamous cell carcinoma (OTSCC).

DESIGN

Initially, bioinformatics analysis was conducted using data from The Cancer Genome Atlas and Gene Expression Omnibus to compare SPINK1 mRNA expression between malignant and adjacent tissues. Subsequently, the impact of differential expression on survival and other clinical variables was examined. Additionally, histology microarray analysis was performed to assess SPINK1 protein expression in 35 cases of malignant and adjacent tissues. Finally, alterations in SPINK1 expression were evaluated to determine its biological phenotypes in OTSCC, including proliferation, apoptosis, invasion, and metastasis.

RESULTS

OTSCC tissues exhibit higher levels of SPINK1 compared to surrounding cancerous tissues. Notably, increased SPINK1 expression correlates with the pathological N stage and independently predicts overall survival among patients with OTSCC.

CONCLUSION

Suppression of SPINK1 inhibited OTSCC cell proliferation, invasion, and motility while promoting apoptosis. These findings suggest that SPINK1 may serve as a prognostic biomarker as well as a potential therapeutic target for managing OTSCC.

Effect of cellular senescence on the response of human peritoneal mesothelial cells to TGF-β

Transforming growth factor β (TGF-β) is implicated in both mesothelial-to-mesenchymal transition (MMT) and cellular senescence of human peritoneal mesothelial cells (HPMCs). We previously showed that senescent HPMCs could spontaneously acquire some phenotypic features of MMT, which in young HPMCs were induced by TGF-β. Here, we used electron microscopy, as well as global gene and protein profiling to assess in detail how exposure to TGF-β impacts on young and senescent HPMCs in vitro. We found that TGF-β induced structural changes consistent with MMT in young, but not in senescent HPMCs. Of all genes and proteins identified reliably in HPMCs across all treatments and states, 4,656 targets represented overlapping genes and proteins. Following exposure to TGF-β, 137 proteins and 46 transcripts were significantly changed in young cells, compared to 225 proteins and only 2 transcripts in senescent cells. Identified differences between young and senescent HPMCs were related predominantly to wound healing, integrin-mediated signalling, production of proteases and extracellular matrix components, and cytoskeleton structure. Thus, the response of senescent HPMCs to TGF-β differs or is less pronounced compared to young cells. As a result, the character and magnitude of the postulated contribution of HPMCs to TGF-β-induced peritoneal remodelling may change with cell senescence.

Prediction of response to promising first-line chemotherapy in ovarian cancer patients with residual peritoneal tumors: practical biomarkers and robust multiplex models

BACKGROUND

Platinum/taxane (TC) chemotherapy with debulking surgery stays the mainstay of the treatment in ovarian cancer patients with peritoneal metastasis, and recently its novel modality, intraperitoneal carboplatin with dose-dense paclitaxel (ddTCip), was shown to have greater therapeutic impact. Nevertheless, the response varies among patients and consequent recurrence, or relapse often occurs. Discovery of therapeutic response predictor to ddTCip and/or TC therapy is eagerly awaited to improve the treatment outcome.

METHODS

Using datasets in 76 participants in our ddTCip study and published databases on patients received TC therapy, we first validated a total of 75 previously suggested markers, sought out more active biomarkers through the association analyses of genome-wide transcriptome and genotyping data with progression-free survival (PFS) and adverse events, and then developed multiplex statistical prediction models for PFS and toxicity by mainly using multiple regression analysis and the classification and regression tree (CART) algorithm.

RESULTS

The association analyses revealed that SPINK1 could be a possible biomarker of ddTCip efficacy, while ABCB1 rs1045642 and ERCC1 rs11615 would be a predictor of hematologic toxicity and peripheral neuropathy, respectively. Multiple regression analyses and CART algorithm finally provided a potent efficacy prediction model using 5 gene expression data and robust multiplex toxicity prediction models-CART models using a total of 4 genotype combinations and multiple regression models using 15 polymorphisms on 12 genes.

CONCLUSION

Biomarkers and multiplex models composed here could work well in the response prediction of ddTCip and/or TC therapy, which might contribute to realize optimal selection of the key therapy.

Early life exposure to vitamin D deficiency impairs molecular mechanisms that regulate liver cholesterol biosynthesis, energy metabolism, inflammation, and detoxification

Introduction

Emerging data suggests liver disease may be initiated during development when there is high genome plasticity and the molecular pathways supporting liver function are being developed.

Methods

Here, we leveraged our Collaborative Cross mouse model of developmental vitamin D deficiency (DVD) to investigate the role of DVD in dysregulating the molecular mechanisms underlying liver disease. We defined the effects on the adult liver transcriptome and metabolome and examined the role of epigenetic dysregulation. Given that the parental origin of the genome (POG) influences response to DVD, we used our established POG model [POG1-(CC011xCC001)F1 and POG2-(CC001xCC011)F1] to identify interindividual differences.

Results

We found that DVD altered the adult liver transcriptome, primarily downregulating genes controlling liver development, response to injury/infection (detoxification & inflammation), cholesterol biosynthesis, and energy production. In concordance with these transcriptional changes, we found that DVD decreased liver cell membrane-associated lipids (including cholesterol) and pentose phosphate pathway metabolites. Each POG also exhibited distinct responses. POG1 exhibited almost 2X more differentially expressed genes (DEGs) with effects indicative of increased energy utilization. This included upregulation of lipid and amino acid metabolism genes and increased intermediate lipid and amino acid metabolites, increased energy cofactors, and decreased energy substrates. POG2 exhibited broader downregulation of cholesterol biosynthesis genes with a metabolomics profile indicative of decreased energy utilization. Although DVD primarily caused loss of liver DNA methylation for both POGs, only one epimutation was shared, and POG2 had 6.5X more differentially methylated genes. Differential methylation was detected at DEGs regulating developmental processes such as amino acid transport (POG1) and cell growth & differentiation (e.g., Wnt & cadherin signaling, POG2).

Conclusions

These findings implicate a novel role for maternal vitamin D in programming essential offspring liver functions that are dysregulated in liver disease. Importantly, impairment of these processes was not rescued by vitamin D treatment at weaning, suggesting these effects require preventative measures. Substantial differences in POG response to DVD demonstrate that the parental genomic context of exposure determines offspring susceptibility.

LARP4B promotes hepatocellular carcinoma progression and impairs sorafenib efficacy by activating SPINK1-mediated EGFR pathway

La-related proteins (LARPs) regulate gene expression by binding to RNAs and exhibit critical effects on disease progression, including tumors. However, the role of LARP4B and its underlying mechanisms in the progression of hepatocellular carcinoma (HCC) remain largely unclear. In this study, we found that LARP4B expression is upregulated and correlates with poor prognosis in patients with HCC. Gain- and loss-of-function assays showed that LARP4B promotes stemness, proliferation, metastasis, and angiogenesis in vitro and in vivo. Furthermore, LARP4B inhibition enhances the antitumor effects of sorafenib and blocks the metastasis-enhancing effects of low sorafenib concentrations in HCC. Mechanistically, LARP4B expression is upregulated by METTL3-mediated N6-methyladenosine (m6A)-IGF2BP3-dependent modification in HCC. RNA- and RNA immunoprecipitation (RIP)- sequencing uncovered that LARP4B upregulates SPINK1 by binding to SPINK1 mRNA via the La motif and maintaining mRNA stability. LARP4B activates the SPINK1-mediated EGFR signaling pathway, which supports stemness, progression and sorafenib resistance in HCC. Additionally, a positive feedback loop with the LARP4B/SPINK1/p-AKT/C/EBP-β axis is responsible for the sorafenib-therapeutic benefit of LARP4B depletion. Overall, this study demonstrated that LARP4B facilitates HCC progression, and LARP4B inhibition provides benefits to sorafenib treatment in HCC, suggesting that LARP4B might be a potential therapeutic target for HCC.

Association of clinicopathologic and molecular factors with the occurrence of positive margins in breast cancer

PURPOSE

To explore the association of clinicopathologic and molecular factors with the occurrence of positive margins after first surgery in breast cancer.

METHODS

The clinical and RNA-Seq data for 951 (75 positive and 876 negative margins) primary breast cancer patients from The Cancer Genome Atlas (TCGA) were used. The role of each clinicopathologic factor for margin prediction and also their impact on survival were evaluated using logistic regression, Fisher's exact test, and Cox proportional hazards regression models. In addition, differential expression analysis on a matched dataset (71 positive and 71 negative margins) was performed using Deseq2 and LASSO regression.

RESULTS

Association studies showed that higher stage, larger tumor size (T), positive lymph nodes (N), and presence of distant metastasis (M) significantly contributed (p ≤ 0.05) to positive surgical margins. In case of surgery, lumpectomy was significantly associated with positive margin compared to mastectomy. Moreover, PAM50 Luminal A subtype had higher chance of positive margin resection compared to Basal-like subtype. Survival models demonstrated that positive margin status along with higher stage, higher TNM, and negative hormone receptor status was significant for disease progression. We also found that margin status might be a surrogate of tumor stage. In addition, 29 genes that could be potential positive margin predictors and 8 pathways were identified from molecular data analysis.

CONCLUSION

The occurrence of positive margins after surgery was associated with various clinical factors, similar to the findings reported in earlier studies. In addition, we found that the PAM50 intrinsic subtype Luminal A has more chance of obtaining positive margins compared to Basal type. As the first effort to pursue molecular understanding of the margin status, a gene panel of 29 genes including 17 protein-coding genes was also identified for potential prediction of the margin status which needs to be validated using a larger sample set.

Leveraging a KRAS-based signature to predict the prognosis and drug sensitivity of colon cancer and identifying SPINK4 as a new biomarker

KRAS is one of the leading mutations reported in colon cancer. However, there are few studies on the application of KRAS related signature in predicting prognosis and drug sensitivity of colon cancer patient. We identified KRAS related differentially expressed genes (DEGs) using The Cancer Genome Atlas (TCGA) database. A signature closely related to overall survival was recognized with Kaplan-Meier survival analysis and univariate cox regression analysis. Then we validated this signature with overall expression score (OE score) algorithm using both scRNA-seq and bulk RNA-seq data. Based on this signature, we performed LASSO cox regression to establish a prognostic model, and corresponding scores were calculated. Differences in genomic alteration, immune microenvironment, drug sensitivity between high- and low-KRD score groups were investigated. A KRAS related signature composed of 80 DEGs in colon cancer were recognized, among which 19 genes were selected to construct a prognostic model. This KRAS related signature was significantly correlated with worse prognosis. Furthermore, patients who scored lower in the prognostic model presented a higher likelihood of responding to chemotherapy, targeted therapy and immunotherapy. Furthermore, among the 19 selected genes in the model, SPINK4 was identified as an independent prognostic biomarker. Further validation in vitro indicated the knockdown of SPINK4 promoted the proliferation and migration of SW48 cells. In conclusion, a novel KRAS related signature was identified and validated based on clinical and genomic information from TCGA and GEO databases. The signature was proved to regulate genomic alteration, immune microenvironment and drug sensitivity in colon cancer, and thus might serve as a predictor for individual prognosis and treatment.

The Carcinogenic Potential of Bisphenol A in the Liver Based on Transcriptomic Studies

Bisphenol A (BPA) is an environmental toxin widely used in the production of polycarbonate plastics. A correlation exists between BPA tissue contamination and the occurrence of pathological conditions, including cancer. First-passage detoxification of high BPA amounts in the liver promotes hepatotoxicity and morphological alterations of this organ, but there is a lack of knowledge about the molecular mechanisms underlying these phenomena. This prompted us to investigate changes in the liver transcriptomics of 3-month-old female mice exposed to BPA (50 mg/kg) in drinking water for 3 months. Five female mice served as controls. The animals were euthanized, the livers were collected, and RNA was extracted to perform RNA-seq analysis. The multistep transcriptomic bioinformatics revealed 120 differentially expressed genes (DEGs) in the BPA-exposed samples. Gene Ontology (GO) annotations indicated that DEGs have been assigned to many biological processes, including "macromolecule modification" and "protein metabolic process". Several of the revealed DEGs have been linked to the pathogenesis of severe metabolic liver disorders and malignant tumors, in particular hepatocellular carcinoma. Data from this study suggest that BPA has a significant impact on gene expression in the liver, which is predictive of the carcinogenic potential of this compound in this organ.

An explainable AI-driven biomarker discovery framework for Non-Small Cell Lung Cancer classification

Non-Small Cell Lung Cancer (NSCLC) exhibits intrinsic heterogeneity at the molecular level that aids in distinguishing between its two prominent subtypes - Lung Adenocarcinoma (LUAD) and Lung Squamous Cell Carcinoma (LUSC). This paper proposes a novel explainable AI (XAI)-based deep learning framework to discover a small set of NSCLC biomarkers. The proposed framework comprises three modules - an autoencoder to shrink the input feature space, a feed-forward neural network to classify NSCLC instances into LUAD and LUSC, and a biomarker discovery module that leverages the combined network comprising the autoencoder and the feed-forward neural network. In the biomarker discovery module, XAI methods uncovered a set of 52 relevant biomarkers for NSCLC subtype classification. To evaluate the classification performance of the discovered biomarkers, multiple machine-learning models are constructed using these biomarkers. Using 10-Fold cross-validation, Multilayer Perceptron achieved an accuracy of 95.74% (±1.27) at 95% confidence interval. Further, using Drug-Gene Interaction Database, we observe that 14 of the discovered biomarkers are druggable. In addition, 28 biomarkers aid the prediction of the survivability of the patients. Out of 52 discovered biomarkers, we find that 45 biomarkers have been reported in previous studies on distinguishing between the two NSCLC subtypes. To the best of our knowledge, the remaining seven biomarkers have not yet been reported for NSCLC subtyping and could be further explored for their contribution to targeted therapy of lung cancer.

The genetic history of Mayotte and Madagascar cattle breeds mirrors the complex pattern of human exchanges in Western Indian Ocean

Despite their central economic and cultural role, the origin of cattle populations living in Indian Ocean islands still remains poorly documented. Here, we unravel the demographic and adaptive histories of the extant Zebus from the Mayotte and Madagascar islands using high-density SNP genotyping data. We found that these populations are very closely related and both display a predominant indicine ancestry. They diverged in the 16th century at the arrival of European people who transformed the trade network in the area. Their common ancestral cattle population originates from an admixture between an admixed African zebu population and an Indian zebu that occurred around the 12th century at the time of the earliest contacts between human African populations of the Swahili corridor and Austronesian people from Southeast Asia in Comoros and Madagascar. A steep increase in the estimated population sizes from the beginning of the 16th to the 17th century coincides with the expansion of the cattle trade. By carrying out genome scans for recent selection in the two cattle populations from Mayotte and Madagascar, we identified sets of candidate genes involved in biological functions (cancer, skin structure, and UV-protection, nervous system and behavior, organ development, metabolism, and immune response) broadly representative of the physiological adaptation to tropical conditions. Overall, the origin of the cattle populations from Western Indian Ocean islands mirrors the complex history of human migrations and trade in this area.

Structural and Biophysical Insights into SPINK1 Bound to Human Cationic Trypsin

(1) The serine protease inhibitor Kazal type 1 (SPINK1) inhibits trypsin activity in zymogen granules of pancreatic acinar cells. Several mutations in the SPINK1 gene are associated with acute recurrent pancreatitis (ARP) and chronic pancreatitis (CP). The most common variant is SPINK1 p.N34S. Although this mutation was identified two decades ago, the mechanism of action has remained elusive. (2) SPINK1 and human cationic trypsin (TRY1) were expressed in E. coli, and inhibitory activities were determined. Crystals of SPINK1-TRY1 complexes were grown by using the hanging-drop method, and phases were solved by molecular replacement. (3) Both SPINK1 variants show similar inhibitory behavior toward TRY1. The crystal structures are almost identical, with minor differences in the mutated loop. Both complexes show an unexpected rotamer conformation of the His63 residue in TRY1, which is a member of the catalytic triad. (4) The SPINK1 p.N34S mutation does not affect the inhibitory behavior or the overall structure of the protein. Therefore, the pathophysiological mechanism of action of the p.N34S variant cannot be explained mechanistically or structurally at the protein level. The observed histidine conformation is part of a mechanism for SPINK1 that can explain the exceptional proteolytic stability of this inhibitor.

Raised SPINK1 levels play a role in angiogenesis and the transendothelial migration of ALL cells

The present study was designed to assess whether raised Serine protease inhibitor Kazal type 1 (SPINK1) expressions modulates angiogenesis. Human umbilical vein endothelial cells (HUVECs) exposed to SPINK1 were noted to exhibit raised expressions of interleukin-8 (IL-8) as well as VCAM-1 and ICAM-1 cell adhesion molecules in a dose-dependent manner. In co-culture system of HUVECs and Acute lymphoblastic leukemia (ALL) cells, SPINK1 exposure also resulted in enhanced endothelial cell motility and ALL cells trans-endothelial migration. High concentrations of SPINK1 caused in vitro cellular reorganization into tubes in Matrigel-cultured HUVECs and induced in vivo vascularization and brain infiltration of NOD/SCID ALL model mice. The further transcriptomic analysis indicated that SPINK1 treatment altered several biological processes of endothelial cells and led to activation of the MAPK pathway. This study is the first to determine the neovascularization effects of raised SPINK1.

SPINKs in Tumors: Potential Therapeutic Targets

The serine protease inhibitor Kazal type (SPINK) family includes SPINK1-14 and is the largest branch in the serine protease inhibitor family. SPINKs play an important role in pancreatic physiology and disease, sperm maturation and capacitation, Nager syndrome, inflammation and the skin barrier. Evidence shows that the unregulated expression of SPINK1, 2, 4, 5, 6, 7, and 13 is closely related to human tumors. Different SPINKs exhibit various regulatory modes in different tumors and can be used as tumor prognostic markers. This article reviews the role of SPINK1, 2, 4, 5, 6, 7, and 13 in different human cancer processes and helps to identify new cancer treatment targets.