Metastatic site influences driver gene function in pancreatic cancer.

Driver gene mutations can increase the metastatic potential of the primary tumor1-3, but their role in sustaining tumor growth at metastatic sites is poorly understood. A paradigm of such mutations is inactivation of SMAD4 - a transcriptional effector of TGFß signaling - which is a hallmark of multiple gastrointestinal malignancies4,5. SMAD4 inactivation mediates TGFß's remarkable anti- to pro-tumorigenic switch during cancer progression and can thus influence both tumor initiation and metastasis6-14. To determine whether metastatic tumors remain dependent on SMAD4 inactivation, we developed a mouse model of pancreatic ductal adenocarcinoma (PDAC) that enables Smad4 depletion in the pre-malignant pancreas and subsequent Smad4 reactivation in established metastases. As expected, Smad4 inactivation facilitated the formation of primary tumors that eventually colonized the liver and lungs. By contrast, Smad4 reactivation in metastatic disease had strikingly opposite effects depending on the tumor's organ of residence: suppression of liver metastases and promotion of lung metastases. Integrative multiomic analysis revealed organ-specific differences in the tumor cells' epigenomic state, whereby the liver and lungs harbored chromatin programs respectively dominated by the KLF and RUNX developmental transcription factors, with Klf4 depletion being sufficient to reverse Smad4's tumor-suppressive activity in liver metastases. Our results show how epigenetic states favored by the organ of residence can influence the function of driver genes in metastatic tumors. This organ-specific gene-chromatin interplay invites consideration of anatomical site in the interpretation of tumor genetics, with implications for the therapeutic targeting of metastatic disease.

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.

Concordance in Oncogenic Alterations Between the Primary Tumor and Advanced Disease: Insights Into the Heterogeneity of Intrahepatic Cholangiocarcinoma.

PURPOSE: Intrahepatic cholangiocarcinoma (ICCA) is characterized by significant phenotypic and clinical heterogeneities and poor response to systemic therapy, potentially related to underlying heterogeneity in oncogenic alterations. We aimed to characterize the genomic heterogeneity between primary tumors and advanced disease in patients with ICCA. METHODS: Biopsy-proven CCA specimens (primary tumor and paired advanced disease [metastatic disease, progressive disease on systemic therapy, or postoperative recurrence]) from two institutions were subjected to targeted next-generation sequencing. Overall concordance (oncogenic driver mutations, copy number alterations, and fusion events) and mutational concordance (only oncogenic mutations) were compared across paired samples. A subgroup analysis was performed on the basis of exposure to systemic therapy. Patients with extrahepatic CCA (ECCA) were included as a comparison group. RESULTS: Sample pairs from 65 patients with ICCA (n = 54) and ECCA (n = 11) were analyzed. The median time between sample collection was 19.6 months (range, 2.7-122.9). For the entire cohort, the overall oncogenic concordance was 49% and the mutational concordance was 62% between primary and advanced disease samples. Subgroup analyses of ICCA and ECCA revealed overall/mutational concordance rates of 47%/58% and 60%/84%, respectively. Oncogenic concordance was similarly low for pairs exposed to systemic therapy between sample collections (n = 50, 53% overall, 68% mutational). In patients treated with targeted therapy for IDH1/2 alterations (n = 6) or FGFR2 fusions (n = 3), there was 100% concordance between the primary and advanced disease specimens. In two patients, FGFR2 (n = 1) and IDH1 (n = 1) alterations were detected de novo in the advanced disease specimens. CONCLUSION: The results reflect a high degree of heterogeneity in ICCA and argue for reassessment of the dominant driver mutations with change in disease status.

A microglia clonal inflammatory disorder in Alzheimer's Disease.

Somatic genetic heterogeneity resulting from post-zygotic DNA mutations is widespread in human tissues and can cause diseases, however few studies have investigated its role in neurodegenerative processes such as Alzheimer's Disease (AD). Here we report the selective enrichment of microglia clones carrying pathogenic variants, that are not present in neuronal, glia/stromal cells, or blood, from patients with AD in comparison to age-matched controls. Notably, microglia-specific AD-associated variants preferentially target the MAPK pathway, including recurrent CBL ring-domain mutations. These variants activate ERK and drive a microglia transcriptional program characterized by a strong neuro-inflammatory response, both in vitro and in patients. Although the natural history of AD-associated microglial clones is difficult to establish in human, microglial expression of a MAPK pathway activating variant was previously shown to cause neurodegeneration in mice, suggesting that AD-associated neuroinflammatory microglial clones may contribute to the neurodegenerative process in patients.

Islands of genomic stability in the face of genetically unstable metastatic cancer.

Introduction: Metastatic cancer affects millions of people worldwide annually and is the leading cause of cancer-related deaths. Most patients with metastatic disease are not eligible for surgical resection, and current therapeutic regimens have varying success rates, some with 5-year survival rates below 5%. Here we test the hypothesis that metastatic cancer can be genetically targeted by exploiting single base substitution mutations unique to individual cells that occur as part of normal aging prior to transformation. These mutations are targetable because ~10% of them form novel tumor-specific "NGG" protospacer adjacent motif (PAM) sites targetable by CRISPR-Cas9. Methods: Whole genome sequencing was performed on five rapid autopsy cases of patient-matched primary tumor, normal and metastatic tissue from pancreatic ductal adenocarcinoma decedents. CRISPR-Cas9 PAM targets were determined by bioinformatic tumor-normal subtraction for each patient and verified in metastatic samples by high-depth capture-based sequencing. Results: We found that 90% of PAM targets were maintained between primary carcinomas and metastases overall. We identified rules that predict PAM loss or retention, where PAMs located in heterozygous regions in the primary tumor can be lost in metastases (private LOH), but PAMs occurring in regions of loss of heterozygosity (LOH) in the primary tumor were universally conserved in metastases. Conclusions: Regions of truncal LOH are strongly retained in the presence of genetic instability, and therefore represent genetic vulnerabilities in pancreatic adenocarcinomas. A CRISPR-based gene therapy approach targeting these regions may be a novel way to genetically target metastatic cancer.

Convergent evolution of BRCA2 reversion mutations under therapeutic pressure by PARP inhibition and platinum chemotherapy.

Reversion mutations that restore wild-type function of the BRCA gene have been described as a key mechanism of resistance to Poly(ADP-ribose) polymerase (PARP) inhibitor therapy in BRCA-associated cancers. Here, we report a case of a patient with metastatic castration-resistant prostate cancer (mCRPC) with a germline BRCA2 mutation who developed acquired resistance to PARP inhibition. Extensive genomic interrogation of cell-free DNA (cfDNA) and tissue at baseline, post-progression, and postmortem revealed ten unique BRCA2 reversion mutations across ten sites. While several of the reversion mutations were private to a specific site, nine out of ten tumors contained at least one mutation, suggesting a powerful clonal selection for reversion mutations in the presence of therapeutic pressure by PARP inhibition. Variable cfDNA shed was seen across tumor sites, emphasizing a potential shortcoming of cfDNA monitoring for PARPi resistance. This report provides a genomic portrait of the temporal and spatial heterogeneity of prostate cancer under the selective pressure of a PARP inhibition and exposes limitations in the current strategies for detection of reversion mutations.

ConDoR: tumor phylogeny inference with a copy-number constrained mutation loss model.

A tumor contains a diverse collection of somatic mutations that reflect its past evolutionary history and that range in scale from single nucleotide variants (SNVs) to large-scale copy-number aberrations (CNAs). However, no current single-cell DNA sequencing (scDNA-seq) technology produces accurate measurements of both SNVs and CNAs, complicating the inference of tumor phylogenies. We introduce a new evolutionary model, the constrained k-Dollo model, that uses SNVs as phylogenetic markers but constrains losses of SNVs according to clusters of cells. We derive an algorithm, ConDoR, that infers phylogenies from targeted scDNA-seq data using this model. We demonstrate the advantages of ConDoR on simulated and real scDNA-seq data.

Subclonal Somatic Copy-Number Alterations Emerge and Dominate in Recurrent Osteosarcoma.

Multiple large-scale genomic profiling efforts have been undertaken in osteosarcoma to define the genomic drivers of tumorigenesis, therapeutic response, and disease recurrence. The spatial and temporal intratumor heterogeneity could also play a role in promoting tumor growth and treatment resistance. We conducted longitudinal whole-genome sequencing of 37 tumor samples from 8 patients with relapsed or refractory osteosarcoma. Each patient had at least one sample from a primary site and a metastatic or relapse site. Subclonal copy-number alterations were identified in all patients except one. In 5 patients, subclones from the primary tumor emerged and dominated at subsequent relapses. MYC gain/amplification was enriched in the treatment-resistant clones in 6 of 7 patients with multiple clones. Amplifications in other potential driver genes, such as CCNE1, RAD21, VEGFA, and IGF1R, were also observed in the resistant copy-number clones. A chromosomal duplication timing analysis revealed that complex genomic rearrangements typically occurred prior to diagnosis, supporting a macroevolutionary model of evolution, where a large number of genomic aberrations are acquired over a short period of time followed by clonal selection, as opposed to ongoing evolution. A mutational signature analysis of recurrent tumors revealed that homologous repair deficiency (HRD)-related SBS3 increases at each time point in patients with recurrent disease, suggesting that HRD continues to be an active mutagenic process after diagnosis. Overall, by examining the clonal relationships between temporally and spatially separated samples from patients with relapsed/refractory osteosarcoma, this study sheds light on the intratumor heterogeneity and potential drivers of treatment resistance in this disease. SIGNIFICANCE: The chemoresistant population in recurrent osteosarcoma is subclonal at diagnosis, emerges at the time of primary resection due to selective pressure from neoadjuvant chemotherapy, and is characterized by unique oncogenic amplifications.

The Ratio of Key Metabolic Transcripts Is a Predictive Biomarker of Breast Cancer Metastasis to the Lung.

Understanding the rewired metabolism underlying organ-specific metastasis in breast cancer could help identify strategies to improve the treatment and prevention of metastatic disease. Here, we used a systems biology approach to compare metabolic fluxes used by parental breast cancer cells and their brain- and lung-homing derivatives. Divergent lineages had distinct, heritable metabolic fluxes. Lung-homing cells maintained high glycolytic flux despite low levels of glycolytic intermediates, constitutively activating a pathway sink into lactate. This strong Warburg effect was associated with a high ratio of lactate dehydrogenase (LDH) to pyruvate dehydrogenase (PDH) expression, which correlated with lung metastasis in patients with breast cancer. Although feature classification models trained on clinical characteristics alone were unable to predict tropism, the LDH/PDH ratio was a significant predictor of metastasis to the lung but not to other organs, independent of other transcriptomic signatures. High lactate efflux was also a trait in lung-homing metastatic pancreatic cancer cells, suggesting that lactate production may be a convergent phenotype in lung metastasis. Together, these analyses highlight the essential role that metabolism plays in organ-specific cancer metastasis and identify a putative biomarker for predicting lung metastasis in patients with breast cancer. SIGNIFICANCE: Lung-homing metastatic breast cancer cells express an elevated ratio of lactate dehydrogenase to pyruvate dehydrogenase, indicating that ratios of specific metabolic gene transcripts have potential as metabolic biomarkers for predicting organ-specific metastasis.

Tissue-specific features of the T cell repertoire after allogeneic hematopoietic cell transplantation in human and mouse.

T cells are the central drivers of many inflammatory diseases, but the repertoire of tissue-resident T cells at sites of pathology in human organs remains poorly understood. We examined the site-specificity of T cell receptor (TCR) repertoires across tissues (5 to 18 tissues per patient) in prospectively collected autopsies of patients with and without graft-versus-host disease (GVHD), a potentially lethal tissue-targeting complication of allogeneic hematopoietic cell transplantation, and in mouse models of GVHD. Anatomic similarity between tissues was a key determinant of TCR repertoire composition within patients, independent of disease or transplant status. The T cells recovered from peripheral blood and spleens in patients and mice captured a limited portion of the TCR repertoire detected in tissues. Whereas few T cell clones were shared across patients, motif-based clustering revealed shared repertoire signatures across patients in a tissue-specific fashion. T cells at disease sites had a tissue-resident phenotype and were of donor origin based on single-cell chimerism analysis. These data demonstrate the complex composition of T cell populations that persist in human tissues at the end stage of an inflammatory disorder after lymphocyte-directed therapy. These findings also underscore the importance of studying T cell in tissues rather than blood for tissue-based pathologies and suggest the tissue-specific nature of both the endogenous and posttransplant T cell landscape.

Cancer cells co-evolve with retrotransposons to mitigate viral mimicry.

Overexpression of repetitive elements is an emerging hallmark of human cancers 1 . Diverse repeats can mimic viruses by replicating within the cancer genome through retrotransposition, or presenting pathogen-associated molecular patterns (PAMPs) to the pattern recognition receptors (PRRs) of the innate immune system 2-5 . Yet, how specific repeats affect tumor evolution and shape the tumor immune microenvironment (TME) in a pro- or anti-tumorigenic manner remains poorly defined. Here, we integrate whole genome and total transcriptome data from a unique autopsy cohort of multiregional samples collected in pancreatic ductal adenocarcinoma (PDAC) patients, into a comprehensive evolutionary analysis. We find that more recently evolved S hort I nterspersed N uclear E lements (SINE), a family of retrotransposable repeats, are more likely to form immunostimulatory double-strand RNAs (dsRNAs). Consequently, younger SINEs are strongly co-regulated with RIG-I like receptor associated type-I interferon genes but anti-correlated with pro-tumorigenic macrophage infiltration. We discover that immunostimulatory SINE expression in tumors is regulated by either L ong I nterspersed N uclear E lements 1 (LINE1/L1) mobility or ADAR1 activity in a TP53 mutation dependent manner. Moreover, L1 retrotransposition activity tracks with tumor evolution and is associated with TP53 mutation status. Altogether, our results suggest pancreatic tumors actively evolve to modulate immunogenic SINE stress and induce pro-tumorigenic inflammation. Our integrative, evolutionary analysis therefore illustrates, for the first time, how dark matter genomic repeats enable tumors to co-evolve with the TME by actively regulating viral mimicry to their selective advantage.

A targetable pathway to eliminate TRA-1-60+/TRA-1-81+ chemoresistant cancer cells.

Chemoresistance is a primary cause of treatment failure in pancreatic cancer. Identifying cell surface markers specifically expressed in chemoresistant cancer cells (CCCs) could facilitate targeted therapies to overcome chemoresistance. We performed an antibody-based screen and found that TRA-1-60 and TRA-1-81, two 'stemness' cell surface markers, are highly enriched in CCCs. Furthermore, TRA-1-60+/TRA-1-81+ cells are chemoresistant compared to TRA-1-60-/TRA-1-81- cells. Transcriptome profiling identified UGT1A10, shown to be both necessary and sufficient to maintain TRA-1-60/TRA-1-81 expression and chemoresistance. From a high-content chemical screen, we identified Cymarin, which downregulates UGT1A10, eliminates TRA-1-60/TRA-1-81 expression, and increases chemosensitivity both in vitro and in vivo. Finally, TRA-1-60/TRA-1-81 expression is highly specific in primary cancer tissue and positively correlated with chemoresistance and short survival, which highlights their potentiality for targeted therapy. Therefore, we discovered a novel CCC surface marker regulated by a pathway that promotes chemoresistance, as well as a leading drug candidate to target this pathway.

Clonal evolution during metastatic spread in high-risk neuroblastoma.

Patients with high-risk neuroblastoma generally present with widely metastatic disease and often relapse despite intensive therapy. As most studies to date focused on diagnosis-relapse pairs, our understanding of the genetic and clonal dynamics of metastatic spread and disease progression remain limited. Here, using genomic profiling of 470 sequential and spatially separated samples from 283 patients, we characterize subtype-specific genetic evolutionary trajectories from diagnosis through progression and end-stage metastatic disease. Clonal tracing timed disease initiation to embryogenesis. Continuous acquisition of structural variants at disease-defining loci (MYCN, TERT, MDM2-CDK4) followed by convergent evolution of mutations targeting shared pathways emerged as the predominant feature of progression. At diagnosis metastatic clones were already established at distant sites where they could stay dormant, only to cause relapses years later and spread via metastasis-to-metastasis and polyclonal seeding after therapy.

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.

Application of high-throughput single-nucleus DNA sequencing in pancreatic cancer.

Despite insights gained by bulk DNA sequencing of cancer it remains challenging to resolve the admixture of normal and tumor cells, and/or of distinct tumor subclones; high-throughput single-cell DNA sequencing circumvents these and brings cancer genomic studies to higher resolution. However, its application has been limited to liquid tumors or a small batch of solid tumors, mainly because of the lack of a scalable workflow to process solid tumor samples. Here we optimize a highly automated nuclei extraction workflow that achieves fast and reliable targeted single-nucleus DNA library preparation of 38 samples from 16 pancreatic ductal adenocarcinoma patients, with an average library yield per sample of 2867 single nuclei. We demonstrate that this workflow not only performs well using low cellularity or low tumor purity samples but reveals genomic evolution patterns of pancreatic ductal adenocarcinoma as well.

Genomic Biomarkers Associated with Response to Induction Chemotherapy in Patients with Localized Pancreatic Ductal Adenocarcinoma.

PURPOSE: There is increasing use of neoadjuvant chemotherapy in the management of localized pancreatic ductal adenocarcinoma (PDAC), yet there are few validated biomarkers to guide therapy selection. We aimed to determine whether somatic genomic biomarkers predict response to induction FOLFIRINOX or gemcitabine/nab-paclitaxel. EXPERIMENTAL DESIGN: This single-institution cohort study included consecutive patients (N = 322) with localized PDAC (2011-2020) who received at least one cycle of FOLFIRINOX (N = 271) or gemcitabine/nab-paclitaxel (N = 51) as initial treatment. We assessed somatic alterations in four driver genes (KRAS, TP53, CDKN2A, and SMAD4) by targeted next-generation sequencing, and determined associations between these alterations and (1) rate of metastatic progression during induction chemotherapy, (2) surgical resection, and (3) complete/major pathologic response. RESULTS: The alteration rates in driver genes KRAS, TP53, CDKN2A, and SMAD4 were 87.0%, 65.5%, 26.7%, and 19.9%, respectively. For patients receiving first-line FOLFIRINOX, SMAD4 alterations were uniquely associated with metastatic progression (30.0% vs. 14.5%; P = 0.009) and decreased rate of surgical resection (37.1% vs. 66.7%; P < 0.001). For patients receiving induction gemcitabine/nab-paclitaxel, alterations in SMAD4 were not associated with metastatic progression (14.3% vs. 16.2%; P = 0.866) nor decreased rate of surgical resection (33.3% vs. 41.9%; P = 0.605). Major pathologic response was rare (6.3%) and not associated with type of chemotherapy regimen. CONCLUSIONS: SMAD4 alterations were associated with more frequent development of metastasis and lower probability of reaching surgical resection during neoadjuvant FOLFIRINOX but not gemcitabine/nab-paclitaxel. Confirmation in a larger, diverse patient cohort will be important before prospective evaluation of SMAD4 as a genomic biomarker to guide treatment selection.

Utility of promoter hypermethylation in malignant risk stratification of intraductal papillary mucinous neoplasms.

BACKGROUND: Intraductal papillary mucinous neoplasms (IPMNs), a type of cystic pancreatic cancer (PC) precursors, are increasingly identified on cross-sectional imaging and present a significant diagnostic challenge. While surgical resection of IPMN-related advanced neoplasia, i.e., IPMN-related high-grade dysplasia or PC, is an essential early PC detection strategy, resection is not recommended for IPMN-low-grade dysplasia (LGD) due to minimal risk of carcinogenesis, and significant procedural risks. Based on their promising results in prior validation studies targeting early detection of classical PC, DNA hypermethylation-based markers may serve as a biomarker for malignant risk stratification of IPMNs. This study investigates our DNA methylation-based PC biomarker panel (ADAMTS1, BNC1, and CACNA1G genes) in differentiating IPMN-advanced neoplasia from IPMN-LGDs. METHODS: Our previously described genome-wide pharmaco-epigenetic method identified multiple genes as potential targets for PC detection. The combination was further optimized and validated for early detection of classical PC in previous case-control studies. These promising genes were evaluated among micro-dissected IPMN tissue (IPMN-LGD: 35, IPMN-advanced neoplasia: 35) through Methylation-Specific PCR. The discriminant capacity of individual and combination of genes were delineated through Receiver Operating Characteristics curve analysis. RESULTS: As compared to IPMN-LGDs, IPMN-advanced neoplasia had higher hypermethylation frequency of candidate genes: ADAMTS1 (60% vs. 14%), BNC1 (66% vs. 3%), and CACGNA1G (25% vs. 0%). We observed Area Under Curve (AUC) values of 0.73 for ADAMTS1, 0.81 for BNC1, and 0.63 for CACNA1G genes. The combination of the BNC1/ CACNA1G genes resulted in an AUC of 0.84, sensitivity of 71%, and specificity of 97%. Combining the methylation status of the BNC1/CACNA1G genes, blood-based CA19-9, and IPMN lesion size enhanced the AUC to 0.92. CONCLUSION: DNA-methylation based biomarkers have shown a high diagnostic specificity and moderate sensitivity for differentiating IPMN-advanced neoplasia from LGDs. Addition of specific methylation targets can improve the accuracy of the methylation biomarker panel and enable the development of noninvasive IPMN stratification biomarkers.