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.
Display omitted
•RAISIN RNA-seq of 1.6 million cells and 6,513 enteric neurons from adult mice and humans•Dozens of neuron subsets that vary by location, circadian phase, age, and species•Putative interactions between neurons and epithelial, stromal, and immune cells•Enteric neurons express risk genes for extraintestinal diseases with gut dysmotility
Comprehensive mapping and comparison of the adult mouse and human enteric nervous system is achieved with development of methods to isolate intact nuclei with ribosome-bound mRNA and label-free profiling of rare cell types. The atlases reveal the neuronal diversity of the ENS; its potential to link the gut, immune system, and brain; and contributions to disease.
Full text
Available for:
GEOZS, IJS, IMTLJ, KILJ, KISLJ, NLZOH, NUK, OILJ, PNG, SAZU, SBCE, SBJE, UILJ, UL, UM, UPCLJ, UPUK, ZAGLJ, ZRSKP
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.
Seventy percent of breast cancers express the estrogen receptor (ER), and agents that target the ER are the mainstay of treatment. However, virtually all people with ER
breast cancer develop ...resistance to ER-directed agents in the metastatic setting. Beyond mutations in the ER itself, which occur in 25-30% of people treated with aromatase inhibitors
, knowledge about clinical resistance mechanisms remains incomplete. We identified activating HER2 mutations in metastatic biopsies from eight patients with ER
metastatic breast cancer who had developed resistance to aromatase inhibitors, tamoxifen or fulvestrant. Examination of treatment-naive primary tumors in five patients showed no evidence of pre-existing mutations in four of five patients, suggesting that these mutations were acquired under the selective pressure of ER-directed therapy. The HER2 mutations and ER mutations were mutually exclusive, suggesting a distinct mechanism of acquired resistance to ER-directed therapies. In vitro analysis confirmed that the HER2 mutations conferred estrogen independence as well as-in contrast to ER mutations-resistance to tamoxifen, fulvestrant and the CDK4 and CDK6 inhibitor palbociclib. Resistance was overcome by combining ER-directed therapy with the irreversible HER2 kinase inhibitor neratinib.
Full text
Available for:
EMUNI, FIS, FZAB, GEOZS, GIS, IJS, IMTLJ, KILJ, KISLJ, MFDPS, NLZOH, NUK, OILJ, PNG, SAZU, SBCE, SBJE, SBMB, SBNM, UKNU, UL, UM, UPUK, VKSCE, ZAGLJ
Non-genetic mechanisms have recently emerged as important drivers of cancer therapy failure
, where some cancer cells can enter a reversible drug-tolerant persister state in response to treatment
. ...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.
Full text
Available for:
GEOZS, IJS, IMTLJ, KISLJ, NLZOH, NUK, OILJ, PNG, SAZU, SBCE, SBMB, UL, UM, UPUK, ZAGLJ
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.
Display omitted
•Distinct CD8+ T cell phenotypes are enriched after immune checkpoint blockade (ICB)•Immune checkpoint ligands are upregulated in macrophages and tumor cells after ICB•Two cancer cell subpopulations are conserved across heterogeneous RCC tumors•Cancer cell programs drive distinct immune interactions and predict patient outcomes
Bi et al. dissect the cancer cell and immune microenvironmental transcriptional programs in immune checkpoint blockade-exposed metastatic renal cell carcinoma, revealing immune subpopulation reprogramming and interactions with distinct cancer cell populations in the context of clinical resistance.
Full text
Available for:
GEOZS, IJS, IMTLJ, KILJ, KISLJ, NLZOH, NUK, OILJ, PNG, SAZU, SBCE, SBJE, UILJ, UL, UM, UPCLJ, UPUK, ZAGLJ, ZRSKP
Immune checkpoint inhibitors (ICIs) produce durable responses in some melanoma patients, but many patients derive no clinical benefit, and the molecular underpinnings of such resistance ...remain elusive. Here, we leveraged single-cell RNA sequencing (scRNA-seq) from 33 melanoma tumors and computational analyses to interrogate malignant cell states that promote immune evasion. We identified a resistance program expressed by malignant cells that is associated with T cell exclusion and immune evasion. The program is expressed prior to immunotherapy, characterizes cold niches in situ, and predicts clinical responses to anti-PD-1 therapy in an independent cohort of 112 melanoma patients. CDK4/6-inhibition represses this program in individual malignant cells, induces senescence, and reduces melanoma tumor outgrowth in mouse models in vivo when given in combination with immunotherapy. Our study provides a high-resolution landscape of ICI-resistant cell states, identifies clinically predictive signatures, and suggests new therapeutic strategies to overcome immunotherapy resistance.
Display omitted
•Single-cell RNA-seq identifies an immune resistance program in malignant cells•Multiple immune resistance mechanisms are co-regulated in the program•The program predicts clinical responses to immunotherapy in melanoma patients•CDK4/6 inhibitors repress the program and may sensitize melanoma to immunotherapy
Single-cell sequencing of checkpoint-inhibitor-resistant melanomas identifies predictive signatures to guide therapeutic approaches to overcome immunotherapy resistance.
Full text
Available for:
GEOZS, IJS, IMTLJ, KILJ, KISLJ, NLZOH, NUK, OILJ, PNG, SAZU, SBCE, SBJE, UILJ, UL, UM, UPCLJ, UPUK, ZAGLJ, ZRSKP
Aneuploidies-whole-chromosome or whole-arm imbalances-are the most prevalent alteration in cancer genomes
. However, it is still debated whether their prevalence is due to selection or ease of ...generation as passenger events
. 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 fitness
. These results provide insight into the driving forces behind aneuploidy and its contribution to tumorigenesis.
Full text
Available for:
IJS, KISLJ, NUK, SBMB, UL, UM, UPUK
Type 2 innate lymphoid cells (ILC2s) both contribute to mucosal homeostasis and initiate pathologic inflammation in allergic asthma. However, the signals that direct ILC2s to promote homeostasis ...versus inflammation are unclear. To identify such molecular cues, we profiled mouse lung-resident ILCs using single-cell RNA sequencing at steady state and after in vivo stimulation with the alarmin cytokines IL-25 and IL-33. ILC2s were transcriptionally heterogeneous after activation, with subpopulations distinguished by expression of proliferative, homeostatic and effector genes. The neuropeptide receptor Nmur1 was preferentially expressed by ILC2s at steady state and after IL-25 stimulation. Neuromedin U (NMU), the ligand of NMUR1, activated ILC2s in vitro, and in vivo co-administration of NMU with IL-25 strongly amplified allergic inflammation. Loss of NMU-NMUR1 signalling reduced ILC2 frequency and effector function, and altered transcriptional programs following allergen challenge in vivo. Thus, NMUR1 signalling promotes inflammatory ILC2 responses, highlighting the importance of neuro-immune crosstalk in allergic inflammation at mucosal surfaces.
Full text
Available for:
IJS, KISLJ, NUK, SBMB, UL, UM, UPUK
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.
Full text
Available for:
FZAB, GEOZS, IJS, IMTLJ, KILJ, KISLJ, MFDPS, NLZOH, NUK, OILJ, PNG, SAZU, SBCE, SBJE, SBMB, SBNM, UKNU, UL, UM, UPUK, VKSCE, ZAGLJ
Regulatory T cells (Tregs) can impair anti-tumor immune responses and are associated with poor prognosis in multiple cancer types. Tregs in human tumors span diverse transcriptional states distinct ...from those of peripheral Tregs, but their contribution to tumor development remains unknown. Here, we use single-cell RNA sequencing (RNA-seq) to longitudinally profile dynamic shifts in the distribution of Tregs in a genetically engineered mouse model of lung adenocarcinoma. In this model, interferon-responsive Tregs are more prevalent early in tumor development, whereas a specialized effector phenotype characterized by enhanced expression of the interleukin-33 receptor ST2 is predominant in advanced disease. Treg-specific deletion of ST2 alters the evolution of effector Treg diversity, increases infiltration of CD8+ T cells into tumors, and decreases tumor burden. Our study shows that ST2 plays a critical role in Treg-mediated immunosuppression in cancer, highlighting potential paths for therapeutic intervention.
Display omitted
•Single-cell profiling of CD4 T cells along tumor development in a mouse model•Treg diversity shifts to a Klrg1+Areg+ (KA) effector phenotype in advanced tumors•Il1rl1 (encoding ST2)+ Tregs have higher expression of KA effector Treg genes•Treg-specific ST2 loss enhances CD8+ T cell infiltration and decreases tumor burden
Li et al. show in a genetic mouse model of lung adenocarcinoma that during tumor development regulatory T cell (Treg) diversity shifts from an interferon-responsive to a ST2-positive, Klrg1+Areg+ effector-like phenotype. Treg-specific deletion of ST2 alters Treg heterogeneity, increases tumor infiltration by CD8+ T cells, and decreases tumor burden.
Full text
Available for:
GEOZS, IJS, IMTLJ, KILJ, KISLJ, NLZOH, NUK, OILJ, PNG, SAZU, SBCE, SBJE, UILJ, UL, UM, UPCLJ, UPUK, ZAGLJ, ZRSKP