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.     

Altered sperm tsRNAs in aged male contribute to anxiety‐like behavior in offspring

Parental age at first pregnancy is increasing worldwide. The offspring of aged father has been associated with higher risk of several neuropsychiatric disorders, such as schizophrenia and autism, but the underlying mechanism remains elusive. Here we report that advanced paternal age in mice alters the profile of transfer RNA‐derived small RNAs (tsRNAs). Injection of sperm tsRNAs from aged male mice into zygotes induced anxiety‐like behaviors in F1 males. RNA sequencing of the cerebral cortex and hippocampus of those F1 male mice altered the gene expression of dopaminergic synapse and neurotrophin. tsRNAs from aged male mice injection also altered the neuropsychiatry‐related gene expression in two‐cell and blastocyst stage embryos. More importantly, the sperm tsRNA profile changes significantly during aging in human. The up‐regulated sperm tsRNA target genes were involved in neurogenesis and nervous system development. These results suggest that aging‐related changes of sperm tsRNA may contribute to the intergenerational transmission of behavioral traits.     

B‐cell capacity for differentiation changes with age

BackgroundAge‐related immune deficiencies are thought to be responsible for increased susceptibility to infection in older adults, with alterations in lymphocyte populations becoming more prevalent over time. The loss of humoral immunity in ageing was attributed to the diminished numbers of B cells and the reduced ability to generate immunoglobulin.AimsTo compare the intrinsic B‐cell capacity for differentiation into mature plasma cells (PCs), between young and old donors, using in vitro assays, providing either effective T‐cell help or activation via TLR engagement.MethodsB cells were isolated from healthy individuals, in younger (30–38 years) and older (60–64 years) donors. An in vitro model system of B‐cell differentiation was used, analysing 5 differentiation markers by flow cytometry, under T‐dependent (TD: CD40/BCR stimulation) or T‐independent (TI: TLR7/BCR activation) conditions. Antibody secretion was measured by ELISA and gene expression using qPCR.ResultsTI and TD differentiation resulted in effective proliferation of B cells followed by their differentiation into PC. B‐cell‐executed TI differentiation was faster, all differentiation marker and genes being expressed earlier than under TD differentiation (day 6), although generating less viable cells and lower antibody levels (day 13). Age‐related differences in B‐cell capacity for differentiation were minimal in TD differentiation. In contrast, in TI differentiation age significantly affected proliferation, viability, differentiation, antibody secretion and gene expression, older donors being more efficient.ConclusionAltogether, B‐cell differentiation into PC appeared similar between age groups when provided with T‐cell help, in contrast to TI differentiation, where multiple age‐related changes suggest better capacities in older donors. These new findings may help explain the emergence of autoantibodies in ageing.     

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.     

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.