Single-cell spatial immune landscapes of primary and metastatic brain tumours.

Single-cell technologies have enabled the characterization of the tumour microenvironment at unprecedented depth and have revealed vast cellular diversity among tumour cells and their niche. Anti-tumour immunity relies on cell-cell relationships within the tumour microenvironment1,2, yet many single-cell studies lack spatial context and rely on dissociated tissues3. Here we applied imaging mass cytometry to characterize the immunological landscape of 139 high-grade glioma and 46 brain metastasis tumours from patients. Single-cell analysis of more than 1.1 million cells across 389 high-dimensional histopathology images enabled the spatial resolution of immune lineages and activation states, revealing differences in immune landscapes between primary tumours and brain metastases from diverse solid cancers. These analyses revealed cellular neighbourhoods associated with survival in patients with glioblastoma, which we leveraged to identify a unique population of myeloperoxidase (MPO)-positive macrophages associated with long-term survival. Our findings provide insight into the biology of primary and metastatic brain tumours, reinforcing the value of integrating spatial resolution to single-cell datasets to dissect the microenvironmental contexture of cancer.

Conserved expression of transposon-derived non-coding transcripts in primate stem cells.

BACKGROUND: A significant portion of expressed non-coding RNAs in human cells is derived from transposable elements (TEs). Moreover, it has been shown that various long non-coding RNAs (lncRNAs), which come from the human endogenous retrovirus subfamily H (HERVH), are not only expressed but required for pluripotency in human embryonic stem cells (hESCs). RESULTS: To identify additional TE-derived functional non-coding transcripts, we generated RNA-seq data from induced pluripotent stem cells (iPSCs) of four primate species (human, chimpanzee, gorilla, and rhesus) and searched for transcripts whose expression was conserved. We observed that about 30% of TE instances expressed in human iPSCs had orthologous TE instances that were also expressed in chimpanzee and gorilla. Notably, our analysis revealed a number of repeat families with highly conserved expression profiles including HERVH but also MER53, which is known to be the source of a placental-specific family of microRNAs (miRNAs). We also identified a number of repeat families from all classes of TEs, including MLT1-type and Tigger families, that contributed a significant amount of sequence to primate lncRNAs whose expression was conserved. CONCLUSIONS: Together, these results describe TE families and TE-derived lncRNAs whose conserved expression patterns can be used to identify what are likely functional TE-derived non-coding transcripts in primate iPSCs.

Transposable elements are major contributors to the origin, diversification, and regulation of vertebrate long noncoding RNAs.

Advances in vertebrate genomics have uncovered thousands of loci encoding long noncoding RNAs (lncRNAs). While progress has been made in elucidating the regulatory functions of lncRNAs, little is known about their origins and evolution. Here we explore the contribution of transposable elements (TEs) to the makeup and regulation of lncRNAs in human, mouse, and zebrafish. Surprisingly, TEs occur in more than two thirds of mature lncRNA transcripts and account for a substantial portion of total lncRNA sequence (~30% in human), whereas they seldom occur in protein-coding transcripts. While TEs contribute less to lncRNA exons than expected, several TE families are strongly enriched in lncRNAs. There is also substantial interspecific variation in the coverage and types of TEs embedded in lncRNAs, partially reflecting differences in the TE landscapes of the genomes surveyed. In human, TE sequences in lncRNAs evolve under greater evolutionary constraint than their non-TE sequences, than their intronic TEs, or than random DNA. Consistent with functional constraint, we found that TEs contribute signals essential for the biogenesis of many lncRNAs, including ~30,000 unique sites for transcription initiation, splicing, or polyadenylation in human. In addition, we identified ~35,000 TEs marked as open chromatin located within 10 kb upstream of lncRNA genes. The density of these marks in one cell type correlate with elevated expression of the downstream lncRNA in the same cell type, suggesting that these TEs contribute to cis-regulation. These global trends are recapitulated in several lncRNAs with established functions. Finally a subset of TEs embedded in lncRNAs are subject to RNA editing and predicted to form secondary structures likely important for function. In conclusion, TEs are nearly ubiquitous in lncRNAs and have played an important role in the lineage-specific diversification of vertebrate lncRNA repertoires.

Fecal microbiota transplantation plus anti-PD-1 immunotherapy in advanced melanoma: a phase I trial.

Fecal microbiota transplantation (FMT) represents a potential strategy to overcome resistance to immune checkpoint inhibitors in patients with refractory melanoma; however, the role of FMT in first-line treatment settings has not been evaluated. We conducted a multicenter phase I trial combining healthy donor FMT with the PD-1 inhibitors nivolumab or pembrolizumab in 20 previously untreated patients with advanced melanoma. The primary end point was safety. No grade 3 adverse events were reported from FMT alone. Five patients (25%) experienced grade 3 immune-related adverse events from combination therapy. Key secondary end points were objective response rate, changes in gut microbiome composition and systemic immune and metabolomics analyses. The objective response rate was 65% (13 of 20), including four (20%) complete responses. Longitudinal microbiome profiling revealed that all patients engrafted strains from their respective donors; however, the acquired similarity between donor and patient microbiomes only increased over time in responders. Responders experienced an enrichment of immunogenic and a loss of deleterious bacteria following FMT. Avatar mouse models confirmed the role of healthy donor feces in increasing anti-PD-1 efficacy. Our results show that FMT from healthy donors is safe in the first-line setting and warrants further investigation in combination with immune checkpoint inhibitors. ClinicalTrials.gov identifier NCT03772899 .

Melanomas with concurrent BRAF non-p.V600 and NF1 loss-of-function mutations are targetable by BRAF/MEK inhibitor combination therapy.

Although combination BRAF/MEK inhibition has produced significant survival benefits for BRAF p.V600 mutant melanomas, targeted therapies approved for BRAF non-p.V600 mutant melanomas remain limited. Through the analysis of 772 cutaneous melanoma exomes, we reveal that BRAF non-p.V600 mutations co-occurs more frequently with NF1 loss, but not with oncogenic NRAS mutations, than expected by chance. We present cell signaling data, which demonstrate that BRAF non-p.V600 mutants can signal as monomers and dimers within an NF1 loss context. Concordantly, BRAF inhibitors that inhibit both monomeric and dimeric BRAF synergize with MEK inhibition to significantly reduce cell viability in vitro and tumor growth in vivo in BRAF non-p.V600 mutant melanomas with co-occurring NF1 loss-of-function mutations. Our data suggest that patients harboring BRAF non-p.V600 mutant melanomas may benefit from current FDA-approved BRAF/MEK inhibitor combination therapy currently reserved for BRAF p.V600 mutant patients.

Spatially mapping the immune landscape of melanoma using imaging mass cytometry.

Melanoma is an immunogenic cancer with a high response rate to immune checkpoint inhibitors (ICIs). It harbors a high mutation burden compared with other cancers and, as a result, has abundant tumor-infiltrating lymphocytes (TILs) within its microenvironment. However, understanding the complex interplay between the stroma, tumor cells, and distinct TIL subsets remains a substantial challenge in immune oncology. To properly study this interplay, quantifying spatial relationships of multiple cell types within the tumor microenvironment is crucial. To address this, we used cytometry time-of-flight (CyTOF) imaging mass cytometry (IMC) to simultaneously quantify the expression of 35 protein markers, characterizing the microenvironment of 5 benign nevi and 67 melanomas. We profiled more than 220,000 individual cells to identify melanoma, lymphocyte subsets, macrophage/monocyte, and stromal cell populations, allowing for in-depth spatial quantification of the melanoma microenvironment. We found that within pretreatment melanomas, the abundance of proliferating antigen-experienced cytotoxic T cells (CD8+CD45RO+Ki67+) and the proximity of antigen-experienced cytotoxic T cells to melanoma cells were associated with positive response to ICIs. Our study highlights the potential of multiplexed single-cell technology to quantify spatial cell-cell interactions within the tumor microenvironment to understand immune therapy responses.

Mutations in the IFNγ-JAK-STAT Pathway Causing Resistance to Immune Checkpoint Inhibitors in Melanoma Increase Sensitivity to Oncolytic Virus Treatment.

PURPOSE: Next-generation sequencing studies and CRISPR-Cas9 screens have established mutations in the IFNγ-JAK-STAT pathway as an immune checkpoint inhibitor (ICI) resistance mechanism in a subset of patients with melanoma. We hypothesized ICI resistance mutations in the IFNγ pathway would simultaneously render melanomas susceptible to oncolytic virus (OV) therapy. EXPERIMENTAL DESIGN: Cytotoxicity experiments were performed with a number of OVs on a matched melanoma cell line pair generated from a baseline biopsy and a progressing lesion with complete JAK2 loss from a patient that relapsed on anti-PD-1 therapy, in melanoma lines following JAK1/2 RNA interference (RNAi) and pharmacologic inhibition and in Jak2 knockout (KO) B16-F10 mouse melanomas. Furthermore, we estimated the frequency of genetic alterations in the IFNγ-JAK-STAT pathway in human melanomas. RESULTS: The melanoma line from an anti-PD-1 progressing lesion was 7- and 22-fold more sensitive to the modified OVs, herpes simplex virus 1 (HSV1-dICP0) and vesicular stomatitis virus (VSV-Δ51), respectively, compared with the line from the baseline biopsy. RNAi, JAK1/2 inhibitor studies, and in vivo studies of Jak2 KOs B16-F10 melanomas revealed a significant increase in VSV-Δ51 sensitivity with JAK/STAT pathway inhibition. Our analysis of The Cancer Genome Atlas data estimated that approximately 11% of ICI-naïve cutaneous melanomas have alterations in IFNγ pathway genes that may confer OV susceptibility. CONCLUSIONS: We provide mechanistic support for the use of OVs as a precision medicine strategy for both salvage therapy in ICI-resistant and first-line treatment in melanomas with IFNγ-JAK-STAT pathway mutations. Our study also supports JAK inhibitor-OV combination therapy for treatment-naïve melanomas without IFN signaling defects.See related commentary by Kaufman, p. 3278.

C3a elicits unique migratory responses in immature low-density neutrophils.

Neutrophils represent the immune system's first line of defense and are rapidly recruited into inflamed tissue. In cancer associated inflammation, phenotypic heterogeneity has been ascribed to this cell type, whereby neutrophils can manifest anti- or pro-metastatic functions depending on the cellular/micro-environmental context. Here, we demonstrate that pro-metastatic immature low-density neutrophils (iLDNs) more efficiently accumulate in the livers of mice bearing metastatic lesions compared with anti-metastatic mature high-density neutrophils (HDNs). Transcriptomic analyses reveal enrichment of a migration signature in iLDNs relative to HDNs. We find that conditioned media derived from liver-metastatic breast cancer cells, but not lung-metastatic variants, specifically induces chemotaxis of iLDNs and not HDNs. Chemotactic responses are due to increased surface expression of C3aR in iLDNs relative to HDNs. In addition, we detect elevated secretion of cancer-cell derived C3a from liver-metastatic versus lung-metastatic breast cancer cells. Perturbation of C3a/C3aR signaling axis with either a small molecule inhibitor, SB290157, or reducing the levels of secreted C3a from liver-metastatic breast cancer cells by short hairpin RNAs, can abrogate the chemotactic response of iLDNs both in vitro and in vivo, respectively. Together, these data reveal novel mechanisms through which iLDNs prefentially accumulate in liver tissue harboring metastases in response to tumor-derived C3a secreted from the liver-aggressive 4T1 breast cancer cells.

Blood transcriptomic biomarker as a surrogate of ischemic brain gene expression.

OBJECTIVES: Blood biomarkers for cerebral tissue ischemia are lacking. The goal was to identify a blood transcriptomic signature jointly identified in the ischemic brain. METHODS: A nonhuman primate model with middle cerebral artery (MCA) territory infarction was used to study gene expression by microarray during acute ischemic cerebral stroke in the brain and the blood. Brain samples were collected in the infarcted and contralateral non-infarcted cortex as well as blood samples before and after occlusion. Gene expression was compared between the two brain locations to find differentially expressed genes. The expressions of these genes were then compared in the blood pre- and post-occlusion. RESULTS: Hierarchical clustering of brain expression data revealed strong independent clustering of ischemic and nonischemic brain samples. The top five enriched, up-regulated gene sets in the brain were TNF α signaling, apoptosis, P53 pathway, hypoxia, and UV response up. A comparison of differentially expressed genes in the brain and blood revealed a significant overlap of gene expression patterns. Stringent analysis of blood expression data from pre- and post-occlusion samples in each monkey identified nine genes highly differentially expressed in both the brain and the blood. Many of these up-regulated genes belong to pathways involved in cell death and DNA damage repair. INTERPRETATION: Common gene expression profile can be identified in the brain and blood and clearly differentiates ischemic from nonischemic conditions. Therefore, specific blood transcriptomic signature may represent a surrogate for brain ischemic gene expression.

Immature Low-Density Neutrophils Exhibit Metabolic Flexibility that Facilitates Breast Cancer Liver Metastasis.

Neutrophils are phenotypically heterogeneous and exert either anti- or pro-metastatic functions. We show that cancer-cell-derived G-CSF is necessary, but not sufficient, to mobilize immature low-density neutrophils (iLDNs) that promote liver metastasis. In contrast, mature high-density neutrophils inhibit the formation of liver metastases. Transcriptomic and metabolomic analyses of high- and low-density neutrophils reveal engagement of numerous metabolic pathways specifically in low-density neutrophils. iLDNs exhibit enhanced global bioenergetic capacity, through their ability to engage mitochondrial-dependent ATP production, and remain capable of executing pro-metastatic neutrophil functions, including NETosis, under nutrient-deprived conditions. We demonstrate that NETosis is an important neutrophil function that promotes breast cancer liver metastasis. iLDNs rely on the catabolism of glutamate and proline to support mitochondrial-dependent metabolism in the absence of glucose, which enables sustained NETosis. These data reveal that distinct pro-metastatic neutrophil populations exhibit a high degree of metabolic flexibility, which facilitates the formation of liver metastases.

In Silico Methods to Identify Exapted Transposable Element Families.

Transposable elements (TEs) have recently been shown to have many regulatory roles within the genome. In this chapter, we will examine two in silico methods for analyzing TEs and identifying families that may have acquired such functions. The first method will look at how the overrepresentation of a repeat family in a set of genomic features can be discovered. The example situation of OCT4 binding sites originating from LTR7 TE sequences will be used to show how this method could be applied. The second method will describe how to determine if a TE family exhibits a cell type-specific expression pattern. As an example, we will look at the expression of HERV-H, an endogenous retrovirus known to act as an lncRNA in embryonic stem cells. We will use this example to demonstrate how RNA-seq data can be used to compare cell type expression of repeats.

The retrovirus HERVH is a long noncoding RNA required for human embryonic stem cell identity.

Human endogenous retrovirus subfamily H (HERVH) is a class of transposable elements expressed preferentially in human embryonic stem cells (hESCs). Here, we report that the long terminal repeats of HERVH function as enhancers and that HERVH is a nuclear long noncoding RNA required to maintain hESC identity. Furthermore, HERVH is associated with OCT4, coactivators and Mediator subunits. Together, these results uncover a new role of species-specific transposable elements in hESCs.