The CD8+ T cell tolerance checkpoint triggers a distinct differentiation state defined by protein translation defects.

Peripheral CD8+ T cell tolerance is a checkpoint in both autoimmune disease and anti-cancer immunity. Despite its importance, the relationship between tolerance-induced states and other CD8+ T cell differentiation states remains unclear. Using flow cytometric phenotyping, single-cell RNA sequencing (scRNA-seq), and chromatin accessibility profiling, we demonstrated that in vivo peripheral tolerance to a self-antigen triggered a fundamentally distinct differentiation state separate from exhaustion, memory, and functional effector cells but analogous to cells defectively primed against tumors. Tolerant cells diverged early and progressively from effector cells, adopting a transcriptionally and epigenetically distinct state within 60 h of antigen encounter. Breaching tolerance required the synergistic actions of strong T cell receptor (TCR) signaling and inflammation, which cooperatively induced gene modules that enhanced protein translation. Weak TCR signaling during bystander infection failed to breach tolerance due to the uncoupling of effector gene expression from protein translation. Thus, tolerance engages a distinct differentiation trajectory enforced by protein translation defects.

Selfish conflict underlies RNA-mediated parent-of-origin effects.

Genomic imprinting-the non-equivalence of maternal and paternal genomes-is a critical process that has evolved independently in many plant and mammalian species1,2. According to kinship theory, imprinting is the inevitable consequence of conflictive selective forces acting on differentially expressed parental alleles3,4. Yet, how these epigenetic differences evolve in the first place is poorly understood3,5,6. Here we report the identification and molecular dissection of a parent-of-origin effect on gene expression that might help to clarify this fundamental question. Toxin-antidote elements (TAs) are selfish elements that spread in populations by poisoning non-carrier individuals7-9. In reciprocal crosses between two Caenorhabditis tropicalis wild isolates, we found that the slow-1/grow-1 TA is specifically inactive when paternally inherited. This parent-of-origin effect stems from transcriptional repression of the slow-1 toxin by the PIWI-interacting RNA (piRNA) host defence pathway. The repression requires PIWI Argonaute and SET-32 histone methyltransferase activities and is transgenerationally inherited via small RNAs. Remarkably, when slow-1/grow-1 is maternally inherited, slow-1 repression is halted by a translation-independent role of its maternal mRNA. That is, slow-1 transcripts loaded into eggs-but not SLOW-1 protein-are necessary and sufficient to counteract piRNA-mediated repression. Our findings show that parent-of-origin effects can evolve by co-option of the piRNA pathway and hinder the spread of selfish genes that require sex for their propagation.

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.

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.

Cholesterol 24-hydroxylase at the choroid plexus contributes to brain immune homeostasis.

The choroid plexus (CP) plays a key role in remotely controlling brain function in health, aging, and disease. Here, we report that CP epithelial cells express the brain-specific cholesterol 24-hydroxylase (CYP46A1) and that its levels are decreased under different mouse and human brain conditions, including amyloidosis, aging, and SARS-CoV-2 infection. Using primary mouse CP cell cultures, we demonstrate that the enzymatic product of CYP46A1, 24(S)-hydroxycholesterol, downregulates inflammatory transcriptomic signatures within the CP, found here to be elevated across multiple neurological conditions. In vitro, the pro-inflammatory cytokine tumor necrosis factor α (TNF-α) downregulates CYP46A1 expression, while overexpression of CYP46A1 or its pharmacological activation in mouse CP organ cultures increases resilience to TNF-α. In vivo, overexpression of CYP46A1 in the CP in transgenic mice with amyloidosis is associated with better cognitive performance and decreased brain inflammation. Our findings suggest that CYP46A1 expression in the CP impacts the role of this niche as a guardian of brain immune homeostasis.

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.

Pumping iron: A multi-omics analysis of two extremophilic algae reveals iron economy management.

Marine algae are responsible for half of the world's primary productivity, but this critical carbon sink is often constrained by insufficient iron. One species of marine algae, Dunaliella tertiolecta, is remarkable for its ability to maintain photosynthesis and thrive in low-iron environments. A related species, Dunaliella salina Bardawil, shares this attribute but is an extremophile found in hypersaline environments. To elucidate how algae manage their iron requirements, we produced high-quality genome assemblies and transcriptomes for both species to serve as a foundation for a comparative multiomics analysis. We identified a host of iron-uptake proteins in both species, including a massive expansion of transferrins and a unique family of siderophore-iron-uptake proteins. Complementing these multiple iron-uptake routes, ferredoxin functions as a large iron reservoir that can be released by induction of flavodoxin. Proteomic analysis revealed reduced investment in the photosynthetic apparatus coupled with remodeling of antenna proteins by dramatic iron-deficiency induction of TIDI1, which is closely related but identifiably distinct from the chlorophyll binding protein, LHCA3. These combinatorial iron scavenging and sparing strategies make Dunaliella unique among photosynthetic organisms.

The transcriptional and regulatory identity of erythropoietin producing cells.

Erythropoietin (Epo) is the master regulator of erythropoiesis and oxygen homeostasis. Despite its physiological importance, the molecular and genomic contexts of the cells responsible for renal Epo production remain unclear, limiting more-effective therapies for anemia. Here, we performed single-cell RNA and transposase-accessible chromatin (ATAC) sequencing of an Epo reporter mouse to molecularly identify Epo-producing cells under hypoxic conditions. Our data indicate that a distinct population of kidney stroma, which we term Norn cells, is the major source of endocrine Epo production in mice. We use these datasets to identify the markers, signaling pathways and transcriptional circuits characteristic of Norn cells. Using single-cell RNA sequencing and RNA in situ hybridization in human kidney tissues, we further provide evidence that this cell population is conserved in humans. These preliminary findings open new avenues to functionally dissect EPO gene regulation in health and disease and may serve as groundwork to improve erythropoiesis-stimulating therapies.

A single-cell census of mouse limb development identifies complex spatiotemporal dynamics of skeleton formation.

Limb development has long served as a model system for coordinated spatial patterning of progenitor cells. Here, we identify a population of naive limb progenitors and show that they differentiate progressively to form the skeleton in a complex, non-consecutive, three-dimensional pattern. Single-cell RNA sequencing of the developing mouse forelimb identified three progenitor states: naive, proximal, and autopodial, as well as Msx1 as a marker for the naive progenitors. In vivo lineage tracing confirmed this role and localized the naive progenitors to the outer margin of the limb, along the anterior-posterior axis. Sequential pulse-chase experiments showed that the progressive transition of Msx1+ naive progenitors into proximal and autopodial progenitors coincides with their differentiation to Sox9+ chondroprogenitors, which occurs along all the forming skeletal segments. Indeed, tracking the spatiotemporal sequence of differentiation showed that the skeleton forms progressively in a complex pattern. These findings suggest an alternative model for limb skeleton development.

Early Infiltration of Innate Immune Cells to the Liver Depletes HNF4α and Promotes Extrahepatic Carcinogenesis.

Multiple studies have identified metabolic changes within the tumor and its microenvironment during carcinogenesis. Yet, the mechanisms by which tumors affect the host metabolism are unclear. We find that systemic inflammation induced by cancer leads to liver infiltration of myeloid cells at early extrahepatic carcinogenesis. The infiltrating immune cells via IL6-pSTAT3 immune-hepatocyte cross-talk cause the depletion of a master metabolic regulator, HNF4α, consequently leading to systemic metabolic changes that promote breast and pancreatic cancer proliferation and a worse outcome. Preserving HNF4α levels maintains liver metabolism and restricts carcinogenesis. Standard liver biochemical tests can identify early metabolic changes and predict patients' outcomes and weight loss. Thus, the tumor induces early metabolic changes in its macroenvironment with diagnostic and potentially therapeutic implications for the host. SIGNIFICANCE: Cancer growth requires a permanent nutrient supply starting from early disease stages. We find that the tumor extends its effect to the host's liver to obtain nutrients and rewires the systemic and tissue-specific metabolism early during carcinogenesis. Preserving liver metabolism restricts tumor growth and improves cancer outcomes. This article is highlighted in the In This Issue feature, p. 1501.

Actin-dependent astrocytic infiltration is a key step for axon defasciculation during remodeling.

Astrocytes are essential for synapse formation, maturation, and plasticity; however, their function during developmental neuronal remodeling is largely unknown. To identify astrocytic molecules required for axon pruning of mushroom body (MB) γ neurons in Drosophila, we profiled astrocytes before (larva) and after (adult) remodeling. Focusing on genes enriched in larval astrocytes, we identified 12 astrocytic genes that are required for axon pruning, including the F-actin regulators Actin-related protein 2/3 complex, subunit 1 (Arpc1) and formin3 (form3). Interestingly, perturbing astrocytic actin dynamics does not affect their gross morphology, migration, or transforming growth factor ß (TGF-ß) secretion. In contrast, actin dynamics is required for astrocyte infiltration into the axon bundle at the onset of pruning. Remarkably, decreasing axonal adhesion facilitates infiltration by Arpc1 knockdown (KD) astrocytes and promotes axon pruning. Conversely, increased axonal adhesion reduces lobe infiltration by wild-type (WT) astrocytes. Together, our findings suggest that actin-dependent astrocytic infiltration is a key step in axon pruning, thus promoting our understanding of neuron-glia interactions during remodeling.

Neighbor-specific gene expression revealed from physically interacting cells during mouse embryonic development.

Development of multicellular organisms is orchestrated by persistent cell-cell communication between neighboring partners. Direct interaction between different cell types can induce molecular signals that dictate lineage specification and cell fate decisions. Current single-cell RNA-seq technology cannot adequately analyze cell-cell contact-dependent gene expression, mainly due to the loss of spatial information. To overcome this obstacle and resolve cell-cell contact-specific gene expression during embryogenesis, we performed RNA sequencing of physically interacting cells (PIC-seq) and assessed them alongside similar single-cell transcriptomes derived from developing mouse embryos between embryonic day (E) 7.5 and E9.5. Analysis of the PIC-seq data identified gene expression signatures that were dependent on the presence of specific neighboring cell types. Our computational predictions, validated experimentally, demonstrated that neural progenitor (NP) cells upregulate Lhx5 and Nkx2-1 genes, when exclusively interacting with definitive endoderm (DE) cells. Moreover, there was a reciprocal impact on the transcriptome of DE cells, as they tend to upregulate Rax and Gsc when in contact with NP cells. Using individual cell transcriptome data, we formulated a means of computationally predicting the impact of one cell type on the transcriptome of its neighboring cell types. We have further developed a distinctive spatial-t-distributed stochastic neighboring embedding to display the pseudospatial distribution of cells in a 2-dimensional space. In summary, we describe an innovative approach to study contact-specific gene regulation during embryogenesis.

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.

Island-specific evolution of a sex-primed autosome in a sexual planarian.

The sexual strain of the planarian Schmidtea mediterranea, indigenous to Tunisia and several Mediterranean islands, is a hermaphrodite1,2. Here we isolate individual chromosomes and use sequencing, Hi-C3,4 and linkage mapping to assemble a chromosome-scale genome reference. The linkage map reveals an extremely low rate of recombination on chromosome 1. We confirm suppression of recombination on chromosome 1 by genotyping individual sperm cells and oocytes. We show that previously identified genomic regions that maintain heterozygosity even after prolonged inbreeding make up essentially all of chromosome 1. Genome sequencing of individuals isolated in the wild indicates that this phenomenon has evolved specifically in populations from Sardinia and Corsica. We find that most known master regulators5-13 of the reproductive system are located on chromosome 1. We used RNA interference14,15 to knock down a gene with haplotype-biased expression, which led to the formation of a more pronounced female mating organ. On the basis of these observations, we propose that chromosome 1 is a sex-primed autosome primed for evolution into a sex chromosome.

Identification of the central intermediate in the extra-embryonic to embryonic endoderm transition through single-cell transcriptomics.

High-resolution maps of embryonic development suggest that acquisition of cell identity is not limited to canonical germ layers but proceeds via alternative routes. Despite evidence that visceral organs are formed via embryonic and extra-embryonic trajectories, the production of organ-specific cell types in vitro focuses on the embryonic one. Here we resolve these differentiation routes using massively parallel single-cell RNA sequencing to generate datasets from FOXA2Venus reporter mouse embryos and embryonic stem cell differentiation towards endoderm. To relate cell types in these datasets, we develop a single-parameter computational approach and identify an intermediate en route from extra-embryonic identity to embryonic endoderm, which we localize spatially in embryos at embryonic day 7.5. While there is little evidence for this cell type in embryonic stem cell differentiation, by following the extra-embryonic trajectory starting with naïve extra-embryonic endoderm stem cells we can generate embryonic gut spheroids. Exploiting developmental plasticity therefore offers alternatives to pluripotent cells and opens alternative avenues for in vitro differentiation.

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 ).

LGR5 expressing skin fibroblasts define a major cellular hub perturbed in scleroderma.

Systemic sclerosis (scleroderma, SSc) is an incurable autoimmune disease with high morbidity and mortality rates. Here, we conducted a population-scale single-cell genomic analysis of skin and blood samples of 56 healthy controls and 97 SSc patients at different stages of the disease. We found immune compartment dysfunction only in a specific subtype of diffuse SSc patients but global dysregulation of the stromal compartment, particularly in a previously undefined subset of LGR5+-scleroderma-associated fibroblasts (ScAFs). ScAFs are perturbed morphologically and molecularly in SSc patients. Single-cell multiome profiling of stromal cells revealed ScAF-specific markers, pathways, regulatory elements, and transcription factors underlining disease development. Systematic analysis of these molecular features with clinical metadata associates specific ScAF targets with disease pathogenesis and SSc clinical traits. Our high-resolution atlas of the sclerodermatous skin spectrum will enable a paradigm shift in the understanding of SSc disease and facilitate the development of biomarkers and therapeutic strategies.

Validation of the society of thoracic surgeons predicted risk of mortality score for long-term survival after cardiac surgery in Israel.

BACKGROUND: Long-term survival is an important metric in assessing procedural value. We previously confirmed that the Society of Thoracic Surgeons predicted risk of mortality score (PROM) accurately predicts 30-day mortality in Israeli patients. The present study investigated the ability of the PROM to reliably predict long-term survival. METHODS: Data on 1279 patients undergoing cardiac surgery were prospectively entered into our database and used to calculate PROM. Long-term mortality was obtained from the Israeli Social Security Database. Patients were stratified into five cohorts according to PROM (A: 0-0.99%, B: 1.0-1.99%, C: 2.0-2.99%, D: 3.0-4.99% and E: ≥ 5.0%). Kaplan-Meier estimates of survival were calculated for each cohort and compared by Wilcoxon signed-rank test. We used C-statistics to assess model discrimination. Cox regression analysis was performed to identify predictors of long-term survival. RESULTS: Follow-up was achieved for 1256 (98%) patients over a mean period of 62 ± 28 months (median 64, range 0-107). Mean survival of the entire cohort was 95 ± 1 (95% CI 93-96) months. Higher PROM was associated with reduced survival: A-104 ± 1 (103-105) months, B-96 ± 2 (93-99) months, C-93 ± 3 (88-98) months, D-89 ± 3 (84-94) months, E-74 ± 3 (68-80) months (p < 0.0001). The Area Under the Curve was 0.76 ± 0.02 indicating excellent model discrimination. Independent predictors of long-term mortality included advanced age, lower ejection fraction, reoperation, diabetes mellitus, dialysis and PROM. CONCLUSIONS: The PROM was a reliable predictor of long-term survival in Israeli patients undergoing cardiac surgery. The PROM might be a useful metric for assessing procedural value and surgical decision-making.