Protective role for kidney TREM2high macrophages in obesity- and diabetes-induced kidney injury.

Diabetic kidney disease (DKD), the most common cause of kidney failure, is a frequent complication of diabetes and obesity, and yet to date, treatments to halt its progression are lacking. We analyze kidney single-cell transcriptomic profiles from DKD patients and two DKD mouse models at multiple time points along disease progression-high-fat diet (HFD)-fed mice aged to 90-100 weeks and BTBR ob/ob mice (a genetic model)-and report an expanding population of macrophages with high expression of triggering receptor expressed on myeloid cells 2 (TREM2) in HFD-fed mice. TREM2high macrophages are enriched in obese and diabetic patients, in contrast to hypertensive patients or healthy controls in an independent validation cohort. Trem2 knockout mice on an HFD have worsening kidney filter damage and increased tubular epithelial cell injury, all signs of worsening DKD. Together, our studies suggest that strategies to enhance kidney TREM2high macrophages may provide therapeutic benefits for DKD.

Long interspersed nuclear elements safeguard neural progenitors from precocious differentiation.

Long interspersed nuclear element-1 (L1 or LINE-1) is a highly abundant mobile genetic element in both humans and mice, comprising almost 20% of each genome. L1s are silenced by several mechanisms, as their uncontrolled expression has the potential to induce genomic instability. However, L1s are paradoxically expressed at high levels in differentiating neural progenitor cells. Using in vitro and in vivo techniques to modulate L1 expression, we report that L1s play a critical role in both human and mouse brain development by regulating the rate of neural differentiation in a reverse-transcription-independent manner.

Junior scientists spotlight social bonds in seminars for diversity, equity, and inclusion in STEM.

Disparities for women and minorities in science, technology, engineering, and math (STEM) careers have continued even amidst mounting evidence for the superior performance of diverse workforces. In response, we launched the Diversity and Science Lecture series, a cross-institutional platform where junior life scientists present their research and comment on diversity, equity, and inclusion in STEM. We characterize speaker representation from 79 profiles and investigate topic noteworthiness via quantitative content analysis of talk transcripts. Nearly every speaker discussed interpersonal support, and three-fifths of speakers commented on race or ethnicity. Other topics, such as sexual and gender minority identity, were less frequently addressed but highly salient to the speakers who mentioned them. We found that significantly co-occurring topics reflected not only conceptual similarity, such as terms for racial identities, but also intersectional significance, such as identifying as a Latina/Hispanic woman or Asian immigrant, and interactions between concerns and identities, including the heightened value of friendship to the LGBTQ community, which we reproduce using transcripts from an independent seminar series. Our approach to scholar profiles and talk transcripts serves as an example for transmuting hundreds of hours of scholarly discourse into rich datasets that can power computational audits of speaker diversity and illuminate speakers' personal and professional priorities.

Cancer aneuploidies are shaped primarily by effects on tumour fitness.

Aneuploidies-whole-chromosome or whole-arm imbalances-are the most prevalent alteration in cancer genomes1,2. However, it is still debated whether their prevalence is due to selection or ease of generation as passenger events1,2. Here we developed a method, BISCUT, that identifies loci subject to fitness advantages or disadvantages by interrogating length distributions of telomere- or centromere-bounded copy-number events. These loci were significantly enriched for known cancer driver genes, including genes not detected through analysis of focal copy-number events, and were often lineage specific. BISCUT identified the helicase-encoding gene WRN as a haploinsufficient tumour-suppressor gene on chromosome 8p, which is supported by several lines of evidence. We also formally quantified the role of selection and mechanical biases in driving aneuploidy, finding that rates of arm-level copy-number alterations are most highly correlated with their effects on cellular fitness1,2. These results provide insight into the driving forces behind aneuploidy and its contribution to tumorigenesis.

Structural and functional properties of mSWI/SNF chromatin remodeling complexes revealed through single-cell perturbation screens.

The mammalian SWI/SNF (mSWI/SNF or BAF) family of chromatin remodeling complexes play critical roles in regulating DNA accessibility and gene expression. The three final-form subcomplexes-cBAF, PBAF, and ncBAF-are distinct in biochemical componentry, chromatin targeting, and roles in disease; however, the contributions of their constituent subunits to gene expression remain incompletely defined. Here, we performed Perturb-seq-based CRISPR-Cas9 knockout screens targeting mSWI/SNF subunits individually and in select combinations, followed by single-cell RNA-seq and SHARE-seq. We uncovered complex-, module-, and subunit-specific contributions to distinct regulatory networks and defined paralog subunit relationships and shifted subcomplex functions upon perturbations. Synergistic, intra-complex genetic interactions between subunits reveal functional redundancy and modularity. Importantly, single-cell subunit perturbation signatures mapped across bulk primary human tumor expression profiles both mirror and predict cBAF loss-of-function status in cancer. Our findings highlight the utility of Perturb-seq to dissect disease-relevant gene regulatory impacts of heterogeneous, multi-component master regulatory complexes.

Single-nucleus and spatial transcriptome profiling of pancreatic cancer identifies multicellular dynamics associated with neoadjuvant treatment.

Pancreatic ductal adenocarcinoma (PDAC) is a highly lethal and treatment-refractory cancer. Molecular stratification in pancreatic cancer remains rudimentary and does not yet inform clinical management or therapeutic development. Here, we construct a high-resolution molecular landscape of the cellular subtypes and spatial communities that compose PDAC using single-nucleus RNA sequencing and whole-transcriptome digital spatial profiling (DSP) of 43 primary PDAC tumor specimens that either received neoadjuvant therapy or were treatment naive. We uncovered recurrent expression programs across malignant cells and fibroblasts, including a newly identified neural-like progenitor malignant cell program that was enriched after chemotherapy and radiotherapy and associated with poor prognosis in independent cohorts. Integrating spatial and cellular profiles revealed three multicellular communities with distinct contributions from malignant, fibroblast and immune subtypes: classical, squamoid-basaloid and treatment enriched. Our refined molecular and cellular taxonomy can provide a framework for stratification in clinical trials and serve as a roadmap for therapeutic targeting of specific cellular phenotypes and multicellular interactions.

Single-nucleus cross-tissue molecular reference maps toward understanding disease gene function.

Understanding gene function and regulation in homeostasis and disease requires knowledge of the cellular and tissue contexts in which genes are expressed. Here, we applied four single-nucleus RNA sequencing methods to eight diverse, archived, frozen tissue types from 16 donors and 25 samples, generating a cross-tissue atlas of 209,126 nuclei profiles, which we integrated across tissues, donors, and laboratory methods with a conditional variational autoencoder. Using the resulting cross-tissue atlas, we highlight shared and tissue-specific features of tissue-resident cell populations; identify cell types that might contribute to neuromuscular, metabolic, and immune components of monogenic diseases and the biological processes involved in their pathology; and determine cell types and gene modules that might underlie disease mechanisms for complex traits analyzed by genome-wide association studies.

Dnmt3a knockout in excitatory neurons impairs postnatal synapse maturation and increases the repressive histone modification H3K27me3.

Two epigenetic pathways of transcriptional repression, DNA methylation and polycomb repressive complex 2 (PRC2), are known to regulate neuronal development and function. However, their respective contributions to brain maturation are unknown. We found that conditional loss of the de novo DNA methyltransferase Dnmt3a in mouse excitatory neurons altered expression of synapse-related genes, stunted synapse maturation, and impaired working memory and social interest. At the genomic level, loss of Dnmt3a abolished postnatal accumulation of CG and non-CG DNA methylation, leaving adult neurons with an unmethylated, fetal-like epigenomic pattern at ~222,000 genomic regions. The PRC2-associated histone modification, H3K27me3, increased at many of these sites. Our data support a dynamic interaction between two fundamental modes of epigenetic repression during postnatal maturation of excitatory neurons, which together confer robustness on neuronal regulation.

Conventional type I dendritic cells maintain a reservoir of proliferative tumor-antigen specific TCF-1+ CD8+ T cells in tumor-draining lymph nodes.

In tumors, a subset of CD8+ T cells expressing the transcription factor TCF-1 drives the response to immune checkpoint blockade. We examined the mechanisms that maintain these cells in an autochthonous model of lung adenocarcinoma. Longitudinal sampling and single-cell sequencing of tumor-antigen specific TCF-1+ CD8+ T cells revealed that while intratumoral TCF-1+ CD8+ T cells acquired dysfunctional features and decreased in number as tumors progressed, TCF-1+ CD8+ T cell frequency in the tumor draining LN (dLN) remained stable. Two discrete intratumoral TCF-1+ CD8+ T cell subsets developed over time-a proliferative SlamF6+ subset and a non-cycling SlamF6- subset. Blocking dLN egress decreased the frequency of intratumoral SlamF6+ TCF-1+ CD8+ T cells. Conventional type I dendritic cell (cDC1) in dLN decreased in number with tumor progression, and Flt3L+anti-CD40 treatment recovered SlamF6+ T cell frequencies and decreased tumor burden. Thus, cDC1s in tumor dLN maintain a reservoir of TCF-1+ CD8+ T cells and their decrease contributes to failed anti-tumor immunity.

Cycling cancer persister cells arise from lineages with distinct programs.

Non-genetic mechanisms have recently emerged as important drivers of cancer therapy failure1, where some cancer cells can enter a reversible drug-tolerant persister state in response to treatment2. Although most cancer persisters remain arrested in the presence of the drug, a rare subset can re-enter the cell cycle under constitutive drug treatment. Little is known about the non-genetic mechanisms that enable cancer persisters to maintain proliferative capacity in the presence of drugs. To study this rare, transiently resistant, proliferative persister population, we developed Watermelon, a high-complexity expressed barcode lentiviral library for simultaneous tracing of each cell's clonal origin and proliferative and transcriptional states. Here we show that cycling and non-cycling persisters arise from different cell lineages with distinct transcriptional and metabolic programs. Upregulation of antioxidant gene programs and a metabolic shift to fatty acid oxidation are associated with persister proliferative capacity across multiple cancer types. Impeding oxidative stress or metabolic reprogramming alters the fraction of cycling persisters. In human tumours, programs associated with cycling persisters are induced in minimal residual disease in response to multiple targeted therapies. The Watermelon system enabled the identification of rare persister lineages that are preferentially poised to proliferate under drug pressure, thus exposing new vulnerabilities that can be targeted to delay or even prevent disease recurrence.

Spatially organized multicellular immune hubs in human colorectal cancer.

Immune responses to cancer are highly variable, with mismatch repair-deficient (MMRd) tumors exhibiting more anti-tumor immunity than mismatch repair-proficient (MMRp) tumors. To understand the rules governing these varied responses, we transcriptionally profiled 371,223 cells from colorectal tumors and adjacent normal tissues of 28 MMRp and 34 MMRd individuals. Analysis of 88 cell subsets and their 204 associated gene expression programs revealed extensive transcriptional and spatial remodeling across tumors. To discover hubs of interacting malignant and immune cells, we identified expression programs in different cell types that co-varied across tumors from affected individuals and used spatial profiling to localize coordinated programs. We discovered a myeloid cell-attracting hub at the tumor-luminal interface associated with tissue damage and an MMRd-enriched immune hub within the tumor, with activated T cells together with malignant and myeloid cells expressing T cell-attracting chemokines. By identifying interacting cellular programs, we reveal the logic underlying spatially organized immune-malignant cell networks.

Tumor and immune reprogramming during immunotherapy in advanced renal cell carcinoma.

Immune checkpoint blockade (ICB) results in durable disease control in a subset of patients with advanced renal cell carcinoma (RCC), but mechanisms driving resistance are poorly understood. We characterize the single-cell transcriptomes of cancer and immune cells from metastatic RCC patients before or after ICB exposure. In responders, subsets of cytotoxic T cells express higher levels of co-inhibitory receptors and effector molecules. Macrophages from treated biopsies shift toward pro-inflammatory states in response to an interferon-rich microenvironment but also upregulate immunosuppressive markers. In cancer cells, we identify bifurcation into two subpopulations differing in angiogenic signaling and upregulation of immunosuppressive programs after ICB. Expression signatures for cancer cell subpopulations and immune evasion are associated with PBRM1 mutation and survival in primary and ICB-treated advanced RCC. Our findings demonstrate that ICB remodels the RCC microenvironment and modifies the interplay between cancer and immune cell populations critical for understanding response and resistance to ICB.

Single-cell meta-analysis of SARS-CoV-2 entry genes across tissues and demographics.

Angiotensin-converting enzyme 2 (ACE2) and accessory proteases (TMPRSS2 and CTSL) are needed for severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) cellular entry, and their expression may shed light on viral tropism and impact across the body. We assessed the cell-type-specific expression of ACE2, TMPRSS2 and CTSL across 107 single-cell RNA-sequencing studies from different tissues. ACE2, TMPRSS2 and CTSL are coexpressed in specific subsets of respiratory epithelial cells in the nasal passages, airways and alveoli, and in cells from other organs associated with coronavirus disease 2019 (COVID-19) transmission or pathology. We performed a meta-analysis of 31 lung single-cell RNA-sequencing studies with 1,320,896 cells from 377 nasal, airway and lung parenchyma samples from 228 individuals. This revealed cell-type-specific associations of age, sex and smoking with expression levels of ACE2, TMPRSS2 and CTSL. Expression of entry factors increased with age and in males, including in airway secretory cells and alveolar type 2 cells. Expression programs shared by ACE2+TMPRSS2+ cells in nasal, lung and gut tissues included genes that may mediate viral entry, key immune functions and epithelial-macrophage cross-talk, such as genes involved in the interleukin-6, interleukin-1, tumor necrosis factor and complement pathways. Cell-type-specific expression patterns may contribute to the pathogenesis of COVID-19, and our work highlights putative molecular pathways for therapeutic intervention.

Blockade of IL-22 signaling reverses erythroid dysfunction in stress-induced anemias.

Patients with myelodysplastic syndromes (MDSs) display severe anemia but the mechanisms underlying this phenotype are incompletely understood. Right open-reading-frame kinase 2 (RIOK2) encodes a protein kinase located at 5q15, a region frequently lost in patients with MDS del(5q). Here we show that hematopoietic cell-specific haploinsufficient deletion of Riok2 (Riok2f/+Vav1cre) led to reduced erythroid precursor frequency leading to anemia. Proteomic analysis of Riok2f/+Vav1cre erythroid precursors suggested immune system activation, and transcriptomic analysis revealed an increase in p53-dependent interleukin (IL)-22 in Riok2f/+Vav1cre CD4+ T cells (TH22). Further, we discovered that the IL-22 receptor, IL-22RA1, was unexpectedly present on erythroid precursors. Blockade of IL-22 signaling alleviated anemia not only in Riok2f/+Vav1cre mice but also in wild-type mice. Serum concentrations of IL-22 were increased in the subset of patients with del(5q) MDS as well as patients with anemia secondary to chronic kidney disease. This work reveals a possible therapeutic opportunity for reversing many stress-induced anemias by targeting IL-22 signaling.

Transcriptional mediators of treatment resistance in lethal prostate cancer.

Metastatic castration-resistant prostate cancer is typically lethal, exhibiting intrinsic or acquired resistance to second-generation androgen-targeting therapies and minimal response to immune checkpoint inhibitors1. Cellular programs driving resistance in both cancer and immune cells remain poorly understood. We present single-cell transcriptomes from 14 patients with advanced prostate cancer, spanning all common metastatic sites. Irrespective of treatment exposure, adenocarcinoma cells pervasively coexpressed multiple androgen receptor isoforms, including truncated isoforms hypothesized to mediate resistance to androgen-targeting therapies2,3. Resistance to enzalutamide was associated with cancer cell-intrinsic epithelial-mesenchymal transition and transforming growth factor-ß signaling. Small cell carcinoma cells exhibited divergent expression programs driven by transcriptional regulators promoting lineage plasticity and HOXB5, HOXB6 and NR1D2 (refs. 4-6). Additionally, a subset of patients had high expression of dysfunction markers on cytotoxic CD8+ T cells undergoing clonal expansion following enzalutamide treatment. Collectively, the transcriptional characterization of cancer and immune cells from human metastatic castration-resistant prostate cancer provides a basis for the development of therapeutic approaches complementing androgen signaling inhibition.

Multimodal pooled Perturb-CITE-seq screens in patient models define mechanisms of cancer immune evasion.

Resistance to immune checkpoint inhibitors (ICIs) is a key challenge in cancer therapy. To elucidate underlying mechanisms, we developed Perturb-CITE-sequencing (Perturb-CITE-seq), enabling pooled clustered regularly interspaced short palindromic repeat (CRISPR)-Cas9 perturbations with single-cell transcriptome and protein readouts. In patient-derived melanoma cells and autologous tumor-infiltrating lymphocyte (TIL) co-cultures, we profiled transcriptomes and 20 proteins in ~218,000 cells under ~750 perturbations associated with cancer cell-intrinsic ICI resistance (ICR). We recover known mechanisms of resistance, including defects in the interferon-γ (IFN-γ)-JAK/STAT and antigen-presentation pathways in RNA, protein and perturbation space, and new ones, including loss/downregulation of CD58. Loss of CD58 conferred immune evasion in multiple co-culture models and was downregulated in tumors of melanoma patients with ICR. CD58 protein expression was not induced by IFN-γ signaling, and CD58 loss conferred immune evasion without compromising major histocompatibility complex (MHC) expression, suggesting that it acts orthogonally to known mechanisms of ICR. This work provides a framework for the deciphering of complex mechanisms by large-scale perturbation screens with multimodal, single-cell readouts, and discovers potentially clinically relevant mechanisms of immune evasion.

Opposing immune and genetic mechanisms shape oncogenic programs in synovial sarcoma.

Synovial sarcoma (SyS) is an aggressive neoplasm driven by the SS18-SSX fusion, and is characterized by low T cell infiltration. Here, we studied the cancer-immune interplay in SyS using an integrative approach that combines single-cell RNA sequencing (scRNA-seq), spatial profiling and genetic and pharmacological perturbations. scRNA-seq of 16,872 cells from 12 human SyS tumors uncovered a malignant subpopulation that marks immune-deprived niches in situ and is predictive of poor clinical outcomes in two independent cohorts. Functional analyses revealed that this malignant cell state is controlled by the SS18-SSX fusion, is repressed by cytokines secreted by macrophages and T cells, and can be synergistically targeted with a combination of HDAC and CDK4/CDK6 inhibitors. This drug combination enhanced malignant-cell immunogenicity in SyS models, leading to induced T cell reactivity and T cell-mediated killing. Our study provides a blueprint for investigating heterogeneity in fusion-driven malignancies and demonstrates an interplay between immune evasion and oncogenic processes that can be co-targeted in SyS and potentially in other malignancies.

Pan-cancer single-cell RNA-seq identifies recurring programs of cellular heterogeneity.

Cultured cell lines are the workhorse of cancer research, but the extent to which they recapitulate the heterogeneity observed among malignant cells in tumors is unclear. Here we used multiplexed single-cell RNA-seq to profile 198 cancer cell lines from 22 cancer types. We identified 12 expression programs that are recurrently heterogeneous within multiple cancer cell lines. These programs are associated with diverse biological processes, including cell cycle, senescence, stress and interferon responses, epithelial-mesenchymal transition and protein metabolism. Most of these programs recapitulate those recently identified as heterogeneous within human tumors. We prioritized specific cell lines as models of cellular heterogeneity and used them to study subpopulations of senescence-related cells, demonstrating their dynamics, regulation and unique drug sensitivities, which were predictive of clinical response. Our work describes the landscape of heterogeneity within diverse cancer cell lines and identifies recurrent patterns of heterogeneity that are shared between tumors and specific cell lines.

The Human and Mouse Enteric Nervous System at Single-Cell Resolution.

The enteric nervous system (ENS) coordinates diverse functions in the intestine but has eluded comprehensive molecular characterization because of the rarity and diversity of cells. Here we develop two methods to profile the ENS of adult mice and humans at single-cell resolution: RAISIN RNA-seq for profiling intact nuclei with ribosome-bound mRNA and MIRACL-seq for label-free enrichment of rare cell types by droplet-based profiling. The 1,187,535 nuclei in our mouse atlas include 5,068 neurons from the ileum and colon, revealing extraordinary neuron diversity. We highlight circadian expression changes in enteric neurons, show that disease-related genes are dysregulated with aging, and identify differences between the ileum and proximal/distal colon. In humans, we profile 436,202 nuclei, recovering 1,445 neurons, and identify conserved and species-specific transcriptional programs and putative neuro-epithelial, neuro-stromal, and neuro-immune interactions. The human ENS expresses risk genes for neuropathic, inflammatory, and extra-intestinal diseases, suggesting neuronal contributions to disease.