Early Injury Landscape in Vein Harvest by Single-Cell and Spatial Transcriptomics.

BACKGROUND: Vein graft failure following cardiovascular bypass surgery results in significant patient morbidity and cost to the healthcare system. Vein graft injury can occur during autogenous vein harvest and preparation, as well as after implantation into the arterial system, leading to the development of intimal hyperplasia, vein graft stenosis, and, ultimately, bypass graft failure. Although previous studies have identified maladaptive pathways that occur shortly after implantation, the specific signaling pathways that occur during vein graft preparation are not well defined and may result in a cumulative impact on vein graft failure. We, therefore, aimed to elucidate the response of the vein conduit wall during harvest and following implantation, probing the key maladaptive pathways driving graft failure with the overarching goal of identifying therapeutic targets for biologic intervention to minimize these natural responses to surgical vein graft injury. METHODS: Employing a novel approach to investigating vascular pathologies, we harnessed both single-nuclei RNA-sequencing and spatial transcriptomics analyses to profile the genomic effects of vein grafts after harvest and distension, then compared these findings to vein grafts obtained 24 hours after carotid-carotid vein bypass implantation in a canine model (n=4). RESULTS: Spatial transcriptomic analysis of canine cephalic vein after initial conduit harvest and distention revealed significant enrichment of pathways (P<0.05) involved in the activation of endothelial cells (ECs), fibroblasts, and vascular smooth muscle cells, namely pathways responsible for cellular proliferation and migration and platelet activation across the intimal and medial layers, cytokine signaling within the adventitial layer, and ECM (extracellular matrix) remodeling throughout the vein wall. Subsequent single-nuclei RNA-sequencing analysis supported these findings and further unveiled distinct EC and fibroblast subpopulations with significant upregulation (P<0.05) of markers related to endothelial injury response and cellular activation of ECs, fibroblasts, and vascular smooth muscle cells. Similarly, in vein grafts obtained 24 hours after arterial bypass, there was an increase in myeloid cell, protomyofibroblast, injury response EC, and mesenchymal-transitioning EC subpopulations with a concomitant decrease in homeostatic ECs and fibroblasts. Among these markers were genes previously implicated in vein graft injury, including VCAN, FBN1, and VEGFC, in addition to novel genes of interest, such as GLIS3 and EPHA3. These genes were further noted to be driving the expression of genes implicated in vascular remodeling and graft failure, such as IL-6, TGFBR1, SMAD4, and ADAMTS9. By integrating the spatial transcriptomics and single-nuclei RNA-sequencing data sets, we highlighted the spatial architecture of the vein graft following distension, wherein activated and mesenchymal-transitioning ECs, myeloid cells, and fibroblasts were notably enriched in the intima and media of distended veins. Finally, intercellular communication network analysis unveiled the critical roles of activated ECs, mesenchymal-transitioning ECs, protomyofibroblasts, and vascular smooth muscle cells in upregulating signaling pathways associated with cellular proliferation (MDK [midkine], PDGF [platelet-derived growth factor], VEGF [vascular endothelial growth factor]), transdifferentiation (Notch), migration (ephrin, semaphorin), ECM remodeling (collagen, laminin, fibronectin), and inflammation (thrombospondin), following distension. CONCLUSIONS: Vein conduit harvest and distension elicit a prompt genomic response facilitated by distinct cellular subpopulations heterogeneously distributed throughout the vein wall. This response was found to be further exacerbated following vein graft implantation, resulting in a cascade of maladaptive gene regulatory networks. Together, these results suggest that distension initiates the upregulation of pathological pathways that may ultimately contribute to bypass graft failure and presents potential early targets warranting investigation for targeted therapies. This work highlights the first applications of single-nuclei and spatial transcriptomic analyses to investigate venous pathologies, underscoring the utility of these methodologies and providing a foundation for future investigations.

Single-cell analysis reveals an antiviral network that controls Zika virus infection in human dendritic cells.

Zika virus (ZIKV) is a mosquito-borne flavivirus that caused an epidemic in the Americas in 2016 and is linked to severe neonatal birth defects, including microcephaly and spontaneous abortion. To better understand the host response to ZIKV infection, we adapted the 10× Genomics Chromium single-cell RNA sequencing (scRNA-seq) assay to simultaneously capture viral RNA and host mRNA. Using this assay, we profiled the antiviral landscape in a population of human monocyte-derived dendritic cells infected with ZIKV at the single-cell level. The bystander cells, which lacked detectable viral RNA, expressed an antiviral state that was enriched for genes coinciding predominantly with a type I interferon (IFN) response. Within the infected cells, viral RNA negatively correlated with type I IFN-dependent and -independent genes (the antiviral module). We modeled the ZIKV-specific antiviral state at the protein level, leveraging experimentally derived protein interaction data. We identified a highly interconnected network between the antiviral module and other host proteins. In this work, we propose a new paradigm for evaluating the antiviral response to a specific virus, combining an unbiased list of genes that highly correlate with viral RNA on a per-cell basis with experimental protein interaction data. IMPORTANCE: Zika virus (ZIKV) remains a public health threat given its potential for re-emergence and the detrimental fetal outcomes associated with infection during pregnancy. Understanding the dynamics between ZIKV and its host is critical to understanding ZIKV pathogenesis. Through ZIKV-inclusive single-cell RNA sequencing (scRNA-seq), we demonstrate on the single-cell level the dynamic interplay between ZIKV and the host: the transcriptional program that restricts viral infection and ZIKV-mediated inhibition of that response. Our ZIKV-inclusive scRNA-seq assay will serve as a useful tool for gaining greater insight into the host response to ZIKV and can be applied more broadly to the flavivirus field.

Single cell analysis reveals an antiviral network that controls Zika virus infection in human dendritic cells.

Zika virus (ZIKV) is a mosquito-borne flavivirus that caused an epidemic in the Americas in 2016 and is linked to severe neonatal birth defects, including microcephaly and spontaneous abortion. To better understand the host response to ZIKV infection, we adapted the 10x Genomics Chromium single cell RNA sequencing (scRNA-seq) assay to simultaneously capture viral RNA and host mRNA. Using this assay, we profiled the antiviral landscape in a population of human moDCs infected with ZIKV at the single cell level. The bystander cells, which lacked detectable viral RNA, expressed an antiviral state that was enriched for genes coinciding predominantly with a type I interferon (IFN) response. Within the infected cells, viral RNA negatively correlated with type I IFN dependent and independent genes (antiviral module). We modeled the ZIKV specific antiviral state at the protein level leveraging experimentally derived protein-interaction data. We identified a highly interconnected network between the antiviral module and other host proteins. In this work, we propose a new paradigm for evaluating the antiviral response to a specific virus, combining an unbiased list of genes that highly correlate with viral RNA on a per cell basis with experimental protein interaction data. Our ZIKV-inclusive scRNA-seq assay will serve as a useful tool to gaining greater insight into the host response to ZIKV and can be applied more broadly to the flavivirus field.

Integrated single-nuclei and spatial transcriptomic analysis reveals propagation of early acute vein harvest and distension injury signaling pathways following arterial implantation.

Background: Vein graft failure (VGF) following cardiovascular bypass surgery results in significant patient morbidity and cost to the healthcare system. Vein graft injury can occur during autogenous vein harvest and preparation, as well as after implantation into the arterial system, leading to the development of intimal hyperplasia, vein graft stenosis, and, ultimately, bypass graft failure. While previous studies have identified maladaptive pathways that occur shortly after implantation, the specific signaling pathways that occur during vein graft preparation are not well defined and may result in a cumulative impact on VGF. We, therefore, aimed to elucidate the response of the vein conduit wall during harvest and following implantation, probing the key maladaptive pathways driving graft failure with the overarching goal of identifying therapeutic targets for biologic intervention to minimize these natural responses to surgical vein graft injury. Methods: Employing a novel approach to investigating vascular pathologies, we harnessed both single-nuclei RNA-sequencing (snRNA-seq) and spatial transcriptomics (ST) analyses to profile the genomic effects of vein grafts after harvest and distension, then compared these findings to vein grafts obtained 24 hours after carotid-cartoid vein bypass implantation in a canine model (n=4). Results: Spatial transcriptomic analysis of canine cephalic vein after initial conduit harvest and distention revealed significant enrichment of pathways (P < 0.05) involved in the activation of endothelial cells (ECs), fibroblasts (FBs), and vascular smooth muscle cells (VSMCs), namely pathways responsible for cellular proliferation and migration and platelet activation across the intimal and medial layers, cytokine signaling within the adventitial layer, and extracellular matrix (ECM) remodeling throughout the vein wall. Subsequent snRNA-seq analysis supported these findings and further unveiled distinct EC and FB subpopulations with significant upregulation (P < 0.00001) of markers related to endothelial injury response and cellular activation of ECs, FBs, and VSMCs. Similarly, in vein grafts obtained 24 hours after arterial bypass, there was an increase in myeloid cell, protomyofibroblast, injury-response EC, and mesenchymal-transitioning EC subpopulations with a concomitant decrease in homeostatic ECs and fibroblasts. Among these markers were genes previously implicated in vein graft injury, including VCAN (versican), FBN1 (fibrillin-1), and VEGFC (vascular endothelial growth factor C), in addition to novel genes of interest such as GLIS3 (GLIS family zinc finger 3) and EPHA3 (ephrin-A3). These genes were further noted to be driving the expression of genes implicated in vascular remodeling and graft failure, such as IL-6, TGFBR1, SMAD4, and ADAMTS9. By integrating the ST and snRNA-seq datasets, we highlighted the spatial architecture of the vein graft following distension, wherein activated and mesenchymal-transitioning ECs, myeloid cells, and FBs were notably enriched in the intima and media of distended veins. Lastly, intercellular communication network analysis unveiled the critical roles of activated ECs, mesenchymal transitioning ECs, protomyofibroblasts, and VSMCs in upregulating signaling pathways associated with cellular proliferation (MDK, PDGF, VEGF), transdifferentiation (Notch), migration (ephrin, semaphorin), ECM remodeling (collagen, laminin, fibronectin), and inflammation (thrombospondin), following distension. Conclusions: Vein conduit harvest and distension elicit a prompt genomic response facilitated by distinct cellular subpopulations heterogeneously distributed throughout the vein wall. This response was found to be further exacerbated following vein graft implantation, resulting in a cascade of maladaptive gene regulatory networks. Together, these results suggest that distension initiates the upregulation of pathological pathways that may ultimately contribute to bypass graft failure and presents potential early targets warranting investigation for targeted therapies. This work highlights the first applications of single-nuclei and spatial transcriptomic analyses to investigate venous pathologies, underscoring the utility of these methodologies and providing a foundation for future investigations.

The Simple prEservatioN of Single cElls method for cryopreservation enables the generation of single-cell immune profiles from whole blood.

Introduction: Current multistep methods utilized for preparing and cryopreserving single-cell suspensions from blood samples for single-cell RNA sequencing (scRNA-seq) are time-consuming, requiring trained personnel and special equipment, so limiting their clinical adoption. We developed a method, Simple prEservatioN of Single cElls (SENSE), for single-step cryopreservation of whole blood (WB) along with granulocyte depletion during single-cell assay, to generate high quality single-cell profiles (SCP). Methods: WB was cryopreserved using the SENSE method and peripheral blood mononuclear cells (PBMCs) were isolated and cryopreserved using the traditional density-gradient method (PBMC method) from the same blood sample (n=6). The SCPs obtained from both methods were processed using a similar pipeline and quality control parameters. Further, entropy calculation, differential gene expression, and cellular communication analysis were performed to compare cell types and subtypes from both methods. Results: Highly viable (86.3 ± 1.51%) single-cell suspensions (22,353 cells) were obtained from the six WB samples cryopreserved using the SENSE method. In-depth characterization of the scRNA-seq datasets from the samples processed with the SENSE method yielded high-quality profiles of lymphoid and myeloid cell types which were in concordance with the profiles obtained with classical multistep PBMC method processed samples. Additionally, the SENSE method cryopreserved samples exhibited significantly higher T-cell enrichment, enabling deeper characterization of T-cell subtypes. Overall, the SENSE and PBMC methods processed samples exhibited transcriptional, and cellular communication network level similarities across cell types with no batch effect except in myeloid lineage cells. Discussion: Comparative analysis of scRNA-seq datasets obtained with the two cryopreservation methods i.e., SENSE and PBMC methods, yielded similar cellular and molecular profiles, confirming the suitability of the former method's incorporation in clinics/labs for cryopreserving and obtaining high-quality single-cells for conducting critical translational research.

Single-cell analysis reveals altered tumor microenvironments of relapse- and remission-associated pediatric acute myeloid leukemia.

Acute myeloid leukemia (AML) microenvironment exhibits cellular and molecular differences among various subtypes. Here, we utilize single-cell RNA sequencing (scRNA-seq) to analyze pediatric AML bone marrow (BM) samples from diagnosis (Dx), end of induction (EOI), and relapse timepoints. Analysis of Dx, EOI scRNA-seq, and TARGET AML RNA-seq datasets reveals an AML blasts-associated 7-gene signature (CLEC11A, PRAME, AZU1, NREP, ARMH1, C1QBP, TRH), which we validate on independent datasets. The analysis reveals distinct clusters of Dx relapse- and continuous complete remission (CCR)-associated AML-blasts with differential expression of genes associated with survival. At Dx, relapse-associated samples have more exhausted T cells while CCR-associated samples have more inflammatory M1 macrophages. Post-therapy EOI residual blasts overexpress fatty acid oxidation, tumor growth, and stemness genes. Also, a post-therapy T-cell cluster associated with relapse samples exhibits downregulation of MHC Class I and T-cell regulatory genes. Altogether, this study deeply characterizes pediatric AML relapse- and CCR-associated samples to provide insights into the BM microenvironment landscape.

Single-cell RNA sequencing distinctly characterizes the wide heterogeneity in pediatric mixed phenotype acute leukemia.

BACKGROUND: Mixed phenotype acute leukemia (MPAL), a rare subgroup of leukemia characterized by blast cells with myeloid and lymphoid lineage features, is difficult to diagnose and treat. A better characterization of MPAL is essential to understand the subtype heterogeneity and how it compares with acute myeloid leukemia (AML) and acute lymphoblastic leukemia (ALL). Therefore, we performed single-cell RNA sequencing (scRNAseq) on pediatric MPAL bone marrow (BM) samples to develop a granular map of the MPAL blasts and microenvironment landscape. METHODS: We analyzed over 40,000 cells from nine pediatric MPAL BM samples to generate a single-cell transcriptomic landscape of B/myeloid (B/My) and T/myeloid (T/My) MPAL. Cells were clustered using unsupervised single-cell methods, and malignant blast and immune clusters were annotated. Differential expression analysis was performed to identify B/My and T/My MPAL blast-specific signatures by comparing transcriptome profiles of MPAL with normal BM, AML, and ALL. Gene set enrichment analysis (GSEA) was performed, and significantly enriched pathways were compared in MPAL subtypes. RESULTS: B/My and T/My MPAL blasts displayed distinct blast signatures. Transcriptomic analysis revealed that B/My MPAL profile overlaps with B-ALL and AML samples. Similarly, T/My MPAL exhibited overlap with T-ALL and AML samples. Genes overexpressed in both MPAL subtypes' blast cells compared to AML, ALL, and healthy BM included MAP2K2 and CD81. Subtype-specific genes included HBEGF for B/My and PTEN for T/My. These marker sets segregated bulk RNA-seq AML, ALL, and MPAL samples based on expression profiles. Analysis comparing T/My MPAL to ETP, near-ETP, and non-ETP T-ALL, showed that T/My MPAL had greater overlap with ETP-ALL cases. Comparisons among MPAL subtypes between adult and pediatric samples showed analogous transcriptomic landscapes of corresponding subtypes. Transcriptomic differences were observed in the MPAL samples based on response to induction chemotherapy, including selective upregulation of the IL-16 pathway in relapsed samples. CONCLUSIONS: We have for the first time described the single-cell transcriptomic landscape of pediatric MPAL and demonstrated that B/My and T/My MPAL have distinct scRNAseq profiles from each other, AML, and ALL. Differences in transcriptomic profiles were seen based on response to therapy, but larger studies will be needed to validate these findings.

Taxane chemotherapy induces stromal injury that leads to breast cancer dormancy escape.

A major cause of cancer recurrence following chemotherapy is cancer dormancy escape. Taxane-based chemotherapy is standard of care in breast cancer treatment aimed at killing proliferating cancer cells. Here, we demonstrate that docetaxel injures stromal cells, which release protumor cytokines, IL-6 and granulocyte colony stimulating factor (G-CSF), that in turn invoke dormant cancer outgrowth both in vitro and in vivo. Single-cell transcriptomics shows a reprogramming of awakened cancer cells including several survival cues such as stemness, chemoresistance in a tumor stromal organoid (TSO) model, as well as an altered tumor microenvironment (TME) with augmented protumor immune signaling in a syngeneic mouse breast cancer model. IL-6 plays a role in cancer cell proliferation, whereas G-CSF mediates tumor immunosuppression. Pathways and differential expression analyses confirmed MEK as the key regulatory molecule in cancer cell outgrowth and survival. Antibody targeting of protumor cytokines (IL-6, G-CSF) or inhibition of cytokine signaling via MEK/ERK pathway using selumetinib prior to docetaxel treatment prevented cancer dormancy outgrowth suggesting a novel therapeutic strategy to prevent cancer recurrence.

Pediatric T-cell acute lymphoblastic leukemia blast signature and MRD associated immune environment changes defined by single cell transcriptomics analysis.

Different driver mutations and/or chromosomal aberrations and dysregulated signaling interactions between leukemia cells and the immune microenvironment have been implicated in the development of T-cell acute lymphoblastic leukemia (T-ALL). To better understand changes in the bone marrow microenvironment and signaling pathways in pediatric T-ALL, bone marrows collected at diagnosis (Dx) and end of induction therapy (EOI) from 11 patients at a single center were profiled by single cell transcriptomics (10 Dx, 5 paired EOI, 1 relapse). T-ALL blasts were identified by comparison with healthy bone marrow cells. T-ALL blast-associated gene signature included SOX4, STMN1, JUN, HES4, CDK6, ARMH1 among the most significantly overexpressed genes, some of which are associated with poor prognosis in children with T-ALL. Transcriptome profiles of the blast cells exhibited significant inter-patient heterogeneity. Post induction therapy expression profiles of the immune cells revealed significant changes. Residual blast cells in MRD+ EOI samples exhibited significant upregulation (P < 0.01) of PD-1 and RhoGDI signaling pathways. Differences in cellular communication were noted in the presence of residual disease in T cell and hematopoietic stem cell compartments in the bone marrow. Together, these studies generate new insights and expand our understanding of the bone marrow landscape in pediatric T-ALL.

The 'omics of obesity in B-cell acute lymphoblastic leukemia.

The obesity pandemic currently affects more than 70 million Americans and more than 650 million individuals worldwide. In addition to increasing susceptibility to pathogenic infections (eg, SARS-CoV-2), obesity promotes the development of many cancer subtypes and increases mortality rates in most cases. We and others have demonstrated that, in the context of B-cell acute lymphoblastic leukemia (B-ALL), adipocytes promote multidrug chemoresistance. Furthermore, others have demonstrated that B-ALL cells exposed to the adipocyte secretome alter their metabolic states to circumvent chemotherapy-mediated cytotoxicity. To better understand how adipocytes impact the function of human B-ALL cells, we used a multi-omic RNA-sequencing (single-cell and bulk transcriptomic) and mass spectroscopy (metabolomic and proteomic) approaches to define adipocyte-induced changes in normal and malignant B cells. These analyses revealed that the adipocyte secretome directly modulates programs in human B-ALL cells associated with metabolism, protection from oxidative stress, increased survival, B-cell development, and drivers of chemoresistance. Single-cell RNA sequencing analysis of mice on low- and high-fat diets revealed that obesity suppresses an immunologically active B-cell subpopulation and that the loss of this transcriptomic signature in patients with B-ALL is associated with poor survival outcomes. Analyses of sera and plasma samples from healthy donors and those with B-ALL revealed that obesity is associated with higher circulating levels of immunoglobulin-associated proteins, which support observations in obese mice of altered immunological homeostasis. In all, our multi-omics approach increases our understanding of pathways that may promote chemoresistance in human B-ALL and highlight a novel B-cell-specific signature in patients associated with survival outcomes.

Cross center single-cell RNA sequencing study of the immune microenvironment in rapid progressing multiple myeloma.

Despite advancements in understanding the pathophysiology of Multiple Myeloma (MM), the cause of rapid progressing disease in a subset of patients is still unclear. MM's progression is facilitated by complex interactions with the surrounding bone marrow (BM) cells, forming a microenvironment that supports tumor growth and drug resistance. Understanding the immune microenvironment is key to identifying factors that promote rapid progression of MM. To accomplish this, we performed a multi-center single-cell RNA sequencing (scRNA-seq) study on 102,207 cells from 48 CD138- BM samples collected at the time of disease diagnosis from 18 patients with either rapid progressing (progression-free survival (PFS) < 18 months) or non-progressing (PFS > 4 years) disease. Comparative analysis of data from three centers demonstrated similar transcriptome profiles and cell type distributions, indicating subtle technical variation in scRNA-seq, opening avenues for an expanded multicenter trial. Rapid progressors depicted significantly higher enrichment of GZMK+ and TIGIT+ exhausted CD8+ T-cells (P = 0.022) along with decreased expression of cytolytic markers (PRF1, GZMB, GNLY). We also observed a significantly higher enrichment of M2 tolerogenic macrophages in rapid progressors and activation of pro-proliferative signaling pathways, such as BAFF, CCL, and IL16. On the other hand, non-progressive patients depicted higher enrichment for immature B Cells (i.e., Pre/Pro B cells), with elevated expression for markers of B cell development (IGLL1, SOX4, DNTT). This multi-center study identifies the enrichment of various pro-tumorigenic cell populations and pathways in those with rapid progressing disease and further validates the robustness of scRNA-seq data generated at different study centers.

Survival Genie, a web platform for survival analysis across pediatric and adult cancers.

The genomics data-driven identification of gene signatures and pathways has been routinely explored for predicting cancer survival and making decisions related to targeted treatments. A large number of packages and tools have been developed to correlate gene expression/mutations to the clinical outcome but lack the ability to perform such analysis based on pathways, gene sets, and gene ratios. Furthermore, in this single-cell omics era, the cluster markers from cancer single-cell transcriptomics studies remain an underutilized prognostic option. Additionally, no bioinformatics online tool evaluates the associations between the enrichment of canonical cell types and survival across cancers. Here we have developed Survival Genie, a web tool to perform survival analysis on single-cell RNA-seq (scRNA-seq) data and a variety of other molecular inputs such as gene sets, genes ratio, tumor-infiltrating immune cells proportion, gene expression profile scores, and tumor mutation burden. For a comprehensive analysis, Survival Genie contains 53 datasets of 27 distinct malignancies from 11 different cancer programs related to adult and pediatric cancers. Users can upload scRNA-seq data or gene sets and select a gene expression partitioning method (i.e., mean, median, quartile, cutp) to determine the effect of expression levels on survival outcomes. The tool provides comprehensive results including box plots of low and high-risk groups, Kaplan-Meier plots with univariate Cox proportional hazards model, and correlation of immune cell enrichment and molecular profile. The analytical options and comprehensive collection of cancer datasets make Survival Genie a unique resource to correlate gene sets, pathways, cellular enrichment, and single-cell signatures to clinical outcomes to assist in developing next-generation prognostic and therapeutic biomarkers. Survival Genie is open-source and available online at https://bbisr.shinyapps.winship.emory.edu/SurvivalGenie/ .

Single cell transcriptomic landscape of diabetic foot ulcers.

Diabetic foot ulceration (DFU) is a devastating complication of diabetes whose pathogenesis remains incompletely understood. Here, we profile 174,962 single cells from the foot, forearm, and peripheral blood mononuclear cells using single-cell RNA sequencing. Our analysis shows enrichment of a unique population of fibroblasts overexpressing MMP1, MMP3, MMP11, HIF1A, CHI3L1, and TNFAIP6 and increased M1 macrophage polarization in the DFU patients with healing wounds. Further, analysis of spatially separated samples from the same patient and spatial transcriptomics reveal preferential localization of these healing associated fibroblasts toward the wound bed as compared to the wound edge or unwounded skin. Spatial transcriptomics also validates our findings of higher abundance of M1 macrophages in healers and M2 macrophages in non-healers. Our analysis provides deep insights into the wound healing microenvironment, identifying cell types that could be critical in promoting DFU healing, and may inform novel therapeutic approaches for DFU treatment.

Comprehensive Characterization of the Multiple Myeloma Immune Microenvironment Using Integrated scRNA-seq, CyTOF, and CITE-seq Analysis.

As part of the Multiple Myeloma Research Foundation (MMRF) immune atlas pilot project, we compared immune cells of multiple myeloma bone marrow samples from 18 patients assessed by single-cell RNA sequencing (scRNA-seq), mass cytometry (CyTOF), and cellular indexing of transcriptomes and epitopes by sequencing (CITE-seq) to understand the concordance of measurements among single-cell techniques. Cell type abundances are relatively consistent across the three approaches, while variations are observed in T cells, macrophages, and monocytes. Concordance and correlation analysis of cell type marker gene expression across different modalities highlighted the importance of choosing cell type marker genes best suited to particular modalities. By integrating data from these three assays, we found International Staging System stage 3 patients exhibited decreased CD4+ T/CD8+ T cells ratio. Moreover, we observed upregulation of RAC2 and PSMB9, in natural killer cells of fast progressors compared with those of nonprogressors, as revealed by both scRNA-seq and CITE-seq RNA measurement. This detailed examination of the immune microenvironment in multiple myeloma using multiple single-cell technologies revealed markers associated with multiple myeloma rapid progression which will be further characterized by the full-scale immune atlas project. Significance: scRNA-seq, CyTOF, and CITE-seq are increasingly used for evaluating cellular heterogeneity. Understanding their concordances is of great interest. To date, this study is the most comprehensive examination of the measurement of the immune microenvironment in multiple myeloma using the three techniques. Moreover, we identified markers predicted to be significantly associated with multiple myeloma rapid progression.

Leukemia vaccine overcomes limitations of checkpoint blockade by evoking clonal T cell responses in a murine acute myeloid leukemia model.

We have developed a personalized vaccine whereby patient derived leukemia cells are fused to autologous dendritic cells, evoking a polyclonal T cell response against shared and neo-antigens. We postulated that the dendritic cell (DC)/AML fusion vaccine would demonstrate synergy with checkpoint blockade by expanding tumor antigen specific lymphocytes that would provide a critical substrate for checkpoint blockade mediated activation. Using an immunocompetent murine leukemia model, we examined the immunologic response and therapeutic efficacy of vaccination in conjunction with checkpoint blockade with respect to leukemia engraftment, disease burden, survival and the induction of tumor specific immunity. Mice treated with checkpoint blockade alone had rapid leukemia progression and demonstrated only a modest extension of survival. Vaccination with DC/AML fusions resulted in the expansion of tumor specific lymphocytes and disease eradication in a subset of animals, while the combination of vaccination and checkpoint blockade induced a fully protective tumor specific immune response in all treated animals. Vaccination followed by checkpoint blockade resulted in upregulation of genes regulating activation and proliferation in memory and effector T cells. Long term survivors exhibited increased T cell clonal diversity and were resistant to subsequent tumor challenge. The combined DC/AML fusion vaccine and checkpoint blockade treatment offers unique synergy inducing the durable activation of leukemia specific immunity, protection from lethal tumor challenge and the selective expansion of tumor reactive clones.

Integrated Skin Transcriptomics and Serum Multiplex Assays Reveal Novel Mechanisms of Wound Healing in Diabetic Foot Ulcers.

Nonhealing diabetic foot ulcers (DFUs) are characterized by low-grade chronic inflammation, both locally and systemically. We prospectively followed a group of patients who either healed or developed nonhealing chronic DFUs. Serum and forearm skin analysis, both at the protein expression and the transcriptomic level, indicated that increased expression of factors such as interferon-γ (IFN-γ), vascular endothelial growth factor, and soluble vascular cell adhesion molecule-1 were associated with DFU healing. Furthermore, foot skin single-cell RNA sequencing analysis showed multiple fibroblast cell clusters and increased inflammation in the dorsal skin of patients with diabetes mellitus (DM) and DFU specimens compared with control subjects. In addition, in myeloid cell DM and DFU upstream regulator analysis, we observed inhibition of interleukin-13 and IFN-γ and dysregulation of biological processes that included cell movement of monocytes, migration of dendritic cells, and chemotaxis of antigen-presenting cells pointing to an impaired migratory profile of immune cells in DM skin. The SLCO2A1 and CYP1A1 genes, which were upregulated at the forearm of nonhealers, were mainly expressed by the vascular endothelial cell cluster almost exclusively in DFU, indicating a potential important role in wound healing. These results from integrated protein and transcriptome analyses identified individual genes and pathways that can potentially be targeted for enhancing DFU healing.

Modulation of lymphocyte-mediated tissue repair by rational design of heterocyclic aryl hydrocarbon receptor agonists.

Aryl hydrocarbon receptor (AHR) is an essential regulator of gut immunity and a promising therapeutic target for inflammatory bowel disease (IBD). Current AHR agonists are inadequate for clinical translation due to low activity, inadequate pharmacokinetics, or toxicity. We synthesized a structurally diverse library and used integrated computational and experimental studies to discover mechanisms governing ligand-receptor interaction and to design potent drug leads PY109 and PY108, which display physiochemical drug-likeness properties, desirable pharmacokinetic profiles, and low toxicity. In a murine model of dextran sulfate sodium-induced colitis, orally administered compounds increase interleukin-22 (IL-22) production and accelerate mucosal healing by modulating mucosal adaptive and innate lymphoid cells. AHR and IL-22 pathway induction was confirmed using RNA sequencing and characterization of the lymphocyte protein-protein interaction network. Significant induction of IL-22 was also observed using human T cells from patients with IBD. Our findings support rationally designed AHR agonists for IBD therapy.

Topical Application of a Mast Cell Stabilizer Improves Impaired Diabetic Wound Healing.

Impaired wound healing in the diabetic foot is a major problem often leading to amputation. Mast cells have been shown to regulate wound healing in diabetes. We developed an indole-carboxamide type mast cell stabilizer, MCS-01, which proved to be an effective mast cell degranulation inhibitor in vitro and can be delivered topically for prolonged periods through controlled release by specifically designed alginate bandages. In diabetic mice, both pre- and post-wounding, topical MCS-01 application accelerated wound healing comparable to that achieved with systemic mast cell stabilization. Moreover, MCS-01 altered the macrophage phenotype, promoting classically activated polarization. Bulk transcriptome analysis from wounds treated with MCS-01 or placebo showed that MCS-01 significantly modulated the mRNA and microRNA profile of diabetic wounds, stimulated upregulation of pathways linked to acute inflammation and immune cell migration, and activated the NF-κB complex along with other master regulators of inflammation. Single-cell RNA sequencing analysis of 6,154 cells from wounded and unwounded mouse skin revealed that MCS-01 primarily altered the gene expression of mast cells, monocytes, and keratinocytes. Taken together, these findings offer insights into the process of diabetic wound healing and suggest topical mast cell stabilization as a potentially successful treatment for diabetic foot ulceration.