Genomic characterization of metastatic patterns from prospective clinical sequencing of 25,000 patients.

Metastatic progression is the main cause of death in cancer patients, whereas the underlying genomic mechanisms driving metastasis remain largely unknown. Here, we assembled MSK-MET, a pan-cancer cohort of over 25,000 patients with metastatic diseases. By analyzing genomic and clinical data from this cohort, we identified associations between genomic alterations and patterns of metastatic dissemination across 50 tumor types. We found that chromosomal instability is strongly correlated with metastatic burden in some tumor types, including prostate adenocarcinoma, lung adenocarcinoma, and HR+/HER2+ breast ductal carcinoma, but not in others, including colorectal cancer and high-grade serous ovarian cancer, where copy-number alteration patterns may be established early in tumor development. We also identified somatic alterations associated with metastatic burden and specific target organs. Our data offer a valuable resource for the investigation of the biological basis for metastatic spread and highlight the complex role of chromosomal instability in cancer progression.

Interrogation of the Microenvironmental Landscape in Brain Tumors Reveals Disease-Specific Alterations of Immune Cells.

Brain malignancies encompass a range of primary and metastatic cancers, including low-grade and high-grade gliomas and brain metastases (BrMs) originating from diverse extracranial tumors. Our understanding of the brain tumor microenvironment (TME) remains limited, and it is unknown whether it is sculpted differentially by primary versus metastatic disease. We therefore comprehensively analyzed the brain TME landscape via flow cytometry, RNA sequencing, protein arrays, culture assays, and spatial tissue characterization. This revealed disease-specific enrichment of immune cells with pronounced differences in proportional abundance of tissue-resident microglia, infiltrating monocyte-derived macrophages, neutrophils, and T cells. These integrated analyses also uncovered multifaceted immune cell activation within brain malignancies entailing converging transcriptional trajectories while maintaining disease- and cell-type-specific programs. Given the interest in developing TME-targeted therapies for brain malignancies, this comprehensive resource of the immune landscape offers insights into possible strategies to overcome tumor-supporting TME properties and instead harness the TME to fight cancer.

The Human Tumor Atlas Network: Charting Tumor Transitions across Space and Time at Single-Cell Resolution.

Crucial transitions in cancer-including tumor initiation, local expansion, metastasis, and therapeutic resistance-involve complex interactions between cells within the dynamic tumor ecosystem. Transformative single-cell genomics technologies and spatial multiplex in situ methods now provide an opportunity to interrogate this complexity at unprecedented resolution. The Human Tumor Atlas Network (HTAN), part of the National Cancer Institute (NCI) Cancer Moonshot Initiative, will establish a clinical, experimental, computational, and organizational framework to generate informative and accessible three-dimensional atlases of cancer transitions for a diverse set of tumor types. This effort complements both ongoing efforts to map healthy organs and previous large-scale cancer genomics approaches focused on bulk sequencing at a single point in time. Generating single-cell, multiparametric, longitudinal atlases and integrating them with clinical outcomes should help identify novel predictive biomarkers and features as well as therapeutically relevant cell types, cell states, and cellular interactions across transitions. The resulting tumor atlases should have a profound impact on our understanding of cancer biology and have the potential to improve cancer detection, prevention, and therapeutic discovery for better precision-medicine treatments of cancer patients and those at risk for cancer.

TGF-ß Tumor Suppression through a Lethal EMT.

TGF-ß signaling can be pro-tumorigenic or tumor suppressive. We investigated this duality in pancreatic ductal adenocarcinoma (PDA), which, with other gastrointestinal cancers, exhibits frequent inactivation of the TGF-ß mediator Smad4. We show that TGF-ß induces an epithelial-mesenchymal transition (EMT), generally considered a pro-tumorigenic event. However, in TGF-ß-sensitive PDA cells, EMT becomes lethal by converting TGF-ß-induced Sox4 from an enforcer of tumorigenesis into a promoter of apoptosis. This is the result of an EMT-linked remodeling of the cellular transcription factor landscape, including the repression of the gastrointestinal lineage-master regulator Klf5. Klf5 cooperates with Sox4 in oncogenesis and prevents Sox4-induced apoptosis. Smad4 is required for EMT but dispensable for Sox4 induction by TGF-ß. TGF-ß-induced Sox4 is thus geared to bolster progenitor identity, whereas simultaneous Smad4-dependent EMT strips Sox4 of an essential partner in oncogenesis. Our work demonstrates that TGF-ß tumor suppression functions through an EMT-mediated disruption of a lineage-specific transcriptional network.

Organoid models of human and mouse ductal pancreatic cancer.

Pancreatic cancer is one of the most lethal malignancies due to its late diagnosis and limited response to treatment. Tractable methods to identify and interrogate pathways involved in pancreatic tumorigenesis are urgently needed. We established organoid models from normal and neoplastic murine and human pancreas tissues. Pancreatic organoids can be rapidly generated from resected tumors and biopsies, survive cryopreservation, and exhibit ductal- and disease-stage-specific characteristics. Orthotopically transplanted neoplastic organoids recapitulate the full spectrum of tumor development by forming early-grade neoplasms that progress to locally invasive and metastatic carcinomas. Due to their ability to be genetically manipulated, organoids are a platform to probe genetic cooperation. Comprehensive transcriptional and proteomic analyses of murine pancreatic organoids revealed genes and pathways altered during disease progression. The confirmation of many of these protein changes in human tissues demonstrates that organoids are a facile model system to discover characteristics of this deadly malignancy.

The nuclear factor ID3 endows macrophages with a potent anti-tumour activity.

Macrophage activation is controlled by a balance between activating and inhibitory receptors1-7, which protect normal tissues from excessive damage during infection8,9 but promote tumour growth and metastasis in cancer7,10. Here we report that the Kupffer cell lineage-determining factor ID3 controls this balance and selectively endows Kupffer cells with the ability to phagocytose live tumour cells and orchestrate the recruitment, proliferation and activation of natural killer and CD8 T lymphoid effector cells in the liver to restrict the growth of a variety of tumours. ID3 shifts the macrophage inhibitory/activating receptor balance to promote the phagocytic and lymphoid response, at least in part by buffering the binding of the transcription factors ELK1 and E2A at the SIRPA locus. Furthermore, loss- and gain-of-function experiments demonstrate that ID3 is sufficient to confer this potent anti-tumour activity to mouse bone-marrow-derived macrophages and human induced pluripotent stem-cell-derived macrophages. Expression of ID3 is therefore necessary and sufficient to endow macrophages with the ability to form an efficient anti-tumour niche, which could be harnessed for cell therapy in cancer.

STING inhibits the reactivation of dormant metastasis in lung adenocarcinoma.

Metastasis frequently develops from disseminated cancer cells that remain dormant after the apparently successful treatment of a primary tumour. These cells fluctuate between an immune-evasive quiescent state and a proliferative state liable to immune-mediated elimination1-6. Little is known about the clearing of reawakened metastatic cells and how this process could be therapeutically activated to eliminate residual disease in patients. Here we use models of indolent lung adenocarcinoma metastasis to identify cancer cell-intrinsic determinants of immune reactivity during exit from dormancy. Genetic screens of tumour-intrinsic immune regulators identified the stimulator of interferon genes (STING) pathway as a suppressor of metastatic outbreak. STING activity increases in metastatic progenitors that re-enter the cell cycle and is dampened by hypermethylation of the STING promoter and enhancer in breakthrough metastases or by chromatin repression in cells re-entering dormancy in response to TGFß. STING expression in cancer cells derived from spontaneous metastases suppresses their outgrowth. Systemic treatment of mice with STING agonists eliminates dormant metastasis and prevents spontaneous outbreaks in a T cell- and natural killer cell-dependent manner-these effects require cancer cell STING function. Thus, STING provides a checkpoint against the progression of dormant metastasis and a therapeutically actionable strategy for the prevention of disease relapse.

Neoantigen quality predicts immunoediting in survivors of pancreatic cancer.

Cancer immunoediting1 is a hallmark of cancer2 that predicts that lymphocytes kill more immunogenic cancer cells to cause less immunogenic clones to dominate a population. Although proven in mice1,3, whether immunoediting occurs naturally in human cancers remains unclear. Here, to address this, we investigate how 70 human pancreatic cancers evolved over 10 years. We find that, despite having more time to accumulate mutations, rare long-term survivors of pancreatic cancer who have stronger T cell activity in primary tumours develop genetically less heterogeneous recurrent tumours with fewer immunogenic mutations (neoantigens). To quantify whether immunoediting underlies these observations, we infer that a neoantigen is immunogenic (high-quality) by two features-'non-selfness'  based on neoantigen similarity to known antigens4,5, and 'selfness'  based on the antigenic distance required for a neoantigen to differentially bind to the MHC or activate a T cell compared with its wild-type peptide. Using these features, we estimate cancer clone fitness as the aggregate cost of T cells recognizing high-quality neoantigens offset by gains from oncogenic mutations. With this model, we predict the clonal evolution of tumours to reveal that long-term survivors of pancreatic cancer develop recurrent tumours with fewer high-quality neoantigens. Thus, we submit evidence that that the human immune system naturally edits neoantigens. Furthermore, we present a model to predict how immune pressure induces cancer cell populations to evolve over time. More broadly, our results argue that the immune system fundamentally surveils host genetic changes to suppress cancer.

Ordered and deterministic cancer genome evolution after p53 loss.

Although p53 inactivation promotes genomic instability1 and presents a route to malignancy for more than half of all human cancers2,3, the patterns through which heterogenous TP53 (encoding human p53) mutant genomes emerge and influence tumorigenesis remain poorly understood. Here, in a mouse model of pancreatic ductal adenocarcinoma that reports sporadic p53 loss of heterozygosity before cancer onset, we find that malignant properties enabled by p53 inactivation are acquired through a predictable pattern of genome evolution. Single-cell sequencing and in situ genotyping of cells from the point of p53 inactivation through progression to frank cancer reveal that this deterministic behaviour involves four sequential phases-Trp53 (encoding mouse p53) loss of heterozygosity, accumulation of deletions, genome doubling, and the emergence of gains and amplifications-each associated with specific histological stages across the premalignant and malignant spectrum. Despite rampant heterogeneity, the deletion events that follow p53 inactivation target functionally relevant pathways that can shape genomic evolution and remain fixed as homogenous events in diverse malignant populations. Thus, loss of p53-the 'guardian of the genome'-is not merely a gateway to genetic chaos but, rather, can enable deterministic patterns of genome evolution that may point to new strategies for the treatment of TP53-mutant tumours.

Computational modeling of pancreatic cancer reveals kinetics of metastasis suggesting optimum treatment strategies.

Pancreatic cancer is a leading cause of cancer-related death, largely due to metastatic dissemination. We investigated pancreatic cancer progression by utilizing a mathematical framework of metastasis formation together with comprehensive data of 228 patients, 101 of whom had autopsies. We found that pancreatic cancer growth is initially exponential. After estimating the rates of pancreatic cancer growth and dissemination, we determined that patients likely harbor metastases at diagnosis and predicted the number and size distribution of metastases as well as patient survival. These findings were validated in an independent database. Finally, we analyzed the effects of different treatment modalities, finding that therapies that efficiently reduce the growth rate of cells earlier in the course of treatment appear to be superior to upfront tumor resection. These predictions can be validated in the clinic. Our interdisciplinary approach provides insights into the dynamics of pancreatic cancer metastasis and identifies optimum therapeutic interventions.

The mutational landscape of human somatic and germline cells.

Over the course of an individual's lifetime, normal human cells accumulate mutations1. Here we compare the mutational landscape in 29 cell types from the soma and germline using multiple samples from the same individuals. Two ubiquitous mutational signatures, SBS1 and SBS5/40, accounted for the majority of acquired mutations in most cell types, but their absolute and relative contributions varied substantially. SBS18, which potentially reflects oxidative damage2, and several additional signatures attributed to exogenous and endogenous exposures contributed mutations to subsets of cell types. The rate of mutation was lowest in spermatogonia, the stem cells from which sperm are generated and from which most genetic variation in the human population is thought to originate. This was due to low rates of ubiquitous mutational processes and may be partially attributable to a low rate of cell division in basal spermatogonia. These results highlight similarities and differences in the maintenance of the germline and soma.

Massive genomic rearrangement acquired in a single catastrophic event during cancer development.

Cancer is driven by somatically acquired point mutations and chromosomal rearrangements, conventionally thought to accumulate gradually over time. Using next-generation sequencing, we characterize a phenomenon, which we term chromothripsis, whereby tens to hundreds of genomic rearrangements occur in a one-off cellular crisis. Rearrangements involving one or a few chromosomes crisscross back and forth across involved regions, generating frequent oscillations between two copy number states. These genomic hallmarks are highly improbable if rearrangements accumulate over time and instead imply that nearly all occur during a single cellular catastrophe. The stamp of chromothripsis can be seen in at least 2%-3% of all cancers, across many subtypes, and is present in ∼25% of bone cancers. We find that one, or indeed more than one, cancer-causing lesion can emerge out of the genomic crisis. This phenomenon has important implications for the origins of genomic remodeling and temporal emergence of cancer.

The mutational landscape of normal human endometrial epithelium.

All normal somatic cells are thought to acquire mutations, but understanding of the rates, patterns, causes and consequences of somatic mutations in normal cells is limited. The uterine endometrium adopts multiple physiological states over a lifetime and is lined by a gland-forming epithelium1,2. Here, using whole-genome sequencing, we show that normal human endometrial glands are clonal cell populations with total mutation burdens that increase at about 29 base substitutions per year and that are many-fold lower than those of endometrial cancers. Normal endometrial glands frequently carry 'driver' mutations in cancer genes, the burden of which increases with age and decreases with parity. Cell clones with drivers often originate during the first decades of life and subsequently progressively colonize the epithelial lining of the endometrium. Our results show that mutational landscapes differ markedly between normal tissues-perhaps shaped by differences in their structure and physiology-and indicate that the procession of neoplastic change that leads to endometrial cancer is initiated early in life.

Precancerous neoplastic cells can move through the pancreatic ductal system.

Most adult carcinomas develop from noninvasive precursor lesions, a progression that is supported by genetic analysis. However, the evolutionary and genetic relationships among co-existing lesions are unclear. Here we analysed the somatic variants of pancreatic cancers and precursor lesions sampled from distinct regions of the same pancreas. After inferring evolutionary relationships, we found that the ancestral cell had initiated and clonally expanded to form one or more lesions, and that subsequent driver gene mutations eventually led to invasive pancreatic cancer. We estimate that this multi-step progression generally spans many years. These new data reframe the step-wise progression model of pancreatic cancer by illustrating that independent, high-grade pancreatic precursor lesions observed in a single pancreas often represent a single neoplasm that has colonized the ductal system, accumulating spatial and genetic divergence over time.

Identification of unique neoantigen qualities in long-term survivors of pancreatic cancer.

Pancreatic ductal adenocarcinoma is a lethal cancer with fewer than 7% of patients surviving past 5 years. T-cell immunity has been linked to the exceptional outcome of the few long-term survivors, yet the relevant antigens remain unknown. Here we use genetic, immunohistochemical and transcriptional immunoprofiling, computational biophysics, and functional assays to identify T-cell antigens in long-term survivors of pancreatic cancer. Using whole-exome sequencing and in silico neoantigen prediction, we found that tumours with both the highest neoantigen number and the most abundant CD8+ T-cell infiltrates, but neither alone, stratified patients with the longest survival. Investigating the specific neoantigen qualities promoting T-cell activation in long-term survivors, we discovered that these individuals were enriched in neoantigen qualities defined by a fitness model, and neoantigens in the tumour antigen MUC16 (also known as CA125). A neoantigen quality fitness model conferring greater immunogenicity to neoantigens with differential presentation and homology to infectious disease-derived peptides identified long-term survivors in two independent datasets, whereas a neoantigen quantity model ascribing greater immunogenicity to increasing neoantigen number alone did not. We detected intratumoural and lasting circulating T-cell reactivity to both high-quality and MUC16 neoantigens in long-term survivors of pancreatic cancer, including clones with specificity to both high-quality neoantigens and predicted cross-reactive microbial epitopes, consistent with neoantigen molecular mimicry. Notably, we observed selective loss of high-quality and MUC16 neoantigenic clones on metastatic progression, suggesting neoantigen immunoediting. Our results identify neoantigens with unique qualities as T-cell targets in pancreatic ductal adenocarcinoma. More broadly, we identify neoantigen quality as a biomarker for immunogenic tumours that may guide the application of immunotherapies.

Unbiased in vivo preclinical evaluation of anticancer drugs identifies effective therapy for the treatment of pancreatic adenocarcinoma.

Pancreatic ductal adenocarcinoma (PDAC) is typically diagnosed at an advanced stage, which limits surgical options and portends a dismal prognosis. Current oncologic PDAC therapies confer marginal benefit and, thus, a significant unmet clinical need exists for new therapeutic strategies. To identify effective PDAC therapies, we leveraged a syngeneic orthotopic PDAC transplant mouse model to perform a large-scale, in vivo screen of 16 single-agent and 41 two-drug targeted therapy combinations in mice. Among 57 drug conditions screened, combined inhibition of heat shock protein (Hsp)-90 and MEK was found to produce robust suppression of tumor growth, leading to an 80% increase in the survival of PDAC-bearing mice with no significant toxicity. Mechanistically, we observed that single-agent MEK inhibition led to compensatory activation of resistance pathways, including components of the PI3K/AKT/mTOR signaling axis, which was overcome with the addition of HSP90 inhibition. The combination of HSP90(i) + MEK(i) was also active in vitro in established human PDAC cell lines and in vivo in patient-derived organoid PDAC transplant models. These findings encourage the clinical development of HSP90(i) + MEK(i) combination therapy and highlight the power of clinically relevant in vivo model systems for identifying cancer therapies.

Cell division rates decrease with age, providing a potential explanation for the age-dependent deceleration in cancer incidence.

A new evaluation of previously published data suggested to us that the accumulation of mutations might slow, rather than increase, as individuals age. To explain this unexpected finding, we hypothesized that normal stem cell division rates might decrease as we age. To test this hypothesis, we evaluated cell division rates in the epithelium of human colonic, duodenal, esophageal, and posterior ethmoid sinonasal tissues. In all 4 tissues, there was a significant decrease in cell division rates with age. In contrast, cell division rates did not decrease in the colon of aged mice, and only small decreases were observed in their small intestine or esophagus. These results have important implications for understanding the relationship between normal stem cells, aging, and cancer. Moreover, they provide a plausible explanation for the enigmatic age-dependent deceleration in cancer incidence in very old humans but not in mice.

Comparison of immune infiltrates in melanoma and pancreatic cancer highlights VISTA as a potential target in pancreatic cancer.

Immune checkpoint therapy (ICT) has transformed cancer treatment in recent years; however, treatment response is not uniform across tumor types. The tumor immune microenvironment plays a critical role in determining response to ICT; therefore, understanding the differential immune infiltration between ICT-sensitive and ICT-resistant tumor types will help to develop effective treatment strategies. We performed a comprehensive analysis of the immune tumor microenvironment of an ICT-sensitive tumor (melanoma, n = 44) and an ICT-resistant tumor (pancreatic cancer, n = 67). We found that a pancreatic tumor has minimal to moderate infiltration of CD3, CD4, and CD8 T cells; however, the immune infiltrates are predominantly present in the stromal area of the tumor and are excluded from tumoral area compared with melanoma, where the immune infiltrates are primarily present in the tumoral area. Metastatic pancreatic ductal adenocarcinomas (PDACs) had a lower infiltration of total T cells compared with resectable primary PDACs, suggesting that metastatic PDACs have poor immunogenicity. Further, a significantly higher number of CD68+ macrophages and VISTA+ cells (also known as V-domain immunoglobulin suppressor of T cell activation) were found in the pancreatic stromal area compared with melanoma. We identified VISTA as a potent inhibitory checkpoint that is predominantly expressed on CD68+ macrophages on PDACs. These data suggest that VISTA may be a relevant immunotherapy target for effective treatment of patients with pancreatic cancer.

Pancreas cancer and BRCA: A critical subset of patients with improving therapeutic outcomes.

BACKGROUND: Patients with germline/somatic BRCA1/BRCA2 mutations (g/sBRCA1/2) comprise a distinct biologic subgroup of pancreas ductal adenocarcinoma (PDAC). METHODS: Institutional databases were queried to identify patients who had PDAC with g/sBRCA1/2. Demographics, clinicopathologic details, genomic data (annotation sBRCA1/2 according to a precision oncology knowledge base for somatic mutations), zygosity, and outcomes were abstracted. Overall survival (OS) was estimated using the Kaplan-Meier method. RESULTS: In total, 136 patients with g/sBRCA1/2 were identified between January 2011 and June 2020. Germline BRCA1/2 (gBRCA1/2) mutation was identified in 116 patients (85%). Oncogenic somatic BRCA1/2 (sBRCA1/2) mutation was present in 20 patients (15%). Seventy-seven patients had biallelic BRCA1/2 mutations (83%), and 16 (17%) had heterozygous mutations. Sixty-five patients with stage IV disease received frontline platinum therapy, and 52 (80%) had a partial response. The median OS for entire cohort was 27.6 months (95% CI, 24.9-34.5 months), and the median OS for patients who had stage IV disease was 23 months (95% CI, 19-26 months). Seventy-one patients received a poly(adenosine diphosphate ribose) polymerase (PARP) inhibitor (PARPi), and 52 received PARPi monotherapy. For maintenance PARPi, 10 patients (36%) had a partial response, 12 (43%) had stable disease, and 6 (21%) had progression of disease as their best response. Six patients (21%) received maintenance PARPi for >2 years. For those with stage IV disease who received frontline platinum, the median OS was 26 months (95% CI, 20-52 months) for biallelic patients (n = 39) and 8.66 months (95% CI, 6.2 months to not reached) for heterozygous patients (n = 4). The median OS for those who received PARPi therapy was 26.5 months (95% CI, 24-53 months) for biallelic patients (n = 25) and 8.66 months (95% CI, 7.23 months to not reached) for heterozygous patients (n = 2). CONCLUSIONS: g/sBRCA1/2 mutations did not appear to have different actionable utility. Platinum and PARPi therapies offer therapeutic benefit, and very durable outcomes are observed in a subset of patients who have g/sBRCA1/2 mutations with biallelic status.

Pancreatic cancer stem cells may define tumor stroma characteristics and recurrence patterns in pancreatic ductal adenocarcinoma.

BACKGROUND: Herein, we investigate the relationship between pancreatic stem cell markers (PCSC markers), CD44, and epithelial-specific antigen (ESA), tumor stroma, and the impact on recurrence outcomes in pancreatic ductal adenocarcinoma (PDAC) patients. METHODS: PDAC patients who underwent surgical resection between 01/2012-06/2014 were identified. CD44 and ESA expression was assessed by immunohistochemistry. Stroma was classified as loose, moderate, and dense based on fibroblast content. Overall survival (OS) and relapse-free survival (RFS) were estimated using the Kaplan-Meier method and compared between subgroups by log-rank test. The association between PCSC markers and stroma type was assessed by Fisher's exact test. RESULTS: N = 93 PDAC patients were identified. The number of PDAC patients with dense, moderate density, and loose stroma was 11 (12%), 51 (54%), and 31 (33%) respectively. PDAC with CD44+/ESA- had highest rate of loose stroma (63%) followed by PDAC CD44+/ESA+ (50%), PDAC CD44-/ESA+ (35%), CD44-/ESA- (9%) (p = 0.0033). Conversely, lack of CD44 and ESA expression was associated with the highest rate of moderate and dense stroma (91% p = 0.0033). No local recurrence was observed in patients with dense stroma and 9 had distant recurrence. The highest rate of cumulative local recurrence was observed in patients with loose stroma. No statistically significant difference in RFS and OS was observed among subgroups (P = 0.089). CONCLUSIONS: These data indicate PCSCs may have an important role in stroma differentiation in PDAC. Our results further suggest that tumor stroma may influence the recurrence pattern in PDAC patients.