Optimizing cancer genome sequencing and analysis

Genomic Mosaicism of the Brain: Origin, Impact, and Utility

Genomic mosaicism describes the phenomenon where some but not all cells within a tissue harbor unique genetic mutations. Traditionally, research focused on the impact of genomic mosaicism on clinical phenotype-motivated by its involvement in cancers and overgrowth syndromes. More recently, we increasingly shifted towards the plethora of neutral mosaic variants that can act as recorders of cellular lineage and environmental exposures. Here, we summarize the current state of the field of genomic mosaicism research with a special emphasis on our current understanding of this phenomenon in brain development and homeostasis. Although the field of genomic mosaicism has a rich history, technological advances in the last decade have changed our approaches and greatly improved our knowledge. We will provide current definitions and an overview of contemporary detection approaches for genomic mosaicism. Finally, we will discuss the impact and utility of genomic mosaicism.

A novel patient-derived meningioma spheroid model as a tool to study and treat epithelial-to-mesenchymal transition (EMT) in meningiomas

Meningiomas are the most common intracranial brain tumours. These tumours are heterogeneous and encompass a wide spectrum of clinical aggressivity. Treatment options are limited to surgery and radiotherapy and have a risk of post-operative morbidities and radiation neurotoxicity, reflecting the need for new therapies. Three-dimensional (3D) patient-derived cell culture models have been shown to closely recapitulate in vivo tumour biology, including microenvironmental interactions and have emerged as a robust tool for drug development. Here, we established a novel easy-to-use 3D patient-derived meningioma spheroid model using a scaffold-free approach. Patient-derived meningioma spheroids were characterised and compared to patient tissues and traditional monolayer cultures by histology, genomics, and transcriptomics studies. Patient-derived meningioma spheroids closely recapitulated morphological and molecular features of matched patient tissues, including patient histology, genomic alterations, and components of the immune microenvironment, such as a CD68 + and CD163 + positive macrophage cell population. Comprehensive transcriptomic profiling revealed an increase in epithelial-to-mesenchymal transition (EMT) in meningioma spheroids compared to traditional monolayer cultures, confirming this model as a tool to elucidate EMT in meningioma. Therefore, as proof of concept study, we developed a treatment strategy to target EMT in meningioma. We found that combination therapy using the MER tyrosine kinase (MERTK) inhibitor UNC2025 and the histone deacetylase (HDAC) inhibitor Trichostatin A (TSA) effectively decreased meningioma spheroid viability and proliferation. Furthermore, we demonstrated this combination therapy significantly increased the expression of the epithelial marker E-cadherin and had a repressive effect on WHO grade 2-derived spheroid invasion, which is suggestive of a partial reversal of EMT in meningioma spheroids.

Ultra-Deep Sequencing Reveals the Mutational Landscape of Classical Hodgkin Lymphoma

The malignant Hodgkin and Reed Sternberg (HRS) cells of classical Hodgkin lymphoma (cHL) are scarce in affected lymph nodes, creating a challenge to detect driver somatic mutations. As an alternative to cell purification techniques, we hypothesized that ultra-deep exome sequencing would allow genomic study of HRS cells, thereby streamlining analysis and avoiding technical pitfalls. To test this, 31 cHL tumor/normal pairs were exome sequenced to approximately 1,000× median depth of coverage. An orthogonal error-corrected sequencing approach verified >95% of the discovered mutations. We identified mutations in genes novel to cHL including: CDH5 and PCDH7, novel stop gain mutations in IL4R, and a novel pattern of recurrent mutations in pathways regulating Hippo signaling. As a further application of our exome sequencing, we attempted to identify expressed somatic single-nucleotide variants (SNV) in single-nuclei RNA sequencing (snRNA-seq) data generated from a patient in our cohort. Our snRNA analysis identified a clear cluster of cells containing a somatic SNV identified in our deep exome data. This cluster has differentially expressed genes that are consistent with genes known to be dysregulated in HRS cells (e.g., PIM1 and PIM3). The cluster also contains cells with an expanded B-cell clonotype further supporting a malignant phenotype. This study provides proof-of-principle that ultra-deep exome sequencing can be utilized to identify recurrent mutations in HRS cells and demonstrates the feasibility of snRNA-seq in the context of cHL. These studies provide the foundation for the further analysis of genomic variants in large cohorts of patients with cHL.

SIGNIFICANCE

Our data demonstrate the utility of ultra-deep exome sequencing in uncovering somatic variants in Hodgkin lymphoma, creating new opportunities to define the genes that are recurrently mutated in this disease. We also show for the first time the successful application of snRNA-seq in Hodgkin lymphoma and describe the expression profile of a putative cluster of HRS cells in a single patient.

Growth exponents reflect evolutionary processes and treatment response in brain metastases

Tumor growth is the result of the interplay of complex biological processes in huge numbers of individual cells living in changing environments. Effective simple mathematical laws have been shown to describe tumor growth in vitro, or simple animal models with bounded-growth dynamics accurately. However, results for the growth of human cancers in patients are scarce. Our study mined a large dataset of 1133 brain metastases (BMs) with longitudinal imaging follow-up to find growth laws for untreated BMs and recurrent treated BMs. Untreated BMs showed high growth exponents, most likely related to the underlying evolutionary dynamics, with experimental tumors in mice resembling accurately the disease. Recurrent BMs growth exponents were smaller, most probably due to a reduction in tumor heterogeneity after treatment, which may limit the tumor evolutionary capabilities. In silico simulations using a stochastic discrete mesoscopic model with basic evolutionary dynamics led to results in line with the observed data.

Variant calling and benchmarking in an era of complete human genome sequences

Genetic variant calling from DNA sequencing has enabled understanding of germline variation in hundreds of thousands of humans. Sequencing technologies and variant-calling methods have advanced rapidly, routinely providing reliable variant calls in most of the human genome. We describe how advances in long reads, deep learning, de novo assembly and pangenomes have expanded access to variant calls in increasingly challenging, repetitive genomic regions, including medically relevant regions, and how new benchmark sets and benchmarking methods illuminate their strengths and limitations. Finally, we explore the possible future of more complete characterization of human genome variation in light of the recent completion of a telomere-to-telomere human genome reference assembly and human pangenomes, and we consider the innovations needed to benchmark their newly accessible repetitive regions and complex variants.

Somatic mosaicism in inborn errors of immunity: Current knowledge, challenges, and future perspectives

Inborn errors of immunity (IEI) are a diverse group of monogenic disorders of the immune system due to germline variants in genes important for the immune response. Over the past decade there has been increasing recognition that acquired somatic variants present in a subset of cells can also lead to immune disorders or 'phenocopies' of IEI. Discovery of somatic mosaicism causing IEI has largely arisen from investigation of seemingly sporadic cases of IEI with predominant symptoms of autoinflammation and/or autoimmunity in which germline disease-causing variants are not detected. Disease-causing somatic mosaicism has been identified in genes that also cause germline IEI, such as FAS, and in genes without significant corresponding germline disease, such as UBA1 and TLR8. There are challenges in detecting low-level somatic variants, and it is likely that the extent of the somatic mosaicism causing IEI is largely uncharted. Here we review the field of somatic mosaicism leading to IEI and discuss challenges and methods for somatic variant detection, including diagnostic approaches for molecular diagnoses of patients.

Clonal interactions in cancer: integrating quantitative models with experimental and clinical data

Tumors consist of different genotypically distinct subpopulations-or subclones-of cells. These subclones can influence neighboring clones in a process called "clonal interaction." Conventionally, research on driver mutations in cancer has focused on their cell-autonomous effects that lead to an increase in fitness of the cells containing the driver. Recently, with the advent of improved experimental and computational technologies for investigating tumor heterogeneity and clonal dynamics, new studies have shown the importance of clonal interactions in cancer initiation, progression, and metastasis. In this review we provide an overview of clonal interactions in cancer, discussing key discoveries from a diverse range of approaches to cancer biology research. We discuss common types of clonal interactions, such as cooperation and competition, its mechanisms, and the overall effect on tumorigenesis, with important implications for tumor heterogeneity, resistance to treatment, and tumor suppression. Quantitative models-in coordination with cell culture and animal model experiments-have played a vital role in investigating the nature of clonal interactions and the complex clonal dynamics they generate. We present mathematical and computational models that can be used to represent clonal interactions and provide examples of the roles they have played in identifying and quantifying the strength of clonal interactions in experimental systems. Clonal interactions have proved difficult to observe in clinical data; however, several very recent quantitative approaches enable their detection. We conclude by discussing ways in which researchers can further integrate quantitative methods with experimental and clinical data to elucidate the critical-and often surprising-roles of clonal interactions in human cancers.

Artificial intelligence assists precision medicine in cancer treatment

Cancer is a major medical problem worldwide. Due to its high heterogeneity, the use of the same drugs or surgical methods in patients with the same tumor may have different curative effects, leading to the need for more accurate treatment methods for tumors and personalized treatments for patients. The precise treatment of tumors is essential, which renders obtaining an in-depth understanding of the changes that tumors undergo urgent, including changes in their genes, proteins and cancer cell phenotypes, in order to develop targeted treatment strategies for patients. Artificial intelligence (AI) based on big data can extract the hidden patterns, important information, and corresponding knowledge behind the enormous amount of data. For example, the ML and deep learning of subsets of AI can be used to mine the deep-level information in genomics, transcriptomics, proteomics, radiomics, digital pathological images, and other data, which can make clinicians synthetically and comprehensively understand tumors. In addition, AI can find new biomarkers from data to assist tumor screening, detection, diagnosis, treatment and prognosis prediction, so as to providing the best treatment for individual patients and improving their clinical outcomes.

From Samples to Germline and Somatic Sequence Variation: A Focus on Next-Generation Sequencing in Melanoma Research

Next-generation sequencing (NGS) applications have flourished in the last decade, permitting the identification of cancer driver genes and profoundly expanding the possibilities of genomic studies of cancer, including melanoma. Here we aimed to present a technical review across many of the methodological approaches brought by the use of NGS applications with a focus on assessing germline and somatic sequence variation. We provide cautionary notes and discuss key technical details involved in library preparation, the most common problems with the samples, and guidance to circumvent them. We also provide an overview of the sequence-based methods for cancer genomics, exposing the pros and cons of targeted sequencing vs. exome or whole-genome sequencing (WGS), the fundamentals of the most common commercial platforms, and a comparison of throughputs and key applications. Details of the steps and the main software involved in the bioinformatics processing of the sequencing results, from preprocessing to variant prioritization and filtering, are also provided in the context of the full spectrum of genetic variation (SNVs, indels, CNVs, structural variation, and gene fusions). Finally, we put the emphasis on selected bioinformatic pipelines behind (a) short-read WGS identification of small germline and somatic variants, (b) detection of gene fusions from transcriptomes, and (c) de novo assembly of genomes from long-read WGS data. Overall, we provide comprehensive guidance across the main methodological procedures involved in obtaining sequencing results for the most common short- and long-read NGS platforms, highlighting key applications in melanoma research.

Inferring parameters of cancer evolution in chronic lymphocytic leukemia

As a cancer develops, its cells accrue new mutations, resulting in a heterogeneous, complex genomic profile. We make use of this heterogeneity to derive simple, analytic estimates of parameters driving carcinogenesis and reconstruct the timeline of selective events following initiation of an individual cancer, where two longitudinal samples are available for sequencing. Using stochastic computer simulations of cancer growth, we show that we can accurately estimate mutation rate, time before and after a driver event occurred, and growth rates of both initiated cancer cells and subsequently appearing subclones. We demonstrate that in order to obtain accurate estimates of mutation rate and timing of events, observed mutation counts should be corrected to account for clonal mutations that occurred after the founding of the tumor, as well as sequencing coverage. Chronic lymphocytic leukemia (CLL), which often does not require treatment for years after diagnosis, presents an optimal system to study the untreated, natural evolution of cancer cell populations. When we apply our methodology to reconstruct the individual evolutionary histories of CLL patients, we find that the parental leukemic clone typically appears within the first fifteen years of life.

FVC as an adaptive and accurate method for filtering variants from popular NGS analysis pipelines

The quality control of variants from whole-genome sequencing data is vital in clinical diagnosis and human genetics research. However, current filtering methods (Frequency, Hard-Filter, VQSR, GARFIELD, and VEF) were developed to be utilized on particular variant callers and have certain limitations. Especially, the number of eliminated true variants far exceeds the number of removed false variants using these methods. Here, we present an adaptive method for quality control on genetic variants from different analysis pipelines, and validate it on the variants generated from four popular variant callers (GATK HaplotypeCaller, Mutect2, Varscan2, and DeepVariant). FVC consistently exhibited the best performance. It removed far more false variants than the current state-of-the-art filtering methods and recalled ~51-99% true variants filtered out by the other methods. Once trained, FVC can be conveniently integrated into a user-specific variant calling pipeline.

FVC is a method for calling specific gene variants from whole genome data, for potential use in clinical diagnosis and human genetics research.

Accurate Identification of Subclones in Tumor Genomes

Understanding intra-tumor heterogeneity is critical for studying tumorigenesis and designing personalized treatments. To decompose the mixed cell population in a tumor, subclones are inferred computationally based on variant allele frequency (VAF) from bulk sequencing data. In this study, we showed that sequencing depth, mean VAF, and variance of VAF of a subclone are confounded. Without considering this effect, current methods require deep-sequencing data (>300x depth) to reliably infer subclones. Here we present a novel algorithm that incorporates depth-variance and mean-variance dependencies in a clustering error model and successfully identifies subclones in tumors sequenced at depths of as low as 30x. We implemented the algorithm as a model-based adaptive grouping of subclones (MAGOS) method. Analyses of computer simulated data and empirical sequencing data showed that MAGOS outperformed existing methods on minimum sequencing depth, decomposition accuracy, and computation efficiency. The most prominent improvements were observed in analyzing tumors sequenced at depths between 30x and 200x, while the performance was comparable between MAGOS and existing methods on deeply sequenced tumors. MAGOS supports analysis of single nucleotide variants and copy number variants from a single sample or multiple samples of a tumor. We applied MAGOS to whole-exome data of late-stage liver cancers and discovered that high subclone count in a tumor was a significant risk factor of poor prognosis. Lastly, our analysis suggested that sequencing multiple samples of the same tumor at standard depth is more cost-effective and robust for subclone characterization than deep sequencing a single sample. MAGOS is available at github (https://github.com/liliulab/magos).

Computational analysis of cancer genome sequencing data

Distilling biologically meaningful information from cancer genome sequencing data requires comprehensive identification of somatic alterations using rigorous computational methods. As the amount and complexity of sequencing data have increased, so has the number of tools for analysing them. Here, we describe the main steps involved in the bioinformatic analysis of cancer genomes, review key algorithmic developments and highlight popular tools and emerging technologies. These tools include those that identify point mutations, copy number alterations, structural variations and mutational signatures in cancer genomes. We also discuss issues in experimental design, the strengths and limitations of sequencing modalities and methodological challenges for the future.

An overview of mutational and copy number signatures in human cancer

The genome of each cell in the human body is constantly under assault from a plethora of exogenous and endogenous processes that can damage DNA. If not successfully repaired, DNA damage generally becomes permanently imprinted in cells, and all their progenies, as somatic mutations. In most cases, the patterns of these somatic mutations contain the tell‐tale signs of the mutagenic processes that have imprinted and are termed mutational signatures. Recent pan‐cancer genomic analyses have elucidated the compendium of mutational signatures for all types of small mutational events, including (1) single base substitutions, (2) doublet base substitutions, and (3) small insertions/deletions. In contrast to small mutational events, where, in most cases, DNA damage is a prerequisite, aneuploidy, which refers to the abnormal number of chromosomes in a cell, usually develops from mistakes during DNA replication. Such mistakes include DNA replication stress, mitotic errors caused by faulty microtubule dynamics, or cohesion defects that contribute to chromosomal breakage and can lead to copy number (CN) alterations (CNAs) or even to structural rearrangements. These aberrations also leave behind genomic scars which can be inferred from sequencing as CN signatures and rearrangement signatures. The analyses of mutational signatures of small mutational events have been extensively reviewed, so we will not comprehensively re‐examine them here. Rather, our focus will be on summarising the existing knowledge for mutational signatures of CNAs. As studying CN signatures is an emerging field, we briefly summarise the utility that mutational signatures of small mutational events have provided in basic science, cancer treatment, and cancer prevention, and we emphasise the future role that CN signatures may play in each of these fields. © 2022 The Authors. The Journal of Pathology published by John Wiley & Sons Ltd on behalf of The Pathological Society of Great Britain and Ireland.

Inferring ongoing cancer evolution from single tumour biopsies using synthetic supervised learning

Variant allele frequencies (VAF) encode ongoing evolution and subclonal selection in growing tumours. However, existing methods that utilize VAF information for cancer evolutionary inference are compressive, slow, or incorrectly specify the underlying cancer evolutionary dynamics. Here, we provide a proof-of-principle synthetic supervised learning method, TumE, that integrates simulated models of cancer evolution with Bayesian neural networks, to infer ongoing selection in bulk-sequenced single tumour biopsies. Analyses in synthetic and patient tumours show that TumE significantly improves both accuracy and inference time per sample when detecting positive selection, deconvoluting selected subclonal populations, and estimating subclone frequency. Importantly, we show how transfer learning can leverage stored knowledge within TumE models for related evolutionary inference tasks-substantially reducing data and computational time for further model development and providing a library of recyclable deep learning models for the cancer evolution community. This extensible framework provides a foundation and future directions for harnessing progressive computational methods for the benefit of cancer genomics and, in turn, the cancer patient.

Cancer Neoantigens: Challenges and Future Directions for Prediction, Prioritization, and Validation

Prioritization of immunogenic neoantigens is key to enhancing cancer immunotherapy through the development of personalized vaccines, adoptive T cell therapy, and the prediction of response to immune checkpoint inhibition. Neoantigens are tumor-specific proteins that allow the immune system to recognize and destroy a tumor. Cancer immunotherapies, such as personalized cancer vaccines, adoptive T cell therapy, and immune checkpoint inhibition, rely on an understanding of the patient-specific neoantigen profile in order to guide personalized therapeutic strategies. Genomic approaches to predicting and prioritizing immunogenic neoantigens are rapidly expanding, raising new opportunities to advance these tools and enhance their clinical relevance. Predicting neoantigens requires acquisition of high-quality samples and sequencing data, followed by variant calling and variant annotation. Subsequently, prioritizing which of these neoantigens may elicit a tumor-specific immune response requires application and integration of tools to predict the expression, processing, binding, and recognition potentials of the neoantigen. Finally, improvement of the computational tools is held in constant tension with the availability of datasets with validated immunogenic neoantigens. The goal of this review article is to summarize the current knowledge and limitations in neoantigen prediction, prioritization, and validation and propose future directions that will improve personalized cancer treatment.

Chromosomal imbalances detected via RNA-sequencing in 28 cancers

MOTIVATION

RNA-sequencing (RNA-seq) of tumor tissue is typically only used to measure gene expression. Here, we present a statistical approach that leverages existing RNA-seq data to also detect somatic copy number alterations (SCNAs), a pervasive phenomenon in human cancers, without a need to sequence the corresponding DNA.

RESULTS

We present an analysis of 4942 participant samples from 28 cancers in The Cancer Genome Atlas (TCGA), demonstrating robust detection of SCNAs from RNA-seq. Using genotype imputation and haplotype information, our RNA-based method had a median sensitivity of 85% to detect SCNAs defined by DNA analysis, at high specificity (∼95%). As an example of translational potential, we successfully replicated SCNA features associated with breast cancer subtypes. Our results credential haplotype-based inference based on RNA-seq to detect SCNAs in clinical and population-based settings.

AVAILABILITY AND IMPLEMENTATION

The analyses presented use the data publicly available from TCGA Research Network (http://cancergenome.nih.gov/). See Methods for details regarding data downloads. hapLOHseq software is freely available under The MIT license and can be downloaded from http://scheet.org/software.html.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Genetic and Transcriptional Contributions to Relapse in Normal Karyotype Acute Myeloid Leukemia

To better understand clonal and transcriptional adaptations after relapse in patients with acute myeloid leukemia (AML), we collected presentation and relapse samples from six normal karyotype AML cases. We performed enhanced whole-genome sequencing to characterize clonal evolution, and deep-coverage single-cell RNA sequencing on the same samples, which yielded 142,642 high-quality cells for analysis. Identifying expressed mutations in individual cells enabled us to discriminate between normal and AML cells, to identify coordinated changes in the genome and transcriptome, and to identify subclone-specific cell states. We quantified the coevolution of genetic and transcriptional heterogeneity during AML progression, and found that transcriptional changes were significantly correlated with genetic changes. However, transcriptional adaptation sometimes occurred independently, suggesting that clonal evolution does not represent all relevant biological changes. In three cases, we identified cells at diagnosis that likely seeded the relapse. Finally, these data revealed a conserved relapse-enriched leukemic cell state bearing markers of stemness, quiescence, and adhesion. SIGNIFICANCE: These data enabled us to identify a relapse-enriched leukemic cell state with distinct transcriptional properties. Detailed case-by-case analyses elucidated the complex ways in which the AML genome, transcriptome, and immune microenvironment interact to evade chemotherapy. These analyses provide a blueprint for evaluating these factors in larger cohorts.This article is highlighted in the In This Issue feature, p. 1.

OUP accepted manuscript

Cell therapy using induced pluripotent stem cell (iPSC) derivatives may result in abnormal tissue generation because the cells undergo numerous cycles of mitosis before clinical application, potentially increasing the accumulation of genetic abnormalities. Therefore, genetic tests may predict abnormal tissue formation after transplantation. Here, we administered iPSC derivatives with or without single-nucleotide variants (SNVs) and deletions in cancer-related genes with various genomic copy number variant (CNV) profiles into immunodeficient mice and examined the relationships between mutations and abnormal tissue formation after transplantation. No positive correlations were found between the presence of SNVs/deletions and the formation of abnormal tissues; the overall predictivity was 29%. However, a copy number higher than 3 was correlated, with an overall predictivity of 86%. Furthermore, we found CNV hotspots at 14q32.33 and 17q12 loci. Thus, CNV analysis may predict abnormal tissue formation after transplantation of iPSC derivatives and reduce the number of tumorigenicity tests.

Measurable Residual Disease Assessment as a Surrogate Marker in New Drug Development in Acute Myeloid Leukemia

ABSTRACT

Response criteria for patients treated for acute myeloid leukemia (AML) based on cytomorphology are inadequate. Many patients achieving a complete remission by such criteria will later relapse. Patients with AML in such remissions who test negative using higher sensitivity measures of residual disease burden (measurable residual disease [MRD]) have on average lower relapse rates and better survival than those testing positive. This association has raised the possibility that these technological advances in measurement of tumor burden could be used to optimize the drug development and regulatory approval processes in AML. The heterogeneous genetic etiology, diverse immunophenotypic profiles, related precursor states and polyclonal architecture however combine to make the development of standardized and validated MRD assessments for AML challenging. Current and future methods to measure residual disease in AML, performance characteristics of testing currently in use, and potential uses for optimized AML MRD tests including as a surrogate endpoint are discussed.

Tracing the evolution of aneuploid cancers by multiregional sequencing with CRUST

Clonal deconvolution of mutational landscapes is crucial to understand the evolutionary dynamics of cancer. Two limiting factors for clonal deconvolution that have remained unresolved are variation in purity and chromosomal copy number across different samples of the same tumor. We developed a semi-supervised algorithm that tracks variant calls through multi-sample spatiotemporal tumor data. While normalizing allele frequencies based on purity, it also adjusts for copy number changes at clonal deconvolution. Absent à priori copy number data, it renders in silico copy number estimations from bulk sequences. Using published and simulated tumor sequences, we reliably segregated clonal/subclonal variants even at a low sequencing depth (~50×). Given at least one pure tumor sample (>70% purity), we could normalize and deconvolve paired samples down to a purity of 40%. This renders a reliable clonal reconstruction well adapted to multi-regionally sampled solid tumors, which are often aneuploid and contaminated by non-cancer cells.

Genomic and Transcriptomic Characterization of Relapsed SCLC Through Rapid Research Autopsy

Introduction

Relapsed SCLC is characterized by therapeutic resistance and high mortality rate. Despite decades of research, mechanisms responsible for therapeutic resistance have remained elusive owing to limited tissues available for molecular studies. Thus, an unmet need remains for molecular characterization of relapsed SCLC to facilitate development of effective therapies.

Methods

We performed whole-exome and transcriptome sequencing of metastatic tumor samples procured from research autopsies of five patients with relapsed SCLC. We implemented bioinformatics tools to infer subclonal phylogeny and identify recurrent genomic alterations. We implemented immune cell signature and single-sample gene set enrichment analyses on tumor and normal transcriptome data from autopsy and additional primary and relapsed SCLC data sets. Furthermore, we evaluated T cell-inflamed gene expression profiles in neuroendocrine (ASCL1, NEUROD1) and non-neuroendocrine (YAP1, POU2F3) SCLC subtypes.

Results

Exome sequencing revealed clonal heterogeneity (intertumor and intratumor) arising from branched evolution and identified resistance-associated truncal and subclonal alterations in relapsed SCLC. Transcriptome analyses further revealed a noninflamed phenotype in neuroendocrine SCLC subtypes (ASCL1, NEUROD1) associated with decreased expression of genes involved in adaptive antitumor immunity whereas non-neuroendocrine subtypes (YAP1, POU2F3) revealed a more inflamed phenotype.

Conclusions

Our results reveal substantial tumor heterogeneity and complex clonal evolution in relapsed SCLC. Furthermore, we report that neuroendocrine SCLC subtypes are immunologically cold, thus explaining decreased responsiveness to immune checkpoint blockade. These results suggest that the mechanisms of innate and acquired therapeutic resistances are subtype-specific in SCLC and highlight the need for continued investigation to bolster therapy selection and development for this cancer.

Toward best practice in cancer mutation detection with whole-genome and whole-exome sequencing

Clinical applications of precision oncology require accurate tests that can distinguish true cancer-specific mutations from errors introduced at each step of next-generation sequencing (NGS). To date, no bulk sequencing study has addressed the effects of cross-site reproducibility, nor the biological, technical and computational factors that influence variant identification. Here we report a systematic interrogation of somatic mutations in paired tumor-normal cell lines to identify factors affecting detection reproducibility and accuracy at six different centers. Using whole-genome sequencing (WGS) and whole-exome sequencing (WES), we evaluated the reproducibility of different sample types with varying input amount and tumor purity, and multiple library construction protocols, followed by processing with nine bioinformatics pipelines. We found that read coverage and callers affected both WGS and WES reproducibility, but WES performance was influenced by insert fragment size, genomic copy content and the global imbalance score (GIV; G > T/C > A). Finally, taking into account library preparation protocol, tumor content, read coverage and bioinformatics processes concomitantly, we recommend actionable practices to improve the reproducibility and accuracy of NGS experiments for cancer mutation detection.

Conifer: clonal tree inference for tumor heterogeneity with single-cell and bulk sequencing data

BACKGROUND

Genetic heterogeneity of a cancer tumor that develops during clonal evolution is one of the reasons for cancer treatment failure, by increasing the chance of drug resistance. Clones are cell populations with different genotypes, resulting from differences in somatic mutations that occur and accumulate during cancer development. An appropriate approach for identifying clones is determining the variant allele frequency of mutations that occurred in the tumor. Although bulk sequencing data can be used to provide that information, the frequencies are not informative enough for identifying different clones with the same prevalence and their evolutionary relationships. On the other hand, single-cell sequencing data provides valuable information about branching events in the evolution of a cancerous tumor. However, the temporal order of mutations may be determined with ambiguities using only single-cell data, while variant allele frequencies from bulk sequencing data can provide beneficial information for inferring the temporal order of mutations with fewer ambiguities.

RESULT

In this study, a new method called Conifer (ClONal tree Inference For hEterogeneity of tumoR) is proposed which combines aggregated variant allele frequency from bulk sequencing data with branching event information from single-cell sequencing data to more accurately identify clones and their evolutionary relationships. It is proven that the accuracy of clone identification and clonal tree inference is increased by using Conifer compared to other existing methods on various sets of simulated data. In addition, it is discussed that the evolutionary tree provided by Conifer on real cancer data sets is highly consistent with information in both bulk and single-cell data.

CONCLUSIONS

In this study, we have provided an accurate and robust method to identify clones of tumor heterogeneity and their evolutionary history by combining single-cell and bulk sequencing data.

Clonal Evolution of Multiple Myeloma-Clinical and Diagnostic Implications

Plasma cell dyscrasias are a heterogeneous group of diseases characterized by the expansion of bone marrow plasma cells. Malignant transformation of plasma cells depends on the continuity of events resulting in a sequence of well-defined disease stages, from monoclonal gammopathy of undetermined significance (MGUS) through smoldering myeloma (SMM) to symptomatic multiple myeloma (MM). Evolution of a pre-malignant cell into a malignant cell, as well as further tumor progression, dissemination, and relapse, require development of multiple driver lesions conferring selective advantage of the dominant clone and allowing subsequent evolution under selective pressure of microenvironment and treatment. This process of natural selection facilitates tumor plasticity leading to the formation of genetically complex and heterogenous tumors that are notoriously difficult to treat. Better understanding of the mechanisms underlying tumor evolution in MM and identification of lesions driving the evolution from the premalignant clone is therefore a key to development of effective treatment and long-term disease control. Here, we review recent advances in clonal evolution patterns and genomic landscape dynamics of MM, focusing on their clinical implications.

Bam-readcount - rapid generation of basepair-resolution sequence metrics

Bam-readcount is a utility for generating low-level information about sequencing data at specific nucleotide positions. Originally designed to help filter genomic mutation calls, the metrics it outputs are useful as input for variant detection tools and for resolving ambiguity between variant callers1,2. In addition, it has found broad applicability in diverse fields including tumor evolution, single-cell genomics, climate change ecology, and tracking community spread of SARS-CoV-2.3-6.

Pooled CRISPR screening in pancreatic cancer cells implicates co-repressor complexes as a cause of multiple drug resistance via regulation of epithelial-to-mesenchymal transition

Background

Pancreatic ductal adenocarcinoma (PDAC) patients suffer poor outcomes, including a five-year survival of below 10%. Poor outcomes result in part from therapeutic resistance that limits the impact of cytotoxic first-line therapy. Novel therapeutic approaches are needed, but currently no targeted therapies exist to treat PDAC.

Methods

To assess cellular resistance mechanisms common to four cytotoxic chemotherapies (gemcitabine, 5-fluorouracil, irinotecan, and oxaliplatin) used to treat PDAC patients, we performed four genome-wide CRISPR activation (CRISPR_act) and CRISPR knock-out (CRISPR_ko) screens in two common PDAC cell lines (Panc-1 and BxPC3). We used pathway analysis to identify gene sets enriched among our hits and conducted RNA-sequencing and chromatin immunoprecipitation-sequencing (ChIP-seq) to characterize top hits from our screen. We used scratch assays to assess changes in cellular migration with HDAC1 overexpression.

Results

Our data revealed activation of ABCG2, a well-described efflux pump, as the most consistent mediator of resistance in each of our screens. CRISPR-mediated activation of genes involved in transcriptional co-repressor complexes also conferred resistance to multiple drugs. Expression of many of these genes, including HDAC1, is associated with reduced survival in PDAC patients. Up-regulation of HDAC1 in vitro increased promoter occupancy and expression of several genes involved in the epithelial-to-mesenchymal transition (EMT). These cells also displayed phenotypic changes in cellular migration consistent with activation of the EMT pathway. The expression changes resulting from HDAC1 activation were also observed with activation of several other co-repressor complex members. Finally, we developed a publicly available analysis tool, PancDS, which integrates gene expression profiles with our screen results to predict drug sensitivity in resected PDAC tumors and cell lines.

Conclusion

Our results provide a comprehensive resource for identifying cellular mechanisms of drug resistance in PDAC, mechanistically implicate HDAC1, and co-repressor complex members broadly, in multi-drug resistance, and provide an analytical tool for predicting treatment response in PDAC tumors and cell lines.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12885-021-08388-1.

MesKit: a tool kit for dissecting cancer evolution of multi-region tumor biopsies through somatic alterations

Background

Multi-region sequencing (MRS) has been widely used to analyze intra-tumor heterogeneity (ITH) and cancer evolution. However, comprehensive analysis of mutational data from MRS is still challenging, necessitating complicated integration of a plethora of computational and statistical approaches.

Findings

Here, we present MesKit, an R/Bioconductor package that can assist in characterizing genetic ITH and tracing the evolutionary history of tumors based on somatic alterations detected by MRS. MesKit provides a wide range of analysis and visualization modules, including ITH evaluation, metastatic route inference, and mutational signature identification. In addition, MesKit implements an auto-layout algorithm to generate phylogenetic trees based on somatic mutations. The application of MesKit for 2 reported MRS datasets of hepatocellular carcinoma and colorectal cancer identified known heterogeneous features and evolutionary patterns, together with potential driver events during cancer evolution.

Conclusions

In summary, MesKit is useful for interpreting ITH and tracing evolutionary trajectory based on MRS data. MesKit is implemented in R and available at https://bioconductor.org/packages/MesKit under the GPL v3 license.

Analytical validation and performance characteristics of a 48-gene next-generation sequencing panel for detecting potentially actionable genomic alterations in myeloid neoplasms

Identification of genomic mutations by molecular testing plays an important role in diagnosis, prognosis, and treatment of myeloid neoplasms. Next-generation sequencing (NGS) is an efficient method for simultaneous detection of clinically significant genomic mutations with high sensitivity. Various NGS based in-house developed and commercial myeloid neoplasm panels have been integrated into routine clinical practice. However, some genes frequently mutated in myeloid malignancies are particularly difficult to sequence with NGS panels (e.g., CEBPA, CARL, and FLT3). We report development and validation of a 48-gene NGS panel that includes genes that are technically challenging for molecular profiling of myeloid neoplasms including acute myeloid leukemia (AML), myelodysplastic syndrome (MDS), and myeloproliferative neoplasms (MPN). Target regions were captured by hybridization with complementary biotinylated DNA baits, and NGS was performed on an Illumina NextSeq500 instrument. A bioinformatics pipeline that was developed in-house was used to detect single nucleotide variations (SNVs), insertions/deletions (indels), and FLT3 internal tandem duplications (FLT3-ITD). An analytical validation study was performed on 184 unique specimens for variants with allele frequencies ≥5%. Variants identified by the 48-gene panel were compared to those identified by a 35-gene hematologic neoplasms panel using an additional 137 unique specimens. The developed assay was applied to a large cohort (n = 2,053) of patients with suspected myeloid neoplasms. Analytical validation yielded 99.6% sensitivity (95% CI: 98.9-99.9%) and 100% specificity (95% CI: 100%). Concordance of variants detected by the 2 tested panels was 100%. Among patients with suspected myeloid neoplasms (n = 2,053), 54.5% patients harbored at least one clinically significant mutation: 77% in AML patients, 48% in MDS, and 45% in MPN. Together, these findings demonstrate that the assay can identify mutations associated with diagnosis, prognosis, and treatment options of myeloid neoplasms even in technically challenging genes.

Moss enables high sensitivity single-nucleotide variant calling from multiple bulk DNA tumor samples

Intra-tumor heterogeneity renders the identification of somatic single-nucleotide variants (SNVs) a challenging problem. In particular, low-frequency SNVs are hard to distinguish from sequencing artifacts. While the increasing availability of multi-sample tumor DNA sequencing data holds the potential for more accurate variant calling, there is a lack of high-sensitivity multi-sample SNV callers that utilize these data. Here we report Moss, a method to identify low-frequency SNVs that recur in multiple sequencing samples from the same tumor. Moss provides any existing single-sample SNV caller the ability to support multiple samples with little additional time overhead. We demonstrate that Moss improves recall while maintaining high precision in a simulated dataset. On multi-sample hepatocellular carcinoma, acute myeloid leukemia and colorectal cancer datasets, Moss identifies new low-frequency variants that meet manual review criteria and are consistent with the tumor's mutational signature profile. In addition, Moss detects the presence of variants in more samples of the same tumor than reported by the single-sample caller. Moss' improved sensitivity in SNV calling will enable more detailed downstream analyses in cancer genomics.

Next Generation Sequencing Technology in the Clinic and Its Challenges

Simple Summary

Precise identification and annotation of mutations are of utmost importance in clinical oncology. Insights of the DNA sequence can provide meaningful knowledge to unravel the underlying genetics of disease. Hence, tailoring of personalized medicine often relies on specific genomic alteration for treatment efficacy. The aim of this review is to highlight that sequencing harbors much more than just four nucleotides. Moreover, the gradual transition from first to second generation sequencing technologies has led to awareness for choosing the most appropriate bioinformatic analytic tools based on the aim, quality and demand for a specific purpose. Thus, the same raw data can lead to various results reflecting the intrinsic features of different datamining pipelines.

Abstract

Data analysis has become a crucial aspect in clinical oncology to interpret output from next-generation sequencing-based testing. NGS being able to resolve billions of sequencing reactions in a few days has consequently increased the demand for tools to handle and analyze such large data sets. Many tools have been developed since the advent of NGS, featuring their own peculiarities. Increased awareness when interpreting alterations in the genome is therefore of utmost importance, as the same data using different tools can provide diverse outcomes. Hence, it is crucial to evaluate and validate bioinformatic pipelines in clinical settings. Moreover, personalized medicine implies treatment targeting efficacy of biological drugs for specific genomic alterations. Here, we focused on different sequencing technologies, features underlying the genome complexity, and bioinformatic tools that can impact the final annotation. Additionally, we discuss the clinical demand and design for implementing NGS.

Genome Sequencing as an Alternative to Cytogenetic Analysis in Myeloid Cancers

BACKGROUND

Genomic analysis is essential for risk stratification in patients with acute myeloid leukemia (AML) or myelodysplastic syndromes (MDS). Whole-genome sequencing is a potential replacement for conventional cytogenetic and sequencing approaches, but its accuracy, feasibility, and clinical utility have not been demonstrated.

METHODS

We used a streamlined whole-genome sequencing approach to obtain genomic profiles for 263 patients with myeloid cancers, including 235 patients who had undergone successful cytogenetic analysis. We adapted sample preparation, sequencing, and analysis to detect mutations for risk stratification using existing European Leukemia Network (ELN) guidelines and to minimize turnaround time. We analyzed the performance of whole-genome sequencing by comparing our results with findings from cytogenetic analysis and targeted sequencing.

RESULTS

Whole-genome sequencing detected all 40 recurrent translocations and 91 copy-number alterations that had been identified by cytogenetic analysis. In addition, we identified new clinically reportable genomic events in 40 of 235 patients (17.0%). Prospective sequencing of samples obtained from 117 consecutive patients was performed in a median of 5 days and provided new genetic information in 29 patients (24.8%), which changed the risk category for 19 patients (16.2%). Standard AML risk groups, as defined by sequencing results instead of cytogenetic analysis, correlated with clinical outcomes. Whole-genome sequencing was also used to stratify patients who had inconclusive results by cytogenetic analysis into risk groups in which clinical outcomes were measurably different.

CONCLUSIONS

In our study, we found that whole-genome sequencing provided rapid and accurate genomic profiling in patients with AML or MDS. Such sequencing also provided a greater diagnostic yield than conventional cytogenetic analysis and more efficient risk stratification on the basis of standard risk categories. (Funded by the Siteman Cancer Research Fund and others.).

Structural variant detection in cancer genomes: computational challenges and perspectives for precision oncology

Cancer is generally characterized by acquired genomic aberrations in a broad spectrum of types and sizes, ranging from single nucleotide variants to structural variants (SVs). At least 30% of cancers have a known pathogenic SV used in diagnosis or treatment stratification. However, research into the role of SVs in cancer has been limited due to difficulties in detection. Biological and computational challenges confound SV detection in cancer samples, including intratumor heterogeneity, polyploidy, and distinguishing tumor-specific SVs from germline and somatic variants present in healthy cells. Classification of tumor-specific SVs is challenging due to inconsistencies in detected breakpoints, derived variant types and biological complexity of some rearrangements. Full-spectrum SV detection with high recall and precision requires integration of multiple algorithms and sequencing technologies to rescue variants that are difficult to resolve through individual methods. Here, we explore current strategies for integrating SV callsets and to enable the use of tumor-specific SVs in precision oncology.

More for less: predicting and maximizing genomic variant discovery via Bayesian nonparametrics

While the cost of sequencing genomes has decreased dramatically in recent years, this expense often remains nontrivial. Under a fixed budget, scientists face a natural trade-off between quantity and quality: spending resources to sequence a greater number of genomes or spending resources to sequence genomes with increased accuracy. Our goal is to find the optimal allocation of resources between quantity and quality. Optimizing resource allocation promises to reveal as many new variations in the genome as possible. We introduce a Bayesian nonparametric methodology to predict the number of new variants in a follow-up study based on a pilot study. When experimental conditions are kept constant between the pilot and follow-up, we find that our prediction is competitive with the best existing methods. Unlike current methods, though, our new method allows practitioners to change experimental conditions between the pilot and the follow-up. We demonstrate how this distinction allows our method to be used for more realistic predictions and for optimal allocation of a fixed budget between quality and quantity. We validate our method on cancer and human genomics data.

A simple method to estimate the in-house limit of detection for genetic mutations with low allele frequencies in whole-exome sequencing analysis by next-generation sequencing

BACKGROUND

Next-generation sequencing (NGS) has profoundly changed the approach to genetic/genomic research. Particularly, the clinical utility of NGS in detecting mutations associated with disease risk has contributed to the development of effective therapeutic strategies. Recently, comprehensive analysis of somatic genetic mutations by NGS has also been used as a new approach for controlling the quality of cell substrates for manufacturing biopharmaceuticals. However, the quality evaluation of cell substrates by NGS largely depends on the limit of detection (LOD) for rare somatic mutations. The purpose of this study was to develop a simple method for evaluating the ability of whole-exome sequencing (WES) by NGS to detect mutations with low allele frequency. To estimate the LOD of WES for low-frequency somatic mutations, we repeatedly and independently performed WES of a reference genomic DNA using the same NGS platform and assay design. LOD was defined as the allele frequency with a relative standard deviation (RSD) value of 30% and was estimated by a moving average curve of the relation between RSD and allele frequency.

RESULTS

Allele frequencies of 20 mutations in the reference material that had been pre-validated by droplet digital PCR (ddPCR) were obtained from 5, 15, 30, or 40 G base pair (Gbp) sequencing data per run. There was a significant association between the allele frequencies measured by WES and those pre-validated by ddPCR, whose p-value decreased as the sequencing data size increased. By this method, the LOD of allele frequency in WES with the sequencing data of 15 Gbp or more was estimated to be between 5 and 10%.

CONCLUSIONS

For properly interpreting the WES data of somatic genetic mutations, it is necessary to have a cutoff threshold of low allele frequencies. The in-house LOD estimated by the simple method shown in this study provides a rationale for setting the cutoff.

A practical guide to cancer subclonal reconstruction from DNA sequencing

Subclonal reconstruction from bulk tumor DNA sequencing has become a pillar of cancer evolution studies, providing insight into the clonality and relative ordering of mutations and mutational processes. We provide an outline of the complex computational approaches used for subclonal reconstruction from single and multiple tumor samples. We identify the underlying assumptions and uncertainties in each step and suggest best practices for analysis and quality assessment. This guide provides a pragmatic resource for the growing user community of subclonal reconstruction methods.

Assessing Limit of Detection in Clinical Sequencing

Clinical reporting of solid tumor sequencing requires reliable assessment of the accuracy and reproducibility of each assay. Somatic mutation variant allele fractions may be below 10% in many samples due to sample heterogeneity, tumor clonality, and/or sample degradation in fixatives such as formalin. The toolkits available to the clinical sequencing community for correlating assay design parameters with assay sensitivity remain limited, and large-scale empirical assessments are often relied upon due to the lack of clear theoretical grounding. To address this uncertainty, a theoretical model was developed for predicting the expected variant calling sensitivity for a given library complexity and sequencing depth. Binomial models were found to be appropriate when assay sensitivity was only limited by library complexity or sequencing depth, but functional scaling for library complexity was necessary when both library complexity and sequencing depth were co-limiting. This model was empirically validated with sequencing experiments by using a series of DNA input amounts and sequencing depths. Based on these findings, a workflow is proposed for determining the limiting factors to sensitivity in different assay designs, and the formulas for these scenarios are presented. The approach described here provides designers of clinical assays with the methods to theoretically predict assay design outcomes a priori, potentially reducing burden in clinical tumor assay design and validation efforts.

Comparative analysis of target gene exon sequencing by cognitive technology using a next generation sequencing platform in patients with lung cancer

Next generation sequencing (NGS) technology is an increasingly important clinical tool for therapeutic decision-making. However, interpretation of NGS data presents challenges at the point of care, due to limitations in understanding the clinical importance of gene variants and efficiently translating results into actionable information for the clinician. The present study compared two approaches for annotating and reporting actionable genes and gene mutations from tumor samples: The traditional approach of manual curation, annotation and reporting using an experienced molecular tumor bioinformationist; and a cloud-based cognitive technology, with the goal to detect gene mutations of potential significance in Chinese patients with lung cancer. Data from 285 gene-targeted exon sequencing previously conducted on 115 patient tissue samples between 2014 and 2016 and subsequently manually annotated and evaluated by the Guangdong Lung Cancer Institute (GLCI) research team were analyzed by the Watson for Genomics (WfG) cognitive genomics technology. A comparative analysis of the annotation results of the two methods was conducted to identify quantitative and qualitative differences in the mutations generated. The complete congruence rate of annotation results between WfG analysis and the GLCI bioinformatician was 43.48%. In 65 (56.52%) samples, WfG analysis identified and interpreted, on average, 1.54 more mutation sites in each sample than the manual GLCI review. These mutation sites were located on 27 genes, including EP300, ARID1A, STK11 and DNMT3A. Mutations in the EP300 gene were most prevalent, and present in 30.77% samples. The Tumor Mutation Burden (TMB) interpreted by WfG analysis (1.82) was significantly higher than the TMB (0.73) interpreted by GLCI review. Compared with manual curation by a bioinformatician, WfG analysis provided comprehensive insights and additional genetic alterations to inform clinical therapeutic strategies for patients with lung cancer. These findings suggest the valuable role of cognitive computing to increase efficiency in the comprehensive detection and interpretation of genetic alterations which may inform opportunities for targeted cancer therapies.

Neoantigen prediction in human breast cancer using RNA sequencing data

Neoantigens have attracted attention as biomarkers or therapeutic targets. However, accurate prediction of neoantigens is still challenging, especially in terms of its accuracy and cost. Variant detection using RNA sequencing (RNA-seq) data has been reported to be a low-accuracy but cost-effective tool, but the feasibility of RNA-seq data for neoantigen prediction has not been fully examined. In the present study, we used whole-exome sequencing (WES) and RNA-seq data of tumor and matched normal samples from six breast cancer patients to evaluate the utility of RNA-seq data instead of WES data in variant calling to detect neoantigen candidates. Somatic variants were called in three protocols using: (i) tumor and normal WES data (DNA method, Dm); (ii) tumor and normal RNA-seq data (RNA method, Rm); and (iii) combination of tumor RNA-seq and normal WES data (Combination method, Cm). We found that the Rm had both high false-positive and high false-negative rates because this method depended greatly on the expression status of normal transcripts. When we compared the results of Dm with those of Cm, only 14% of the neoantigen candidates detected in Dm were identified in Cm, but the majority of the missed candidates lacked coverage or variant allele reads in the tumor RNA. In contrast, about 70% of the neoepitope candidates with higher expression and rich mutant transcripts could be detected in Cm. Our results showed that Cm could be an efficient and a cost-effective approach to predict highly expressed neoantigens in tumor samples.

Emerging epigenetic therapeutics for myeloid leukemia: modulating demethylase activity with ascorbate

The past decade has seen a proliferation of drugs that target epigenetic pathways. Many of these drugs were developed to treat acute myeloid leukemia, a condition in which dysregulation of the epigenetic landscape is well established. While these drugs have shown promise, critical issues persist. Specifically, patients with the same mutations respond quite differently to treatment. This is true even with highly specific drugs that are designed to target the underlying oncogenic driver mutations. Furthermore, patients who do respond may eventually develop resistance. There is now evidence that epigenetic heterogeneity contributes, in part, to these issues. Cancer cells also have a remarkable capacity to ‘rewire’ themselves at the epigenetic level in response to drug treatment, and thereby maintain expression of key oncogenes. This epigenetic plasticity is a promising new target for drug development. It is therefore important to consider combination therapy in cases in which both driver mutations and epigenetic plasticity are targeted. Using ascorbate as an example of an emerging epigenetic therapeutic, we review the evidence for its potential use in both of these modes. We provide an overview of 2-oxoglutarate dependent dioxygenases with DNA, histone and RNA demethylase activity, focusing on those which require ascorbate as a cofactor. We also evaluate their role in the development and maintenance of acute myeloid leukemia. Using this information, we highlight situations in which the use of ascorbate to restore 2-oxoglutarate dependent dioxygenase activity could prove beneficial, in contrast to contexts in which targeted inhibition of specific enzymes might be preferred. Finally, we discuss how these insights could be incorporated into the rational design of future clinical trials.

Beyond regulations at DNA levels: A review of epigenetic therapeutics targeting cancer stem cells

In the past few years, the paramount role of cancer stem cells (CSCs), in terms of cancer initiation, proliferation, metastasis, invasion and chemoresistance, has been revealed by accumulating studies. However, this level of cellular plasticity cannot be entirely explained by genetic mutations. Research on epigenetic modifications as a complementary explanation for the properties of CSCs has been increasing over the past several years. Notably, therapeutic strategies are currently being developed in an effort to reverse aberrant epigenetic alterations using specific chemical inhibitors. In this review, we summarize the current understanding of CSCs and their role in cancer progression, and provide an overview of epigenetic alterations seen in CSCs. Importantly, we focus on primary cancer therapies that target the epigenetic modification of CSCs by the use of specific chemical inhibitors, such as histone deacetylase (HDAC) inhibitors, DNA methyltransferase (DNMT) inhibitors and microRNA‐based (miRNA‐based) therapeutics.

Subclonal reconstruction of tumors by using machine learning and population genetics

Most cancer genomic data are generated from bulk samples composed of mixtures of cancer subpopulations, as well as normal cells. Subclonal reconstruction methods based on machine learning aim to separate those subpopulations in a sample and infer their evolutionary history. However, current approaches are entirely data driven and agnostic to evolutionary theory. We demonstrate that systematic errors occur in the analysis if evolution is not accounted for, and this is exacerbated with multi-sampling of the same tumor. We present a novel approach for model-based tumor subclonal reconstruction, called MOBSTER, which combines machine learning with theoretical population genetics. Using public whole-genome sequencing data from 2,606 samples from different cohorts, new data and synthetic validation, we show that this method is more robust and accurate than current techniques in single-sample, multiregion and longitudinal data. This approach minimizes the confounding factors of nonevolutionary methods, thus leading to more accurate recovery of the evolutionary history of human cancers.

HPV-EM: an accurate HPV detection and genotyping EM algorithm

Accurate HPV genotyping is crucial in facilitating epidemiology studies, vaccine trials, and HPV-related cancer research. Contemporary HPV genotyping assays only detect < 25% of all known HPV genotypes and are not accurate for low-risk or mixed HPV genotypes. Current genomic HPV genotyping algorithms use a simple read-alignment and filtering strategy that has difficulty handling repeats and homology sequences. Therefore, we have developed an optimized expectation–maximization algorithm, designated HPV-EM, to address the ambiguities caused by repetitive sequencing reads. HPV-EM achieved 97–100% accuracy when benchmarked using cell line data and TCGA cervical cancer data. We also validated HPV-EM using DNA tiling data on an institutional cervical cancer cohort (96.5% accuracy). Using HPV-EM, we demonstrated HPV genotypic differences in recurrence and patient outcomes in cervical and head and neck cancers.

Mutation of CFAP57, a protein required for the asymmetric targeting of a subset of inner dynein arms in Chlamydomonas, causes primary ciliary dyskinesia

Primary ciliary dyskinesia (PCD) is characterized by chronic airway disease, reduced fertility, and randomization of the left/right body axis. It is caused by defects of motile cilia and sperm flagella. We screened a cohort of affected individuals that lack an obvious axonemal defect for pathogenic variants using whole exome capture, next generation sequencing, and bioinformatic analysis assuming an autosomal recessive trait. We identified one subject with an apparently homozygous nonsense variant [(c.1762C>T), p.(Arg588*)] in the uncharacterized CFAP57 gene. Interestingly, the variant results in the skipping of exon 11 (58 amino acids), which may be due to disruption of an exonic splicing enhancer. In normal human nasal epithelial cells, CFAP57 localizes throughout the ciliary axoneme. Nasal cells from the PCD patient express a shorter, mutant version of CFAP57 and the protein is not incorporated into the axoneme. The missing 58 amino acids include portions of WD repeats that may be important for loading onto the intraflagellar transport (IFT) complexes for transport or docking onto the axoneme. A reduced beat frequency and an alteration in ciliary waveform was observed. Knockdown of CFAP57 in human tracheobronchial epithelial cells (hTECs) recapitulates these findings. Phylogenetic analysis showed that CFAP57 is highly conserved in organisms that assemble motile cilia. CFAP57 is allelic with the BOP2/IDA8/FAP57 gene identified previously in Chlamydomonas reinhardtii. Two independent, insertional fap57 Chlamydomonas mutant strains show reduced swimming velocity and altered waveforms. Tandem mass tag (TMT) mass spectroscopy shows that FAP57 is missing, and the "g" inner dyneins (DHC7 and DHC3) and the "d" inner dynein (DHC2) are reduced, but the FAP57 paralog FBB7 is increased. Together, our data identify a homozygous variant in CFAP57 that causes PCD that is likely due to a defect in the inner dynein arm assembly process.

SomaticCombiner: improving the performance of somatic variant calling based on evaluation tests and a consensus approach

It is challenging to identify somatic variants from high-throughput sequence reads due to tumor heterogeneity, sub-clonality, and sequencing artifacts. In this study, we evaluated the performance of eight primary somatic variant callers and multiple ensemble methods using both real and synthetic whole-genome sequencing, whole-exome sequencing, and deep targeted sequencing datasets with the NA12878 cell line. The test results showed that a simple consensus approach can significantly improve performance even with a limited number of callers and is more robust and stable than machine learning based ensemble approaches. To fully exploit the multi-callers, we also developed a software package, SomaticCombiner, that can combine multiple callers and integrates a new variant allelic frequency (VAF) adaptive majority voting approach, which can maintain sensitive detection for variants with low VAFs.

Applying next-generation sequencing to unravel the mutational landscape in viral quasispecies

Next-generation sequencing (NGS) has revolutionized the scale and depth of biomedical sciences. Because of its unique ability for the detection of sub-clonal variants within genetically diverse populations, NGS has been successfully applied to analyze and quantify the exceptionally-high diversity within viral quasispecies, and many low-frequency drug- or vaccine-resistant mutations of therapeutic importance have been discovered. Although many works have intensively discussed the latest NGS approaches and applications in general, none of them has focused on applying NGS in viral quasispecies studies, mostly due to the limited ability of current NGS technologies to accurately detect and quantify rare viral variants. Here, we summarize several error-correction strategies that have been developed to enhance the detection accuracy of minority variants. We also discuss critical considerations for preparing a sequencing library from viral RNAs and for analyzing NGS data to unravel the mutational landscape.

The clonal evolution of metastatic colorectal cancer

Tumor heterogeneity and evolution drive treatment resistance in metastatic colorectal cancer (mCRC). Patient-derived xenografts (PDXs) can model mCRC biology; however, their ability to accurately mimic human tumor heterogeneity is unclear. Current genomic studies in mCRC have limited scope and lack matched PDXs. Therefore, the landscape of tumor heterogeneity and its impact on the evolution of metastasis and PDXs remain undefined. We performed whole-genome, deep exome, and targeted validation sequencing of multiple primary regions, matched distant metastases, and PDXs from 11 patients with mCRC. We observed intricate clonal heterogeneity and evolution affecting metastasis dissemination and PDX clonal selection. Metastasis formation followed both monoclonal and polyclonal seeding models. In four cases, metastasis-seeding clones were not identified in any primary region, consistent with a metastasis-seeding-metastasis model. PDXs underrepresented the subclonal heterogeneity of parental tumors. These suggest that single sample tumor sequencing and current PDX models may be insufficient to guide precision medicine.