Integrative molecular and spatial analysis reveals evolutionary dynamics and tumor-immune interplay of in situ and invasive acral melanoma.

In acral melanoma (AM), progression from in situ (AMis) to invasive AM (iAM) leads to significantly reduced survival. However, evolutionary dynamics during this process remain elusive. Here, we report integrative molecular and spatial characterization of 147 AMs using genomics, bulk and single-cell transcriptomics, and spatial transcriptomics and proteomics. Vertical invasion from AMis to iAM displays an early and monoclonal seeding pattern. The subsequent regional expansion of iAM exhibits two distinct patterns, clonal expansion and subclonal diversification. Notably, molecular subtyping reveals an aggressive iAM subset featured with subclonal diversification, increased epithelial-mesenchymal transition (EMT), and spatial enrichment of APOE+/CD163+ macrophages. In vitro and ex vivo experiments further demonstrate that APOE+CD163+ macrophages promote tumor EMT via IGF1-IGF1R interaction. Adnexal involvement can predict AMis with higher invasive potential whereas APOE and CD163 serve as prognostic biomarkers for iAM. Altogether, our results provide implications for the early detection and treatment of AM.

Deep whole-genome analysis of 494 hepatocellular carcinomas.

Over half of hepatocellular carcinoma (HCC) cases diagnosed worldwide are in China1-3. However, whole-genome analysis of hepatitis B virus (HBV)-associated HCC in Chinese individuals is limited4-8, with current analyses of HCC mainly from non-HBV-enriched populations9,10. Here we initiated the Chinese Liver Cancer Atlas (CLCA) project and performed deep whole-genome sequencing (average depth, 120×) of 494 HCC tumours. We identified 6 coding and 28 non-coding previously undescribed driver candidates. Five previously undescribed mutational signatures were found, including aristolochic-acid-associated indel and doublet base signatures, and a single-base-substitution signature that we termed SBS_H8. Pentanucleotide context analysis and experimental validation confirmed that SBS_H8 was distinct to the aristolochic-acid-associated SBS22. Notably, HBV integrations could take the form of extrachromosomal circular DNA, resulting in elevated copy numbers and gene expression. Our high-depth data also enabled us to characterize subclonal clustered alterations, including chromothripsis, chromoplexy and kataegis, suggesting that these catastrophic events could also occur in late stages of hepatocarcinogenesis. Pathway analysis of all classes of alterations further linked non-coding mutations to dysregulation of liver metabolism. Finally, we performed in vitro and in vivo assays to show that fibrinogen alpha chain (FGA), determined as both a candidate coding and non-coding driver, regulates HCC progression and metastasis. Our CLCA study depicts a detailed genomic landscape and evolutionary history of HCC in Chinese individuals, providing important clinical implications.

A TGF-ß-dominant chemoresistant phenotype of hepatoblastoma associated with aflatoxin exposure in children.

BACKGROUND AND AIMS: Hepatoblastoma (HB) is the most common liver cancer in children, posing a serious threat to children's health. Chemoresistance is the leading cause of mortality in patients with HB. A more explicit definition of the features of chemotherapy resistance in HB represents a fundamental urgent need. APPROACH AND RESULTS: We performed an integrative analysis including single-cell RNA sequencing, whole-exome sequencing, and bulk RNA sequencing in 180 HB samples, to reveal genomic features, transcriptomic profiles, and the immune microenvironment of HB. Multicolor immunohistochemistry staining and in vitro experiments were performed for validation. Here, we reported four HB transcriptional subtypes primarily defined by differential expression of transcription factors. Among them, the S2A subtype, characterized by strong expression of progenitor ( MYCN , MIXL1 ) and mesenchymal transcription factors ( TWIST1 , TBX5 ), was defined as a new chemoresistant subtype. The S2A subtype showed increased TGF-ß cancer-associated fibroblast and an immunosuppressive microenvironment induced by the upregulated TGF-ß of HB. Interestingly, the S2A subtype enriched SBS24 signature and significantly higher serum aflatoxin B1-albumin (AFB1-ALB) level in comparison with other subtypes. Functional assays indicated that aflatoxin promotes HB to upregulate TGF-ß. Furthermore, clinical prognostic analysis showed that serum AFB1-ALB is a potential indicator of HB chemoresistance and prognosis. CONCLUSIONS: Our studies offer new insights into the relationship between aflatoxin and HB chemoresistance and provide important implications for its diagnosis and treatment.

SCAR: Single-cell and Spatially-resolved Cancer Resources.

Advances in sequencing and imaging technologies offer a unique opportunity to unravel cell heterogeneity and develop new immunotherapy strategies for cancer research. There is an urgent need for a resource that effectively integrates a vast amount of transcriptomic profiling data to comprehensively explore cancer tissue heterogeneity and the tumor microenvironment. In this context, we developed the Single-cell and Spatially-resolved Cancer Resources (SCAR) database, a combined tumor spatial and single-cell transcriptomic platform, which is freely accessible at http://8.142.154.29/SCAR2023 or http://scaratlas.com. SCAR contains spatial transcriptomic data from 21 tumor tissues and single-cell transcriptomic data from 11 301 352 cells encompassing 395 cancer subtypes and covering a wide variety of tissues, organoids, and cell lines. This resource offers diverse functional modules to address key cancer research questions at multiple levels, including the screening of tumor cell types, metabolic features, cell communication and gene expression patterns within the tumor microenvironment. Moreover, SCAR enables the analysis of biomarker expression patterns and cell developmental trajectories. SCAR also provides a comprehensive analysis of multi-dimensional datasets based on 34 state-of-the-art omics techniques, serving as an essential tool for in-depth mining and understanding of cell heterogeneity and spatial location. The implications of this resource extend to both cancer biology research and cancer immunotherapy development.

Th1 cells inducing IFNγ response improves immunotherapy efficacy in gastric cancer.

Objective: Cancer immunotherapy has made remarkable advances in recent years, but its effectiveness in treating gastric cancer is often limited by the complexity of the tumor microenvironment and the lack of effective biomarkers. This study aimed to identify effective biomarkers for immunotherapy treatment by characterizing the tumor microenvironment. Methods: We retrieved the RNA-seq data from gastric cancer patients treated with the programmed death 1 (PD-1) blockade pembrolizumab. Differentially expressed genes associated with clinical outcomes were identified and further analyzed using gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis. Gene signature scores were calculated by single sample Gene Set Enrichment Analysis (ssGSEA). The infiltration levels of immune cells were quantified using the xCell website. Cell type enrichment analysis was performed to compare treatment response and non-response groups, and regression analysis was used to investigate the relationship between interferon gamma (IFNγ) immune response and immune cell infiltration. Biomarkers were identified using least absolute shrinkage and selection operator (LASSO) analysis. Results: Compared to normal tissues, cytokine activity and interleukin-6 production were highly activated in gastric tumors. Responders to pembrolizumab showed significantly up-regulated expression of IFNγ response-related genes. Cell type enrichment analysis revealed that Th1 cells were significantly enriched in the tumor microenvironment of responders. Regression analysis indicated that Th1 cells induced IFNγ response more efficiently than other cell types. Using signatures of Th1 cells, stromal cells and IFNγ response, a set of eight genes were identified that effectively predicted the efficacy of immunotherapy treatment and patient prognosis. Conclusions: Th1 cells promote therapeutic efficacy of PD-1 blockade by promoting IFNγ immune response in gastric cancer. The identified biomarkers have the potential to improve the effectiveness of immunotherapy treatment for gastric cancer patients.

Intratumoral erythroblastic islands restrain anti-tumor immunity in hepatoblastoma.

Erythroblastic islands (EBIs) are the specialized structures for erythropoiesis, but they have never been found functional in tumors. As the most common pediatric liver malignancy, hepatoblastoma (HB) requires more effective and safer therapies to prevent progression and the lifelong impact of complications on young children. However, developing such therapies is impeded by a lack of comprehensive understanding of the tumor microenvironment. By single-cell RNA sequencing of 13 treatment-naive HB patients, we discover an immune landscape characterized by aberrant accumulation of EBIs, formed by VCAM1+ macrophages and erythroid cells, which is inversely correlated with survival of HB. Erythroid cells inhibit the function of dendritic cells (DCs) via the LGALS9/TIM3 axis, leading to impaired anti-tumor T cell immune responses. Encouragingly, TIM3 blockades relieve the inhibitory effect of erythroid cells on DCs. Our study provides an immune evasion mechanism mediated by intratumoral EBIs and proposes TIM3 as a promising therapeutic target for HB.

Liver tumour immune microenvironment subtypes and neutrophil heterogeneity.

The heterogeneity of the tumour immune microenvironment (TIME), organized by various immune and stromal cells, is a major contributing factor of tumour metastasis, relapse and drug resistance1-3, but how different TIME subtypes are connected to the clinical relevance in liver cancer remains unclear. Here we performed single-cell RNA-sequencing (scRNA-seq) analysis of 189 samples collected from 124 patients and 8 mice with liver cancer. With more than 1 million cells analysed, we stratified patients into five TIME subtypes, including immune activation, immune suppression mediated by myeloid or stromal cells, immune exclusion and immune residence phenotypes. Different TIME subtypes were spatially organized and associated with chemokine networks and genomic features. Notably, tumour-associated neutrophil (TAN) populations enriched in the myeloid-cell-enriched subtype were associated with an unfavourable prognosis. Through in vitro induction of TANs and ex vivo analyses of patient TANs, we showed that CCL4+ TANs can recruit macrophages and that PD-L1+ TANs can suppress T cell cytotoxicity. Furthermore, scRNA-seq analysis of mouse neutrophil subsets revealed that they are largely conserved with those of humans. In vivo neutrophil depletion in mouse models attenuated tumour progression, confirming the pro-tumour phenotypes of TANs. With this detailed cellular heterogeneity landscape of liver cancer, our study illustrates diverse TIME subtypes, highlights immunosuppressive functions of TANs and sheds light on potential immunotherapies targeting TANs.

Single-Cell DNA Sequencing Reveals Punctuated and Gradual Clonal Evolution in Hepatocellular Carcinoma.

BACKGROUND & AIMS: Copy number alterations (CNAs), elicited by genome instability, are a major source of intratumor heterogeneity. How CNAs evolve in hepatocellular carcinoma (HCC) remains unknown. METHODS: We performed single-cell DNA sequencing (scDNA-seq) on 1275 cells isolated from 10 patients with HCC, ploidy-resolved scDNA-seq on 356 cells from 1 additional patient, and single-cell RNA sequencing on 27,344 cells from 3 additional patients. Three statistical fitting models were compared to investigate the CNA accumulation pattern. RESULTS: Cells in the tumor were categorized into the following 3 subpopulations: euploid, pseudoeuploid, and aneuploid. Our scDNA-seq analysis revealed that CNA accumulation followed a dual-phase copy number evolution model, that is, a punctuated phase followed by a gradual phase. Patients who exhibited prolonged gradual phase showed higher intratumor heterogeneity and worse disease-free survival. Integrating bulk RNA sequencing of 17 patients with HCC, published datasets of 1196 liver tumors, and immunohistochemical staining of 202 HCC tumors, we found that high expression of CAD, a gene involved in pyrimidine synthesis, was correlated with rapid tumorigenesis and reduced survival. The dual-phase copy number evolution model was validated by our single-cell RNA sequencing data and published scDNA-seq datasets of other cancer types. Furthermore, ploidy-resolved scDNA-seq revealed the common clonal origin of diploid- and polyploid-aneuploid cells, suggesting that polyploid tumor cells were generated by whole genome doubling of diploid tumor cells. CONCLUSIONS: Our work revealed a novel dual-phase copy number evolution model, showed HCC with longer gradual phase was more severe, identified CAD as a promising biomarker for early recurrence of HCC, and supported the diploid origin of polyploid HCC.

Genomic and transcriptomic profiling of hepatoid adenocarcinoma of the stomach.

Hepatoid adenocarcinoma of the stomach (HAS), a rare subtype of gastric cancer (GC), has a low incidence but a high mortality rate. Little is known about the molecular features of HAS. Here we applied whole-exome sequencing (WES) on 58 tumours and the matched normal controls from 54 HAS patients, transcriptome sequencing on 30 HAS tumours, and single-cell RNA sequencing (scRNA-seq) on one HAS tumour. Our results reveal that the adenocarcinomatous component and hepatocellular-like component of the same HAS tumour originate monoclonally, and HAS is likely to initiate from pluripotent precursor cells. HAS has high stemness and high methionine cycle activity compared to classical GC. Two genes in the methionine cycle, MAT2A, and AHCY are potential targets for HAS treatments. We provide the first integrative genomic profiles of HAS, which may facilitate its diagnosis, prognosis, and treatment.

IDH Mutation Subgroup Status Associates with Intratumor Heterogeneity and the Tumor Microenvironment in Intrahepatic Cholangiocarcinoma.

Intrahepatic cholangiocarcinoma (ICC) is highly heterogeneous. Here, the authors perform exome sequencing and bulk RNA sequencing on 73 tumor regions from 14 ICC patients to portray the multi-faceted intratumor heterogeneity (ITH) landscape of ICC. The authors show that ITH is highly concordant across genomic, transcriptomic, and immune levels. Comparison of these data to 8 published datasets reveals significantly higher degrees of ITH in ICC than hepatocellular carcinoma. Remarkably, the authors find that high-ITH tumors highly overlap with the IDH (isocitrate dehydrogenase)-mutant subgroup (IDH-SG), comprising of IDH-mutated tumors and IDH-like tumors, that is, those IDH-wildtype tumors that exhibit similar molecular profiles to the IDH-mutated ones. Furthermore, IDH-SG exhibits less T cell infiltration and lower T cell cytotoxicity, indicating a colder tumor microenvironment (TME). The higher ITH and colder TME of IDH-SG are successfully validated by single-cell RNA sequencing on 17 503 cells from 4 patients. Collectively, the study shows that IDH mutant subgroup status, rather than IDH mutation alone, is associated with ITH and the TME of ICC tumors. The results highlight that IDH-like patients may also benefit from IDH targeted therapies and provide important implications for the diagnosis and treatment of ICC.

Mapping the spreading routes of lymphatic metastases in human colorectal cancer.

Lymphatic metastases are closely associated with tumor relapse and reduced survival in colorectal cancer (CRC). How tumor cells disseminate within the lymphatic network remains largely unknown. Here, we analyze the subclonal structure of 94 tumor samples, covering the primary tumors, lymph node metastases (LNMs), and liver metastases from 10 CRC patients. We portray a high-resolution lymphatic metastatic map for CRC by dividing LNMs into paracolic, intermediate, and central subgroups. Among the 61 metastatic routes identified, 38 (62.3%) are initiated from the primary tumors, 22 (36.1%) from LNMs, and 1 from liver metastasis (1.6%). In 5 patients, we find 6 LNMs that reseed 2 or more LNMs. We summarize 3 diverse modes of metastasis in CRC and show that skip spreading of tumor cells within the lymphatic network is common. Our study sheds light on the complicated metastatic pattern in CRC and has great clinical implications.

Somatic mutation landscape reveals differential variability of cell-of-origin for primary liver cancer.

Primary liver tissue cancer types are renowned to display a consistent increase in global disease burden and mortality, thus needing more effective diagnostics and treatments. Yet, integrative research efforts to identify cell-of-origin for these cancers by utilizing human specimen data were poorly established. To this end, we analyzed previously published whole-genome sequencing data for 384 tumor and progenitor tissues along with 423 publicly available normal tissue epigenomic features and single cell RNA-seq data from human livers to assess correlation patterns and extended this information to conduct in-silico prediction of the cell-of-origin for primary liver cancer subtypes. Despite mixed histological features, the cell-of-origin for mixed hepatocellular carcinoma/intrahepatic cholangiocarcinoma subtype was predominantly predicted to be hepatocytic origin. Individual sample-level predictions also revealed hepatocytes as one of the major predicted cell-of-origin for intrahepatic cholangiocarcinoma, thus implying trans-differentiation process during cancer progression. Additional analyses on the whole genome sequencing data of hepatic progenitor cells suggest these cells may not be a direct cell-of-origin for liver cancers. These results provide novel insights on the nature and potential contributors of cell-of-origins for primary liver cancers.

Genomic and Transcriptomic Profiling of Combined Hepatocellular and Intrahepatic Cholangiocarcinoma Reveals Distinct Molecular Subtypes.

We performed genomic and transcriptomic sequencing of 133 combined hepatocellular and intrahepatic cholangiocarcinoma (cHCC-ICC) cases, including separate, combined, and mixed subtypes. Integrative comparison of cHCC-ICC with hepatocellular carcinoma and intrahepatic cholangiocarcinoma revealed that combined and mixed type cHCC-ICCs are distinct subtypes with different clinical and molecular features. Integrating laser microdissection, cancer cell fraction analysis, and single nucleus sequencing, we revealed both mono- and multiclonal origins in the separate type cHCC-ICCs, whereas combined and mixed type cHCC-ICCs were all monoclonal origin. Notably, cHCC-ICCs showed significantly higher expression of Nestin, suggesting Nestin may serve as a biomarker for diagnosing cHCC-ICC. Our results provide important biological and clinical insights into cHCC-ICC.

Parity-time symmetry in periodically curved optical waveguides.

The coupling strength between two parity-time (PT) symmetric resonators determines whether the PT phase is broken or not. Here we investigate the scenario that two optical waveguides are spatially curved so that they switch periodically between unbroken and broken PT phases. We show that the existence of locally broken PT phase does not necessarily render a broken phase to waves propagating inside. Criteria are proposed to characterize the collective dynamics of wave near the Brillouin zone (BZ) edge, toward the cases of a totally broken phase, a partially broken phase, or a totally unbroken phase. We also discuss the characteristics of two special kinds of exceptional points (EPs) at the BZ edge, and show that their field patterns are displaced by half a period with each other. Full-wave numerical simulation proves our analysis. Potential applications especially these associated with EPs are discussed. This study helps us to understand how the locally PT-symmetric related eigenstate influences the globally collective dynamics of wave in spatially periodic configuration.

Unidirectional Excitation of Radiative-Loss-Free Surface Plasmon Polaritons in PT-Symmetric Systems.

We investigate the excitation and propagation of surface plasmon polaritons (SPPs) at a geometrically flat metal-dielectric interface with a parity-time (PT) symmetric modulation on the permittivity ϵ(x) of the dielectric medium. We show that two striking effects can be simultaneously achieved thanks to the nonreciprocal nature of the Bloch modes in the system. First, SPPs can be unidirectionally excited when light is normally incident on the interface. Secondly, the backscattering of SPPs into the far field is suppressed, producing a radiative-loss-free effect on the unidirectional SPPs. As a result, the lifetime and propagation distance of SPPs can be significantly improved. These results show that PT symmetry can be employed as a new approach to designing transformative nanoscale optical devices, such as low-loss plasmonic routers and isolators for efficient optical computation, communication, and information processing.

Coherent-trapped helical mode in parity-time symmetric metamaterials.

Coaxial optical subwavelength elements support helical modes Lm with different topological indexes m. Here we propose to couple the two bright L±1 modes with the dark one L0 via a parity-time (PT) symmetric perturbation. We show that the cascading coupled configuration is similar to a three-level atomic system, and supports a special hybridized mode Lc via a classic analog of coherent-population-trapping effect. Resonant frequency of Lc is independent of the PT-symmetric perturbation. Populations in L±1 can be manipulated by tuning the PT-symmetric perturbation, and no population is trapped in L0. Since the L±1 modes are associated with optical waves of opposite circular polarizations, the polarization of transmitted wave is independent of the polarization of incidence but solely determined by the PT-symmetric perturbation. Such an effect can be utilized to manipulate the polarization state of light. Numerical simulation in a well-designed coaxial metamaterial verifies our analysis.

Genomic signatures of pancreatic adenosquamous carcinoma (PASC).

In pancreatic cancer, pancreatic adenosquamous carcinoma (PASC) containing both ductal adenocarcinoma and squamous carcinoma in the same tumour represents ∼4% of the total incidence. To date, the genomic features of this mixed tumour are still unknown. We analysed the genomes of 17 PASCs and 34 pancreatic ductal adenocarcinomas (PDACs), and showed that PASC carried highly enriched TP53 mutations and 3p loss as compared with PDAC. We also showed that adenomatous and squamous components of PASC harboured comparable genomic alterations, suggesting that the two cellular components develop from the same progenitor cancer cells. Our study has updated genomic knowledge to help with understanding mixed cancers of the pancreas. Copyright © 2017 Pathological Society of Great Britain and Ireland. Published by John Wiley & Sons, Ltd.

Localization and oscillation of optical beams in Moiré lattices.

We study the propagation of optical beams in two-dimensional Moiré lattices, and demonstrate position-dependent beam dynamics when a quasi-Bragg condition is satisfied. We show that when the optical beam is incident to a peak of the lattice envelop, an optical Zitterbewegung is obtained. If the optical beam is incident to a node of the envelop, a field localization effect takes place. The localized beam oscillates with a much larger spatial period than that of the optical Zitterbewegung. Variation of the oscillation period versus the split in periods is discussed. The position-dependent beam dynamics are explained by the excitation of proper bandedge eigenmodes of the Moiré lattice, and can be engineered via tuning the periods of the two superimposed Bragg lattices.