Non-canonical odor coding in the mosquito.

Aedes aegypti mosquitoes are a persistent human foe, transmitting arboviruses including dengue when they feed on human blood. Mosquitoes are intensely attracted to body odor and carbon dioxide, which they detect using ionotropic chemosensory receptors encoded by three large multi-gene families. Genetic mutations that disrupt the olfactory system have modest effects on human attraction, suggesting redundancy in odor coding. The canonical view is that olfactory sensory neurons each express a single chemosensory receptor that defines its ligand selectivity. We discovered that Ae. aegypti uses a different organizational principle, with many neurons co-expressing multiple chemosensory receptor genes. In vivo electrophysiology demonstrates that the broad ligand-sensitivity of mosquito olfactory neurons depends on this non-canonical co-expression. The redundancy afforded by an olfactory system in which neurons co-express multiple chemosensory receptors may increase the robustness of the mosquito olfactory system and explain our long-standing inability to disrupt the detection of humans by mosquitoes.

Cell-type-specific effects of autism-associated 15q duplication syndrome in the human brain.

Recurrent copy-number variation represents one of the most well-established genetic drivers in neurodevelopmental disorders, including autism spectrum disorder. Duplication of 15q11-q13 (dup15q) is a well-described neurodevelopmental syndrome that increases the risk of autism more than 40-fold. However, the effects of this duplication on gene expression and chromatin accessibility in specific cell types in the human brain remain unknown. To identify the cell-type-specific transcriptional and epigenetic effects of dup15q in the human frontal cortex, we conducted single-nucleus RNA sequencing and multi-omic sequencing on dup15q-affected individuals (n = 6) as well as individuals with non-dup15q autism (n = 7) and neurotypical control individuals (n = 7). Cell-type-specific differential expression analysis identified significantly regulated genes, critical biological pathways, and differentially accessible genomic regions. Although there was overall increased gene expression across the duplicated genomic region, cellular identity represented an important factor mediating gene-expression changes. As compared to other cell types, neuronal subtypes showed greater upregulation of gene expression across a critical region within the duplication. Genes that fell within the duplicated region and had high baseline expression in control individuals showed only modest changes in dup15q, regardless of cell type. Of note, dup15q and autism had largely distinct signatures of chromatin accessibility but shared the majority of transcriptional regulatory motifs, suggesting convergent biological pathways. However, the transcriptional binding-factor motifs implicated in each condition implicated distinct biological mechanisms: neuronal JUN and FOS networks in autism vs. an inflammatory transcriptional network in dup15q microglia. This work provides a cell-type-specific analysis of how dup15q changes gene expression and chromatin accessibility in the human brain, and it finds evidence of marked cell-type-specific effects of this genetic driver. These findings have implications for guiding therapeutic development in dup15q syndrome, as well as understanding the functional effects of copy-number variants more broadly in neurodevelopmental disorders.

Benchmarking mapping algorithms for cell-type annotating in mouse brain by integrating single-nucleus RNA-seq and Stereo-seq data.

Limited gene capture efficiency and spot size of spatial transcriptome (ST) data pose significant challenges in cell-type characterization. The heterogeneity and complexity of cell composition in the mammalian brain make it more challenging to accurately annotate ST data from brain. Many algorithms attempt to characterize subtypes of neuron by integrating ST data with single-nucleus RNA sequencing (snRNA-seq) or single-cell RNA sequencing. However, assessing the accuracy of these algorithms on Stereo-seq ST data remains unresolved. Here, we benchmarked 9 mapping algorithms using 10 ST datasets from four mouse brain regions in two different resolutions and 24 pseudo-ST datasets from snRNA-seq. Both actual ST data and pseudo-ST data were mapped using snRNA-seq datasets from the corresponding brain regions as reference data. After comparing the performance across different areas and resolutions of the mouse brain, we have reached the conclusion that both robust cell-type decomposition and SpatialDWLS demonstrated superior robustness and accuracy in cell-type annotation. Testing with publicly available snRNA-seq data from another sequencing platform in the cortex region further validated our conclusions. Altogether, we developed a workflow for assessing suitability of mapping algorithm that fits for ST datasets, which can improve the efficiency and accuracy of spatial data annotation.

Glial dysregulation in the human brain in fragile X-associated tremor/ataxia syndrome.

Short trinucleotide expansions at the FMR1 locus are associated with the late-onset condition fragile X-associated tremor/ataxia syndrome (FXTAS), which shows very different clinical and pathological features from fragile X syndrome (associated with longer expansions), with no clear molecular explanation for these marked differences. One prevailing theory posits that the shorter, premutation expansion uniquely causes extreme neurotoxic increases in FMR1 mRNA (i.e., four to eightfold increases), but evidence to support this hypothesis is largely derived from analysis of peripheral blood. We applied single-nucleus RNA sequencing to postmortem frontal cortex and cerebellum from 7 individuals with premutation and matched controls (n = 6) to assess cell type-specific molecular neuropathology. We found only modest upregulation (~1.3-fold) of FMR1 in some glial populations associated with premutation expansions. In premutation cases, we also identified decreased astrocyte proportions in the cortex. Differential expression and gene ontology analysis demonstrated altered neuroregulatory roles of glia. Using network analyses, we identified cell type-specific and region-specific patterns of FMR1 protein target gene dysregulation unique to premutation cases, with notable network dysregulation in the cortical oligodendrocyte lineage. We used pseudotime trajectory analysis to determine how oligodendrocyte development was altered and identified differences in early gene expression in oligodendrocyte trajectories in premutation cases specifically, implicating early cortical glial developmental perturbations. These findings challenge dogma regarding extremely elevated FMR1 increases in FXTAS and implicate glial dysregulation as a critical facet of premutation pathophysiology, representing potential unique therapeutic targets directly derived from the human condition.

Distinct myofibre domains of the human myotendinous junction revealed by single-nucleus RNA sequencing.

The myotendinous junction (MTJ) is a specialized domain of the multinucleated myofibre that is faced with the challenge of maintaining robust cell-matrix contact with the tendon under high mechanical stress and strain. Here, we profiled 24,124 nuclei in semitendinosus muscle-tendon samples from three healthy males by using single-nucleus RNA sequencing (snRNA-seq), alongside spatial transcriptomics, to gain insight into the genes characterizing this specialization in humans. We identified a cluster of MTJ myonuclei represented by 47 enriched transcripts, of which the presence of ABI3BP, ABLIM1, ADAMTSL1, BICD1, CPM, FHOD3, FRAS1 and FREM2 was confirmed at the MTJ at the protein level in immunofluorescence assays. Four distinct subclusters of MTJ myonuclei were apparent, comprising two COL22A1-expressing subclusters and two subclusters lacking COL22A1 expression but with differing fibre type profiles characterized by expression of either MYH7 or MYH1 and/or MYH2. Our findings reveal distinct myonuclei profiles of the human MTJ, which represents a weak link in the musculoskeletal system that is selectively affected in pathological conditions ranging from muscle strains to muscular dystrophies.

Transcriptional and functional motifs defining renal function revealed by single-nucleus RNA sequencing.

Recent advances in single-cell sequencing provide a unique opportunity to gain novel insights into the diversity, lineage, and functions of cell types constituting a tissue/organ. Here, we performed a single-nucleus study of the adult Drosophila renal system, consisting of Malpighian tubules and nephrocytes, which shares similarities with the mammalian kidney. We identified 11 distinct clusters representing renal stem cells, stellate cells, regionally specific principal cells, garland nephrocyte cells, and pericardial nephrocytes. Characterization of the transcription factors specific to each cluster identified fruitless (fru) as playing a role in stem cell regeneration and Hepatocyte nuclear factor 4 (Hnf4) in regulating glycogen and triglyceride metabolism. In addition, we identified a number of genes, including Rho guanine nucleotide exchange factor at 64C (RhoGEF64c), Frequenin 2 (Frq2), Prip, and CG1093 that are involved in regulating the unusual star shape of stellate cells. Importantly, the single-nucleus dataset allows visualization of the expression at the organ level of genes involved in ion transport and junctional permeability, providing a systems-level view of the organization and physiological roles of the tubules. Finally, a cross-species analysis allowed us to match the fly kidney cell types to mouse kidney cell types and planarian protonephridia, knowledge that will help the generation of kidney disease models. Altogether, our study provides a comprehensive resource for studying the fly kidney.

Keratinocytes drive the epithelial hyperplasia key to sea lice resistance in coho salmon.

BACKGROUND: Salmonid species have followed markedly divergent evolutionary trajectories in their interactions with sea lice. While sea lice parasitism poses significant economic, environmental, and animal welfare challenges for Atlantic salmon (Salmo salar) aquaculture, coho salmon (Oncorhynchus kisutch) exhibit near-complete resistance to sea lice, achieved through a potent epithelial hyperplasia response leading to rapid louse detachment. The molecular mechanisms underlying these divergent responses to sea lice are unknown. RESULTS: We characterized the cellular and molecular responses of Atlantic salmon and coho salmon to sea lice using single-nuclei RNA sequencing. Juvenile fish were exposed to copepodid sea lice (Lepeophtheirus salmonis), and lice-attached pelvic fin and skin samples were collected 12 h, 24 h, 36 h, 48 h, and 60 h after exposure, along with control samples. Comparative analysis of control and treatment samples revealed an immune and wound-healing response that was common to both species, but attenuated in Atlantic salmon, potentially reflecting greater sea louse immunomodulation. Our results revealed unique but complementary roles of three layers of keratinocytes in the epithelial hyperplasia response leading to rapid sea lice rejection in coho salmon. Our results suggest that basal keratinocytes direct the expansion and mobility of intermediate and, especially, superficial keratinocytes, which eventually encapsulate the parasite. CONCLUSIONS: Our results highlight the key role of keratinocytes in coho salmon's sea lice resistance and the diverged biological response of the two salmonid host species when interacting with this parasite. This study has identified key pathways and candidate genes that could be manipulated using various biotechnological solutions to improve Atlantic salmon sea lice resistance.

Single-nucleus RNA sequencing and lipidomics reveal characteristics of transcriptional and lipid composition in porcine longissimus dorsi muscle.

BACKGROUND: Global per capita meat consumption continues to rise, especially pork. Meat quality is influenced by the content of intramuscular fat (IMF) as a key factor. The longissimus dorsi muscle of Dahe pigs (DHM, IMF: 7.98% ± 1.96%) and Dahe black pigs (DHBM, IMF: 3.30% ± 0.64%) was studied to explore cellular heterogeneity and differentially expressed genes (DEGs) associated with IMF deposition using single-nucleus RNA sequencing (snRNA-seq). The lipid composition was then analyzed using non-targeted lipidomics. RESULTS: A total of seven cell subpopulations were identified, including myocytes, fibroblast/fibro/adipogenic progenitors (FAPs), satellite cells, endothelial cells, macrophages, pericytes, and adipocytes. Among them, FAPs and adipocytes were more focused because they could be associated with lipid deposition. 1623 DEGs in the FAPs subpopulation of DHBM were up-regulated compared with DHM, while 1535 were down-regulated. These DEGs enriched in the glycolysis/gluconeogenesis pathway. 109 DEGs were up-regulated and 806 were down-regulated in the adipocyte subpopulation of DHBM compared with DHM, which were mainly enriched in the PPAR signaling pathway and fatty acid (FA) biosynthesis. The expression level of PPARG, ABP4, LEP, and ACSL1 genes in DHM was higher than that in DHBM. Lipidomics reveals porcine lipid composition characteristics of muscle tissue. A total of 41 lipid classes and 2699 lipid species were identified in DHM and DHBM groups. The top ten relative peak areas of lipid classes in DHM and DHBM were triglyceride (TG), phosphatidylcholine (PC), phosphatidylethanolamine (PE), phosphatidylserine (PS), diglyceride (DG), cardiolipin (CL), ceramides (Cer), Simple Glc series (Hex1Cer), sphingomyelin (phSM), and phosphatidylinositol (PI). The relative peak areas of 35 lipid species in DHM were lower than DHBM, and 28 lipid species that were higher. There was a significant increase in the TG fatty acyl chains C6:0, C17:0, and C11:4, and a significant decrease in C16:0, C18:1, C18:2, and C22:4 in DHBM (p < 0.05). CONCLUSIONS: C16:0 FA may downregulate the expression level of PPARG gene, which leads to the downregulation of fat metabolism-related genes such as ACSL, PLIN2, and FABP4 in DHBM compared with DHM. This may be the reason that the lipid deposition ability of Dahe pigs is stronger than that of Dahe black pigs, which need further investigation.

snRNA-seq of adipose tissues reveals the potential cellular and molecular mechanisms of cold and disease resistance in Mongolian cattle.

BACKGROUND: Mongolian cattle are local breeds in northern China with excellent adaptability to harsh environmental conditions. Adipose tissues play essential roles in tolerance to cold and disease, but the associated cellular and molecular mechanisms are unclear. METHODS: Single-nucleus RNA sequencing (snRNA-seq) was performed on the adipose tissues from the subcutaneous (SAT), greater omentum (OAT) and perirenal (PAT) of 3 healthy cattle. The adipogenic trajectory was analyzed, and the functional roles of gene of interest were verified in vitro. RESULTS: There were different cell subpopulations in adipose tissues. The lipid-deposition adipocytes identified by the PTGER3 marker exhibited outstanding characteristics in SAT. In PAT and OAT, aldosterone was expressed to provide clues for the differential brown adipocytes. Among the DEGs by comparing OAT with SAT and PAT with OAT, C3 was significantly expressed in most of the cell populations in SAT. G0S2, LIPE, LPIN1, PTGER3 and RGCC took part in the adipogenic trajectory from preadipocyte commitment to mature adipocytes. S100A4 expression affected Ca2+ signaling and the expression of UCP1 ~ 3, FABP4 and PTGER3. CONCLUSION: The cell heterogeneity and genes expressed in adipose tissues of Mongolian cattle not only determine the endocrine and energy storage, but contribute to adapt to cold and disease resistance.

Cell-specific AHR-driven differential gene expression in the mouse liver cell following acute TCDD exposure.

2,3,7,8-Tetrachlorodibenzo-p-dioxin (TCDD) is a persistent environmental contaminant that disrupts hepatic function leading to steatotic liver disease (SLD)-like pathologies, such as steatosis, steatohepatitis, and fibrosis. These effects are mediated by the aryl hydrocarbon receptor following changes in gene expression. Although diverse cell types are involved, initial cell-specific changes in gene expression have not been reported. In this study, differential gene expression in hepatic cell types was examined in male C57BL/6 mice gavaged with 30 µg/kg of TCDD using single-nuclei RNA-sequencing. Ten liver cell types were identified with the proportions of most cell types remaining unchanged, except for neutrophils which increased at 72 h. Gene expression suggests TCDD induced genes related to oxidative stress in hepatocytes as early as 2 h. Lipid homeostasis was disrupted in hepatocytes, macrophages, B cells, and T cells, characterized by the induction of genes associated with lipid transport, steroid hormone biosynthesis, and the suppression of ß-oxidation, while linoleic acid metabolism was altered in hepatic stellate cells (HSCs), B cells, portal fibroblasts, and plasmacytoid dendritic cells. Pro-fibrogenic processes were also enriched, including the induction retinol metabolism genes in HSCs and the early induction of anti-fibrolysis genes in hepatocytes, endothelial cells, HSCs, and macrophages. Hepatocytes also had gene expression changes consistent with hepatocellular carcinoma. Collectively, these findings underscore the effects of TCDD in initiating SLD-like phenotypes and identified cell-specific gene expression changes related to oxidative stress, steatosis, fibrosis, cell proliferation and the development of HCC.

Dentate gyrus granule cells are a locus of pathology in Scn8a developmental encephalopathy.

Gain-of-function mutations in SCN8A cause developmental and epileptic encephalopathy (DEE), a disorder characterized by early-onset refractory seizures, deficits in motor and intellectual functions, and increased risk of sudden unexpected death in epilepsy. Altered activity of neurons in the corticohippocampal circuit has been reported in mouse models of DEE. We examined the effect of chronic seizures on gene expression in the hippocampus by single-nucleus RNA sequencing in mice expressing the patient mutation SCN8A-p.Asn1768Asp (N1768D). One hundred and eighty four differentially expressed genes were identified in dentate gyrus granule cells, many more than in other cell types. Electrophysiological recording from dentate gyrus granule cells demonstrated an elevated firing rate. Targeted reduction of Scn8a expression in the dentate gyrus by viral delivery of an shRNA resulted in doubling of median survival time from 4 months to 8 months, whereas delivery of shRNA to the CA1 and CA3 regions did not result in lengthened survival. These data indicate that granule cells of the dentate gyrus are a specific locus of pathology in SCN8A-DEE.

Unveiling the Ovarian Cell Characteristics and Molecular Mechanism of Prolificacy in Goats via Single-Nucleus Transcriptomics Data Analysis.

Increases in litter size, which are influenced by ovulation, are responsible for between 74% and 96% of the economic value of genetic progress, which influences selection. For the selection and breeding of highly prolific goats, genetic mechanisms underlying variations in litter size should be elucidated. Here, we used single-nucleus RNA sequencing to analyze 44,605 single nuclei from the ovaries of polytocous and monotocous goats during the follicular phase. Utilizing known reference marker genes, we identified 10 ovarian cell types characterized by distinct gene expression profiles, transcription factor networks, and reciprocal interaction signatures. An in-depth analysis of the granulosa cells revealed three subtypes exhibiting distinct gene expression patterns and dynamic regulatory mechanisms. Further investigation of cell-type-specific prolificacy-associated transcriptional changes elucidated that "downregulation of apoptosis", "increased anabolism", and "upstream responsiveness to hormonal stimulation" are associated with prolificacy. This study provides a comprehensive understanding of the cell-type-specific mechanisms and regulatory networks in the goat ovary, providing insights into the molecular mechanisms underlying goat prolificacy. These findings establish a vital foundation for furthering understanding of the molecular mechanisms governing folliculogenesis and for improving the litter size in goats via molecular design breeding.

Decoding brain memory formation by single-cell RNA sequencing.

To understand how distinct memories are formed and stored in the brain is an important and fundamental question in neuroscience and computational biology. A population of neurons, termed engram cells, represents the physiological manifestation of a specific memory trace and is characterized by dynamic changes in gene expression, which in turn alters the synaptic connectivity and excitability of these cells. Recent applications of single-cell RNA sequencing (scRNA-seq) and single-nucleus RNA sequencing (snRNA-seq) are promising approaches for delineating the dynamic expression profiles in these subsets of neurons, and thus understanding memory-specific genes, their combinatorial patterns and regulatory networks. The aim of this article is to review and discuss the experimental and computational procedures of sc/snRNA-seq, new studies of molecular mechanisms of memory aided by sc/snRNA-seq in human brain diseases and related mouse models, and computational challenges in understanding the regulatory mechanisms underlying long-term memory formation.

Impaired GABAergic regulation and developmental immaturity in interneurons derived from the medial ganglionic eminence in the tuberous sclerosis complex.

GABAergic interneurons play a critical role in maintaining neural circuit balance, excitation-inhibition regulation, and cognitive function modulation. In tuberous sclerosis complex (TSC), GABAergic neuron dysfunction contributes to disrupted network activity and associated neurological symptoms, assumingly in a cell type-specific manner. This GABAergic centric study focuses on identifying specific interneuron subpopulations within TSC, emphasizing the unique characteristics of medial ganglionic eminence (MGE)- and caudal ganglionic eminence (CGE)-derived interneurons. Using single-nuclei RNA sequencing in TSC patient material, we identify somatostatin-expressing (SST+) interneurons as a unique and immature subpopulation in TSC. The disrupted maturation of SST+ interneurons may undergo an incomplete switch from excitatory to inhibitory GABAergic signaling during development, resulting in reduced inhibitory properties. Notably, this study reveals markers of immaturity specifically in SST+ interneurons, including an abnormal NKCC1/KCC2 ratio, indicating an imbalance in chloride homeostasis crucial for the postsynaptic consequences of GABAergic signaling as well as the downregulation of GABAA receptor subunits, GABRA1, and upregulation of GABRA2. Further exploration of SST+ interneurons revealed altered localization patterns of SST+ interneurons in TSC brain tissue, concentrated in deeper cortical layers, possibly linked to cortical dyslamination. In the epilepsy context, our research underscores the diverse cell type-specific roles of GABAergic interneurons in shaping seizures, advocating for precise therapeutic considerations. Moreover, this study illuminates the potential contribution of SST+ interneurons to TSC pathophysiology, offering insights for targeted therapeutic interventions.

Xenografted human iPSC-derived neurons with the familial Alzheimer's disease APPV717I mutation reveal dysregulated transcriptome signatures linked to synaptic function and implicate LINGO2 as a disease signaling mediator.

Alzheimer's disease (AD) is the most common cause of dementia, and disease mechanisms are still not fully understood. Here, we explored pathological changes in human induced pluripotent stem cell (iPSC)-derived neurons carrying the familial AD APPV717I mutation after cell injection into the mouse forebrain. APPV717I mutant iPSCs and isogenic controls were differentiated into neurons revealing enhanced Aß42 production, elevated phospho-tau, and impaired neurite outgrowth in APPV717I neurons. Two months after transplantation, APPV717I and control neural cells showed robust engraftment but at 12 months post-injection, APPV717I grafts were smaller and demonstrated impaired neurite outgrowth compared to controls, while plaque and tangle pathology were not seen. Single-nucleus RNA-sequencing of micro-dissected grafts, performed 2 months after cell injection, identified significantly altered transcriptome signatures in APPV717I iPSC-derived neurons pointing towards dysregulated synaptic function and axon guidance. Interestingly, APPV717I neurons showed an increased expression of genes, many of which are also upregulated in postmortem neurons of AD patients including the transmembrane protein LINGO2. Downregulation of LINGO2 in cultured APPV717I neurons rescued neurite outgrowth deficits and reversed key AD-associated transcriptional changes related but not limited to synaptic function, apoptosis and cellular senescence. These results provide important insights into transcriptional dysregulation in xenografted APPV717I neurons linked to synaptic function, and they indicate that LINGO2 may represent a potential therapeutic target in AD.

Single-nucleus RNA transcriptome profiling reveals murine adipose tissue endothelial cell proliferation gene networks involved in obesity development.

Endothelial cells play an important role in the metabolism of adipose tissue (AT). This study aimed to analyze the changes that adipose tissue in AT endothelial cells undergo during the development of obesity, using single-nucleus RNA sequence (snRNA-seq). Mouse paraepididymal AT cells were subjected to snRNA-seq with the 10X Genomics platform. The cell types were then clustered using t-distributed stochastic neighbor embedding and unbiased computational informatics analyses. Protein-protein interactions network was established using the STRING database and visualized using Cytoscape. The dataset was subjected to differential gene enrichment analysis. In total, 21,333 cells acquired from 24 mouse paraepididymal AT samples were analyzed using snRNA-seq. This study identified 18 distinct clusters and annotated macrophages, fibroblasts, epithelial cells, T cells, endothelial cells, stem cells, neutrophil cells, and neutrophil cell types based on representative markers. Cluster 12 was defined as endothelial cells. The proportion of endothelial cells decreased with the development of obesity. Inflammatory factors, such as Vegfa and Prdm16 were upregulated in the medium obesity group but downregulated in the obesity group. Genes, such as Prox1, Erg, Flt4, Kdr, Flt1, and Pecam1 promoted the proliferation of AT endothelial cells and maintained the internal environment of AT. This study established a reference model and general framework for studying the mechanisms, biomarkers, and therapeutic targets of endothelial cell dysfunction-related diseases at the single-cell level.

Disease-Associated Neurotoxic Astrocyte Markers in Alzheimer Disease Based on Integrative Single-Nucleus RNA Sequencing.

Alzheimer disease (AD) is an irreversible neurodegenerative disease, and astrocytes play a key role in its onset and progression. The aim of this study is to analyze the characteristics of neurotoxic astrocytes and identify novel molecular targets for slowing down the progression of AD. Single-nucleus RNA sequencing (snRNA-seq) data were analyzed from various AD cohorts comprising about 210,654 cells from 53 brain tissue. By integrating snRNA-seq data with bulk RNA-seq data, crucial astrocyte types and genes associated with the prognosis of patients with AD were identified. The expression of neurotoxic astrocyte markers was validated using 5 × FAD and wild-type (WT) mouse models, combined with experiments such as western blot, quantitative real-time PCR (qRT-PCR), and immunofluorescence. A group of neurotoxic astrocytes closely related to AD pathology was identified, which were involved in inflammatory responses and pathways related to neuron survival. Combining snRNA and bulk tissue data, ZEP36L, AEBP1, WWTR1, PHYHD1, DST and RASL12 were identified as toxic astrocyte markers closely related to disease severity, significantly elevated in brain tissues of 5 × FAD mice and primary astrocytes treated with Aß. Among them, WWTR1 was significantly increased in astrocytes of 5 × FAD mice, driving astrocyte inflammatory responses, and has been identified as an important marker of neurotoxic astrocytes. snRNA-seq analysis reveals the biological functions of neurotoxic astrocytes. Six genes related to AD pathology were identified and validated, among which WWTR1 may be a novel marker of neurotoxic astrocytes.

Decoding the transcriptome of muscular dystrophy due to Ptrf deficiency using single-nucleus RNA sequencing.

Lacking PTRF (polymerase I and transcript release factor), an essential caveolae component, causes a secondary deficiency of caveolins resulting in muscular dystrophy. The transcriptome responses of different types of muscle fibers and mononuclear cells in skeletal muscle to muscular dystrophy caused by Ptrf deletion have not been explored. Here, we created muscular dystrophy mice by Ptrf knockout and applied single-nucleus RNA sequencing (snRNA-seq) to unveil the transcriptional changes of the skeletal muscle at single-nucleus resolution. 11 613 muscle nuclei (WT, 5838; Ptrf KO, 5775) were classified into 12 clusters corresponding to 11 nuclear types. Trajectory analysis revealed the potential transition between type IIb_1 and IIb_2 myonuclei upon muscular dystrophy. Functional enrichment analysis indicated that apoptotic signaling and enzyme-linked receptor protein signaling pathway were significantly enriched in type IIb_1 and IIb_2 myonuclei of Ptrf KO, respectively. The muscle structure development and the PI3K-AKT signaling pathway were significantly enriched in type IIa and IIx myonuclei of Ptrf KO. Meanwhile, metabolic pathway analysis showed a decrease in overall metabolic pathway activity of myonuclei subtypes upon muscular dystrophy, with the most decrease in type IIb_1 myonuclei. Gene regulatory network analysis found that the activity of Mef2c, Mef2d, Myf5, and Pax3 regulons was enhanced in type II myonuclei of Ptrf KO, especially in type IIb_2 myonuclei. In addition, we investigated the transcriptome changes in adipocytes and found that muscular dystrophy enhanced the lipid metabolic capacity of adipocytes. Our findings provide a valuable resource for exploring the molecular mechanism of muscular dystrophy due to Ptrf deficiency.

Empagliflozin and liraglutide ameliorate HFpEF in mice via augmenting the Erbb4 signaling pathway.

Heart failure with preserved ejection fraction (HFpEF) is closely associated with metabolic derangement. Sodium glucose cotransporter-2 inhibitors (SGLT2i) and glucagon-like peptide-1 receptor agonists (GLP-1RA) exert anti-HFpEF effects, but the underlying mechanisms remain unclear. In this study, we explored the anti-HFpEF effects of empagliflozin and liraglutide and the underlying molecular mechanisms in a mouse model of HFpEF. This model was established by high-fat diet (HFD) feeding plus Nω-nitro-L-arginine methyl ester (L-NAME) treatment. The mice were treated with empagliflozin (20 mg·kg-1·d-1, i.g.) or liraglutide (0.3 mg·kg-1·d-1, i.p.) or their combination for 4 weeks. At the end of the experimental protocol, cardiac function was measured using ultrasound, then mice were euthanized and heart, liver, and kidney tissues were collected. Nuclei were isolated from frozen mouse ventricular tissue for single-nucleus RNA-sequencing (snRNA-seq). We showed that administration of empagliflozin or liraglutide alone or in combination significantly improved diastolic function, ameliorated cardiomyocyte hypertrophy and cardiac fibrosis, as well as exercise tolerance but no synergism was observed in the combination group. Furthermore, empagliflozin and/or liraglutide lowered body weight, improved glucose metabolism, lowered blood pressure, and improved liver and kidney function. After the withdrawal of empagliflozin or liraglutide for 1 week, these beneficial effects tended to diminish. The snRNA-seq analysis revealed a subcluster of myocytes, in which Erbb4 expression was down-regulated under HFpEF conditions, and restored by empagliflozin or liraglutide. Pseudo-time trajectory analysis and cell-to-cell communication studies confirmed that the Erbb4 pathway was a prominent pathway essential for both drug actions. In the HFpEF mouse model, both empagliflozin and liraglutide reversed Erbb4 down-regulation. In rat h9c2 cells, we showed that palmitic acid- or high glucose-induced changes in PKCα and/or ERK1/2 phosphorylation at least in part through Erbb4. Collectively, the single-cell atlas reveals the anti-HFpEF mechanism of empagliflozin and liraglutide, suggesting that Erbb4 pathway represents a new therapeutic target for HFpEF. Effects and mechanisms of action of empagliflozin and liraglutide in HFpEF mice. HFpEF was induced with a high-fat diet and L-NAME for 15 weeks, and treatment with empagliflozin and liraglutide improved the HFpEF phenotype. Single nucleus RNA sequencing (snRNA-seq) was used to reveal the underlying mechanism of action of empagliflozin and liraglutide.

Distinct subpopulations of D1 medium spiny neurons exhibit unique transcriptional responsiveness to cocaine.

Drugs of abuse increase extracellular concentrations of dopamine in the nucleus accumbens (NAc), resulting in transcriptional alterations that drive long-lasting cellular and behavioral adaptations. While decades of research have focused on the transcriptional mechanisms by which drugs of abuse influence neuronal physiology and function, few studies have comprehensively defined NAc cell type heterogeneity in transcriptional responses to drugs of abuse. Here, we used single nucleus RNA-seq (snRNA-seq) to characterize the transcriptome of over 39,000 NAc cells from male and female adult Sprague-Dawley rats following acute or repeated cocaine experience. This dataset identified 16 transcriptionally distinct cell populations, including two populations of medium spiny neurons (MSNs) that express the Drd1 dopamine receptor (D1-MSNs). Critically, while both populations expressed classic marker genes of D1-MSNs, only one population exhibited a robust transcriptional response to cocaine. Validation of population-selective transcripts using RNA in situ hybridization revealed distinct spatial compartmentalization of these D1-MSN populations within the NAc. Finally, analysis of published NAc snRNA-seq datasets from non-human primates and humans demonstrated conservation of MSN subtypes across rat and higher order mammals, and further highlighted cell type-specific transcriptional differences across the NAc and broader striatum. These results highlight the utility in using snRNA-seq to characterize both cell type heterogeneity and cell type-specific responses to cocaine and provides a useful resource for cross-species comparisons of NAc cell composition.