Citations

Abstract Background: Although CAR T-cell therapy has achieved remarkable response rates in adult patients (pts) with relapsed and refractory B-cell non-Hodgkin’s lymphoma (NHL), 20-30% of pts still have primary refractory disease, and 30-40% of pts relapse, often within 6 months of therapy. Given that durable clinical response can be achieved without CAR T cell persistence, we hypothesized that durable remission is correlated with the expansion of tumor-reactive T-cell receptor (TCR) clones following CAR T-cell infusion. To assess this, we examined the TCR repertoire of adult pts with B-NHL who received CD19 CAR T-cell (axi-cel) therapy to identify correlates of clinical response. Methods: Single-cell (sc) RNA sequencing (seq) and scTCR-seq were performed on peripheral blood mononuclear cells from B-NHL pts before lymphodepletion chemotherapy (pre), at peak CAR-T expansion (peak), and one month following CAR-T (M1). scRNAseq and scTCR-seq data were processed using CellRanger v7.0.1, which were then analyzed using Immunarch (v0.9.0) for TCR clonality and sample diversity. Diversity was assessed using the Gini coefficient to indicate clone size inequality (0 = maximum diversity). TCR diversity and phenotype of expanded and emerged TCR clones post CAR-T were compared across clinical response categories: CR (complete remission for at least 6 months), PD1 (primary refractory), and PD2 (relapsed after partial or complete response). Results: We analyzed 107,035 T cells from 77 serial samples of 32 pts (16 CR, 4 PD1, 12 PD2) and 5 healthy controls. Compared to pts with PD1 and PD2, pts with CR had a larger Gini coefficient for TCR clones at M1 relative to Pre and Peak. The Gini coefficient was the largest for tumor recirculating CD8 T effector memory (Tem) and activated CD8 peripheral memory (Tpm) cells at M1 for pts with CR. Examining the repertoire of emerged and expanded TCR clones at Peak and M1 from Pre, pts with PD2, compared to pts with CR, had a higher percentage of activated CD8 T peripheral memory (Tpm) and T regulatory cells (Treg) among the emerged clones and activated CD8 Tpm among the expanded clones. Patients with PD1, relative to patients with CR, had more proliferating and non-activated CD8 Tpm among the expanded and emerged clones. Conclusion: NHL pts with CR after CAR-T had more robust clonal expansion of T cells, particularly tumor recirculating CD8 Tem and activated CD8 Tpm, compared to pts with relapsed or refractory disease. In contrast, pts with relapsed or refractory disease had an expansion of TCR clones with a larger proportion of CD8 Tpm and Tregs. Endogenous tumor-reactive TCR expansion may help predict favorable clinical response following CAR T-cell therapy. Citation Format: Aiman J. Faruqi, Panwen Wang, Radhika Bansal, Andre De Menezes Silva Corraes, Henan Zhang, Zuoyi Shao, Kevin Regan, Eider Moreno Cortes, Arushi Khurana, Nora Bennani, Yucai Wang, Paul Hampel, Jonas Paludo, Patrick Johnston, Hemant Murthy, Madiha Iqbal, Stephen Ansell, Javier Munoz, Allison Rosenthal, Mohamed Kharfan Dabaja, Januario E. Castro, Saad Kenderian, Haidong Dong, Ying Li, Yi Lin. T-cell receptor repertoire changes associated with clinical response in patients with B-cell non-Hodgkin’s lymphoma receiving CD19 CAR-T therapy [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2024; Part 1 (Regular Abstracts); 2024 Apr 5-10; San Diego, CA. Philadelphia (PA): AACR; Cancer Res 2024;84(6_Suppl):Abstract nr 53.

Abstract Background: Although targeting myeloid cells has become one of the next generation targets for cancer immunotherapy, there is a lack of therapeutic tools to specifically target these cells. Our group has developed a novel humanized PD-L1-targeting antibody, H1A, that induces PD-L1 degradation via blockade of PD-L1 and CMTM6, therefore limiting PD-L1 intrinsic signaling while also preventing interactions with PD-1. In the immune cell compartment, CMTM6 expression is restricted to myeloid cells, making these cells the primary target of H1A. Recently, we unexpectedly found that targeting intrinsic PD-L1 in myeloid cells via H1A resulted in increased release of CCL2 from myeloid cells. In this study, we examined the role of CCL2/CCR2 signaling in CD8 T cell function using human peripheral blood-derived immune cells and evaluated the therapeutic potential of H1A antibody using our humanized PD-L1/PD-1 mouse model. Methods: Human PBMCs derived from healthy donors were treated with H1A or Atezolizumab, a blocking PD-L1 antibody, followed by mass cytometry and flow cytometry to evaluate immune cell phenotype and functional states. PD-L1 expressing myeloid cells isolated from human PBMCs were treated with H1A or Atezolizumab and assessed for CCL2 secretion using multiplex immunoassays. CCL2 production and transcription were confirmed by flow cytometry and single-cell RNA-sequencing. CCR2 blockade was used on human PBMCs to evaluate the effects of CCL2/CCR2 signaling on T cell function alongside H1A treatment. Finally, in vivo antitumor activity of H1A and Atezolizumab were compared in two humanized PD-L1 tumor models. Results: H1A reduced the frequency of regulatory T cells and increased effector T cells, in whole PBMCs upon global T cell activation with anti-CD3/CD28. H1A enhanced the secretion and production of CCL2 from myeloid cells, which was confirmed at a transcript level. CCR2 blockade diminished CD8 T cells’ ability to differentiate into effector cells in H1A treated PBMCs. H1A demonstrated comparable therapeutic effects with Atezolizumab in treatment of immunogenic and ICI-responsive mouse tumor MC38, which is CD8 T cell-dependent, but showed superior therapeutic effects in treatment of ICI-insensitive mouse tumor E0771 compared to Atezolizumab. Both treatments elicited durable memory specific immune responses upon E0771 tumor rechallenge. Conclusion: H1A demonstrated a superior antitumor activity compared to Atezolizumab in our preclinical models via a new mechanism by which H1A caused degradation of intrinsic PD-L1 in myeloid cells. This causes more release of CCL2 to improve T cell function via CCR2 upon T cell activation. Thus, H1A is emerging as a new immunotherapy tool to target PD-L1/CMTM6 in myeloid cells for treatment of cancers that are refractory to current immune checkpoint inhibitor therapy. Citation Format: Michelle Alexandria Hsu, Xin Liu, Ying Li, Whitney Barham, Kevin Pavelko, Jake Hirdler, Fabrice Lucien, Haidong Dong. Targeting PD-L1/CMTM6 in myeloid cells results in elevated release of CCL2 to improve CD8 T cell-mediated antitumor response [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2024; Part 1 (Regular Abstracts); 2024 Apr 5-10; San Diego, CA. Philadelphia (PA): AACR; Cancer Res 2024;84(6_Suppl):Abstract nr 2664.

Abstract Background: Assessing treatment response in gliomas with high fidelity remains challenging. Cell-free DNA (cfDNA) is increasingly recognized as a liquid biomarker in plasma and cerebrospinal fluid (CSF), showing promise in optimizing therapeutic regimens for various cancers. While CSF has been explored as a source of cfDNA for gliomas, understanding the dynamics of cfDNA abundance and sequencing changes during treatment and recurrence is essential. Leveraging our biobank of longitudinal CSF specimens obtained through CSF access devices, we aimed to determine the potential utility of CSF cfDNA as a monitoring tool for treatment response in glioma patients. Materials and methods: Longitudinal CSF specimens were collected through Ommaya reservoirs (NCT04692337) or ventriculoperitoneal shunts (NCT04692324) from patients with brain tumors, including gliomas. cfDNA was extracted from 1-5 mL of CSF based on available sample volumes and quantified using Qubit. Next-Generation sequencing using PredicineCARE or PredicineSCORE low-pass whole genome sequencing (LPWGS) was performed based on the available cfDNA. Results: CSF collected before versus after resection showed an increase in quantified cfDNA (2.97x, range: 1.58-5.26x), indicating the impact of increased parenchymal disruption and closer contact with CSF. Subsequently, CSF cfDNA decreased during chemoradiation and adjuvant treatment and increased at recurrence, even in patients co-enrolled in immunotherapy trials for IL-7 agonists and pembrolizumab. Notably, copy number burden (CNB) increased at recurrence in three patients with known disease progression and decreased throughout immunotherapy treatment in one patient despite concerns of pseudoprogression. In a patient with a hypermutated glioblastoma, CSF cfDNA revealed over 59 genomic alterations, including MAP2K1, KIT, and PDGFRA. Another patient with a progressing EGFR-amplified GBM showed over 200 new variants in the final sample, with EGFR copy number increasing from 2 to over 30. Conclusion: Analysis of cfDNA throughout a patient’s disease course, encompassing resection and treatment with standard-of-care and experimental therapies, may aid in disease monitoring and treatment response. Further evaluation in a larger patient cohort is needed to determine the sensitivity of changes in mutations, CNB, and cfDNA quantity for glioma disease burden. Citation Format: Riviere-Cazaux Cecile, Xiaoxi Dong, Chao Dai, Wei Mo, Pan Du, Shidong Jia, Terence Burns. Monitoring glioma treatment response through longitudinal analysis of cell-free DNA in cerebrospinal fluid: Insights from a comprehensive study using next-generation sequencing [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2024; Part 1 (Regular Abstracts); 2024 Apr 5-10; San Diego, CA. Philadelphia (PA): AACR; Cancer Res 2024;84(6_Suppl):Abstract nr 2294.

Abstract Introduction: Single-cell RNA-sequencing (scRNA-seq) pairing with single-cell T-Cell-Receptor-sequencing (scTCR-seq) enables profiling of gene expression and TCR repertoires in individual T cells, elevating our understanding of T-cell-mediated immunity. Here, we present immunopipe, a comprehensive and flexible pipeline (https://github.com/pwwang/immunopipe) for analyzing the paired data. Besides the command-line tool, it offers a user-friendly web interface that allows users with varying programming expertise to configure, initiate, and monitor the pipeline. By combining extensive functionality with ease of use, immunopipe empowers users to unlock valuable insights from the data, which notably advances our understanding of immune response and fosters the development of innovative immunotherapies. Results: Immunopipe is a comprehensive analytical pipeline for scRNA-seq and scTCR-seq data, offering various modules leveraging existing tools and novel methods to ensure data quality and enable downstream analytics. It performs quality control for scRNA-seq data, clusters cells based on gene expression to facilitate the identification of cell populations and provides insights into T-cell heterogeneity and functional diversity through T-cell selection and clustering. Enrichment analysis identifies markers and pathways for T-cell characterization. Our automated T-cell annotation tool annotates cell types accurately and hierarchically. For scTCR-seq data, immunopipe analyzes TCR clones, including clonality, diversity, gene usage, and repertoire overlap. Clone residency analysis examines the presence of clones across samples. Integrative analyses combine scRNA-seq and scTCR-seq data, providing a comprehensive understanding of the cellular landscape. Immunopipe incorporates diverse cell group information to analyze T-cell populations, which enables granular perspectives of the data for comprehensive insights. Immunopipe offers a flexible platform, allowing users to customize the analyses as needed. It supports both core and optional analyses through configuration, where users have control over the adjustable parameters. The pipeline relies on cell grouping information that can be referenced from metadata or modified to generate new variables for analysis and data filtering. Immunopipe can be executed locally or on scheduler platforms, and a docker image is provided for compatibility and simplified installation. It also offers a web interface that facilitates configuration, execution monitoring, and access to results. Conclusion: Immunopipe is a flexible and user-friendly pipeline for analyzing scRNA-seq and scTCR-seq data. It empowers researchers to gain valuable insights that could enhance the understanding of immune responses and facilitate the development of innovative immunotherapies. Citation Format: Panwen Wang, Haidong Dong, Yue Yu, Shuwen Zhang, Zhifu Sun, Jean-Pierre A. Kocher, Junwen Wang, Lin Yi, Ying Li. Immunopipe: A comprehensive and flexible scRNA-seq and scTCR-seq data analysis pipeline [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2024; Part 1 (Regular Abstracts); 2024 Apr 5-10; San Diego, CA. Philadelphia (PA): AACR; Cancer Res 2024;84(6_Suppl):Abstract nr 4959.

Abstract Metastasis-ablative radiotherapy is a new therapeutic approach that prolongs survival in patients with advanced cancers. Tumour irradiation can elicit an immunogenic cell death resulting in the expansion of tumour-reactive T cells and eradication of both irradiated and distant non-irradiated metastases, referred to as the “abscopal effect”. This abscopal response, nevertheless, is observed in a minority of patients suggesting uncharacterized mechanisms that prevent robust systemic anti-tumor immune response. Tumour-derived extracellular vesicles (TdEVs) have recently emerged as mediators of antitumor immunity but their immunomodulatory role in response to radiotherapy remains largely unknown. Recent clinical studies by our lab have demonstrated that circulating levels of tumour-derived extracellular vesicles (tdEVs) are elevated in response to radiotherapy and high tdEV concentrations are associated with reduced levels of peripheral cytotoxic CD8+ T cells and higher risk of disease recurrence in advanced prostate cancer. Based on our findings, we hypothesize that radiotherapy induces the release of tumor-derived EVs that suppress anti-tumor immunity. Prostate cancer cells (PC3 and DU145) were irradiated with clinically used doses of ionizing radiation (20Gy x 1 or 10Gy x 3 fractions) and tdEVs were co-cultured with human CD8+ T cells to measure their proliferation and cytotoxic activity. Increased levels of tdEVs were released by tumour cells treated with radiotherapy. A higher frequency of impaired T cell cytotoxic activity and a reduction in CD8 T cell proliferation was found in the 10Gy x 3 radiation regimen. Proteomic analysis and western blot showed an increased expression of the immune checkpoint molecule B7-H3 on the tumour cell surface and tdEVs in response to radiotherapy. While coculture of CD8 T cells with B7-H3 depleted tdEVs led to recovery of T-cell cytotoxic ability, overexpression of the B7-H3 in tdEVs inhibited cytotoxic activity. This study provides a rationale for future studies in which we will test the immune-mediated effects of radiotherapy-induced tdEVs in immunocompetent mouse models and investigate B7-H3 blockade as a therapeutic approach to boost abscopal response in the treatment of advanced prostate cancer. Citation Format: Chisom J. Uzendu, Yohan Kim, Roxane R. Lavoie, Haidong Dong, Sean Park, Fabrice Lucien-Matteoni. Impact of radiotherapy on extracellular vesicle mediated anti-tumor immunity in oligometastatic prostate cancer [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2024; Part 1 (Regular Abstracts); 2024 Apr 5-10; San Diego, CA. Philadelphia (PA): AACR; Cancer Res 2024;84(6_Suppl):Abstract nr 4995.

Abstract Background Cholangiocarcinoma (CCA) is a heterogeneous malignancy. Platinum-based chemotherapy (Gemcitabine and Cisplatin) has been the standard treatment for over a decade and is the base for the current immunotherapy combinations. Response rates are typically between 20-30%. Herein, we present outcomes of preclinical efficacy trials on CCA PDX models in correlation with baseline multiomics profiling of CCA patient-derived xenograft (PDX) models to identify molecular features associated with response. Methods Tumor samples from twenty-eight CCA PDX were used to isolate RNA and protein to be included in multiomics analysis (Phosphoproteomics, Proteomics, and Transcriptomics). A subset of PDX models were treated with various doses of Gemcitabine and Cisplatin. The Median Efficacy Index (MEI) was calculated as the ratio of the difference in median of delta in tumor volume at the last day of the study between the treatment and the control arms divided by the median of delta in tumor volume of the control at the last day of the treatment (positive= growth, negative=inhibition, and 0=no efficacy). A Hierarchical All-Against-All algorithm was utilized to compare multiomics features to the MEI and pathway analysis completed for identified molecular features. Results After conducting preclinical efficacy trials, the five tested CCA PDX models were classified into sensitive, moderately sensitive, and resistant models based on the MEI on each treatment arm. The top up-regulated genes were identified for Gemcitabine response (SRPK2, GLK2, and RIOK2) and Cisplatin response (NTRK1, PRPF4B, and SRPK2). The top up-regulated genes related to Gemcitabine resistance were NTRK1, PIP5K1C, and CSNK2A2, while those related to Cisplatin resistance were SRPK1, PRPF4B, and CSNK2A1. Conclusion Cholangiocarcinoma multiomics signatures can predict therapeutic response to standard-of-care chemotherapy in preclinical models. Novel mechanisms of Gemcitabine and Cisplatin resistance were identified related to NTRK1, SRPK1 and SRPK2. Resistant signatures will be challenged with the appropriate targeted therapies on the resistant CCA PDX models. Citation Format: Amro M. Abdelrahman, Danielle M. Carlson, Erik Jessen, Dong-Gi Mun, Isaac Lynch, Alessandro Fogliati, Stella Konadu Adjei Antwi, Aushinie Abeynayake, Akhilesh Pandey, Mark J. Truty, Gregory J. Gores, Rory L. Smoot. Pretreatment multiomics in conjunction with preclinical trials on patient-derived xenograft models of cholangiocarcinoma can identify chemotherapy response signatures [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2024; Part 1 (Regular Abstracts); 2024 Apr 5-10; San Diego, CA. Philadelphia (PA): AACR; Cancer Res 2024;84(6_Suppl):Abstract nr 6912.

Abstract Background: Hepatocellular carcinoma (HCC) is one of the leading causes of cancer-related mortality worldwide. In sub-Saharan Africa, HCC is the second most common cancer in men and the sixth in women. African patients present with HCC at a young age and at advanced stage of disease. Improved surveillance can increase early HCC detection. The GALAD score has good sensitivity in the diagnosis of HCC in Asia, Europe and North America, however, it has not been validated in an African cohort. The aim of this study was to assess the performance of the GALAD score in the diagnosis of HCC in Ghanaian patients. Methods: Clinical data from patients with cirrhosis (n=93) or HCC (n=78) from out-patient hepatology clinics at 3 teaching hospitals in Ghana were abstracted, and serum samples analyzed. A logistic regression model predicting HCC status based on GALAD score was constructed to obtain the Receiver Operating Characteristics (ROC) curve for GALAD. The area under the curve (AUC) with 95% confidence interval was calculated and used to assess the performance of the GALAD score in surveillance for HCC. Results: The median AFP (5010.2 vs. 3.4), L3% (28.3 vs. 0.3) and DCP (507.6 vs. 0.2) levels were all significantly higher among HCC cases than the controls (p<0.0001). The median GALAD score was also higher among HCC patients (8.0 vs. -4.1, p<0.0001). The AUC of the GALAD score in the diagnosis of HCC was 0.86 (95% CI 0.79 - 0.92). At a cut-off value of -0.37, the GALAD score had a sensitivity of 0.81 and a specificity of 0.86. In a re-fit of the logistic regression model, the most significant factors in the GALAD model in our cohort were AFP-L3% and log(DCP). The AUC of the new GALAD model after model selection was 0.89. Conclusion: The GALAD score can be used to detect HCC in Ghana and has the potential to improve HCC surveillance in low-resource settings in sub-Saharan Africa. Citation Format: Yvonne Ayerki Nartey, Ju Dong Yang, Tyler Zemla, Joshua Ayawin, Shadrack Osei Asibey, Mohamed El-Kassas, Sally Afua Bampoh, Amoako Duah, Adwoa Agyei-Nkansah, Yaw Asante Awuku, Mary Yeboah Afihene, Hiroyuki Yamada, Jun Yin, Amelie Plymoth, Lewis Rowland Roberts. Accuracy of the GALAD serologic model in the diagnosis of hepatocellular carcinoma in Ghanaian patients [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2024; Part 1 (Regular Abstracts); 2024 Apr 5-10; San Diego, CA. Philadelphia (PA): AACR; Cancer Res 2024;84(6_Suppl):Abstract nr 4802.

Abstract Background The electronic health record (EHR) presents new opportunities for the timely identification of patients at high risk of critical illness and the implementation of preventive strategies. This study aims to externally validate an EHR‐based Elders Risk Assessment (ERA) score to identify older patients at high risk of future critical illness during a primary care visit. Methods This historical cohort study included patients aged ≥65 years who had primary care visits at Mayo Clinic Rochester, MN, between July 2019 and December 2021. The ERA score at the time of the primary care visit was used to predict critical illness, defined as death or ICU admission within 1 year of the visit. Results A total of 12,885 patients were included in the analysis. The median age at the time of the primary care visit was 75 years, with 44.6% being male. 93.7% of participants were White, and 64.2% were married. The median (25th, 75th percentile) ERA score was 4 (0, 9). 11.3% of study participants were admitted to the ICU or died within 1 year of the visit. The ERA score predicted critical illness within 1 year of a primary care visit with an area under the receiver operating characteristic curve of 0.84 (95% CI 0.83–0.85), which indicates good discrimination. An ERA score of 9 was identified as optimal for implementing and testing potential preventive strategies, with the odds ratio of having the primary outcome in patients with ERA score ≥9 being 11.33 (95%CI 9.98–12.87). Conclusions This simple EHR‐based risk assessment model can predict critical illness within 1 year of primary care visits in older patients. The findings of this study can serve as a basis for testing and implementation of preventive strategies to promote the well‐being of older adults at risk of critical illness and its consequences.

AbstractBackgroundClot based assays used for lupus anticoagulant (LAC) detection are typically interpreted in a qualitative fashion and may not reflect LAC potency. In this cross‐sectional study, we describe a method for quantifying the LAC titer using serial (dependent) two‐fold dilutions in normal pooled plasma.MethodsSerial dilutions of 51 residual plasma samples from 50 patients were tested using the Russell's viper venom screening time (DRVVT) and activated partial thromboplastin screening time (APTT) methodologies. The measured clotting times and the corresponding dilution factors were then used to derive a four‐parameter logistic model. The LAC titer for each patient was interpolated as the sample dilution that corresponds to the upper reference interval limit of the corresponding assay.ResultsCalculated APTT and DRVVT LAC titers displayed a strong linear correlation (R2 = 0.84) between each other, but not with the degree of prolongation of the APTT/DRVVT screening time in the neat undiluted samples. Using data driven partitioning, patients could be grouped into low (<10) or high (≥10) DRVVT LAC titer. There were no significant differences in anticardiolipin (aCL) or anti‐beta 2 glycoprotein 1 (aB2GPI) antibody levels or prevalence of thromboembolic events between low and high LAC titer groups. In contrast, antiphosphatidylserine/prothrombin (aPS/PT) IgM antibody levels, but not IgG, were significantly higher in the high LAC titer group.ConclusionsThe degree of prolongation of the APTT/DRVVT screening time is not correlated with the LAC titer. Only aPS/PT IgM antibodies levels were strongly correlated with the LAC titers. Additional studies are warranted to determine clinical implications of high LAC titers.

Background: NAFLD is highly prevalent with limited treatment options. Bile acids (BAs) increase in the systemic circulation and liver during NAFLD progression. Changes in plasma membrane localization and zonal distribution of BA transporters can influence transport function and BA homeostasis. However, a thorough characterization of how NAFLD influences these factors is currently lacking. This study aimed to evaluate the impact of NAFLD and the accompanying histologic features on the functional capacity of key hepatocyte BA transporters across zonal regions in human liver biopsies. Methods: A novel machine learning image classification approach was used to quantify relative zonal abundance and plasma membrane localization of BA transporters (bile salt export pump [BSEP], sodium-taurocholate cotransporting polypeptide, organic anion transporting polypeptide [OATP] 1B1 and OATP1B3) in non-diseased (n = 10), NAFL (n = 9), and NASH (n = 11) liver biopsies. Based on these data, membrane-localized zonal abundance (MZA) measures were developed to estimate transporter functional capacity. Results: NAFLD diagnosis and histologic scoring were associated with changes in transporter membrane localization and zonation. Increased periportal BSEP MZA (mean proportional difference compared to non-diseased liver of 0.090) and decreased pericentral BSEP MZA (−0.065) were observed with NASH and also in biopsies with higher histologic scores. Compared to Non-diseased Liver, periportal OATP1B3 MZA was increased in NAFL (0.041) and NASH (0.047). Grade 2 steatosis (mean proportional difference of 0.043 when compared to grade 0) and grade 1 lobular inflammation (0.043) were associated with increased periportal OATP1B3 MZA . Conclusions: These findings provide novel mechanistic insight into specific transporter alterations that impact BA homeostasis in NAFLD. Changes in BSEP MZA likely contribute to altered BA disposition and pericentral microcholestasis previously reported in some patients with NAFLD. BSEP MZA assessment could inform future development and optimization of NASH-related pharmacotherapies.