A spatial cell atlas of neuroblastoma reveals developmental, epigenetic and spatial axis of tumor heterogeneity

Neuroblastoma is a pediatric cancer arising from the developing sympathoadrenal lineage with complex inter- and intra-tumoral heterogeneity. To chart this complexity, we generated a comprehensive cell atlas of 55 neuroblastoma patient tumors, collected from two pediatric cancer institutions, spanning a range of clinical, genetic, and histologic features. Our atlas combines single-cell/nucleus RNA-seq (sc/scRNA-seq), bulk RNA-seq, whole exome sequencing, DNA methylation profiling, spatial transcriptomics, and two spatial proteomic methods. Sc/snRNA-seq revealed three malignant cell states with features of sympathoadrenal lineage development. All of the neuroblastomas had malignant cells that resembled sympathoblasts and the more differentiated adrenergic cells. A subset of tumors had malignant cells in a mesenchymal cell state with molecular features of Schwann cell precursors. DNA methylation profiles defined four groupings of patients, which differ in the degree of malignant cell heterogeneity and clinical outcomes. Using spatial proteomics, we found that neuroblastomas are spatially compartmentalized, with malignant tumor cells sequestered away from immune cells. Finally, we identify spatially restricted signaling patterns in immune cells from spatial transcriptomics. To facilitate the visualization and analysis of our atlas as a resource for further research in neuroblastoma, single cell, and spatial-omics, all data are shared through the Human Tumor Atlas Network Data Commons at www.humantumoratlas.org.

A single-nucleus and spatial transcriptomic atlas of the COVID-19 liver reveals topological, functional, and regenerative organ disruption in patients

The molecular underpinnings of organ dysfunction in acute COVID-19 and its potential long-term sequelae are under intense investigation. To shed light on these in the context of liver function, we performed single-nucleus RNA-seq and spatial transcriptomic profiling of livers from 17 COVID-19 decedents. We identified hepatocytes positive for SARS-CoV-2 RNA with an expression phenotype resembling infected lung epithelial cells. Integrated analysis and comparisons with healthy controls revealed extensive changes in the cellular composition and expression states in COVID-19 liver, reflecting hepatocellular injury, ductular reaction, pathologic vascular expansion, and fibrogenesis. We also observed Kupffer cell proliferation and erythrocyte progenitors for the first time in a human liver single-cell atlas, resembling similar responses in liver injury in mice and in sepsis, respectively. Despite the absence of a clinical acute liver injury phenotype, endothelial cell composition was dramatically impacted in COVID-19, concomitantly with extensive alterations and profibrogenic activation of reactive cholangiocytes and mesenchymal cells. Our atlas provides novel insights into liver physiology and pathology in COVID-19 and forms a foundational resource for its investigation and understanding.

A single-nucleus and spatial transcriptomic atlas of the COVID-19 liver reveals topological, functional, and regenerative organ disruption in patients

The molecular underpinnings of organ dysfunction in acute COVID-19 and its potential long-term sequelae are under intense investigation. To shed light on these in the context of liver function, we performed single-nucleus RNA-seq and spatial transcriptomic profiling of livers from 17 COVID-19 decedents. We identified hepatocytes positive for SARS-CoV-2 RNA with an expression phenotype resembling infected lung epithelial cells. Integrated analysis and comparisons with healthy controls revealed extensive changes in the cellular composition and expression states in COVID-19 liver, reflecting hepatocellular injury, ductular reaction, pathologic vascular expansion, and fibrogenesis. We also observed Kupffer cell proliferation and erythrocyte progenitors for the first time in a human liver single-cell atlas, resembling similar responses in liver injury in mice and in sepsis, respectively. Despite the absence of a clinical acute liver injury phenotype, endothelial cell composition was dramatically impacted in COVID-19, concomitantly with extensive alterations and profibrogenic activation of reactive cholangiocytes and mesenchymal cells. Our atlas provides novel insights into liver physiology and pathology in COVID-19 and forms a foundational resource for its investigation and understanding.

Mouse fetal growth restriction through parental and fetal immune gene variation and intercellular communications cascade

Fetal growth restriction (FGR) affects 5-10% of pregnancies, and can have serious consequences for both mother and child. Prevention and treatment are limited because FGR pathogenesis is poorly understood. Genetic studies implicate KIR and HLA genes in FGR, however, linkage disequilibrium, genetic influence from both parents, and challenges with investigating human pregnancies make the risk alleles and their functional effects difficult to map. Here, we demonstrate that the interaction between the maternal KIR2DL1, expressed on uterine natural killer (NK) cells, and the paternally inherited HLA-C*0501, expressed on fetal trophoblast cells, leads to FGR in a humanized mouse model. We show that the KIR2DL1 and C*0501 interaction leads to pathogenic uterine arterial remodeling and modulation of uterine NK cell function. This initial effect cascades to altered transcriptional expression and intercellular communication at the maternal-fetal interface. These findings provide mechanistic insight into specific FGR risk alleles, and provide avenues of prevention and treatment.

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.

Increased energy demand from anabolic-catabolic processes drives β-lactam antibiotic lethality

β-Lactam antibiotics disrupt the assembly of peptidoglycan (PG) within the bacterial cell wall by inhibiting the enzymatic activity of penicillin-binding proteins (PBPs). It was recently shown that β-lactam treatment initializes a futile cycle of PG synthesis and degradation, highlighting major gaps in our understanding of the lethal effects of PBP inhibition by β-lactam antibiotics. Here, we assess the downstream metabolic consequences of treatment of Escherichia coli with the β-lactam mecillinam and show that lethality from PBP2 inhibition is a specific consequence of toxic metabolic shifts induced by energy demand from multiple catabolic and anabolic processes, including accelerated protein synthesis downstream of PG futile cycling. Resource allocation into these processes is coincident with alterations in ATP synthesis and utilization, as well as a broadly dysregulated cellular redox environment. These results indicate that the disruption of normal anabolic-catabolic homeostasis by PBP inhibition is an essential factor for β-lactam antibiotic lethality.

Skin-resident innate lymphoid cells converge on a pathogenic effector state

Tissue-resident innate lymphoid cells (ILCs) help sustain barrier function and respond to local signals. ILCs are traditionally classified as ILC1, ILC2 or ILC3 on the basis of their expression of specific transcription factors and cytokines¹. In the skin, disease-specific production of ILC3-associated cytokines interleukin (IL)-17 and IL-22 in response to IL-23 signalling contributes to dermal inflammation in psoriasis. However, it is not known whether this response is initiated by pre-committed ILCs or by cell-state transitions. Here we show that the induction of psoriasis in mice by IL-23 or imiquimod reconfigures a spectrum of skin ILCs, which converge on a pathogenic ILC3-like state. Tissue-resident ILCs were necessary and sufficient, in the absence of circulatory ILCs, to drive pathology. Single-cell RNA-sequencing (scRNA-seq) profiles of skin ILCs along a time course of psoriatic inflammation formed a dense transcriptional continuum-even at steady state-reflecting fluid ILC states, including a naive or quiescent-like state and an ILC2 effector state. Upon disease induction, the continuum shifted rapidly to span a mixed, ILC3-like subset also expressing cytokines characteristic of ILC2s, which we inferred as arising through multiple trajectories. We confirmed the transition potential of quiescent-like and ILC2 states using in vitro experiments, single-cell assay for transposase-accessible chromatin using sequencing (scATAC-seq) and in vivo fate mapping. Our results highlight the range and flexibility of skin ILC responses, suggesting that immune activities primed in healthy tissues dynamically adapt to provocations and, left unchecked, drive pathological remodelling.

A Distinct Transcriptional Program in Human CAR T Cells Bearing the 4-1BB Signaling Domain Revealed by scRNA-Seq

T cells engineered to express chimeric antigen receptors (CARs) targeting CD19 have produced impressive outcomes for the treatment of B cell malignancies, but different products vary in kinetics, persistence, and toxicity profiles based on the co-stimulatory domains included in the CAR. In this study, we performed transcriptional profiling of bulk CAR T cell populations and single cells to characterize the transcriptional states of human T cells transduced with CD3ζ, 4-1BB-CD3ζ (BBζ), or CD28-CD3ζ (28ζ) co-stimulatory domains at rest and after activation by triggering their CAR or their endogenous T cell receptor (TCR). We identified a transcriptional signature common across CARs with the CD3ζ signaling domain, as well as a distinct program associated with the 4-1BB co-stimulatory domain at rest and after activation. CAR T cells bearing BBζ had increased expression of human leukocyte antigen (HLA) class II genes, ENPP2, and interleukin (IL)-21 axis genes, and decreased PD1 compared to 28ζ CAR T cells. Similar to previous studies, we also found BBζ CAR CD8 T cells to be enriched in a central memory cell phenotype and fatty acid metabolism genes. Our data uncovered transcriptional signatures related to costimulatory domains and demonstrated that signaling domains included in CARs uniquely shape the transcriptional programs of T cells.

Emergence of a High-Plasticity Cell State during Lung Cancer Evolution

Tumor evolution from a single cell into a malignant, heterogeneous tissue remains poorly understood. Here, we profile single-cell transcriptomes of genetically engineered mouse lung tumors at seven stages, from pre-neoplastic hyperplasia to adenocarcinoma. The diversity of transcriptional states increases over time and is reproducible across tumors and mice. Cancer cells progressively adopt alternate lineage identities, computationally predicted to be mediated through a common transitional, high-plasticity cell state (HPCS). Accordingly, HPCS cells prospectively isolated from mouse tumors and human patient-derived xenografts display high capacity for differentiation and proliferation. The HPCS program is associated with poor survival across human cancers and demonstrates chemoresistance in mice. Our study reveals a central principle underpinning intra-tumoral heterogeneity and motivates therapeutic targeting of the HPCS.

A single-cell landscape of high-grade serous ovarian cancer

Malignant abdominal fluid (ascites) frequently develops in women with advanced high-grade serous ovarian cancer (HGSOC) and is associated with drug resistance and a poor prognosis1. To comprehensively characterize the HGSOC ascites ecosystem, we used single-cell RNA sequencing to profile ~11,000 cells from 22 ascites specimens from 11 patients with HGSOC. We found significant inter-patient variability in the composition and functional programs of ascites cells, including immunomodulatory fibroblast sub-populations and dichotomous macrophage populations. We found that the previously described immunoreactive and mesenchymal subtypes of HGSOC, which have prognostic implications, reflect the abundance of immune infiltrates and fibroblasts rather than distinct subsets of malignant cells2. Malignant cell variability was partly explained by heterogeneous copy number alteration patterns or expression of a stemness program. Malignant cells shared expression of inflammatory programs that were largely recapitulated in single-cell RNA sequencing of ~35,000 cells from additionally collected samples, including three ascites, two primary HGSOC tumors and three patient ascites-derived xenograft models. Inhibition of the JAK/STAT pathway, which was expressed in both malignant cells and cancer-associated fibroblasts, had potent anti-tumor activity in primary short-term cultures and patient-derived xenograft models. Our work contributes to resolving the HSGOC landscape3-5 and provides a resource for the development of novel therapeutic approaches.

Author Correction: A single-cell and single-nucleus RNA-Seq toolbox for fresh and frozen human tumors

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

A single-cell and single-nucleus RNA-Seq toolbox for fresh and frozen human tumors

Single-cell genomics is essential to chart tumor ecosystems. Although single-cell RNA-Seq (scRNA-Seq) profiles RNA from cells dissociated from fresh tumors, single-nucleus RNA-Seq (snRNA-Seq) is needed to profile frozen or hard-to-dissociate tumors. Each requires customization to different tissue and tumor types, posing a barrier to adoption. Here, we have developed a systematic toolbox for profiling fresh and frozen clinical tumor samples using scRNA-Seq and snRNA-Seq, respectively. We analyzed 216,490 cells and nuclei from 40 samples across 23 specimens spanning eight tumor types of varying tissue and sample characteristics. We evaluated protocols by cell and nucleus quality, recovery rate and cellular composition. scRNA-Seq and snRNA-Seq from matched samples recovered the same cell types, but at different proportions. Our work provides guidance for studies in a broad range of tumors, including criteria for testing and selecting methods from the toolbox for other tumors, thus paving the way for charting tumor atlases.

Transcriptional Atlas of Intestinal Immune Cells Reveals that Neuropeptide α-CGRP Modulates Group 2 Innate Lymphoid Cell Responses

Signaling abnormalities in immune responses in the small intestine can trigger chronic type 2 inflammation involving interaction of multiple immune cell types. To systematically characterize this response, we analyzed 58,067 immune cells from the mouse small intestine by single-cell RNA sequencing (scRNA-seq) at steady state and after induction of a type 2 inflammatory reaction to ovalbumin (OVA). Computational analysis revealed broad shifts in both cell-type composition and cell programs in response to the inflammation, especially in group 2 innate lymphoid cells (ILC2s). Inflammation induced the expression of exon 5 of Calca, which encodes the alpha-calcitonin gene-related peptide (α-CGRP), in intestinal KLRG1⁺ ILC2s. α-CGRP antagonized KLRG1⁺ ILC2s proliferation but promoted IL-5 expression. Genetic perturbation of α-CGRP increased the proportion of intestinal KLRG1⁺ ILC2s. Our work highlights a model where α-CGRP-mediated neuronal signaling is critical for suppressing ILC2 expansion and maintaining homeostasis of the type 2 immune machinery.

Author Correction: Design, execution, and analysis of CRISPR-Cas9-based deletions and genetic interaction networks in the fungal pathogen Candida albicans

The version of this paper originally published contained reference errors. The sentence "To dissect complex genetic interactions in C. albicans, a CRISPR-Cas9-based Gene Drive Array (GDA) was developed" incorrectly cited ref. 13, and should have cited ref. 14. In addition, the reference included as ref. 13 in the original paper was incorrect, and should have been the following: Shapiro, R. S., Chavez, A. & Collins, J. J. CRISPR-based genomic tools for the manipulation of genetically intractable microorganisms. Nat. Rev. Microbiol. 16, 333-339 (2018). This reference should have been cited after the sentence "Recent innovations in CRISPR-Cas9-based genome editing have facilitated such genetic interaction analyses." The original reference 13 (Gerami-Nejad, M., Zacchi, L. F., McClellan, M., Matter, K. & Berman, J. Shuttle vectors for facile gap repair cloning and integration into a neutral locus in Candida albicans. Microbiology 159, 565-579 (2013)) should have been cited later in the paper, and is now in the reference list as ref. 27. As a result, original references 27-33 have been renumbered in the reference list and in the text. These changes have been made in the PDF and HTML versions of the protocol.

Design, execution, and analysis of CRISPR-Cas9-based deletions and genetic interaction networks in the fungal pathogen Candida albicans

The study of fungal pathogens is of immediate importance, yet progress is hindered by the technical challenges of genetic manipulation. For Candida species, their inability to maintain plasmids, unusual codon usage, and inefficient homologous recombination are among the obstacles limiting efficient genetic manipulation. New advances in genomic biotechnologies-particularly CRISPR-based tools-have revolutionized genome editing for many fungal species. Here, we present a protocol for CRISPR-Cas9-based manipulation in Candida albicans using a modified gene-drive-based strategy that takes ~1 month to complete. We detail the generation of Candida-optimized Cas9-based plasmids for gene deletion, an efficient transformation protocol using C. albicans haploids, and an optimized mating strategy to generate homozygous single- and double-gene diploid mutants. We further describe protocols for quantifying cell growth and analysis pipelines to calculate fitness and genetic interaction scores for genetic mutants. This protocol overcomes previous limitations associated with genetic manipulation in C. albicans and advances researchers' ability to perform genetic analysis in this pathogen; the protocol also has broad applicability to other mating-competent microorganisms.

A CRISPR-Cas9-based gene drive platform for genetic interaction analysis in Candida albicans

Candida albicans is the leading cause of fungal infections; yet, complex genetic interaction analysis remains cumbersome in this diploid pathogen. Here, we developed a CRISPR-Cas9-based 'gene drive array' platform to facilitate efficient genetic analysis in C. albicans. In our system, a modified DNA donor molecule acts as a selfish genetic element, replaces the targeted site and propagates to replace additional wild-type loci. Using mating-competent C. albicans haploids, each carrying a different gene drive disabling a gene of interest, we are able to create diploid strains that are homozygous double-deletion mutants. We generate double-gene deletion libraries to demonstrate this technology, targeting antifungal efflux and biofilm adhesion factors. We screen these libraries to identify virulence regulators and determine how genetic networks shift under diverse conditions. This platform transforms our ability to perform genetic interaction analysis in C. albicans and is readily extended to other fungal pathogens.

Carbon Sources Tune Antibiotic Susceptibility in Pseudomonas aeruginosa via Tricarboxylic Acid Cycle Control

Metabolically dormant bacteria present a critical challenge to effective antimicrobial therapy because these bacteria are genetically susceptible to antibiotic treatment but phenotypically tolerant. Such tolerance has been attributed to impaired drug uptake, which can be reversed by metabolic stimulation. Here, we evaluate the effects of central carbon metabolite stimulations on aminoglycoside sensitivity in the pathogen Pseudomonas aeruginosa. We identify fumarate as a tobramycin potentiator that activates cellular respiration and generates a proton motive force by stimulating the tricarboxylic acid (TCA) cycle. In contrast, we find that glyoxylate induces phenotypic tolerance by inhibiting cellular respiration with acetyl-coenzyme A diversion through the glyoxylate shunt, despite drug import. Collectively, this work demonstrates that TCA cycle activity is important for both aminoglycoside uptake and downstream lethality and identifies a potential strategy for potentiating aminoglycoside treatment of P. aeruginosa infections.

Bactericidal Antibiotics Induce Toxic Metabolic Perturbations that Lead to Cellular Damage

Understanding how antibiotics impact bacterial metabolism may provide insight into their mechanisms of action and could lead to enhanced therapeutic methodologies. Here, we profiled the metabolome of Escherichia coli after treatment with three different classes of bactericidal antibiotics (?-lactams, aminoglycosides, quinolones). These treatments induced a similar set of metabolic changes after 30 min that then diverged into more distinct profiles at later time points. The most striking changes corresponded to elevated concentrations of central carbon metabolites, active breakdown of the nucleotide pool, reduced lipid levels, and evidence of an elevated redox state. We examined potential end-target consequences of these metabolic perturbations and found that antibiotic-treated cells exhibited cytotoxic changes indicative of oxidative stress, including higher levels of protein carbonylation, malondialdehyde adducts, nucleotide oxidation, and double-strand DNA breaks. This work shows that bactericidal antibiotics induce a complex set of metabolic changes that are correlated with the buildup of toxic metabolic by-products.