Fibrosis is a major complication of Crohn's disease (CD) marked by excess deposition of extracellular matrix, leading to stricturing and functional impairment. As mechanistic characterization and therapeutic options are lacking, we paired single-cell and spatial transcriptomics in 61 samples from 21 patients with CD and 10 patients without inflammatory bowel disease (IBD). Intestinal strictures were characterized by increased immune cells, including IgG+ plasma cells, CCR7-hi CD4+ T cells and inflammatory fibroblasts. Spatial transcriptomics showed that key subsets colocalize within diseased tissues and identified additional populations such as interstitial cells of Cajal and enteric neurons. Furthermore, we mapped gene expression onto intestinal biogeography, finding that known genetic risk loci are enriched within discrete spatial modules, defined by the presence of inflammatory fibroblasts and lymphoid follicles. Altogether, our datasets chart the key transcriptomic and cellular networks in stricturing CD and highlight the spatial organization of multicellular genetic risk factors.

Antibody-secreting plasma cells (PCs) are generated in secondary lymphoid organs but are reported to reside in an emerging range of anatomical sites. Analysis of the transcriptome of different tissue-resident (Tr)PC populations revealed that they each have their own transcriptional signature indicative of functional adaptation to the host tissue environment. In contrast to expectation, all TrPCs were extremely long-lived, regardless of their organ of residence, with longevity influenced by intrinsic factors like the immunoglobulin isotype. Analysis at single-cell resolution revealed that the bone marrow is unique in housing a compendium of PCs generated all over the body that retain aspects of the transcriptional program indicative of their tissue of origin. This study reveals that extreme longevity is an intrinsic property of TrPCs whose transcriptome is imprinted by signals received both at the site of induction and within the tissue of residence.

Genome-wide association studies (GWASs) are a valuable tool for understanding the biology of complex human traits and diseases, but associated variants rarely point directly to causal genes. In the present study, we introduce a new method, polygenic priority score (PoPS), that learns trait-relevant gene features, such as cell-type-specific expression, to prioritize genes at GWAS loci. Using a large evaluation set of genes with fine-mapped coding variants, we show that PoPS and the closest gene individually outperform other gene prioritization methods, but observe the best overall performance by combining PoPS with orthogonal methods. Using this combined approach, we prioritize 10,642 unique gene–trait pairs across 113 complex traits and diseases with high precision, finding not only well-established gene–trait relationships but nominating new genes at unresolved loci, such as LGR4 for estimated glomerular filtration rate and CCR7 for deep vein thrombosis. Overall, we demonstrate that PoPS provides a powerful addition to the gene prioritization toolbox. Polygenic Priority Score (PoPS) prioritizes candidate effector genes at complex trait loci by integrating genome-wide association summary statistics with other data types. Combining PoPS with methods that leverage local genetic signals further improves the performance.

We introduce BacDrop, a highly scalable technology for bacterial single-cell RNA sequencing that has overcome many challenges hindering the development of scRNA-seq in bacteria. BacDrop can be applied to thousands to millions of cells from both gram-negative and gram-positive species. It features universal ribosomal RNA depletion and combinatorial barcodes that enable multiplexing and massively parallel sequencing. We applied BacDrop to study Klebsiella pneumoniae clinical isolates and to elucidate their heterogeneous responses to antibiotic stress. In an unperturbed population presumed to be homogeneous, we found within-population heterogeneity largely driven by the expression of mobile genetic elements that promote the evolution of antibiotic resistance. Under antibiotic perturbation, BacDrop revealed transcriptionally distinct subpopulations associated with different phenotypic outcomes including antibiotic persistence. BacDrop thus can capture cellular states that cannot be detected by bulk RNA-seq, which will unlock new microbiological insights into bacterial responses to perturbations and larger bacterial communities such as the microbiome.

Immune checkpoint blockade (ICB) has revolutionized the treatment of many cancers but has been ineffective for the treatment of microsatellite stable (MSS) PDAC. The lack of efficacy of immunotherapies in PDAC is due to: 1) a desmoplastic tumor microenvironment (TME); 2) the presence of suppressive cells, including myeloid derived suppressor cells and regulatory T cells; and 3) the lack of antigen-presenting dendritic cells (DCs) that are important in priming an effective immune response to generate functionally effective tumor antigen-specific T cells. We recently completed a pilot study of dual ICB (Ipilumamab and Nivolumab) with radiation therapy (SBRT 8Gy for 3 fractions) in a cohort of 25 metastatic PDAC patients that had progressed on conventional chemotherapy; this combination conferred an impressive 18% ORR and 29% disease control rate measured on non-irradiated lesions (historical 0% ORR with ICB in PDAC). This led to a phase 2 study in 30 metastatic PDAC patients using this dual modality treatment paradigm. To understand the role of radiation and ICB in altering the PDAC tumor microenvironment we performed single-cell RNA-sequencing and TCR-sequencing (>180k cells), and single-nucleus RNA-sequencing (>300k cells) on 36 tumor biopsies (23 pre-treatment, 13 paired on-treatment between day 10 and 21) from patients undergoing treatment in our phase 2 study. Tumor tissue was taken from distinct tissue sites, including primary tumors in the pancreas, and liver and abdominal wall metastases. We identified distinct tumor cell state distributions within different tissues, and a redistribution of cells from basal/mesenchymal states to classical states after radiation. We identified several state-specific interferon stimulated gene programs thus cataloging distinct responses of epithelial cells with different transcriptional states. Importantly, we found a redistribution of T cells states towards proliferating and exhausted T cells with unique clonality after radiation. Additionally, the myeloid compartment after radiation was enriched for C1QC+ and MHCII+ macrophage subsets, as well as infiltrating CD16/CD16 monocytes and CD14 monocytes, each showing induction of unique sets of interferon stimulated genes (ISGs). We next sought to better understand immunotherapy resistance mechanisms within these PDAC patients despite finding strong ISG induction in several subsets. We analyzed covarying gene programs and identified multicellular communities of cells before and after radiation that underlie interaction networks associated with radiation. Together our data provides the most comprehensive single-cell atlas of paired biopsies to study tumor and immune cell states in the context of radiation and ICB response. Citation Format: Arnav Mehta, Aparna Parikh, Milan Parikh, Ryan Park, Moshe Sade-Feldman, Lynn Bi, Nicole Carzo, Tarin M. Grillo, Islam Baiev, Olanike Asupoto, Irena Gushterova, Tom LaSalle, Anna Gonye, Emily Blaum, Sebastien Vigneau, Ronan Chaligne, Ana Lako, Thomas Lila, David Nelson, Caroline Porter, Orr Ashenberg, Karthik Jagadesh, William L. Hwang, Christopher Smillie, David P. Ryan, David T. Ting, Theodore Hong, Dana Pe'er, Nir Hacohen. Dissecting the reorganization of pancreatic tumor microenvironments after radiation and immunotherapy reveals insights into immunotherapy resistance [abstract]. In: Proceedings of the AACR Special Conference on Pancreatic Cancer; 2022 Sep 13-16; Boston, MA. Philadelphia (PA): AACR; Cancer Res 2022;82(22 Suppl):Abstract nr C012.

BACKGROUND & AIMS:Monogenic forms of inflammatory bowel disease (IBD) illustrate the essential roles of individual genes in pathways and networks safeguarding immune tolerance and gut homeostasis. METHODS:To build a taxonomy model, we assessed 165 disorders. Genes were prioritized based on penetrance of IBD and disease phenotypes were integrated with multi-omics datasets. Monogenic IBD genes were classified by (1) overlapping syndromic features, (2) response to hematopoietic stem cell transplantation, (3) bulk RNA-sequencing of 32 tissues, (4) single-cell RNA-sequencing of >50 cell subsets from the intestine of healthy individuals and patients with IBD (pediatric and adult), and (5) proteomes of 43 immune subsets. The model was validated by addition of newly identified monogenic IBD defects. As a proof-of-concept, we explore the intersection between immunometabolism and antimicrobial activity for a group of disorders (G6PC3/SLC37A4). RESULTS:Our quantitative integrated taxonomy defines the cellular landscape of monogenic IBD gene expression across 102 genes with high and moderate penetrance (81 in the model set and 21 genes in the validation set). We illustrate distinct cellular networks, highlight expression profiles across understudied cell types (e.g., CD8+ T cells, neutrophils, epithelial subsets, and endothelial cells) and define genotype-phenotype associations (perianal disease and defective antimicrobial activity). We illustrate processes and pathways shared across cellular compartments and phenotypic groups and highlight cellular immunometabolism with mammalian target of rapamycin activation as one of the converging pathways. There is an overlap of genes and enriched cell-specific expression between monogenic and polygenic IBD. CONCLUSION:Our taxonomy integrates genetic, clinical and multi-omic data; providing a basis for genomic diagnostics and testable hypotheses for disease functions and treatment responses.

We introduce BacDrop, a highly scalable technology for bacterial single-cell RNA sequencing that has overcome many challenges hindering the development of scRNA-seq in bacteria. BacDrop can be applied to thousands or millions of cells from both gram-negative and gram-positive species. It features universal ribosomal RNA depletion and combinatorial barcodes that enable multiplexing and massively parallel sequencing. We applied BacDrop to study Klebsiella pneumoniae clinical isolates and to elucidate their heterogeneous responses to antibiotic stress. In an unperturbed population presumed to be homogenous, we find within- population heterogeneity largely driven by the expression of mobile genetic elements that promote the evolution of antibiotic resistance. Under antibiotic perturbation, BacDrop revealed transcriptionally distinct subpopulations associated with different phenotypic outcomes including antibiotic persistence. BacDrop thus can capture cellular states that cannot be detected by bulk RNA-seq, which will unlock new microbiological insights into bacterial responses to perturbations and larger bacterial communities such as the microbiome.

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.

Mechanisms underlying severe coronavirus disease 2019 (COVID-19) disease remain poorly understood. We analyze several thousand plasma proteins longitudinally in 306 COVID-19 patients and 78 symptomatic controls, uncovering immune and non-immune proteins linked to COVID-19. Deconvolution of our plasma proteome data using published scRNA-seq datasets reveals contributions from circulating immune and tissue cells. Sixteen percent of patients display reduced inflammation yet comparably poor outcomes. Comparison of patients who died to severely ill survivors identifies dynamic immune-cell-derived and tissue-associated proteins associated with survival, including exocrine pancreatic proteases. Using derived tissue-specific and cell-type-specific intracellular death signatures, cellular angiotensin-converting enzyme 2 (ACE2) expression, and our data, we infer whether organ damage resulted from direct or indirect effects of infection. We propose a model in which interactions among myeloid, epithelial, and T cells drive tissue damage. These datasets provide important insights and a rich resource for analysis of mechanisms of severe COVID-19 disease.

The intestine is a site of direct encounter with the external environment and must consequently balance barrier defense with nutrient uptake. To investigate how nutrient uptake is regulated in the small intestine, we tested the effect of diets with different macronutrient compositions on epithelial gene expression. We found that enzymes and transporters required for carbohydrate digestion and absorption were regulated by carbohydrate availability. The "on-demand" induction of this machinery required γδ T cells, which regulated this program through the suppression of interleukin-22 production by type 3 innate lymphoid cells. Nutrient availability altered the tissue localization and transcriptome of γδ T cells. Additionally, transcriptional responses to diet involved cellular remodeling of the epithelial compartment. Thus, this work identifies a role for γδ T cells in nutrient sensing.

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. An integrated analysis of over 100 single-cell and single-nucleus transcriptomics studies illustrates severe acute respiratory syndrome coronavirus 2 viral entry gene coexpression patterns across different human tissues, and shows association of age, smoking status and sex with viral entry gene expression in respiratory cell populations.

Gut dysmotility is associated with constipation, diarrhea, and functional gastrointestinal disorders like irritable bowel syndrome (IBS), although its molecular underpinnings are poorly characterized. We studied stool frequency (defined by the number of bowel movements per day, based on questionnaire data) as a proxy for gut motility in a GWAS meta-analysis including 167,875 individuals from UK Biobank and four smaller population based cohorts. We identify 14 loci associated with stool frequency (p % 5.0 3 10-8). Gene set and pathway analyses detected enrichment for genes involved in neurotransmitter/neuropeptide signaling and preferentially expressed in enteric motor neurons controlling peristalsis. PheWAS identified pleiotropic associations with dysmotility syndromes and the response to their pharmacological treatment. The genetic architecture of stool frequency correlates with that of IBS, and UK Biobank participants from the top 1% of stool frequency polygenic score distribution were associated with 53 higher risk of IBS with diarrhea. These findings pave the way for the identification of actionable pathological mechanisms in IBS and the dysmotility syndromes.

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.

ABSTRACT The COVID-19 pandemic, caused by the novel coronavirus SARS-CoV-2, creates an urgent need for identifying molecular mechanisms that mediate viral entry, propagation, and tissue pathology. Cell membrane bound angiotensin-converting enzyme 2 (ACE2) and associated proteases, transmembrane protease serine 2 (TMPRSS2) and Cathepsin L (CTSL), were previously identified as mediators of SARS-CoV2 cellular entry. Here, we assess the cell type-specific RNA expression of ACE2 , TMPRSS2 , and CTSL through an integrated analysis of 107 single-cell and single-nucleus RNA-Seq studies, including 22 lung and airways datasets (16 unpublished), and 85 datasets from other diverse organs. Joint expression of ACE2 and the accessory proteases identifies specific subsets of respiratory epithelial cells as putative targets of viral infection in the nasal passages, airways, and alveoli. Cells that co-express ACE2 and proteases are also identified in cells from other organs, some of which have been associated with COVID-19 transmission or pathology, including gut enterocytes, corneal epithelial cells, cardiomyocytes, heart pericytes, olfactory sustentacular cells, and renal epithelial cells. Performing the first meta-analyses of scRNA-seq studies, we analyzed 1,176,683 cells from 282 nasal, airway, and lung parenchyma samples from 164 donors spanning fetal, childhood, adult, and elderly age groups, associate increased levels of ACE2 , TMPRSS2 , and CTSL in specific cell types with increasing age, male gender, and smoking, all of which are epidemiologically linked to COVID-19 susceptibility and outcomes. Notably, there was a particularly low expression of ACE2 in the few young pediatric samples in the analysis. Further analysis reveals a gene expression program shared by ACE2 + TMPRSS2 + cells in nasal, lung and gut tissues, including genes that may mediate viral entry, subtend key immune functions, and mediate epithelial-macrophage cross-talk. Amongst these are IL6, its receptor and co-receptor, IL1R , TNF response pathways, and complement genes. Cell type specificity in the lung and airways and smoking effects were conserved in mice. Our analyses suggest that differences in the cell type-specific expression of mediators of SARS-CoV-2 viral entry may be responsible for aspects of COVID-19 epidemiology and clinical course, and point to putative molecular pathways involved in disease susceptibility and pathogenesis.

While horizontal gene transfer is prevalent across the biosphere, the regulatory features that enable expression and functionalization of foreign DNA remain poorly understood. Here, we combine high-throughput promoter activity measurements and large-scale genomic analysis of regulatory regions to investigate the cross-compatibility of regulatory elements (REs) in bacteria. Functional characterization of thousands of natural REs in three distinct bacterial species revealed distinct expression patterns according to RE and recipient phylogeny. Host capacity to activate foreign promoters was proportional to their genomic GC content, while many low GC regulatory elements were both broadly active and had more transcription start sites across hosts. The difference in expression capabilities could be explained by the influence of the host GC content on the stringency of the AT-rich canonical σ70 motif necessary for transcription initiation. We further confirm the generalizability of this model and find widespread GC content adaptation of the σ70 motif in a set of 1,545 genomes from all major bacterial phyla. Our analysis identifies a key mechanism by which the strength of the AT-rich σ70 motif relative to a host’s genomic GC content governs the capacity for expression of acquired DNA. These findings shed light on regulatory adaptation in the context of evolving genomic composition.

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

COVID-19 has caused over 1 million deaths globally, yet the cellular mechanisms underlying severe disease remain poorly understood. By analyzing several thousand plasma proteins in 306 COVID-19 patients and 78 symptomatic controls over serial timepoints using two complementary approaches, we uncover COVID-19 host immune and non-immune proteins not previously linked to this disease. Integration of plasma proteomics with nine published scRNAseq datasets shows that SARS-CoV-2 infection upregulates monocyte/macrophage, plasmablast, and T cell effector proteins. By comparing patients who died to severely ill patients who survived, we identify dynamic immunomodulatory and tissue-associated proteins associated with survival, providing insights into which host responses are beneficial and which are detrimental to survival. We identify intracellular death signatures from specific tissues and cell types, and by associating these with angiotensin converting enzyme 2 (ACE2) expression, we map tissue damage associated with severe disease and propose which damage results from direct viral infection rather than from indirect effects of illness. We find that disease severity in lung tissue is driven by myeloid cell phenotypes and cell-cell interactions with lung epithelial cells and T cells. Based on these results, we propose a model of immune and epithelial cell interactions that drive cell-type specific and tissue-specific damage in severe COVID-19.

ABSTRACTObjectiveGut dysmotility is associated with constipation, diarrhea, and functional gastrointestinal disorders (FGID) like irritable bowel syndrome (IBS). Its molecular underpinnings, and their anomalies in FGID disorders are poorly characterized, hence we sought to gain mechanistic insight through a large-scale genetic investigation.DesignWe used stool frequency (STL-FRQ) as a (surrogate) quantitative trait to study the genetics of gut motility, exploiting questionnaire and genotype data from UK Biobank and four smaller population-based cohorts (LifeLines-Deep, Genes for Good, Flemish Gut Flora Project and PopCol), in a GWAS meta-analysis spanning 8,817,117 high-quality SNP markers and 167,875 individuals of European descent.ResultsWe identify 13 genome-wide significant loci (P≤5.0×10−8) harboring prioritized genes that are: i) involved in sensory perception and neurotransmitter/neuropeptide signaling; ii) enriched for their expression in enteric motor neurons associated with the control of peristalsis (P=7.0×10−8) iii) previously linked to other traits and conditions, including GI motility and dysmotility syndromes, and the response to their pharmacological treatment. The genetic architecture of STL-FRQ most strongly correlates with that of IBS (rg=0.42; P=1.1×10−3). In UK Biobank, the risk of IBS with diarrhea was 4x higher in individuals from the top 1% of the distribution of polygenic scores (PGS) computed based on STL-FRQ GWAS summary statistics (ORs=4.14; P=1.2×10−97).ConclusionWe identify loci harboring genes with a plausible role in GI motility, possibly acting via neurotransmission and similar pathways in specialized enteric neurons. The demonstrated relevance of these findings to IBS warrants further study for the identification of actionable pathomechanisms in the dysmotility syndromes.

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.