Bioinformatics for precision oncology

Utility analyses of AVITI sequencing chemistry

BACKGROUND

DNA sequencing is a critical tool in modern biology. Over the last two decades, it has been revolutionized by the advent of massively parallel sequencing, leading to significant advances in the genome and transcriptome sequencing of various organisms. Nevertheless, challenges with accuracy, lack of competitive options and prohibitive costs associated with high throughput parallel short-read sequencing persist.

RESULTS

Here, we conduct a comparative analysis using matched DNA and RNA short-reads assays between Element Biosciences' AVITI and Illumina's NextSeq 550 chemistries. Similar comparisons were evaluated for synthetic long-read sequencing for RNA and targeted single-cell transcripts between the AVITI and Illumina's NovaSeq 6000. For both DNA and RNA short-read applications, the study found that the AVITI produced significantly higher per sequence quality scores. For PCR-free DNA libraries, we observed an average 89.7% lower experimentally determined error rate when using the AVITI chemistry, compared to the NextSeq 550. For short-read RNA quantification, AVITI platform had an average of 32.5% lower error rate than that for NextSeq 550. With regards to synthetic long-read mRNA and targeted synthetic long read single cell mRNA sequencing, both platforms' respective chemistries performed comparably in quantification of genes and isoforms. The AVITI displayed a marginally lower error rate for long reads, with fewer chemistry-specific errors and a higher mutation detection rate.

CONCLUSION

These results point to the potential of the AVITI platform as a competitive candidate in high-throughput short read sequencing analyses when juxtaposed with the Illumina NextSeq 550.

Target-Driven Tissue-Agnostic Drug Approvals-A New Path of Drug Development

The regulatory approvals of tumor-agnostic therapies have led to the re-evaluation of the drug development process. The conventional models of drug development are histology-based. On the other hand, the tumor-agnostic drug development of a new drug (or combination) focuses on targeting a common genomic biomarker in multiple cancers, regardless of histology. The basket-like clinical trials with multiple cohorts allow clinicians to evaluate pan-cancer efficacy and toxicity. There are currently eight tumor agnostic approvals granted by the Food and Drug Administration (FDA). This includes two immune checkpoint inhibitors, and five targeted therapy agents. Pembrolizumab is an anti-programmed cell death protein-1 (PD-1) antibody that was the first FDA-approved tumor-agnostic treatment for unresectable or metastatic microsatellite instability-high (MSI-H) or deficient mismatch repair (dMMR) solid tumors in 2017. It was later approved for tumor mutational burden-high (TMB-H) solid tumors, although the TMB cut-off used is still debated. Subsequently, in 2021, another anti-PD-1 antibody, dostarlimab, was also approved for dMMR solid tumors in the refractory setting. Patients with fusion-positive cancers are typically difficult to treat due to their rare prevalence and distribution. Gene rearrangements or fusions are present in a variety of tumors. Neurotrophic tyrosine kinase (NTRK) fusions are present in a range of pediatric and adult solid tumors in varying frequency. Larotrectinib and entrectinib were approved for neurotrophic tyrosine kinase (NTRK) fusion-positive cancers. Similarly, selpercatinib was approved for rearranged during transfection (RET) fusion-positive solid tumors. The FDA approved the first combination therapy of dabrafenib, a B-Raf proto-oncogene serine/threonine kinase (BRAF) inhibitor, plus trametinib, a mitogen-activated protein kinase (MEK) inhibitor for patients 6 months or older with unresectable or metastatic tumors (except colorectal cancer) carrying a BRAFV600E mutation. The most recent FDA tumor-agnostic approval is of fam-trastuzumab deruxtecan-nxki (T-Dxd) for HER2-positive solid tumors. It is important to identify and expeditiously develop drugs that have the potential to provide clinical benefit across tumor types.

Utility Analyses of AVITI Sequencing Chemistry

Background

DNA sequencing is a critical tool in modern biology. Over the last two decades, it has been revolutionized by the advent of massively parallel sequencing, leading to significant advances in the genome and transcriptome sequencing of various organisms. Nevertheless, challenges with accuracy, lack of competitive options and prohibitive costs associated with high throughput parallel short-read sequencing persist.

Results

Here, we conduct a comparative analysis using matched DNA and RNA short-reads assays between Element Biosciences' AVITI and Illumina's NextSeq 550 chemistries. Similar comparisons were evaluated for synthetic long-read sequencing for RNA and targeted single-cell transcripts between the AVITI and Illumina's NovaSeq 6000. For both DNA and RNA short-read applications, the study found that the AVITI produced significantly higher per sequence quality scores. For PCR-free DNA libraries, we observed an average 89.7% lower experimentally determined error rate when using the AVITI chemistry, compared to the NextSeq 550. For short-read RNA quantification, AVITI platform had an average of 32.5% lower error rate than that for NextSeq 550. With regards to synthetic long-read mRNA and targeted synthetic long read single cell mRNA sequencing, both platforms' respective chemistries performed comparably in quantification of genes and isoforms. The AVITI displayed a marginally lower error rate for long reads, with fewer chemistry-specific errors and a higher mutation detection rate.

Conclusion

These results point to the potential of the AVITI platform as a competitive candidate in high-throughput short read sequencing analyses when juxtaposed with the Illumina NextSeq 550.

Evaluation of single-sample network inference methods for precision oncology

A major challenge in precision oncology is to detect targetable cancer vulnerabilities in individual patients. Modeling high-throughput omics data in biological networks allows identifying key molecules and processes of tumorigenesis. Traditionally, network inference methods rely on many samples to contain sufficient information for learning, resulting in aggregate networks. However, to implement patient-tailored approaches in precision oncology, we need to interpret omics data at the level of individual patients. Several single-sample network inference methods have been developed that infer biological networks for an individual sample from bulk RNA-seq data. However, only a limited comparison of these methods has been made and many methods rely on 'normal tissue' samples as reference, which are not always available. Here, we conducted an evaluation of the single-sample network inference methods SSN, LIONESS, SWEET, iENA, CSN and SSPGI using transcriptomic profiles of lung and brain cancer cell lines from the CCLE database. The methods constructed functional gene networks with distinct network characteristics. Hub gene analyses revealed different degrees of subtype-specificity across methods. Single-sample networks were able to distinguish between tumor subtypes, as exemplified by node strength clustering, enrichment of known subtype-specific driver genes among hubs and differential node strength. We also showed that single-sample networks correlated better to other omics data from the same cell line as compared to aggregate networks. We conclude that single-sample network inference methods can reflect sample-specific biology when 'normal tissue' samples are absent and we point out peculiarities of each method.

Advances in next-generation sequencing and emerging technologies for hematologic malignancies

Innovations in molecular diagnostics have often evolved through the study of hematologic malignancies. Examples include the pioneering characterization of the Philadelphia chromosome by cytogenetics in the 1970s, the implementation of polymerase chain reaction for high-sensitivity detection and monitoring of mutations and, most recently, targeted next- generation sequencing to drive the prognostic and therapeutic assessment of leukemia. Hematologists and hematopath- ologists have continued to advance in the past decade with new innovations improving the type, amount, and quality of data generated for each molecule of nucleic acid. In this review article, we touch on these new developments and discuss their implications for diagnostics in hematopoietic malignancies. We review advances in sequencing platforms and library preparation chemistry that can lead to faster turnaround times, novel sequencing techniques, the development of mobile laboratories with implications for worldwide benefits, the current status of sample types, improvements to quality and reference materials, bioinformatic pipelines, and the integration of machine learning and artificial intelligence into mol- ecular diagnostic tools for hematologic malignancies.

Comprehensive analysis of cuproptosis and copper homeostasis genotyping and related immune land scape in lung adenocarcinoma

Cuproptosis is a manner of cell death which is related to the homeostasis of copper ions in the cellular environment and is expected to open a new direction of anti-tumor therapy. However, the studies on cuproptosis and copper homeostasis in lung adenocarcinoma (LUAD) are still limited. In this study, we identified new cuproptosis and copper homeostasis related genes (CHRGs) which were effective in stratifying genotyping clusters with survival differences based on transcriptomic data obtained from The Cancer Genome Atlas (TCGA) and Gene Expression Omnibus (GEO). Weighted Gene Co-expression Network Analysis (WGCNA) further expands the screening boundary of CHRGs, and finally we established a 10-CHRGs-based prognostic signature using lasso-penalized cox regression method, which were validated in GSE30219. Comprehensive bioinformatics analysis revealed these genes are potential regulators of modulating immunotherapy efficacy, drug resistance, tumor microenvironment infiltration, and tumor mutation patterns. Lastly, the scRNA-seq datasets GSE183219 and GSE203360 offers the evidences that CHRGs signature are mainly distributed in cancer epithelial cells, real time quantitative polymerase chain reaction (RT-qPCR) also confirmed the differential expression of these genes between normal lung cell line and lung adenocarcinoma cell lines. Collectively, our findings revealed new cuproptosis and copper homeostasis related genotyping clusters and genes which may play important roles in predicting prognosis, influencing tumor microenvironment and drug efficacy in LUAD patients.

Bioinformatics in Early Cancer Detection

Bioinformatics is a pretty recent branch of biology that encompasses the use of algebraic, analytic, and computing approaches to the processing and interpretation of biological information. A wide term, "bioinformatics" refers to the use of digital technology to study biological processes using high-dimensional data collected from many resources. The design and testing of the software tools required to evaluate the information are the core of bioinformatics research, which is conducted in great portions in silico and typically involves the synthesis of new learning from available data. Early diagnosis of cancer results in improved prognosis, but at the same time, it is difficult to conform to diagnosis at a very early stage. The use of DNA microarrays and proteomics studies for large-scale gene expression research has advanced technology, thus elevating the significance of bioinformatics tools. In today's research, wet experimentation and the application of bioinformatics analytics go side by side. Molecular profiling of tumor biopsies is becoming more and more crucial to both cancer research and the treatment of cancer.

Exploring the Role of Enhancer-Mediated Transcriptional Regulation in Precision Biology

The emergence of precision biology has been driven by the development of advanced technologies and techniques in high-resolution biological research systems. Enhancer-mediated transcriptional regulation, a complex network of gene expression and regulation in eukaryotes, has attracted significant attention as a promising avenue for investigating the underlying mechanisms of biological processes and diseases. To address biological problems with precision, large amounts of data, functional information, and research on the mechanisms of action of biological molecules is required to address biological problems with precision. Enhancers, including typical enhancers and super enhancers, play a crucial role in gene expression and regulation within this network. The identification and targeting of disease-associated enhancers hold the potential to advance precision medicine. In this review, we present the concepts, progress, importance, and challenges in precision biology, transcription regulation, and enhancers. Furthermore, we propose a model of transcriptional regulation for multi-enhancers and provide examples of their mechanisms in mammalian cells, thereby enhancing our understanding of how enhancers achieve precise regulation of gene expression in life processes. Precision biology holds promise in providing new tools and platforms for discovering insights into gene expression and disease occurrence, ultimately benefiting individuals and society as a whole.

An immune-related prognostic model predicts neoplasm-immunity interactions for metastatic nasopharyngeal carcinoma

Background

The prognosis of nasopharyngeal carcinoma (NPC) has been recognized to improve immensely owing to radiotherapy combined with chemotherapy. However, patients with metastatic NPC have a poor prognosis. Immunotherapy has dramatically prolonged the survival of patients with NPC. Hence, further research on immune-related biomarkers is imperative to establish the prognosis of metastatic NPC.

Methods

10 NPC RNA expression profiles were generated from patients with or without distant metastasis after chemoradiotherapy from the Fujian Cancer Hospital. The differential immune-related genes were identified and validated by immunohistochemistry analysis. The method of least absolute shrinkage and selection operator (LASSO)was used to further establish the immune-related prognostic model in an external GEO database (GSE102349, n=88). The immune microenvironment and signal pathways were evaluated in multiple dimensions at the transcriptome and single-cell levels.

Results

1328 differential genes were identified, out of which 520 were upregulated and 808 were downregulated. Notably, most of the immune genes and pathways were down-regulated in the metastasis group. A prognostic immune model involving nine hub genes. Patients in low-risk group were characterized by survival advantage, hot immune phenotype and benefit from immunotherapy. Compared with immune cells, malignant cell exhibited the most active levels of risk score by ssGSEA. Accordingly, intercellular communications including LT, CD70, CD40 and SPP1, and the like, between high-risk and low-risk were explored by the R package “Cellchat”.

Conclusion

We have constructed a model based on immunity of metastatic NPC and determined its prognostic value. The model identified the level of immune cell infiltration, cell-cell communication, along with potential immunotherapy for metastatic NPC.

Comprehensive pan-cancer analysis identifies centromere associated protein E as a novel prognostic and immunological biomarker in human tumors

Centromere-associated protein E (CENP-E), a core component of the kinetochore, mediates chromosome congression and spindle microtubule capture during mitosis. Partial experimental evidence has illustrated the carcinogenic effects of CENPE in tumors, but the corresponding pan-cancer analysis of CENPE still lacking. Based on public databases, including the Cancer Genome Atlas (TCGA), Genotype-Tissue Expression (GTEx), and Human Protein Atlas (HPA), we take an array of bioinformatics methods to investigate the potential oncogenic roles of CENPE. Then, we validated CENPE, cell cycle-related proteins, and immune checkpoint molecule findings expression in clinical colon cancer samples by western blot. Our results showed that CENPE was up-regulated in almost all tumors, and the expression level of CENPE was associated with worse overall survival (OS) and disease-specific survival (DSS) in patients. The strong relationship between CENPE with gene mutation and MMR has also been validated. Moreover, CENPE gene expression was positively correlated with immune checkpoint molecular, and reversely correlated with infiltration levels of most immune cells. In the human colon cancer tissues, the expression of CENPE, cell cycle-related proteins, and immune checkpoint molecules were significantly higher than in the adjacent normal tissues. Our results indicated that CENPE can function as an oncogene in various cancers, and may be regarded as a promising prognostic and diagnostic biomarker in cancer treatment.

Data-driven support to decision-making in molecular tumour boards for lymphoma: A design science approach

Background

The increasing amount of molecular data and knowledge about genomic alterations from next-generation sequencing processes together allow for a greater understanding of individual patients, thereby advancing precision medicine. Molecular tumour boards feature multidisciplinary teams of clinical experts who meet to discuss complex individual cancer cases. Preparing the meetings is a manual and time-consuming process.

Purpose

To design a clinical decision support system to improve the multimodal data interpretation in molecular tumour board meetings for lymphoma patients at Karolinska University Hospital, Stockholm, Sweden. We investigated user needs and system requirements, explored the employment of artificial intelligence, and evaluated the proposed design with primary stakeholders.

Methods

Design science methodology was used to form and evaluate the proposed artefact. Requirements elicitation was done through a scoping review followed by five semi-structured interviews. We used UML Use Case diagrams to model user interaction and UML Activity diagrams to inform the proposed flow of control in the system. Additionally, we modelled the current and future workflow for MTB meetings and its proposed machine learning pipeline. Interactive sessions with end-users validated the initial requirements based on a fictive patient scenario which helped further refine the system.

Results

The analysis showed that an interactive secure Web-based information system supporting the preparation of the meeting, multidisciplinary discussions, and clinical decision-making could address the identified requirements. Integrating artificial intelligence via continual learning and multimodal data fusion were identified as crucial elements that could provide accurate diagnosis and treatment recommendations.

Impact

Our work is of methodological importance in that using artificial intelligence for molecular tumour boards is novel. We provide a consolidated proof-of-concept system that could support the end-to-end clinical decision-making process and positively and immediately impact patients.

Conclusion

Augmenting a digital decision support system for molecular tumour boards with retrospective patient material is promising. This generates realistic and constructive material for human learning, and also digital data for continual learning by data-driven artificial intelligence approaches. The latter makes the future system adaptable to human bias, improving adequacy and decision quality over time and over tasks, while building and maintaining a digital log.

Spatial Transcriptomic Analysis Reveals Associations between Genes and Cellular Topology in Breast and Prostate Cancers

BACKGROUND

Cancer is the leading cause of death worldwide with breast and prostate cancer the most common among women and men, respectively. Gene expression and image features are independently prognostic of patient survival; but until the advent of spatial transcriptomics (ST), it was not possible to determine how gene expression of cells was tied to their spatial relationships (i.e., topology).

METHODS

We identify topology-associated genes (TAGs) that correlate with 700 image topological features (ITFs) in breast and prostate cancer ST samples. Genes and image topological features are independently clustered and correlated with each other. Themes among genes correlated with ITFs are investigated by functional enrichment analysis.

RESULTS

Overall, topology-associated genes (TAG) corresponding to extracellular matrix (ECM) and Collagen Type I Trimer gene ontology terms are common to both prostate and breast cancer. In breast cancer specifically, we identify the ZAG-PIP Complex as a TAG. In prostate cancer, we identify distinct TAGs that are enriched for GI dysmotility and the IgA immunoglobulin complex. We identified TAGs in every ST slide regardless of cancer type.

CONCLUSIONS

These TAGs are enriched for ontology terms, illustrating the biological relevance to our image topology features and their potential utility in diagnostic and prognostic models.

Bioinformatics: From NGS Data to Biological Complexity in Variant Detection and Oncological Clinical Practice

The use of next-generation sequencing (NGS) techniques for variant detection has become increasingly important in clinical research and in clinical practice in oncology. Many cancer patients are currently being treated in clinical practice or in clinical trials with drugs directed against specific genomic alterations. In this scenario, the development of reliable and reproducible bioinformatics tools is essential to derive information on the molecular characteristics of each patient’s tumor from the NGS data. The development of bioinformatics pipelines based on the use of machine learning and statistical methods is even more relevant for the determination of complex biomarkers. In this review, we describe some important technologies, computational algorithms and models that can be applied to NGS data from Whole Genome to Targeted Sequencing, to address the problem of finding complex cancer-associated biomarkers. In addition, we explore the future perspectives and challenges faced by bioinformatics for precision medicine both at a molecular and clinical level, with a focus on an emerging complex biomarker such as homologous recombination deficiency (HRD).

J-SPACE: a Julia package for the simulation of spatial models of cancer evolution and of sequencing experiments

Background

The combined effects of biological variability and measurement-related errors on cancer sequencing data remain largely unexplored. However, the spatio-temporal simulation of multi-cellular systems provides a powerful instrument to address this issue. In particular, efficient algorithmic frameworks are needed to overcome the harsh trade-off between scalability and expressivity, so to allow one to simulate both realistic cancer evolution scenarios and the related sequencing experiments, which can then be used to benchmark downstream bioinformatics methods.

Result

We introduce a Julia package for SPAtial Cancer Evolution (J-SPACE), which allows one to model and simulate a broad set of experimental scenarios, phenomenological rules and sequencing settings.Specifically, J-SPACE simulates the spatial dynamics of cells as a continuous-time multi-type birth-death stochastic process on a arbitrary graph, employing different rules of interaction and an optimised Gillespie algorithm. The evolutionary dynamics of genomic alterations (single-nucleotide variants and indels) is simulated either under the Infinite Sites Assumption or several different substitution models, including one based on mutational signatures. After mimicking the spatial sampling of tumour cells, J-SPACE returns the related phylogenetic model, and allows one to generate synthetic reads from several Next-Generation Sequencing (NGS) platforms, via the ART read simulator. The results are finally returned in standard FASTA, FASTQ, SAM, ALN and Newick file formats.

Conclusion

J-SPACE is designed to efficiently simulate the heterogeneous behaviour of a large number of cancer cells and produces a rich set of outputs. Our framework is useful to investigate the emergent spatial dynamics of cancer subpopulations, as well as to assess the impact of incomplete sampling and of experiment-specific errors. Importantly, the output of J-SPACE is designed to allow the performance assessment of downstream bioinformatics pipelines processing NGS data. J-SPACE is freely available at: https://github.com/BIMIB-DISCo/J-Space.jl.

Application of non-negative matrix factorization in oncology: one approach for establishing precision medicine

The increase in the expectations of artificial intelligence (AI) technology has led to machine learning technology being actively used in the medical field. Non-negative matrix factorization (NMF) is a machine learning technique used for image analysis, speech recognition, and language processing; recently, it is being applied to medical research. Precision medicine, wherein important information is extracted from large-scale medical data to provide optimal medical care for every individual, is considered important in medical policies globally, and the application of machine learning techniques to this end is being handled in several ways. NMF is also introduced differently because of the characteristics of its algorithms. In this review, the importance of NMF in the field of medicine, with a focus on the field of oncology, is described by explaining the mathematical science of NMF and the characteristics of the algorithm, providing examples of how NMF can be used to establish precision medicine, and presenting the challenges of NMF. Finally, the direction regarding the effective use of NMF in the field of oncology is also discussed.

Sequential chemotherapy and icotinib as first-line treatment for advanced epidermal growth factor receptor-mutated non-small cell lung cancer

BACKGROUND

Icotinib could have potential effect and tolerability when used sequentially with chemotherapy for advanced epidermal growth factor receptor (EGFR)-mutated non-small cell lung cancer (NSCLC).

AIM

To evaluate the efficacy and safety of chemotherapy followed by icotinib maintenance therapy as first-line treatment for advanced EGFR-mutated NSCLC.

METHODS

This multicenter, open-label, pilot randomized controlled trial enrolled 68 EGFR-mutated stage IIIB/IV NSCLC patients randomized 2:3 to the icotinib alone and chemotherapy + icotinib groups.

RESULTS

The median progression-free survival in the icotinib alone and chemotherapy + icotinib groups was 8.0 mo (95%CI: 3.84-11.63) and 13.4 mo (95%CI: 10.18-16.33), respectively (P = 0.0249). No significant differences were found in the curative effect when considering different cycles of chemotherapy or chemotherapy regimen (all P > 0.05).

CONCLUSION

A sequential combination of chemotherapy and EGFR-tyrosine kinase inhibitor is feasible for stage IV EGFR-mutated NSCLC patients.

5-Methylcytosine RNA Methyltransferases-Related Long Non-coding RNA to Develop and Validate Biochemical Recurrence Signature in Prostate Cancer

The effects of 5-methylcytosine in RNA (m5C) in various human cancers have been increasingly studied recently; however, the m5C regulator signature in prostate cancer (PCa) has not been well established yet. In this study, we identified and characterized a series of m5C-related long non-coding RNAs (lncRNAs) in PCa. Univariate Cox regression analysis and least absolute shrinkage and selector operation (LASSO) regression analysis were implemented to construct a m5C-related lncRNA prognostic signature. Consequently, a prognostic m5C-lnc model was established, including 17 lncRNAs: MAFG-AS1, AC012510.1, AC012065.3, AL117332.1, AC132192.2, AP001160.2, AC129510.1, AC084018.2, UBXN10-AS1, AC138956.2, ZNF32-AS2, AC017100.1, AC004943.2, SP2-AS1, Z93930.2, AP001486.2, and LINC01135. The high m5C-lnc score calculated by the model significantly relates to poor biochemical recurrence (BCR)-free survival (p < 0.0001). Receiver operating characteristic (ROC) curves and a decision curve analysis (DCA) further validated the accuracy of the prognostic model. Subsequently, a predictive nomogram combining the prognostic model with clinical features was created, and it exhibited promising predictive efficacy for BCR risk stratification. Next, the competing endogenous RNA (ceRNA) network and lncRNA–protein interaction network were established to explore the potential functions of these 17 lncRNAs mechanically. In addition, functional enrichment analysis revealed that these lncRNAs are involved in many cellular metabolic pathways. Lastly, MAFG-AS1 was selected for experimental validation; it was upregulated in PCa and probably promoted PCa proliferation and invasion in vitro. These results offer some insights into the m5C's effects on PCa and reveal a predictive model with the potential clinical value to improve the prognosis of patients with PCa.

Knowledge bases and software support for variant interpretation in precision oncology

Precision oncology is a rapidly evolving interdisciplinary medical specialty. Comprehensive cancer panels are becoming increasingly available at pathology departments worldwide, creating the urgent need for scalable cancer variant annotation and molecularly informed treatment recommendations. A wealth of mainly academia-driven knowledge bases calls for software tools supporting the multi-step diagnostic process. We derive a comprehensive list of knowledge bases relevant for variant interpretation by a review of existing literature followed by a survey among medical experts from university hospitals in Germany. In addition, we review cancer variant interpretation tools, which integrate multiple knowledge bases. We categorize the knowledge bases along the diagnostic process in precision oncology and analyze programmatic access options as well as the integration of knowledge bases into software tools. The most commonly used knowledge bases provide good programmatic access options and have been integrated into a range of software tools. For the wider set of knowledge bases, access options vary across different parts of the diagnostic process. Programmatic access is limited for information regarding clinical classifications of variants and for therapy recommendations. The main issue for databases used for biological classification of pathogenic variants and pathway context information is the lack of standardized interfaces. There is no single cancer variant interpretation tool that integrates all identified knowledge bases. Specialized tools are available and need to be further developed for different steps in the diagnostic process.

Statistical Process Control Charts for Monitoring Next-Generation Sequencing and Bioinformatics Turnaround in Precision Medicine Initiatives

Purpose

Precision oncology, such as next generation sequencing (NGS) molecular analysis and bioinformatics are used to guide targeted therapies. The laboratory turnaround time (TAT) is a key performance indicator of laboratory performance. This study aims to formally apply statistical process control (SPC) methods such as CUSUM and EWMA to a precision medicine programme to analyze the learning curves of NGS and bioinformatics processes.

Patients and Methods

Trends in NGS and bioinformatics TAT were analyzed using simple regression models with TAT as the dependent variable and chronologically-ordered case number as the independent variable. The M-estimator "robust" regression and negative binomial regression were chosen to serve as sensitivity analyses to each other. Next, two popular statistical process control (SPC) approaches which are CUSUM and EWMA were utilized and the CUSUM log-likelihood ratio (LLR) charts were also generated. All statistical analyses were done in Stata version 16.0 (StataCorp), and nominal P < 0.05 was considered to be statistically significant.

Results

A total of 365 patients underwent successful molecular profiling. Both the robust linear model and negative binomial model showed statistically significant reductions in TAT with accumulating experience. The EWMA and CUSUM charts of overall TAT largely corresponded except that the EWMA chart consistently decreased while the CUSUM analyses indicated improvement only after a nadir at the 82^nd case. CUSUM analysis found that the bioinformatics team took a lower number of cases (54 cases) to overcome the learning curve compared to the NGS team (85 cases).

Conclusion

As NGS and bioinformatics lead precision oncology into the forefront of cancer management, characterizing the TAT of NGS and bioinformatics processes improves the timeliness of data output by potentially spotlighting problems early for rectification, thereby improving care delivery.

A gentle introduction to understanding preclinical data for cancer pharmaco-omic modeling

A central goal of precision oncology is to administer an optimal drug treatment to each cancer patient. A common preclinical approach to tackle this problem has been to characterize the tumors of patients at the molecular and drug response levels, and employ the resulting datasets for predictive in silico modeling (mostly using machine learning). Understanding how and why the different variants of these datasets are generated is an important component of this process. This review focuses on providing such introduction aimed at scientists with little previous exposure to this research area.

Machine learning in the prediction of cancer therapy

Resistance to therapy remains a major cause of cancer treatment failures, resulting in many cancer-related deaths. Resistance can occur at any time during the treatment, even at the beginning. The current treatment plan is dependent mainly on cancer subtypes and the presence of genetic mutations. Evidently, the presence of a genetic mutation does not always predict the therapeutic response and can vary for different cancer subtypes. Therefore, there is an unmet need for predictive models to match a cancer patient with a specific drug or drug combination. Recent advancements in predictive models using artificial intelligence have shown great promise in preclinical settings. However, despite massive improvements in computational power, building clinically useable models remains challenging due to a lack of clinically meaningful pharmacogenomic data. In this review, we provide an overview of recent advancements in therapeutic response prediction using machine learning, which is the most widely used branch of artificial intelligence. We describe the basics of machine learning algorithms, illustrate their use, and highlight the current challenges in therapy response prediction for clinical practice.

[Molecular tumor boards - insights and perspectives]

The rapid development of molecular technologies and targeted therapies has fostered the implementation of specialized tumor conferences, known as molecular tumor boards (MTBs). MTBs become particularly important when treatment recommendations are needed based on molecular alterations beyond the approved targeted therapies. While an MTB's goals are based on individualized diagnostics and therapies of tumor patients using innovative technologies and biomarkers, the procedures of MTBs are still quite heterogeneous. This applies to the primary inclusion criteria for tumor patients, the composition of MTBs, the applied diagnostic tests and their assessment and reporting, the evaluation of their clinical value and implementation in a therapeutic strategy, and the associated quality assurance measurements as well as knowledge-gaining, economical, legal, and ethical aspects.This article provides an overview of the spectrum of MTBs, their challenges, and the potential for individualized cancer medicine.

Conserved pan-cancer microenvironment subtypes predict response to immunotherapy

The clinical use of molecular targeted therapy is rapidly evolving but has primarily focused on genomic alterations. Transcriptomic analysis offers an opportunity to dissect the complexity of tumors, including the tumor microenvironment (TME), a crucial mediator of cancer progression and therapeutic outcome. TME classification by transcriptomic analysis of >10,000 cancer patients identifies four distinct TME subtypes conserved across 20 different cancers. The TME subtypes correlate with patient response to immunotherapy in multiple cancers, with patients possessing immune-favorable TME subtypes benefiting the most from immunotherapy. Thus, the TME subtypes act as a generalized immunotherapy biomarker across many cancer types due to the inclusion of malignant and microenvironment components. A visual tool integrating transcriptomic and genomic data provides a global tumor portrait, describing the tumor framework, mutational load, immune composition, anti-tumor immunity, and immunosuppressive escape mechanisms. Integrative analyses plus visualization may aid in biomarker discovery and the personalization of therapeutic regimens.

Comprehensive tumor molecular profile analysis in clinical practice

BACKGROUND

Tumor molecular profile analysis by Next Generation Sequencing technology is currently widely applied in clinical practice and has enabled the detection of predictive biomarkers of response to targeted treatment. In parallel with targeted therapies, immunotherapies are also evolving, revolutionizing cancer therapy, with Programmed Death-ligand 1 (PD-L1), Microsatellite instability (MSI), and Tumor Mutational Burden (TMB) analysis being the biomarkers employed most commonly.

METHODS

In the present study, tumor molecular profile analysis was performed using a 161 gene NGS panel, containing the majority of clinically significant genes for cancer treatment selection. A variety of tumor types have been analyzed, including aggressive and hard to treat cancers such as pancreatic cancer. Besides, the clinical utility of immunotherapy biomarkers (TMB, MSI, PD-L1), was also studied.

RESULTS

Molecular profile analysis was conducted in 610 cancer patients, while in 393 of them a at least one biomarker for immunotherapy response was requested. An actionable alteration was detected in 77.87% of the patients. 54.75% of them received information related to on-label or off-label treatment (Tiers 1A.1, 1A.2, 2B, and 2C.1) and 21.31% received a variant that could be used for clinical trial inclusion. The addition to immunotherapy biomarker to targeted biomarkers' analysis in 191 cases increased the number of patients with an on-label treatment recommendation by 22.92%, while an option for on-label or off-label treatment was provided in 71.35% of the cases.

CONCLUSIONS

Tumor molecular profile analysis using NGS is a first-tier method for a variety of tumor types and provides important information for decision making in the treatment of cancer patients. Importantly, simultaneous analysis for targeted therapy and immunotherapy biomarkers could lead to better tumor characterization and offer actionable information in the majority of patients. Furthermore, our data suggest that one in two patients may be eligible for on-label ICI treatment based on biomarker analysis. However, appropriate interpretation of results from such analysis is essential for implementation in clinical practice and accurate refinement of treatment strategy.

Identification of core genes and clinical outcomes in tumors originated from endoderm (gastric cancer and lung carcinoma) via bioinformatics analysis

ABSTRACT

During last decade, bioinformatics analysis has provided an effective way to study the relationship between various genes and biological processes. In this study, we aimed to identify potential core candidate genes and underlying mechanisms of progression of lung and gastric carcinomas which both originated from endoderm. The expression profiles, GSE54129 (gastric carcinoma) and GSE27262 (lung carcinoma), were collected from GEO database. One hundred eleven patients with gastric carcinoma and 21 health people were included in this research. Meanwhile, there were 25 lung carcinoma patients. Then, 75 differentially expressed genes were selected via GEO2R online tool and Venn software, including 31 up-regulated genes and 44 down-regulated genes. Next, we used Database for Annotation, Visualization, and Integrated Discovery and Metascpe software to analyze Kyoto Encyclopedia of Gene and Genome pathway and gene ontology. Furthermore, Cytoscape software and MCODE App were performed to construct complex of these differentially expressed genes . Twenty core genes were identified, which mainly enriched in extracellular matrix-receptor interaction, focal adhesion, and PI3K-Akt pathway (P < .01). Finally, the significant difference of gene expression between cancer tissues and normal tissues in both lung and gastric carcinomas was examined by Gene Expression Profiling Interactive Analysis database. Twelve candidate genes with positive statistical significance (P < .01), COMP CTHRC1 COL1A1 SPP1 COL11A1 COL10A1 CXCL13 CLDN3 CLDN1 matrix metalloproteinases 7 ADAM12 PLAU, were picked out to further analysis. The Kaplan-Meier plotter website was applied to examine relationship among these genes and clinical outcomes. We found 4 genes (ADAM12, SPP1, COL1A1, COL11A1) were significantly associated with poor prognosis in both lung and gastric carcinoma patients (P  < .05). In conclusion, these candidate genes may be potential therapeutic targets for cancer treatment.

Artificial Intelligence in Decision-Making for Colorectal Cancer Treatment Strategy: An Observational Study of Implementing Watson for Oncology in a 250-Case Cohort

Background

Personalized and novel evidence-based clinical treatment strategy consulting for colorectal cancer has been available through various artificial intelligence (AI) supporting systems such as Watson for Oncology (WFO) from IBM. However, the potential effects of this supporting tool in cancer care have not been thoroughly explored in real-world studies. This research aims to investigate the concordance between treatment recommendations for colorectal cancer patients made by WFO and a multidisciplinary team (MDT) at a major comprehensive gastrointestinal cancer center.

Methods

In this prospective study, both WFO and the blinded MDT's treatment recommendations were provided concurrently for enrolled colorectal cancers of stages II to IV between March 2017 and January 2018 at Shanghai Minimally Invasive Surgery Center. Concordance was achieved if the cancer team's decisions were listed in the "recommended" or "for consideration" classification in WFO. A review was carried out after 100 cases for all non-concordant patients to explain the inconsistency, and corresponding feedback was given to WFO's database. The concordance of the subsequent cases was analyzed to evaluate both the performance and learning ability of WFO.

Results

Overall, 250 patients met the inclusion criteria and were recruited in the study. Eighty-one were diagnosed with colon cancer and 189 with rectal cancer. The concordances for colon cancer, rectal cancer, or overall were all 91%. The overall rates were 83, 94, and 88% in subgroups of stages II, III, and IV. When categorized by treatment strategy, concordances were 97, 93, 89, 87, and 100% for neoadjuvant, surgery, adjuvant, first line, and second line treatment groups, respectively. After analyzing the main factors causing discordance, relative updates were made in the database accordingly, which led to the concordance curve rising in most groups compared with the initial rates.

Conclusion

Clinical recommendations made by WFO and the cancer team were highly matched for colorectal cancer. Patient age, cancer stage, and the consideration of previous therapy details had a significant influence on concordance. Addressing these perspectives will facilitate the use of the cancer decision-support systems to help oncologists achieve the promise of precision medicine.

Dihydroartemisinin inhibits the tumorigenesis and invasion of gastric cancer by regulating STAT1/KDR/MMP9 and P53/BCL2L1/CASP3/7 pathways

Dihydroartemisinin (DHA), an effective antimalarial drug, has been widely investigated as an anti-tumor agent. Although previous studies have indicated the potential therapeutic effects of DHA on multiple malignancies, its detailed molecular mechanisms in gastric cancer (GC) are still undocumented. In the present study, we applied network pharmacology and bioinformatics (gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses) to obtain the collective targets of DHA and GC and analyzed their involvement in constructing a protein-protein interaction (PPI) network. The top 10% hub targets in this network were identified, and TCGA database was utilized for the single gene analysis of their correlation with the prognosis of GC. CCK8, EdU, Transwell, and flow cytometry analyses were conducted, and subcutaneous xenograft tumor models were constructed to assess the effects of DHA on the tumorigenesis and invasion of GC. Furthermore, the targets of DHA were verified by molecular docking, quantitative real-time PCR (qPCR) and western blot analyses in GC cells. The results indicated that the common targets of DHA and GC were enriched in multiple cancer-related pathways including KDR, STAT1 and apoptosis signaling pathways, where the core genes included KDR, MMP9, STAT1, TP53, CASP3/7 and BCL2L1. The lowered expression of KDR and increased expression of TP53 and CASP7 harbored a favorable survival for patients with GC patients. CASP7 showed a positive correlation with CASP3 but a negative correlation with KDR and could be regarded as an independent protective factor for overall survival in GC. Moreover, DHA treatment induced cell apoptosis and suppressed the cell proliferation, DNA synthesis, cycle progression and invasive capabilities both in vitro and in vivo. DHA also upregulated p53, CASP3, and cleaved-CASP3 and downregulated BCL2L1, MMP9, KDR, p-KDR, STAT1 and p-STAT1 in GC cell lines. In conclusion, DHA could suppress the tumorigenesis and invasion of GC by regulating STAT1/KDR/MMP9 and p53/BCL2L1/CASP3/7 pathways. Our findings might provide a novel approach for the treatment of GC.

Bioinformatics workflows for clinical applications in precision oncology

High-throughput molecular profiling of tumors is a fundamental aspect of precision oncology, enabling the identification of genomic alterations that can be targeted therapeutically. In this context, a patient is matched to a specific drug or therapy based on the tumor's underlying genetic driver events rather than the histologic classification. This approach requires extensive bioinformatics methodology and workflows, including raw sequencing data processing and quality control, variant calling and annotation, integration of different molecular data types, visualization and finally reporting the data to physicians, cancer researchers and pharmacologists in a format that is readily interpretable for clinical decision making. This review comprises a broad overview of these bioinformatics aspects and discusses the multiple analytical, technical and interpretational challenges that remain to efficiently translate molecular findings into personalized treatment recommendations.

Screening and identification of hub genes in bladder cancer by bioinformatics analysis and KIF11 is a potential prognostic biomarker

Bladder cancer (BC) is the ninth most common lethal malignancy worldwide. Great efforts have been devoted to clarify the pathogenesis of BC, but the underlying molecular mechanisms remain unclear. To screen for the genes associated with the progression and carcinogenesis of BC, three datasets were obtained from the Gene Expression Omnibus. A total of 37 tumor and 16 non-cancerous samples were analyzed to identify differentially expressed genes (DEGs). Subsequently, 141 genes were identified, including 55 upregulated and 86 downregulated genes. The protein-protein interaction network was established using the Search Tool for Retrieval of Interacting Genes database. Hub gene identification and module analysis were performed using Cytoscape software. Hierarchical clustering of hub genes was conducted using the University of California, Santa Cruz Cancer Genomics Browser. Among the hub genes, kinesin family member 11 (KIF11) was identified as one of the most significant prognostic biomarkers among all the candidates. The Kaplan Meier Plotter database was used for survival analysis of KIF11. The expression profile of KIF11 was analyzed using the ONCOMINE database. The expression levels of KIF11 in BC samples and bladder cells were measured using reverse transcription-quantitative pCR, immunohistochemistry and western blotting. In summary, KIF11 was significantly upregulated in BC and might act as a potential prognostic biomarker. The present identification of DEGs and hub genes in BC may provide novel insight for investigating the molecular mechanisms of BC.

Establishment of a Molecular Tumor Board (MTB) and Uptake of Recommendations in a Community Setting

In the precision medicine era, molecular testing in advanced cancer is foundational to patient management. Molecular tumor boards (MTBs) can be effective in processing comprehensive genomic profiling (CGP) results and providing expert recommendations. We assessed an MTB and its role in a community setting. This retrospective analysis included patients with MTB recommendations at a community-based oncology practice January 2015 to December 2018; exclusions were death within 60 days of the MTB and/or no metastatic disease. Potentially actionable genomic alterations from CGP (immunohistochemistry, in-situ hybridization, next-generation sequencing) were reviewed bi-weekly by MTB practice experts, pathologists, genetic counselors, and other support staff, and clinical care recommendations were provided. Subsequent chart reviews determined implementation rates of recommendations. In 613 patients, the most common cancers were lung (23%), breast (19%), and colorectal (17%); others included ovarian, endometrial, bladder, and melanoma. Patients received 837 actionable recommendations: standard therapy (37%), clinical trial (31%), germline testing and genetic counseling (17%), off-label therapy (10%), subspecialty multidisciplinary tumor board review (2%), and advice for classifying tumor of unknown origin (2%). Of these recommendations, 36% to 78% were followed by the treating physician. For clinical trial recommendations (n = 262), 13% of patients enrolled in a clinical trial. The median time between CPG result availability and MTB presentation was 12 days. A community oncology-based comprehensive and high-throughput MTB provided useful clinical guidance in various treatment domains within an acceptable timeframe for patients with cancer in a large community setting.

Differential impact of cognitive computing augmented by real world evidence on novice and expert oncologists

INTRODUCTION

Cognitive computing point-of-care decision support tools which ingest patient attributes from electronic health records and display treatment options based on expert training and medical literature, supplemented by real world evidence (RWE), might prove useful to expert and novice oncologists. The concordance of augmented intelligence systems with best medical practices and potential influences on physician behavior remain unknown.

METHODS

Electronic health records from 88 breast cancer patients evaluated at a USA tertiary care center were presented to subspecialist experts and oncologists focusing on other disease states with and without reviewing the IBM Watson for Oncology with Cota RWE platform.

RESULTS

The cognitive computing "recommended" option was concordant with selection by breast cancer experts in 78.5% and "for consideration" option was selected in 9.4%, yielding agreements in 87.9%. Fifty-nine percent of non-concordant responses were generated from 8% of cases. In the Cota observational database 69.3% of matched controls were treated with "recommended," 11.4% "for consideration", and 19.3% "not recommended." Without guidance from Watson for Oncology (WfO)/Cota RWE, novice oncologists chose 75.5% recommended/for consideration treatments which improved to 95.3% with WfO/Cota RWE. The novices were more likely than experts to choose a non-recommended option (P < .01) without WfO/Cota RWE and changed decisions in 39% cases.

CONCLUSIONS

Watson for Oncology with Cota RWE options were largely concordant with disease expert judged best oncology practices, and was able to improve treatment decisions among breast cancer novices. The observation that nearly a fifth of patients with similar disease characteristics received non-recommended options in a real world database highlights a need for decision support.

Review: Precision medicine and driver mutations: Computational methods, functional assays and conformational principles for interpreting cancer drivers

At the root of the so-called precision medicine or precision oncology, which is our focus here, is the hypothesis that cancer treatment would be considerably better if therapies were guided by a tumor’s genomic alterations. This hypothesis has sparked major initiatives focusing on whole-genome and/or exome sequencing, creation of large databases, and developing tools for their statistical analyses—all aspiring to identify actionable alterations, and thus molecular targets, in a patient. At the center of the massive amount of collected sequence data is their interpretations that largely rest on statistical analysis and phenotypic observations. Statistics is vital, because it guides identification of cancer-driving alterations. However, statistics of mutations do not identify a change in protein conformation; therefore, it may not define sufficiently accurate actionable mutations, neglecting those that are rare. Among the many thematic overviews of precision oncology, this review innovates by further comprehensively including precision pharmacology, and within this framework, articulating its protein structural landscape and consequences to cellular signaling pathways. It provides the underlying physicochemical basis, thereby also opening the door to a broader community.

SwissMTB: establishing comprehensive molecular cancer diagnostics in Swiss clinics

Background

Molecular precision oncology is an emerging practice to improve cancer therapy by decreasing the risk of choosing treatments that lack efficacy or cause adverse events. However, the challenges of integrating molecular profiling into routine clinical care are manifold. From a computational perspective these include the importance of a short analysis turnaround time, the interpretation of complex drug-gene and gene-gene interactions, and the necessity of standardized high-quality workflows. In addition, difficulties faced when integrating molecular diagnostics into clinical practice are ethical concerns, legal requirements, and limited availability of treatment options beyond standard of care as well as the overall lack of awareness of their existence.

Methods

To the best of our knowledge, we are the first group in Switzerland that established a workflow for personalized diagnostics based on comprehensive high-throughput sequencing of tumors at the clinic. Our workflow, named SwissMTB (Swiss Molecular Tumor Board), links genetic tumor alterations and gene expression to therapeutic options and clinical trial opportunities. The resulting treatment recommendations are summarized in a clinical report and discussed in a molecular tumor board at the clinic to support therapy decisions.

Results

Here we present results from an observational pilot study including 22 late-stage cancer patients. In this study we were able to identify actionable variants and corresponding therapies for 19 patients. Half of the patients were analyzed retrospectively. In two patients we identified resistance-associated variants explaining lack of therapy response. For five out of eleven patients analyzed before treatment the SwissMTB diagnostic influenced treatment decision.

Conclusions

SwissMTB enables the analysis and clinical interpretation of large numbers of potentially actionable molecular targets. Thus, our workflow paves the way towards a more frequent use of comprehensive molecular diagnostics in Swiss hospitals.

Electronic supplementary material

The online version of this article (10.1186/s12911-018-0680-0) contains supplementary material, which is available to authorized users.

Personalized Clinical Decision Making Through Implementation of a Molecular Tumor Board: A German Single-Center Experience

PURPOSE

Dramatic advances in our understanding of the molecular pathophysiology of cancer, along with a rapidly expanding portfolio of molecular targeted drugs, have led to a paradigm shift toward personalized, biomarker-driven cancer treatment. Here, we report the 2-year experience of the Comprehensive Cancer Center Freiburg Molecular Tumor Board (MTB), one of the first interdisciplinary molecular tumor conferences established in Europe. The role of the MTB is to recommend personalized therapy for patients with cancer beyond standard-of-care treatment.

METHODS

This retrospective case series includes 198 patients discussed from March 2015 through February 2017. The MTB guided individual molecular diagnostics, assessed evidence of actionability of molecular alterations, and provided therapy recommendations, including approved and off-label treatments as well as available matched clinical trials.

RESULTS

The majority of patients had metastatic solid tumors (73.7%), mostly progressive (77.3%) after a mean of 2.0 lines of standard treatment. Diagnostic recommendations resulted in 867 molecular diagnostic tests for 172 patients (five per case), including exome analysis in 36 cases (18.2%). With a median turnaround time of 28 days, treatment recommendations were given to 104 patients (52.5%). These included single-agent targeted therapies (42.3%), checkpoint inhibitors (37.5%), and combination therapies (18.3%). Treatment recommendations were implemented in 33 of 104 patients (31.7%), of whom 19 (57.6%) showed stable disease or partial response, including 14 patients (7.1% of the entire population) receiving off-label treatments.

CONCLUSION

Personalized extended molecular-guided patient care is effective for a small but clinically meaningful proportion of patients in challenging clinical situations. Limited access to targeted drugs, lack of trials, and submission at late disease stage prevents broader applicability, whereas genome-wide analyses are not a strict requirement for predictive molecular testing.

The overlap between regeneration and fibrosis in injured skeletal muscle is regulated by phosphatidylinositol 3-kinase/Akt signaling pathway - A bioinformatic analysis based on lncRNA microarray

Injured skeletal muscle would go through a sequence of the pathological phases of degeneration, myogenesis and fibrosis. Growing evidence indicated that fibrotic and myogenic phases might overlap within the injured skeletal muscle in the early time after injury. However, the mechanism underlying this overlapping remains unclear. Here, we performed an lncRNA microarray to identify the activated pathways in mice muscle seven days after contusion. KEGG analysis indicated that phosphatidylinositol 3-kinase/Akt (PI3K/Akt) signaling cascade was predicted to be activated by lncRNAs. The top genes targeted by lncRNAs in PI3K/Akt signaling were subunits of laminin, collagen 5, and collagen 6, which participated in either myogenic or fibrotic process. Reverse transcriptase-polymerase chain reaction analysis and immunohistochemical stain further confirmed the prediction in silico. These results suggested that the overlap might be related to an activated PI3K/Akt pathway by lncRNA regulation.