Time-resolved single-cell transcriptomics defines immune trajectories in glioblastoma.

Deciphering the cell-state transitions underlying immune adaptation across time is fundamental for advancing biology. Empirical in vivo genomic technologies that capture cellular dynamics are currently lacking. We present Zman-seq, a single-cell technology recording transcriptomic dynamics across time by introducing time stamps into circulating immune cells, tracking them in tissues for days. Applying Zman-seq resolved cell-state and molecular trajectories of the dysfunctional immune microenvironment in glioblastoma. Within 24 hours of tumor infiltration, cytotoxic natural killer cells transitioned to a dysfunctional program regulated by TGFB1 signaling. Infiltrating monocytes differentiated into immunosuppressive macrophages, characterized by the upregulation of suppressive myeloid checkpoints Trem2, Il18bp, and Arg1, over 36 to 48 hours. Treatment with an antagonistic anti-TREM2 antibody reshaped the tumor microenvironment by redirecting the monocyte trajectory toward pro-inflammatory macrophages. Zman-seq is a broadly applicable technology, enabling empirical measurements of differentiation trajectories, which can enhance the development of more efficacious immunotherapies.

Coupled scRNA-Seq and Intracellular Protein Activity Reveal an Immunosuppressive Role of TREM2 in Cancer.

Cell function and activity are regulated through integration of signaling, epigenetic, transcriptional, and metabolic pathways. Here, we introduce INs-seq, an integrated technology for massively parallel recording of single-cell RNA sequencing (scRNA-seq) and intracellular protein activity. We demonstrate the broad utility of INs-seq for discovering new immune subsets by profiling different intracellular signatures of immune signaling, transcription factor combinations, and metabolic activity. Comprehensive mapping of Arginase 1-expressing cells within tumor models, a metabolic immune signature of suppressive activity, discovers novel Arg1+ Trem2+ regulatory myeloid (Mreg) cells and identifies markers, metabolic activity, and pathways associated with these cells. Genetic ablation of Trem2 in mice inhibits accumulation of intra-tumoral Mreg cells, leading to a marked decrease in dysfunctional CD8+ T cells and reduced tumor growth. This study establishes INs-seq as a broadly applicable technology for elucidating integrated transcriptional and intra-cellular maps and identifies the molecular signature of myeloid suppressive cells in tumors.

Dysfunctional CD8 T Cells Form a Proliferative, Dynamically Regulated Compartment within Human Melanoma.

Tumor immune cell compositions play a major role in response to immunotherapy, but the heterogeneity and dynamics of immune infiltrates in human cancer lesions remain poorly characterized. Here, we identify conserved intratumoral CD4 and CD8 T cell behaviors in scRNA-seq data from 25 melanoma patients. We discover a large population of CD8 T cells showing continuous progression from an early effector "transitional" into a dysfunctional T cell state. CD8 T cells that express a complete cytotoxic gene set are rare, and TCR sharing data suggest their independence from the transitional and dysfunctional cell states. Notably, we demonstrate that dysfunctional T cells are the major intratumoral proliferating immune cell compartment and that the intensity of the dysfunctional signature is associated with tumor reactivity. Our data demonstrate that CD8 T cells previously defined as exhausted are in fact a highly proliferating, clonal, and dynamically differentiating cell population within the human tumor microenvironment.

Dissecting Immune Circuits by Linking CRISPR-Pooled Screens with Single-Cell RNA-Seq.

In multicellular organisms, dedicated regulatory circuits control cell type diversity and responses. The crosstalk and redundancies within these circuits and substantial cellular heterogeneity pose a major research challenge. Here, we present CRISP-seq, an integrated method for massively parallel single-cell RNA sequencing (RNA-seq) and clustered regularly interspaced short palindromic repeats (CRISPR)-pooled screens. We show that profiling the genomic perturbation and transcriptome in the same cell enables us to simultaneously elucidate the function of multiple factors and their interactions. We applied CRISP-seq to probe regulatory circuits of innate immunity. By sampling tens of thousands of perturbed cells in vitro and in mice, we identified interactions and redundancies between developmental and signaling-dependent factors. These include opposing effects of Cebpb and Irf8 in regulating the monocyte/macrophage versus dendritic cell lineages and differential functions for Rela and Stat1/2 in monocyte versus dendritic cell responses to pathogens. This study establishes CRISP-seq as a broadly applicable, comprehensive, and unbiased approach for elucidating mammalian regulatory circuits.

Author Correction: Genome-wide SWAp-Tag yeast libraries for proteome exploration.

The version of Supplementary Table 1 originally published online with this article contained incorrect localization annotations for one plate. This error has been corrected in the online Supplementary Information.

Evolution of intron splicing towards optimized gene expression is based on various Cis- and Trans-molecular mechanisms.

Splicing expands, reshapes, and regulates the transcriptome of eukaryotic organisms. Despite its importance, key questions remain unanswered, including the following: Can splicing evolve when organisms adapt to new challenges? How does evolution optimize inefficiency of introns' splicing and of the splicing machinery? To explore these questions, we evolved yeast cells that were engineered to contain an inefficiently spliced intron inside a gene whose protein product was under selection for an increased expression level. We identified a combination of mutations in Cis (within the gene of interest) and in Trans (in mRNA-maturation machinery). Surprisingly, the mutations in Cis resided outside of known intronic functional sites and improved the intron's splicing efficiency potentially by easing tight mRNA structures. One of these mutations hampered a protein's domain that was not under selection, demonstrating the evolutionary flexibility of multi-domain proteins as one domain functionality was improved at the expense of the other domain. The Trans adaptations resided in two proteins, Npl3 and Gbp2, that bind pre-mRNAs and are central to their maturation. Interestingly, these mutations either increased or decreased the affinity of these proteins to mRNA, presumably allowing faster spliceosome recruitment or increased time before degradation of the pre-mRNAs, respectively. Altogether, our work reveals various mechanistic pathways toward optimizations of intron splicing to ultimately adapt gene expression patterns to novel demands.

Genome-wide SWAp-Tag yeast libraries for proteome exploration.

Yeast libraries revolutionized the systematic study of cell biology. To extensively increase the number of such libraries, we used our previously devised SWAp-Tag (SWAT) approach to construct a genome-wide library of ~5,500 strains carrying the SWAT NOP1promoter-GFP module at the N terminus of proteins. In addition, we created six diverse libraries that restored the native regulation, created an overexpression library with a Cherry tag, or enabled protein complementation assays from two fragments of an enzyme or fluorophore. We developed methods utilizing these SWAT collections to systematically characterize the yeast proteome for protein abundance, localization, topology, and interactions.

Melanoma-derived extracellular vesicles instigate proinflammatory signaling in the metastatic microenvironment.

The major cause of melanoma mortality is metastasis to distant organs, including lungs and brain. Reciprocal interactions of metastasizing tumor cells with stromal cells in secondary sites play a critical role in all stages of tumorigenesis and metastasis. Changes in the metastatic microenvironment were shown to precede clinically relevant metastases, and may occur prior to the arrival of disseminated tumor cells to the distant organ, thus creating a hospitable "premetastatic niche." Exosomes secreted by tumor cells were demonstrated to play an important role in the preparation of a hospitable metastatic niche. However, the functional role of melanoma-derived exosomes on metastatic niche formation, and the downstream pathways activated in stromal cells at the metastatic niche are largely unresolved. Here we show that extracellular vesicles (EVs) secreted by metastatic melanoma cells that spontaneously metastasize to lungs and to brain, activate proinflammatory signaling in lung fibroblasts and in astrocytes. Interestingly, unlike paracrine signaling by melanoma cells, EVs secreted by metastatic melanoma cells instigated a proinflammatory gene signature in lung fibroblasts but did not activate wound-healing functions, suggesting that tumor cell-secreted EVs activate distinct CAF characteristics and tumor-promoting functions. Moreover, melanoma-secreted EVs also activated proinflammatory signaling in astrocytes, indicating that EV-mediated reprogramming of stromal cells is a general mechanism of modulating the metastatic niche in multiple distant organs. Thus, our study demonstrates that melanoma-derived EVs reprogram tumor-promoting functions in stromal cells in a distinct manner, implicating a central role for tumor-derived EV signaling in promoting the formation of an inflammatory metastatic niche.

Pex35 is a regulator of peroxisome abundance.

Peroxisomes are cellular organelles with vital functions in lipid, amino acid and redox metabolism. The cellular formation and dynamics of peroxisomes are governed by PEX genes; however, the regulation of peroxisome abundance is still poorly understood. Here, we use a high-content microscopy screen in Saccharomyces cerevisiae to identify new regulators of peroxisome size and abundance. Our screen led to the identification of a previously uncharacterized gene, which we term PEX35, which affects peroxisome abundance. PEX35 encodes a peroxisomal membrane protein, a remote homolog to several curvature-generating human proteins. We systematically characterized the genetic and physical interactome as well as the metabolome of mutants in PEX35, and we found that Pex35 functionally interacts with the vesicle-budding-inducer Arf1. Our results highlight the functional interaction between peroxisomes and the secretory pathway.

One library to make them all: streamlining the creation of yeast libraries via a SWAp-Tag strategy.

The yeast Saccharomyces cerevisiae is ideal for systematic studies relying on collections of modified strains (libraries). Despite the significance of yeast libraries and the immense variety of available tags and regulatory elements, only a few such libraries exist, as their construction is extremely expensive and laborious. To overcome these limitations, we developed a SWAp-Tag (SWAT) method that enables one parental library to be modified easily and efficiently to give rise to an endless variety of libraries of choice. To showcase the versatility of the SWAT approach, we constructed and investigated a library of ∼1,800 strains carrying SWAT-GFP modules at the amino termini of endomembrane proteins and then used it to create two new libraries (mCherry and seamless GFP). Our work demonstrates how the SWAT method allows fast and effortless creation of yeast libraries, opening the door to new ways of systematically studying cell biology.

Characterization of proteome dynamics during growth in oleate reveals a new peroxisome-targeting receptor.

To optimally perform the diversity of metabolic functions that occur within peroxisomes, cells must dynamically regulate peroxisome size, number and content in response to the cell state and the environment. Except for transcriptional regulation little is known about the mechanisms used to perform this complicated feat. Focusing on the yeast Saccharomyces cerevisiae, we used complementary high-content screens to follow changes in localization of most proteins during growth in oleate. We found extensive changes in cellular architecture and identified several proteins that colocalized with peroxisomes that had not previously been considered peroxisomal proteins. One of the newly identified peroxisomal proteins, Ymr018w, is a protein with an unknown function that is similar to the yeast and human peroxisomal targeting receptor Pex5. We demonstrate that Ymr018w is a new peroxisomal-targeting receptor that targets a subset of matrix proteins to peroxisomes. We, therefore, renamed Ymr018w, Pex9, and suggest that Pex9 is a condition-specific targeting receptor that enables the dynamic rewiring of peroxisomes in response to metabolic needs. Moreover, we suggest that Pex5-like receptors might also exist in vertebrates.

Accurate, model-based tuning of synthetic gene expression using introns in S. cerevisiae.

Introns are key regulators of eukaryotic gene expression and present a potentially powerful tool for the design of synthetic eukaryotic gene expression systems. However, intronic control over gene expression is governed by a multitude of complex, incompletely understood, regulatory mechanisms. Despite this lack of detailed mechanistic understanding, here we show how a relatively simple model enables accurate and predictable tuning of synthetic gene expression system in yeast using several predictive intron features such as transcript folding and sequence motifs. Using only natural Saccharomyces cerevisiae introns as regulators, we demonstrate fine and accurate control over gene expression spanning a 100 fold expression range. These results broaden the engineering toolbox of synthetic gene expression systems and provide a framework in which precise and robust tuning of gene expression is accomplished.

Yeast phospholipid biosynthesis is linked to mRNA localization.

Regulation of the localization of mRNAs and local translation are universal features in eukaryotes and contribute to cellular asymmetry and differentiation. In Saccharomyces cerevisiae, localization of mRNAs that encode membrane proteins requires the She protein machinery, including the RNA-binding protein She2p, as well as movement of the cortical endoplasmic reticulum (cER) to the yeast bud. In a screen for ER-specific proteins necessary for the directional transport of WSC2 and EAR1 mRNAs, we have identified enzymes that are involved in phospholipid metabolism. Loss of the phospholipid methyltransferase Cho2p, which showed the strongest impact on mRNA localization, disturbs mRNA localization, as well as ER morphology and segregation, owing to an increase in the amount of cellular phosphatidylethanolamine (PtdEtn). Mislocalized mRNPs containing She2p colocalize with aggregated cER structures, suggestive of the entrapment of mRNA and She2p by the elevated PtdEtn level. This was confirmed by the elevated binding of She2p to PtdEtn-containing liposomes. These findings underscore the importance of ER membrane integrity in mRNA transport.

Spanning high-dimensional expression space using ribosome-binding site combinatorics.

Protein levels are a dominant factor shaping natural and synthetic biological systems. Although proper functioning of metabolic pathways relies on precise control of enzyme levels, the experimental ability to balance the levels of many genes in parallel is a major outstanding challenge. Here, we introduce a rapid and modular method to span the expression space of several proteins in parallel. By combinatorially pairing genes with a compact set of ribosome-binding sites, we modulate protein abundance by several orders of magnitude. We demonstrate our strategy by using a synthetic operon containing fluorescent proteins to span a 3D color space. Using the same approach, we modulate a recombinant carotenoid biosynthesis pathway in Escherichia coli to reveal a diversity of phenotypes, each characterized by a distinct carotenoid accumulation profile. In a single combinatorial assembly, we achieve a yield of the industrially valuable compound astaxanthin 4-fold higher than previously reported. The methodology presented here provides an efficient tool for exploring a high-dimensional expression space to locate desirable phenotypes.

Primers-4-Yeast: a comprehensive web tool for planning primers for Saccharomyces cerevisiae.

The budding yeast Saccharomyces cerevisiae is a key model organism of functional genomics, due to its ease and speed of genetic manipulations. In fact, in this yeast, the requirement for homologous sequences for recombination purposes is so small that 40 base pairs (bp) are sufficient. Hence, an enormous variety of genetic manipulations can be performed by simply planning primers with the correct homology, using a defined set of transformation plasmids. Although designing primers for yeast transformations and for the verification of their correct insertion is a common task in all yeast laboratories, primer planning is usually done manually and a tool that would enable easy, automated primer planning for the yeast research community is still lacking. Here we introduce Primers-4-Yeast, a web tool that allows primers to be designed in batches for S. cerevisiae gene-targeting transformations, and for the validation of correct insertions. This novel tool enables fast, automated, accurate primer planning for large sets of genes, introduces consistency in primer planning and is therefore suggested to serve as a standard in yeast research. Primers-4-Yeast is available at: http://www.weizmann.ac.il/Primers-4-Yeast

Modeling T cell temporal response to cancer immunotherapy rationalizes development of combinatorial treatment protocols.

Successful immunotherapy relies on triggering complex responses involving T cell dynamics in tumors and the periphery. Characterizing these responses remains challenging using static human single-cell atlases or mouse models. To address this, we developed a framework for in vivo tracking of tumor-specific CD8+ T cells over time and at single-cell resolution. Our tools facilitate the modeling of gene program dynamics in the tumor microenvironment (TME) and the tumor-draining lymph node (tdLN). Using this approach, we characterize two modes of anti-programmed cell death protein 1 (PD-1) activity, decoupling induced differentiation of tumor-specific activated precursor cells from conventional type 1 dendritic cell (cDC1)-dependent proliferation and recruitment to the TME. We demonstrate that combining anti-PD-1 therapy with anti-4-1BB agonist enhances the recruitment and proliferation of activated precursors, resulting in tumor control. These data suggest that effective response to anti-PD-1 therapy is dependent on sufficient influx of activated precursor CD8+ cells to the TME and highlight the importance of understanding system-level dynamics in optimizing immunotherapies.

Spatial and Temporal Mapping of Breast Cancer Lung Metastases Identify TREM2 Macrophages as Regulators of the Metastatic Boundary.

Cancer mortality primarily stems from metastatic recurrence, emphasizing the urgent need for developing effective metastasis-targeted immunotherapies. To better understand the cellular and molecular events shaping metastatic niches, we used a spontaneous breast cancer lung metastasis model to create a single-cell atlas spanning different metastatic stages and regions. We found that premetastatic lungs are infiltrated by inflammatory neutrophils and monocytes, followed by the accumulation of suppressive macrophages with the emergence of metastases. Spatial profiling revealed that metastasis-associated immune cells were present in the metastasis core, with the exception of TREM2+ regulatory macrophages uniquely enriched at the metastatic invasive margin, consistent across both murine models and human patient samples. These regulatory macrophages (Mreg) contribute to the formation of an immune-suppressive niche, cloaking tumor cells from immune surveillance. Our study provides a compendium of immune cell dynamics across metastatic stages and niches, informing the development of metastasis-targeting immunotherapies. SIGNIFICANCE: Temporal and spatial single-cell analysis of metastasis stages revealed new players in modulating immune surveillance and suppression. Our study highlights distinct populations of TREM2 macrophages as modulators of the microenvironment in metastasis, and as the key immune determinant defining metastatic niches, pointing to myeloid checkpoints to improve therapeutic strategies. This article is featured in Selected Articles from This Issue, p. 2489.

DestVI identifies continuums of cell types in spatial transcriptomics data.

Most spatial transcriptomics technologies are limited by their resolution, with spot sizes larger than that of a single cell. Although joint analysis with single-cell RNA sequencing can alleviate this problem, current methods are limited to assessing discrete cell types, revealing the proportion of cell types inside each spot. To identify continuous variation of the transcriptome within cells of the same type, we developed Deconvolution of Spatial Transcriptomics profiles using Variational Inference (DestVI). Using simulations, we demonstrate that DestVI outperforms existing methods for estimating gene expression for every cell type inside every spot. Applied to a study of infected lymph nodes and of a mouse tumor model, DestVI provides high-resolution, accurate spatial characterization of the cellular organization of these tissues and identifies cell-type-specific changes in gene expression between different tissue regions or between conditions. DestVI is available as part of the open-source software package scvi-tools ( https://scvi-tools.org ).

Anti-CTLA-4 antibodies drive myeloid activation and reprogram the tumor microenvironment through FcγR engagement and type I interferon signaling.

Despite the clinical success of checkpoint inhibitors, a substantial gap still exists in our understanding of their mechanism of action. While antibodies to cytotoxic T lymphocyte-associated protein-4 (CTLA-4) were developed to block inhibitory signals in T cells, several recent studies have demonstrated that Fcγ receptor (FcγR)-dependent depletion of regulatory T cells (Treg) is critical for antitumor activity. Here, using single-cell RNA sequencing, we dissect the impact of anti-CTLA-4-blocking, Treg cell-depleting and FcR-engaging activity on the immune response within tumors. We observed a rapid remodeling of the innate immune landscape as early as 24 h after treatment. Using genetic Treg cell ablation models, we show that immune remodeling was not driven solely by Treg cell depletion or CTLA-4 blockade but mainly through FcγR engagement, downstream activation of type I interferon signaling and reduction of suppressive macrophages. Our findings indicate that FcγR engagement and innate immune remodeling are involved in successful anti-CTLA-4 treatment, supporting the development of optimized immunotherapy agents bearing these features.