Wang, Junwen

Abstract The role of the bone marrow (BM) microenvironment in regulating the antitumor immune response in Waldenstrom macroglobulinemia (WM) remains limited. Here we transcriptionally and phenotypically profiled non-malignant (CD19− CD138−) BM cells from WM patients with a focus on myeloid derived suppressive cells (MDSCs) to provide a deeper understanding of their role in WM. We found that HLA-DRlowCD11b+CD33+ MDSCs were significantly increased in WM patients as compared to normal controls, with an expansion of predominantly polymorphonuclear (PMN)-MDSCs. Single-cell immunogenomic profiling of WM MDSCs identified an immune-suppressive gene signature with upregulated inflammatory pathways associated with interferon and TNF signaling. Gene signatures associated with an inflammatory and immune suppressive environment were predominately expressed in PMN-MDSCs. In vitro WM PMN-MDSCs demonstrated robust T-cell suppression, and their viability and expansion was notably enhanced by G-CSF and TNFα. Furthermore, BM malignant B-cells attracted PMN-MDSCs to a greater degree than monocytic MDSCs. Collectively, these data suggest that malignant WM B cells actively recruit PMN-MDSCs which promote an immunosuppressive BM microenvironment through a direct T cell inhibition, while release of G-CSF/TNFα in the microenvironment further promotes PMN-MDSC expansion and in turn immune suppression. Targeting PMN-MDSCs may therefore represent a potential therapeutic strategy in patients with WM.

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.

Background: Alzheimer’s disease (AD) is the most common type of age-related dementia, affecting 6.2 million people aged 65 or older according to CDC data. It is commonly agreed that discovering an effective AD diagnosis biomarker could have enormous public health benefits, potentially preventing or delaying up to 40% of dementia cases. Tau neurofibrillary tangles are the primary driver of downstream neurodegeneration and subsequent cognitive impairment in AD, resulting in structural deformations such as hippocampal atrophy that can be observed in magnetic resonance imaging (MRI) scans. Objective: To build a surface-based model to 1) detect differences between APOE subgroups in patterns of tau deposition and hippocampal atrophy, and 2) use the extracted surface-based features to predict cognitive decline. Methods: Using data obtained from different institutions, we develop a surface-based federated Chow test model to study the synergistic effects of APOE, a previously reported significant risk factor of AD, and tau on hippocampal surface morphometry. Results: We illustrate that the APOE-specific morphometry features correlate with AD progression and better predict future AD conversion than other MRI biomarkers. For example, a strong association between atrophy and abnormal tau was identified in hippocampal subregion cornu ammonis 1 (CA1 subfield) and subiculum in e4 homozygote cohort. Conclusion: Our model allows for identifying MRI biomarkers for AD and cognitive decline prediction and may uncover a corner of the neural mechanism of the influence of APOE and tau deposition on hippocampal morphology.

Abstract High grade glioma (HGG) represents a group of devastating diseases with dismal prognosis. Surgical resection of the contrast enhancing (CE) region of HGG remains the mainstay of treatment, but recurrence inevitably arises from the unresected non-contrast enhancing (NE) region, surgically inaccessible due to cancer cell invasion into healthy brain tissue. Due to its critical role in recurrence, understanding of the NE region is central to the improvement of clinical outcomes. We reveal the biological characteristics of this region through image localized multi-regional sampling. We linked microenvironmental characteristics measured by multi-parametric MRI to genomic mutations and transcriptional phenotypes using mixed effect modeling which allowed us to control for individualized patient effects. We first confirmed that T2 is a significant indicator of IDH mutation status in the NE region, being the first description of such a relationship in a HGG cohort. We found the combination of EGFR amplification and CDKN2A homozygous loss was associated with a significantly lower mean diffusivity (MD) compared to double wild type tumors in the NE region, indicating the presence of greater cellular packing and proliferation in EGFR amplification/CDKN2A loss regions. Finally, using single cell pathway based tumor classifications, we showed that nK2, a DSC-MRI metric representing cell size heterogeneity, correlated positively with neuronal signature and negatively with glycolytic/plurimetabolic signature within the NE tumor, indicating that glycolytic/plurimetobolic tumors possessed a high amount of cell size heterogeneity compared to neuronal samples. This hypothesis was supported using digital reference object (DRO) modeling which confirmed that cell size and heterogeneity drove the differential nK2 signal between neuronal and glycolytic/plurimetabolic samples. We identified immune cell infiltrate as one possible mechanism of increased cell size heterogeneity using transcriptomic signature analysis which found more immune cell signatures within glycolytic/plurimetobolic tumors compared to neuronal. Collectively this study demonstrates the central role of multi-parametric MRI as a non-invasive measure of tumor biology and a tool for understanding the clinically critical NE region which can then inform new therapies targeting this region of HGG recurrence. Citation Format: Matthew Flick, Taylor Weiskittel, Kevin Meng-Lin, Fulvio D'Angelo, Francesca Caruso, Shannon Ensign, Mylan Blomquist, Luija Wang, Christopher Sereduk, Gustavo De Leon, Ashley Nespodzany, Javier Urcuyo, Ashlynn Gonzalez, Lee Curtin, Kyle Singleton, Aliya Anil, Natenael Simmineh, Erika Lewis, Teresa Noviello, Reyna Patel, Panwen Wang, Junwen Wang, Jennifer Eschbacher, Andrea Hawkins-Daarud, Pamela Jackson, Kris Smith, Peter Nakaji, Bernard Bendok, Richard Zimmerman, Chandan Krishna, Devi Patra, Naresh Patel, Mark Lyons, Matthew Neal, Kliment Donev, Maciej Mrugala, Alyx Porter, Scott Beeman, Yuxiang Zhou, Leslie Baxter, Christopher Plaisier, Jing Li, Hu Li, Anna Lasorella, Chad Quarles, Kristin Swanson, Michele Ceccarelli, Antonio Iavarone, Nhan Tran, Leland Hu. Multi-parametric MRI maps regional heterogeneity of high grade glioma phenotypes. [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2023; Part 1 (Regular and Invited Abstracts); 2023 Apr 14-19; Orlando, FL. Philadelphia (PA): AACR; Cancer Res 2023;83(7_Suppl):Abstract nr 5621.

Abstract TIGIT (T cell immunoreceptor with Ig and ITIM domains) is a co-inhibitory receptor, highly expressed by Tregs, and contributes to the suppressive function of Tregs by limiting pro-inflammatory Th1 and Th17 cell responses. Herein, we report that subsets of human regulatory B cells (Bregs) express surface TIGIT and can thus play a critical role in immune regulation. Using data from human B cell surface arrays and RNA-seq, we identified two major TIGIT+ human Bregs subsets, exhibiting CD24hiCD27+CD39hiIgD+IgMhiCD1chi marginal zone (MZ)-like and CD24hiCD27+CD39+IgD−IgM+CD1c+ memory Bregs. We further demonstrated that TIGIT+ Bregs were also capable of suppressing immune responses by their ability to express not only IL-10 and PD-L1, but also granzyme B and TGFb1. In addition, TIGIT+ Bregs can also suppress immune responses by altering dendritic cell functions. In line with these observations, liver allograft recipients with donor-specific alloantibody (DSA) had significantly lower numbers of TIGIT+ Bregs than those without DSA. Moreover, the frequency of TIGIT+ Bregs correlated with that of CD4+CD25+CD127lo Tregs, while it inversely correlated with the frequency of follicular helper T cells. We thus concluded that TIGIT+ Bregs play critical roles in immune regulation in liver allograft recipients.

AbstractUnderstanding intratumor heterogeneity is critical for studying tumorigenesis and designing personalized treatments. To decompose the mixed cell population in a tumor, subclones are inferred computationally based on variant allele frequency (VAF) from bulk sequencing data. In this study, we showed that sequencing depth, mean VAF, and variance of VAF of a subclone are confounded. Without considering this effect, current methods require deep-sequencing data (>300× depth) to reliably infer subclones. Here, we present a novel algorithm that incorporates depth-variance and mean-variance dependencies in a clustering error model and successfully identifies subclones in tumors sequenced at depths of as low as 30×. We implemented the algorithm as a model-based adaptive grouping of subclones (MAGOS) method. Analyses of computer simulated data and empirical sequencing data showed that MAGOS outperformed existing methods on minimum sequencing depth, decomposition accuracy, and computation efficiency. The most prominent improvements were observed in analyzing tumors sequenced at depths between 30× and 200×, whereas the performance was comparable between MAGOS and existing methods on deeply sequenced tumors. MAGOS supports analysis of single-nucleotide variants and copy number variants from a single sample or multiple samples of a tumor. We applied MAGOS to whole-exome data of late-stage liver cancers and discovered that high subclone count in a tumor was a significant risk factor of poor prognosis. Lastly, our analysis suggested that sequencing multiple samples of the same tumor at standard depth is more cost-effective and robust for subclone characterization than deep sequencing a single sample. MAGOS is available at github (https://github.com/liliulab/magos).

Abstract Purpose: Regulatory T-cells (Treg) are essential to Tregs homeostasis and modulate the antitumor immune response in patients with lymphoma. However, the biology and prognostic impact of Tregs in splenic marginal zone lymphoma (SMZL) have not been studied. Experimental Design: Biopsy specimens from 24 patients with SMZL and 12 reactive spleens (rSP) from individuals without lymphoma were analyzed by using CITE-seq (cellular indexing of transcriptomes and epitopes by sequencing), CyTOF (mass cytometry) analysis, and flow cytometry to explore the phenotype, transcriptomic profile, and clinical significance of intratumoral Tregs and their subsets. The biological characteristics and cell signaling pathways of intratumoral Treg subsets were confirmed by in vitro functional assays. Results: We found that Tregs are more abundant in SMZL patients' spleens than rSP, and Tregs from patients with SMZL and rSP can be separated into CD161+Treg and CD26+Treg subsets. CD161+Tregs are increased in SMZL but have dysregulated immune function. We found that CD161+Treg and CD26+Tregs have unique gene expression and phenotypic profiles and are differentially correlated with patient outcomes. Specifically, increased CD161+Tregs are significantly associated with a favorable prognosis in patients with SMZL, whereas CD26+Tregs are associated with a poor prognosis. Furthermore, activation of the IL2/STAT5 pathway contributes to the induction of CD26+Tregs and can be reversed by STAT5 inhibition. Conclusions: IL2/STAT5-mediated expansion of CD26+Tregs contributes to a poor clinical outcome in SMZL and may represent a therapeutic opportunity in this disease.