Autoreactive CD8+ T cells are restrained by an exhaustion-like program that is maintained by LAG3.

Impaired chronic viral and tumor clearance has been attributed to CD8+ T cell exhaustion, a differentiation state in which T cells have reduced and altered effector function that can be partially reversed upon blockade of inhibitory receptors. The role of the exhaustion program and transcriptional networks that control CD8+ T cell function and fate in autoimmunity is not clear. Here we show that intra-islet CD8+ T cells phenotypically, transcriptionally, epigenetically and metabolically possess features of canonically exhausted T cells, yet maintain important differences. This 'restrained' phenotype can be perturbed and disease accelerated by CD8+ T cell-restricted deletion of the inhibitory receptor lymphocyte activating gene 3 (LAG3). Mechanistically, LAG3-deficient CD8+ T cells have enhanced effector-like functions, trafficking to the islets, and have a diminished exhausted phenotype, highlighting a physiological role for an exhaustion program in limiting autoimmunity and implicating LAG3 as a target for autoimmune therapy.

Microbiota-specific T follicular helper cells drive tertiary lymphoid structures and anti-tumor immunity against colorectal cancer.

The composition of the intestinal microbiota is associated with both the development of tumors and the efficacy of anti-tumor immunity. Here, we examined the impact of microbiota-specific T cells in anti-colorectal cancer (CRC) immunity. Introduction of Helicobacter hepaticus (Hhep) in a mouse model of CRC did not alter the microbial landscape but increased tumor infiltration by cytotoxic lymphocytes and inhibited tumor growth. Anti-tumor immunity was independent of CD8+ T cells but dependent upon CD4+ T cells, B cells, and natural killer (NK) cells. Hhep colonization induced Hhep-specific T follicular helper (Tfh) cells, increased the number of colon Tfh cells, and supported the maturation of Hhep+ tumor-adjacent tertiary lymphoid structures. Tfh cells were necessary for Hhep-mediated tumor control and immune infiltration, and adoptive transfer of Hhep-specific CD4+ T cells to Tfh cell-deficient Bcl6fl/flCd4Cre mice restored anti-tumor immunity. Thus, introduction of immunogenic intestinal bacteria can promote Tfh-associated anti-tumor immunity in the colon, suggesting therapeutic approaches for the treatment of CRC.

Immune Landscape of Viral- and Carcinogen-Driven Head and Neck Cancer.

Head and neck squamous cell carcinoma (HNSCC) arises through exposure to environmental carcinogens or malignant transformation by human papillomavirus (HPV). Here, we assessed the transcriptional profiles of 131,224 single cells from peripheral and intra-tumoral immune populations from patients with HPV- and HPV+ HNSCC and healthy donors. Immune cells within tumors of HPV- and HPV+ HNSCC displayed a spectrum of transcriptional signatures, with helper CD4+ T cells and B cells being relatively divergent and CD8+ T cells and CD4+ regulatory T cells being relatively similar. Transcriptional results were contextualized through multispectral immunofluorescence analyses and evaluating putative cell-cell communication based on spatial proximity. These analyses defined a gene expression signature associated with CD4+ T follicular helper cells that is associated with longer progression-free survival in HNSCC patients. The datasets and analytical approaches herein provide a resource for the further study of the impact of immune cells on viral- and carcinogen-induced cancers.

Single-cell analysis of HIV-1 transcriptional activity reveals expression of proviruses in expanded clones during ART.

Little is known about the fraction of human immunodeficiency virus type 1 (HIV-1) proviruses that express unspliced viral RNA in vivo or about the levels of HIV RNA expression within single infected cells. We developed a sensitive cell-associated HIV RNA and DNA single-genome sequencing (CARD-SGS) method to investigate fractional proviral expression of HIV RNA (1.3-kb fragment of p6, protease, and reverse transcriptase) and the levels of HIV RNA in single HIV-infected cells from blood samples obtained from individuals with viremia or individuals on long-term suppressive antiretroviral therapy (ART). Spiking experiments show that the CARD-SGS method can detect a single cell expressing HIV RNA. Applying CARD-SGS to blood mononuclear cells in six samples from four HIV-infected donors (one with viremia and not on ART and three with viremia suppressed on ART) revealed that an average of 7% of proviruses (range: 2-18%) expressed HIV RNA. Levels of expression varied from one to 62 HIV RNA molecules per cell (median of 1). CARD-SGS also revealed the frequent expression of identical HIV RNA sequences across multiple single cells and across multiple time points in donors on suppressive ART consistent with constitutive expression of HIV RNA in infected cell clones. Defective proviruses were found to express HIV RNA at levels similar to those proviruses that had no obvious defects. CARD-SGS is a useful tool to characterize fractional proviral expression in single infected cells that persist despite ART and to assess the impact of experimental interventions on proviral populations and their expression.

HIV-1 Virion Production from Single Inducible Proviruses following T-Cell Activation Ex Vivo.

Quantifying induced virion production from single proviruses is important for assessing the effects of HIV-1 latency reversal agents. Limiting dilution ex vivo cultures of resting CD4(+) T cells from 14 HIV-positive volunteers revealed that virion production after T-cell activation from individual proviruses varies by 10,000-fold to 100,000-fold. High-producing proviruses were associated with increases in cell-associated HIV-1 DNA levels, suggesting that reactivated proviruses proliferate. Single-cell analyses are needed to investigate differences in proviral expansion and virus production following latency reversal.

Structured inquiry-based learning: Drosophila GAL4 enhancer trap characterization in an undergraduate laboratory course.

We have developed and tested two linked but separable structured inquiry exercises using a set of Drosophila melanogaster GAL4 enhancer trap strains for an upper-level undergraduate laboratory methods course at Bucknell University. In the first, students learn to perform inverse PCR to identify the genomic location of the GAL4 insertion, using FlyBase to identify flanking sequences and the primary literature to synthesize current knowledge regarding the nearest gene. In the second, we cross each GAL4 strain to a UAS-CD8-GFP reporter strain, and students perform whole mount CNS dissection, immunohistochemistry, confocal imaging, and analysis of developmental expression patterns. We have found these exercises to be very effective in teaching the uses and limitations of PCR and antibody-based techniques as well as critical reading of the primary literature and scientific writing. Students appreciate the opportunity to apply what they learn by generating novel data of use to the wider research community.

Quantification of HIV-1 latency reversal in resting CD4+ T cells from patients on suppressive antiretroviral therapy.

Reversal of proviral latency is being pursued as a curative strategy for HIV-1 infection. Recent clinical studies of in vivo administration of the histone deacetylase inhibitor suberoylanilide hydroxamic acid (SAHA; vorinostat) show increases in unspliced cellular HIV-1 RNA levels in resting CD4(+) T cells. A critical unknown, however, is the proportion of latent proviruses that can be transcriptionally reactivated by SAHA or T-cell activation. In this study, we quantified the fraction of HIV-1 proviruses in resting CD4(+) T cells from patients on suppressive antiretroviral therapy that were reactivated ex vivo with SAHA or antibodies to CD3/CD28. At concentrations of SAHA achieved clinically, only 0.079% of proviruses in resting CD4(+) T cells were reactivated to produce virions, compared with 1.5% of proviruses in cells treated with anti-CD3/CD28 antibodies after correcting for spontaneous virion production in the medium control. A significant positive correlation (ρ = 0.67, P < 0.001) was found between levels of virions in the supernatant and unspliced cellular HIV-1 RNA following anti-CD3/CD28 treatment, but not following SAHA treatment (ρ = 0.21, P = 0.99). These results reveal that the majority of HIV-1 proviruses are not reactivated by current therapeutic approaches and that more effective means of reversing proviral latency will likely be required to deplete HIV-1 reservoirs.

Novel Assays for Measurement of Total Cell-Associated HIV-1 DNA and RNA.

Although a number of PCR-based quantitative assays for measuring HIV-1 persistence during suppressive antiretroviral therapy (ART) have been reported, a simple, sensitive, reproducible method is needed for application to large clinical trials. We developed novel quantitative PCR assays for cell-associated (CA) HIV-1 DNA and RNA, targeting a highly conserved region in HIV-1pol, with sensitivities of 3 to 5 copies/1 million cells. We evaluated the performance characteristics of the assays using peripheral blood mononuclear cells (PBMCs) from 5 viremic patients and 20 patients receiving effective ART. Total and resting CD4(+)T cells were isolated from a subset of patients and tested for comparison with PBMCs. The estimated standard deviations including interassay variability and intra-assay variability of the assays were modest, i.e., 0.15 and 0.10 log10copies/10(6)PBMCs, respectively, for CA HIV-1 DNA and 0.40 and 0.19 log10copies/10(6)PBMCs for CA HIV-1 RNA. Testing of longitudinally obtained PBMC samples showed little variation for either viremic patients (median fold differences of 0.80 and 0.88 for CA HIV-1 DNA and RNA, respectively) or virologically suppressed patients (median fold differences of 1.14 and 0.97, respectively). CA HIV-1 DNA and RNA levels were strongly correlated (r= 0.77 to 1;P= 0.0001 to 0.037) for assays performed using PBMCs from different sources (phlebotomy versus leukapheresis) or using total or resting CD4(+)T cells purified by either bead selection or flow cytometric sorting. Their sensitivity, reproducibility, and broad applicability to small numbers of mononuclear cells make these assays useful for observational and interventional studies that examine longitudinal changes in the numbers of HIV-1-infected cells and their levels of transcription.

Simian Immunodeficiency Virus SIVsab Infection of Rhesus Macaques as a Model of Complete Immunological Suppression with Persistent Reservoirs of Replication-Competent Virus: Implications for Cure Research.

Simian immunodeficiency virus SIVsab infection is completely controlled in rhesus macaques (RMs) through functional immune responses. We report that in SIVsab-infected RMs, (i) viral replication is controlled to <0 to 3 copies/ml, (ii) about one-third of the virus strains in reservoirs are replication incompetent, and (iii) rebounding virus after CD8(+) cell depletion is replication competent and genetically similar to the original virus stock, suggesting early reservoir seeding. This model permits assessment of strategies aimed at depleting the reservoir without multidrug antiretroviral therapy.

Improved single-copy assays for quantification of persistent HIV-1 viremia in patients on suppressive antiretroviral therapy.

A quantitative real-time PCR (qRT-PCR) assay with single-copy sensitivity targeting HIV-1 gag RNA (the gag single-copy assay [gSCA]) has been used widely to quantify plasma viremia below the limit of detection of clinical assays in patients on effective antiretroviral therapy (ART), but viral RNA in 15 to 30% of samples amplifies inefficiently because of primer/probe mismatches. We sought to develop improved single-copy assays with increased sensitivity by improving nucleic acid recovery, designing qRT-PCR primers and a probe for a highly conserved region of integrase in the HIV-1 pol gene (the integrase single-copy assay [iSCA]), and increasing the plasma volume tested (Mega-iSCA). We evaluated gSCA versus iSCA in paired plasma samples from 10 consecutive patients with viremia of >1,000 copies/ml and 25 consecutive patients on suppressive ART. Three of 10 viremic samples amplified inefficiently with gSCA compared to the Roche Cobas Ampliprep/TaqMan 2.0, whereas all 10 samples amplified efficiently with iSCA. Among 25 samples from patients on suppressive ART, 8 of 12 samples that were negative for HIV-1 RNA by gSCA had detectable HIV-1 RNA by iSCA, and iSCA detected 3-fold or higher HIV-1 RNA levels compared to gSCA in 10 of 25 samples. Large-volume plasma samples (>20 ml) from 7 patients were assayed using Mega-iSCA, and HIV-1 RNA was quantifiable in 6, including 4 of 5 that were negative by standard-volume iSCA. These improved assays with superior sensitivity will be useful for evaluating whether in vivo interventions can reduce plasma viremia and for assessing relationships between residual viremia and other virologic parameters, including the inducible proviral reservoir.

Three distinct phases of HIV-1 RNA decay in treatment-naive patients receiving raltegravir-based antiretroviral therapy: ACTG A5248.

OBJECTIVE: The goal of this study was to define viral kinetics after initiation of raltegravir (RAL)-based antiretroviral therapy (ART). METHODS: ART-naive patients received RAL, tenofovir disoproxil fumarate, and emtricitabine for 72 weeks. Human immunodeficiency virus type 1 (HIV-1) RNA were measured by ultrasensitive and single-copy assays, and first (d1)-, second (d2)-, and, third (d3)-phase decay rates were estimated by mixed-effects models. Decay data were compared to historical estimates for efavirenz (EFV)- and ritonavir/lopinavir (LPV/r)-based regimens. RESULTS: Bi- and tri-exponential models for ultrasensitive assay (n = 38) and single-copy assay (n = 8) data, respectively, provided the best fits over 8 and 72 weeks. The median d1 with ultrasensitive data was 0.563/day (interquartile range [IQR], 0.501-0.610/day), significantly slower than d1 for EFV-based regimens [P < .001]). The median duration of d1 was 15.1 days, transitioning to d2 at an HIV-1 RNA of 91 copies/mL, indicating a longer duration of d1 and a d2 transition at lower viremia levels than with EFV. Median patient-specific decay estimates with the single-copy assay were 0.607/day (IQR, 0.582-0.653) for d1, 0.070/day (IQR, 0.042-0.079) for d2, and 0.0016/day (IQR, 0.0005-0.0022) for d3; the median d1 duration was 16.1 days, transitioning to d2 at 69 copies/mL. d3 transition occurred at 110 days, at 2.6 copies/mL, similar to values for LPV/r-based regimens. CONCLUSIONS: Models using single-copy assay data revealed 3 phases of decay with RAL-containing ART, with a longer duration of first-phase decay consistent with RAL-mediated blockade of productive infection from preintegration complexes.

Immunocompetent murine model of Ewing sarcoma reveals role for TGFß inhibition to enhance immune infiltrates in Ewing tumors during radiation.

Ewing sarcoma (ES) is an aggressive cancer diagnosed in adolescents and young adults. The fusion oncoprotein (EWSR1::FLI1) that drives Ewing sarcoma is known to downregulate TGFBR2 expression (part of the TGFß receptor). Because TGFBR2 is downregulated, it was thought that TGFß likely plays an inconsequential role in Ewing biology. However, the expression of TGFß in the Ewing tumor immune microenvironment (TIME) and functional impact of TGFß in the TIME remains largely unknown given the historical lack of immunocompetent preclinical models. Here, we use single-cell RNAseq analysis of human Ewing tumors to show that immune cells, such as NK cells, are the largest source of TGFß production in human Ewing tumors. We develop a humanized (immunocompetent) mouse model of ES and demonstrate distinct TME signatures and metastatic potential in these models as compared to tumors developed in immunodeficient mice. Using this humanized model, we study the effect of TGFß inhibition on the Ewing TME during radiation therapy, a treatment that both enhances TGFß activation and is used to treat aggressive ES. Utilizing a trivalent ligand TGFß TRAP to inhibit TGFß, we demonstrate that in combination with radiation, TGFß inhibition both increases ES immune cell infiltration and decreases lung metastatic burden in vivo . The culmination of these data demonstrates the value of humanized models to address immunobiologic preclinical questions in Ewing sarcoma and suggests TGFß inhibition as a promising intervention during radiation therapy to promote metastatic tumor control.

All-trans retinoic acid induces durable tumor immunity in IDH-mutant gliomas by rescuing transcriptional repression of the CRBP1-retinoic acid axis.

Diffuse gliomas are epigenetically dysregulated, immunologically cold, and fatal tumors characterized by mutations in isocitrate dehydrogenase (IDH). Although IDH mutations yield a uniquely immunosuppressive tumor microenvironment, the regulatory mechanisms that drive the immune landscape of IDH mutant (IDHm) gliomas remain unknown. Here, we reveal that transcriptional repression of retinoic acid (RA) pathway signaling impairs both innate and adaptive immune surveillance in IDHm glioma through epigenetic silencing of retinol binding protein 1 (RBP1) and induces a profound anti-inflammatory landscape marked by loss of inflammatory cell states and infiltration of suppressive myeloid phenotypes. Restorative retinoic acid therapy in murine glioma models promotes clonal CD4 + T cell expansion and induces tumor regression in IDHm, but not IDH wildtype (IDHwt), gliomas. Our findings provide a mechanistic rationale for RA immunotherapy in IDHm glioma and is the basis for an ongoing investigator-initiated, single-center clinical trial investigating all-trans retinoic acid (ATRA) in recurrent IDHm human subjects.

Integrated BATF transcriptional network regulates suppressive intratumoral regulatory T cells.

Human regulatory T cells (Tregs) are crucial regulators of tissue repair, autoimmune diseases, and cancer. However, it is challenging to inhibit the suppressive function of Tregs for cancer therapy without affecting immune homeostasis. Identifying pathways that may distinguish tumor-restricted Tregs is important, yet the transcriptional programs that control intratumoral Treg gene expression, and that are distinct from Tregs in healthy tissues, remain largely unknown. We profiled single-cell transcriptomes of CD4+ T cells in tumors and peripheral blood from patients with head and neck squamous cell carcinomas (HNSCC) and those in nontumor tonsil tissues and peripheral blood from healthy donors. We identified a subpopulation of activated Tregs expressing multiple tumor necrosis factor receptor (TNFR) genes (TNFR+ Tregs) that is highly enriched in the tumor microenvironment (TME) compared with nontumor tissue and the periphery. TNFR+ Tregs are associated with worse prognosis in HNSCC and across multiple solid tumor types. Mechanistically, the transcription factor BATF is a central component of a gene regulatory network that governs key aspects of TNFR+ Tregs. CRISPR-Cas9-mediated BATF knockout in human activated Tregs in conjunction with bulk RNA sequencing, immunophenotyping, and in vitro functional assays corroborated the central role of BATF in limiting excessive activation and promoting the survival of human activated Tregs. Last, we identified a suite of surface molecules reflective of the BATF-driven transcriptional network on intratumoral Tregs in patients with HNSCC. These findings uncover a primary transcriptional regulator of highly suppressive intratumoral Tregs, highlighting potential opportunities for therapeutic intervention in cancer without affecting immune homeostasis.

Immune landscape in invasive ductal and lobular breast cancer reveals a divergent macrophage-driven microenvironment.

T cell-centric immunotherapies have shown modest clinical benefit thus far for estrogen receptor-positive (ER+) breast cancer. Despite accounting for 70% of all breast cancers, relatively little is known about the immunobiology of ER+ breast cancer in women with invasive ductal carcinoma (IDC) and invasive lobular carcinoma (ILC). To investigate this, we performed phenotypic, transcriptional and functional analyses for a cohort of treatment-naive IDC (n = 94) and ILC (n = 87) tumors. We show that macrophages, and not T cells, are the predominant immune cells infiltrating the tumor bed and the most transcriptionally diverse cell subset between IDC and ILC. Analysis of cellular neighborhoods revealed an interplay between macrophages and T cells associated with longer disease-free survival in IDC but not ILC. Our datasets provide a rich resource for further interrogation into immune cell dynamics in ER+ IDC and ILC and highlight macrophages as a potential target for ER+ breast cancer.

New tools for quantifying HIV-1 reservoirs: plasma RNA single copy assays and beyond.

Quantification of plasma HIV-1 RNA below the limit of FDA-approved assays by a single copy quantitative PCR assays (SCA) has provided significant insights into HIV-1 persistence despite potent antiretroviral therapy as well as a means to assess the impact of therapeutic strategies, such as treatment intensification, on residual viremia. In this review, we discuss insights gained from plasma HIV-1 RNA SCA and highlight the need for additional assays to characterize better the cellular and tissue reservoirs of HIV-1. Accurate, reproducible, and sensitive assays to quantify HIV-1 reservoirs, before and after therapeutic interventions, are essential tools in the quest for a cure of HIV-1 infection.

Molecular Pathways and Mechanisms of LAG3 in Cancer Therapy.

Immunotherapy targeting coinhibitory receptors has been highly successful in treating a wide variety of malignancies; however, only a subset of patients exhibits durable responses. The first FDA-approved immunotherapeutics targeting coinhibitory receptors PD1 and CTLA4, alone or in combination, significantly improved survival but were also accompanied by substantial toxicity in combination. The third FDA-approved immune checkpoint inhibitor targets LAG3, a coinhibitory receptor expressed on activated CD4+ and CD8+ T cells, especially in settings of long-term antigenic stimulation, such as chronic viral infection or cancer. Mechanistically, LAG3 expression limits both the expansion of activated T cells and the size of the memory pool, suggesting that LAG3 may be a promising target for immunotherapy. Importantly, the mechanism(s) by which LAG3 contributes to CD8+ T-cell exhaustion may be distinct from those governed by PD1, indicating that the combination of anti-LAG3 and anti-PD1 may synergistically enhance antitumor immunity. Clinical studies evaluating the role of anti-LAG3 in combination with anti-PD1 are underway, and recent phase III trial results in metastatic melanoma demonstrate both the efficacy and safety of this combination. Further ongoing clinical trials are evaluating this combination across multiple tumor types and the adjuvant setting, with accompanying translational and biomarker-focused studies designed to elucidate the molecular pathways that lead to improved antitumor T-cell responses following dual blockade of PD1 and LAG3. Overall, LAG3 plays an important role in limiting T-cell activation and has now become part of the repertoire of combinatorial immunotherapeutics available for the treatment of metastatic melanoma.

NKG7 Is Required for Optimal Antitumor T-cell Immunity.

Tumor antigen-specific CD8+ T cells play a critical role in antitumor immunity. Clinical trials reinvigorating the immune system via immune checkpoint blockade (ICB) have shown remarkable clinical promise. Numerous studies have identified an association between NKG7 expression and patient outcome across different malignancies. However, aside from these correlative observations, very little is known about NKG7 and its role in antitumor immunity. Herein, we utilized single-cell RNA sequencing (scRNA-seq) datasets, NKG7-deficient mice, NKG7-reporter mice, and mouse tumor models to investigate the role of NKG7 in neoantigen-mediated tumor rejection and ICB immunotherapy. scRNA-seq of tumors from patients with metastatic melanoma or head and neck squamous cell carcinoma revealed that NKG7 expression is highly associated with cytotoxicity and specifically expressed by CD8+ T cells and natural killer (NK) cells. Furthermore, we identified a key role for NKG7 in controlling intratumor T-cell accumulation and activation. NKG7 was upregulated on intratumor antigen-specific CD8+ T cells and NK cells and required for the accumulation of T cells in the tumor microenvironment. Accordingly, neoantigen-expressing mouse tumors grew faster in Nkg7-deficient mice. Strikingly, efficacy of single or combination ICB was significantly reduced in Nkg7-deficient mice.See related article by Wen et al., p. 162.

Mobilizing phospholipids on tumor plasma membrane implicates phosphatidylserine externalization blockade for cancer immunotherapy.

In "healthy" tumor cells, phosphatidylserine (PS) is predominately localized in the inner plasma membrane leaflet. During apoptosis, PS relocates to the outer leaflet. Herein, we established PSout tumor models with tumor cells lacking PS flippase component CDC50A, constantly exposing PS but alive. PSout tumors developed bigger than wild-type (WT) tumors, featuring M2 polarized tumor-associated macrophages (TAMs) and fewer tumor-antigen-specific T cells. The PS receptor TIM-3 is responsible for PS recognition. Employing an opposite tumor model, PSin, with tumor cells lacking the PS scramblase Xkr8 and unable to expose PS during otherwise normal apoptosis, we find that the accumulated apoptotic tumor cells produce and release cyclic GAMP (cGAMP) to immune cells to activate the STING pathway, leading to TAM M1 polarization, suppressed interleukin (IL)-10 secretion, and natural killer (NK) cell cytotoxicity. Silencing Xkr8 in vivo by either short hairpin RNA (shRNA) or small interfering RNA (siRNA) to achieve a PS externalization blockade provides robust therapeutic anti-tumor efficiency.

Ewing Sarcoma and Osteosarcoma Have Distinct Immune Signatures and Intercellular Communication Networks.

PURPOSE: Ewing sarcoma and osteosarcoma are primary bone sarcomas occurring most commonly in adolescents. Metastatic and relapsed disease are associated with dismal prognosis. Although effective for some soft tissue sarcomas, current immunotherapeutic approaches for the treatment of bone sarcomas have been largely ineffective, necessitating a deeper understanding of bone sarcoma immunobiology. EXPERIMENTAL DESIGN: Multiplex immunofluorescence analysis of immune infiltration in relapsed versus primary disease was conducted. To better understand immune states and drivers of immune infiltration, especially during disease progression, we performed single-cell RNA sequencing (scRNAseq) of immune populations from paired blood and bone sarcoma tumor samples. RESULTS: Our multiplex immunofluorescence analysis revealed increased immune infiltration in relapsed versus primary disease in both Ewing sarcoma and osteosarcoma. scRNAseq analyses revealed terminally exhausted CD8+ T cells expressing co-inhibitory receptors in osteosarcoma and an effector T-cell subpopulation in Ewing sarcoma. In addition, distinct subsets of CD14+CD16+ macrophages were present in Ewing sarcoma and osteosarcoma. To determine pathways driving tumor immune infiltration, we conducted intercellular communication analyses and uncovered shared mechanisms of immune infiltration driven by CD14+CD16+ macrophages and unique pathways of immune infiltration driven by CXCL10 and CXCL12 in osteosarcoma. CONCLUSIONS: Our study provides preclinical rationale for future investigation of specific immunotherapeutic targets upon relapse and provides an invaluable resource of immunologic data from bone sarcomas.