Algorithmically Reconstructed Molecular Pathways as the New Generation of Prognostic Molecular Biomarkers in Human Solid Cancers.

Individual gene expression and molecular pathway activation profiles were shown to be effective biomarkers in many cancers. Here, we used the human interactome model to algorithmically build 7470 molecular pathways centered around individual gene products. We assessed their associations with tumor type and survival in comparison with the previous generation of molecular pathway biomarkers (3022 "classical" pathways) and with the RNA transcripts or proteomic profiles of individual genes, for 8141 and 1117 samples, respectively. For all analytes in RNA and proteomic data, respectively, we found a total of 7441 and 7343 potential biomarker associations for gene-centric pathways, 3020 and 2950 for classical pathways, and 24,349 and 6742 for individual genes. Overall, the percentage of RNA biomarkers was statistically significantly higher for both types of pathways than for individual genes (p < 0.05). In turn, both types of pathways showed comparable performance. The percentage of cancer-type-specific biomarkers was comparable between proteomic and transcriptomic levels, but the proportion of survival biomarkers was dramatically lower for proteomic data. Thus, we conclude that pathway activation level is the advanced type of biomarker for RNA and proteomic data, and momentary algorithmic computer building of pathways is a new credible alternative to time-consuming hypothesis-driven manual pathway curation and reconstruction.

Uniformly shaped harmonization combines human transcriptomic data from different platforms while retaining their biological properties and differential gene expression patterns.

Introduction: Co-normalization of RNA profiles obtained using different experimental platforms and protocols opens avenue for comprehensive comparison of relevant features like differentially expressed genes associated with disease. Currently, most of bioinformatic tools enable normalization in a flexible format that depends on the individual datasets under analysis. Thus, the output data of such normalizations will be poorly compatible with each other. Recently we proposed a new approach to gene expression data normalization termed Shambhala which returns harmonized data in a uniform shape, where every expression profile is transformed into a pre-defined universal format. We previously showed that following shambhalization of human RNA profiles, overall tissue-specific clustering features are strongly retained while platform-specific clustering is dramatically reduced. Methods: Here, we tested Shambhala performance in retention of fold-change gene expression features and other functional characteristics of gene clusters such as pathway activation levels and predicted cancer drug activity scores. Results: Using 6,793 cancer and 11,135 normal tissue gene expression profiles from the literature and experimental datasets, we applied twelve performance criteria for different versions of Shambhala and other methods of transcriptomic harmonization with flexible output data format. Such criteria dealt with the biological type classifiers, hierarchical clustering, correlation/regression properties, stability of drug efficiency scores, and data quality for using machine learning classifiers. Discussion: Shambhala-2 harmonizer demonstrated the best results with the close to 1 correlation and linear regression coefficients for the comparison of training vs validation datasets and more than two times lesser instability for calculation of drug efficiency scores compared to other methods.

Macrophage - tumor cell interaction beyond cytokines.

Tumor cells communication with tumor associated macrophages is a highly important factor of tumor malignant potential development. For a long time, studies of this interaction were focused on a cytokine- and other soluble factors -mediated processes. Discovery of exosomes and regulatory RNAs as their cargo opened a broad field of research. Non-coding RNAs (ncRNAs) were demonstrated to contribute significantly to the development of macrophage phenotype, not only by regulating expression of certain genes, but also by providing for feedback loops of macrophage activation. Being a usual cargo of macrophage- or tumor cell-derived exosomes ncRNAs provide an important mechanism of tumor-stromal cell interaction that contributes significantly to the pathogenesis of various types of tumors. Despite the volume of ongoing research there are still many gaps that must be filled before the practical use of ncRNAs will be possible. In this review we discuss the role of regulatory RNAs in the development of macrophage phenotype. Further we review recent studies supporting the hypothesis that macrophages may affect the properties of tumor cells and vice versa tumor cells influence macrophage phenotype by miRNA and lncRNA transported between these cells by exosomes. We suggest that this mechanism of tumor cell - macrophage interaction is highly promising for the development of novel diagnostic and therapeutic strategies, though many problems are still to be solved.

Transcriptome-Based Traits of Radioresistant Sublines of Non-Small Cell Lung Cancer Cells.

Radioresistance is a major obstacle for the successful therapy of many cancers, including non-small cell lung cancer (NSCLC). To elucidate the mechanism of radioresistance of NSCLC cells and to identify key molecules conferring radioresistance, the radioresistant subclones of p53 wild-type A549 and p53-deficient H1299 cell cultures were established. The transcriptional changes between parental and radioresistant NSCLC cells were investigated by RNA-seq. In total, expression levels of 36,596 genes were measured. Changes in the activation of intracellular molecular pathways of cells surviving irradiation relative to parental cells were quantified using the Oncobox bioinformatics platform. Following 30 rounds of 2 Gy irradiation, a total of 322 genes were differentially expressed between p53 wild-type radioresistant A549IR and parental A549 cells. For the p53-deficient (H1299) NSCLC cells, the parental and irradiated populations differed in the expression of 1628 genes and 1616 pathways. The expression of genes associated with radioresistance reflects the complex biological processes involved in clinical cancer cell eradication and might serve as a potential biomarker and therapeutic target for NSCLC treatment.

Pan-cancer antagonistic inhibition pattern of ATM-driven G2/M checkpoint pathway vs other DNA repair pathways.

DNA repair mechanisms keep genome integrity and limit tumor-associated alterations and heterogeneity, but on the other hand they promote tumor survival after radiation and genotoxic chemotherapies. We screened pathway activation levels of 38 DNA repair pathways in nine human cancer types (gliomas, breast, colorectal, lung, thyroid, cervical, kidney, gastric, and pancreatic cancers). We took RNAseq profiles of the experimental 51 normal and 408 tumor samples, and from The Cancer Genome Atlas and Clinical Proteomic Tumor Analysis Consortium databases - of 500/407 normal and 5752/646 tumor samples, and also 573 normal and 984 tumor proteomic profiles from Proteomic Data Commons portal. For all the samplings we observed a congruent trend that all cancer types showed inhibition of G2/M arrest checkpoint pathway compared to the normal samples, and relatively low activities of p53-mediated pathways. In contrast, other DNA repair pathways were upregulated in most of the cancer types. The G2/M checkpoint pathway was statistically significantly downregulated compared to the other DNA repair pathways, and this inhibition was strongly impacted by antagonistic regulation of (i) promitotic genes CCNB and CDK1, and (ii) GADD45 genes promoting G2/M arrest. At the DNA level, we found that ATM, TP53, and CDKN1A genes accumulated loss of function mutations, and cyclin B complex genes - transforming mutations. These findings suggest importance of activation for most of DNA repair pathways in cancer progression, with remarkable exceptions of G2/M checkpoint and p53-related pathways which are downregulated and neutrally activated, respectively.

Personalized targeted therapy prescription in colorectal cancer using algorithmic analysis of RNA sequencing data.

BACKGROUND: Overall survival of advanced colorectal cancer (CRC) patients remains poor, and gene expression analysis could potentially complement detection of clinically relevant mutations to personalize CRC treatments. METHODS: We performed RNA sequencing of formalin-fixed, paraffin-embedded (FFPE) cancer tissue samples of 23 CRC patients and interpreted the data obtained using bioinformatic method Oncobox for expression-based rating of targeted therapeutics. Oncobox ranks cancer drugs according to the efficiency score calculated using target genes expression and molecular pathway activation data. The patients had primary and metastatic CRC with metastases in liver, peritoneum, brain, adrenal gland, lymph nodes and ovary. Two patients had mutations in NRAS, seven others had mutated KRAS gene. Patients were treated by aflibercept, bevacizumab, bortezomib, cabozantinib, cetuximab, crizotinib, denosumab, panitumumab and regorafenib as monotherapy or in combination with chemotherapy, and information on the success of totally 39 lines of therapy was collected. RESULTS: Oncobox drug efficiency score was effective biomarker that could predict treatment outcomes in the experimental cohort (AUC 0.77 for all lines of therapy and 0.91 for the first line after tumor sampling). Separately for bevacizumab, it was effective in the experimental cohort (AUC 0.87) and in 3 independent literature CRC datasets, n = 107 (AUC 0.84-0.94). It also predicted progression-free survival in univariate (Hazard ratio 0.14) and multivariate (Hazard ratio 0.066) analyses. Difference in AUC scores evidences importance of using recent biosamples for the prediction quality. CONCLUSION: Our results suggest that RNA sequencing analysis of tumor FFPE materials may be helpful for personalizing prescriptions of targeted therapeutics in CRC.

Experimentally Deduced Criteria for Detection of Clinically Relevant Fusion 3' Oncogenes from FFPE Bulk RNA Sequencing Data.

Drugs targeting receptor tyrosine kinase (RTK) oncogenic fusion proteins demonstrate impressive anti-cancer activities. The fusion presence in the cancer is the respective drug prescription biomarker, but their identification is challenging as both the breakpoint and the exact fusion partners are unknown. RNAseq offers the advantage of finding both fusion parts by screening sequencing reads. Paraffin (FFPE) tissue blocks are the most common way of storing cancer biomaterials in biobanks. However, finding RTK fusions in FFPE samples is challenging as RNA fragments are short and their artifact ligation may appear in sequencing libraries. Here, we annotated RNAseq reads of 764 experimental FFPE solid cancer samples, 96 leukemia samples, and 2 cell lines, and identified 36 putative clinically relevant RTK fusions with junctions corresponding to exon borders of the fusion partners. Where possible, putative fusions were validated by RT-PCR (confirmed for 10/25 fusions tested). For the confirmed 3'RTK fusions, we observed the following distinguishing features. Both moieties were in-frame, and the tyrosine kinase domain was preserved. RTK exon coverage by RNAseq reads upstream of the junction site were lower than downstream. Finally, most of the true fusions were present by more than one RNAseq read. This provides the basis for automatic annotation of 3'RTK fusions using FFPE RNAseq profiles.

Human Blood Serum Inhibits Ductal Carcinoma Cells BT474 Growth and Modulates Effect of HER2 Inhibition.

Trastuzumab, a HER2-targeted antibody, is widely used for targeted therapy of HER2-positive breast cancer (BC) patients; yet, not all of them respond to this treatment. We investigated here whether trastuzumab activity on the growth of HER2-overexpressing BT474 cells may interfere with human peripheral blood endogenous factors. Among 33 individual BC patient blood samples supplemented to the media, BT474 sensitivity to trastuzumab varied up to 14 times. In the absence of trastuzumab, human peripheral blood serum samples could inhibit growth of BT474, and this effect varied ~10 times for 50 individual samples. In turn, the epidermal growth factor (EGF) suppressed the trastuzumab effect on BT474 cell growth. Trastuzumab treatment increased the proportion of BT474 cells in the G0/G1 phases of cell cycle, while simultaneous addition of EGF decreased it, yet not to the control level. We used RNA sequencing profiling of gene expression to elucidate the molecular mechanisms involved in EGF- and human-sera-mediated attenuation of the trastuzumab effect on BT474 cell growth. Bioinformatic analysis of the molecular profiles suggested that trastuzumab acts similarly to the inhibition of PI3K/Akt/mTOR signaling axis, and the mechanism of EGF suppression of trastuzumab activity may be associated with parallel activation of PKC and transcriptional factors ETV1-ETV5.

Macrophage Phenotype in Combination with Tumor Microbiome Composition Predicts RCC Patients' Survival: A Pilot Study.

The identification of new prognostic markers of renal cell carcinoma (RCC) is an urgent problem in oncourology. To investigate the potential prognostic significance of tumor microbiome and stromal inflammatory markers, we studied a cohort of 66 patients with RCC (23 clear cell RCC, 19 papillary RCC and 24 chromophobe RCC). The microbiome was analyzed in tumor and normal tissue by 16S rRNA amplicon sequencing. Characterization of the tumor stroma was performed using immunohistochemistry. A significant difference in alpha diversity was demonstrated between normal kidney tissue and all types of RCC. Further, we demonstrated that the bacterial burden was higher in adjacent normal tissue than in a tumor. For the first time, we demonstrated a significant correlation between bacterial burden and the content of PU.1+ macrophages and CD66b+ neutrophils in kidney tumors. Tumors with high content of PU.1+ cells and CD66b+ cells in the stroma were characterized by a lower bacterial burden. In the tumors with high bacterial burden, the number of PU.1+ cells and CD66b+ was associated with a poor prognosis. The identified associations indicate the great prognostic potential of a combined tumor microbiome and stromal cell analysis.

Subtype of Neuroblastoma Cells with High KIT Expression Are Dependent on KIT and Its Knockdown Induces Compensatory Activation of Pro-Survival Signaling.

Neuroblastoma (NB) is a pediatric cancer with high clinical and molecular heterogeneity, and patients with high-risk tumors have limited treatment options. Receptor tyrosine kinase KIT has been identified as a potential marker of high-risk NB and a promising target for NB treatment. We investigated 19,145 tumor RNA expression and molecular pathway activation profiles for 20 cancer types and detected relatively high levels of KIT expression in NB. Increased KIT expression was associated with activation of cell survival pathways, downregulated apoptosis induction, and cell cycle checkpoint control pathways. KIT knockdown with shRNA encoded by lentiviral vectors in SH-SY5Y cells led to reduced cell proliferation and apoptosis induction up to 50%. Our data suggest that apoptosis induction was caused by mitotic catastrophe, and there was a 2-fold decrease in percentage of G2-M cell cycle phase after KIT knockdown. We found that KIT knockdown in NB cells leads to strong upregulation of other pro-survival growth factor signaling cascades such as EPO, NGF, IL-6, and IGF-1 pathways. NGF, IGF-1 and EPO were able to increase cell proliferation in KIT-depleted cells in an ERK1/2-dependent manner. Overall, we show that KIT is a promising therapeutic target in NB, although such therapy efficiency could be impeded by growth factor signaling activation.

Identification of cell type-specific correlations between ERK activity and cell viability upon treatment with ERK1/2 inhibitors.

Increased MAPK signaling is a hallmark of various cancers and is a central regulator of cell survival. Direct ERK1/2 inhibition is considered a promising approach to avoid ERK1/2 reactivation caused by upstream kinases BRAF, MEK1/2, and KRAS, as well as by receptor tyrosine kinase inhibitors, but the dynamics and selectivity of ERK1/2 inhibitors are much less studied compared with BRAF or MEK inhibitors. Using ERK1/2 and downstream kinase ELK1 reporter cell lines of lung cancer (H1299; NRASQ61K), colon cancer (HCT-116; KRASG13D), neuroblastoma (SH-SY5Y), and leukemia (U937), we examined the relationship between ERK inhibition and drug-induced toxicity for five ERK inhibitors: SCH772984, ravoxertinib, LY3214996, ulixertinib, and VX-11e, as well as one MEK inhibitor, PD0325901. Comparing cell viability and ERK inhibition revealed different ERK dependencies for these cell lines. We identify several drugs, such as SCH772984 and VX-11e, which induce excessive toxicity not directly related to ERK1/2 inhibition in specific cell lines. We also show that PD0325901, LY3214996, and ulixertinib are prone to ERK1/2 reactivation over time. We distinguished two types of ERK1/2 reactivation: the first could be reversed by adding a fresh dose of inhibitors, while the second persists even after additional treatments. We also showed that cells that became resistant to the MEK1/2 inhibitor PD0325901 due to ERK1/2 reactivation remained sensitive to ERK1/2 inhibitor ulixertinib. Our data indicate that correlation of ERK inhibition with drug-induced toxicity in multiple cell lines may help to find more selective and effective ERK1/2 inhibitors.

Next-Generation Grade and Survival Expression Biomarkers of Human Gliomas Based on Algorithmically Reconstructed Molecular Pathways.

In gliomas, expression of certain marker genes is strongly associated with survival and tumor type and often exceeds histological assessments. Using a human interactome model, we algorithmically reconstructed 7494 new-type molecular pathways that are centered each on an individual protein. Each single-gene expression and gene-centric pathway activation was tested as a survival and tumor grade biomarker in gliomas and their diagnostic subgroups (IDH mutant or wild type, IDH mutant with 1p/19q co-deletion, MGMT promoter methylated or unmethylated), including the three major molecular subtypes of glioblastoma (proneural, mesenchymal, classical). We used three datasets from The Cancer Genome Atlas and the Chinese Glioma Genome Atlas, which in total include 527 glioblastoma and 1097 low grade glioma profiles. We identified 2724 such gene and 2418 pathway survival biomarkers out of total 17,717 genes and 7494 pathways analyzed. We then assessed tumor grade and molecular subtype biomarkers and with the threshold of AUC > 0.7 identified 1322/982 gene biomarkers and 472/537 pathway biomarkers. This suggests roughly two times greater efficacy of the reconstructed pathway approach compared to gene biomarkers. Thus, we conclude that activation levels of algorithmically reconstructed gene-centric pathways are a potent class of new-generation diagnostic and prognostic biomarkers for gliomas.

OncoboxPD: human 51 672 molecular pathways database with tools for activity calculating and visualization.

OncoboxPD (Oncobox pathway databank) available at https://open.oncobox.com is the collection of 51 672 uniformly processed human molecular pathways. Superposition of all pathways formed interactome graph of protein-protein interactions and metabolic reactions containing 361 654 interactions and 64 095 molecular participants. Pathways are uniformly classified by biological processes, and each pathway node is algorithmically functionally annotated by specific activator/repressor role. This enables online calculation of statistically supported pathway activation levels (PALs) with the built-in bioinformatic tool using custom RNA/protein expression profiles. Each pathway can be visualized as static or dynamic graph, where vertices are molecules participating in a pathway and edges are interactions or reactions between them. Differentially expressed nodes in a pathway can be visualized in two-color mode with user-defined color scale. For every comparison, OncoboxPD also generates a graph summarizing top up- and downregulated pathways.

Transcriptomic Portraits and Molecular Pathway Activation Features of Adult Spinal Intramedullary Astrocytomas.

In this study, we report 31 spinal intramedullary astrocytoma (SIA) RNA sequencing (RNA-seq) profiles for 25 adult patients with documented clinical annotations. To our knowledge, this is the first clinically annotated RNA-seq dataset of spinal astrocytomas derived from the intradural intramedullary compartment. We compared these tumor profiles with the previous healthy central nervous system (CNS) RNA-seq data for spinal cord and brain and identified SIA-specific gene sets and molecular pathways. Our findings suggest a trend for SIA-upregulated pathways governing interactions with the immune cells and downregulated pathways for the neuronal functioning in the context of normal CNS activity. In two patient tumor biosamples, we identified diagnostic KIAA1549-BRAF fusion oncogenes, and we also found 16 new SIA-associated fusion transcripts. In addition, we bioinformatically simulated activities of targeted cancer drugs in SIA samples and predicted that several tyrosine kinase inhibitory drugs and thalidomide analogs could be potentially effective as second-line treatment agents to aid in the prevention of SIA recurrence and progression.

A Review of Recent Studies on the Effects of SARS-CoV-2 Infection and SARS-CoV-2 Vaccines on Male Reproductive Health.

According to world statistics, men are more susceptible to the coronavirus disease 2019 (COVID-19) than are women. Considering the interconnection between infections and male infertility, investigation of the potential impact of COVID-19 on men's reproductive health is now a particularly relevant topic. Published data indicate decreased sperm quality and orchitis development in patients with COVID-19, including reduced sperm count, decreased sperm motility, and elevated DNA fragmentation index. Although mass vaccination against COVID-19 is currently being carried out worldwide using available authorized vaccines, the effect of these vaccines on men's reproductive health has not yet been investigated. There is currently no evidence that SARS-CoV-2 can be transmitted in semen, but available data suggest that it can infect spermatogonia, spermatids, Leydig cells, and Sertoli cells. Therefore, SARS-CoV-2 orchitis and reduced male fertility may be long-term complications of COVID-19, which requires further investigation. Currently, there is also no evidence that vaccines against SARS-CoV-2 have any pathological effects on spermatogenesis or male reproductive health. Thus, further studies are needed to determine the effects of COVID-19 and COVID-19 vaccines on men's reproductive health, which will help to optimize the management and rehabilitation of these patients. This review aims to discuss recent studies on the impact of the COVID-19 and COVID-19 vaccines on men's reproductive health. The article addresses various issues such as the effect of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) on testosterone biosynthesis, semen parameters, testicular tissue, and epididymis.

Better Agreement of Human Transcriptomic and Proteomic Cancer Expression Data at the Molecular Pathway Activation Level.

Previously, we have shown that the aggregation of RNA-level gene expression profiles into quantitative molecular pathway activation metrics results in lesser batch effects and better agreement between different experimental platforms. Here, we investigate whether pathway level of data analysis provides any advantage when comparing transcriptomic and proteomic data. We compare the paired proteomic and transcriptomic gene expression and pathway activation profiles obtained for the same human cancer biosamples in The Cancer Genome Atlas (TCGA) and the NCI Clinical Proteomic Tumor Analysis Consortium (CPTAC) projects, for a total of 755 samples of glioblastoma, breast, liver, lung, ovarian, pancreatic, and uterine cancers. In a CPTAC assay, expression levels of 15,112 protein-coding genes were profiled using the Thermo QE series of mass spectrometers. In TCGA, RNA expression levels of the same genes were obtained using the Illumina HiSeq 4000 engine for the same biosamples. At the gene level, absolute gene expression values are compared, whereas pathway-grade comparisons are made between the pathway activation levels (PALs) calculated using average sample-normalized transcriptomic and proteomic profiles. We observed remarkably different average correlations between the primary RNA- and protein expression data for different cancer types: Spearman Rho between 0.017 (p = 1.7 × 10−13) and 0.27 (p < 2.2 × 10−16). However, at the pathway level we detected overall statistically significantly higher correlations: averaged Rho between 0.022 (p < 2.2 × 10−16) and 0.56 (p < 2.2 × 10−16). Thus, we conclude that data analysis at the PAL-level yields results of a greater similarity when comparing high-throughput RNA and protein expression profiles.

EGFR Pathway-Based Gene Signatures of Druggable Gene Mutations in Melanoma, Breast, Lung, and Thyroid Cancers.

EGFR, BRAF, PIK3CA, and KRAS genes play major roles in EGFR pathway, and accommodate activating mutations that predict response to many targeted therapeutics. However, connections between these mutations and EGFR pathway expression patterns remain unexplored. Here, we investigated transcriptomic associations with these activating mutations in three ways. First, we compared expressions of these genes in the mutant and wild type tumors, respectively, using RNA sequencing profiles from The Cancer Genome Atlas project database (n = 3660). Second, mutations were associated with the activation level of EGFR pathway. Third, they were associated with the gene signatures of differentially expressed genes from these pathways between the mutant and wild type tumors. We found that the upregulated EGFR pathway was linked with mutations in the BRAF (thyroid cancer, melanoma) and PIK3CA (breast cancer) genes. Gene signatures were associated with BRAF (thyroid cancer, melanoma), EGFR (squamous cell lung cancer), KRAS (colorectal cancer), and PIK3CA (breast cancer) mutations. However, only for the BRAF gene signature in the thyroid cancer we observed strong biomarker diagnostic capacity with AUC > 0.7 (0.809). Next, we validated this signature on the independent literature-based dataset (n = 127, fresh-frozen tissue samples, AUC 0.912), and on the experimental dataset (n = 42, formalin fixed, paraffin embedded tissue samples, AUC 0.822). Our results suggest that the RNA sequencing profiles can be used for robust identification of the replacement of Valine at position 600 with Glutamic acid in the BRAF gene in the papillary subtype of thyroid cancer, and evidence that the specific gene expression levels could provide information about the driver carcinogenic mutations.

RNA Sequencing Data for FFPE Tumor Blocks Can Be Used for Robust Estimation of Tumor Mutation Burden in Individual Biosamples.

Tumor mutation burden (TMB) is a well-known efficacy predictor for checkpoint inhibitor immunotherapies. Currently, TMB assessment relies on DNA sequencing data. Gene expression profiling by RNA sequencing (RNAseq) is another type of analysis that can inform clinical decision-making and including TMB estimation may strongly benefit this approach, especially for the formalin-fixed, paraffin-embedded (FFPE) tissue samples. Here, we for the first time compared TMB levels deduced from whole exome sequencing (WES) and RNAseq profiles of the same FFPE biosamples in single-sample mode. We took TCGA project data with mean sequencing depth 23 million gene-mapped reads (MGMRs) and found 0.46 (Pearson)-0.59 (Spearman) correlation with standard mutation calling pipelines. This was converted into low (<10) and high (>10) TMB per megabase classifier with area under the curve (AUC) 0.757, and application of machine learning increased AUC till 0.854. We then compared 73 experimental pairs of WES and RNAseq profiles with lower (mean 11 MGMRs) and higher (mean 68 MGMRs) RNA sequencing depths. For higher depth, we observed ~1 AUC for the high/low TMB classifier and 0.85 (Pearson)-0.95 (Spearman) correlation with standard mutation calling pipelines. For the lower depth, the AUC was below the high-quality threshold of 0.7. Thus, we conclude that using RNA sequencing of tumor materials from FFPE blocks with enough coverage can afford for high-quality discrimination of tumors with high and low TMB levels in a single-sample mode.

RNA Sequencing for Personalized Treatment of Metastatic Leiomyosarcoma: Case Report.

Uterine leiomyosarcoma (UL) is a rare malignant tumor that develops from the uterine smooth muscle tissue. Due to the low frequency and lack of sufficient data from clinical trials there is currently no effective treatment that is routinely accepted for UL. Here we report a case of a 65-years-old female patient with metastatic UL, who progressed on ifosfamide and doxorubicin therapy and developed severe hypertensive crisis after administration of second line pazopanib, which lead to treatment termination. Rapid progression of the tumor stressed the need for the alternative treatment options. We performed RNA sequencing and whole exome sequencing profiling of the patient's biopsy and applied Oncobox bioinformatic algorithm to prioritize targeted therapeutics. No clinically relevant mutations associated with drug efficiencies were found, but the Oncobox transcriptome analysis predicted regorafenib as the most effective targeted treatment option. Regorafenib administration resulted in a complete metabolic response which lasted for 10 months. In addition, RNA sequencing analysis revealed a novel cancer fusion transcript of YWHAE gene with fusion partner JAZF1. Several chimeric transcripts for YWHAE and JAZF1 genes were previously found in uterine neoplasms and some of them were associated with tumor prognosis. However, their combination was detected in this study for the first time. Taken together, these findings evidence that RNA sequencing may complement analysis of clinically relevant mutations and enhance management of oncological patients by suggesting putative treatment options.

Growth factor signaling predicts therapy resistance mechanisms and defines neuroblastoma subtypes.

Neuroblastoma (NB) has a low frequency of recurrent mutations compared to other cancers, which hinders the development of targeted therapies and novel risk stratification strategies. Multikinase inhibitors have shown potential in treating high-risk NB, but their efficacy is likely impaired by the cancer cells' ability to adapt to these drugs through the employment of alternative signaling pathways. Based on the expression of 48 growth factor-related genes in 1189 NB tumors, we have developed a model for NB patient survival prediction. This model discriminates between stage 4 NB tumors with favorable outcomes (>80% overall survival) and very poor outcomes (<10%) independently from MYCN-amplification status. Using signaling pathway analysis and gene set enrichment methods in 60 NB patients with known therapy response, we identified signaling pathways, including EPO, NGF, and HGF, upregulated in patients with no or partial response. In a therapeutic setting, we showed that among six selected growth factors, EPO, and NGF showed the most pronounced protective effects in vitro against several promising anti-NB multikinase inhibitors: imatinib, dasatinib, crizotinib, cabozantinib, and axitinib. Mechanistically kinase inhibitors potentiated NB cells to stronger ERK activation by EPO and NGF. The protective action of these growth factors strongly correlated with ERK activation and was ERK-dependent. ERK inhibitors combined with anticancer drugs, especially with dasatinib, showed a synergistic effect on NB cell death. Consideration of growth factor signaling activity benefits NB outcome prediction and tailoring therapy regimens to treat NB.