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.

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.

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.

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.

Single nuclei transcriptomics delineates complex immune and kidney cell interactions contributing to kidney allograft fibrosis.

Chronic allograft dysfunction (CAD), characterized histologically by interstitial fibrosis and tubular atrophy, is the major cause of kidney allograft loss. Here, using single nuclei RNA sequencing and transcriptome analysis, we identified the origin, functional heterogeneity, and regulation of fibrosis-forming cells in kidney allografts with CAD. A robust technique was used to isolate individual nuclei from kidney allograft biopsies and successfully profiled 23,980 nuclei from five kidney transplant recipients with CAD and 17,913 nuclei from three patients with normal allograft function. Our analysis revealed two distinct states of fibrosis in CAD; low and high extracellular matrix (ECM) with distinct kidney cell subclusters, immune cell types, and transcriptional profiles. Imaging mass cytometry analysis confirmed increased ECM deposition at the protein level. Proximal tubular cells transitioned to an injured mixed tubular (MT1) phenotype comprised of activated fibroblasts and myofibroblast markers, generated provisional ECM which recruited inflammatory cells, and served as the main driver of fibrosis. MT1 cells in the high ECM state achieved replicative repair evidenced by dedifferentiation and nephrogenic transcriptional signatures. MT1 in the low ECM state showed decreased apoptosis, decreased cycling tubular cells, and severe metabolic dysfunction, limiting the potential for repair. Activated B, T and plasma cells were increased in the high ECM state, while macrophage subtypes were increased in the low ECM state. Intercellular communication between kidney parenchymal cells and donor-derived macrophages, detected several years post-transplantation, played a key role in injury propagation. Thus, our study identified novel molecular targets for interventions aimed to ameliorate or prevent allograft fibrogenesis in kidney transplant recipients.

Dissociation protocols used for sarcoma tissues bias the transcriptome observed in single-cell and single-nucleus RNA sequencing.

BACKGROUND: Single-cell RNA-seq has emerged as an innovative technology used to study complex tissues and characterize cell types, states, and lineages at a single-cell level. Classification of bulk tumors by their individual cellular constituents has also created new opportunities to generate single-cell atlases for many organs, cancers, and developmental models. Despite the tremendous promise of this technology, recent evidence studying epithelial tissues and diverse carcinomas suggests the methods used for tissue processing, cell disaggregation, and preservation can significantly bias gene expression and alter the observed cell types. To determine whether sarcomas - tumors of mesenchymal origin - are subject to the same technical artifacts, we profiled patient-derived tumor explants (PDXs) propagated from three aggressive subtypes: osteosarcoma (OS), Ewing sarcoma (ES), desmoplastic small round cell tumor (DSRCT). Given the rarity of these sarcoma subtypes, we explored whether single-nuclei RNA-seq from more widely available archival frozen specimens could accurately be identified by gene expression signatures linked to tissue phenotype or pathognomonic fusion proteins. RESULTS: We systematically assessed dissociation methods across different sarcoma subtypes. We compared gene expression from single-cell and single-nucleus RNA-sequencing of 125,831 whole-cells and nuclei from ES, DSRCT, and OS PDXs. We detected warm dissociation artifacts in single-cell samples and gene length bias in single-nucleus samples. Classic sarcoma gene signatures were observed regardless of the dissociation method. In addition, we showed that dissociation method biases could be computationally corrected. CONCLUSIONS: We highlighted transcriptional biases, including warm dissociation and gene-length biases, introduced by the dissociation method for various sarcoma subtypes. This work is the first to characterize how the dissociation methods used for sc/snRNA-seq may affect the interpretation of the molecular features in sarcoma PDXs.

Patch-seq: Advances and Biological Applications.

Multimodal analysis of gene-expression patterns, electrophysiological properties, and morphological phenotypes at the single-cell/single-nucleus level has been arduous because of the diversity and complexity of neurons. The emergence of Patch-sequencing (Patch-seq) directly links transcriptomics, morphology, and electrophysiology, taking neuroscience research to a multimodal era. In this review, we summarized the development of Patch-seq and recent applications in the cortex, hippocampus, and other nervous systems. Through generating multimodal cell type atlases, targeting specific cell populations, and correlating transcriptomic data with phenotypic information, Patch-seq has provided new insight into outstanding questions in neuroscience. We highlight the challenges and opportunities of Patch-seq in neuroscience and hope to shed new light on future neuroscience research.

Single-Nucleus RNA-Seq: Open the Era of Great Navigation for FFPE Tissue.

Single-cell sequencing (scRNA-seq) has revolutionized our ability to explore heterogeneity and genetic variations at the single-cell level, opening up new avenues for understanding disease mechanisms and cell-cell interactions. Single-nucleus RNA-sequencing (snRNA-seq) is emerging as a promising solution to scRNA-seq due to its reduced ionized transcription bias and compatibility with richer samples. This approach will provide an exciting opportunity for in-depth exploration of billions of formalin-fixed paraffin-embedded (FFPE) tissues. Recent advancements in single-cell/nucleus gene expression workflows tailored for FFPE tissues have demonstrated their feasibility and provided crucial guidance for future studies utilizing FFPE specimens. In this review, we provide a broad overview of the nuclear preparation strategies, the latest technologies of snRNA-seq applicable to FFPE samples. Finally, the limitations and potential technical developments of snRNA-seq in FFPE samples are summarized. The development of snRNA-seq technologies for FFPE samples will lay a foundation for transcriptomic studies of valuable samples in clinical medicine and human sample banks.