Growing evidence of transcriptional and metabolomic differentiation induced many studies which analyze such differentiation in context of outcome of disease progression, treatment or influence of many different factors affecting cellular and tissue metabolism. Particularly, cancer researchers are looking for new biomarkers that can serve as a diagnostic/prognostic factor and its further corresponding relationship regarding clinical effects. As a result of the increasing interest in use of dichotomization of continuous variables involving clinical or epidemiological data (gene expression, biomarkers, biochemical parameters, etc.) there is a large demand for cutoff point determination tools with simultaneous lack of software offering stratification of patients based on continuous and binary variables. Therefore, we developed "Evaluate Cutpoints" application offering wide set of statistical and graphical methods for cutpoint optimization enabling stratification of population into two or three groups. Application is based on R language including algorithms of packages such as survival, survMisc, OptimalCutpoints, maxstat, Rolr, ggplot2, GGally and plotly offering Kaplan-Meier plots and ROC curves with cutoff point determination. All capabilities of Evaluate Cutpoints were illustrated with example analysis of estrogen, progesterone and human epidermal growth factor 2 receptors in breast cancer cohort. Through ROC curve the cutoff points were established for expression of ESR1, PGR and ERBB2 in correlation with their immunohistochemical status (cutoff: 1301.253, 243.35, 11,434.438, respectively; sensitivity: 94%, 85%, 64%, respectively; specificity: 93%, 86%, 91%, respectively). Through disease-free survival analysis we divided patients into two and three groups regarding expression of ESR1, PGR and ERBB2. Example algorithm cutp showed that lowered expression of ESR1 and ERBB2 was more favorable (HR = 2.07, p = 0.0412; HR = 2.79, p = 0.0777, respectively), whereas heightened PGR expression was correlated with better prognosis (HR = 0.192, p = 0.0115). This work presents application Evaluate Cutpoints that is freely available to download at http://wnbikp.umed.lodz.pl/Evaluate-Cutpoints/. Currently, many softwares are used to split continuous variables such as Cutoff Finder and X-Tile, which offer distinct algorithms. Unlike them, Evaluate Cutpoints allows not only dichotomization of populations into groups according to continuous variables and binary variables, but also stratification into three groups as well as manual selection of cutoff point thus preventing potential loss of information.

Lung malignancies comprise lethal and aggressive tumours that remain the leading cancer-related death cause worldwide. Regarding histological classification, lung squamous cell carcinoma (LUSC) and adenocarcinoma (LUAD) account for the majority of cases. Surgical resection and various combinations of chemo- and radiation therapies are the golden standards in the treatment of lung cancers, although the five-year survival rate remains very poor. Notch, Hedgehog, Wnt and Erbb signalling are evolutionarily conserved pathways regulating pivotal cellular processes such as differentiation, proliferation, and angiogenesis during embryogenesis and post-natal life. However, to date, there is no study comprehensively revealing signalling networks of these four pathways in LUSC and LUAD. Therefore, the aim of the present study was the investigation profiles of downstream target genes of pathways that differ between LUSC and LUAD biology. Our results showed a few co-expression modules, identified through weighted gene co-expression network analysis (WGCNA), which significantly differentiated downstream signaling of Notch, ErbB, Hedgehog, and Wnt in LUSC and LUAD. Among co-expressed genes essential regulators of the cell cycle, DNA damage response, apoptosis, and proliferation have been found. Most of them were upregulated in LUSC compared to LUAD. In conclusion, identified downstream networks revealed distinct biological mechanisms underlying cancer development and progression in LUSC and LUAD that may diversify the clinical outcome of the disease.

Notch signaling is an evolutionarily conserved pathway regulating normal embryonic development and homeostasis in a wide variety of tissues. It is also critically involved in carcinogenesis, as well as cancer progression. Activation of the Notch pathway members can be either oncogenic or suppressive, depending on tissue context. The present study is a comprehensive overview, extended with a bioinformatics analysis of TCGA cohorts, including breast, bladder, cervical, colon, kidney, lung, ovary, prostate and rectum carcinomas. We performed global expression profiling of the Notch pathway core components and downstream targets. For this purpose, we implemented the Uniform Manifold Approximation and Projection algorithm to reduce the dimensions. Furthermore, we determined the optimal cutpoint using Evaluate Cutpoint software to established disease-free and overall survival with respect to particular Notch members. Our results demonstrated separation between tumors and their corresponding normal tissue, as well as between tumors in general. The differentiation of the Notch pathway, at its various stages, in terms of expression and survival resulted in distinct profiles of biological processes such as proliferation, adhesion, apoptosis and epithelial to mesenchymal transition. In conclusion, whether oncogenic or suppressive, Notch signaling is proven to be associated with various types of malignancies, and thus may be of interest as a potential therapeutic target.

Prostate carcinoma (PRAD) is one of the most frequently diagnosed malignancies amongst men worldwide. It is well-known that androgen receptor (AR) plays a pivotal role in a vast majority of prostate tumors. However, recent evidence emerged stating that estrogen receptors (ERs) may also contribute to prostate tumor development. Moreover, progression and aggressiveness of prostate cancer may be associated with differential expression genes of epithelial-to-mesenchymal transition (EMT). Therefore we aimed to assess the significance of receptors status as well as EMT marker genes expression among PRAD patients in accordance to their age and Gleason score. We analyzed TCGA gene expression profiles of 497 prostate tumor samples according to 43 genes involved in EMT and 3 hormone receptor genes (AR, ESR1, ESR2) as well as clinical characteristic of cancer patients. Then patients were divided into four groups according to their age and 5 groups according to Gleason score. Next, we evaluated PRAD samples according to relationship between the set of variables in different combinations and compared differential expression in subsequent groups of patients. The analysis was applied using R packages: FactoMineR, gplots, RColorBrewer and NMF. MFA analysis resulted in distinct grouping of PRAD patients into four age categories according to expression level of AR, ESR1 and ESR2 with the most distinct group of age less than 50 years old. Further investigations indicated opposite expression profiles of EMT markers between different age groups as well as strong association of EMT gene expression with Gleason score. We found that depending on age of prostate cancer patients and Gleason score EMT genes with distinctly altered expression are: KRT18, KRT19, MUC1 and COL4A1, CTNNB1, SNAI2, ZEB1 and MMP3. Our major observation is that prostate cancer from patients under 50 years old compared to older ones has entirely different EMT gene expression profiles showing potentially more aggressive invasive phenotype, despite Gleason score classification.

Most of the currently used cancer immunotherapies inhibit the programmed cell death protein 1 (PD1)–programmed cell death 1 ligand 1 (PDL1) axis of T‐cells. However, dendritic cells (DCs) controlled by natural killer (NK) cells via the FMS‐related tyrosine kinase 3 (FLT3) axis are necessary for activation of T‐cells. The aim of the study was to evaluate FLT3 as a prognostic factor and determine its role in immune infiltration (with emphasis on NK cells and DCs). Using The Cancer Genome Atlas (TCGA) database, we performed bioinformatic analysis of the gene expression datasets of 501 lung squamous cell carcinoma (LUSC) and 515 lung adenocarcinoma (LUAD) patient who had corresponding clinical data [analysis was performed in R (version 4.2.0)]. Disease‐free survival (DFS) differed between the FLT3‐low and FLT3‐high expression groups, respectively, in LUSC (61.0 vs 71.3 months P = 0.075) and LUAD (32.7 vs 47.5 months P = 0.045). A tumor microenvironment (TME) with high immune infiltration and rich in T‐cell exhaustion markers was observed in the FLT3‐high group. We showed overexpression of NK cell and DC gene signatures in the FLT3‐high expression group as well as overexpression of key effector genes of the cyclic GMP‐AMP synthase (cGAS)–stimulator of interferon genes protein (STING) pathway, which is crucial in response to radiotherapy. High expression of FLT3 in the TME was associated with immune cell infiltration (especially of NK cells and DCs), increased expression of T‐cell exhaustion markers and expression of effector genes of the cGAS‐STING pathway, which may consequently increase susceptibility to immunotherapy and radiotherapy. High FLT3 expression correlated with prolonged DFS in the LUSC and LUAD cohorts.

The brain is the most functionally organized structure of all organs. It manages behavior, perception and higher cognitive functions. The WWOX gene is non-classical tumor suppressor gene, which has been shown to have an impact on proliferation, apoptosis and migration processes. Moreover, genetic aberrations in WWOX induce severe neuropathological phenotypes in humans and rodents. The aim of the present study was to investigate in detail the impact of WWOX on human neural progenitor cell (hNPC) maintenance and how depletion of WWOX disturbs signaling pathways playing a pivotal role in neuronal differentiation and central nervous system (CNS) organogenesis. hNPC with a silenced WWOX gene exhibited lowered mitochondrial redox potential, enhanced adhesion to fibronectin and extracellular matrix protein mixture, downregulation of MMP2/9 expression and impaired 3D growth. Global transcriptome analysis using cap analysis of gene expression (CAGE) found that WWOX downregulation significantly changes the expression of multiple genes engaged in cytoskeleton organization, adhesion, cell signaling and chromatin remodeling. The massive changes in gene expression caused by WWOX silencing may strongly affect the differentiation and migration of neurons in organogenesis, brain injury, cancerogenesis or neurodifferentiation. WWOX gene appears to be an important regulator of neural tissue architecture and function.