Tracing Oncogene Rearrangements in the Mutational History of Lung Adenocarcinoma.

Mutational processes giving rise to lung adenocarcinomas (LADCs) in non-smokers remain elusive. We analyzed 138 LADC whole genomes, including 83 cases with minimal contribution of smoking-associated mutational signature. Genomic rearrangements were not correlated with smoking-associated mutations and frequently served as driver events of smoking-signature-low LADCs. Complex genomic rearrangements, including chromothripsis and chromoplexy, generated 74% of known fusion oncogenes, including EML4-ALK, CD74-ROS1, and KIF5B-RET. Unlike other collateral rearrangements, these fusion-oncogene-associated rearrangements were frequently copy-number-balanced, representing a genomic signature of early oncogenesis. Analysis of mutation timing revealed that fusions and point mutations of canonical oncogenes were often acquired in the early decades of life. During a long latency, cancer-related genes were disrupted or amplified by complex rearrangements. The genomic landscape was different between subgroups-EGFR-mutant LADCs had frequent whole-genome duplications with p53 mutations, whereas fusion-oncogene-driven LADCs had frequent SETD2 mutations. Our study highlights LADC oncogenesis driven by endogenous mutational processes.

Characterizing Mutational Signatures in Human Cancer Cell Lines Reveals Episodic APOBEC Mutagenesis.

Multiple signatures of somatic mutations have been identified in cancer genomes. Exome sequences of 1,001 human cancer cell lines and 577 xenografts revealed most common mutational signatures, indicating past activity of the underlying processes, usually in appropriate cancer types. To investigate ongoing patterns of mutational-signature generation, cell lines were cultured for extended periods and subsequently DNA sequenced. Signatures of discontinued exposures, including tobacco smoke and ultraviolet light, were not generated in vitro. Signatures of normal and defective DNA repair and replication continued to be generated at roughly stable mutation rates. Signatures of APOBEC cytidine deaminase DNA-editing exhibited substantial fluctuations in mutation rate over time with episodic bursts of mutations. The initiating factors for the bursts are unclear, although retrotransposon mobilization may contribute. The examined cell lines constitute a resource of live experimental models of mutational processes, which potentially retain patterns of activity and regulation operative in primary human cancers.

Tumor hypoxia represses γδ T cell-mediated antitumor immunity against brain tumors.

The anatomic location and immunologic characteristics of brain tumors result in strong lymphocyte suppression. Consequently, conventional immunotherapies targeting CD8 T cells are ineffective against brain tumors. Tumor cells escape immunosurveillance by various mechanisms and tumor cell metabolism can affect the metabolic states and functions of tumor-infiltrating lymphocytes. Here, we discovered that brain tumor cells had a particularly high demand for oxygen, which affected γδ T cell-mediated antitumor immune responses but not those of conventional T cells. Specifically, tumor hypoxia activated the γδ T cell protein kinase A pathway at a transcriptional level, resulting in repression of the activatory receptor NKG2D. Alleviating tumor hypoxia reinvigorated NKG2D expression and the antitumor function of γδ T cells. These results reveal a hypoxia-mediated mechanism through which brain tumors and γδ T cells interact and emphasize the importance of γδ T cells for antitumor immunity against brain tumors.

Widespread somatic L1 retrotransposition in normal colorectal epithelium.

Throughout an individual's lifetime, genomic alterations accumulate in somatic cells1-11. However, the mutational landscape induced by retrotransposition of long interspersed nuclear element-1 (L1), a widespread mobile element in the human genome12-14, is poorly understood in normal cells. Here we explored the whole-genome sequences of 899 single-cell clones established from three different cell types collected from 28 individuals. We identified 1,708 somatic L1 retrotransposition events that were enriched in colorectal epithelium and showed a positive relationship with age. Fingerprinting of source elements showed 34 retrotransposition-competent L1s. Multidimensional analysis demonstrated that (1) somatic L1 retrotranspositions occur from early embryogenesis at a substantial rate, (2) epigenetic on/off of a source element is preferentially determined in the early organogenesis stage, (3) retrotransposition-competent L1s with a lower population allele frequency have higher retrotransposition activity and (4) only a small fraction of L1 transcripts in the cytoplasm are finally retrotransposed in somatic cells. Analysis of matched cancers further suggested that somatic L1 retrotransposition rate is substantially increased during colorectal tumourigenesis. In summary, this study illustrates L1 retrotransposition-induced somatic mosaicism in normal cells and provides insights into the genomic and epigenomic regulation of transposable elements over the human lifetime.

Spatial genomics maps the structure, nature and evolution of cancer clones.

Genome sequencing of cancers often reveals mosaics of different subclones present in the same tumour1-3. Although these are believed to arise according to the principles of somatic evolution, the exact spatial growth patterns and underlying mechanisms remain elusive4,5. Here, to address this need, we developed a workflow that generates detailed quantitative maps of genetic subclone composition across whole-tumour sections. These provide the basis for studying clonal growth patterns, and the histological characteristics, microanatomy and microenvironmental composition of each clone. The approach rests on whole-genome sequencing, followed by highly multiplexed base-specific in situ sequencing, single-cell resolved transcriptomics and dedicated algorithms to link these layers. Applying the base-specific in situ sequencing workflow to eight tissue sections from two multifocal primary breast cancers revealed intricate subclonal growth patterns that were validated by microdissection. In a case of ductal carcinoma in situ, polyclonal neoplastic expansions occurred at the macroscopic scale but segregated within microanatomical structures. Across the stages of ductal carcinoma in situ, invasive cancer and lymph node metastasis, subclone territories are shown to exhibit distinct transcriptional and histological features and cellular microenvironments. These results provide examples of the benefits afforded by spatial genomics for deciphering the mechanisms underlying cancer evolution and microenvironmental ecology.

Clonal dynamics in early human embryogenesis inferred from somatic mutation.

Cellular dynamics and fate decision in early human embryogenesis remain largely unknown owing to the challenges of performing studies in human embryos1. Here, we explored whole-genomes of 334 single-cell colonies and targeted deep sequences of 379 bulk tissues obtained from various anatomical locations of seven recently deceased adult human donors. Using somatic mutations as an intrinsic barcode, we reconstructed early cellular phylogenies that demonstrate (1) an endogenous mutational rate that is higher in the first cell division but decreases to approximately one per cell per cell division later in life; (2) universal unequal contribution of early cells to embryo proper, resulting from early cellular bottlenecks that stochastically set aside epiblast cells within the embryo; (3) examples of varying degrees of early clonal imbalances between tissues on the left and right sides of the body, different germ layers and specific anatomical parts and organs; (4) emergence of a few ancestral cells that will substantially contribute to adult cell pools in blood and liver; and (5) presence of mitochondrial DNA heteroplasmy in the fertilized egg. Our approach also provides insights into the age-related mutational processes and loss of sex chromosomes in normal somatic cells. In sum, this study provides a foundation for future studies to complete cellular phylogenies in human embryogenesis.

Genomic features of renal cell carcinoma developed during end-stage renal disease and dialysis.

Patients with end-stage renal disease (ESRD) or receiving dialysis have a much higher risk for renal cell carcinoma (RCC), but carcinogenic mechanisms and genomic features remain little explored and undefined. This study's goal was to identify the genomic features of ESRD RCC and characterize them for associations with tumor histology and dialysis exposure. In this study, we obtained 33 RCCs, with various histological subtypes, that developed in ESRD patients receiving dialysis and performed whole-genome sequencing and transcriptome analyses. Driver events, copy-number alteration (CNA) analysis and mutational signature profiling were performed using an analysis pipeline that integrated data from germline and somatic SNVs, Indels and structural variants as well as CNAs, while transcriptome data were analyzed for differentially expressed genes and through gene set enrichment analysis. ESRD related clear cell RCCs' driver genes and mutations mirrored those in sporadic ccRCCs. Longer dialysis periods significantly correlated with a rare mutational signature SBS23, whose etiology is unknown, and increased mitochondrial copy number. All acquired cystic disease (ACD)-RCCs, which developed specifically in ESRD patients, showed chromosome 16q amplification. Gene expression analysis suggests similarity between certain ACD-RCCs and papillary RCCs and in TCGA papillary RCCs with chromosome 16 gain identified enrichment for genes related to DNA repair, as well as pathways related to reactive oxygen species, oxidative phosphorylation and targets of Myc. This analysis suggests that ESRD or dialysis could induce types of cellular stress that impact some specific types of genomic damage leading to oncogenesis.

Estimation of intrafamilial DNA contamination in family trio genome sequencing using deviation from Mendelian inheritance.

With the increasing number of sequencing projects involving families, quality control tools optimized for family genome sequencing are needed. However, accurately quantifying contamination in a DNA mixture is particularly difficult when genetically related family members are the sources. We developed TrioMix, a maximum likelihood estimation (MLE) framework based on Mendel's law of inheritance, to quantify DNA mixture between family members in genome sequencing data of parent-offspring trios. TrioMix can accurately deconvolute any intrafamilial DNA contamination, including parent-offspring, sibling-sibling, parent-parent, and even multiple familial sources. In addition, TrioMix can be applied to detect genomic abnormalities that deviate from Mendelian inheritance patterns, such as uniparental disomy (UPD) and chimerism. A genome-wide depth and variant allele frequency plot generated by TrioMix facilitates tracing the origin of Mendelian inheritance deviations. We showed that TrioMix could accurately deconvolute genomes in both simulated and real data sets.

Author Correction: Landscape of somatic mutations in 560 breast cancer whole-genome sequences.

In the Methods section of this Article, 'greater than' should have been 'less than' in the sentence 'Putative regions of clustered rearrangements were identified as having an average inter-rearrangement distance that was at least 10 times greater than the whole-genome average for the individual sample. '. The Article has not been corrected.

Clinical and molecular characteristics of Korean patients with Kabuki syndrome.

INTRODUCTION: Kabuki syndrome (KS) is a rare disorder characterized by typical facial features, skeletal anomalies, fetal fingertip pad persistence, postnatal growth retardation, and intellectual disabilities. Heterozygous variants of the KMT2D and KDM6A genes are major genetic causes of KS. This study aimed to report the clinical and genetic characteristics of KS. METHODS: This study included 28 Korean patients (14 boys and 14 girls) with KS through molecular genetic testing, including direct Sanger sequencing, whole-exome sequencing, or whole-genome sequencing. RESULTS: The median age at clinical diagnosis was 18.5 months (IQR 7-58 months), and the median follow-up duration was 80.5 months (IQR 48-112 months). Molecular genetic testing identified different pathogenic variants of the KMT2D (n = 23) and KDM6A (n = 3) genes, including 15 novel variants. Patients showed typical facial features (100%), such as long palpebral fissure and eversion of the lower eyelid; intellectual disability/developmental delay (96%); short stature (79%); and congenital cardiac anomalies (75%). Although 71% experienced failure to thrive in infancy, 54% of patients showed a tendency toward overweight/obesity in early childhood. Patients with KDM6A variants demonstrated severe genotype-phenotype correlation. CONCLUSION: This study enhances the understanding of the clinical and genetic characteristics of KS.

Ultra-Low Level Somatic Mutations and Structural Variations in Focal Cortical Dysplasia Type II.

OBJECTIVE: Brain somatic mutations in mTOR pathway genes are a major genetic etiology of focal cortical dysplasia type II (FCDII). Despite a greater ability to detect low-level somatic mutations in the brain by deep sequencing and analytics, about 40% of cases remain genetically unexplained. METHODS: We included 2 independent cohorts consisting of 21 patients with mutation-negative FCDII without apparent mutations on conventional deep sequencing of bulk brain. To find ultra-low level somatic variants or structural variants, we isolated cells exhibiting phosphorylation of the S6 ribosomal protein (p-S6) in frozen brain tissues using fluorescence-activated cell sorting (FACS). We then performed deep whole-genome sequencing (WGS; >90×) in p-S6+ cells in a cohort of 11 patients with mutation-negative. Then, we simplified the method to whole-genome amplification and target gene sequencing of p-S6+ cells in independent cohort of 10 patients with mutation-negative followed by low-read depth WGS (10×). RESULTS: We found that 28.6% (6 of 21) of mutation-negative FCDII carries ultra-low level somatic mutations (less than 0.2% of variant allele frequency [VAF]) in mTOR pathway genes. Our method showed ~34 times increase of the average mutational burden in FACS mediated enrichment of p-S6+ cells (average VAF = 5.84%) than in bulky brain tissues (average VAF = 0.17%). We found that 19% (4 of 21) carried germline structural variations in GATOR1 complex undetectable in whole exome or targeted gene sequencing. CONCLUSIONS: Our method facilitates the detection of ultra-low level somatic mutations, in specifically p-S6+ cells, and germline structural variations and increases the genetic diagnostic rate up to ~80% for the entire FCDII cohort. ANN NEUROL 2023;93:1082-1093.

Sox17 Deficiency Promotes Pulmonary Arterial Hypertension via HGF/c-Met Signaling.

BACKGROUND: In large-scale genomic studies, Sox17, an endothelial-specific transcription factor, has been suggested as a putative causal gene of pulmonary arterial hypertension (PAH); however, its role and molecular mechanisms remain to be elucidated. We investigated the functional impacts and acting mechanisms of impaired Sox17 (SRY-related HMG-box17) pathway in PAH and explored its potential as a therapeutic target. METHODS: In adult mice, Sox17 deletion in pulmonary endothelial cells (ECs) induced PAH under hypoxia with high penetrance and severity, but not under normoxia. RESULTS: Key features of PAH, such as hypermuscularization, EC hyperplasia, and inflammation in lung arterioles, right ventricular hypertrophy, and elevated pulmonary arterial pressure, persisted even after long rest in normoxia. Mechanistically, transcriptomic profiling predicted that the combination of Sox17 deficiency and hypoxia activated c-Met signaling in lung ECs. HGF (hepatocyte grow factor), a ligand of c-Met, was upregulated in Sox17-deficient lung ECs. Pharmacologic inhibition of HGF/c-Met signaling attenuated and reversed the features of PAH in both preventive and therapeutic settings. Similar to findings in animal models, Sox17 levels in lung ECs were repressed in 26.7% of PAH patients (4 of 15), while those were robust in all 14 non-PAH controls. HGF levels in pulmonary arterioles were increased in 86.7% of patients with PAH (13 of 15), but none of the controls showed that pattern. CONCLUSIONS: The downregulation of Sox17 levels in pulmonary arterioles increases the susceptibility to PAH, particularly when exposed to hypoxia. Our findings suggest the reactive upregulation of HGF/c-Met signaling as a novel druggable target for PAH treatment.

Human glioblastoma arises from subventricular zone cells with low-level driver mutations.

Glioblastoma (GBM) is a devastating and incurable brain tumour, with a median overall survival of fifteen months1,2. Identifying the cell of origin that harbours mutations that drive GBM could provide a fundamental basis for understanding disease progression and developing new treatments. Given that the accumulation of somatic mutations has been implicated in gliomagenesis, studies have suggested that neural stem cells (NSCs), with their self-renewal and proliferative capacities, in the subventricular zone (SVZ) of the adult human brain may be the cells from which GBM originates3-5. However, there is a lack of direct genetic evidence from human patients with GBM4,6-10. Here we describe direct molecular genetic evidence from patient brain tissue and genome-edited mouse models that show astrocyte-like NSCs in the SVZ to be the cell of origin that contains the driver mutations of human GBM. First, we performed deep sequencing of triple-matched tissues, consisting of (i) normal SVZ tissue away from the tumour mass, (ii) tumour tissue, and (iii) normal cortical tissue (or blood), from 28 patients with isocitrate dehydrogenase (IDH) wild-type GBM or other types of brain tumour. We found that normal SVZ tissue away from the tumour in 56.3% of patients with wild-type IDH GBM contained low-level GBM driver mutations (down to approximately 1% of the mutational burden) that were observed at high levels in their matching tumours. Moreover, by single-cell sequencing and laser microdissection analysis of patient brain tissue and genome editing of a mouse model, we found that astrocyte-like NSCs that carry driver mutations migrate from the SVZ and lead to the development of high-grade malignant gliomas in distant brain regions. Together, our results show that NSCs in human SVZ tissue are the cells of origin that contain the driver mutations of GBM.

Non-coding de novo mutations in chromatin interactions are implicated in autism spectrum disorder.

Three-dimensional chromatin interactions regulate gene expressions. The significance of de novo mutations (DNMs) in chromatin interactions remains poorly understood for autism spectrum disorder (ASD). We generated 813 whole-genome sequences from 242 Korean simplex families to detect DNMs, and identified target genes which were putatively affected by non-coding DNMs in chromatin interactions. Non-coding DNMs in chromatin interactions were significantly involved in transcriptional dysregulations related to ASD risk. Correspondingly, target genes showed spatiotemporal expressions relevant to ASD in developing brains and enrichment in biological pathways implicated in ASD, such as histone modification. Regarding clinical features of ASD, non-coding DNMs in chromatin interactions particularly contributed to low intelligence quotient levels in ASD probands. We further validated our findings using two replication cohorts, Simons Simplex Collection (SSC) and MSSNG, and showed the consistent enrichment of non-coding DNM-disrupted chromatin interactions in ASD probands. Generating human induced pluripotent stem cells in two ASD families, we were able to demonstrate that non-coding DNMs in chromatin interactions alter the expression of target genes at the stage of early neural development. Taken together, our findings indicate that non-coding DNMs in ASD probands lead to early neurodevelopmental disruption implicated in ASD risk via chromatin interactions.

Somatic mutations reveal asymmetric cellular dynamics in the early human embryo.

Somatic cells acquire mutations throughout the course of an individual's life. Mutations occurring early in embryogenesis are often present in a substantial proportion of, but not all, cells in postnatal humans and thus have particular characteristics and effects. Depending on their location in the genome and the proportion of cells they are present in, these mosaic mutations can cause a wide range of genetic disease syndromes and predispose carriers to cancer. They have a high chance of being transmitted to offspring as de novo germline mutations and, in principle, can provide insights into early human embryonic cell lineages and their contributions to adult tissues. Although it is known that gross chromosomal abnormalities are remarkably common in early human embryos, our understanding of early embryonic somatic mutations is very limited. Here we use whole-genome sequences of normal blood from 241 adults to identify 163 early embryonic mutations. We estimate that approximately three base substitution mutations occur per cell per cell-doubling event in early human embryogenesis and these are mainly attributable to two known mutational signatures. We used the mutations to reconstruct developmental lineages of adult cells and demonstrate that the two daughter cells of many early embryonic cell-doubling events contribute asymmetrically to adult blood at an approximately 2:1 ratio. This study therefore provides insights into the mutation rates, mutational processes and developmental outcomes of cell dynamics that operate during early human embryogenesis.

Dll4 Suppresses Transcytosis for Arterial Blood-Retinal Barrier Homeostasis.

RATIONALE: Central nervous system has low vascular permeability by organizing tight junction (TJ) and limiting endothelial transcytosis. While TJ has long been considered to be responsible for vascular barrier in central nervous system, suppressed transcytosis in endothelial cells is now emerging as a complementary mechanism. Whether transcytosis regulation is independent of TJ and its dysregulation dominantly causes diseases associated with edema remain elusive. Dll4 signaling is important for various vascular contexts, but its role in the maintenance of vascular barrier in central nervous system remains unknown. OBJECTIVE: To find a TJ-independent regulatory mechanism selective for transcytosis and identify its dysregulation as a cause of pathological leakage. METHODS AND RESULTS: We studied transcytosis in the adult mouse retina with low vascular permeability and employed a hypertension-induced retinal edema model for its pathological implication. Both antibody-based and genetic inactivation of Dll4 or Notch1 induce hyperpermeability by increasing transcytosis without junctional destabilization in arterial endothelial cells, leading to nonhemorrhagic leakage predominantly in the superficial retinal layer. Endothelial Sox17 deletion represses Dll4 in retinal arteries, phenocopying Dll4 blocking-driven vascular leakage. Ang II (angiotensin II)-induced hypertension represses arterial Sox17 and Dll4, followed by transcytosis-driven retinal edema, which is rescued by a gain of Notch activity. Transcriptomic profiling of retinal endothelial cells suggests that Dll4 blocking activates SREBP1 (sterol regulatory element-binding protein 1)-mediated lipogenic transcription and enriches gene sets favorable for caveolae formation. Profiling also predicts the activation of VEGF (vascular endothelial growth factor) signaling by Dll4 blockade. Inhibition of SREBP1 or VEGF-VEGFR2 (VEGF receptor 2) signaling attenuates both Dll4 blockade-driven and hypertension-induced retinal leakage. CONCLUSIONS: In the retina, Sox17-Dll4-SREBP1 signaling axis controls transcytosis independently of TJ in superficial arteries among heterogeneous regulations for the whole vessels. Uncontrolled transcytosis via dysregulated Dll4 underlies pathological leakage in hypertensive retina and could be a therapeutic target for treating hypertension-associated retinal edema.

Landscape of somatic mutations in 560 breast cancer whole-genome sequences.

We analysed whole-genome sequences of 560 breast cancers to advance understanding of the driver mutations conferring clonal advantage and the mutational processes generating somatic mutations. We found that 93 protein-coding cancer genes carried probable driver mutations. Some non-coding regions exhibited high mutation frequencies, but most have distinctive structural features probably causing elevated mutation rates and do not contain driver mutations. Mutational signature analysis was extended to genome rearrangements and revealed twelve base substitution and six rearrangement signatures. Three rearrangement signatures, characterized by tandem duplications or deletions, appear associated with defective homologous-recombination-based DNA repair: one with deficient BRCA1 function, another with deficient BRCA1 or BRCA2 function, the cause of the third is unknown. This analysis of all classes of somatic mutation across exons, introns and intergenic regions highlights the repertoire of cancer genes and mutational processes operating, and progresses towards a comprehensive account of the somatic genetic basis of breast cancer.

A fusion of CD63-BCAR4 identified in lung adenocarcinoma promotes tumorigenicity and metastasis.

BACKGROUND: Recently, fusion variants of the breast cancer anti-oestrogen-resistance 4 (BCAR4) gene were recurrently discovered in lung adenocarcinoma from the genome-wide studies. However, the functional characterisation of BCAR4 fusion has not been investigated. METHODS: Based on the analysis of RNA-sequencing data, we identified a fusion transcript of CD63-BCAR4 in a Korean patient with lung adenocarcinoma who did not harbour any known activating mutations in EGFR and KRAS genes. To investigate the oncogenic effect of CD63-BCAR4, in vitro and in vivo animal experiments were performed. RESULTS: In vitro experiments showed strongly enhanced cell migration and proliferation by the exogenous expression of CD63-BCAR4 protein in bronchial epithelial cells. Cell migration was notably reduced after knockdown of BCAR4 fusion by small-interfering RNA. The tumorigenic and metastatic capability of the CD63-BCAR4 fusion was confirmed by using the mouse xenograft model. Fusion-overexpressed cells result in metastasis to the liver and lung as well as the primary tumours after subcutaneous injection into mice. Cyclin D1, MMP1, Slug and mesenchymal markers were significantly increased after CD63-BCAR4 overexpression in the in vitro and in vivo experiments. CONCLUSIONS: Taken together, our results suggest a newly identified fusion gene, CD63-BCAR4 as a potential novel oncogene in lung adenocarcinoma.

4-1BB Delineates Distinct Activation Status of Exhausted Tumor-Infiltrating CD8+ T Cells in Hepatocellular Carcinoma.

BACKGROUND AND AIMS: Targeting costimulatory receptors with agonistic antibodies is a promising cancer immunotherapy option. We aimed to investigate costimulatory receptor expression, particularly 4-1BB (CD137 or tumor necrosis factor receptor superfamily member 9), on tumor-infiltrating CD8+ T cells (CD8+ tumor-infiltrating lymphocytes [TILs]) and its association with distinct T-cell activation features among exhausted CD8+ TILs in hepatocellular carcinoma (HCC). APPROACH AND RESULTS: Tumor tissues, adjacent nontumor tissues, and peripheral blood were collected from HCC patients undergoing surgical resection (n = 79). Lymphocytes were isolated and used for multicolor flow cytometry, RNA-sequencing, and in vitro functional restoration assays. Among the examined costimulatory receptors, 4-1BB was most prominently expressed on CD8+ TILs. 4-1BB expression was almost exclusively detected on CD8+ T cells in the tumor-especially on programmed death 1 (PD-1)high cells and not PD-1int and PD-1neg cells. Compared to PD-1int and 4-1BBneg PD-1high CD8+ TILs, 4-1BBpos PD-1high CD8+ TILs exhibited higher levels of tumor reactivity and T-cell activation markers and significant enrichment for T-cell activation gene signatures. Per-patient analysis revealed positive correlations between percentages of 4-1BBpos cells among CD8+ TILs and levels of parameters of tumor reactivity and T-cell activation. Among highly exhausted PD-1high CD8+ TILs, 4-1BBpos cells harbored higher proportions of cells with proliferative and reinvigoration potential. Our 4-1BB-related gene signature predicted survival outcomes of HCC patients in the The Cancer Genome Atlas cohort. 4-1BB agonistic antibodies enhanced the function of CD8+ TILs and further enhanced the anti-PD-1-mediated reinvigoration of CD8+ TILs, especially in cases showing high levels of T-cell activation. CONCLUSION: 4-1BB expression on CD8+ TILs represents a distinct activation state among highly exhausted CD8+ T cells in HCC. 4-1BB costimulation with agonistic antibodies may be a promising strategy for treating HCCs exhibiting prominent T-cell activation.

Implication of CD69+ CD103+ tissue-resident-like CD8+ T cells as a potential immunotherapeutic target for cholangiocarcinoma.

BACKGROUND: The heterogeneous immune landscapes of intrahepatic cholangiocarcinoma (ICC) remain largely unknown. Here we aimed to investigate the implications of tissue-resident memory (TRM)-related features of tumour-infiltrating CD8+ T cells (CD8+ TILs) from ICC patients. METHODS: From ICC patients, we obtained blood samples and ICC surgical specimens (n = 33). We performed multicolour flow cytometry, multiplexed immunohistochemistry and RNA sequencing. RESULTS: When compared to peripheral CD8+ T cells, the CD8+ TILs included significantly higher proportions of the CD69+ CD103- and CD69+ CD103+ TRM-like subsets (P < .001 for both). Relative to CD69- and CD69+ CD103- cells, the CD69+ CD103+ CD8+ TILs harboured higher levels of T-cell markers representing tumour specificity (ie CD39), proliferation (ie Ki-67) and T-cell activation (ie HLA-DR and CD38) (all P < .001). Moreover, compared to the stroma, the tumour margin and core density each had a significantly higher density of CD103+ CD8+ TILs (P < .001 for both). ICCs with high proportions of CD69+ CD103+ cells displayed higher levels of parameters associated with response to immune checkpoint inhibitors (ICIs)-including number of CD8+ TIL infiltrates (P = .019), PD-L1 expression in the tumour (P = .046) and expression of the T cell-inflamed gene signature (P < .001). ICCs with lower proportions of CD69+ CD103+ CD8+ TILs exhibited significant enrichment of genes related to the Wnt/ß-catenin (P < .001) and TGF-ß pathways (P = .002). CONCLUSION: CD69+ CD103+ TRM-like CD8+ TILs represent prominent tumour-specific immune responses and hold promise as a potential therapeutic target in ICC patients. Differential TRM-related features of ICCs may help develop future immunotherapeutic strategies such as maximizing TRM responses or inhibiting pathways contributing to immune evasion.