MesKit: a tool kit for dissecting cancer evolution of multi-region tumor biopsies through somatic alterations

Multiregion WES of metastatic pancreatic neuroendocrine tumors revealed heterogeneity in genomic alterations, immune microenvironment and evolutionary patterns

Pancreatic neuroendocrine tumors (PanNETs), though uncommon, have a high likelihood of spreading to other body parts. Previously, the genetic diversity and evolutionary patterns in metastatic PanNETs were not well understood. To investigate this, we performed multiregion sampling whole-exome sequencing (MRS-WES) on samples from 10 patients who had not received prior treatment for metastatic PanNETs. This included 29 primary tumor samples, 31 lymph node metastases, and 15 liver metastases. We used the MSK-MET dataset for survival analysis and validation of our findings. Our research indicates that mutations in the MEN1/DAXX genes might trigger the early stages of PanNET development. We categorized the patients based on the presence (MEN1/DAXXmut, n = 7) or absence (MEN1/DAXXwild, n = 3) of these mutations. Notable differences were observed between the two groups in terms of genetic alterations and clinically relevant mutations, confirmed using the MSK-MET dataset. Notably, patients with mutations in MEN1/DAXX/ATRX genes had a significantly longer median overall survival compared to those without these mutations (median not reached vs. 43.63 months, p = 0.047). Multiplex immunohistochemistry (mIHC) analysis showed a more prominent immunosuppressive environment in metastatic tumors, especially in patients with MEN1/DAXX mutations. These findings imply that MEN1/DAXX mutations lead PanNETs through a unique evolutionary path. The disease's progression pattern indicates that PanNETs can spread early, even before clinical detection, highlighting the importance of identifying biomarkers related to metastasis to guide personalized treatment strategies.

A functional personalised oncology approach against metastatic colorectal cancer in matched patient derived organoids

Globally, colorectal cancer (CRC) is the third most frequently occurring cancer. Progression on to an advanced metastatic malignancy (metCRC) is often indicative of poor prognosis, as the 5-year survival rates of patients decline rapidly. Despite the availability of many systemic therapies for the management of metCRC, the long-term efficacies of these regimens are often hindered by the emergence of treatment resistance due to intratumoral and intertumoral heterogeneity. Furthermore, not all systemic therapies have associated biomarkers that can accurately predict patient responses. Hence, a functional personalised oncology (FPO) approach can enable the identification of patient-specific combinatorial vulnerabilities and synergistic combinations as effective treatment strategies. To this end, we established a panel of CRC patient-derived organoids (PDOs) as clinically relevant biological systems, of which three pairs of matched metCRC PDOs were derived from the primary sites (ptCRC) and metastatic lesions (mCRC). Histological and genomic characterisation of these PDOs demonstrated the preservation of histopathological and genetic features found in the parental tumours. Subsequent application of the phenotypic-analytical drug combination interrogation platform, Quadratic Phenotypic Optimisation Platform, in these pairs of PDOs identified patient-specific drug sensitivity profiles to epigenetic-based combination therapies. Most notably, matched PDOs from one patient exhibited differential sensitivity patterns to the rationally designed drug combinations despite being genetically similar. These findings collectively highlight the limitations of current genomic-driven precision medicine in guiding treatment strategies for metCRC patients. Instead, it suggests that epigenomic profiling and application of FPO could complement the identification of novel combinatorial vulnerabilities to target synchronous ptCRC and mCRC.

Molecular cartography uncovers evolutionary and microenvironmental dynamics in sporadic colorectal tumors

Colorectal cancer exhibits dynamic cellular and genetic heterogeneity during progression from precursor lesions toward malignancy. Analysis of spatial multi-omic data from 31 human colorectal specimens enabled phylogeographic mapping of tumor evolution that revealed individualized progression trajectories and accompanying microenvironmental and clonal alterations. Phylogeographic mapping ordered genetic events, classified tumors by their evolutionary dynamics, and placed clonal regions along global pseudotemporal progression trajectories encompassing the chromosomal instability (CIN+) and hypermutated (HM) pathways. Integrated single-cell and spatial transcriptomic data revealed recurring epithelial programs and infiltrating immune states along progression pseudotime. We discovered an immune exclusion signature (IEX), consisting of extracellular matrix regulators DDR1, TGFBI, PAK4, and DPEP1, that charts with CIN+ tumor progression, is associated with reduced cytotoxic cell infiltration, and shows prognostic value in independent cohorts. This spatial multi-omic atlas provides insights into colorectal tumor-microenvironment co-evolution, serving as a resource for stratification and targeted treatments.

Genetic analysis of low-grade adenosquamous carcinoma of the breast progressing to high-grade metaplastic carcinoma

PURPOSE

Low-grade adenosquamous carcinoma (LGASC) is a rare type of metaplastic carcinoma of the breast (MBC) with an indolent clinical course. A few LGASC cases with high-grade transformation have been reported; however, the genetics underlying malignant progression of LGASC remain unclear.

METHODS

We performed whole-genome sequencing analysis on five MBCs from four patients, including one case with matching primary LGASC and a lymph node metastatic tumor consisting of high-grade MBC with a predominant metaplastic squamous cell carcinoma component (MSC) that progressed from LGASC and three cases of independent de novo MSC.

RESULTS

Unlike de novo MSC, LGASC and its associated MSC showed no TP53 mutation and tended to contain fewer structural variants than de novo MSC. Both LGASC and its associated MSC harbored the common GNAS c.C2530T:p.Arg844Cys mutation, which was more frequently detected in the cancer cell fraction of MSC. MSC associated with LGASC showed additional pathogenic deletions of multiple tumor-suppressor genes, such as KMT2D and BTG1. Copy number analysis revealed potential 18q loss of heterozygosity in both LGASC and associated MSC. The frequency of SMAD4::DCC fusion due to deletions increased with progression to MSC; however, chimeric proteins were not detected. SMAD4 protein expression was already decreased at the LGASC stage due to unknown mechanisms.

CONCLUSION

Not only LGASC but also its associated high-grade MBC may be genetically different from de novo high-grade MBC. Progression from LGASC to high-grade MBC may involve the concentration of driver mutations caused by clonal selection and inactivation of tumor-suppressor genes.

Strategic Approach to Heterogeneity Analysis of Cutaneous Adnexal Carcinomas Using Computational Pathology and Genomics

Cutaneous adnexal tumors are neoplasms that arise from skin appendages. Their morphologic diversity and phenotypic variability with rare progression to malignancy make them difficult to diagnose and classify, and there is currently no established treatment strategy. To overcome these difficulties, this study investigated the transcription factor SOX9 expression, morphology, and genetics of skin adnexal tumors for understanding their biology, especially their histogenesis. We showed that cutaneous adnexal tumors and their nontumor counterparts of skin and appendages exhibit expression patterns similar to that of SOX9. Its expression intensity and pattern, as well as histopathologic evaluation of tumors, were analyzed using digital images of 69 normal skin adnexal 9-type organs and 185 skin adnexal 29-type tumors as references. It was possible to distinguish basal cell carcinoma from squamous cell carcinoma, sebaceous carcinoma, and pilomatrixoma with significant differences, along with porocarcinoma from squamous cell carcinoma. Furthermore, unsupervised machine learning "computational pathology" was used to derive a multiregion whole-exome sequencing fusion method termed "genocomputed pathology." The genocomputed pathology of three representable adnexal carcinomas (porocarcinoma, hidradenocarcinoma, and spiradenocarcinoma) was evaluated for total nine cases. We showed that there was more heterogeneity than expected within the tumors as well as the coexistence of components lacking driver fusion genes. The presence or absence of potential driver genes, such as PIK3CA, YAP1, and PTEN, in each region was identified, highlighting a therapeutic strategy for cutaneous adnexal carcinoma encompassing heterogeneous tumors.

Precise microdissection of gastric mixed adeno-neuroendocrine carcinoma dissects its genomic landscape and evolutionary clonal origins

Gastric mixed adenoneuroendocrine carcinoma (MANEC) is a clinically aggressive and heterogeneous tumor composed of adenocarcinoma (ACA) and neuroendocrine carcinoma (NEC). The genomic properties and evolutionary clonal origins of MANEC remain unclear. We conduct whole-exome and multiregional sequencing on 101 samples from 33 patients to elucidate their evolutionary paths. We identify four significantly mutated genes, TP53, RB1, APC, and CTNNB1. MANEC resembles chromosomal instability stomach adenocarcinoma in that whole-genome doubling in MANEC is predominant and occurs earlier than most copy-number losses. All tumors are of monoclonal origin, and NEC components show more aggressive genomic properties than their ACA counterparts. The phylogenetic trees show two tumor divergence patterns, including sequential and parallel divergence. Furthermore, ACA-to-NEC rather than NEC-to-ACA transition is confirmed by immunohistochemistry on 6 biomarkers in ACA- and NEC-dominant regions. These results provide insights into the clonal origin and tumor differentiation of MANEC.

Genomic Patterns of Malignant Peripheral Nerve Sheath Tumor (MPNST) Evolution Correlate with Clinical Outcome and Are Detectable in Cell-Free DNA

Malignant peripheral nerve sheath tumor (MPNST), an aggressive soft-tissue sarcoma, occurs in people with neurofibromatosis type 1 (NF1) and sporadically. Whole-genome and multiregional exome sequencing, transcriptomic, and methylation profiling of 95 tumor samples revealed the order of genomic events in tumor evolution. Following biallelic inactivation of NF1, loss of CDKN2A or TP53 with or without inactivation of polycomb repressive complex 2 (PRC2) leads to extensive somatic copy-number aberrations (SCNA). Distinct pathways of tumor evolution are associated with inactivation of PRC2 genes and H3K27 trimethylation (H3K27me3) status. Tumors with H3K27me3 loss evolve through extensive chromosomal losses followed by whole-genome doubling and chromosome 8 amplification, and show lower levels of immune cell infiltration. Retention of H3K27me3 leads to extensive genomic instability, but an immune cell-rich phenotype. Specific SCNAs detected in both tumor samples and cell-free DNA (cfDNA) act as a surrogate for H3K27me3 loss and immune infiltration, and predict prognosis.

SIGNIFICANCE

MPNST is the most common cause of death and morbidity for individuals with NF1, a relatively common tumor predisposition syndrome. Our results suggest that somatic copy-number and methylation profiling of tumor or cfDNA could serve as a biomarker for early diagnosis and to stratify patients into prognostic and treatment-related subgroups. This article is highlighted in the In This Issue feature, p. 517.

Induction immune-checkpoint inhibitors for resectable oncogene-mutant NSCLC: A multicenter pooled analysis

Despite limited efficacy of immunotherapy for advanced non-small-cell lung cancer (NSCLC) with driver mutations, whether neoadjuvant immunotherapy could be clinically valuable in those patients warrants further investigation. We utilized 40 oncogene-mutant NSCLC treated with induction immunotherapy from a large consecutive multicenter cohort. Overall response rate was 62.5% while 2 patients had disease progression. Of 39 patients that received surgery, R0 resection rate was 97.4%. The major pathological response (MPR) rate was 37.5% and the pathological complete response (pCR) rate was 12.5%. Pre-treatment PD-L1 expression was not a predictive biomarker in these patients. Median disease-free survival for all oncogenic mutation and EGFR mutation was 28.5 months. Indirect comparison through integrating CTONG1103 cohort showed neoadjuvant immunotherapy plus chemotherapy yielded the most superior efficacy among erlotinib and chemotherapy for resectable EGFR-mutant NSCLC. No MPR patients were identified with neoadjuvant immunotherapy plus chemotherapy for uncommon EGFR insertion or point mutations. Our results indicated the potential clinical feasibility of neoadjuvant immunotherapy for resectable localized oncogene-mutant NSCLC especially for EGFR-mutant NSCLC.

Deciphering clonal dynamics and metastatic routines in a rare patient of synchronous triple-primary tumors and multiple metastases with MPTevol

Multiple primary tumor (MPT) is a special and rare cancer type, defined as more than two primary tumors presenting at the diagnosis in a single patient. The molecular characteristics and tumorigenesis of MPT remain unclear due to insufficient approaches. Here, we present MPTevol, a practical computational framework for comprehensively exploring the MPT from multiregion sequencing (MRS) experiments. To verify the utility of MPTevol, we performed whole-exome MRS for 33 samples of a rare patient with triple-primary tumors and three metastatic sites and systematically investigated clonal dynamics and metastatic routines. MPTevol assists in comparing genomic profiles across samples, detecting clonal evolutionary history and metastatic routines and quantifying the metastatic history. All triple-primary tumors were independent origins and their genomic characteristics were consistent with corresponding sporadic tumors, strongly supporting their independent tumorigenesis. We further showed two independent early monoclonal seeding events for the metastases in the ovary and uterus. We revealed that two ovarian metastases were disseminated from the same subclone of the primary tumor through undergoing whole-genome doubling processes, suggesting metastases-to-metastases seeding occurred when tumors had similar microenvironments. Surprisingly, according to the metastasis timing model of MPTevol, we found that primary tumors of about 0.058–0.124 cm diameter have been disseminating to distant organs, which is much earlier than conventional clinical views. We developed MPT-specialized analysis framework MPTevol and demonstrated its utility in explicitly resolving clonal evolutionary history and metastatic seeding routines with a rare MPT case. MPTevol is implemented in R and is available at https://github.com/qingjian1991/MPTevol under the GPL v3 license.

Characterization of Genetic Heterogeneity in Recurrent Metastases of Renal Cell Carcinoma

Metastatic renal cell carcinoma (RCC) exhibits poor prognosis. Better knowledge of distant metastases is crucial to foster personalized treatment strategies. Here, we aimed to investigate the genetic landscape of metastases, including synchronous and/or recurrent metastases to elucidate potential drug target genes and clinically relevant mutations in a real-world setting of patients. We assessed 81 metastases from 56 RCC patients, including synchronous and/or recurrent metastases of 19 patients. Samples were analysed through next-generation sequencing with a high coverage (~1000× mean coverage). We therefore established a novel sequencing panel comprising 32 genes with impact on RCC development. We observed a high frequency of mutations in known RCC driver genes (e.g., >40% carriers of VHL and PBRM1 mutations) in metastases irrespective of the metastatic site. The somatic mutational composition was significantly associated with cancer-specific survival (p(logrank) = 0.03). Moreover, we identified in 34 patients at least one drug target gene as well as clinically relevant mutations listed in the VICC Meta-Knowledgebase in 7%. In addition to significantly higher mutational burden in recurrent metastases compared to earlier ones, synchronous and/or recurrent metastases of individual patients, even after a time-period >2 yrs, shared a high proportion of somatic events. Our data demonstrate the importance of somatic profiling in metastases for precision medicine in RCC.

MesKit: a tool kit for dissecting cancer evolution of multi-region tumor biopsies through somatic alterations

Background

Multi-region sequencing (MRS) has been widely used to analyze intra-tumor heterogeneity (ITH) and cancer evolution. However, comprehensive analysis of mutational data from MRS is still challenging, necessitating complicated integration of a plethora of computational and statistical approaches.

Findings

Here, we present MesKit, an R/Bioconductor package that can assist in characterizing genetic ITH and tracing the evolutionary history of tumors based on somatic alterations detected by MRS. MesKit provides a wide range of analysis and visualization modules, including ITH evaluation, metastatic route inference, and mutational signature identification. In addition, MesKit implements an auto-layout algorithm to generate phylogenetic trees based on somatic mutations. The application of MesKit for 2 reported MRS datasets of hepatocellular carcinoma and colorectal cancer identified known heterogeneous features and evolutionary patterns, together with potential driver events during cancer evolution.

Conclusions

In summary, MesKit is useful for interpreting ITH and tracing evolutionary trajectory based on MRS data. MesKit is implemented in R and available at https://bioconductor.org/packages/MesKit under the GPL v3 license.