Single-Nucleus and In Situ RNA-Sequencing Reveal Cell Topographies in the Human Pancreas.

BACKGROUND & AIMS: Molecular evidence of cellular heterogeneity in the human exocrine pancreas has not been yet established because of the local concentration and cascade of hydrolytic enzymes that can rapidly degrade cells and RNA upon pancreatic resection. We sought to better understand the heterogeneity and cellular composition of the pancreas in neonates and adults in healthy and diseased conditions using single-cell sequencing approaches. METHODS: We innovated single-nucleus RNA-sequencing protocols and profiled more than 120,000 cells from pancreata of adult and neonatal human donors. We validated the single-nucleus findings using RNA fluorescence in situ hybridization, in situ sequencing, and computational approaches. RESULTS: We created the first comprehensive atlas of human pancreas cells including epithelial and nonepithelial constituents, and uncovered 3 distinct acinar cell types, with possible implications for homeostatic and inflammatory processes of the pancreas. The comparison with neonatal single-nucleus sequencing data showed a different cellular composition of the endocrine tissue, highlighting the tissue dynamics occurring during development. By applying spatial cartography, involving cell proximity mapping through in situ sequencing, we found evidence of specific cell type neighborhoods, dynamic topographies in the endocrine and exocrine pancreas, and principles of morphologic organization of the organ. Furthermore, similar analyses in chronic pancreatitis biopsy samples showed the presence of acinar-REG+ cells, a reciprocal association between macrophages and activated stellate cells, and a new potential role of tuft cells in this disease. CONCLUSIONS: Our human pancreas cell atlas can be interrogated to understand pancreatic cell biology and provides a crucial reference set for comparisons with diseased tissue samples to map the cellular foundations of pancreatic diseases.

Proteomic and transcriptomic profiling of brainstem, cerebellum and olfactory tissues in early- and late-phase COVID-19.

Neurological symptoms, including cognitive impairment and fatigue, can occur in both the acute infection phase of coronavirus disease 2019 (COVID-19) and at later stages, yet the mechanisms that contribute to this remain unclear. Here we profiled single-nucleus transcriptomes and proteomes of brainstem tissue from deceased individuals at various stages of COVID-19. We detected an inflammatory type I interferon response in acute COVID-19 cases, which resolves in the late disease phase. Integrating single-nucleus RNA sequencing and spatial transcriptomics, we could localize two patterns of reaction to severe systemic inflammation, one neuronal with a direct focus on cranial nerve nuclei and a separate diffuse pattern affecting the whole brainstem. The latter reflects a bystander effect of the respiratory infection that spreads throughout the vascular unit and alters the transcriptional state of mainly oligodendrocytes, microglia and astrocytes, while alterations of the brainstem nuclei could reflect the connection of the immune system and the central nervous system via, for example, the vagus nerve. Our results indicate that even without persistence of severe acute respiratory syndrome coronavirus 2 in the central nervous system, local immune reactions are prevailing, potentially causing functional disturbances that contribute to neurological complications of COVID-19.

resVAE ensemble: Unsupervised identification of gene sets in multi-modal single-cell sequencing data using deep ensembles.

Feature identification and manual inspection is currently still an integral part of biological data analysis in single-cell sequencing. Features such as expressed genes and open chromatin status are selectively studied in specific contexts, cell states or experimental conditions. While conventional analysis methods construct a relatively static view on gene candidates, artificial neural networks have been used to model their interactions after hierarchical gene regulatory networks. However, it is challenging to identify consistent features in this modeling process due to the inherently stochastic nature of these methods. Therefore, we propose using ensembles of autoencoders and subsequent rank aggregation to extract consensus features in a less biased manner. Here, we performed sequencing data analyses of different modalities either independently or simultaneously as well as with other analysis tools. Our resVAE ensemble method can successfully complement and find additional unbiased biological insights with minimal data processing or feature selection steps while giving a measurement of confidence, especially for models using stochastic or approximation algorithms. In addition, our method can also work with overlapping clustering identity assignment suitable for transitionary cell types or cell fates in comparison to most conventional tools.

Single-cell analysis of patient-derived PDAC organoids reveals cell state heterogeneity and a conserved developmental hierarchy.

Pancreatic ductal adenocarcinoma (PDAC) is projected to be the second leading cause of cancer mortality by 2030. Bulk transcriptomic analyses have distinguished 'classical' from 'basal-like' tumors with more aggressive clinical behavior. We derive PDAC organoids from 18 primary tumors and two matched liver metastases, and show that 'classical' and 'basal-like' cells coexist in individual organoids. By single-cell transcriptome analysis of PDAC organoids and primary PDAC, we identify distinct tumor cell states shared across patients, including a cycling progenitor cell state and a differentiated secretory state. Cell states are connected by a differentiation hierarchy, with 'classical' cells concentrated at the endpoint. In an imaging-based drug screen, expression of 'classical' subtype genes correlates with better drug response. Our results thus uncover a functional hierarchy of PDAC cell states linked to transcriptional tumor subtypes, and support the use of PDAC organoids as a clinically relevant model for in vitro studies of tumor heterogeneity.

Functional States in Tumor-Initiating Cell Differentiation in Human Colorectal Cancer.

Intra-tumor heterogeneity of tumor-initiating cell (TIC) activity drives colorectal cancer (CRC) progression and therapy resistance. Here, we used single-cell RNA-sequencing of patient-derived CRC models to decipher distinct cell subpopulations based on their transcriptional profiles. Cell type-specific expression modules of stem-like, transit amplifying-like, and differentiated CRC cells resemble differentiation states of normal intestinal epithelial cells. Strikingly, identified subpopulations differ in proliferative activity and metabolic state. In summary, we here show at single-cell resolution that transcriptional heterogeneity identifies functional states during TIC differentiation. Furthermore, identified expression signatures are linked to patient prognosis. Targeting transcriptional states associated to cancer cell differentiation might unravel novel vulnerabilities in human CRC.