Citations

We describe the design and results from the BraTS 2023 Intracranial Meningioma Segmentation Challenge. The BraTS Meningioma Challenge differed from prior BraTS Glioma challenges in that it focused on meningiomas, which are typically benign extra-axial tumors with diverse radiologic and anatomical presentation and a propensity for multiplicity. Nine participating teams each developed deep-learning automated segmentation models using image data from the largest multi-institutional systematically expert annotated multilabel multi-sequence meningioma MRI dataset to date, which included 1000 training set cases, 141 validation set cases, and 283 hidden test set cases. Each case included T2, FLAIR, T1, and T1Gd brain MRI sequences with associated tumor compartment labels delineating enhancing tumor, non-enhancing tumor, and surrounding non-enhancing FLAIR hyperintensity. Participant automated segmentation models were evaluated and ranked based on a scoring system evaluating lesion-wise metrics including dice similarity coefficient (DSC) and 95% Hausdorff Distance. The top ranked team had a lesion-wise median dice similarity coefficient (DSC) of 0.976, 0.976, and 0.964 for enhancing tumor, tumor core, and whole tumor, respectively and a corresponding average DSC of 0.899, 0.904, and 0.871, respectively. These results serve as state-of-the-art benchmarks for future pre-operative meningioma automated segmentation algorithms. Additionally, we found that 1286 of 1424 cases (90.3%) had at least 1 compartment voxel abutting the edge of the skull-stripped image edge, which requires further investigation into optimal pre-processing face anonymization steps.

Abstract Background Inter‐fractional and intra‐fractional positioning and motion monitoring based on implanted fiducials through kilovoltage (kV) imaging is a cost‐effective approach to enhance intact prostate radiotherapy treatment accuracy. However, comprehensive studies for three‐dimensional (3D) position extraction and the corrective action threshold are lacking. Purpose To develop and verify a fiducial motion monitoring system with 3D position information based on triggered kV images taken during treatment delivery, and to comprehensively evaluate the recommended thresholds for corrective action during treatment. Methods An in‐house fiducial triangulation algorithm was developed to monitor fiducial positions using subsequent kV triggered images. The precision of fiducial triangulation and motion detection was validated on a pelvis phantom with four fiducial inserts. A retrospective study was conducted on prostate cancer patients who received either moderately hypofractionated radiotherapy (MHRT, n = 16) or stereotactic body radiotherapy (SBRT, n = 12), categorized by endo‐rectal balloon (ERB) use. Intra‐fractional positions relative to isocenter were computed, and these positions were analyzed in comparison to the original marker locations of the patient's treatment plan. A linear regression fit was used per fraction to determine correlation coefficients between motions in the left–right (LR), superior–inferior (SI), and anterior–posterior (AP) directions. Average fraction time (AFT) of treatment was reported based on the observed average time spent on steps in the workflow and the percentage of fiducials caught outside of the current two‐dimensional (2D) tolerance threshold of 5 mm. The observed 3D offsets were used to estimate AFT for various 3D tolerance thresholds. Results Overall, our triangulation method proved to be very precise for static cases, where phantom measurements revealed a fiducial position triangulation precision of mm for stationary targets, but had a spread of mm for targets with 1 mm of motion. Employing triangulation, true motion exceeding 5 mm was detected above 4.0 mm in magnitude 90.9% of the time, a noticeable improvement in comparison to the 28%–65% successful detection rate reported by 2D projection methods. Most detection errors were attributed to depth disparity. Similar to other reports, correlation coefficients for intra‐fractional motion generally indicated no LR/SI, no LR/AP, and weak positive SI/AP correlations for MHRT and SBRT patients. Offsets beyond a 3D tolerance threshold of 5 mm were observed at a rate of 4.25%–5.25%, while the 5 mm 2D out‐of‐tolerance catch‐rate was 1.6%. The AFT was 8.1–8.2 min using the 2D tolerance threshold of 5 mm. In comparison, the estimated AFT for the proposed 3D tolerance monitoring of offsets beyond 2–5 mm was slightly higher at 8.2–10.9 min due to the higher amount of out‐of‐tolerance instances for the higher precision intra‐fractional motion management. Conclusions Our study showcases a promising subsequent kV‐based triangulation method for intra‐fractional prostate motion monitoring. Acquiring 3D motion information results in higher out‐of‐tolerance catch‐rates particularly in the depth dimension of the kV images, which is perpendicular to the treatment beam. Failure to properly observe and catch these offsets would result in sub‐optimal conformality and accuracy of the dose delivery.

Abstract The aging process is universal, and it is characterized by a progressive deterioration and decrease in physiological function leading to decline on the organismal level. Nevertheless, a number of genetic and nongenetic interventions have been described, which successfully extend healthspan and lifespan in different species. Furthermore, a number of clinical trials have been evaluating the feasibility of different interventions to promote human health. The goal of the annual Biological Sciences Section of the Gerontological Society of America meeting was to share current knowledge of different topics in aging research and provide a vision of the future of aging research. The meeting gathered international experts in diverse areas of aging research including basic biology, demography, and clinical and translational studies. Specific topics included metabolism, inflammaging, epigenetic clocks, frailty, senescence, neuroscience, stem cells, reproductive aging, inter-organelle crosstalk, comparative transcriptomics of longevity, circadian clock, metabolomics, and biodemography.

ABSTRACT The clinical relevance of TP53 mutations ( TP53 MUT ) in myeloproliferative neoplasms (MPN) and their prognostic interaction with MPN subtype designation has not been systematically studied. In the current study, 114 patients with MPN harboring TP53 MUT (VAF ≥ 2%) were evaluated for overall survival (OS), calculated from the time of TP53 MUT detection: chronic phase myelofibrosis (MF‐CP; N = 61); blast‐phase (MPN‐BP; N = 31) or accelerated‐phase (MPN‐AP; N = 16) MPN, and polycythemia vera/essential thrombocythemia (PV/ET; N = 6). Sixty‐five (57%) patients harbored International Consensus Classification (ICC)‐defined multihit TP53 MUT and 56 (49%) monosomal/complex karyotype (MK/CK). Majority of MPN‐BP (90%) and MPN‐AP (81%) while 39% of MF‐CP and none of PV/ET patients harbored multihit TP53 MUT . OS in MPN‐BP and MPN‐AP was equally dismal (median 6 vs. 4.5 months, respectively; p = 1.0), regardless of multihit configuration ( p = 0.44), while OS in TP53 MUT MPN‐BP/AP ( N = 47; median 4 months) was inferior to that of a separate cohort ( N = 49) with TP53 wild‐type MPN‐BP/AP (median 11 months; p < 0.01). OS in MF‐CP was significantly shorter with multihit versus non‐multihit TP53 MUT (median 10 vs. 35 months; HR 2.9; p < 0.01), independent of other MF‐relevant genetic risk factors, including ASXL1/SRSF2/U2AF1 mutations. Multihit TP53 MUT was also associated with inferior survival following allogeneic stem cell transplant (ASCT): median 9 months versus “not reached” in patients with ( N = 9) versus without ( N = 8) multihit TP53 MUT ( p < 0.01). The presence of multihit or non‐multihit TP53 MUT in MPN‐BP/AP or multihit TP53 MUT in MF‐CP is associated with exceptionally poor prognosis and justifies inclusion into the ICC category of “myeloid neoplasms with mutated TP53 .” By contrast, detection of non‐multihit TP53 MUT , by itself, might not endanger short‐term survival in MF‐CP, PV, or ET.

AbstractBackgroundConventional approaches for emergent or expedited palliative radiotherapy (RT) involve the application of cumbersome vendor‐provided solutions and/or multiple patient appointments to complete the RT workflow within a compressed timeframe.PurposeThis report delineates the clinical development of an in‐house, semi‐automated Cone‐beam computed tomography (CBCT)‐based simulation‐free platform for expedited palliative RT on conventional linacs, intended to supplant existing techniques employed at this institution.MethodsThe internal software, termed SimFree Wizard (SFW), was engineered utilizing a C#‐based application programming interface integrated within the treatment planning system (TPS). Generated scripts were compiled as stand‐alone executables, with a graphical user interface (GUI) customized via an integrated development environment. The platform was conceived as a framework for accelerated CBCT‐based RT, bypassing the requirement for standard simulation imaging. SFW employs full automation where feasible to minimize user intervention, supplemented by graphical instructions for tasks requiring manual execution. During development, relevant temporal metrics from 10 end‐to‐end tests for palliative spine RT were quantified. User feedback was solicited via a simple questionnaire assessing the overall platform usability. Automated in‐house secondary verification software was developed for validation of the TPS‐calculated monitor units (MUs).ResultsThe mean duration for workflow execution was 41:42 ± 3:18 [mm:ss] (range ∼37–46 min). SFW satisfactorily generated simple, multi‐field CBCT‐based 3D treatment plans within seconds following delineation of the desired treatment area. User feedback indicated enhanced usability compared to previously employed solutions. Validation of the secondary verification software demonstrated accurate results for palliative spine RT and other simple cases wherein the dose calculation point resides in a predominantly homogenous medium.ConclusionA novel in‐house solution for expedited CBCT‐based RT was successfully developed, facilitating completion of the entire workflow within approximately 1‐hour or less for simple palliative/emergent scenarios. Overall, this application is expected to improve the quality and safety of palliative RT while greatly reducing workflow duration, thereby improving access to palliative care.

Background: Therapy-related myeloid neoplasms (t-MN) occur after exposure to DNA-damaging agentsincluding cytotoxic chemo and radiation therapy and are associated with dismal outcomes. TP53 and PPM1D (protein phosphatase Mn2+/Mg2+-dependent 1D) mutations (MT) disrupt the DNA damage response and repair (DDR) pathway, are often co-mutated (25-35%), and although enriched in t-MN (20%), can also be seen in age-related clonal hematopoiesis (CH) (5-10%). While progressive genomic instability with biallelic inactivation has been described as a mechanism of MN progression with TP53 MT, similar mechanisms remain to be elucidated for PPM1D MT. We conducted this study to assess the clinical phenotype and outcomes of TP53 and PPM1D co-mutated (PPM1DTP53m) versus PPM1D MT and wild type TP53 (PPM1DTP53wt) CH and MN. Methods: After approval from the Mayo Clinic IRB, the CH and MN databases were queried and 119 pts with PPM1D MT (CH= 26, clonal cytopenia of unknown significance [CCUS]= 64, myelodysplastic syndrome [MDS]=16, myeloproliferative neoplasm [MPN]=8, acute myeloid leukemia [AML]= 5) were identified. We analyzed baseline characteristics, karyotypic changes, somatic mutations, and survival outcomes between the PPM1DTP53m (n= 22) and PPM1DTP53wt (n= 97) groups (gp). Results: Baseline characteristics The median age in years (yrs) was 68.5 (range [R], 43-89) and 71 (R,30-92) in the PPM1DTP53m and PPM1DTP53wt gp (p= 0.65). 22% vs 18% (p=0.78), 58% vs 32% (p=0.03), 9% vs 32% (p=0.01), 7% vs 4.5% (p=&gt;0.99) and 3% vs 14% (p=0.07) pts met criteria for CH, CCUS, MDS, MPN and AML in the PPM1DTP53wt and PPM1DTP53m gp, respectively. Similarly, 19% vs 18% (p=&gt;0.99), 47% vs 32% (p=0.23), 7% vs 36% (p=0.001), 2% vs 0% (p=&gt;0.99) and 2% vs 14% (p=0.04) of pts met criteria for t-CH, t-CCUS, t-MDS, t-MPN and t-AML in the PPM1DTP53wt and PPM1DTP53m gp, respectively. All PPM1D variants were truncating, located in exon 6 of the gene, with median variant allele frequency (VAF) (8% [R,1-64] vs 5.5% [R,2-24], p= 0.10) and number of PPM1D MT (1 [R,1-5] vs 1 [R, 1-3], p=&gt;0.99) being similar between the two gp. While proportion of pts with t-MN were similar between the two gps (77% vs 86%, p= 0.40), cytogenetic (CG) abnormalities (67% vs 31%, p= 0.006), including monosomies (28% vs 2%, p= &lt;0.001) and complex CG (53% vs 4%, p= &lt;0.001) were more frequent in PPM1DTP53m gp compared to PPM1DTP53wt gp, respectively. The median TP53 MT VAF was 14% (R, 1-14); 50% had multi-hit TP53 MT. Somatic mutations in PPM1D MT myeloid neoplasm with and without TP53 MT Somatic co-mutations excluding TP53 were seen in 71% of pts; 73.5% in PPM1DTP53wt and 59% in PPM1DTP53m gp (p= 0.20). Apart from CBL (0% vs 9%,p=0.03), MT involving DNMT3A (20% vs 14%, p= 0.56), TET2 (15% vs. 14%, p= &gt;0.99), IDH1 (6% vs 19%, p= 0.07), splicing factors (8% vs 14%, p=0.08), ASXL1 (11% vs 4.5%, p=0.46), and JAK2 (9% vs 0%,p=0.20) were comparable between PPM1DTP53wt and PPM1DTP53m gps, respectively. The proportion of pts with ≥ 2 somatic co-mutations excluding TP53 MT were numerically higher in PPM1DTP53wt gp (43%) compared to PPM1DTP53m gp (23%) (p= 0.09). In the PPM1DTP53wt gp, other somatic driver mutations included RUNX1 (n= 7 [7%]; CH [n=1], CCUS [n=2], MDS [n=3], MPN [n=1]), NRAS (n=1 [1%]; CCUS), BCOR (n=1 [1%]; CCUS), and PTPN11 (n=1 [1%]; CCUS). Leukemic transformation and survival With a median follow up of 6.43 months (mo) (R, 0.23-36.4), median overall survival (mOS) in mo among high-risk MDS/AML (16.4 vs 15.0, p= 0.66) and CHIP/CCUS pts (35.79 vs 26.28, p= 0.93) was similar between PPM1DTP53wt and PPM1DTP53m gps, respectively. CHIP/CCUS pts harboring multiple vs single PPM1D MT had an inferior mOS (23.9 vs 36.7 mo, p= 0.05). Three (2.5%) MDS pts underwent AML transformation; 1 (1%) and 2 (10%) in PPM1DTP53wt and PPM1DTP53m gp, respectively (p= 0.07). Conclusions In our study, 18% of pts had concurrent TP53 and PPM1D MT (82% t-MN), with the PPM1DTP53wt genotype being more frequent in CCUS, while the PPM1DTP53m genotype was more frequent in MN. The PPM1DTP53m genotype had a higher frequency of CG abnormalities, including complex CG. Survival outcome was inferior and comparable in PPM1D MT MN with or without TP53 MT. While 73.5 % of the PPM1DTP53wt gp had somatic driver MT and 31% had karyotypic abnormalities, 22% had no additional MT or CG changes (all CH/CCUS), underscoring the need to understand how PPM1D MT contribute to hematopoietic dysfunction and MN.

Background: The majority of patients with chronic myelomonocytic leukemia (CMML) harbor somatic mutations with the most frequent being TET2 (56%), SRSF2 (42%), and ASXL1 (41%) (Leukemia 2020;34:1407). The prognostic relevance of mutations in CMML has been studied extensively with consistent observations suggesting favorable survival impact from TET2 and unfavorable from DNMT3A, ASXL1, NRAS, and SETBP1 mutations. Among the latter, DNMT3A mutation stands out as being the strongest predictor of inferior overall (OS) and leukemia-free (LFS) survival (AJH 2017;92:56). The current study was inspired by observations from exploratory analyses of a large Mayo Clinic CMML cohort. Methods: The current retrospective study was conducted under an institutional review board approved minimum risk protocol that allowed retrospective collection and analysis of data from patient records (Mayo Clinic, Rochester, MN; Scottsdale, AZ; and Jacksonville, FL). Mutations were screened by multi-gene next-generation sequencing (NGS). Diagnostic criteria were according to the International Consensus Classification (Blood 2022; 140:1200). Conventional statistical methods were employed for comparisons of groups and survival analysis (JMP Pro 17.0.0 software SAS Institute, Cary, NC, USA). Results: Exploratory analysis of a well-characterized cohort of 398 Mayo Clinic patients with CMML suggested an unusually long median survival of 96 months for 7 patients harboring PHF6 mutation vs. 27 months for those without the mutation. A subsequent institutional database search focusing on PHF6-mutated CMML identified 20 more cases, expanding the total number of study patients to 418, including 27 PHF6-mutated cases. PHF6 mutations clustered with (co-mutation frequency) TET2 (93%) and less likely to be co-mutated with TP53 (0%), PTPTN11 (0%) and SF3B1 (0%) mutations. PHF6-mutated patients often presented with thrombocytopenia (platelet count &lt;100 x 10(9)/L; 85% vs. 51% in PHF6 wild-type cases; p&lt;0.01). Median time from initial presentation with thrombocytopenia to CMML diagnosis was 4 years (0-14 years). In order to isolate the prognostic impact of PHF6 mutation without associated DNMT3A mutation and independent of TET2 mutation, we considered 4 mutational categories: i) DNMT3A-mutated (N=19; 2 co-expressed PHF6 and 7 TET2 mutations), ii) PHF6-mutated with wild-type DNMT3A (N=25; 23 co-expressed TET2 mutation), iii) TET2-mutated without DNMT3A or PHF6 co-expression (N=169), and iv) wild-type for DNMT3A, TET2, and PHF6 (N=205). Compared to the other groups, PHF6-mutated patients with wild-type DNMT3A were less likely to display leukocytosis (p&lt;0.01) or red cell transfusion need at presentation (p&lt;0.01). At a median follow-up of 18 months, 302 deaths, 78 leukemic transformations, and 44 allogenic stem cell transplants were documented. Univariate analysis for OS disclosed superiority in favor of PHF6-mutated patients with wild-type DNMT3A (median OS 90 months) vs. the other three groups with median OS ranging from 10 months for DNMT3A-mutated patients (HR 20.7, 95% CI 8.1-52.7; p&lt;0.01) to 40 months for TET2-mutated cases with wild-type DNMT3A (HR 3.5, 95% CI 1.7-7.6; p&lt;0.01). Results were similar for multivariable analysis of LFS, where, remarkably, none (0%) of the 25 PHF6-mutated patients with wild-type DNMT3A experienced leukemic transformation vs. 15-30% risk of leukemic transformation in the other three groups (p&lt;0.01). Multivariable analysis for both OS and LFS, which included other risk factors such as red cell transfusion need, leukocytosis (leukocyte count ≥13 x 109/L) and karyotype, confirmed the independent and favorable impact of PHF6 mutation/wild-type DNMT3A (p&lt;0.01 in both instances). Finally, all-inclusive multivariable analysis of karyotype and mutations, in DNMT3A wild-type patients, confirmed the highly favorable and independent prognostic impact of PHF6 mutation (HR 0.18, 95% CI 0.09-0.37). Conclusion: The current study proposes a 2-gene genomic signature (PHF6 mutated and DNMT3A wild-type) to identify a unique subset of CMML patients who often present with thrombocytopenia and whose clinical course is characterized by a very low risk of leukemic transformation, resulting in superior overall survival. Additional studies are warranted to validate these findings and include the information in contemporary risk models.

Background: TP53-mutated (TP53mut) myeloid neoplasms (MN) have particularly poor survival. Discrepancies in defining TP53mut MN have impact on diagnostic work up, therapeutic decisions, and clinical trial enrolment. Methods: MN cases harboring TP53mut (variant allele frequency, VAF ≥2%) were classified based on bone marrow (BM) and peripheral blood (PB) blast % using the 4th revision of the WHO classification (WHO-4R). Biallelic loss was defined as (i) ≥2 mutations each with VAF ≥2%; (ii) single TP53mut VAF ≥2% plus 17p loss detected on karyotype, FISH, or SNP; or (iii) single TP53mut VAF ≥50%, whereas monoallelic loss was defined as single TP53mut VAF &lt;50% without 17p loss. Primary endpoint was overall survival (OS) from diagnosis. Conditional Inference Trees (CIT) analysis was used to determine the true effect of a predictor. Results: We identified 580 cases including WHO-4R defined acute myeloid leukemia (AML, BM/PB blasts ≥20%, n=219, 37.8%), myelodysplastic syndrome (MDS) with excess blasts-2 (EB-2, BM 10-19% or PB 5%-19%, n=92, 15.9%), MDS with excess blasts-1 (EB-1, BM 5-9% or PB 2-4%, n=75, 12.9%), and the remaining were MDS with low blasts (-LB, BM &lt;5% and PB ≤1%, n=194, 33.4%). Of the 578 cases evaluable for allelic status, 437 (75.6%) and 141 (24.4%) had biallelic and monoallelic TP53 loss, respectively. Higher frequency of biallelic loss was observed with increasing blasts. 63.2%, 73.3%, 78% and 86.3% of cases with MDS-LB, MDS-EB1, MDS-EB2, and AML had bi-allelic TP53 loss (P&lt;0.0001). In a multivariable analysis, WHO-4R categories, VAF ≥10%, and biallelic loss were independent predictors of poor survival. CIT analysis further stratified the cohort into 4 groups with distinct survival: (1) MDS-LB (n=194, 33.4%); (2) MDS-EB1/EB2 and AML with VAF &lt;10% (n=38, 6.6%); (3) MDS-EB1/EB2 with VAF ≥10% (n=149, 25.7%); and (4) AML with VAF ≥10% (n=199, 34.3%). MDS-LB cases could be further stratified based on allelic status and complex karyotype (CK): OS was poor for cases with biallelic TP53 loss and monoallelic TP53mut with CK compared to monoallelic TP53mut without CK, who had the most favorable survival (12.7 vs. 13 vs. 34.8 months, P&lt;0.0001).MDS-EB1/EB2/AML with VAF&lt;10% could also be stratified by CK: MDS-EB1/EB2/AML with VAF&lt;10% cases without CK had the favorable survival compared to their counterparts with CK and MDS-EB1/EB2 with VAF ≥10% cases (26.2 vs. 5.6 vs. 6.3 months, P=0.003).MDS-EB1/EB2 with VAF ≥10% was a homogenous subgroup: Biological characteristics including the proportion of cases with biallelic TP53 loss, CK (P&gt;0.9), monosomal karyotype (MK, P&gt;0.9), del 5q (P=0.6), del 7q (P=0.4), and del 17p (P=0.2) were comparable between MDS-EB1 and -EB2 with VAF ≥10%. Importantly, median OS of MDS-EB1 and EB2 with VAF ≥10% was comparable (9.6 vs. 7.2 months, P=0.12). OS in this group was comparable regardless of the blast %, allelic status, or presence of CK.AML with VAF ≥10% had extremely poor survival: Median OS was 3.9 months. Survival of AML was equally poor irrespective of the mono- or biallelic (3.6 vs. 3.9, P=0.89) and non-CK vs. CK (1.8 vs. 3.9 months, P=0.09). Combined, our model would acknowledge poor survival of 514 (89.1%) cases. In contrast, retrospective application of the 5th edition of the WHO classification (WHO-5) and International Consensus Classification (ICC) would classify only 211 (36.4%) and 468 (80.7%) cases as biallelic/multi-hit TP53mut MN, respectively. Finally, the Harrell's c-index of the proposed model, WHO-5, and ICC was 0.66, 0.571, and 0.617, respectively. Conclusion: We provide evidence for subsequent refinements in classifying MN harboring TP53mut. Our salient findings include: (i) acknowledgement of CK as ‘biallelic equivalent’ in MDS-LB; as well as EB1/EB2/AML with VAF 2 to &lt;10% cases; and (ii) demonstration of poor survival of MDS-EB1 and -EB2 with VAF ≥10% regardless of the allelic status. In conclusion, our proposed model acknowledged poor survival of a higher proportion of cases as TP53mut MN compared to the current models with a higher c- index.

Multiple myeloma (MM), an incurable hematologic malignancy, presents significant challenges in managing its recurrence. The bone marrow tumor microenvironment (BM-TME) is critical to MM progression. However, classical immune checkpoint blockade (ICB) therapies, like anti-PD1, have shown limited efficacy in MM. Therefore, elucidating the compositional features and regulatory mechanisms of the TIME in MM is crucial for developing novel immunotherapeutic targets. Tissue-resident memory CD8+ T cells (CD8+ TRMs) reside indefinitely at multiple anatomical locations like the skin, intestine, and bone marrow, functioning as frontline responders against reencountered pathogens and tumors. Whether CD8+ TRMs within BM-TME are involved in the MM cell immune escape and pathogenesis has not been systematically elucidated. Firstly, we leverage single-cell RNA sequencing (scRNA-seq) to dissect the BM CD8+ T cells across patients at different stages of MM, and identify that CD69+CD8+ T subsets were BM residents, and the proportion of BM CD8+ TRMs between health donors and MM patients was comparable. Moreover, our findings indicated that BM CD8+ TRMs exhibit exhausted features and enrichment of NK genes with MM progresses. Importantly, we observed a considerable elevation of the NK gene KLRB1 and its encoded protein CD161 in BM CD8+ TRMs in MM patients. Immunofluorescence labeling of BM biopsy samples from MM patients confirmed the infiltration of the BM-TME with CD3+CD8+CD69+CD161+ cells. Then, we found that KLRB1 expression in BM CD8+ TRMs was higher than that of other classical inhibitory receptor genes (e.g., PDCD1, CTLA4, LAG3), and CLEC2D (encodes the ligand of CD161) was also markedly expressed in patients' MM cells, exhibiting considerably higher than classical immunosuppressive receptor ligand genes (e.g., PDL1, CD80, CD86). Moreover, patients with higher CLEC2D expression had unfavorable overall survival. We further found that the CD161 expression was higher in CD69+CD8+ TRMs compared to circulating CD69-CD8+ T cells, and the percentage of BM CD161+CD8+ TRMs in MM patients was significantly higher than in healthy donors, while there was no difference for PD-1, LAG3, and TIM3 expression. Notably, only a small fraction of BM CD8+ TRMs have mucosal-associated invariant T cell phenotypes (identified by CD161+TCR-Vα7.2+). Additionally, the CLEC2D expression of BM cells was analyzed, and the highest level of CLEC2D expression was observed in B cells from MM patients, followed by plasma cells. Importantly, CLEC2D expression in MM patients was significantly higher than in plasma cells from healthy donors. Next, to investigate whether the blockade of CD161 could restore the function of CD8+ TRMs, the primary BM cells from MM patients were isolated and stimulated ex vivo. Compared to the anti-PD-1 and isotype group, CD161 blocking presents enhanced CD8+ TRMs proliferation and cytokine secretion levels (TNF-α, IFN-γ, IL2), increasing MM cells apoptosis. Chimeric antigen receptor T (CAR-T) cell therapy has achieved impressive success in treating relapsed/refractory MM. How to prevent CAR-T cell exhaustion and enhance function remains a significant challenge. Therefore, we next set out to investigate the impact of CD161 on CAR-T cells. Our results demonstrated that CD161 was highly expressed in early peripheral blood (PB) CAR-T cells after infusion in MM patients but decreased with time. Moreover, the CD161 expression of CAR-T cells was considerably higher than endogenous T cells. Notably, BM-derived CAR-T cells have higher CD161 expression than PB counterparts for the same patient, suggesting CD161 may have BM-specific expression enrichment. After co-culturing MM and CAR-T cells in vitro for 24 hours, we found that blocking CD161 led to an augmented proportion of CAR-T cells, increased IFN-γ secretion, and a notable elevation in tumor cell apoptosis compared to the anti-PD-1 and isotype groups. Taken together, we demonstrated that CD8+ TRMs constitute a vital anti-tumor cell subset in the MM BM-TME, which exhibits a high expression of the inhibitory CD161 receptor, binding to the CLEC2D expressed by MM cells, thereby mediating the dysfunction of CD8+ TRMs. Blocking the CD161-CLEC2D pathway is expected to become a novel target for MM immunotherapy. The combination of CD161 inhibitors and CAR-T cell therapy can potentially enhance the anti-myeloma function of CAR-T cells.