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.     

Cellular and molecular atlas of the placenta from a COVID‐19 pregnant woman infected at midgestation highlights the defective impacts on foetal health

ObjectivesThe impacts of the current COVID‐19 pandemic on maternal and foetal health are enormous and of serious concern. However, the influence of SARS‐CoV‐2 infection at early‐to‐mid gestation on maternal and foetal health remains unclear.Materials and methodsHere, we report the follow‐up study of a pregnant woman of her whole infective course of SARS‐CoV‐2, from asymptomatic infection at gestational week 20 to mild and then severe illness state, and finally cured at Week 24. Following caesarean section due to incomplete uterine rupture at Week 28, histological examinations on the placenta and foetal tissues as well as single‐cell RNA sequencing (scRNA‐seq) for the placenta were performed.ResultsCompared with the gestational age‐matched control placentas, the placenta from this COVID‐19 case exhibited more syncytial knots and lowered expression of syncytiotrophoblast‐related genes. The scRNA‐seq analysis demonstrated impaired trophoblast differentiation, activation of antiviral and inflammatory CD8 T cells, as well as the tight association of increased inflammatory responses in the placenta with complement over‐activation in macrophages. In addition, levels of several inflammatory factors increased in the placenta and foetal blood.ConclusionThese findings illustrate a systematic cellular and molecular signature of placental insufficiency and immune activation at the maternal–foetal interface that may be attributed to SARS‐CoV‐2 infection at the midgestation stage, which highly suggests the extensive care for maternal and foetal outcomes in pregnant women suffering from COVID‐19.     

Neural G0: a quiescent‐like state found in neuroepithelial‐derived cells and glioma

Single‐cell RNA sequencing has emerged as a powerful tool for resolving cellular states associated with normal and maligned developmental processes. Here, we used scRNA‐seq to examine the cell cycle states of expanding human neural stem cells (hNSCs). From these data, we constructed a cell cycle classifier that identifies traditional cell cycle phases and a putative quiescent‐like state in neuroepithelial‐derived cell types during mammalian neurogenesis and in gliomas. The Neural G0 markers are enriched with quiescent NSC genes and other neurodevelopmental markers found in non‐dividing neural progenitors. Putative glioblastoma stem‐like cells were significantly enriched in the Neural G0 cell population. Neural G0 cell populations and gene expression are significantly associated with less aggressive tumors and extended patient survival for gliomas. Genetic screens to identify modulators of Neural G0 revealed that knockout of genes associated with the Hippo/Yap and p53 pathways diminished Neural G0 in vitro, resulting in faster G1 transit, down‐regulation of quiescence‐associated markers, and loss of Neural G0 gene expression. Thus, Neural G0 represents a dynamic quiescent‐like state found in neuroepithelial‐derived cells and gliomas.     

Single‐cell transcriptional profiling of human carotid plaques reveals a subpopulation of endothelial cells associated with stroke incidences

The differences in plaque histology between symptomatic and asymptomatic patients have been widely accepted. Whether there is a heterogeneity of cells between symptomatic and asymptomatic plaques remains largely unclear. To reveal the potential heterogeneity within different plaques, which may contribute to different stroke incidences, we obtained the scRNA‐seq data from symptomatic and asymptomatic patients and identified eight cell types present in plaques. Further analysis of endothelial cells (ECs) revealed three distinct EC subpopulations appeared to be endowed with specific biological functions such as antigen processing and presentation, cell adhesion, and smooth muscle cell proliferation. Of note, the differentially expressed genes of the EC 2 subpopulation showed that the genes involved in cell adhesion were up‐regulated in asymptomatic plaques compared to symptomatic plaques. Integrating the data of intraplaque haemorrhage and plaque stability, the 5th top‐enriched biological process was cell adhesion in the stable or non‐haemorrhaged plaques compared to unstable plaques or haemorrhaged plaques. Among these cell adhesion‐related genes, the intersection gene AOC3 may play a vital role in plaque haemorrhage and plaque stability. Targeting cell adhesion and the specialized genes may provide potential new therapeutic directions to prevent asymptomatic patients from stroke.     

Cytochrome P450 26A1 regulates the clusters and killing activity of NK cells during the peri‐implantation period

Cytochrome P450 26A1 (CYP26A1) plays a vital role in early pregnancy in mice. Our previous studies have found that CYP26A1 affects embryo implantation by modulating natural killer (NK) cells, and that there is a novel population of CYP26A1⁺ NK cells in the uteri of pregnant mice. The aim of this study was to investigate the effects of CYP26A1 on the subsets and killing activity of NK cells. Through single‐cell RNA sequencing (scRNA‐seq), we identified four NK cell subsets in the uterus, namely, conventional NK (cNK), tissue‐resident NK (trNK) 1 and 2, and proliferating trNK (trNKp). The two most variable subpopulations after uterine knockdown of CYP26A1 were trNKp and trNK2 cells. CYP26A1 knockdown significantly downregulated the expression of the NK cell function‐related genes Cd44, Cd160, Vegfc, and Slamf6 in trNK2 cells, and Klra17 and Ogn in trNKp cells. Both RNA‐seq and cytotoxicity assays confirmed that CYP26A1⁺ NK cells had low cytotoxicity. These results indicate that CYP26A1 may affect the immune microenvironment at the maternal‐foetal interface by regulating the activity of NK cells.     

Analyses of erythropoiesis from embryonic stem cell‐CD34+ and cord blood‐CD34+ cells reveal mechanisms for defective expansion and enucleation of embryomic stem cell‐erythroid cells

Red blood cells (RBCs) generated ex vivo have the potential to be used for transfusion. Human embryonic stem cells (ES) and induced pluripotent stem cells (iPS) possess unlimited self‐renewal capacity and are the preferred cell sources to be used for ex vivo RBC generation. However, their applications are hindered by the facts that the expansion of ES/iPS‐derived erythroid cells is limited and the enucleation of ES/iPS‐derived erythroblasts is low compared to that derived from cord blood (CB) or peripheral blood (PB). To address this, we sought to investigate the underlying mechanisms by comparing the in vitro erythropoiesis profiles of CB CD34⁺ and ES CD34⁺ cells. We found that the limited expansion of ES CD34⁺ cell‐derived erythroid cells was associated with defective cell cycle of erythroid progenitors. In exploring the cellular and molecular mechanisms for the impaired enucleation of ES CD34⁺ cell‐derived orthochromatic erythroblasts (ES‐ortho), we found the chromatin of ES‐ortho was less condensed than that of CB CD34⁺ cell‐derived orthochromatic erythroblasts (CB‐ortho). At the molecular level, both RNA‐seq and ATAC‐seq analyses revealed that pathways involved in chromatin modification were down‐regulated in ES‐ortho. Additionally, the expression levels of molecules known to play important role in chromatin condensation or/and enucleation were significantly lower in ES‐ortho compared to that in CB‐ortho. Together, our findings have uncovered mechanisms for the limited expansion and impaired enucleation of ES CD34⁺ cell‐derived erythroid cells and may help to improve ex vivo RBC production from stem cells.     

Senescence‐associated β‐galactosidase reveals the abundance of senescent CD8+ T cells in aging humans

Aging leads to a progressive functional decline of the immune system, rendering the elderly increasingly susceptible to disease and infection. The degree to which immune cell senescence contributes to this decline remains unclear, however, since markers that label immune cells with classical features of cellular senescence accurately and comprehensively have not been identified. Using a second‐generation fluorogenic substrate for β‐galactosidase and multi‐parameter flow cytometry, we demonstrate here that peripheral blood mononuclear cells (PBMCs) isolated from healthy humans increasingly display cells with high senescence‐associated β‐galactosidase (SA‐βGal) activity with advancing donor age. The greatest age‐associated increases were observed in CD8+ T‐cell populations, in which the fraction of cells with high SA‐βGal activity reached average levels of 64% in donors in their 60s. CD8+ T cells with high SA‐βGal activity, but not those with low SA‐βGal activity, were found to exhibit features of telomere dysfunction‐induced senescence and p16‐mediated senescence, were impaired in their ability to proliferate, developed in various T‐cell differentiation states, and had a gene expression signature consistent with the senescence state previously observed in human fibroblasts. Based on these results, we propose that senescent CD8+ T cells with classical features of cellular senescence accumulate to levels that are significantly higher than previously reported and additionally provide a simple yet robust method for the isolation and characterization of senescent CD8+ T cells with predictive potential for biological age.     

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.