CellRegMap: a statistical framework for mapping context‐specific regulatory variants using scRNA‐seq

Single‐cell RNA sequencing (scRNA‐seq) enables characterizing the cellular heterogeneity in human tissues. Recent technological advances have enabled the first population‐scale scRNA‐seq studies in hundreds of individuals, allowing to assay genetic effects with single‐cell resolution. However, existing strategies to analyze these data remain based on principles established for the genetic analysis of bulk RNA‐seq. In particular, current methods depend on a priori definitions of discrete cell types, and hence cannot assess allelic effects across subtle cell types and cell states. To address this, we propose the Cell Regulatory Map (CellRegMap), a statistical framework to test for and quantify genetic effects on gene expression in individual cells. CellRegMap provides a principled approach to identify and characterize genotype–context interactions of known eQTL variants using scRNA‐seq data. This model‐based approach resolves allelic effects across cellular contexts of different granularity, including genetic effects specific to cell subtypes and continuous cell transitions. We validate CellRegMap using simulated data and apply it to previously identified eQTL from two recent studies of differentiating iPSCs, where we uncover hundreds of eQTL displaying heterogeneity of genetic effects across cellular contexts. Finally, we identify fine‐grained genetic regulation in neuronal subtypes for eQTL that are colocalized with human disease variants.     

Single‐cell dynamics of chromatin activity during cell lineage differentiation in Caenorhabditis elegans embryos

Elucidating the chromatin dynamics that orchestrate embryogenesis is a fundamental question in developmental biology. Here, we exploit position effects on expression as an indicator of chromatin activity and infer the chromatin activity landscape in every lineaged cell during Caenorhabditis elegans early embryogenesis. Systems‐level analyses reveal that chromatin activity distinguishes cellular states and correlates with fate patterning in the early embryos. As cell lineage unfolds, chromatin activity diversifies in a lineage‐dependent manner, with switch‐like changes accompanying anterior–posterior fate asymmetry and characteristic landscapes being established in different cell lineages. Upon tissue differentiation, cellular chromatin from distinct lineages converges according to tissue types but retains stable memories of lineage history, contributing to intra‐tissue cell heterogeneity. However, the chromatin landscapes of cells organized in a left–right symmetric pattern are predetermined to be analogous in early progenitors so as to pre‐set equivalent states. Finally, genome‐wide analysis identifies many regions exhibiting concordant chromatin activity changes that mediate the co‐regulation of functionally related genes during differentiation. Collectively, our study reveals the developmental and genomic dynamics of chromatin activity at the single‐cell level.     

Aging‐related cell type‐specific pathophysiologic immune responses that exacerbate disease severity in aged COVID‐19 patients

Coronavirus disease 2019 (COVID‐19) is especially severe in aged patients, defined as 65 years or older, for reasons that are currently unknown. To investigate the underlying basis for this vulnerability, we performed multimodal data analyses on immunity, inflammation, and COVID‐19 incidence and severity as a function of age. Our analysis leveraged age‐specific COVID‐19 mortality and laboratory testing from a large COVID‐19 registry, along with epidemiological data of ~3.4 million individuals, large‐scale deep immune cell profiling data, and single‐cell RNA‐sequencing data from aged COVID‐19 patients across diverse populations. We found that decreased lymphocyte count and elevated inflammatory markers (C‐reactive protein, D‐dimer, and neutrophil–lymphocyte ratio) are significantly associated with age‐specific COVID‐19 severities. We identified the reduced abundance of naïve CD8 T cells with decreased expression of antiviral defense genes (i.e., IFITM3 and TRIM22) in aged severe COVID‐19 patients. Older individuals with severe COVID‐19 displayed type I and II interferon deficiencies, which is correlated with SARS‐CoV‐2 viral load. Elevated expression of SARS‐CoV‐2 entry factors and reduced expression of antiviral defense genes (LY6E and IFNAR1) in the secretory cells are associated with critical COVID‐19 in aged individuals. Mechanistically, we identified strong TGF‐beta‐mediated immune–epithelial cell interactions (i.e., secretory‐non‐resident macrophages) in aged individuals with critical COVID‐19. Taken together, our findings point to immuno‐inflammatory factors that could be targeted therapeutically to reduce morbidity and mortality in aged COVID‐19 patients.     

An organoid‐derived bronchioalveolar model for SARS‐CoV‐2 infection of human alveolar type II‐like cells

Severe acute respiratory syndrome coronavirus 2 (SARS‐CoV‐2) causes coronavirus disease 2019 (COVID‐19), which may result in acute respiratory distress syndrome (ARDS), multiorgan failure, and death. The alveolar epithelium is a major target of the virus, but representative models to study virus host interactions in more detail are currently lacking. Here, we describe a human 2D air–liquid interface culture system which was characterized by confocal and electron microscopy and single‐cell mRNA expression analysis. In this model, alveolar cells, but also basal cells and rare neuroendocrine cells, are grown from 3D self‐renewing fetal lung bud tip organoids. These cultures were readily infected by SARS‐CoV‐2 with mainly surfactant protein C‐positive alveolar type II‐like cells being targeted. Consequently, significant viral titers were detected and mRNA expression analysis revealed induction of type I/III interferon response program. Treatment of these cultures with a low dose of interferon lambda 1 reduced viral replication. Hence, these cultures represent an experimental model for SARS‐CoV‐2 infection and can be applied for drug screens.     

Single‐cell transcriptomics of suprachiasmatic nuclei reveal a Prokineticin‐driven circadian network

Circadian rhythms in mammals are governed by the hypothalamic suprachiasmatic nucleus (SCN), in which 20,000 clock cells are connected together into a powerful time‐keeping network. In the absence of network‐level cellular interactions, the SCN fails as a clock. The topology and specific roles of its distinct cell populations (nodes) that direct network functions are, however, not understood. To characterise its component cells and network structure, we conducted single‐cell sequencing of SCN organotypic slices and identified eleven distinct neuronal sub‐populations across circadian day and night. We defined neuropeptidergic signalling axes between these nodes, and built neuropeptide‐specific network topologies. This revealed their temporal plasticity, being up‐regulated in circadian day. Through intersectional genetics and real‐time imaging, we interrogated the contribution of the Prok2‐ProkR2 neuropeptidergic axis to network‐wide time‐keeping. We showed that Prok2‐ProkR2 signalling acts as a key regulator of SCN period and rhythmicity and contributes to defining the network‐level properties that underpin robust circadian co‐ordination. These results highlight the diverse and distinct contributions of neuropeptide‐modulated communication of temporal information across the SCN.     

Understanding of mouse and human bladder at single‐cell resolution: integrated analysis of trajectory and cell‐cell interactive networks based on multiple scRNA‐seq datasets

ObjectivesTo elaborately decipher the mouse and human bladders at single‐cell levels.Materials and MethodsWe collected more than 50,000 cells from multiple datasets and created, up to date, the largest integrated bladder datasets. Pseudotime trajectory of urothelium and interstitial cells, as well as dynamic cell‐cell interactions, was investigated. Biological activity scores and different roles of signaling pathways between certain cell clusters were also identified.ResultsThe glucose score was significantly high in most urothelial cells, while the score of H3 acetylation was roughly equally distributed across all cell types. Several genes via a pseudotime pattern in mouse (Car3, Dkk2, Tnc, etc.) and human (FBLN1, S100A10, etc.) were discovered. S100A6, TMSB4X, and typical uroplakin genes seemed as shared pseudotime genes for urothelial cells in both human and mouse datasets. In combinational mouse (n = 16,688) and human (n = 22,080) bladders, we verified 1,330 and 1,449 interactive ligand‐receptor pairs, respectively. The distinct incoming and outgoing signaling was significantly associated with specific cell types. Collagen was the strongest signal from fibroblasts to urothelial basal cells in mouse, while laminin pathway for urothelial basal cells to smooth muscle cells (SMCs) in human. Fibronectin 1 pathway was intensely sent by myofibroblasts, received by urothelial cells, and almost exclusively mediated by SMCs in mouse bladder. Interestingly, the cell cluster of SMCs 2 was the dominant sender and mediator for Notch signaling in the human bladder, while SMCs 1 was not. The expression of integrin superfamily (the most common communicative pairs) was depicted, and their co‐expression patterns were located in certain cell types (eg, Itgb1 and Itgb4 in mouse and human basal cells).ConclusionsThis study provides a complete interpretation of the normal bladder at single‐cell levels, offering an in‐depth resource and foundation for future research.     

Nonproductive exposure of PBMCs to SARS‐CoV‐2 induces cell‐intrinsic innate immune responses

Cell‐intrinsic responses mounted in PBMCs during mild and severe COVID‐19 differ quantitatively and qualitatively. Whether they are triggered by signals emitted by productively infected cells of the respiratory tract or result from physical interaction with virus particles remains unclear. Here, we analyzed susceptibility and expression profiles of PBMCs from healthy donors upon ex vivo exposure to SARS‐CoV and SARS‐CoV‐2. In line with the absence of detectable ACE2 receptor expression, human PBMCs were refractory to productive infection. RT–PCR experiments and single‐cell RNA sequencing revealed JAK/STAT‐dependent induction of interferon‐stimulated genes (ISGs) but not proinflammatory cytokines. This SARS‐CoV‐2‐specific response was most pronounced in monocytes. SARS‐CoV‐2‐RNA‐positive monocytes displayed a lower ISG signature as compared to bystander cells of the identical culture. This suggests a preferential invasion of cells with a low ISG baseline profile or delivery of a SARS‐CoV‐2‐specific sensing antagonist upon efficient particle internalization. Together, nonproductive physical interaction of PBMCs with SARS‐CoV‐2‐ and, to a much lesser extent, SARS‐CoV particles stimulate JAK/STAT‐dependent, monocyte‐accentuated innate immune responses that resemble those detected in vivo in patients with mild COVID‐19.