Inferring mutational timing and reconstructing tumour evolutionary histories

A war on many fronts: cross disciplinary approaches for novel cancer treatment strategies

Cancer is a disease characterized by uncontrolled cellular growth where cancer cells take advantage of surrounding cellular populations to obtain resources and promote invasion. Carcinomas are the most common type of cancer accounting for almost 90% of cancer cases. One of the major subtypes of carcinomas are adenocarcinomas, which originate from glandular cells that line certain internal organs. Cancers such as breast, prostate, lung, pancreas, colon, esophageal, kidney are often adenocarcinomas. Current treatment strategies include surgery, chemotherapy, radiation, targeted therapy, and more recently immunotherapy. However, patients with adenocarcinomas often develop resistance or recur after the first line of treatment. Understanding how networks of tumor cells interact with each other and the tumor microenvironment is crucial to avoid recurrence, resistance, and high-dose therapy toxicities. In this review, we explore how mathematical modeling tools from different disciplines can aid in the development of effective and personalized cancer treatment strategies. Here, we describe how concepts from the disciplines of ecology and evolution, economics, and control engineering have been applied to mathematically model cancer dynamics and enhance treatment strategies.

Evolutionary mode and timing of dissemination of high-grade serous carcinomas

Dissemination within the peritoneal cavity is a main determinant of poor patient outcomes from high-grade serous carcinomas (HGSCs). The dissemination process is poorly understood from a cancer evolutionary perspective. We reconstructed the evolutionary trajectories across a median of 5 tumor sites and regions from each of 23 patients based on deep whole-exome sequencing. Polyclonal cancer origin was detected in 1 patient. Ovarian tumors had more complex subclonal architectures than other intraperitoneal tumors in each patient, which indicated that tumors developed earlier in the ovaries. Three common modes of dissemination were identified, including monoclonal or polyclonal dissemination of monophyletic (linear) or polyphyletic (branched) subclones. Mutation profiles of initial or disseminated clones varied greatly among cancers, but recurrent mutations were found in 7 cancer-critical genes, including TP53, BRCA1, BRCA2, and DNMT3A, and in the PI3K/AKT1 pathway. Disseminated clones developed late in the evolutionary trajectory models of most cancers, in particular in cancers with DNA damage repair deficiency. Polyclonal dissemination was predicted to occur predominantly as a single and rapid wave, but chemotherapy exposure was associated with higher genomic diversity of disseminated clones. In conclusion, we described three common evolutionary dissemination modes across HGSCs and proposed factors associated with dissemination diversity.

Mutational landscape of paired primary and synchronous metastatic lymph node in chemotherapy naive gallbladder cancer

BACKGROUND

Comprehensive genomic analysis of paired primary tumors and their metastatic lesions may provide new insights into the biology of metastatic processes and therefore guide the development of novel strategies for intervention. To date, our knowledge of the genetic divergence and phylogenetic relationships in gallbladder cancer (GBC) is limited.

METHODS

We performed whole exome sequencing for 5 patients with primary tumor, metastatic lymph node (LNM) and corresponding normal tissue. Mutations, mutation signatures and copy number variations were analyzed with state-of-art bioinformatics methods. Phylogenetic tree was also generated to infer metastatic pattern.

RESULTS

Five driver mutations were detected in these patients. Among which, TP53 was the only shared mutation between primary tumor and LNM. Although tumor mutational burden was comparable between primary tumor and LNM, higher mutation burden was observed in LNM of one patient. Copy number variations (CNVs) burden was higher in LNM than their primary tumor. Phylogenetic analysis indicated both linear and parallel progression of metastasis exist in these patients. TP53 mutation and CNVs were homogenously between primary tumor and LNM.

CONCLUSIONS

High consistence of genetic landscape were shown between primary tumor and LNM in GBC. However, heterogenicity still exist between primary tumor and LNM in particular patients in term of driver mutation, TMB and CNV burden. Phylogenetic analysis indicated both Linear and parallel progression of metastasis were exist among these patients.

CloneSig can jointly infer intra-tumor heterogeneity and mutational signature activity in bulk tumor sequencing data

Systematic DNA sequencing of cancer samples has highlighted the importance of two aspects of cancer genomics: intra-tumor heterogeneity (ITH) and mutational processes. These two aspects may not always be independent, as different mutational processes could be involved in different stages or regions of the tumor, but existing computational approaches to study them largely ignore this potential dependency. Here, we present CloneSig, a computational method to jointly infer ITH and mutational processes in a tumor from bulk-sequencing data. Extensive simulations show that CloneSig outperforms current methods for ITH inference and detection of mutational processes when the distribution of mutational signatures changes between clones. Applied to a large cohort of 8,951 tumors with whole-exome sequencing data from The Cancer Genome Atlas, and on a pan-cancer dataset of 2,632 whole-genome sequencing tumor samples from the Pan-Cancer Analysis of Whole Genomes initiative, CloneSig obtains results overall coherent with previous studies.

Mathematical modeling of multiple pathways in colorectal carcinogenesis using dynamical systems with Kronecker structure

Like many other types of cancer, colorectal cancer (CRC) develops through multiple pathways of carcinogenesis. This is also true for colorectal carcinogenesis in Lynch syndrome (LS), the most common inherited CRC syndrome. However, a comprehensive understanding of the distribution of these pathways of carcinogenesis, which allows for tailored clinical treatment and even prevention, is still lacking. We suggest a linear dynamical system modeling the evolution of different pathways of colorectal carcinogenesis based on the involved driver mutations. The model consists of different components accounting for independent and dependent mutational processes. We define the driver gene mutation graphs and combine them using the Cartesian graph product. This leads to matrix components built by the Kronecker sum and product of the adjacency matrices of the gene mutation graphs enabling a thorough mathematical analysis and medical interpretation. Using the Kronecker structure, we developed a mathematical model which we applied exemplarily to the three pathways of colorectal carcinogenesis in LS. Beside a pathogenic germline variant in one of the DNA mismatch repair (MMR) genes, driver mutations in APC, CTNNB1, KRAS and TP53 are considered. We exemplarily incorporate mutational dependencies, such as increased point mutation rates after MMR deficiency, and based on recent experimental data, biallelic somatic CTNNB1 mutations as common drivers of LS-associated CRCs. With the model and parameter choice, we obtained simulation results that are in concordance with clinical observations. These include the evolution of MMR-deficient crypts as early precursors in LS carcinogenesis and the influence of variants in MMR genes thereon. The proportions of MMR-deficient and MMR-proficient APC-inactivated crypts as first measure for the distribution among the pathways in LS-associated colorectal carcinogenesis are compatible with clinical observations. The approach provides a modular framework for modeling multiple pathways of carcinogenesis yielding promising results in concordance with clinical observations in LS CRCs.

Role of Alpelisib in the Treatment of PIK3CA-Mutated Breast Cancer: Patient Selection and Clinical Perspectives

The PI3K/AKT/mTOR pathway has long been known to play a major role in the growth and survival of cancer cells. Breast tumors often harbor PIK3CA gene alterations, which therefore constitute a rational drug target. However, it has taken many years to demonstrate clinically-relevant efficacy of PI3K inhibition and eventually attain regulatory approvals. As data on PI3K inhibitors continue to mature, this review updates and summarizes the current state of the science, including the prognostic role of PIK3CA alterations in breast cancer; the evolution of PI3K inhibitors; the clinical utility of the first-in-class oral selective PI3Kα inhibitor, alpelisib; PIK3CA mutation detection techniques; and adverse effect management. PIK3CA-mutated breast carcinomas predict survival benefit from PI3K inhibitor therapy. The pan-PI3K inhibitor, buparlisib and the beta-isoform-sparing PI3K inhibitor, taselisib, met efficacy endpoints in clinical trials, but pictilisib did not; moreover, poor tolerability of these three drugs abrogated further clinical trials. Alpelisib is better tolerated, with a more manageable toxicity profile; the principal adverse events, hyperglycemia, rash and diarrhea, can be mitigated by intensive monitoring and timely intervention, thereby enabling patients to remain adherent to clinically beneficial treatment. Alpelisib plus endocrine therapy shows promising efficacy for treating postmenopausal women with HR+/HER2- advanced breast cancer. Available evidence supporting using alpelisib after disease progression on first-line endocrine therapy with or without CDK4/6 inhibitors justifies PIK3CA mutation testing upon diagnosing HR+/HER2- advanced breast cancer, which can be done using either tumor tissue or circulating tumor DNA. With appropriate toxicity management and patient selection using validated testing methods, all eligible patients can potentially benefit from this new treatment. Further clinical trials to assess combinations of hormone therapy with PI3K, AKT, mTOR, or CDK 4/6 inhibitors, or studies in men and women with other breast subtypes are ongoing.

Molecular Pap Smear: Validation of HPV Genotype and Host Methylation Profiles of ADCY8, CDH8, and ZNF582 as a Predictor of Cervical Cytopathology

Primary high-risk Human Papillomavirus (hrHPV) screening has recently become an accepted standalone or co-test with conventional cytology. Unfortunately, hrHPV singularly lacks specificity for cytopathological grade. However, mechanisms and markers of evolving virus-host interactions at the epigenome level may be harnessed as a better predictor of carcinogenesis. This study aimed to validate and expand the clinical performance of a multiparametric biomarker panel, referred to as the "Molecular Pap smear" based, on HPV genotype and ADCY8, CDH8 and ZNF582 CpG-methylation as a predictive classifier of cervical cytology. This prospective, cross-sectional study used an independent cohort of residual liquid-based cytology for HPV genotyping and epigenetic analysis. Extracted DNA underwent parallel PCR using 3 primer sets for HPV DNA amplification. HPV-infected samples were genotyped by Sanger sequencing. Promoter methylation levels of 3 tumor suppressor genes were quantified by bisulfite-pyrosequencing of genomic DNA on the newest high-resolution PyroMark Q48 platform. Logistic model performance was compared, and model parameters were used to predict and classify binary cytological outcomes. A total of 883 samples were analyzed. HPV DNA positivity correlated with worsening grade: 125/237 (53%) NILM; 136/235 (58%) ASCUS; 222/229 (97%) LSIL; and 157/182 (86%) HSIL samples. The proportion of carcinogenic HPV-types in PCR-positive sequenceable samples correlated with worsening grade: NILM 34/98 (35%); ASCUS 50/113 (44%); LSIL 92/214 (43%); HSIL 129/152 (85%). Additionally, ADCY8, CDH8, and ZNF582 methylation levels increased in direct correlation with worsening grade. Overall, the multi-marker modeling parameters predicted binarized cytological outcomes better than HPV-type alone with significantly higher area under the receiver operator curve (AUC)s, respectively: NILM vs. > NILM (AUC 0.728 vs. 0.709); NILM/ASCUS vs. LSIL/HSIL (AUC 0.805 vs. 0.776); and Collapse

Key Words

• DNA methylation
• carcinogenesis
• epigenetic modification
• evolution
• host-pathogen interactions
• human papillomavirus infection
• pap smear
• pyrosequencing

Collapse

MESH Headings Collapse

Grants Collapse

Affiliation(s)

Jane Shen-Gunther Gynecologic Oncology & Clinical Investigation, Department of Clinical Investigation, Brooke Army Medical Center, Fort Sam Houston, TX, United States Department of Molecular Medicine, Cancer Therapy and Research Center, University of Texas Health Science Center at San Antonio, San Antonio, TX, United States
Qingqing Xia Department of Clinical Investigation, Brooke Army Medical Center, Fort Sam Houston, TX, United States
Winfred Stacey Department of Clinical Investigation, Brooke Army Medical Center, Fort Sam Houston, TX, United States
Heisy B. Asusta Department of Obstetrics and Gynecology, Brooke Army Medical Center, Fort Sam Houston, TX, United States

Collapse

Chromosome Instability; Implications in Cancer Development, Progression, and Clinical Outcomes

Chromosome instability (CIN) refers to an ongoing rate of chromosomal changes and is a driver of genetic, cell-to-cell heterogeneity. It is an aberrant phenotype that is intimately associated with cancer development and progression. The presence, extent, and level of CIN has tremendous implications for the clinical management and outcomes of those living with cancer. Despite its relevance in cancer, there is still extensive misuse of the term CIN, and this has adversely impacted our ability to identify and characterize the molecular determinants of CIN. Though several decades of genetic research have provided insight into CIN, the molecular determinants remain largely unknown, which severely limits its clinical potential. In this review, we provide a definition of CIN, describe the two main types, and discuss how it differs from aneuploidy. We subsequently detail its impact on cancer development and progression, and describe how it influences metastatic potential with reference to cancer prognosis and outcomes. Finally, we end with a discussion of how CIN induces genetic heterogeneity to influence the use and efficacy of several precision medicine strategies, including patient and risk stratification, as well as its impact on the acquisition of drug resistance and disease recurrence.

Effects of ordered mutations on dynamics in signaling networks

BACKGROUND

Many previous clinical studies have found that accumulated sequential mutations are statistically related to tumorigenesis. However, they are limited in fully elucidating the significance of the ordered-mutation because they did not focus on the network dynamics. Therefore, there is a pressing need to investigate the dynamics characteristics induced by ordered-mutations.

METHODS

To quantify the ordered-mutation-inducing dynamics, we defined the mutation-sensitivity and the order-specificity that represent if the network is sensitive against a double knockout mutation and if mutation-sensitivity is specific to the mutation order, respectively, using a Boolean network model.

RESULTS

Through intensive investigations, we found that a signaling network is more sensitive when a double-mutation occurs in the direction order inducing a longer path and a smaller number of paths than in the reverse order. In addition, feedback loops involving a gene pair decreased both the mutation-sensitivity and the order-specificity. Next, we investigated relationships of functionally important genes with ordered-mutation-inducing dynamics. The network is more sensitive to mutations subject to drug-targets, whereas it is less specific to the mutation order. Both the sensitivity and specificity are increased when different-drug-targeted genes are mutated. Further, we found that tumor suppressors can efficiently suppress the amplification of oncogenes when the former are mutated earlier than the latter.

CONCLUSION

Taken together, our results help to understand the importance of the order of mutations with respect to the dynamical effects in complex biological systems.

Molecular Biology and Evolution of Cancer: From Discovery to Action

Cancer progression is an evolutionary process. During this process, evolving cancer cell populations encounter restrictive ecological niches within the body, such as the primary tumor, circulatory system, and diverse metastatic sites. Efforts to prevent or delay cancer evolution-and progression-require a deep understanding of the underlying molecular evolutionary processes. Herein we discuss a suite of concepts and tools from evolutionary and ecological theory that can inform cancer biology in new and meaningful ways. We also highlight current challenges to applying these concepts, and propose ways in which incorporating these concepts could identify new therapeutic modes and vulnerabilities in cancer.

Colon cancer stemness as a reversible epigenetic state: Implications for anticancer therapies

The recent discovery of cancer cell plasticity, i.e. their ability to reprogram into cancer stem cells (CSCs) either naturally or under chemotherapy and/or radiotherapy, has changed, once again, the way we consider cancer treatment. If cancer stemness is a reversible epigenetic state rather than a genetic identity, opportunities will arise for therapeutic strategies that remodel epigenetic landscapes of CSCs. However, the systematic use of DNA methyltransferase and histone deacetylase inhibitors, alone or in combination, in advanced solid tumors including colorectal cancers, regardless of their molecular subtypes, does not seem to be the best strategy. In this review, we first summarize the knowledge researchers have gathered on the epigenetic signatures of CSCs with the difficulty of isolating rare populations of cells. We raise questions about the relevant use of currently available epigenetic inhibitors (epidrugs) while the expression of numerous cancer stem cell markers are often repressed by epigenetic mechanisms. These markers include the three cluster of differentiation CD133, CD44 and CD166 that have been extensively used for the isolation of colon CSCs.and . Finally, we describe current treatment strategies using epidrugs, and we hypothesize that, using correlation tools comparing associations of relevant CSC markers with chromatin modifier expression, we could identify better candidates for epienzyme targeting.

Predicting clone genotypes from tumor bulk sequencing of multiple samples

Motivation

Analyses of data generated from bulk sequencing of tumors have revealed extensive genomic heterogeneity within patients. Many computational methods have been developed to enable the inference of genotypes of tumor cell populations (clones) from bulk sequencing data. However, the relative and absolute accuracy of available computational methods in estimating clone counts and clone genotypes is not yet known.

Results

We have assessed the performance of nine methods, including eight previously-published and one new method (CloneFinder), by analyzing computer simulated datasets. CloneFinder, LICHeE, CITUP and cloneHD inferred clone genotypes with low error (<5% per clone) for a majority of datasets in which the tumor samples contained evolutionarily-related clones. Computational methods did not perform well for datasets in which tumor samples contained mixtures of clones from different clonal lineages. Generally, the number of clones was underestimated by cloneHD and overestimated by PhyloWGS, and BayClone2, Canopy and Clomial required prior information regarding the number of clones. AncesTree and Canopy did not produce results for a large number of datasets. Overall, the deconvolution of clone genotypes from single nucleotide variant (SNV) frequency differences among tumor samples remains challenging, so there is a need to develop more accurate computational methods and robust software for clone genotype inference.

Availability and implementation

CloneFinder is implemented in Python and is available from https://github.com/gstecher/CloneFinderAPI.

Supplementary information

Supplementary data are available at Bioinformatics online.

Cell-lineage level-targeted sequencing to identify acute myeloid leukemia with myelodysplasia-related changes

Acute myeloid leukemia (AML) is a clonal myeloid neoplasm that typically arises de novo; however, some cases evolve from a preleukemic state, such as myelodysplastic syndrome (MDS). Such secondary AMLs and those with typical MDS-related clinical features are known as AMLs with myelodysplasia-related changes (AML-MRC). Because patients with AML-MRC have poor prognosis, more accurate diagnostic approaches are required. In this study, we performed targeted sequencing of 54 genes in 3 cell populations (granulocyte, blast, and T-cell fractions) using samples from 13 patients with MDS, 16 patients with clinically diagnosed AML-MRC, 4 patients with suspected AML-MRC but clinically diagnosed as AML not otherwise specified (AML-NOS), and 11 patients with de novo AML. We found that overlapping mutations, defined as those shared at least by the blast and granulocyte fractions, were significantly enriched in patients with MDS and AML-MRC, including those with suspected AML-MRC, indicating a substantial history of clonal hematopoiesis. In contrast, blast-specific nonoverlapping mutations were significantly enriched in patients with de novo AML. Furthermore, the presence of overlapping mutations, excluding DNMT3A, TET2, and ASXL1, effectively segregated patients with MDS and AML-MRC or suspected AML-MRC from patients with de novo AML. Additionally, the presence of ≥3 mutations in the blast fraction was useful for distinguishing patients with AML-MRC from those with MDS. In conclusion, our approach is useful for classifying clinically diagnosable AML-MRC and identifying clinically diagnosed AML-NOS as latent AML-MRC. Additional prospective studies are needed to confirm the utility of this approach.

Modeling breast cancer progression to bone: how driver mutation order and metabolism matter

BACKGROUND

Not all the mutations are equally important for the development of metastasis. What about their order? The survival of cancer cells from the primary tumour site to the secondary seeding sites depends on the occurrence of very few driver mutations promoting oncogenic cell behaviours. Usually these driver mutations are among the most effective clinically actionable target markers. The quantitative evaluation of the effects of a mutation across primary and secondary sites is an important challenging problem that can lead to better predictability of cancer progression trajectory.

RESULTS

We introduce a quantitative model in the framework of Cellular Automata to investigate the effects of metabolic mutations and mutation order on cancer stemness and tumour cell migration from breast, blood to bone metastasised sites. Our approach models three types of mutations: driver, the order of which is relevant for the dynamics, metabolic which support cancer growth and are estimated from existing databases, and non-driver mutations. We integrate the model with bioinformatics analysis on a cancer mutation database that shows metabolism-modifying alterations constitute an important class of key cancer mutations.

CONCLUSIONS

Our work provides a quantitative basis of how the order of driver mutations and the number of mutations altering metabolic processis matter for different cancer clones through their progression in breast, blood and bone compartments. This work is innovative because of multi compartment analysis and could impact proliferation of therapy-resistant clonal populations and patient survival. Mathematical modelling of the order of mutations is presented in terms of operators in an accessible way to the broad community of researchers in cancer models so to inspire further developments of this useful (and underused in biomedical models) methodology. We believe our results and the theoretical framework could also suggest experiments to measure the overall personalised cancer mutational signature.

Timing somatic events in the evolution of cancer

Cancer arises through the accumulation of somatic mutations over time. An understanding of the sequence of events during this process should allow both earlier diagnosis and better prediction of cancer progression. However, the pathways of tumor evolution have not yet been comprehensively characterized. With the advent of whole genome sequencing, it is now possible to infer the evolutionary history of single tumors from the snapshot of their genome taken at diagnosis, giving new insights into the biology of tumorigenesis.

Cancer-mutation network and the number and specificity of driver mutations

Cancer genomics has produced extensive information on cancer-associated genes, but the number and specificity of cancer-driver mutations remains a matter of debate. We constructed a bipartite network in which 7,665 tumors from 30 cancer types are connected via shared mutations in 198 previously identified cancer genes. We show that about 27% of the tumors can be assigned to statistically supported modules, most of which encompass one or two cancer types. The rest of the tumors belong to a diffuse network component suggesting lower gene specificity of driver mutations. Linear regression of the mutational loads in cancer genes was used to estimate the number of drivers required for the onset of different cancers. The mean number of drivers in known cancer genes is approximately two, with a range of one to five. Cancers that are associated with modules had more drivers than those from the diffuse network component, suggesting that unidentified and/or interchangeable drivers exist in the latter.

Sensitivity to sequencing depth in single-cell cancer genomics

BACKGROUND

Querying cancer genomes at single-cell resolution is expected to provide a powerful framework to understand in detail the dynamics of cancer evolution. However, given the high costs currently associated with single-cell sequencing, together with the inevitable technical noise arising from single-cell genome amplification, cost-effective strategies that maximize the quality of single-cell data are critically needed. Taking advantage of previously published single-cell whole-genome and whole-exome cancer datasets, we studied the impact of sequencing depth and sampling effort towards single-cell variant detection.

METHODS

Five single-cell whole-genome and whole-exome cancer datasets were independently downscaled to 25, 10, 5, and 1× sequencing depth. For each depth level, ten technical replicates were generated, resulting in a total of 6280 single-cell BAM files. The sensitivity of variant detection, including structural and driver mutations, genotyping, clonal inference, and phylogenetic reconstruction to sequencing depth was evaluated using recent tools specifically designed for single-cell data.

RESULTS

Altogether, our results suggest that for relatively large sample sizes (25 or more cells) sequencing single tumor cells at depths > 5× does not drastically improve somatic variant discovery, characterization of clonal genotypes, or estimation of single-cell phylogenies.

CONCLUSIONS

We suggest that sequencing multiple individual tumor cells at a modest depth represents an effective alternative to explore the mutational landscape and clonal evolutionary patterns of cancer genomes.

NSAID use and somatic exomic mutations in Barrett's esophagus

Background

Use of aspirin and other non-steroidal anti-inflammatory drugs (NSAIDs) has been shown to protect against tetraploidy, aneuploidy, and chromosomal alterations in the metaplastic condition Barrett’s esophagus (BE) and to lower the incidence and mortality of esophageal adenocarcinoma (EA). The esophagus is exposed to both intrinsic and extrinsic mutagens resulting from gastric reflux, chronic inflammation, and exposure to environmental carcinogens such as those found in cigarettes. Here we test the hypothesis that NSAID use inhibits accumulation of point mutations/indels during somatic genomic evolution in BE.

Methods

Whole exome sequences were generated from 82 purified epithelial biopsies and paired blood samples from a cross-sectional study of 41 NSAID users and 41 non-users matched by sex, age, smoking, and continuous time using or not using NSAIDs.

Results

NSAID use reduced overall frequency of point mutations across the spectrum of mutation types, lowered the frequency of mutations even when adjusted for both TP53 mutation and smoking status, and decreased the prevalence of clones with high variant allele frequency. Never smokers who consistently used NSAIDs had fewer point mutations in signature 17, which is commonly found in EA. NSAID users had, on average, a 50% reduction in functional gene mutations in nine cancer-associated pathways and also had less diversity in pathway mutational burden compared to non-users.

Conclusions

These results indicate NSAID use functions to limit overall mutations on which selection can act and supports a model in which specific mutant cell populations survive or expand better in the absence of NSAIDs.

Electronic supplementary material

The online version of this article (10.1186/s13073-018-0520-y) contains supplementary material, which is available to authorized users.

Pathological Bases and Clinical Impact of Intratumor Heterogeneity in Clear Cell Renal Cell Carcinoma

PURPOSE OF REVIEW

Intratumor heterogeneity is an inherent event in tumor development that is receiving much attention in the last years since it is responsible for most failures of current targeted therapies. The purpose of this review is to offer clinicians an updated insight of the multiple manifestations of a complex event that impacts significantly patient's life.

RECENT FINDINGS

Clear cell renal cell carcinoma is the most common renal tumor and a paradigmatic example of a heterogeneous neoplasm. Next-generation sequencing has demonstrated that intratumor heterogeneity encompasses genetic, epigenetic, and microenvironmental variability. Currently accepted protocols of tumor sampling seem insufficient in unveiling intratumor heterogeneity with reliability and need to be updated. This variability challenges the precise morphological diagnosis, its molecular characterization, and the selection of optimal personalized therapies in clear cell renal cell carcinoma, a neoplasm traditionally considered chemo- and radio-resistant. We review the state of the art of the different approaches to intratumor heterogeneity in clear cell renal cell carcinomas, from the simple morphology to the most sophisticated massive sequencing tools.

Integrative genomic and transcriptomic analysis of leiomyosarcoma

Leiomyosarcoma (LMS) is an aggressive mesenchymal malignancy with few therapeutic options. The mechanisms underlying LMS development, including clinically actionable genetic vulnerabilities, are largely unknown. Here we show, using whole-exome and transcriptome sequencing, that LMS tumors are characterized by substantial mutational heterogeneity, near-universal inactivation of TP53 and RB1, widespread DNA copy number alterations including chromothripsis, and frequent whole-genome duplication. Furthermore, we detect alternative telomere lengthening in 78% of cases and identify recurrent alterations in telomere maintenance genes such as ATRX, RBL2, and SP100, providing insight into the genetic basis of this mechanism. Finally, most tumors display hallmarks of "BRCAness", including alterations in homologous recombination DNA repair genes, multiple structural rearrangements, and enrichment of specific mutational signatures, and cultured LMS cells are sensitive towards olaparib and cisplatin. This comprehensive study of LMS genomics has uncovered key biological features that may inform future experimental research and enable the design of novel therapies.

Single-Cell mRNA Sequencing in Cancer Research: Integrating the Genomic Fingerprint

Critical cancer mutations are often regional and mosaic, confounding the efficacy of targeted therapeutics. Single cell mRNA sequencing (scRNA-seq) has enabled unprecedented studies of intra-tumor heterogeneity and its role in cancer progression, metastasis, and treatment resistance. When coupled with DNA sequencing, scRNA-seq allows one to infer the in vivo impact of genomic alterations on gene expression. This combination can be used to reliably distinguish neoplastic from non-neoplastic cells, to correlate paracrine-signaling pathways between neoplastic cells and stroma, and to map expression signatures to inferred clones and phylogenies. Here we review recent advances in scRNA-seq, with a special focus on cancer. We discuss the challenges and prospects of combining scRNA-seq with DNA sequencing to assess intra-tumor heterogeneity.

The Role of Somatic L1 Retrotransposition in Human Cancers

The human LINE-1 (or L1) element is a non-LTR retrotransposon that is mobilized through an RNA intermediate by an L1-encoded reverse transcriptase and other L1-encoded proteins. L1 elements remain actively mobile today and continue to mutagenize human genomes. Importantly, when new insertions disrupt gene function, they can cause diseases. Historically, L1s were thought to be active in the germline but silenced in adult somatic tissues. However, recent studies now show that L1 is active in at least some somatic tissues, including epithelial cancers. In this review, we provide an overview of these recent developments, and examine evidence that somatic L1 retrotransposition can initiate and drive tumorigenesis in humans. Recent studies have: (i) cataloged somatic L1 activity in many epithelial tumor types; (ii) identified specific full-length L1 source elements that give rise to somatic L1 insertions; and (iii) determined that L1 promoter hypomethylation likely plays an early role in the derepression of L1s in somatic tissues. A central challenge moving forward is to determine the extent to which L1 driver mutations can promote tumor initiation, evolution, and metastasis in humans.

DNA Mutations May Not Be the Cause of Cancer

Cancer is the most challenging disease of our time with increasing numbers of new cases each year, worldwide. Great achievements have been reached in cancer research through deep sequencing which helped define druggable targets. However, the still-evolving targeted therapy suffers resistance suggesting that DNA mutations considered as drivers may not have a role in tumor initiation. The present work discusses the role of DNA mutations as drivers and passengers in cancer initiation and development. First, it is important to discern the role of these DNA mutations as initiating events causing cancer or as contributors crucial for the development of a tumor once it has initiated. Second, breast cancer shown here illustrates how identification of DNA mutations in cancerous cells has influenced our approach for anti-cancer drug design. The cancer trilogy we have reached and described as: initial drug; resistance/recurrence; drug/treatment combinations, calls for a paradigm shift. To design more effective cancer drugs with durable and positive outcome, future cancer research needs to move beyond the sequencing era and explore changes which are taking place in cancer cells at levels other than the DNA. Evolutionary constraints may be acting as a barrier to preserve the human species from being transformed and, for that matter, all multi-cellular species which can incur cancer. Furthermore, mutations in the DNA do occur and for a multitude of reasons but without necessarily causing cancer. New directions will draw themselves when more focus is given to the event responsible for the switch of a cell from normalcy to malignancy. Until then, targeted therapy will certainly continue to improve the outcome of patients; however, it is unlikely to eradicate breast cancer depicted here.

The evolution of tumour phylogenetics: principles and practice

Rapid advances in high-throughput sequencing and a growing realization of the importance of evolutionary theory to cancer genomics have led to a proliferation of phylogenetic studies of tumour progression. These studies have yielded not only new insights but also a plethora of experimental approaches, sometimes reaching conflicting or poorly supported conclusions. Here, we consider this body of work in light of the key computational principles underpinning phylogenetic inference, with the goal of providing practical guidance on the design and analysis of scientifically rigorous tumour phylogeny studies. We survey the range of methods and tools available to the researcher, their key applications, and the various unsolved problems, closing with a perspective on the prospects and broader implications of this field.

Building a lineage from single cells: genetic techniques for cell lineage tracking

Resolving lineage relationships between cells in an organism is a fundamental interest of developmental biology. Furthermore, investigating lineage can drive understanding of pathological states, including cancer, as well as understanding of developmental pathways that are amenable to manipulation by directed differentiation. Although lineage tracking through the injection of retroviral libraries has long been the state of the art, a recent explosion of methodological advances in exogenous labelling and single-cell sequencing have enabled lineage tracking at larger scales, in more detail, and in a wider range of species than was previously considered possible. In this Review, we discuss these techniques for cell lineage tracking, with attention both to those that trace lineage forwards from experimental labelling, and those that trace backwards across the life history of an organism.

Reconstructing metastatic seeding patterns of human cancers

Reconstructing the evolutionary history of metastases is critical for understanding their basic biological principles and has profound clinical implications. Genome-wide sequencing data has enabled modern phylogenomic methods to accurately dissect subclones and their phylogenies from noisy and impure bulk tumour samples at unprecedented depth. However, existing methods are not designed to infer metastatic seeding patterns. Here we develop a tool, called Treeomics, to reconstruct the phylogeny of metastases and map subclones to their anatomic locations. Treeomics infers comprehensive seeding patterns for pancreatic, ovarian, and prostate cancers. Moreover, Treeomics correctly disambiguates true seeding patterns from sequencing artifacts; 7% of variants were misclassified by conventional statistical methods. These artifacts can skew phylogenies by creating illusory tumour heterogeneity among distinct samples. In silico benchmarking on simulated tumour phylogenies across a wide range of sample purities (15–95%) and sequencing depths (25-800 × ) demonstrates the accuracy of Treeomics compared with existing methods.

Tumours frequently metastasize to multiple anatomical sites and understanding how these different metastases evolve may be important for therapy. Here, the authors develop a method—Treeomics—that can construct phylogenies from multiple metastases from next-generation sequencing data.

FISHtrees 3.0: Tumor Phylogenetics Using a Ploidy Probe

Advances in fluorescence in situ hybridization (FISH) make it feasible to detect multiple copy-number changes in hundreds of cells of solid tumors. Studies using FISH, sequencing, and other technologies have revealed substantial intra-tumor heterogeneity. The evolution of subclones in tumors may be modeled by phylogenies. Tumors often harbor aneuploid or polyploid cell populations. Using a FISH probe to estimate changes in ploidy can guide the creation of trees that model changes in ploidy and individual gene copy-number variations. We present FISHtrees 3.0, which implements a ploidy-based tree building method based on mixed integer linear programming (MILP). The ploidy-based modeling in FISHtrees includes a new formulation of the problem of merging trees for changes of a single gene into trees modeling changes in multiple genes and the ploidy. When multiple samples are collected from each patient, varying over time or tumor regions, it is useful to evaluate similarities in tumor progression among the samples. Therefore, we further implemented in FISHtrees 3.0 a new method to build consensus graphs for multiple samples. We validate FISHtrees 3.0 on a simulated data and on FISH data from paired cases of cervical primary and metastatic tumors and on paired breast ductal carcinoma in situ (DCIS) and invasive ductal carcinoma (IDC). Tests on simulated data show improved accuracy of the ploidy-based approach relative to prior ploidyless methods. Tests on real data further demonstrate novel insights these methods offer into tumor progression processes. Trees for DCIS samples are significantly less complex than trees for paired IDC samples. Consensus graphs show substantial divergence among most paired samples from both sets. Low consensus between DCIS and IDC trees may help explain the difficulty in finding biomarkers that predict which DCIS cases are at most risk to progress to IDC. The FISHtrees software is available at ftp://ftp.ncbi.nih.gov/pub/FISHtrees.

The role of Six1 signaling in paclitaxel-dependent apoptosis in MCF-7 cell line

The resistance of cancer cells to chemotherapeutic agents represents the main problem in cancer treatment. Despite intensive research, mechanisms of resistance have not yet been fully elucidated. Six1 signaling has an important role in the expansion of progenitor cell populations during early embryogenesis. Six1 gene overexpression has been strongly associated with aggressiveness, invasiveness, and poor prognosis of different cancers. In this study, we investigated the role of Six1 signaling in resistance of MCF-7 breast cancer cells to taxanes. We first established in vitro paclitaxel-resistant MCF-7 breast cancer cells. Morphological modifications in paclitaxel-resistant cells were examined via light microscopic images and fluorescence-activated cell sorting analysis. Applying quantitative real-time polymerase chain reaction, we measured Six1, B-cell lymphoma/leukemia(BCL-2), BAX, and P53 mRNA expression levels in both non-resistant and resistant cells. Resistant cells were developed from the parent MCF-7 cells by applying increasing concentrations of paclitaxel up to 64 nM. The inhibitory concentration 50% value in resistant cells increased from 3.5 ± 0.03 to 511 ± 10.22 nM (p = 0.015). In paclitaxel-resistant cells, there was a significant increase in Six1 and BCL-2 mRNA levels (p = 0.0007) with a marked decrease in pro-apoptotic Bax mRNA expression level (p = 0.03); however, there was no significant change in P53 expression (p = 0.025). Our results suggest that identifying cancer patients with high Six1 expression and then inhibition of Six1 signaling can improve the efficiency of chemotherapeutic agents in the induction of apoptosis.

Defining order and timing of mutations during cancer progression: the TO-DAG probabilistic graphical model

Somatic mutations arise and accumulate both during tumor genesis and progression. However, the order in which mutations occur is an open question and the inference of the temporal ordering at the gene level could potentially impact on patient treatment. Thus, exploiting recent observations suggesting that the occurrence of mutations is a non-memoryless process, we developed a computational approach to infer timed oncogenetic directed acyclic graphs (TO-DAGs) from human tumor mutation data. Such graphs represent the path and the waiting times of alterations during tumor evolution. The probability of occurrence of each alteration in a path is the probability that the alteration occurs when all alterations prior to it have occurred. The waiting time between an alteration and the subsequent is modeled as a stochastic function of the conditional probability of the event given the occurrence of the previous one. TO-DAG performances have been evaluated both on synthetic data and on somatic non-silent mutations from prostate cancer and melanoma patients and then compared with those of current well-established approaches. TO-DAG shows high performance scores on synthetic data and recognizes mutations in gatekeeper tumor suppressor genes as trigger for several downstream mutational events in the human tumor data.