Stromal changes in the aged lung induce an emergence from melanoma dormancy.

Disseminated cancer cells from primary tumours can seed in distal tissues, but may take several years to form overt metastases, a phenomenon that is termed tumour dormancy. Despite its importance in metastasis and residual disease, few studies have been able to successfully characterize dormancy within melanoma. Here we show that the aged lung microenvironment facilitates a permissive niche for efficient outgrowth of dormant disseminated cancer cells-in contrast to the aged skin, in which age-related changes suppress melanoma growth but drive dissemination. These microenvironmental complexities can be explained by the phenotype switching model, which argues that melanoma cells switch between a proliferative cell state and a slower-cycling, invasive state1-3. It was previously shown that dermal fibroblasts promote phenotype switching in melanoma during ageing4-8. We now identify WNT5A as an activator of dormancy in melanoma disseminated cancer cells within the lung, which initially enables the efficient dissemination and seeding of melanoma cells in metastatic niches. Age-induced reprogramming of lung fibroblasts increases their secretion of the soluble WNT antagonist sFRP1, which inhibits WNT5A in melanoma cells and thereby enables efficient metastatic outgrowth. We also identify the tyrosine kinase receptors AXL and MER as promoting a dormancy-to-reactivation axis within melanoma cells. Overall, we find that age-induced changes in distal metastatic microenvironments promote the efficient reactivation of dormant melanoma cells in the lung.

Profiling of syngeneic mouse HCC tumor models as a framework to understand anti-PD-1 sensitive tumor microenvironments.

BACKGROUND AND AIMS: The treatment of hepatocellular carcinoma (HCC) has been transformed by the use of immune checkpoint inhibitors. However, most patients with HCC do not benefit from treatment with immunotherapy. There is an urgent need to understand the mechanisms that underlie response or resistance to immunotherapy for patients with HCC. The use of syngeneic mouse models that closely recapitulate the heterogeneity of human HCC will provide opportunities to examine the complex interactions between cancer cells and nonmalignant cells in the tumor microenvironment. APPROACH AND RESULTS: We leverage a multifaceted approach that includes imaging mass cytometry and suspension cytometry by time of flight to profile the tumor microenvironments of the Hep53.4, Hepa 1-6, RIL-175, and TIBx (derivative of TIB-75) syngeneic mouse HCC models. The immune tumor microenvironments vary across these four models, and various immunosuppressive pathways exist at baseline in orthotopic liver tumors derived from these models. For instance, TIBx, which is resistant to anti-programmed cell death protein 1 therapy, contains a high proportion of "M2-like" tumor-associated macrophages with the potential to diminish antitumor immunity. Investigation of The Cancer Genome Atlas reveals that the baseline immunologic profiles of Hep53.4, RIL-175, and TIBx are broadly representative of human HCCs; however, Hepa 1-6 does not recapitulate the immune tumor microenvironment of the vast majority of human HCCs. CONCLUSIONS: There is a wide diversity in the immune tumor microenvironments in preclinical models and in human HCC, highlighting the need to use multiple syngeneic HCC models to improve the understanding of how to treat HCC through immune modulation.

The estrogen receptor-alpha S118P variant does not affect breast cancer incidence or response to endocrine therapies.

PURPOSE: Estrogen receptor-alpha (ER) is a therapeutic target of ER-positive (ER+) breast cancers. Although ER signaling is complex, many mediators of this pathway have been identified. Specifically, phosphorylation of ER at serine 118 affects responses to estrogen and therapeutic ligands and has been correlated with clinical outcomes in ER+ breast cancer patients. We hypothesized that a newly described germline variant (S118P) at this residue would drive cellular changes consistent with breast cancer development and/or hormone resistance. METHODS: Isogenic human breast epithelial cell line models harboring ER S118P were developed via genome editing and characterized to determine the functional effects of this variant. We also examined the frequency of ER S118P in a case-control study (N = 536) of women with and without breast cancer with a familial risk. RESULTS: In heterozygous knock-in models, the S118P variant demonstrated no significant change in proliferation, migration, MAP Kinase pathway signaling, or response to the endocrine therapies tamoxifen and fulvestrant. Further, there was no difference in the prevalence of S118P between women with and without cancer relative to population registry databases. CONCLUSIONS: This study suggests that the ER S118P variant does not affect risk for breast cancer or hormone therapy resistance. Germline screening and modification of treatments for patients harboring this variant are likely not warranted.

NDRG1 links p53 with proliferation-mediated centrosome homeostasis and genome stability.

The tumor protein 53 (TP53) tumor suppressor gene is the most frequently somatically altered gene in human cancers. Here we show expression of N-Myc down-regulated gene 1 (NDRG1) is induced by p53 during physiologic low proliferative states, and mediates centrosome homeostasis, thus maintaining genome stability. When placed in physiologic low-proliferating conditions, human TP53 null cells fail to increase expression of NDRG1 compared with isogenic wild-type controls and TP53 R248W knockin cells. Overexpression and RNA interference studies demonstrate that NDRG1 regulates centrosome number and amplification. Mechanistically, NDRG1 physically associates with γ-tubulin, a key component of the centrosome, with reduced association in p53 null cells. Strikingly, TP53 homozygous loss was mutually exclusive of NDRG1 overexpression in over 96% of human cancers, supporting the broad applicability of these results. Our study elucidates a mechanism of how TP53 loss leads to abnormal centrosome numbers and genomic instability mediated by NDRG1.

HER2 missense mutations have distinct effects on oncogenic signaling and migration.

Recurrent human epidermal growth factor receptor 2 (HER2) missense mutations have been reported in human cancers. These mutations occur primarily in the absence of HER2 gene amplification such that most HER2-mutant tumors are classified as "negative" by FISH or immunohistochemistry assays. It remains unclear whether nonamplified HER2 missense mutations are oncogenic and whether they are targets for HER2-directed therapies that are currently approved for the treatment of HER2 gene-amplified breast cancers. Here we functionally characterize HER2 kinase and extracellular domain mutations through gene editing of the endogenous loci in HER2 nonamplified human breast epithelial cells. In in vitro and in vivo assays, the majority of HER2 missense mutations do not impart detectable oncogenic changes. However, the HER2 V777L mutation increased biochemical pathway activation and, in the context of a PIK3CA mutation, enhanced migratory features in vitro. However, the V777L mutation did not alter in vivo tumorigenicity or sensitivity to HER2-directed therapies in proliferation assays. Our results suggest the oncogenicity and potential targeting of HER2 missense mutations should be considered in the context of cooperating genetic alterations and provide previously unidentified insights into functional analysis of HER2 mutations and strategies to target them.

PIK3CA mutations and TP53 alterations cooperate to increase cancerous phenotypes and tumor heterogeneity.

BACKGROUND/PURPOSE: The combined contributions of oncogenes and tumor suppressor genes toward carcinogenesis remain poorly understood. Elucidation of cancer gene cooperativity can provide new insights leading to more effective use of therapies. EXPERIMENTAL DESIGN/METHODS: We used somatic cell genome editing to introduce singly and in combination PIK3CA mutations (E545K or H1047R) with TP53 alterations (R248W or knockout), to assess any enhanced cancerous phenotypes. The non-tumorigenic human breast epithelial cell line, MCF10A, was used as the parental cell line, and resultant cells were assessed via various in vitro assays, growth as xenografts, and drug sensitivity assays using targeted agents and chemotherapies. RESULTS: Compared to single-gene-targeted cells and parental controls, cells with both a PIK3CA mutation and TP53 alteration had increased cancerous phenotypes including cell proliferation, soft agar colony formation, aberrant morphology in acinar formation assays, and genomic heterogeneity. Cells also displayed varying sensitivities to anti-neoplastic drugs, although all cells with PIK3CA mutations showed a relative increased sensitivity to paclitaxel. All cell lines remained non-tumorigenic. CONCLUSIONS: This cell line panel provides a resource for further elucidating cooperative genetic mediators of carcinogenesis and response to therapies.

Single-Nucleotide Polymorphism Leading to False Allelic Fraction by Droplet Digital PCR.

BACKGROUND: Molecular-based diagnostics have great utility for cancer detection. We have used droplet digital PCR (ddPCR) as a platform for identifying mutations in circulating plasma tumor DNA (ptDNA). We present the unexpected finding of a spurious mutant allele fraction that was discovered to be artifactual because of the presence of a single-nucleotide polymorphism (SNP) in a patient sample. DESIGN AND METHODS: Probe and primer combinations for the K700 and V701 loci of the SF3B1 spliceosome gene were designed for ddPCR to identify the percentage of mutant and wild-type alleles. Clinical samples from patients with cancer with known SF3B1 mutations were collected and tested to evaluate the assays' ability to detect SF3B1 mutations in ptDNA. RESULTS: Patient samples showed SF3B1 K700E mutations within the ptDNA of 4 patients with acute leukemia and 3 with myelodysplastic syndrome who were known to harbor this mutation. A blood sample from a patient with lung cancer with a known SF3B1 V701F mutation was also analyzed and this mutation was successfully identified in ptDNA. However, 1 of the patients with a K700E mutation was found to have a mutational burden of 98%. After careful analysis of this locus by Sanger sequencing and ddPCR, this patient was found to have an SNP (R702R), which prevented binding of the ddPCR wild-type probe to its cognate allele. CONCLUSIONS: These results further support that ddPCR-based assays may be valuable companion diagnostics for the identification and monitoring of patients with cancer, but the results also emphasize the need to identify SNPs at loci that are being analyzed.

MACROD2 overexpression mediates estrogen independent growth and tamoxifen resistance in breast cancers.

Tamoxifen is effective for treating estrogen receptor-alpha (ER) positive breast cancers. However, few molecular mediators of tamoxifen resistance have been elucidated. Here we describe a previously unidentified gene, MACROD2 that confers tamoxifen resistance and estrogen independent growth. We found MACROD2 is amplified and overexpressed in metastatic tamoxifen-resistant tumors. Transgene overexpression of MACROD2 in breast cancer cell lines results in tamoxifen resistance, whereas RNAi-mediated gene knock down reverses this phenotype. MACROD2 overexpression also leads to estrogen independent growth in xenograft assays. Mechanistically, MACROD2 increases p300 binding to estrogen response elements in a subset of ER regulated genes. Primary breast cancers and matched metastases demonstrate MACROD2 expression can change with disease evolution, and increased expression and amplification of MACROD2 in primary tumors is associated with worse overall survival. These studies establish MACROD2 as a key mediator of estrogen independent growth and tamoxifen resistance, as well as a potential novel target for diagnostics and therapy.

A phosphoproteomic screen demonstrates differential dependence on HER3 for MAP kinase pathway activation by distinct PIK3CA mutations.

The PIK3CA gene encodes for the p110 alpha isoform of PI3 kinase and is one of the most frequently mutated oncogenes in human cancers. However, the mechanisms by which PIK3CA mutations activate cell signaling are not fully understood. Here we used a phosphoproteomic approach to compare differential phosphorylation patterns between human breast epithelial cells and two isogenic somatic cell knock in derivatives, each harboring a distinct PIK3CA mutation. We demonstrated differential phosphorylation patterns between isogenic cell lines containing a PIK3CA helical domain mutation (E545K) compared to cells with a PIK3CA kinase domain mutation (H1047R). In particular, the receptor tyrosine kinase, HER3, showed increased phosphorylation at tyrosine 1328 in H1047R cells versus E545K cells. Genetic studies using shRNA demonstrated that H1047R cells have a profound decrease in growth factor independent proliferation upon HER3 knock down, but this effect was attenuated in E545K cells. In addition, HER3 knock down led to reductions in both PI3 kinase and MAP kinase pathway activation in H1047R cells, but in E545K cells only PI3 kinase pathway diminution was observed. These studies demonstrate the power of using paired isogenic cell lines for proteomic analysis to gain new insights into oncogenic signal transduction pathways.

Fibroblasts in the Aged Pancreas Drive Pancreatic Cancer Progression.

Pancreatic cancer is more prevalent in older individuals and often carries a poorer prognosis for them. The relationship between the microenvironment and pancreatic cancer is multifactorial, and age-related changes in nonmalignant cells in the tumor microenvironment may play a key role in promoting cancer aggressiveness. Because fibroblasts have profound impacts on pancreatic cancer progression, we investigated whether age-related changes in pancreatic fibroblasts influence cancer growth and metastasis. Proteomics analysis revealed that aged fibroblasts secrete different factors than young fibroblasts, including increased growth/differentiation factor 15 (GDF-15). Treating young mice with GDF-15 enhanced tumor growth, whereas aged GDF-15 knockout mice showed reduced tumor growth. GDF-15 activated AKT, rendering tumors sensitive to AKT inhibition in an aged but not young microenvironment. These data provide evidence for how aging alters pancreatic fibroblasts and promotes tumor progression, providing potential therapeutic targets and avenues for studying pancreatic cancer while accounting for the effects of aging. SIGNIFICANCE: Aged pancreatic fibroblasts secrete GDF-15 and activate AKT signaling to promote pancreatic cancer growth, highlighting the critical role of aging-mediated changes in the pancreatic cancer microenvironment in driving tumor progression. See related commentary by Isaacson et al., p. 1185.

ERK hyperactivation serves as a unified mechanism of escape in intrinsic and acquired CDK4/6 inhibitor resistance in acral lentiginous melanoma.

Patients with metastatic acral lentiginous melanoma (ALM) suffer worse outcomes relative to patients with other forms of cutaneous melanoma (CM), and do not benefit as well to approved melanoma therapies. Identification of cyclin-dependent kinase 4 and 6 (CDK4/6) pathway gene alterations in >60% of ALMs has led to clinical trials of the CDK4/6 inhibitor (CDK4i/6i) palbociclib for ALM; however, median progression free survival with CDK4i/6i treatment was only 2.2 months, suggesting existence of resistance mechanisms. Therapy resistance in ALM remains poorly understood; here we report hyperactivation of MAPK signaling and elevated cyclin D1 expression serve as a mechanism of intrinsic early/adaptive CDK4i/6i resistance. ALM cells that have acquired CDK4i/6i resistance following chronic treatment exposure also exhibit hyperactivation of the MAPK pathway. MEK and/or ERK inhibition increases CDK4i/6i efficacy against therapy naïve and CDK4i/6i-resistant AM cells in xenograft and patient-derived xenograft (PDX) models and promotes a defective DNA repair, cell cycle arrested and apoptotic program. Notably, gene alterations poorly correlate with protein expression of cell cycle proteins in ALM or efficacy of CDK4i/6i, urging additional strategies when stratifying patients for CDK4i/6i trial inclusion. Concurrent targeting of the MAPK pathway and CDK4/6 represents a new approach for patients with metastatic ALM to improve outcomes.

Oncogene-induced matrix reorganization controls CD8+ T cell function in the soft-tissue sarcoma microenvironment.

CD8+ T cell dysfunction impedes antitumor immunity in solid cancers, but the underlying mechanisms are diverse and poorly understood. Extracellular matrix (ECM) composition has been linked to impaired T cell migration and enhanced tumor progression; however, impacts of individual ECM molecules on T cell function in the tumor microenvironment (TME) are only beginning to be elucidated. Upstream regulators of aberrant ECM deposition and organization in solid tumors are equally ill-defined. Therefore, we investigated how ECM composition modulates CD8+ T cell function in undifferentiated pleomorphic sarcoma (UPS), an immunologically active desmoplastic tumor. Using an autochthonous murine model of UPS and data from multiple human patient cohorts, we discovered a multifaceted mechanism wherein the transcriptional coactivator YAP1 promotes collagen VI (COLVI) deposition in the UPS TME. In turn, COLVI induces CD8+ T cell dysfunction and immune evasion by remodeling fibrillar collagen and inhibiting T cell autophagic flux. Unexpectedly, collagen I (COLI) opposed COLVI in this setting, promoting CD8+ T cell function and acting as a tumor suppressor. Thus, CD8+ T cell responses in sarcoma depend on oncogene-mediated ECM composition and remodeling.

Transfer Learning Reveals Cancer-Associated Fibroblasts Are Associated with Epithelial-Mesenchymal Transition and Inflammation in Cancer Cells in Pancreatic Ductal Adenocarcinoma.

Pancreatic ductal adenocarcinoma (PDAC) is an aggressive malignancy characterized by an immunosuppressive tumor microenvironment enriched with cancer-associated fibroblasts (CAF). This study used a convergence approach to identify tumor cell and CAF interactions through the integration of single-cell data from human tumors with human organoid coculture experiments. Analysis of a comprehensive atlas of PDAC single-cell RNA sequencing data indicated that CAF density is associated with increased inflammation and epithelial-mesenchymal transition (EMT) in epithelial cells. Transfer learning using transcriptional data from patient-derived organoid and CAF cocultures provided in silico validation of CAF induction of inflammatory and EMT epithelial cell states. Further experimental validation in cocultures demonstrated integrin beta 1 (ITGB1) and vascular endothelial factor A (VEGFA) interactions with neuropilin-1 mediating CAF-epithelial cell cross-talk. Together, this study introduces transfer learning from human single-cell data to organoid coculture analyses for experimental validation of discoveries of cell-cell cross-talk and identifies fibroblast-mediated regulation of EMT and inflammation. SIGNIFICANCE: Adaptation of transfer learning to relate human single-cell RNA sequencing data to organoid-CAF cocultures facilitates discovery of human pancreatic cancer intercellular interactions and uncovers cross-talk between CAFs and tumor cells through VEGFA and ITGB1.

Immune landscape of tertiary lymphoid structures in hepatocellular carcinoma (HCC) treated with neoadjuvant immune checkpoint blockade.

Neoadjuvant immunotherapy is thought to produce long-term remissions through induction of antitumor immune responses before removal of the primary tumor. Tertiary lymphoid structures (TLS), germinal center-like structures that can arise within tumors, may contribute to the establishment of immunological memory in this setting, but understanding of their role remains limited. Here, we investigated the contribution of TLS to antitumor immunity in hepatocellular carcinoma (HCC) treated with neoadjuvant immunotherapy. We found that neoadjuvant immunotherapy induced the formation of TLS, which were associated with superior pathologic response, improved relapse free survival, and expansion of the intratumoral T and B cell repertoire. While TLS in viable tumor displayed a highly active mature morphology, in areas of tumor regression we identified an involuted TLS morphology, which was characterized by dispersion of the B cell follicle and persistence of a T cell zone enriched for ongoing antigen presentation and T cell-mature dendritic cell interactions. Involuted TLS showed increased expression of T cell memory markers and expansion of CD8+ cytotoxic and tissue resident memory clonotypes. Collectively, these data reveal the circumstances of TLS dissolution and suggest a functional role for late-stage TLS as sites of T cell memory formation after elimination of viable tumor.

Lenalidomide modulates IL-8 and anti-prostate antibody levels in men with biochemically recurrent prostate cancer.

BACKGROUND: We retrospectively explored changes in immunological parameters in men with biochemically recurrent prostate cancer treated with either 5 or 25 mg of lenalidomide in a randomized phase 2 trial, and determined whether those changes correlated with disease progression. METHODS: Cytokine levels were compared for each patient at baseline and after 6 months of treatment with lenalidomide. Regression models for correlated data were used to assess associations of cytokine levels with lenalidomide treatment effect. Estimates were obtained using generalized estimating equations. Changes in circulating anti-prostate antibodies were evaluated using a high-throughput immunoblot technique. RESULTS: Treatment with lenalidomide was associated with global changes in immunoreactivity to a number of prostate-associated antigens, as well as with changes in circulating levels of the T(H) 2 cytokines IL-4, IL-5, IL-10, and IL-13. Disease progression in treated patients was associated with an increase in circulating IL-8 levels, while IL-8 levels decreased significantly in non-progressors. CONCLUSIONS: Lenalidomide demonstrates immunomodulatory properties in patients with biochemically recurrent prostate cancer. The induction of novel anti-prostate antibodies is a potential mechanism for lenalidomide response. Changes in serum IL-8 levels may serve as a potential biomarker in treated patients. These hypotheses require formal testing in future prospective trials.

Shared genetic and epigenetic changes link aging and cancer.

Aging is a universal biological process that increases the risk of multiple diseases including cancer. Growing evidence shows that alterations in the genome and epigenome, driven by similar mechanisms, are found in both aged cells and cancer cells. In this review, we detail the genetic and epigenetic changes associated with normal aging and the mechanisms responsible for these changes. By highlighting genetic and epigenetic alterations in the context of tumorigenesis, cancer progression, and the aging tumor microenvironment, we examine the possible impacts of the normal aging process on malignant transformation. Finally, we examine the implications of age-related genetic and epigenetic alterations in both tumors and patients for the treatment of cancer.

Clinical, Genomic, and Transcriptomic Data Profiling of Biliary Tract Cancer Reveals Subtype-Specific Immune Signatures.

PURPOSE: Biliary tract cancers (BTCs) are aggressive cancers that carry a poor prognosis. An enhanced understanding of the immune landscape of anatomically and molecularly defined subsets of BTC may improve patient selection for immunotherapy and inform immune-based combination treatment strategies. METHODS: We analyzed deidentified clinical, genomic, and transcriptomic data from the Tempus database to determine the mutational frequency and mutational clustering across the three major BTC subtypes (intrahepatic cholangiocarcinoma [IHC], extrahepatic cholangiocarcinoma, and gallbladder cancer). We subsequently determined the relationship between specific molecular alterations and anatomical subsets and features of the BTC immune microenvironment. RESULTS: We analyzed 454 samples of BTC, of which the most commonly detected alterations were TP53 (42.5%), CDKN2A (23.4%), ARID1A (19.6%), BAP1 (15.5%), KRAS (15%), CDKN2B (14.2%), PBRM1 (11.7%), IDH1 (11.7%), TERT (8.4%), KMT2C (10.4%) and LRP1B (8.4%), and FGFR2 fusions (8.7%). Potentially actionable molecular alterations were identified in 30.5% of BTCs including 39.1% of IHC. Integrative cluster analysis revealed four distinct molecular clusters, with cluster 4 predominately associated with FGFR2 rearrangements and BAP1 mutations in IHC. Immune-related biomarkers indicative of an inflamed tumor-immune microenvironment were elevated in gallbladder cancers and in cluster 1, which was enriched for TP53, KRAS, and ATM mutations. Multiple common driver genes, including TP53, FGFR2, IDH1, TERT, BRAF, and BAP1, were individually associated with unique BTC immune microenvironments. CONCLUSION: BTC subtypes exhibit diverse DNA alterations, RNA inflammatory signatures, and immune biomarkers. The association between specific BTC anatomical subsets, molecular alterations, and immunophenotypes highlights new opportunities for therapeutic development.

NOTCH1 PEST domain variants are responsive to standard of care treatments despite distinct transformative properties in a breast cancer model.

Activating variants in the PEST region of NOTCH1 have been associated with aggressive phenotypes in human cancers, including triple-negative breast cancer (TNBC). Previous studies suggested that PEST domain variants in TNBC patients resulted in increased cell proliferation, invasiveness, and decreased overall survival. In this study, we assess the phenotypic transformation of activating NOTCH1 variants and their response to standard of care therapies. AAV-mediated gene targeting was used to isogenically incorporate 3 NOTCH1 variants, including a novel TNBC frameshift variant, in two non-tumorigenic breast epithelial cell lines, MCF10A and hTERT-IMEC. Two different variants at the NOTCH1 A2241 site (A2441fs and A2441T) both demonstrated increased transformative properties when compared to a non-transformative PEST domain variant (S2523L). These phenotypic changes include proliferation, migration, anchorage-independent growth, and MAPK pathway activation. In contrast to previous studies, activating NOTCH1 variants did not display sensitivity to a gamma secretase inhibitor (GSI) or resistance to chemotherapies. This study demonstrates distinct transformative phenotypes are specific to a given variant within NOTCH1 and these phenotypes do not correlate with sensitivities or resistance to chemotherapies or GSIs. Although previous studies have suggested NOTCH1 variants may be prognostic for TNBC, our study does not demonstrate prognostic ability of these variants and suggests further characterization would be required for clinical applications.

Hierarchical tumor heterogeneity mediated by cell contact between distinct genetic subclones.

Intratumor heterogeneity is an important mediator of poor outcomes in many cancers, including breast cancer. Genetic subclones frequently contribute to this heterogeneity; however, their growth dynamics and interactions remain poorly understood. PIK3CA and HER2 alterations are known to coexist in breast and other cancers. Herein, we present data that describe the ability of oncogenic PIK3CA mutant cells to induce the proliferation of quiescent HER2 mutant cells through a cell contact-mediated mechanism. Interestingly, the HER2 cells proliferated to become the major subclone over PIK3CA counterparts both in vitro and in vivo. Furthermore, this phenotype was observed in both hormone receptor-positive and -negative cell lines, and was dependent on the expression of fibronectin from mutant PIK3CA cells. Analysis of human tumors demonstrated similar HER2:PIK3CA clonal dynamics and fibronectin expression. Our study provides insight into nonrandom subclonal architecture of heterogenous tumors, which may aid the understanding of tumor evolution and inform future strategies for personalized medicine.