Co-expression of immune checkpoints in glioblastoma revealed by single-nucleus RNA sequencing and spatial transcriptomics

Glioblastoma (GBM) is one of the most malignant tumors of the central nervous system. The pattern of immune checkpoint expression in GBM remains largely unknown. We performed snRNA-Seq and spatial transcriptomic (ST) analyses on untreated GBM samples. 8 major cell types were found in both tumor and adjacent normal tissues, with variations in infiltration grade. Neoplastic cells_6 was identified in malignant cells with high expression of invasion and proliferator-related genes, and analyzed its interactions with microglia, MDM cells and T cells. Significant alterations in ligand-receptor interactions were observed, particularly between Neoplastic cells_6 and microglia, and found prominent expression of VISTA/VSIG3, suggesting a potential mechanism for evading immune system attacks. High expression of TIM-3, VISTA, PSGL-1 and VSIG-3 with similar expression patterns in GBM, may have potential as therapeutic targets. The prognostic value of VISTA expression was cross-validated in 180 glioma patients, and it was observed that patients with high VISTA expression had a poorer prognosis. In addition, multimodal cross analysis integrated SnRNA-seq and ST, revealing complex intracellular communication and mapping the GBM tumor microenvironment. This study reveals novel molecular characteristics of GBM, co-expression of immune checkpoints, and potential therapeutic targets, contributing to improving the understanding and treatment of GBM.

Single-nucleus RNA and multiomics in situ pairwise sequencing reveals cellular heterogeneity of the abnormal ligamentum teres in patients with developmental dysplasia of the hip

Developmental dysplasia of the hip (DDH) is the most common hip deformity in pediatric orthopedics. One of the common pathological changes in DDH is the thickening and hypertrophy of the ligamentum teres. However, the underlying pathogenic mechanism responsible for these changes remains unclear. This study represents the first time that the heterogeneity of cell subsets in the abnormal ligamentum teres of patients with DDH has been resolved at the single-cell and spatial levels by snRNA-Seq and MiP-Seq. Through gene set enrichment and intercellular communication network analyses, we found that receptor-like cells and ligament stem cells may play an essential role in the pathological changes resulting in ligamentum teres thickening and hypertrophy. Eight ligand-receptor pairs related to the ECM-receptor pathway were observed to be closely associated with DDH. Further, using the Monocle R package, we predicted a differentiation trajectory of pericytes into two branches, leading to junctional ligament stem cells or fibroblasts. The expression of extracellular matrix-related genes along pseudotemporal trajectories was also investigated. Using MiP-Seq, we determined the expression distribution of marker genes specific to different cell types within the ligamentum teres, as well as differentially expressed DDH-associated genes at the spatial level.

MEF2C Hypofunction in GABAergic Cells Alters Sociability and Prefrontal Cortex Inhibitory Synaptic Transmission in a Sex-Dependent Manner

Background

Heterozygous mutations or deletions of MEF2C cause a neurodevelopmental disorder termed MEF2C haploinsufficiency syndrome (MCHS), characterized by autism spectrum disorder and neurological symptoms. In mice, global Mef2c heterozygosity has produced multiple MCHS-like phenotypes. MEF2C is highly expressed in multiple cell types of the developing brain, including GABAergic (gamma-aminobutyric acidergic) inhibitory neurons, but the influence of MEF2C hypofunction in GABAergic neurons on MCHS-like phenotypes remains unclear.

Methods

We employed GABAergic cell type-specific manipulations to study mouse Mef2c heterozygosity in a battery of MCHS-like behaviors. We also performed electroencephalography, single-cell transcriptomics, and patch-clamp electrophysiology and optogenetics to assess the impact of Mef2c haploinsufficiency on gene expression and prefrontal cortex microcircuits.

Results

Mef2c heterozygosity in developing GABAergic cells produced female-specific deficits in social preference and altered approach-avoidance behavior. In female, but not male, mice, we observed that Mef2c heterozygosity in developing GABAergic cells produced 1) differentially expressed genes in multiple cell types, including parvalbumin-expressing GABAergic neurons, 2) baseline and social-related frontocortical network activity alterations, and 3) reductions in parvalbumin cell intrinsic excitability and inhibitory synaptic transmission onto deep-layer pyramidal neurons.

Conclusions

MEF2C hypofunction in female, but not male, developing GABAergic cells is important for typical sociability and approach-avoidance behaviors and normal parvalbumin inhibitory neuron function in the prefrontal cortex of mice. While there is no apparent sex bias in autism spectrum disorder symptoms of MCHS, our findings suggest that GABAergic cell-specific dysfunction in females with MCHS may contribute disproportionately to sociability symptoms.

Single-Nucleus Sequencing of Silkworm Larval Brain Reveals the Key Role of Lysozyme in the Antiviral Immune Response in Brain Hemocytes

INTRODUCTION

The brain is considered as an immune-privileged organ, yet innate immune reactions can occur in the central nervous system of vertebrates and invertebrates. Silkworm (Bombyx mori) is an economically important insect and a lepidopteran model species. The diversity of cell types in the silkworm brain, and how these cell subsets produce an immune response to virus infection, remains largely unknown.

METHODS

Single-nucleus RNA sequencing (snRNA-seq), bioinformatics analysis, RNAi, and other methods were mainly used to analyze the cell types and gene functions of the silkworm brain.

RESULTS

We used snRNA-seq to identify 19 distinct clusters representing Kenyon cell, glial cell, olfactory projection neuron, optic lobes neuron, hemocyte-like cell, and muscle cell types in the B. mori nucleopolyhedrovirus (BmNPV)-infected and BmNPV-uninfected silkworm larvae brain at the late stage of infection. Further, we found that the cell subset that exerts an antiviral function in the silkworm larvae brain corresponds to hemocytes. Specifically, antimicrobial peptides were significantly induced by BmNPV infection in the hemocytes, especially lysozyme, exerting antiviral effects.

CONCLUSION

Our single-cell dataset reveals the diversity of silkworm larvae brain cells, and the transcriptome analysis provides insights into the immune response following virus infection at the single-cell level.

Integration of single-nucleus RNA sequencing and network disturbance to elucidate crosstalk between multicomponent drugs and trigeminal ganglia cells in migraine

ETHNOPHARMACOLOGICAL RELEVANCE

Migraine is caused by hyperactivity of the trigeminovascular system, where trigeminal ganglia (TG) plays an important role. TG is composed of multiple neuronal and non-neuronal cell types, which is related to "neuro-inflammation-vascular" disorder in migraine. Tou Tong Ning capsule (TTNC), a CFDA-approved traditional Chinese medicine for treating migraine, has the characteristics of "multicomponents, multitargets, multipathways".

AIM OF THE STUDY

To clarify the mechanism of TTNC and elucidate crosstalk between multicomponent drugs and neuronal and non-neuronal functions and cells in migraine.

MATERIALS AND METHODS

We integrated single-nucleus RNA sequencing and a quantitative evaluation algorithm of the disturbance of multitarget drugs on the disease network and explored the specific pathology of migraine and corresponding compounds. A cerebrovascular smooth muscle spasmolytic activity experiment was carried out to verify the results of the bioinformatics analysis.

RESULTS

TTNC exhibited its regulation activities in neuronal and non-neuronal aspects based on drugs attack to four subnetworks and cell specific networks, which explored the MoA of TTNC in comprehensive and refined perspectives. Compared to neuronal regulation, TTNC showed more significant attack score on non-neuronal biological function (smooth muscle and vessel). And TTNC compound clusters C1, C6 and C7, targeting non-neuronal function and cells, had larger group area than C10, C4 and C6 for neuronal function and cell, which implied that TTNC may mainly regulate the non-neuronal function, e.g., vessel smooth muscle contraction. Contraction of cerebrovascular smooth muscle of mice ex vivo confirmed the vasodilation activity of TTNC and active compounds from C1, C6, C9 (Emodin, Luteolin and Levistilide A). Literature mining confirmed the vasospasmodolytic activity and neuroprotective effect of TTNC.

CONCLUSIONS

The study found that TTNC may primarily alleviate non-neuronal functional disorders in migraine by relaxing cerebral vascular smooth muscle cell spasm to alleviate migraine. Integrating single-nucleus RNA sequencing data and network disturbance tools provides a new strategy for the pharmacological mechanism of multicomponent drugs through cell subtyping.

Unraveling the intercellular communication disruption and key pathways in Alzheimer's disease: an integrative study of single-nucleus transcriptomes and genetic association

BACKGROUND

Recently, single-nucleus RNA-seq (snRNA-seq) analyses have revealed important cellular and functional features of Alzheimer's disease (AD), a prevalent neurodegenerative disease. However, our knowledge regarding intercellular communication mediated by dysregulated ligand-receptor (LR) interactions remains very limited in AD brains.

METHODS

We systematically assessed the intercellular communication networks by using a discovery snRNA-seq dataset comprising 69,499 nuclei from 48 human postmortem prefrontal cortex (PFC) samples. We replicated the findings using an independent snRNA-seq dataset of 56,440 nuclei from 18 PFC samples. By integrating genetic signals from AD genome-wide association studies (GWAS) summary statistics and whole genome sequencing (WGS) data, we prioritized AD-associated Gene Ontology (GO) terms containing dysregulated LR interactions. We further explored drug repurposing for the prioritized LR pairs using the Therapeutic Targets Database.

RESULTS

We identified 190 dysregulated LR interactions across six major cell types in AD PFC, of which 107 pairs were replicated. Among the replicated LR signals, we found globally downregulated communications in the astrocytes-to-neurons signaling axis, characterized, for instance, by the downregulation of APOE-related and Calmodulin (CALM)-related LR interactions and their potential regulatory connections to target genes. Pathway analyses revealed 44 GO terms significantly linked to AD, highlighting Biological Processes such as 'amyloid precursor protein processing' and 'ion transmembrane transport,' among others. We prioritized several drug repurposing candidates, such as cromoglicate, targeting the identified dysregulated LR pairs.

CONCLUSIONS

Our integrative analysis identified key dysregulated LR interactions in a cell type-specific manner and the associated GO terms in AD, offering novel insights into potential therapeutic targets involved in disrupted cell-cell communication in AD.

Single-nucleus sequencing of silkworm larval midgut reveals the immune escape strategy of BmNPV in the midgut during the late stage of infection

The midgut is an important barrier against microorganism invasion and proliferation, yet is the first tissue encountered when a baculovirus naturally invades the host. However, only limited knowledge is available how different midgut cell types contribute to the immune response and the clearance or promotion of viral infection. Here, single-nucleus RNA sequencing (snRNA seq) was employed to analyze the responses of various cell subpopulations in the silkworm larval midgut to B. mori nucleopolyhedrovirus (BmNPV) infection. We identified 22 distinct clusters representing enteroendocrine cells (EEs), enterocytes (ECs), intestinal stem cells (ISCs), Goblet cell-like and muscle cell types in the BmNPV-infected and uninfected silkworm larvae midgut at 72 h post infection. Further, our results revealed that the strategies for immune escape of BmNPV in the midgut at the late stage of infection include (1) inhibiting the response of antiviral pathways; (2) inhibiting the expression of antiviral host factors; (3) stimulating expression levels of genes promoting BmNPV replication. These findings suggest that the midgut, as the first line of defense against the invasion of the baculovirus, has dual characteristics of "resistance" and "tolerance". Our single-cell dataset reveals the diversity of silkworm larval midgut cells, and the transcriptome analysis provides insights into the interaction between host and virus infection at the single-cell level.

Single-nucleus transcriptome sequencing reveals hepatic cell atlas in pigs

BACKGROUND

As the largest substantive organ of animals, the liver plays an essential role in the physiological processes of digestive metabolism and immune defense. However, the cellular composition of the pig liver remains poorly understood. This investigation used single-nucleus RNA sequencing technology to identify cell types from liver tissues of pigs, providing a theoretical basis for further investigating liver cell types in pigs.

RESULTS

The analysis revealed 13 cells clusters which were further identified 7 cell types including endothelial cells, T cells, hepatocytes, Kupffer cells, stellate cells, B cells, and cholangiocytes. The dominant cell types were endothelial cells, T cells and hepatocytes in the liver tissue of Dahe pigs and Dahe black pigs, which accounts for about 85.76% and 82.74%, respectively. The number of endothelial cells was higher in the liver tissue of Dahe pigs compared to Dahe black pigs, while the opposite tendency was observed for T cells. Moreover, functional enrichment analysis demonstrated that the differentially expressed genes in pig hepatic endothelial cells were significantly enriched in the protein processing in endoplasmic reticulum, MAPK signaling pathway, and FoxO signaling pathway. Functional enrichment analysis demonstrated that the differentially expressed genes in pig hepatic T cells were significantly enriched in the thyroid hormone signaling pathway, B cell receptor signaling pathway, and focal adhesion. Functional enrichment analysis demonstrated that the differentially expressed genes in pig hepatic hepatocytes were significantly enriched in the metabolic pathways.

CONCLUSIONS

In summary, this study provides a comprehensive cell atlas of porcine hepatic tissue. The number, gene expression level and functional characteristics of each cell type in pig liver tissue varied between breeds.

Neuronal SNCA transcription during Lewy body formation

Misfolded α-synuclein (α-syn) is believed to contribute to neurodegeneration in Lewy body disease (LBD) based on considerable evidence including a gene-dosage effect observed in relation to point mutations and multiplication of SNCA in familial Parkinson's disease. A contradictory concept proposes early loss of the physiological α-syn as the major driver of neurodegeneration. There is a paucity of data on SNCA transcripts in various α-syn immunoreactive cytopathologies. Here, the total cell body, nuclear, and cytoplasmic area density of SNCA transcripts in neurons without and with various α-syn immunoreactive cytopathologies in the substantia nigra and amygdala in autopsy cases of LBD (n = 5) were evaluated using RNAscope combined with immunofluorescence for disease-associated α-syn. Single-nucleus RNA sequencing was performed to elucidate cell-type specific SNCA expression in non-diseased frontal cortex (n = 3). SNCA transcripts were observed in the neuronal nucleus and cytoplasm in neurons without α-syn, those containing punctate α-syn immunoreactivity, irregular-shaped compact inclusion, and brainstem-type and cortical-type LBs. However, SNCA transcripts were only rarely found in the α-syn immunoreactive LB areas. The total cell body SNCA transcript area densities in neurons with punctate α-syn immunoreactivity were preserved but were significantly reduced in neurons with compact α-syn inclusions both in the substantia nigra and amygdala. This reduction was also observed in the cytoplasm but not in the nucleus. Only single SNCA transcripts were detected in astrocytes with or without disease-associated α-syn immunoreactivity in the amygdala. Single-nucleus RNA sequencing revealed that excitatory and inhibitory neurons, oligodendrocyte progenitor cells, oligodendrocytes, and homeostatic microglia expressed SNCA transcripts, while expression was largely absent in astrocytes and microglia. The preserved cellular SNCA expression in the more abundant non-Lewy body type α-syn cytopathologies might provide a pool for local protein production that can aggregate and serve as a seed for misfolded α-syn. Successful segregation of disease-associated α-syn is associated with the exhaustion of SNCA production in the terminal cytopathology, the Lewy body. Our observations inform therapy development focusing on targeting SNCA transcription in LBD.

Novel insights into molecular signatures and pathogenic cell populations shared by systemic lupus erythematosus and vascular dementia

Although vascular dementia (VD) and systemic lupus erythematosus (SLE) may share immune-mediated pathophysiologic processes, the underlying mechanisms are unclear. This study investigated shared gene signatures in SLE versus VD, as well as their potential molecular mechanisms. Bulk RNA sequencing (RNAseq) and single-cell or single-nucleus RNAseq (sc/snRNAseq) datasets from SLE blood samples and VD brain samples were obtained from Gene Expression Omnibus. The identification of genes associated with both SLE and VD was performed using the weighted gene co-expression network analysis (WGCNA) and machine learning algorithms. For the sc/snRNAseq data, an unbiased clustering pipeline based on Seurat and CellChat was used to determine the cellular landscape profile and examine intracellular communication, respectively. The results were subsequently validated using a mice model of SLE with cognitive dysfunction (female MRL/lpr mice). WGCNA and machine learning identified C1QA, LY96, CD163, and MS4A4A as key genes for SLE and VD. sc/snRNAseq analyses revealed that CD163 and MS4A4A were upregulated in mononuclear phagocytes (MPs) from SLE and VD samples and were associated with monocyte-macrophage differentiation. Intriguingly, LGALS9-associated molecular pathway, as the only signaling pathway common between SLE and VD via CellChat analysis, exhibited significant upregulation in cortical microglia of MRL/lpr mice. Our analyses identified C1QA, LY96, CD163, and MS4A4A as potential biomarkers for SLE and VD. Moreover, the upregulation of CD163/MS4A4A and activation of LGALS9 signaling in MPs may contribute to the pathogenesis of VD with SLE. These findings offer novel insight into the mechanisms underlying VD in SLE patients.

Uncovering the prominent role of satellite cells in paravertebral muscle development and aging by single-nucleus RNA sequencing

To uncover the role of satellite cells (SCs) in paravertebral muscle development and aging, we constructed a single-nucleus transcriptomic atlas of mouse paravertebral muscle across seven timepoints spanning the embryo (day 16.5) to old (month 24) stages. Eight cell types, including SCs, fast muscle cells, and slow muscle cells, were identified. An energy metabolism-related gene set, TCA CYCLE IN SENESCENCE, was enriched in SCs. Forty-two skeletal muscle disease-related genes were highly expressed in SCs and exhibited similar expression patterns. Among them, Pdha1 was the core gene in the TCA CYCLE IN SENESCENCE; Pgam2, Sod1, and Suclg1 are transcription factors closely associated with skeletal muscle energy metabolism. Transcription factor enrichment analysis of the 42 genes revealed that Myod1 and Mef2a were also highly expressed in SCs, which regulated Pdha1 expression and were associated with skeletal muscle development. These findings hint that energy metabolism may be pivotal in SCs development and aging. Three ligand-receptor pairs of extracellular matrix (ECM)-receptor interactions, Lamc1-Dag1, Lama2-Dag1, and Hspg2-Dag1, may play a vital role in SCs interactions with slow/fast muscle cells and SCs self-renewal. Finally, we built the first database of a skeletal muscle single-cell transcriptome, the Musculoskeletal Cell Atlas (http://www.mskca.tech), which lists 630,040 skeletal muscle cells and provides interactive visualization, a useful resource for revealing skeletal muscle cellular heterogeneity during development and aging. Our study could provide new targets and ideas for developing drugs to inhibit skeletal muscle aging and treat skeletal muscle diseases.

Unraveling the intercellular communication disruption and key pathways in Alzheimer's disease: An integrative study of single-nucleus transcriptomes and genetic association

Background

Recently, single-nucleus RNA-seq (snRNA-seq) analyses have revealed important cellular and functional features of Alzheimer's disease (AD), a prevalent neurodegenerative disease. However, our knowledge regarding intercellular communication mediated by dysregulated ligand-receptor (LR) interactions remains very limited in AD brains.

Methods

We systematically assessed the intercellular communication networks by using a discovery snRNA-seq dataset comprising 69,499 nuclei from 48 human postmortem prefrontal cortex (PFC) samples. We replicated the findings using an independent snRNA-seq dataset of 56,440 nuclei from 18 PFC samples. By integrating genetic signals from AD genome-wide association studies (GWAS) summary statistics and whole genome sequencing (WGS) data, we prioritized AD-associated Gene Ontology (GO) terms containing dysregulated LR interactions. We further explored drug repurposing for the prioritized LR pairs using the Therapeutic Targets Database.

Results

We identified 316 dysregulated LR interactions across six major cell types in AD PFC, of which 210 pairs were replicated. Among the replicated LR signals, we found globally downregulated communications in astrocytes-to-neurons signaling axis, characterized, for instance, by the downregulation of APOE-related and Calmodulin (CALM)-related LR interactions and their potential regulatory connections to target genes. Pathway analyses revealed 60 GO terms significantly linked to AD, highlighting Biological Processes such as 'amyloid precursor protein processing' and 'ion transmembrane transport', among others. We prioritized several drug repurposing candidates, such as cromoglicate, targeting the identified dysregulated LR pairs.

Conclusions

Our integrative analysis identified key dysregulated LR interactions in a cell type-specific manner and the associated GO terms in AD, offering novel insights into potential therapeutic targets involved in disrupted cell-cell communication in AD.

Updated skin transcriptomic atlas depicted by reciprocal contribution of single-nucleus RNA sequencing and single-cell RNA sequencing

BACKGROUND

Single-cell RNA sequencing (scRNA-seq) has advanced our understanding of skin biology, but its utility is restricted by the requirement of fresh samples, inadequate dissociation-induced cell loss or death, and activation during tissue digestion. Single-nucleus RNA sequencing (snRNA-seq) can use frozen, hard-to-dissociate materials, which might be a promising method to circumvent the limitations of scRNA-seq for the skin tissue.

OBJECTIVE

To profile skin cells using snRNA-seq in parallel with scRNA-seq.

METHODS

We performed snRNA-seq in parallel with scRNA-seq for the bisected skin sample of one person and integrated previously published scRNA-seq data for analysis. We comparatively analyzed the differences in cell proportions and gene expression between the two methods. The differentiation trajectories of keratinocytes and fibroblasts were analyzed by Slingshot analysis.

RESULTS

snRNA-seq was less susceptible to contamination from mitochondrial and ribosomal RNA, and exhibited a greater capacity to detect transcription factors. snRNA-seq identified more spatially and functionally relevant keratinocyte clusters that constitute cell trajectories with expected differentiation dynamics. Novel markers, e.g., LYPD3, EMP2, and CSTB, were revealed for different differentiation stages of keratinocytes, and NFIB and GRHL1 were identified as transcription factors involving in the proliferation and functional differentiation of keratinocytes. Fibroblasts were found in a state of activation in scRNA-seq. And scRNA-seq detected a greater number of immune cells.

CONCLUSIONS

We generated an updated atlas of the skin transcriptome based on the reciprocal contribution of scRNA-seq and snRNA-seq.

Single-nucleus and Spatially Resolved Intratumor Subtype Heterogeneity in Bladder Cancer

Background

Current bulk transcriptomic classification systems for bladder cancer do not consider the level of intratumor subtype heterogeneity.

Objective

To investigate the extent and possible clinical impact of intratumor subtype heterogeneity across early and more advanced stages of bladder cancer.

Design setting and participants

We performed single-nucleus RNA sequencing (RNA-seq) of 48 bladder tumors and additional spatial transcriptomics for four of these tumors. Total bulk RNA-seq and spatial proteomics data were available from the same tumors for comparison, along with detailed clinical follow-up of the patients.

Outcome measurements and statistical analysis

The primary outcome was progression-free survival for non-muscle-invasive bladder cancer. Cox regression analysis, log-rank tests, Wilcoxon rank-sum tests, Spearman correlation, and Pearson correlation were used for statistical analysis.

Results and limitations

We found that the tumors exhibited varying levels of intratumor subtype heterogeneity and that the level of subtype heterogeneity can be estimated from both single-nucleus and bulk RNA-seq data, with high concordance between the two. We found that a higher class 2a weight estimated from bulk RNA-seq data is associated with worse outcome for patients with molecular high-risk class 2a tumors. The sparsity of the data generated using the DroNc-seq sequencing protocol is a limitation.

Conclusions

Our results indicate that discrete subtype assignments from bulk RNA-seq data may lack biological granularity and that continuous class scores may improve clinical risk stratification of patients with bladder cancer.

Patient summary

We found that several molecular subtypes can exist within a single bladder tumor and that continuous subtype scores can be used to identify a subgroup of patients with poor outcomes. Use of these subtype scores may improve risk stratification for patients with bladder cancer, which can help in making decisions on treatment.

Single-nucleus RNA sequencing of human pancreatic islets identifies novel gene sets and distinguishes β-cell subpopulations with dynamic transcriptome profiles

BACKGROUND

Single-cell RNA sequencing (scRNA-seq) provides valuable insights into human islet cell types and their corresponding stable gene expression profiles. However, this approach requires cell dissociation that complicates its utility in vivo. On the other hand, single-nucleus RNA sequencing (snRNA-seq) has compatibility with frozen samples, elimination of dissociation-induced transcriptional stress responses, and affords enhanced information from intronic sequences that can be leveraged to identify pre-mRNA transcripts.

METHODS

We obtained nuclear preparations from fresh human islet cells and generated snRNA-seq datasets. We compared these datasets to scRNA-seq output obtained from human islet cells from the same donor. We employed snRNA-seq to obtain the transcriptomic profile of human islets engrafted in immunodeficient mice. In both analyses, we included the intronic reads in the snRNA-seq data with the GRCh38-2020-A library.

RESULTS

First, snRNA-seq analysis shows that the top four differentially and selectively expressed genes in human islet endocrine cells in vitro and in vivo are not the canonical genes but a new set of non-canonical gene markers including ZNF385D, TRPM3, LRFN2, PLUT (β-cells); PTPRT, FAP, PDK4, LOXL4 (α-cells); LRFN5, ADARB2, ERBB4, KCNT2 (δ-cells); and CACNA2D3, THSD7A, CNTNAP5, RBFOX3 (γ-cells). Second, by integrating information from scRNA-seq and snRNA-seq of human islet cells, we distinguish three β-cell sub-clusters: an INS pre-mRNA cluster (β3), an intermediate INS mRNA cluster (β2), and an INS mRNA-rich cluster (β1). These display distinct gene expression patterns representing different biological dynamic states both in vitro and in vivo. Interestingly, the INS mRNA-rich cluster (β1) becomes the predominant sub-cluster in vivo.

CONCLUSIONS

In summary, snRNA-seq and pre-mRNA analysis of human islet cells can accurately identify human islet cell populations, subpopulations, and their dynamic transcriptome profile in vivo.

Recent advances in deciphering hippocampus complexity using single-cell transcriptomics

Single-cell and single-nucleus RNA sequencing (scRNA-seq and snRNA-seq) technologies have emerged as revolutionary and powerful tools, which have helped in achieving significant progress in biomedical research over the last decade. scRNA-seq and snRNA-seq resolve heterogeneous cell populations from different tissues and help reveal the function and dynamics at the single-cell level. The hippocampus is an essential component for cognitive functions, including learning, memory, and emotion regulation. However, the molecular mechanisms underlying the activity of hippocampus have not been fully elucidated. The development of scRNA-seq and snRNA-seq technologies provides strong support for attaining an in-depth understanding of hippocampal cell types and gene expression regulation from the single-cell transcriptome profiling perspective. This review summarizes the applications of scRNA-seq and snRNA-seq in the hippocampus to further expand our knowledge of the molecular mechanisms related to hippocampal development, health, and diseases.

Single-nucleus RNA sequencing and deep tissue proteomics reveal distinct tumour microenvironment in stage-I and II cervical cancer

BACKGROUND

Cervical cancer (CC) is the 3rd most common cancer in women and the 4th leading cause of deaths in gynaecological malignancies, yet the exact progression of CC is inconclusive, mainly due to the high complexity of the changing tumour microenvironment (TME) at different stages of tumorigenesis. Importantly, a detailed comparative single-nucleus transcriptomic analysis of tumour microenvironment (TME) of CC patients at different stages is lacking.

METHODS

In this study, a total of 42,928 and 29,200 nuclei isolated from the tumour tissues of stage-I and II CC patients and subjected to single-nucleus RNA sequencing (snRNA-seq) analysis. The cell heterogeneity and functions were comparatively investigated using bioinformatic tools. In addition, label-free quantitative mass spectrometry based proteomic analysis was carried out. The proteome profiles of stage-I and II CC patients were compared, and an integrative analysis with the snRNA-seq was performed.

RESULTS

Compared with the stage-I CC (CCI) patients, the immune response relevant signalling pathways were largely suppressed in various immune cells of the stage-II CC (CCII) patients, yet the signalling associated with cell and tissue development was enriched, as well as metabolism for energy production suggested by the upregulation of genes associated with mitochondria. This was consistent with the quantitative proteomic analysis that showed the dominance of proteins promoting cell growth and intercellular matrix development in the TME of CCII group. The interferon-α and γ responses appeared the most activated pathways in many cell populations of the CCI patients. Several collagens, such as COL12A1, COL5A1, COL4A1 and COL4A2, were found significantly upregulated in the CCII group, suggesting their roles in diagnosing CC progression. A novel transcript AC244205.1 was detected as the most upregulated gene in CCII patients, and its possible mechanistic role in CC may be investigated further.

CONCLUSIONS

Our study provides important resources for decoding the progression of CC and set the foundation for developing novel approaches for diagnosing CC and tackling the immunosuppressive TME.

Identification of resident progenitors labeled with Top2a responsible for proximal tubular regeneration in ischemia reperfusion-induced acute kidney injury

BACKGROUND

Acute kidney injury is a common fatal disease with complex etiology and limited treatment methods. Proximal tubules (PTs) are the most vulnerable segment. Not only in injured kidneys but also in normal kidneys, shedding of PTs often happens. However, the source cells and mechanism of their regeneration remain unclear.

METHODS

ScRNA and snRNA sequencing data of acute injured or normal kidney were downloaded from GEO database to identify the candidate biomarker of progenitor of proximal tubules. SLICE algorithm and CytoTRACE analyses were employed to evaluate the stemness of progenitors. Then the repairing trajectory was constructed through pseudotime analyses. SCENIC algorithm was used to detect cell-type-specific regulon. With spatial transcriptome data, the location of progenitors was simulated. Neonatal/ adult/ aged mice and preconditioning AKI mice model and deconvolution of 2 RNA-seq data were employed for validation.

RESULTS

Through cluster identification, PT cluster expressed Top2a specifically was identified to increase significantly during AKI. With relatively strong stemness, the Top2a-labeled PT cluster tended to be the origin of the repairing trajectory. Moreover, the cluster was regulated by Pbx3-based regulon and possessed great segmental heterogeneity. Changes of Top2a between neonatal and aged mice and among AKI models validated the progenitor role of Top2a-labeled cluster.

CONCLUSIONS

Our study provided transcriptomic evidence that resident proximal tubular progenitors labeled with Top2a participated in regeneration. Considering the segmental heterogeneity, we find that there is a group of reserve progenitor cells in each tubular segment. When AKI occurs, the reserve progenitors of each tubular segment proliferate and replenish first, and PT-progenitors, a cluster with no obvious PT markers replenish each subpopulation of the reserve cells.

Single-nucleus transcriptomic profiling of multiple organs in a rhesus macaque model of SARS-CoV-2 infection

Infection with severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) causes diverse clinical manifestations and tissue injuries in multiple organs. However, cellular and molecular understanding of SARS-CoV-2 infection-associated pathology and immune defense features in different organs remains incomplete. Here, we profiled approximately 77 000 single-nucleus transcriptomes of the lung, liver, kidney, and cerebral cortex in rhesus macaques ( Macaca mulatta) infected with SARS-CoV-2 and healthy controls. Integrated analysis of the multi-organ dataset suggested that the liver harbored the strongest global transcriptional alterations. We observed prominent impairment in lung epithelial cells, especially in AT2 and ciliated cells, and evident signs of fibrosis in fibroblasts. These lung injury characteristics are similar to those reported in patients with coronavirus disease 2019 (COVID-19). Furthermore, we found suppressed MHC class I/II molecular activity in the lung, inflammatory response in the liver, and activation of the kynurenine pathway, which induced the development of an immunosuppressive microenvironment. Analysis of the kidney dataset highlighted tropism of tubule cells to SARS-CoV-2, and we found membranous nephropathy (an autoimmune disease) caused by podocyte dysregulation. In addition, we identified the pathological states of astrocytes and oligodendrocytes in the cerebral cortex, providing molecular insights into COVID-19-related neurological implications. Overall, our multi-organ single-nucleus transcriptomic survey of SARS-CoV-2-infected rhesus macaques broadens our understanding of disease features and antiviral immune defects caused by SARS-CoV-2 infection, which may facilitate the development of therapeutic interventions for COVID-19.

Single-cell and spatially resolved analysis uncovers cell heterogeneity of breast cancer

The heterogeneity and the complex cellular architecture have a crucial effect on breast cancer progression and response to treatment. However, deciphering the neoplastic subtypes and their spatial organization is still challenging. Here, we combine single-nucleus RNA sequencing (snRNA-seq) with a microarray-based spatial transcriptomics (ST) to identify cell populations and their spatial distribution in breast cancer tissues. Malignant cells are clustered into distinct subpopulations. These cell clusters not only have diverse features, origins and functions, but also emerge to the crosstalk within subtypes. Furthermore, we find that these subclusters are mapped in distinct tissue regions, where discrepant enrichment of stromal cell types are observed. We also inferred the abundance of these tumorous subpopulations by deconvolution of large breast cancer RNA-seq cohorts, revealing differential association with patient survival and therapeutic response. Our study provides a novel insight for the cellular architecture of breast cancer and potential therapeutic strategies.

Cellular diversity and gene expression profiles in the male and female brain of Aedes aegypti

BACKGROUND

Aedes aegypti is a medically-important mosquito vector that transmits arboviruses including yellow fever, dengue, chikungunya, and Zika viruses to humans. The mosquito exhibits typical sexually dimorphic behaviors such as courtship, mating, host seeking, bloodfeeding, and oviposition. All these behaviors are mainly regulated by the brain; however, little is known about the function and neuron composition of the mosquito brain. In this study, we generated an initial atlas of the adult male and female brain of Ae. aegypti using 10xGenomics based single-nucleus RNA sequencing.

RESULTS

We identified 35 brain cell clusters in male and female brains, and 15 of those clusters were assigned to known cell types. Identified cell types include glia (astrocytes), Kenyon cells, (ventral) projection neurons, monoaminergic neurons, medulla neurons, and proximal medulla neurons. In addition, the cell type compositions of male and female brains were compared to each other showing that they were quantitatively distinct, as 17 out of 35 cell clusters varied significantly in their cell type proportions. Overall, the transcriptomes from each cell cluster looked very similar between the male and female brain as only up to 25 genes were differentially expressed in these clusters. The sex determination factor Nix was highly expressed in neurons and glia of the male brain, whereas doublesex (dsx) was expressed in all neuron and glia cell clusters of the male and female brain.

CONCLUSIONS

An initial cell atlas of the brain of the mosquito Ae. aegypti has been generated showing that the cellular compositions of the male and female brains of this hematophagous insect differ significantly from each other. Although some of the rare brain cell types have not been detected in our single biological replicate, this study provides an important basis for the further development of a complete brain cell atlas as well as a better understanding of the neurobiology of the brains of male and female mosquitoes and their sexually dimorphic behaviors.

Hippocampal neuropathology in suicide: Gaps in our knowledge and opportunities for a breakthrough

Suicide is a major global hazard. There is a need for increasing suicide awareness and effective and evidence-based interventions, targeting both suicidal ideation and conduct. However, anti-suicide pharmacological effects are unsatisfactory. The human hippocampus is vulnerable to neuropsychiatric damages and subsequently releases psychobiological signals. Human hippocampal studies of suicide completers have shown mechanistic changes in neurobiology, which, however, could not reflect the neuropathological 'fingerprints' of fatal suicide ideations and suicide attempts. In this review, we provide several leading theories of suicide, including the serotoninergic system, Wnt pathway and brain-derived neurotrophic factor/tropomyosin receptor kinase B signalling, and discuss the evidence for their roles in suicide and treatment. Moreover, the cognitive dysfunctions associated with suicide risk are discussed, as well as the novel evidence on cognitive therapies that decrease suicidal ideation. We highlight the need to apply multi-omics techniques (including single-nucleus RNA sequencing and mass spectrometry histochemistry) on hippocampal samples from donors who died by suicide or legal euthanasia, to clarify the aetiology of suicide and propose novel therapeutic strategies.

Single-Cell RNAseq Analysis of lncRNAs

Mammalian genomes are pervasively transcribed and a small fraction of RNAs produced codify for proteins. The importance of noncoding RNAs for the maintenance of cell functions is well known (e.g., rRNAs, tRNAs), but only recently it was first demonstrated the involvement of microRNAs (miRNAs) in posttranscriptional regulation and then the activity of long noncoding RNAs (lncRNAs) in the regulation of miRNAs, DNA structure and protein function. LncRNAs have an expression more cell specific than other RNAs and basing on their subcellular localization exert different functions. In this book chapter we consider different protocols to evaluate the expression of lncRNAs at the single cell level using genome-wide approaches. We considered the skeletal muscle as example because the most abundant tissue in mammals involved in the regulation of metabolism and body movement. We firstly described how to isolate the smallest complete contractile system responsible for muscle metabolic and contractile traits (myofibers). We considered how to separate long and short RNAs to allow the sequencing of the full-length transcript using the SMART technique for the retrotranscription. Because of myofibers are multinucleated cells and because of it is better to perform single cell sequencing on fresh tissues we described the single-nucleus sequencing that can be applied to frozen tissues. The chapter concludes with a description of bioinformatics approaches to evaluate differential expression from single-cell or single-nucleus RNA sequencing.

Distinct amyloid-β and tau-associated microglia profiles in Alzheimer's disease

Alzheimer’s disease (AD) is the most prevalent form of dementia and is characterized by abnormal extracellular aggregates of amyloid-β and intraneuronal hyperphosphorylated tau tangles and neuropil threads. Microglia, the tissue-resident macrophages of the central nervous system (CNS), are important for CNS homeostasis and implicated in AD pathology. In amyloid mouse models, a phagocytic/activated microglia phenotype has been identified. How increasing levels of amyloid-β and tau pathology affect human microglia transcriptional profiles is unknown. Here, we performed snRNAseq on 482,472 nuclei from non-demented control brains and AD brains containing only amyloid-β plaques or both amyloid-β plaques and tau pathology. Within the microglia population, distinct expression profiles were identified of which two were AD pathology-associated. The phagocytic/activated AD1-microglia population abundance strongly correlated with tissue amyloid-β load and localized to amyloid-β plaques. The AD2-microglia abundance strongly correlated with tissue phospho-tau load and these microglia were more abundant in samples with overt tau pathology. This full characterization of human disease-associated microglia phenotypes provides new insights in the pathophysiological role of microglia in AD and offers new targets for microglia-state-specific therapeutic strategies.

Heterogeneity of midgut cells and their differential responses to blood meal ingestion by the mosquito, Aedes aegypti

Mosquitoes are the most notorious hematophagous insects and due to their blood feeding behavior and genetic compatibility, numerous mosquito species are highly efficient vectors for certain human pathogenic parasites and viruses. The mosquito midgut is the principal organ of blood meal digestion and nutrient absorption. It is also the initial site of infection with blood meal acquired parasites and viruses. We conducted an analysis based on single-nucleus RNA sequencing (snRNA-Seq) to assess the cellular diversity of the midgut and how individual cells respond to blood meal ingestion to facilitate its digestion. Our study revealed the presence of 20 distinguishable cell-type clusters in the female midgut of Aedes aegypti. The identified cell types included intestinal stem cells (ISC), enteroblasts (EB), differentiating EB (dEB), enteroendocrine cells (EE), enterocytes (EC), EC-like cells, cardia cells, and visceral muscle (VM) cells. Blood meal ingestion dramatically changed the overall midgut cell type composition, profoundly increasing the proportions of ISC and three EC/EC-like clusters. In addition, transcriptional profiles of all cell types were strongly affected while genes involved in various metabolic processes were significantly upregulated. Our study provides a basis for further physiological and molecular studies on blood digestion, nutrient absorption, and cellular homeostasis in the mosquito midgut.