Citations

The translation of AI-generated brain metastases (BM) segmentation into clinical practice relies heavily on diverse, high-quality annotated medical imaging datasets. The BraTS-METS 2023 challenge has gained momentum for testing and benchmarking algorithms using rigorously annotated internationally compiled real-world datasets. This study presents the results of the segmentation challenge and characterizes the challenging cases that impacted the performance of the winning algorithms. Untreated brain metastases on standard anatomic MRI sequences (T1, T2, FLAIR, T1PG) from eight contributed international datasets were annotated in stepwise method: published UNET algorithms, student, neuroradiologist, final approver neuroradiologist. Segmentations were ranked based on lesion-wise Dice and Hausdorff distance (HD95) scores. False positives (FP) and false negatives (FN) were rigorously penalized, receiving a score of 0 for Dice and a fixed penalty of 374 for HD95. The mean scores for the teams were calculated. Eight datasets comprising 1303 studies were annotated, with 402 studies (3076 lesions) released on Synapse as publicly available datasets to challenge competitors. Additionally, 31 studies (139 lesions) were held out for validation, and 59 studies (218 lesions) were used for testing. Segmentation accuracy was measured as rank across subjects, with the winning team achieving a LesionWise mean score of 7.9. The Dice score for the winning team was 0.65 ± 0.25. Common errors among the leading teams included false negatives for small lesions and misregistration of masks in space. The Dice scores and lesion detection rates of all algorithms diminished with decreasing tumor size, particularly for tumors smaller than 100 mm3. In conclusion, algorithms for BM segmentation require further refinement to balance high sensitivity in lesion detection with the minimization of false positives and negatives. The BraTS-METS 2023 challenge successfully curated well- annotated, diverse datasets and identified common errors, facilitating the translation of BM segmentation across varied environments and providing the tools for future development of personalized volumetric reports to patients undergoing BM treatment.

Abstract Monoclonal antibodies targeting CD38 are a therapeutic mainstay in multiple myeloma (MM). Although they have contributed to improved outcomes, most patients still experience disease relapse, and little is known about tumor-intrinsic mechanisms of resistance to these drugs. Antigen escape has been implicated as a mechanism of tumor-cell evasion in immunotherapy. Yet, it is unknown whether MM cells can develop permanent resistance to anti-CD38 antibodies by acquiring genomic events leading to biallelic disruption of the CD38 gene locus. Here, we analyzed whole-genome and whole-exome sequencing data from patients 701 newly diagnosed MM, 67 patients at relapse with naivety to anti-CD38 antibodies, and 50 patients collected at relapse after anti-CD38 antibodies. We report a loss of CD38 in 10 of 50 patients (20%) after CD38 therapy, 3 of whom exhibited a loss of both copies. Two of these cases showed convergent evolution in which distinct subclones independently acquired similar advantageous variants. Functional studies on missense mutations involved in biallelic CD38 events revealed that 2 variants, L153H and C275Y, decreased binding affinity and antibody-dependent cellular cytotoxicity of the commercial antibodies daratumumab and isatuximab. However, a third mutation, R140G, conferred selective resistance to daratumumab, while retaining sensitivity to isatuximab. Clinically, patients with MM are often rechallenged with CD38 antibodies after disease progression and these data suggest that next-generation sequencing may play a role in subsequent treatment selection for a subset of patients.

PURPOSE The combination of the Bruton tyrosine kinase inhibitor ibrutinib with bendamustine-rituximab for first-line treatment of mantle cell lymphoma (MCL) prolonged progression-free survival (PFS), but without improvement in overall survival (OS), likely because of toxicity. Acalabrutinib was shown to be efficacious and less toxic than ibrutinib in a head-to-head trial in chronic lymphocytic leukemia and therefore might lead to better outcomes in MCL. METHODS Patients 65 years and older with previously untreated MCL received acalabrutinib (100 mg twice daily) or placebo (until disease progression or unacceptable toxicity), plus six cycles of bendamustine (90 mg/m 2 once daily; days 1 and 2) and rituximab (375 mg/m 2 as a single dose; day 1) followed by rituximab maintenance in responding patients for 2 years. Crossover to acalabrutinib at disease progression was permitted. The primary end point was PFS per the independent review committee; overall response rate and OS were secondary end points. RESULTS In total, 598 patients were randomly assigned, with 299 in each arm. At a median follow-up of 49.8 months using the reverse Kaplan-Meier method, the median PFS was 66.4 months in the acalabrutinib arm and 49.6 months in the placebo arm (hazard ratio [HR], 0.73 [95% CI, 0.57 to 0.94]; P = .0160). Benefit was seen across all subgroups, including those with high-risk features. Overall response/complete response rates were 91.0%/66.6% and 88.0%/53.5% in the acalabrutinib and placebo arms, respectively. OS was not significantly different (HR, 0.86 [95% CI, 0.65 to 1.13]; P = .27). Grade 3 or greater adverse events were reported in 88.9% and 88.2% in the acalabrutinib and placebo arms, respectively. CONCLUSION The combination of acalabrutinib with bendamustine-rituximab significantly improved PFS. Clinical benefit of acalabrutinib with bendamustine-rituximab was achieved with manageable toxicity.

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 Increasing evidence suggests that Alzheimer’s disease (AD) pathogenesis strongly correlates with neuroinflammation. Peripheral monocytes are crucial components of the human immune system that may play a role in neuroinflammation, but their contribution to AD pathogenesis is largely understudied partially due to the lack of appropriate human models. Here, we present human cortical organoid microphysiological systems (hCO-MPSs) for modeling dynamic AD neuroinflammation mediated by monocytes. By incorporating 3D printed devices into an existing cortical organoid protocol, 96 hCO-MPSs can be established with significantly reduced necrosis and hypoxia as well as enhanced viability within a commonly used 96 well plate, and each hCO-MPS consists of a doughnut-shaped hCO and a 3D printed device per well. Using this approach, monocytes from AD patients exhibit higher infiltration, decreased amyloid-beta (Aβ) clearance, and stronger inflammatory responses compared to monocytes from age-matched control donors. Moreover, pro-inflammatory effects such as elevated astrocyte activation and neuronal apoptosis were observed to be induced by AD monocytes. Furthermore, the significant increase in the expression of IL1B and CCL3, both at the transcriptional and protein levels, indicated the pivotal role of these cytokine and chemokine in monocyte-mediated AD neuroinflammation. Our findings provide insight for understanding monocytes’ role in AD pathogenesis, and the user-friendly MPS models we present are compatible with existing laboratory settings, highlighting their potential for modeling neuroinflammation and developing new therapeutics for various neuroinflammatory diseases.

Introduction Monoclonal antibodies (MoAbs) directed against CD38 are a mainstay in multiple myeloma (MM). Downregulation of CD38 expression is a resistance mechanism of antigenic escape, however, cell surface expression may return. While this tumor cell plasticity is likely driven by epigenetic mechanisms, it is unknown if MM cells can develop permanent resistance to anti-CD38 MoAb by acquiring genomic events leading to biallelic disruption of CD38 thus precluding re-treatment. Here, we report biallelic loss of CD38 as a recurrent mechanism of resistance, often via convergent evolution where distinct subclones are advantaged towards the same selective pressure. Methods We interrogated two sets of patients (Mayo Clinic; n = 31, Calgary University; n =19) with relapse following anti-CD38 MoAbs. Two additional cases of interest were selected from Heidelberg University. Whole genome sequencing (WGS, 60-100x) or whole exome sequencing (WXS), flow cytometry (flow), and bulk RNAseq were used to validate the effect of genomic events on CD38. Missense mutations predicted to affect CD38 MoAb binding were validated with site-directed mutagenesis of a CD38 plasmid, transduction into K562 MM cells, and CD38 MoAb binding assays. Results The prevalence of biallelic inactivation of CD38 at anti-CD38 MoAb relapse was 6% (3/50). An additional selected biallelic event was noted in the Heidelberg set. Each case is described in detail: MM-19 (Calgary) received daratumumab (Dara), lenalidomide, and dexamethasone (Dex), and then had a second exposure to Dara, following relapse. At this second relapse, WGS was performed on two independent extramedullary disease localizations in the liver. Both samples had unique, large deletions (>3 Mb) on chromosome 4 coupled with additional unique focal deletions involving CD38. Flow confirmed the presence of CD138+/CD38- plasma cell populations in both sites. A second case of convergent evolution was seen in a patient (Mayo) after seventh-line Dara/Len/Dex. WGS revealed a common deletion encompassing CD38. On the remaining allele, one subclone had a structural variation-mediated loss of CD38(cancer cell fraction; CCF 40%) while another harbored an L153H missense mutation (CCF 60%). Flow confirmed presence of two CD138+ clones; one with diminished CD38 expression. L153H was functionally validated to abrogate Dara binding, while maintaining CD38 expression. The third sample was collected after fifth-line Dara and sixth-line Dara/Bortezomib/Dex (Heidelberg) and represents another case of convergent evolution with a common loss of chromosome 4 and two subclones detected by WGS (P98Lfs*12, CCF 39% resulting in biallelic inactivation; R140G, CCF 61%). Bulk RNAseq revealed, in comparison to NDMM from CoMMpass, CD38 expression below the first decile. Notably, the R140G variant reduced binding affinity to Dara but not to Isatuximab (Isa). A fourth case (Mayo) relapsing after 6 months of Dara-based therapy similarly had a focal deletion of CD38 with an L18Sfs*16 frameshift deletion on the remaining allele associated with a 4.35-fold reduction in CD38 expression by RNAseq compared to pre-treatment. Interestingly, in 3/4 cases the CD38-negative clone emerged after a second exposure to anti-CD38 MoAb, suggesting emergence with the first exposure and then selection during the second. One other monoallelic event (C275Y; Heidelberg) was detected post-Isa. Copy number of the alternate allele is not known given targeted sequencing methodology, but the variant was seen to reduce binding affinity to both Dara and Isa. The prevalence of monoallelic loss of CD38 in newly diagnosed MM from MMRF CoMMpass was 7% (50/701) and two patients had mutations in CD38, but there were no cases of biallelic loss suggesting these events are driven by treatment pressure. Conclusion Distinction of the mode of CD38 loss is valuable for predicting outcome of re-challenge to anti-CD38 MoAbs. We see here an estimated 5-10% of anti-CD38 MoAb relapses are characterized by genomic antigen escape. Results here are in line with antigen escape to BCMA- and GPRC5D-directed immunotherapies (Lee et al, Nat Med 2023), suggesting that acquired antigen escape resulting from prolonged exposure to the selective pressures of immunotherapy is a recurrent mechanism of resistance in MM. Moreover, in select cases, non-overlapping resistance mutations provide rationale for switching agents at second exposure.

Introduction: Multiple myeloma (MM) with t(11;14)(q13;q32), resulting in CCND1;IGH rearrangement, remains the only subset amenable to targeted therapy with the BCL2 inhibitor, Venetoclax. Nevertheless, not all patients with this translocation respond to treatment and some progress early after initial response. The mechanisms behind this clinical heterogeneity are largely unknown. Methods: To investigate mechanisms of heterogeneity in response to Venetoclax, we interrogated 44 whole genome sequencing (WGS) and 24 whole exome sequencing data from a training cohort including 33 patients with t(11;14) MM treated with Venetoclax at the Mayo Clinic. In a validation cohort, we extracted single nucleotide variants (SNV) and indels from 17 MM patients with t(11;14) with available RNAseq data treated with Venetoclax at Emory University. Results: In the training cohort the median number of lines of therapy was 4 and Venetoclax drug regimens included combinations with dexamethasone, proteasome inhibitors, anti-CD38 agents, and immunomodulatory agents. After a median follow up of 30 months, the median progression free survival (PFS) was 9 months, and 13 patients (39%) achieved a complete remission (CR) to therapy, while 9 (27%) were refractory (i.e. no response or progressed within the first 60 days). For our driver discovery analysis, we integrated SNV, indels, copy number variants (CNV) and structural variants (SV). Deletion 17p (del17p) or TP53 mutation, MCL1 gain/amp, and hyper-APOBEC did not correlate with PFS in our cohort (p=0.59, p=0.71, and p=0.05, respectively). In contrast, mutations involving the MAPK pathway (NRAS, KRAS, BRAF, FGFR3) displayed the strongest association with short PFS, which was retained after adjusting for del17p/TP53 mutations, and 1q gain/amp (HR: 4.76, CI: 1.8-12.9, p=0.002). The association between suboptimal response and resistance to Venetoclax in MM patients with MAPK pathway mutations (n=7) was also observed in our validation cohort from Emory University (p=0.046). This is in line with what has been reported in functional studies in acute myeloid leukemia, where RAS pathway promotes MCL1 under the selective pressure of Venetoclax (Zhang et al. Signal Transduct Target Ther 2022). While MAPK mutations explained most cases of refractory patients (86%), the mechanisms through which initially responding patients subsequently progressed remain unclear. To investigate this, we interrogated 10 patients with pre- and post-treatment samples collected at disease progression (PD). We identified genomic events acquired/selected at PD of genes involved in the anti-apoptotic MCL1/BCL2 and MAPK pathways in 6/10 cases. One patient was refractory to Venetoclax and had double minutes involving BCL2, which caused multiple focal amplifications (>30 copies), likely reflecting extrachromosomal DNA. Another patient progressed with an acquired, focal deletion on BCL2L11. In two patients, the progressing clones harbored mutations in MAPK pathways (one BRAF and the other WNK2). Two patients who relapsed after sustained CR had major branching evolution where a previously undetectable focal deletion of PMAIP1 (NOXA) was seen at PD. Importantly, one of these deletions was a biallelic event. Venetoclax sensitivity assays conducted on PMAIP1 knockout in the MM t(11;14) cell line, KMS12PE, demonstrated significantly reduced cell death across all drug concentrations (p<0.05), validating the loss of NOXA as a mechanism of resistance to Venetoclax in MM. Finally, similarly with what was previously reported in chronic lymphocyte leukemia treated with Venetoclax (e.g., Herling et al Nat Comm 2018), CDKN2A/B monoallelic and biallelic focal deletions were selected post-Venetoclax, at PD in 4 patients while TP53 SNV/deletions were selected from subclonal to clonal in 2 patients at PD. Conclusion: In this study, we demonstrate that mutations involving the MAPK pathway identify a population of MM patients with t(11;14)(CCND1;IGH) who are unlikely to respond to Venetoclax. We also show that genomic events involving the MAPK and BCL2/MCL1 pathways are primarily responsible for disease progression after initial response. These findings highlight the critical role of comprehensive genomic profiling through WGS in identifying mechanisms of resistance to targeted therapies in MM, aiding patient selection for Venetoclax, and guiding the development of novel drug combinations.

Background: Multiple myeloma (MM) is an incurable plasma cell malignancy with several recurrent primary cytogenetic abnormalities including t(11,14)(q13;q32), resulting in increased expression of CCND1. Although 40% of patients with t(11;14) MM respond favorably to the BCL-2 inhibitor venetoclax (ven), 60% fail to achieve an objective response. Previous studies have suggested upregulation of PMAIP1, encoding the pro-apoptotic NOXA protein, in association with ven sensitivity (Leblay, Blood, 2024). Here we sought to evaluate the role for NOXA in ven response in t(11;14) MM. Methods: We analyzed whole exome, whole genome and RNA sequencing data of CD138+ tumor cells from 1,359 patients with MM in the CoMMpass and Mayo Clinic cohorts. The histone mark H3K27ac was assessed in t(11;14)+ KMS12PE and U266B1 human myeloma cell lines using ChIP-seq. Gene expression and protein levels of PMAIP1 (NOXA) were also detected by qRT-PCR and Western blot, respectively. Live cell detection was evaluated by Annexin V and PI using flow cytometry and by YoYo3 dye using IncuCyte following treatment with ven. Results: The t(11;14) MM subtype was associated with a 1.6-5.5-fold increased expression of PMAIP1 compared to non-t(11,14) MM in the CoMMpass and Mayo Clinic cohorts (p<0.0001). These findings are consistent with a significant correlation between CCND1 and PMAIP1 expression in both cohorts (p<0.001). To evaluate whether cyclin D1 is directly associated with increasing expression of PMAIP1 in t(11;14) MM, we generated a cyclin D1 knock out (KO) version of the U266B1 line using CRISPR/cas9. Two independent U266B1 KO cell lines showed reduced expression of PMAIP1 and NOXA protein compared to the WT U266B1 consistent with a role for cyclin D1 in regulating PMAIP1 expression. In addition, the KO lines had a 4.8-fold reduced H3K27ac signal at the PMAIP1 locus in comparison to the WT control suggesting a role for cyclin D1 in promoting H3K27ac. We previously observed increased expression of CCND1 in association with the G allele compared to the A allele of rs9344 in both the CoMMpass and Mayo Clinic cohorts (Leukemia, in press). We next asked whether the G risk allele of rs9344 within CCND1 impacted gene expression and H3k27ac levels at the PMAIP1 locus. To study the allele specific effect of CCND1, we used the CRISPR/cas9 edited KMS12PE WT (rs9344 AA) and Knock In (KI) (rs9344 GG) cell lines. In the KI cells, the G allele was accompanied by an ~2-fold increased expression of PMAIP1 and NOXA protein detected by western blot (<0.05). Increased PMAIP1 expression in KI lines was also accompanied by a 2.2-fold increased H3K27ac further suggesting a role for the G allele of rs9344 in promoting H3K27ac and gene expression at the PMAIP1 locus. Finally, we evaluated the consequence of PMAIP1 expression levels in relation to ven response. Increased expression of PMAIP1 and elevated levels of NOXA protein in the KMS12PE KI lines was associated with increased sensitivity to ven resulting in ~50% increase in cell death of KI vs. WT lines (p<0.02). Conclusion: We demonstrate a role for cyclin D1 in regulating PMAIP1 expression through H3K27ac resulting in increased NOXA protein. Levels of NOXA are associated with the G allele of rs9344 and increased venetoclax response.

Introduction: Cyclin D1, encoded by the CCND1 gene, is primarily known as a major regulator of cell cycle progression. Emerging studies have demonstrated a role for cyclin D1 in regulating gene transcription. In mantle cell lymphoma with t(11;14), cyclin D1 strongly binds to the promoters of highly transcribed genes and its overexpression led to genome-wide transcriptional downregulation. However, the transcriptional role for cyclin D1 in t(11;14) MM remains undefined. This study aims to elucidate the global chromatin and transcriptional changes induced by cyclin D1 in t(11;14) MM. Methods: Knockout (KO) of CCND1 in U266B1, a MM cell line carrying t(11;14), was performed by the Synthego (Redwood City, CA) using CRISPR/Cas9, which generated a frameshift insertion of base pair (bp) “A” after 69,641,478 bp within exon 1. Two independent KO clones were used for downstream analysis. ChIP-seq and ATAC-seq were performed for KO (n=2) and wild-type clones (WT, n=2). ChIP-seq used antibodies against H3K4me3 (promoter) and H3K4me1 and H3K27ac (enhancer). Libraries were sequenced to 2x51 bp on the Illumina NextSeq 2000 platform. Total RNA libraries (in triplicate) were sequenced to 2x150 bp. ChIP-seq and ATAC-seq reads were mapped to hg38 with BWA, and peaks were identified by MACS2 at FDR 0.01. RNA-seq reads were mapped to hg38 with STAR and gene raw counts were estimated using featureCounts. Differential peaks and differentially expressed genes between KO and WT were identified using edgeR at FDR 0.05 and fold-change 2. Only peaks present in at least 2 samples and genes having reads-per-million value 0.1 in at least 1 sample were included, excluding those from chrY. Pathway analysis was performed with the Enrichr package. Results: Differential expression analysis identified 2-fold more genes that were down-regulated than up-regulated in KO compared to WT (467 vs. 233). Pathway analysis revealed an enrichment of TNF-alpha signaling via NF-kappaB, inflammatory response, and apoptosis in association with down-regulated genes, while up-regulated genes are enriched in myogenesis, KRAS signaling, and epithelial mesenchymal transition pathway. In parallel, differential chromatin analysis identified 2.5% (1,161) of the peaks from H3K27ac and 0.5-1% (359-610) of the peaks for the other 3 marks that showed significant changes of signal levels between KO and WT. For the promoter mark H3K4me3, there are largely a comparable number of peaks showing increased (191 peaks) and decreased (168 peaks) signals, respectively, in KO compared to WT. However, for ATAC-seq and H3K4me1, about 70% of the differential peaks showed an increased signal in KO, in contrast with H3K27ac for which 60% of differential peaks showed decreased signal in KO. The preferential loss of H3K27ac occupancy is consistent with the observation of a higher proportion of genes being down-regulated in KO. These results indicated that knockout of CCND1 led to a global increase of chromatin accessibility, but a reduction of H3K27ac and gene expression. Further annotation revealed that, for H3K27ac and H3K4me1/ATAC-seq, 92% and 97% of the differential peaks were located >2kb away from the transcription start sites (TSSs) of expressed genes. Even for H3K4me3, 88% of the differential peaks are >2kb away from TSSs. The results suggest that the changes of chromatin upon CCND1 knockout predominantly occurred at distal regulatory regions. For up- and down-regulated genes, we estimated ChIP-seq and ATAC-seq signals in the TSS 2kb regions, revealing that up-regulated genes were overall associated with increased signals. Further analysis found that, for ATAC-seq and H3K4me3, 44% and 59% of the differential peaks located in the promoter regions were associated with differentially expressed genes, which dropped to ~20% for the two enhancer marks. These data indicate that changes in chromatin state impact the expression of the associated genes. Conclusion: Knockout of CCND1 impacts local chromatin state, particularly at enhancer regions, and the transcriptional program. In promoter-proximal regulatory regions, the changes of chromatin accessibility and H3K4me3 are strongly correlated with the expression changes of nearby genes. Further studies will identify the genes targeted by the differential enhancers and their possible roles in MM pathogenesis mediated by cyclin D1 following the t(11;14) event.

Acute pancreatitis (AP) is among the most common hospital gastrointestinal diagnoses; understanding the mechanisms underlying the severity of AP is critical for development of new treatment options for this disease. Here, we evaluate the biological function of phosphofructo-2-kinase/fructose-2,6-biphosphatase 3 (PFKFB3) in AP pathogenesis in 2 independent genetically engineered mouse models of AP. PFKFB3 was elevated in AP and severe AP (SAP), and KO of Pfkfb3 abrogated the severity of alcoholic SAP (FAEE-SAP). Using a combination of genetic, pharmacological, and molecular studies, we defined the interaction of PFKFB3 with inositol 1,4,5-trisphosphate receptor (IP3R) as a key event mediating this phenomenon. Further analysis demonstrated that the interaction between PFKFB3 and IP3R promotes FAEE-SAP severity by altering intracellular calcium homeostasis in acinar cells. Together, our results support a PFKFB3-driven mechanism controlling AP pathobiology and define this enzyme as a therapeutic target to ameliorate the severity of this condition.