Deciphering intratumor heterogeneity using cancer genome analysis

Cellular and Genetic Background of Osteosarcoma

Osteosarcoma describes a tumor of mesenchymal origin with an annual incidence rate of four to five people per million. Even though chemotherapy treatment has shown success in non-metastatic osteosarcoma, metastatic disease still has a low survival rate of 20%. A targeted therapy approach is limited due to high heterogeneity of tumors, and different underlying mutations. In this review, we will summarize new advances obtained by new technologies, such as next generation sequencing and single-cell sequencing. These new techniques have enabled better assessment of cell populations within osteosarcoma, as well as an understanding of the molecular pathogenesis. We also discuss the presence and properties of osteosarcoma stem cells-the cell population within the tumor that is responsible for metastasis, recurrence, and drug resistance.

Intratumor genetic heterogeneity and clonal evolution to decode endometrial cancer progression

Analyzing different tumor regions by next generation sequencing allows the assessment of intratumor genetic heterogeneity (ITGH), a phenomenon that has been studied widely in some tumor types but has been less well explored in endometrial carcinoma (EC). In this study, we sought to characterize the spatial and temporal heterogeneity of 9 different ECs using whole-exome sequencing, and by performing targeted sequencing validation of the 42 primary tumor regions and 30 metastatic samples analyzed. In addition, copy number alterations of serous carcinomas were assessed by comparative genomic hybridization arrays. From the somatic mutations, identified by whole-exome sequencing, 532 were validated by targeted sequencing. Based on these data, the phylogenetic tree reconstructed for each case allowed us to establish the tumors’ evolution and correlate this to tumor progression, prognosis, and the presence of recurrent disease. Moreover, we studied the genetic landscape of an ambiguous EC and the molecular profile obtained was used to guide the selection of a potential personalized therapy for this patient, which was subsequently validated by preclinical testing in patient-derived xenograft models. Overall, our study reveals the impact of analyzing different tumor regions to decipher the ITGH in ECs, which could help make the best treatment decision.

Single-cell sequencing technology in tumor research

Tumor heterogeneity is a key characteristic of malignant tumors and a significant obstacle in cancer treatment and research. Although bulk tissue sequencing has wide coverage and high accuracy, it can only represent the dominant cell signal information of each sample, while masking the unique gene expression of rare cells; therefore it cannot represent genes that are unstable within a subgroup, but unchanged in a majority of cells. With the progress of genomic technology, the emergence of single-cell sequencing (SCS) has effectively solved the above problem. Genetic, transcriptomic and epigenetic sequencing at the single-cell level provides an important basis for us to correctly classify the cell subsets of heterogeneous tumor populations and to reveal the process of complex changes in tumor cells at the molecular level. Single-cell sequencing technology has been applied to the field of cancer, revealing exciting discoveries in the potential mechanisms of tumor driver gene mutation, clonal evolution, invasion and metastasis. It also provides favorable conditions for developing new tumor biomarkers and providing more accurate and individualized targeted tumor therapy. Herein, we review the steps and methods of single-cell sequencing and highlight the application of SCS in tumor diagnosis and clinical treatment.

Genotoxins exaggerate the stressed state of aneuploid embryonic stem cells via activation of autophagy

The cellular consequences of aneuploidy are largely dependent on the cell types examined. Aneuploid yeasts and mouse embryonic fibroblasts exhibit cell proliferation defects and can be selectively inhibited by compounds that cause proteotoxic or energy stress. By contrast, most aneuploid pluripotent stem cells proliferate rapidly and reach higher saturation densities. The responses of aneuploid pluripotent stem cells to the stress-inducing compounds remain uncharacterized. Here, we tested the response of aneuploid embryonic stem cells to several compounds that caused proteotoxic, energy and genotoxic stress using previously established mouse embryonic stem cell lines trisomic for chromosome 6, 8, 11, or 15. Not all trisomic embryonic stem cells were selectively inhibited by compounds that cause proteotoxic or energy stress. However, most of these cells exhibited increased sensitivity to genotoxins. They displayed elevated DNA damage response as characterized by increased γH2A.X foci under genotoxic stress. Further investigations indicated that elevated autophagy levels might contribute to the increased cytotoxic effects of genotoxins on trisomic embryonic stem cells. Our study laid the foundation for eliminating aneuploidy that might be an effective approach for controlling cancer progression.

Clonal selection confers distinct evolutionary trajectories in BRAF-driven cancers

Molecular determinants governing the evolution of tumor subclones toward phylogenetic branches or fixation remain unknown. Using sequencing data, we model the propagation and selection of clones expressing distinct categories of BRAF mutations to estimate their evolutionary trajectories. We show that strongly activating BRAF mutations demonstrate hard sweep dynamics, whereas mutations with less pronounced activation of the BRAF signaling pathway confer soft sweeps or are subclonal. We use clonal reconstructions to estimate the strength of "driver" selection in individual tumors. Using tumors cells and human-derived murine xenografts, we show that tumor sweep dynamics can significantly affect responses to targeted inhibitors of BRAF/MEK or DNA damaging agents. Our study uncovers patterns of distinct BRAF clonal evolutionary dynamics and nominates therapeutic strategies based on the identity of the BRAF mutation and its clonal composition.

BRCA Genes: The Role in Genome Stability, Cancer Stemness and Therapy Resistance

Carcinogenesis is a multistep process, and tumors frequently harbor multiple mutations regulating genome integrity, cell division and death. The integrity of cellular genome is closely controlled by the mechanisms of DNA damage signaling and DNA repair. The association of breast cancer susceptibility genes BRCA1 and BRCA2 with breast and ovarian cancer development was first demonstrated over 20 years ago. Since then the germline mutations within these genes were linked to genomic instability and increased risk of many other cancer types. Genomic instability is an engine of the oncogenic transformation of non-tumorigenic cells into tumor-initiating cells and further tumor evolution. In this review we discuss the biological functions of BRCA1 and BRCA2 genes and the role of BRCA mutations in tumor initiation, regulation of cancer stemness, therapy resistance and tumor progression.

Intratumor heterogeneity inferred from targeted deep sequencing as a prognostic indicator

Tumor genetic heterogeneity may underlie poor clinical outcomes because diverse subclones could be comprised of metastatic and drug resistant cells. Targeted deep sequencing has been used widely as a diagnostic tool to identify actionable mutations in cancer patients. In this study, we evaluated the clinical utility of estimating tumor heterogeneity using targeted panel sequencing data. We investigated the prognostic impact of a tumor heterogeneity (TH) index on clinical outcomes, using mutational profiles from targeted deep sequencing data acquired from 1,352 patients across 8 cancer types. The TH index tended to be increased in high pathological stage disease in several cancer types, indicating clonal expansion of cancer cells as tumor progression proceeds. In colorectal cancer patients, TH index values also correlated significantly with clinical prognosis. Integration of the TH index with genomic and clinical features could improve the power of risk prediction for clinical outcomes. In conclusion, deep sequencing to determine the TH index could serve as a promising prognostic indicator in cancer patients.

The sensitivity and specificity of serum glycan-based biomarkers for cancer detection

Most of clinically used serum biomarkers for cancer detection were established in early 1980s when the Nobel Prize in physiology or medicine was awarded for the "discovery of the principle for the production of monoclonal antibodies." Using this "Nobel" technology, various monoclonal antibodies were obtained when different types of cancer cells were injected into mice and the ligands on the cancer cell surface were characterized. Both aberrant glycan structures and aberrant glycan-associated glycoproteins were revealed as a common feature of cancer cell surfaces through the specific interactions with the monoclonal antibodies. These results indicate that the biosynthesis of the environment-sensitive glycan structures goes awry in cancer cells, which is beyond genetic mutations. Later on, the glycan-related biomarkers were detected in the sera of cancer patients and then developed into serum biomarkers, such as CA125, CA153, CA195, CA199, CA242, CA27.29, CA50, and CA724, which are still in clinical use as of today. During the past 30 years, even with the advancement of different OMICS technologies not limited to genomics, epigenomics, proteomics, glycomics, lipidomics, and metabolomics, very few serum biomarkers have been introduced into clinical practice. The reason is that most of the newly discovered cancer biomarkers are inferior in terms of sensitivity and specificity to these biomarkers. We will summarize the reported sensitivity and specificity of currently used cancer biomarkers, especially the glycan-related biomarkers, in the forms of tables and radar plots and discuss the pros and cons of currently used cancer biomarkers.

Performance validation of an amplicon-based targeted next-generation sequencing assay and mutation profiling of 648 Chinese colorectal cancer patients

Next-generation sequencing (NGS) has become a promising approach for tumor somatic mutation detection. However, stringent validation is required for its application on clinical specimens, especially for low-quality formalin-fixed paraffin-embedded (FFPE) tissues. Here, we validated the performance of an amplicon-based targeted NGS assay, OncoAim™ DNA panel, on both commercial reference FFPE samples and clinical FFPE samples of Chinese colorectal cancer (CRC) patients. Then we profiled the mutation spectrum of 648 Chinese CRC patients in a multicenter study to explore its clinical utility. This NGS assay achieved 100% test specificity and 95-100% test sensitivity for variants with mutant allele frequency (MAF) ≥ 5% when median read depth ≥ 500×. The orthogonal methods including amplification refractory mutation system (ARMS)-PCR and Sanger sequencing validated that NGS generated three false negatives (FNs) but no false positives (FPs) among 516 clinical samples for KRAS aberration detection. Genomic profiling of Chinese CRC patients with this assay revealed that 63.3% of the tumors harbored clinically actionable alterations. Besides the commonly mutated genes including TP53 (52.82%), KRAS (46.68%), APC (24.09%), PIK3CA (18.94%), SMAD4 (9.47%), BRAF (6.15%), FBXW7 (5.32%), and NRAS (4.15%), other less frequently mutated genes were also identified. Statistically significant association of specific mutated genes with certain clinicopathological features was detected, e.g., both BRAF and PIK3CA were more prevalent in right-side CRC (p < 0.001 and p = 0.002, respectively). We concluded this targeted NGS assay is qualified for clinical practice, and our findings could help the diagnosis and prognosis of Chinese CRC patients.

Sample-Index Misassignment Impacts Tumour Exome Sequencing

Sample pooling enabled by dedicated indexes is a common strategy for cost-effective and robust high-throughput sequencing. Index misassignment leading to mutual contamination between pooled samples has however been described as a general problem of the latest Illumina sequencing instruments utilizing exclusion amplification. Using real-life data from multiple tumour sequencing projects, we demonstrate that index misassignment can induce artefactual variant calls closely resembling true, high-quality somatic variants. These artefactual calls potentially impact cancer applications utilizing low allelic frequencies, such as in clonal analysis of tumours. We discuss the available countermeasures with an emphasis on improved library indexing methods, and provide software that can assist in the identification of variants that may be consequences of index misassignment.

Novel putative drivers revealed by targeted exome sequencing of advanced solid tumors

Next generation sequencing (NGS) is becoming increasingly integrated into oncological practice and clinical research. NGS methods have also provided evidence for clonal evolution of cancers during disease progression and treatment. The number of variants associated with response to specific therapeutic agents keeps increasing. However, the identification of novel driver mutations as opposed to passenger (phenotypically silent or clinically irrelevant) mutations remains a major challenge. We conducted targeted exome sequencing of advanced solid tumors from 44 pre-treated patients with solid tumors including breast, colorectal and lung carcinomas, neuroendocrine tumors, sarcomas and others. We catalogued established driver mutations and putative new drivers as predicted by two distinct algorithms. The established drivers we detected were consistent with published observations. However, we also detected a significant number of mutations with driver potential never described before in each tumor type we studied. These putative drivers belong to key cell fate regulatory networks, including potentially druggable pathways. Should our observations be confirmed, they would support the hypothesis that new driver mutations are selected by treatment in clinically aggressive tumors, and indicate a need for longitudinal genomic testing of solid tumors to inform second line cancer treatment.

Somatic evolutionary timings of driver mutations

BACKGROUND

A unified analysis of DNA sequences from hundreds of tumors concluded that the driver mutations primarily occur in the earliest stages of cancer formation, with relatively few driver mutation events detected in the late-arising subclones. However, emerging evidence from the sequencing of multiple tumors and tumor regions per individual suggests that late-arising subclones with additional driver mutations are underestimated in single-sample analyses.

METHODS

To test whether driver mutations generally map to early tumor development, we examined multi-regional tumor sequencing data from 101 individuals reported in 11 published studies. Following previous studies, we annotated mutations as early-arising when all tumors/regions had those mutations (ubiquitous). We then inferred the fraction of mutations occurring early and compared it with late-arising mutations that were found in only single tumors/regions.

RESULTS

While a large fraction of driver mutations in tumors occurred relatively early in cancers, later driver mutations occurred at least as frequently as the early drivers in a substantial number of patients. This result was robust to many different approaches to annotate driver mutations. The relative frequency of early and late driver mutations varied among patients of the same cancer type and in different cancer types. We found that previous reports of the preponderance of early driver mutations were primarily informed by analysis of single tumor variant allele profiles, with which it is challenging to clearly distinguish between early and late drivers.

CONCLUSIONS

The origin and preponderance of new driver mutations are not limited to early stages of tumor evolution, with different tumors and regions showing distinct driver mutations and, consequently, distinct characteristics. Therefore, tumors with extensive intratumor heterogeneity appear to have many newly acquired drivers.

The diverse roles of glutathione-associated cell resistance against hypericin photodynamic therapy

The diverse responses of different cancers to treatments such as photodynamic therapy of cancer (PDT) have fueled a growing need for reliable predictive markers for treatment outcome. In the present work we have studied the differential response of two phenotypically and genotypically different breast adenocarcinoma cell lines, MCF7 and MDA-MB-231, to hypericin PDT (HYP-PDT). MDA-MB-231 cells were 70% more sensitive to HYP PDT than MCF7 cells at LD50. MCF7 were found to express a substantially higher level of glutathione peroxidase (GPX4) than MDA-MB-231, while MDA-MB-231 differentially expressed glutathione-S-transferase (GSTP1), mainly used for xenobiotic detoxification. Eighty % reduction of intracellular glutathione (GSH) by buthionine sulfoximine (BSO), largely enhanced the sensitivity of the GSTP1 expressing MDA-MB-231 cells to HYP-PDT, but not in MCF7 cells. Further inhibition of the GSH reduction however by carmustine (BCNU) resulted in an enhanced sensitivity of MCF7 to HYP-PDT. HYP loading studies suggested that HYP can be a substrate of GSTP for GSH conjugation as BSO enhanced the cellular HYP accumulation by 20% in MDA-MB-231 cells, but not in MCF7 cells. Studies in solutions showed that L-cysteine can bind the GSTP substrate CDNB in the absence of GSTP. This means that the GSTP-lacking MCF7 may use L-cysteine for xenobiotic detoxification, especially during GSH synthesis inhibition, which leads to L-cysteine build-up. This was confirmed by the lowered accumulation of HYP in both cell lines in the presence of BSO and the L-cysteine source NAC. NAC reduced the sensitivity of MCF7, but not MDA-MB-231, cells to HYP PDT which is in accordance with the antioxidant effects of L-cysteine and its potential as a GSTP substrate. As a conclusion we have herein shown that the different GSH based cell defense mechanisms can be utilized as predictive markers for the outcome of PDT and as a guide for selecting optimal combination strategies.

Towards personalized tumor markers

The cancer biomarker discovery pipeline is progressing slowly. The difficulties of finding novel and effective biomarkers for diagnosis and management of cancer patients are well-known. We speculate that it is unlikely to discover new serological biomarkers characterized by high sensitivity and specificity. This projection is supported by recent findings that cancers are genetically highly heterogeneous. Here, we propose a new way of improving the landscape of cancer biomarker research. There are currently hundreds, if not thousands, of described biomarkers which perform at high specificity (> 90%), but at relatively low sensitivity (< 30%). We call these “rare tumor markers.” Borrowing from the principles of precision medicine, we advocate that among these low sensitivity markers, some may be useful to specific patients. We suggest screening new patients for hundreds to thousands of cancer biomarkers to identify a few that are informative, and then use them clinically. This is similar to what we currently do with genomics to identify personalized therapies. We further suggest that this approach may explain as to why some biomarkers are elevated in only a small group of patients. It is likely that these differences in expression are linked to specific genomic alterations, which could then be found with genomic sequencing.

Whole tumor RNA-sequencing and deconvolution reveal a clinically-prognostic PTEN/PI3K-regulated glioma transcriptional signature

The concept that solid tumors are maintained by a productive interplay between neoplastic and non-neoplastic elements has gained traction with the demonstration that stromal fibroblasts and immune system cells dictate cancer development and progression. While less studied, brain tumor (glioma) biology is likewise influenced by non-neoplastic immune system cells (macrophages and microglia) which interact with neoplastic glioma cells to create a unique physiological state (glioma ecosystem) distinct from that found in the normal tissue. To explore this neoplastic ground state, we leveraged several preclinical mouse models of neurofibromatosis type 1 (NF1) optic glioma, a low-grade astrocytoma whose formation and maintenance requires productive interactions between non-neoplastic and neoplastic cells, and employed whole tumor RNA-sequencing and mathematical deconvolution strategies to characterize this low-grade glioma ecosystem as an aggregate of cellular and acellular elements. Using this approach, we demonstrate that optic gliomas generated by altering the germline Nf1 gene mutation, the glioma cell of origin, or the presence of co-existing genetic alterations represent molecularly-distinct tumors. However, these optic glioma tumors share a 25-gene core signature, not found in normal optic nerve, that is normalized by microglia inhibition (minocycline), but not conventional (carboplatin) or molecularly-targeted (rapamycin) chemotherapy. Lastly, we identify a genetic signature conferred by Pten reduction and corrected by PI3K inhibition. This signature predicts progression-free survival in patients with either low-grade or high-grade glioma. Collectively, these findings support the concept that gliomas are composite ecological systems whose biology and response to therapy may be best defined by examining the tumor as a whole.

Cancer stem cells in osteosarcoma

Osteosarcoma is the most common primary bone tumour in children and adolescents and advanced osteosarcoma patients with evidence of metastasis share a poor prognosis. Osteosarcoma frequently gains resistance to standard therapies highlighting the need for improved treatment regimens and identification of novel therapeutic targets. Cancer stem cells (CSC) represent a sub-type of tumour cells attributed to critical steps in cancer including tumour propagation, therapy resistance, recurrence and in some cases metastasis. Recent published work demonstrates evidence of cancer stem cell phenotypes in osteosarcoma with links to drug resistance and tumorigenesis. In this review we will discuss the commonly used isolation techniques for cancer stem cells in osteosarcoma as well as the identified biochemical and molecular markers.

Rnd3 in Cancer: A Review of the Evidence for Tumor Promoter or Suppressor

Rho-GTPases are members of the Ras superfamily of small GTPases and are general modulators of important cellular processes in tumor biology such as migration and proliferation. Among these proteins, Rnd3/RhoE, an atypical Rho-GTPase devoid of GTP hydrolytic activity, has recently been studied for its putative role in tumorigenesis. Indeed, Rnd3 is implicated in processes, such as proliferation and migration, whose deregulation is linked to cancer development and metastasis. The aim of this review is to provide an overview of the data surrounding Rnd3 deregulation in cancers, its origin, and consequences. Presented here is a comprehensive account of the expression status and biological output obtained in prostate, liver, stomach, colon, lung, and brain cancers as well as in melanoma and squamous cell carcinoma. Although there appears to be no general consensus about Rnd3 expression in cancers as this protein is differently altered according to the tumor context, these alterations overwhelmingly favor a protumorigenic role. Thus, depending on the tumor type, it may behave either as a tumor suppressor or as a tumor promoter. Importantly, the deregulation of Rnd3, in most cases, is linked to patient poor outcome.

IMPLICATIONS

Rnd3 has prognostic marker potential as exemplified in lung cancers and Rnd3 or Rnd3-associated signaling pathways may represent a new putative therapeutic target. Mol Cancer Res; 14(11); 1033-44. ©2016 AACR.