1
|
Borisov N, Tkachev V, Simonov A, Sorokin M, Kim E, Kuzmin D, Karademir-Yilmaz B, Buzdin A. Uniformly shaped harmonization combines human transcriptomic data from different platforms while retaining their biological properties and differential gene expression patterns. Front Mol Biosci 2023; 10:1237129. [PMID: 37745690 PMCID: PMC10511763 DOI: 10.3389/fmolb.2023.1237129] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/08/2023] [Accepted: 08/28/2023] [Indexed: 09/26/2023] Open
Abstract
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.
Collapse
Affiliation(s)
- Nicolas Borisov
- Omicsway Corp, Walnut, CA, United States
- Moscow Institute of Physics and Technology, Dolgoprudny, Russia
| | | | - Alexander Simonov
- Moscow Institute of Physics and Technology, Dolgoprudny, Russia
- Oncobox Ltd., Moscow, Russia
| | - Maxim Sorokin
- Moscow Institute of Physics and Technology, Dolgoprudny, Russia
- Oncobox Ltd., Moscow, Russia
- World-Class Research Center “Digital Biodesign and Personalized Healthcare”, Sechenov First Moscow State Medical University, Moscow, Russia
| | - Ella Kim
- Clinic for Neurosurgery, Laboratory of Experimental Neurooncology, Johannes Gutenberg University Medical Centre, Mainz, Germany
| | - Denis Kuzmin
- Moscow Institute of Physics and Technology, Dolgoprudny, Russia
| | - Betul Karademir-Yilmaz
- Department of Biochemistry, School of Medicine/Genetic and Metabolic Diseases Research and Investigation Center (GEMHAM) Marmara University, Istanbul, Türkiye
| | - Anton Buzdin
- Moscow Institute of Physics and Technology, Dolgoprudny, Russia
- World-Class Research Center “Digital Biodesign and Personalized Healthcare”, Sechenov First Moscow State Medical University, Moscow, Russia
- Shemyakin-Ovchinnikov Institute of Bioorganic Chemistry, Moscow, Russia
- PathoBiology Group, European Organization for Research and Treatment of Cancer (EORTC), Brussels, Belgium
| |
Collapse
|
2
|
Sorokin M, Zolotovskaia M, Nikitin D, Suntsova M, Poddubskaya E, Glusker A, Garazha A, Moisseev A, Li X, Sekacheva M, Naskhletashvili D, Seryakov A, Wang Y, Buzdin A. Personalized targeted therapy prescription in colorectal cancer using algorithmic analysis of RNA sequencing data. BMC Cancer 2022; 22:1113. [PMCID: PMC9623986 DOI: 10.1186/s12885-022-10177-3] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/27/2021] [Accepted: 09/26/2022] [Indexed: 12/24/2022] Open
Abstract
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. Supplementary Information The online version contains supplementary material available at 10.1186/s12885-022-10177-3.
Collapse
Affiliation(s)
- Maxim Sorokin
- grid.448878.f0000 0001 2288 8774I.M. Sechenov First Moscow State Medical University, 119991 Moscow, Russia ,grid.18763.3b0000000092721542Moscow Institute of Physics and Technology, 141701 Moscow Region, Russia ,OmicsWay Corp, 91789 Walnut, CA USA
| | - Marianna Zolotovskaia
- grid.18763.3b0000000092721542Moscow Institute of Physics and Technology, 141701 Moscow Region, Russia
| | - Daniil Nikitin
- OmicsWay Corp, 91789 Walnut, CA USA ,grid.418853.30000 0004 0440 1573Shemyakin-Ovchinnikov Institute of Bioorganic Chemistry, 117997 Moscow, Russia
| | - Maria Suntsova
- grid.448878.f0000 0001 2288 8774World-Class Research Center “Digital biodesign and personalized healthcare”, Sechenov First Moscow State Medical University, 119991 Moscow, Russia
| | - Elena Poddubskaya
- grid.448878.f0000 0001 2288 8774I.M. Sechenov First Moscow State Medical University, 119991 Moscow, Russia ,Clinical Center Vitamed, 121309 Moscow, Russia
| | - Alexander Glusker
- grid.448878.f0000 0001 2288 8774I.M. Sechenov First Moscow State Medical University, 119991 Moscow, Russia
| | | | - Alexey Moisseev
- grid.448878.f0000 0001 2288 8774I.M. Sechenov First Moscow State Medical University, 119991 Moscow, Russia
| | - Xinmin Li
- grid.19006.3e0000 0000 9632 6718Department of Pathology and Laboratory Medicine, University of California, 90095 Los Angeles, CA USA
| | - Marina Sekacheva
- grid.448878.f0000 0001 2288 8774World-Class Research Center “Digital biodesign and personalized healthcare”, Sechenov First Moscow State Medical University, 119991 Moscow, Russia
| | - David Naskhletashvili
- grid.466904.90000 0000 9092 133XN.N. Blokhin Russian Cancer Research Center, 115478 Moscow, Russia
| | | | - Ye Wang
- grid.410645.20000 0001 0455 0905Core Laboratory, The Affiliated Qingdao Central Hospital of Qingdao University, Qingdao, China
| | - Anton Buzdin
- grid.18763.3b0000000092721542Moscow Institute of Physics and Technology, 141701 Moscow Region, Russia ,grid.418853.30000 0004 0440 1573Shemyakin-Ovchinnikov Institute of Bioorganic Chemistry, 117997 Moscow, Russia ,grid.448878.f0000 0001 2288 8774World-Class Research Center “Digital biodesign and personalized healthcare”, Sechenov First Moscow State Medical University, 119991 Moscow, Russia
| |
Collapse
|
3
|
Transcriptomic Harmonization as the Way for Suppressing Cross-Platform Bias and Batch Effect. Biomedicines 2022; 10:biomedicines10092318. [PMID: 36140419 PMCID: PMC9496268 DOI: 10.3390/biomedicines10092318] [Citation(s) in RCA: 4] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/22/2022] [Revised: 09/14/2022] [Accepted: 09/16/2022] [Indexed: 11/16/2022] Open
Abstract
(1) Background: Emergence of methods interrogating gene expression at high throughput gave birth to quantitative transcriptomics, but also posed a question of inter-comparison of expression profiles obtained using different equipment and protocols and/or in different series of experiments. Addressing this issue is challenging, because all of the above variables can dramatically influence gene expression signals and, therefore, cause a plethora of peculiar features in the transcriptomic profiles. Millions of transcriptomic profiles were obtained and deposited in public databases of which the usefulness is however strongly limited due to the inter-comparison issues; (2) Methods: Dozens of methods and software packages that can be generally classified as either flexible or predefined format harmonizers have been proposed, but none has become to the date the gold standard for unification of this type of Big Data; (3) Results: However, recent developments evidence that platform/protocol/batch bias can be efficiently reduced not only for the comparisons of limited transcriptomic datasets. Instead, instruments were proposed for transforming gene expression profiles into the universal, uniformly shaped format that can support multiple inter-comparisons for reasonable calculation costs. This forms a basement for universal indexing of all or most of all types of RNA sequencing and microarray hybridization profiles; (4) Conclusions: In this paper, we attempted to overview the landscape of modern approaches and methods in transcriptomic harmonization and focused on the practical aspects of their application.
Collapse
|
4
|
Borisov N, Sorokin M, Zolotovskaya M, Borisov C, Buzdin A. Shambhala-2: A Protocol for Uniformly Shaped Harmonization of Gene Expression Profiles of Various Formats. Curr Protoc 2022; 2:e444. [PMID: 35617464 DOI: 10.1002/cpz1.444] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/11/2022]
Abstract
Uniformly shaped harmonization of gene expression profiles is central for the simultaneous comparison of multiple gene expression datasets. It is expected to operate with the gene expression data obtained using various experimental methods and equipment, and to return harmonized profiles in a uniform shape. Such uniformly shaped expression profiles from different initial datasets can be further compared directly. However, current harmonization techniques have strong limitations that prevent their broad use for bioinformatic applications. They can either operate with only up to two datasets/platforms or return data in a dynamic format that will be different for every comparison under analysis. This also does not allow for adding new data to the previously harmonized dataset(s), which complicates the analysis and increases calculation costs. We propose here a new method termed Shambhala-2 that can transform multi-platform expression data into a universal format that is identical for all harmonizations made using this technique. Shambhala-2 is based on sample-by-sample cubic conversion of the initial expression dataset into a preselected shape of the reference definitive dataset. Using 8390 samples of 12 healthy human tissue types and 4086 samples of colorectal, kidney, and lung cancer tissues, we verified Shambhala-2's capacity in restoring tissue-specific expression patterns for seven microarray and three RNA sequencing platforms. Shambhala-2 performed well for all tested combinations of RNAseq and microarray profiles, and retained gene-expression ranks, as evidenced by high correlations between different single- or aggregated gene expression metrics in pre- and post-Shambhalized samples, including preserving cancer-specific gene expression and pathway activation features. © 2022 Wiley Periodicals LLC. Basic Protocol: Shambhala-2 harmonizer Alternate Protocol 1: Linear Shambhala/Shambhala-1 Alternate Protocol 2: Alternative (flexible-format and uniformly shaped) normalization methods Support Protocol 1: Watermelon multisection (WM) Support Protocol 2: Calculation of cancer-to-normal log-fold-change (LFC) and pathway activation level (PAL).
Collapse
Affiliation(s)
- Nicolas Borisov
- Omicsway Corp., Walnut, California.,Moscow Institute of Physics and Technology, Dolgoprudny, Moscow Region, Russia
| | - Maksim Sorokin
- Omicsway Corp., Walnut, California.,Moscow Institute of Physics and Technology, Dolgoprudny, Moscow Region, Russia.,I.M. Sechenov First Moscow State Medical University, Moscow, Russia
| | - Marianna Zolotovskaya
- Moscow Institute of Physics and Technology, Dolgoprudny, Moscow Region, Russia.,Oncobox Ltd., Moscow, Russia
| | | | - Anton Buzdin
- Moscow Institute of Physics and Technology, Dolgoprudny, Moscow Region, Russia.,Shemyakin-Ovchinnikov Institute of Bioorganic Chemistry, Moscow, Russia.,World-Class Research Center "Digital biodesign and personalized healthcare", Sechenov First Moscow State Medical University, Moscow, Russia.,PathoBiology Group, European Organization for Research and Treatment of Cancer (EORTC), Brussels, Belgium
| |
Collapse
|
5
|
Seryakov A, Magomedova Z, Suntsova M, Prokofieva A, Rabushko E, Glusker A, Makovskaia L, Zolotovskaia M, Buzdin A, Sorokin M. RNA Sequencing for Personalized Treatment of Metastatic Leiomyosarcoma: Case Report. Front Oncol 2021; 11:666001. [PMID: 34527573 PMCID: PMC8435728 DOI: 10.3389/fonc.2021.666001] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/15/2021] [Accepted: 08/11/2021] [Indexed: 01/14/2023] Open
Abstract
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.
Collapse
Affiliation(s)
| | - Zaynab Magomedova
- The Laboratory of Clinical and Genomic Bioinformatics, I.M. Sechenov First Moscow State Medical University, Moscow, Russia
| | - Maria Suntsova
- World-Class Research Center “Digital Biodesign and Personalized Healthcare”, Sechenov First Moscow State Medical University, Moscow, Russia
| | - Anastasia Prokofieva
- The Laboratory of Clinical and Genomic Bioinformatics, I.M. Sechenov First Moscow State Medical University, Moscow, Russia
| | - Elizaveta Rabushko
- The Laboratory of Clinical and Genomic Bioinformatics, I.M. Sechenov First Moscow State Medical University, Moscow, Russia
| | - Alexander Glusker
- The Laboratory of Clinical and Genomic Bioinformatics, I.M. Sechenov First Moscow State Medical University, Moscow, Russia
| | - Lyudmila Makovskaia
- Faculty of Fundamental Medicine, Lomonosov Moscow State University, Moscow, Russia
| | - Marianna Zolotovskaia
- Laboratory of Translational Genomic Bioinformatics, Moscow Institute of Physics and Technology, Dolgoprudny, Russia
| | - Anton Buzdin
- World-Class Research Center “Digital Biodesign and Personalized Healthcare”, Sechenov First Moscow State Medical University, Moscow, Russia
- Laboratory of Translational Genomic Bioinformatics, Moscow Institute of Physics and Technology, Dolgoprudny, Russia
- Shemyakin-Ovchinnikov Institute of Bioorganic Chemistry, Russian Academy of Sciences, Moscow, Russia
- OmicsWay Corp, Walnut, CA, United States
| | - Maxim Sorokin
- The Laboratory of Clinical and Genomic Bioinformatics, I.M. Sechenov First Moscow State Medical University, Moscow, Russia
- Laboratory of Translational Genomic Bioinformatics, Moscow Institute of Physics and Technology, Dolgoprudny, Russia
- OmicsWay Corp, Walnut, CA, United States
| |
Collapse
|
6
|
Using proteomic and transcriptomic data to assess activation of intracellular molecular pathways. ADVANCES IN PROTEIN CHEMISTRY AND STRUCTURAL BIOLOGY 2021; 127:1-53. [PMID: 34340765 DOI: 10.1016/bs.apcsb.2021.02.005] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/13/2022]
Abstract
Analysis of molecular pathway activation is the recent instrument that helps to quantize activities of various intracellular signaling, structural, DNA synthesis and repair, and biochemical processes. This may have a deep impact in fundamental research, bioindustry, and medicine. Unlike gene ontology analyses and numerous qualitative methods that can establish whether a pathway is affected in principle, the quantitative approach has the advantage of exactly measuring the extent of a pathway up/downregulation. This results in emergence of a new generation of molecular biomarkers-pathway activation levels, which reflect concentration changes of all measurable pathway components. The input data can be the high-throughput proteomic or transcriptomic profiles, and the output numbers take both positive and negative values and positively reflect overall pathway activation. Due to their nature, the pathway activation levels are more robust biomarkers compared to the individual gene products/protein levels. Here, we review the current knowledge of the quantitative gene expression interrogation methods and their applications for the molecular pathway quantization. We consider enclosed bioinformatic algorithms and their applications for solving real-world problems. Besides a plethora of applications in basic life sciences, the quantitative pathway analysis can improve molecular design and clinical investigations in pharmaceutical industry, can help finding new active biotechnological components and can significantly contribute to the progressive evolution of personalized medicine. In addition to the theoretical principles and concepts, we also propose publicly available software for the use of large-scale protein/RNA expression data to assess the human pathway activation levels.
Collapse
|
7
|
Sorokin M, Ignatev K, Poddubskaya E, Vladimirova U, Gaifullin N, Lantsov D, Garazha A, Allina D, Suntsova M, Barbara V, Buzdin A. RNA Sequencing in Comparison to Immunohistochemistry for Measuring Cancer Biomarkers in Breast Cancer and Lung Cancer Specimens. Biomedicines 2020; 8:E114. [PMID: 32397474 PMCID: PMC7277916 DOI: 10.3390/biomedicines8050114] [Citation(s) in RCA: 20] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/29/2020] [Revised: 05/02/2020] [Accepted: 05/07/2020] [Indexed: 12/11/2022] Open
Abstract
RNA sequencing is considered the gold standard for high-throughput profiling of gene expression at the transcriptional level. Its increasing importance in cancer research and molecular diagnostics is reflected in the growing number of its mentions in scientific literature and clinical trial reports. However, the use of different reagents and protocols for RNA sequencing often produces incompatible results. Recently, we published the Oncobox Atlas of RNA sequencing profiles for normal human tissues obtained from healthy donors killed in road accidents. This is a database of molecular profiles obtained using uniform protocol and reagents settings that can be broadly used in biomedicine for data normalization in pathology, including cancer. Here, we publish new original 39 breast cancer (BC) and 19 lung cancer (LC) RNA sequencing profiles obtained for formalin-fixed paraffin-embedded (FFPE) tissue samples, fully compatible with the Oncobox Atlas. We performed the first correlation study of RNA sequencing and immunohistochemistry-measured expression profiles for the clinically actionable biomarker genes in FFPE cancer tissue samples. We demonstrated high (Spearman's rho 0.65-0.798) and statistically significant (p < 0.00004) correlations between the RNA sequencing (Oncobox protocol) and immunohistochemical measurements for HER2/ERBB2, ER/ESR1 and PGR genes in BC, and for PDL1 gene in LC; AUC: 0.963 for HER2, 0.921 for ESR1, 0.912 for PGR, and 0.922 for PDL1. To our knowledge, this is the first validation that total RNA sequencing of archived FFPE materials provides a reliable estimation of marker protein levels. These results show that in the future, RNA sequencing can complement immunohistochemistry for reliable measurements of the expression biomarkers in FFPE cancer samples.
Collapse
Affiliation(s)
- Maxim Sorokin
- Institute of Personalized Medicine, I.M. Sechenov First Moscow State Medical University, 119048 Moscow, Russia; (M.S.); (E.P.); (D.A.); (M.S.)
- Omicsway Corp., Walnut, CA 91789, USA;
- Shemyakin-Ovchinnikov Institute of Bioorganic Chemistry, 117997 Moscow, Russia;
| | - Kirill Ignatev
- Karelia Republic Oncological Hospital, 185000 Petrozavodsk, Russia;
| | - Elena Poddubskaya
- Institute of Personalized Medicine, I.M. Sechenov First Moscow State Medical University, 119048 Moscow, Russia; (M.S.); (E.P.); (D.A.); (M.S.)
- Vitamed Oncological Clinical Center, 121309 Moscow, Russia
| | - Uliana Vladimirova
- Shemyakin-Ovchinnikov Institute of Bioorganic Chemistry, 117997 Moscow, Russia;
| | - Nurshat Gaifullin
- Faculty of Fundamental Medicine, Lomonosov Moscow State University, 119991 Moscow, Russia;
| | - Dmitriy Lantsov
- Kaluga Regional Oncological Hospital, 248007 Kaluga, Russia;
| | | | - Daria Allina
- Institute of Personalized Medicine, I.M. Sechenov First Moscow State Medical University, 119048 Moscow, Russia; (M.S.); (E.P.); (D.A.); (M.S.)
| | - Maria Suntsova
- Institute of Personalized Medicine, I.M. Sechenov First Moscow State Medical University, 119048 Moscow, Russia; (M.S.); (E.P.); (D.A.); (M.S.)
| | - Victoria Barbara
- Oncological Dispensary of the Republic of Karelia, 185002 Petrozavodsk, Russia;
| | - Anton Buzdin
- Institute of Personalized Medicine, I.M. Sechenov First Moscow State Medical University, 119048 Moscow, Russia; (M.S.); (E.P.); (D.A.); (M.S.)
- Omicsway Corp., Walnut, CA 91789, USA;
- Shemyakin-Ovchinnikov Institute of Bioorganic Chemistry, 117997 Moscow, Russia;
- Moscow Institute of Physics and Technology, 141701 Moscow, Russia
| |
Collapse
|
8
|
Sorokin M, Poddubskaya E, Baranova M, Glusker A, Kogoniya L, Markarova E, Allina D, Suntsova M, Tkachev V, Garazha A, Sekacheva M, Buzdin A. RNA sequencing profiles and diagnostic signatures linked with response to ramucirumab in gastric cancer. Cold Spring Harb Mol Case Stud 2020; 6:a004945. [PMID: 32060041 PMCID: PMC7133748 DOI: 10.1101/mcs.a004945] [Citation(s) in RCA: 20] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/06/2019] [Accepted: 02/03/2020] [Indexed: 02/06/2023] Open
Abstract
Gastric cancer (GC) is the fifth-ranked cancer type by associated mortality. The proportion of early diagnosis is low, and most patients are diagnosed at the advanced stages. First-line therapy standardly includes fluoropyrimidines and platinum compounds with trastuzumab for HER2-positive cases. For recurrent disease, there are several alternative options including ramucirumab, a monoclonal therapeutic antibody that inhibits VEGF-mediated tumor angiogenesis by binding with VEGFR2, alone or in combination with other cancer drugs. However, overall response rate following ramucirumab or its combinations is 30%-80% of the patients, suggesting that personalization of drug prescription is needed to increase efficacy of treatment. We report here original tumor RNA sequencing profiles for 15 advanced GC patients linked with data on clinical response to ramucirumab or its combinations. Three genes showed differential expression in the tumors for responders versus nonresponders: CHRM3, LRFN1, and TEX15 Of them, CHRM3 was up-regulated in the responders. Using the bioinformatic platform Oncobox we simulated ramucirumab efficiency and compared output model results with actual tumor response data. An agreement was observed between predicted and real clinical outcomes (AUC ≥ 0.7). These results suggest that RNA sequencing may be used to personalize the prescription of ramucirumab for GC and indicate potential molecular mechanisms underlying ramucirumab resistance. The RNA sequencing profiles obtained here are fully compatible with the previously published Oncobox Atlas of Normal Tissue Expression (ANTE) data.
Collapse
Affiliation(s)
- Maxim Sorokin
- I.M. Sechenov First Moscow State Medical University, Moscow, 119991, Russia
- Omicsway Corp., Walnut, California 91789, USA
- Shemyakin-Ovchinnikov Institute of Bioorganic Chemistry, Moscow, 117997, Russia
| | - Elena Poddubskaya
- I.M. Sechenov First Moscow State Medical University, Moscow, 119991, Russia
| | - Madina Baranova
- N.N. Blokhin Russian Cancer Research Center, Moscow, 115478, Russia
- Clinical Center Vitamed, Moscow, 121309, Russia
| | - Alex Glusker
- I.M. Sechenov First Moscow State Medical University, Moscow, 119991, Russia
| | - Lali Kogoniya
- M.F. Vladimirsky Moscow Regional Research Clinical Institute, Moscow, 129110, Russia
| | - Ekaterina Markarova
- M.F. Vladimirsky Moscow Regional Research Clinical Institute, Moscow, 129110, Russia
| | - Daria Allina
- I.M. Sechenov First Moscow State Medical University, Moscow, 119991, Russia
| | - Maria Suntsova
- Shemyakin-Ovchinnikov Institute of Bioorganic Chemistry, Moscow, 117997, Russia
| | | | | | - Marina Sekacheva
- I.M. Sechenov First Moscow State Medical University, Moscow, 119991, Russia
| | - Anton Buzdin
- I.M. Sechenov First Moscow State Medical University, Moscow, 119991, Russia
- Omicsway Corp., Walnut, California 91789, USA
- Shemyakin-Ovchinnikov Institute of Bioorganic Chemistry, Moscow, 117997, Russia
- Moscow Institute of Physics and Technology, Moscow Region, 141701, Russia
| |
Collapse
|
9
|
Moisseev A, Albert E, Lubarsky D, Schroeder D, Clark J. Transcriptomic and Genomic Testing to Guide Individualized Treatment in Chemoresistant Gastric Cancer Case. Biomedicines 2020; 8:biomedicines8030067. [PMID: 32210001 PMCID: PMC7148467 DOI: 10.3390/biomedicines8030067] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/25/2020] [Revised: 03/20/2020] [Accepted: 03/20/2020] [Indexed: 12/20/2022] Open
Abstract
Gastric cancer is globally the fifth leading cause of cancer death. We present a case report describing the unique genomic characteristics of an Epstein–Barr virus-negative gastric cancer with esophageal invasion and regional lymph node metastasis. Genomic tests were performed first with the stomach biopsy using platforms FoundationOne, OncoDNA, and Oncopanel at Dana Farber Institute. Following neoadjuvant chemotherapy, residual tumor was resected and the stomach and esophageal residual tumor samples were compared with the initial biopsy by whole exome sequencing and molecular pathway analysis platform Oncobox. Copy number variation profiling perfectly matched the whole exome sequencing results. A moderate agreement was seen between the diagnostic platforms in finding mutations in the initial biopsy. Final data indicate somatic activating mutation Q546K in PIK3CA gene, somatic frameshifts in PIH1D1 and FBXW7 genes, stop-gain in TP53BP1, and a few somatic mutations of unknown significance. RNA sequencing analysis revealed upregulated expressions of MMP7, MMP9, BIRC5, and PD-L1 genes and strongly differential regulation of several molecular pathways linked with the mutations identified. According to test results, the patient received immunotherapy with anti-PD1 therapy and is now free of disease for 2 years. Our data suggest that matched tumor and normal tissue analyses have a considerable advantage over tumor biopsy-only genomic tests in stomach cancer.
Collapse
Affiliation(s)
- Alexey Moisseev
- Institute for personalized medicine, I.M. Sechenov First Moscow State Medical University, 119048 Moscow, Russia;
- Correspondence: ; Tel.: +7(926)1443639
| | - Eugene Albert
- Institute for personalized medicine, I.M. Sechenov First Moscow State Medical University, 119048 Moscow, Russia;
| | | | - David Schroeder
- Wellesley Internal Medicine, 372 Washington St Ste 2, Wellesley Hills, MA 02481, USA;
| | - Jeffrey Clark
- Department of Hematology and Oncology, Massachusetts General Hospital, 55 Fruit Street Boston, MA 02114, USA;
| |
Collapse
|
10
|
Buzdin A, Sorokin M, Garazha A, Glusker A, Aleshin A, Poddubskaya E, Sekacheva M, Kim E, Gaifullin N, Giese A, Seryakov A, Rumiantsev P, Moshkovskii S, Moiseev A. RNA sequencing for research and diagnostics in clinical oncology. Semin Cancer Biol 2019; 60:311-323. [PMID: 31412295 DOI: 10.1016/j.semcancer.2019.07.010] [Citation(s) in RCA: 32] [Impact Index Per Article: 6.4] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/30/2019] [Accepted: 07/16/2019] [Indexed: 12/26/2022]
Abstract
Molecular diagnostics is becoming one of the major drivers of personalized oncology. With hundreds of different approved anticancer drugs and regimens of their administration, selecting the proper treatment for a patient is at least nontrivial task. This is especially sound for the cases of recurrent and metastatic cancers where the standard lines of therapy failed. Recent trials demonstrated that mutation assays have a strong limitation in personalized selection of therapeutics, consequently, most of the drugs cannot be ranked and only a small percentage of patients can benefit from the screening. Other approaches are, therefore, needed to address a problem of finding proper targeted therapies. The analysis of RNA expression (transcriptomic) profiles presents a reasonable solution because transcriptomics stands a few steps closer to tumor phenotype than the genome analysis. Several recent studies pioneered using transcriptomics for practical oncology and showed truly encouraging clinical results. The possibility of directly measuring of expression levels of molecular drugs' targets and profiling activation of the relevant molecular pathways enables personalized prioritizing for all types of molecular-targeted therapies. RNA sequencing is the most robust tool for the high throughput quantitative transcriptomics. Its use, potentials, and limitations for the clinical oncology will be reviewed here along with the technical aspects such as optimal types of biosamples, RNA sequencing profile normalization, quality controls and several levels of data analysis.
Collapse
Affiliation(s)
- Anton Buzdin
- I.M. Sechenov First Moscow State Medical University, Moscow, Russia; Omicsway Corp., Walnut, CA, USA; Shemyakin-Ovchinnikov Institute of Bioorganic Chemistry, Moscow, Russia.
| | - Maxim Sorokin
- I.M. Sechenov First Moscow State Medical University, Moscow, Russia; Omicsway Corp., Walnut, CA, USA; Shemyakin-Ovchinnikov Institute of Bioorganic Chemistry, Moscow, Russia
| | | | | | - Alex Aleshin
- Stanford University School of Medicine, Stanford, 94305, CA, USA
| | - Elena Poddubskaya
- I.M. Sechenov First Moscow State Medical University, Moscow, Russia; Vitamed Oncological Clinics, Moscow, Russia
| | - Marina Sekacheva
- I.M. Sechenov First Moscow State Medical University, Moscow, Russia
| | - Ella Kim
- Johannes Gutenberg University Mainz, Mainz, Germany
| | - Nurshat Gaifullin
- Lomonosov Moscow State University, Faculty of Medicine, Moscow, Russia
| | | | | | | | - Sergey Moshkovskii
- Institute of Biomedical Chemistry, Moscow, 119121, Russia; Pirogov Russian National Research Medical University (RNRMU), Moscow, 117997, Russia
| | - Alexey Moiseev
- I.M. Sechenov First Moscow State Medical University, Moscow, Russia
| |
Collapse
|