Single-cell transcriptomics of 20 mouse organs creates a Tabula Muris

Adipose stem cells are sexually dimorphic cells with dual roles as preadipocytes and resident fibroblasts

Cell identities are defined by intrinsic transcriptional networks and spatio-temporal environmental factors. Here, we explored multiple factors that contribute to the identity of adipose stem cells, including anatomic location, microvascular neighborhood, and sex. Our data suggest that adipose stem cells serve a dual role as adipocyte precursors and fibroblast-like cells that shape the adipose tissue's extracellular matrix in an organotypic manner. We further find that adipose stem cells display sexual dimorphism regarding genes involved in estrogen signaling, homeobox transcription factor expression and the renin-angiotensin-aldosterone system. These differences could be attributed to sex hormone effects, developmental origin, or both. Finally, our data demonstrate that adipose stem cells are distinct from mural cells, and that the state of commitment to adipogenic differentiation is linked to their anatomic position in the microvascular niche. Our work supports the importance of sex and microvascular function in adipose tissue physiology.

A novel STING variant triggers endothelial toxicity and SAVI disease

Gain-of-function mutations in STING cause STING-associated vasculopathy with onset in infancy (SAVI) characterized by early-onset systemic inflammation, skin vasculopathy, and interstitial lung disease. Here, we report and characterize a novel STING variant (F269S) identified in a SAVI patient. Single-cell transcriptomics of patient bone marrow revealed spontaneous activation of interferon (IFN) and inflammatory pathways across cell types and a striking prevalence of circulating naïve T cells was observed. Inducible STING F269S expression conferred enhanced signaling through ligand-independent translocation of the protein to the Golgi, protecting cells from viral infections but preventing their efficient immune priming. Additionally, endothelial cell activation was promoted and further exacerbated by cytokine secretion by SAVI immune cells, resulting in inflammation and endothelial damage. Our findings identify STING F269S mutation as a novel pathogenic variant causing SAVI, highlight the importance of the crosstalk between endothelial and immune cells in the context of lung disease, and contribute to a better understanding of how aberrant STING activation can cause pathology.

Tumour-derived small extracellular vesicles act as a barrier to therapeutic nanoparticle delivery

Nanoparticles are promising for drug delivery applications, with several clinically approved products. However, attaining high nanoparticle accumulation in solid tumours remains challenging. Here we show that tumour cell-derived small extracellular vesicles (sEVs) block nanoparticle delivery to tumours, unveiling another barrier to nanoparticle-based tumour therapy. Tumour cells secrete large amounts of sEVs in the tumour microenvironment, which then bind to nanoparticles entering tumour tissue and traffic them to liver Kupffer cells for degradation. Knockdown of Rab27a, a gene that controls sEV secretion, decreases sEV levels and improves nanoparticle accumulation in tumour tissue. The therapeutic efficacy of messenger RNAs encoding tumour suppressing and proinflammatory proteins is greatly improved when co-encapsulated with Rab27a small interfering RNA in lipid nanoparticles. Together, our results demonstrate that tumour cell-derived sEVs act as a defence system against nanoparticle tumour delivery and that this system may be a potential target for improving nanoparticle-based tumour therapies.

scDiffusion: conditional generation of high-quality single-cell data using diffusion model

MOTIVATION

Single-cell RNA sequencing (scRNA-seq) data are important for studying the laws of life at single-cell level. However, it is still challenging to obtain enough high-quality scRNA-seq data. To mitigate the limited availability of data, generative models have been proposed to computationally generate synthetic scRNA-seq data. Nevertheless, the data generated with current models are not very realistic yet, especially when we need to generate data with controlled conditions. In the meantime, diffusion models have shown their power in generating data with high fidelity, providing a new opportunity for scRNA-seq generation.

RESULTS

In this study, we developed scDiffusion, a generative model combining the diffusion model and foundation model to generate high-quality scRNA-seq data with controlled conditions. We designed multiple classifiers to guide the diffusion process simultaneously, enabling scDiffusion to generate data under multiple condition combinations. We also proposed a new control strategy called Gradient Interpolation. This strategy allows the model to generate continuous trajectories of cell development from a given cell state. Experiments showed that scDiffusion could generate single-cell gene expression data closely resembling real scRNA-seq data. Also, scDiffusion can conditionally produce data on specific cell types including rare cell types. Furthermore, we could use the multiple-condition generation of scDiffusion to generate cell type that was out of the training data. Leveraging the Gradient Interpolation strategy, we generated a continuous developmental trajectory of mouse embryonic cells. These experiments demonstrate that scDiffusion is a powerful tool for augmenting the real scRNA-seq data and can provide insights into cell fate research.

AVAILABILITY AND IMPLEMENTATION

scDiffusion is openly available at the GitHub repository https://github.com/EperLuo/scDiffusion or Zenodo https://zenodo.org/doi/10.5281/zenodo.13268742.

TLNRD1 is a CCM complex component and regulates endothelial barrier integrity

We previously identified talin rod domain-containing protein 1 (TLNRD1) as a potent actin-bundling protein in vitro. Here, we report that TLNRD1 is expressed in the vasculature in vivo. Its depletion leads to vascular abnormalities in vivo and modulation of endothelial cell monolayer integrity in vitro. We demonstrate that TLNRD1 is a component of the cerebral cavernous malformations (CCM) complex through its direct interaction with CCM2, which is mediated by a hydrophobic C-terminal helix in CCM2 that attaches to a hydrophobic groove on the four-helix domain of TLNRD1. Disruption of this binding interface leads to CCM2 and TLNRD1 accumulation in the nucleus and actin fibers. Our findings indicate that CCM2 controls TLNRD1 localization to the cytoplasm and inhibits its actin-bundling activity and that the CCM2-TLNRD1 interaction impacts endothelial actin stress fiber and focal adhesion formation. Based on these results, we propose a new pathway by which the CCM complex modulates the actin cytoskeleton and vascular integrity.

Trefoil Factor 2 Expressed by the Murine Pancreatic Acinar Cells Is Required for the Development of Islets and for β-Cell Function During Aging

Exocrine-to-endocrine cross talk in the pancreas is crucial to maintain β-cell function. However, the molecular mechanisms underlying this cross talk are largely undefined. Trefoil factor 2 (Tff2) is a secreted factor known to promote the proliferation of β-cells in vitro, but its physiological role in vivo in the pancreas is unknown. Also, it remains unclear which pancreatic cell type expresses Tff2 protein. We therefore created a mouse model with a conditional knockout of Tff2 in the murine pancreas. We find that the Tff2 protein is preferentially expressed in acinar but not ductal or endocrine cells. Tff2 deficiency in the pancreas reduces β-cell mass on embryonic day 16.5. However, homozygous mutant mice are born without a reduction of β-cells and with acinar Tff3 compensation by day 7. When mice are aged to 1 year, both male and female homozygous and male heterozygous mutants develop impaired glucose tolerance without affected insulin sensitivity. Perifusion analysis reveals that the second phase of glucose-stimulated insulin secretion from islets is reduced in aged homozygous mutant compared with controls. Collectively, these results demonstrate a previously unknown role of Tff2 as an exocrine acinar cell-derived protein required for maintaining functional endocrine β-cells in mice.

ARTICLE HIGHLIGHTS

ZBTB46 coordinates angiogenesis and immunity to control tumor outcome

Tumor angiogenesis and immunity show an inverse correlation in cancer progression and outcome1. Here, we report that ZBTB46, a repressive transcription factor and a widely accepted marker for classical dendritic cells (DCs)2,3, controls both tumor angiogenesis and immunity. Zbtb46 was downregulated in both DCs and endothelial cells by tumor-derived factors to facilitate robust tumor growth. Zbtb46 downregulation led to a hallmark pro-tumor microenvironment (TME), including dysfunctional vasculature and immunosuppressive conditions. Analysis of human cancer data revealed a similar association of low ZBTB46 expression with an immunosuppressive TME and a worse prognosis. In contrast, enforced Zbtb46 expression led to TME changes to restrict tumor growth. Mechanistically, Zbtb46-deficient endothelial cells were highly angiogenic, and Zbtb46-deficient bone marrow progenitors upregulated Cebpb and diverted the DC program to immunosuppressive myeloid lineage output, potentially explaining the myeloid lineage skewing phenomenon in cancer4. Conversely, enforced Zbtb46 expression normalized tumor vessels and, by suppressing Cebpb, skewed bone marrow precursors toward immunostimulatory myeloid lineage output, leading to an immune-hot TME. Remarkably, Zbtb46 mRNA treatment synergized with anti-PD1 immunotherapy to improve tumor management in preclinical models. These findings identify ZBTB46 as a critical factor for angiogenesis and for myeloid lineage skewing in cancer and suggest that maintaining its expression could have therapeutic benefits.

Rbm24 modulates neuronal RNA splicing to restrict cognitive dysfunction

Synaptic dysfunction is associated with early neurodegenerative changes and cognitive deficits. Neuronal cell-specific alternative splicing (AS) programs exclusively encode unique neuron- and synapse-specific proteins. However, it remains unclear whether splicing disturbances in neurons influence the pathogenesis of cognitive impairment. Here, we observed that RNA-binding motif protein 24 (RBM24) expression was decreased in Alzheimer's disease (AD) patients. Using conditional RBM24 knockout mice, we demonstrated that deletion of RBM24 in the brain resulted in learning and memory impairment. Electrophysiological recordings from hippocampal slices from mice lacking RBM24 revealed multiple defects in excitatory synaptic function and plasticity. Furthermore, RNA sequencing and splicing analysis showed that RBM24 regulates a network of genes related to cognitive function. Deletion of RBM24 disrupted the AS of synapse-associated genes, including GluR2 and Prrt1, the major disease genes involved in cognitive impairment and memory loss, leading to cognitive dysfunction. Together, our results suggest that the regulation of mRNA splicing by RBM24 is a key process involved in maintaining normal synaptic function and provide novel mechanistic insights into the pathogenesis of AD.

Single-cell analysis of chromatin and expression reveals age- and sex-associated alterations in the human heart

Sex differences and age-related changes in the human heart at the tissue, cell, and molecular level have been well-documented and many may be relevant for cardiovascular disease. However, how molecular programs within individual cell types vary across individuals by age and sex remains poorly characterized. To better understand this variation, we performed single-nucleus combinatorial indexing (sci) ATAC- and RNA-Seq in human heart samples from nine donors. We identify hundreds of differentially expressed genes by age and sex and find epigenetic signatures of variation in ATAC-Seq data in this discovery cohort. We then scale up our single-cell RNA-Seq analysis by combining our data with five recently published single nucleus RNA-Seq datasets of healthy adult hearts. We find variation such as metabolic alterations by sex and immune changes by age in differential expression tests, as well as alterations in abundance of cardiomyocytes by sex and neurons with age. In addition, we compare our adult-derived ATAC-Seq profiles to analogous fetal cell types to identify putative developmental-stage-specific regulatory factors. Finally, we train predictive models of cell-type-specific RNA expression levels utilizing ATAC-Seq profiles to link distal regulatory sequences to promoters, quantifying the predictive value of a simple TF-to-expression regulatory grammar and identifying cell-type-specific TFs. Our analysis represents the largest single-cell analysis of cardiac variation by age and sex to date and provides a resource for further study of healthy cardiac variation and transcriptional regulation at single-cell resolution.

Endothelin 3/EDNRB signaling induces thermogenic differentiation of white adipose tissue

Thermogenic adipose tissue, consisting of brown and beige fat, regulates nutrient utilization and energy metabolism. Human brown fat is relatively scarce and decreases with obesity and aging. Hence, inducing thermogenic differentiation of white fat offers an attractive way to enhance whole-body metabolic capacity. Here, we show the role of endothelin 3 (EDN3) and endothelin receptor type B (EDNRB) in promoting the browning of white adipose tissue (WAT). EDNRB overexpression stimulates thermogenic differentiation of human white preadipocytes through cAMP-EPAC1-ERK activation. In mice, cold induces the expression of EDN3 and EDNRB in WAT. Deletion of EDNRB in adipose progenitor cells impairs cold-induced beige adipocyte formation in WAT, leading to excessive weight gain, glucose intolerance, and insulin resistance upon high-fat feeding. Injection of EDN3 into WAT promotes browning and improved whole-body glucose metabolism. The findings shed light on the mechanism of WAT browning and offer potential therapeutics for obesity and metabolic disorders.

Profiling cell identity and tissue architecture with single-cell and spatial transcriptomics

Single-cell transcriptomics has broadened our understanding of cellular diversity and gene expression dynamics in healthy and diseased tissues. Recently, spatial transcriptomics has emerged as a tool to contextualize single cells in multicellular neighbourhoods and to identify spatially recurrent phenotypes, or ecotypes. These technologies have generated vast datasets with targeted-transcriptome and whole-transcriptome profiles of hundreds to millions of cells. Such data have provided new insights into developmental hierarchies, cellular plasticity and diverse tissue microenvironments, and spurred a burst of innovation in computational methods for single-cell analysis. In this Review, we discuss recent advancements, ongoing challenges and prospects in identifying and characterizing cell states and multicellular neighbourhoods. We discuss recent progress in sample processing, data integration, identification of subtle cell states, trajectory modelling, deconvolution and spatial analysis. Furthermore, we discuss the increasing application of deep learning, including foundation models, in analysing single-cell and spatial transcriptomics data. Finally, we discuss recent applications of these tools in the fields of stem cell biology, immunology, and tumour biology, and the future of single-cell and spatial transcriptomics in biological research and its translation to the clinic.

Transcriptomic analysis unveils bona fide molecular signatures of microglia under conditions of homeostasis and viral encephalitis

Microglia serve as a front-line defense against neuroinvasive viral infection, however, determination of their actual transcriptional profiles under conditions of health and disease is challenging. Here, we used various experimental approaches to delineate the transcriptional landscape of microglia during viral infection. Intriguingly, multiple activation genes were found to be artificially induced in sorted microglia and we demonstrated that shear stress encountered during cell sorting was one of the key inducers. Post-hoc analysis revealed that publicly available large-scale single-cell RNA sequencing datasets were significantly tainted by aberrant signatures that are associated with cell sorting. By exploiting the ribosomal tagging approach, we developed a strategy to enrich microglia-specific transcripts by comparing immunoprecipitated RNA with total RNA. Such enriched transcripts were instrumental in defining bona fide signatures of microglia under conditions of health and virus infection. These unified microglial signatures may serve as a benchmark to retrospectively assess ex vivo artefacts from available atlases. Leveraging the microglial translatome, we found enrichment of genes implicated in T-cell activation and cytokine production during the course of VSV infection. These data linked microglia with T-cell re-stimulation and further underscored that microglia are involved in shaping antiviral T-cell responses in the brain. Collectively, our study defines the transcriptional landscape of microglia under steady state and during viral encephalitis and highlights cellular interactions between microglia and T cells that contribute to the control of virus dissemination.

Single-cell chromatin accessibility and transposable element landscapes reveal shared features of tissue-residing immune cells

Tissue adaptation is required for regulatory T (Treg) cell function within organs. Whether this program shares aspects with other tissue-localized immune populations is unclear. Here, we analyzed single-cell chromatin accessibility data, including the transposable element (TE) landscape of CD45+ immune cells from colon, skin, adipose tissue, and spleen. We identified features of organ-specific tissue adaptation across different immune cells. Focusing on tissue Treg cells, we found conservation of the Treg tissue adaptation program in other tissue-localized immune cells, such as amphiregulin-producing T helper (Th)17 cells. Accessible TEs can act as regulatory elements, but their contribution to tissue adaptation is not understood. TE landscape analysis revealed an enrichment of specific transcription factor binding motifs in TE regions within accessible chromatin peaks. TEs, specifically from the LTR family, were located in enhancer regions and associated with tissue adaptation. These findings broaden our understanding of immune tissue residency and provide an important step toward organ-specific immune interventions.

Integrative clinical and preclinical studies identify FerroTerminator1 as a potent therapeutic drug for MASH

The complex etiological factors associated with metabolic dysfunction-associated fatty liver disease (MAFLD), including perturbed iron homeostasis, and the unclear nature by which they contribute to disease progression have resulted in a limited number of effective therapeutic interventions. Here, we report that patients with metabolic dysfunction-associated steatohepatitis (MASH), a pathological subtype of MAFLD, exhibit excess hepatic iron and that it has a strong positive correlation with disease progression. FerroTerminator1 (FOT1) effectively reverses liver injury across multiple MASH models without notable toxic side effects compared with clinically approved iron chelators. Mechanistically, our multi-omics analyses reveal that FOT1 concurrently inhibits hepatic iron accumulation and c-Myc-Acsl4-triggered ferroptosis in various MASH models. Furthermore, MAFLD cohort studies suggest that serum ferritin levels might serve as a predictive biomarker for FOT1-based therapy in MASH. These findings provide compelling evidence to support FOT1 as a promising novel therapeutic option for all stages of MAFLD and for future clinical trials.

Liver Extracellular Vesicles and Particles Enriched β-Sitosterol Effectively Promote Liver Regeneration in Mice

Background

The liver's regenerative capacity allows it to repair itself after injury. Extracellular vesicles and particles (EVPs) in the liver's interstitial space are crucial for signal transduction, metabolism, and immune regulation. Understanding the role and mechanism of liver-derived EVPs in regeneration is significant, particularly after partial hepatectomy, where the mechanisms remain unclear.

Methods

A 70% hepatectomy model was established in mice, and EVPs were isolated and characterized using electron microscopy, nanocharacterization, and Western blot analysis. Combined metabolomic and transcriptomic analyses revealed β-sitosterol enrichment in EVPs and activation of the Hedgehog signaling pathway during regeneration. The role of β-sitosterol in EVPs on the Hedgehog pathway and its targets were identified using qRT-PCR, Western blot analysis. The regulation of carnitine synthesis by this pathway was determined using a dual luciferase assay. The effect of a β-sitosterol diet on liver regeneration was verified in mice.

Results

After 70% hepatectomy, the liver successfully regenerated without liver failure or death. At 24 hours post-surgery, tissue staining showed transient regeneration-associated steatosis (TRAS), with increased Ki67 positivity at 48 hours. EVPs displayed a spherical lipid bilayer structure with particle sizes of 70-130 nm. CD9, CD63, and CD81 in liver-derived EVPs were confirmed. Transcriptomic and metabolomic analyses showed EVPs supplementation significantly promoted carnitine synthesis and fatty acid oxidation. Tissue staining confirmed accelerated TRAS resolution and enhanced liver regeneration with EVP supplementation. Mass spectrometry identified β-sitosterol in EVPs, which binds to Smo protein, activating the Hedgehog pathway. This led to the nuclear transport of Gli3, stimulating Setd5 transcription and inducing carnitine synthesis, thereby accelerating fatty acid oxidation. Mice with increased β-sitosterol intake showed faster TRAS resolution and liver regeneration compared to controls.

Conclusion

Liver-derived EVPs promote regeneration after partial hepatectomy. β-sitosterol from EVPs accelerates fatty acid oxidation and promotes liver regeneration by activating Hedgehog signaling pathway.

Transcriptomic Signatures in Lung Allografts and Their Therapeutic Implications

Ex vivo lung perfusion (EVLP) is a well-established method of lung preservation in clinical transplantation. Transcriptomic analyses of cells and tissues uncover gene expression patterns which reveal granular molecular pathways and cellular programs under various conditions. Coupling EVLP and transcriptomics may provide insights into lung allograft physiology at a molecular level with the potential to develop targeted therapies to enhance or repair the donor lung. This review examines the current landscape of transcriptional analysis of lung allografts in the context of state-of-the-art therapeutics that have been developed to optimize lung allograft function.

Physical modeling of embryonic transcriptomes identifies collective modes of gene expression

Starting from one totipotent cell, complex multicellular organisms form through a series of differentiation and morphogenetic events, culminating in a multitude of cell types arranged in a functional and intricate spatial pattern. To do so, cells coordinate with each other, resulting in dynamics which follow a precise developmental trajectory, constraining the space of possible embryo-to-embryo variation. Using recent single-cell sequencing data of early ascidian embryos, we leverage natural variation together with modeling and inference techniques from statistical physics to investigate development at the level of a complete interconnected embryo - an embryonic transcriptome. After developing a robust and biophysically motivated approach to identifying distinct transcriptomic states or cell types, a statistical analysis reveals correlations within embryos and across cell types demonstrating the presence of collective variation. From these intra-embryo correlations, we infer minimal networks of cell-cell interactions, which reveal the collective modes of gene expression. Our work demonstrates how the existence and nature of spatial interactions along with the collective modes of expression that they give rise to can be inferred from single-cell gene expression measurements, opening up a wider range of biological questions that can be addressed using sequencing-based modalities.

A unified model for interpretable latent embedding of multi-sample, multi-condition single-cell data

Single-cell analysis across multiple samples and conditions requires quantitative modeling of the interplay between the continuum of cell states and the technical and biological sources of sample-to-sample variability. We introduce GEDI, a generative model that identifies latent space variations in multi-sample, multi-condition single-cell datasets and attributes them to sample-level covariates. GEDI enables cross-sample cell state mapping on par with state-of-the-art integration methods, cluster-free differential gene expression analysis along the continuum of cell states, and machine learning-based prediction of sample characteristics from single-cell data. GEDI can also incorporate gene-level prior knowledge to infer pathway and regulatory network activities in single cells. Finally, GEDI extends all these concepts to previously unexplored modalities that require joint consideration of dual measurements, such as the joint analysis of exon inclusion/exclusion reads to model alternative cassette exon splicing, or spliced/unspliced reads to model the mRNA stability landscapes of single cells.

Functional specialisation and coordination of myonuclei

Myofibres serve as the functional unit for locomotion, with the sarcomere as fundamental subunit. Running the entire length of this structure are hundreds of myonuclei, located at the periphery of the myofibre, juxtaposed to the plasma membrane. Myonuclear specialisation and clustering at the centre and ends of the fibre are known to be essential for muscle contraction, yet the molecular basis of this regionalisation has remained unclear. While the 'myonuclear domain hypothesis' helped explain how myonuclei can independently govern large cytoplasmic territories, novel technologies have provided granularity on the diverse transcriptional programs running simultaneously within the syncytia and added a new perspective on how myonuclei communicate. Building upon this, we explore the critical cellular and molecular sources of transcriptional and functional heterogeneity within myofibres, discussing the impact of intrinsic and extrinsic factors on myonuclear programs. This knowledge provides new insights for understanding muscle development, repair, and disease, but also opens avenues for the development of novel and precise therapeutic approaches.

GSK3 regulation Wnt/β-catenin signaling affects adipogenesis in bovine skeletal muscle fibro/adipogenic progenitors

Clarifying the cellular origin and regulatory mechanisms of intramuscular fat (IMF) deposition is crucial for improving beef quality. Here, we used single-nucleus RNA sequencing to analyze the structure and heterogeneity of skeletal muscle cell populations in different developmental stages of Yanbian cattle and identified eight cell types in two developmental stages of calves and adults. Among them, fibro/adipogenic progenitors (FAPs) expressing CD29 (ITGA7)pos and CD56 (NCAM1)neg surface markers were committed to IMF deposition in beef cattle and expressed major Wnt ligands and receptors. LY2090314/XAV-939 was used to activate/inhibit Wnt/β-catenin signal. The results showed that the blockade of Glycogen Synthase Kinase 3 (GSK3) by LY2090314 promoted the stabilization of β-catenin and reduced the expression of genes related adipogenic differentiation (e.g., PPARγ and C/EBPα) in bovine FAPs, confirming the anti-adipogenic effect of GSK3. XAV-939 inhibition of the Wnt/β-catenin pathway promoted the lipid accumulation capacity of FAPs. Furthermore, we found that blocking GSK3 enhanced the paracrine effects of FAPs-MuSCs and increased myotube formation in muscle satellite cells (MuSCs). Overall, our results outline a single-cell atlas of skeletal muscle development in Yanbian cattle, revealed the role of Wnt/GSK3/β-catenin signaling in FAPs adipogenesis, and provide a theoretical basis for further regulation of bovine IMF deposition.

Toward universal cell embeddings: integrating single-cell RNA-seq datasets across species with SATURN

Analysis of single-cell datasets generated from diverse organisms offers unprecedented opportunities to unravel fundamental evolutionary processes of conservation and diversification of cell types. However, interspecies genomic differences limit the joint analysis of cross-species datasets to homologous genes. Here we present SATURN, a deep learning method for learning universal cell embeddings that encodes genes' biological properties using protein language models. By coupling protein embeddings from language models with RNA expression, SATURN integrates datasets profiled from different species regardless of their genomic similarity. SATURN can detect functionally related genes coexpressed across species, redefining differential expression for cross-species analysis. Applying SATURN to three species whole-organism atlases and frog and zebrafish embryogenesis datasets, we show that SATURN can effectively transfer annotations across species, even when they are evolutionarily remote. We also demonstrate that SATURN can be used to find potentially divergent gene functions between glaucoma-associated genes in humans and four other species.

Endothelial cell diversity: the many facets of the crystal

Endothelial cells (ECs) form the inner lining of blood vessels and play crucial roles in angiogenesis. While it has been known for a long time that there are considerable differences among ECs from lymphatic and blood vessels, as well as among arteries, veins and capillaries, the full repertoire of endothelial diversity is only beginning to be elucidated. It has become apparent that the role of ECs is not just limited to their exchange functions. Indeed, a multitude of organ-specific functions, including release of growth factors, regulation of immune functions, have been linked to ECs. Recent years have seen a surge into the identification of spatiotemporal molecular and functional heterogeneity of ECs, supported by technologies such as single-cell RNA sequencing (scRNA-seq), lineage tracing and intersectional genetics. Together, these techniques have spurred the generation of epigenomic, transcriptomic and proteomic signatures of ECs. It is now clear that ECs across organs and in different vascular beds, but even within the same vessel, have unique molecular identities and employ specialized molecular mechanisms to fulfil highly specialized needs. Here, we focus on the molecular heterogeneity of the endothelium in different organs and pathological conditions.

The anion exchanger slc26a3 regulates colonic mucus expansion during steady state and in response to prostaglandin E2, while Cftr regulates de novo mucus release in response to carbamylcholine

The intestinal epithelium is covered by mucus that protects the tissue from the luminal content. Studies have shown that anion secretion via the cystic fibrosis conductance regulator (Cftr) regulates mucus formation in the small intestine. However, mechanisms regulating mucus formation in the colon are less understood. The aim of this study was to explore the role of anion transport in the regulation of mucus formation during steady state and in response to carbamylcholine (CCh) and prostaglandin E2 (PGE2). The broad-spectrum anion transport inhibitor 4,4'-diisothiocyanatostilbene-2,2'-disulfonate (DIDS), CftrdF508 (CF) mice, and the slc26a3 inhibitor SLC26A3-IN-2 were used to inhibit anion transport. In the distal colon, steady-state mucus expansion was reduced by SLC26A3-IN-2 and normal in CF mice. PGE2 stimulated mucus expansion without de novo mucus release in wild type (WT) and CF colon via slc26a3 sensitive mechanisms, while CCh induced de novo mucus secretion in WT but not in CF colon. However, when added simultaneously, CCh and PGE2 stimulated de novo mucus secretion in the CF colon via DIDS-sensitive pathways. A similar response was observed in CF ileum that responded to CCh and PGE2 with DIDS-sensitive de novo mucus secretion. In conclusion, this study suggests that slc26a3 regulates colonic mucus expansion, while Cftr regulates CCh-induced de novo mucus secretion from ileal and distal colon crypts. Furthermore, these findings demonstrate that in the absence of a functional Cftr channel, parallel stimulation with CCh and PGE2 activates additional anion transport processes that help release mucus from intestinal goblet cells.

In vivo interaction screening reveals liver-derived constraints to metastasis

It is estimated that only 0.02% of disseminated tumour cells are able to seed overt metastases1. While this suggests the presence of environmental constraints to metastatic seeding, the landscape of host factors controlling this process remains largely unclear. Here, combining transposon technology2 and fluorescence niche labelling3, we developed an in vivo CRISPR activation screen to systematically investigate the interactions between hepatocytes and metastatic cells. We identify plexin B2 as a critical host-derived regulator of liver colonization in colorectal and pancreatic cancer and melanoma syngeneic mouse models. We dissect a mechanism through which plexin B2 interacts with class IV semaphorins on tumour cells, leading to KLF4 upregulation and thereby promoting the acquisition of epithelial traits. Our results highlight the essential role of signals from the liver parenchyma for the seeding of disseminated tumour cells before the establishment of a growth-promoting niche. Our findings further suggest that epithelialization is required for the adaptation of CRC metastases to their new tissue environment. Blocking the plexin-B2-semaphorin axis abolishes metastatic colonization of the liver and therefore represents a therapeutic strategy for the prevention of hepatic metastases. Finally, our screening approach, which evaluates host-derived extrinsic signals rather than tumour-intrinsic factors for their ability to promote metastatic seeding, is broadly applicable and lays a framework for the screening of environmental constraints to metastasis in other organs and cancer types.

Data acquisition approaches for single cell proteomics

Single-cell proteomics (SCP) aims to characterize the proteome of individual cells, providing insights into complex biological systems. It reveals subtle differences in distinct cellular populations that bulk proteome analysis may overlook, which is essential for understanding disease mechanisms and developing targeted therapies. Mass spectrometry (MS) methods in SCP allow the identification and quantification of thousands of proteins from individual cells. Two major challenges in SCP are the limited material in single-cell samples necessitating highly sensitive analytical techniques and the efficient processing of samples, as each biological sample requires thousands of single cell measurements. This review discusses MS advancements to mitigate these challenges using data-dependent acquisition (DDA) and data-independent acquisition (DIA). Additionally, we examine the use of short liquid chromatography gradients and sample multiplexing methods that increase the sample throughput and scalability of SCP experiments. We believe these methods will pave the way for improving our understanding of cellular heterogeneity and its implications for systems biology.

Intestinal retinol saturase is implicated in the development of obesity and epithelial homeostasis upon injury

Retinol saturase (RetSat) is an oxidoreductase involved in lipid metabolism and the cellular sensitivity to peroxides. RetSat is highly expressed in metabolic organs like the liver and adipose tissue and its global loss in mice increases body weight and adiposity. The regulation of RetSat expression and its function in the intestine are unexplored. Here, we show that RetSat is present in different segments of the digestive system, localizes to intestinal epithelial cells, and is upregulated by feeding mice high-fat diet (HFD). Intestine-specific RetSat deletion in adult mice did not affect nutrient absorption and energy homeostasis basally, but lowered body weight gain and fat mass of HFD-fed mice, potentially via increasing locomotor activity. Moreover, jejunal expression of genes related to β-oxidation and cholesterol efflux was decreased, and colonic cholesterol content was reduced upon RetSat deletion. In colitis, which we show to downregulate intestinal RetSat expression in humans and mice, RetSat ablation improved epithelial architecture of the murine colon. Thus, intestinal RetSat expression is regulated by dietary interventions and inflammation, and its loss reduces weight gain upon HFD feeding and alleviates epithelial damage upon injury.NEW & NOTEWORTHY Retinol saturase (RetSat) is an oxidoreductase with unknown function in the intestine. We found that RetSat localizes in intestinal epithelial cells and that its deletion reduced weight gain and fat mass in obese mice. In colitis, which decreased intestinal RetSat expression in humans and mice, RetSat ablation improved the epithelial architecture of the murine colon, presumably by decreasing ROS production, thus rendering RetSat a novel target for metabolic and inflammatory bowel disease.

Single-cell atlas of the human brain vasculature across development, adulthood and disease

A broad range of brain pathologies critically relies on the vasculature, and cerebrovascular disease is a leading cause of death worldwide. However, the cellular and molecular architecture of the human brain vasculature remains incompletely understood1. Here we performed single-cell RNA sequencing analysis of 606,380 freshly isolated endothelial cells, perivascular cells and other tissue-derived cells from 117 samples, from 68 human fetuses and adult patients to construct a molecular atlas of the developing fetal, adult control and diseased human brain vasculature. We identify extensive molecular heterogeneity of the vasculature of healthy fetal and adult human brains and across five vascular-dependent central nervous system (CNS) pathologies, including brain tumours and brain vascular malformations. We identify alteration of arteriovenous differentiation and reactivated fetal as well as conserved dysregulated genes and pathways in the diseased vasculature. Pathological endothelial cells display a loss of CNS-specific properties and reveal an upregulation of MHC class II molecules, indicating atypical features of CNS endothelial cells. Cell-cell interaction analyses predict substantial endothelial-to-perivascular cell ligand-receptor cross-talk, including immune-related and angiogenic pathways, thereby revealing a central role for the endothelium within brain neurovascular unit signalling networks. Our single-cell brain atlas provides insights into the molecular architecture and heterogeneity of the developing, adult/control and diseased human brain vasculature and serves as a powerful reference for future studies.

Histone β-hydroxybutyrylation is critical in reversal of sarcopenia

Sarcopenia, a leading cause for global disability and mortality, is an age-related muscular disorder, characterized by accelerated muscle mass loss and functional decline. It is known that caloric restriction (CR), ketogenic diet or endurance exercise lessen sarcopenia and elevate circulating β-hydroxybutyrate (β-HB) levels. Whether the elevated β-HB is essential to the reversal of sarcopenia, however, remains unclear. Here we show in both Caenorhabditis elegans and mouse models that an increase of β-HB reverse myofiber atrophy and improves motor functions at advanced ages. β-HB-induced histone lysine β-hydroxybutyrylation (Kbhb) is indispensable for the reversal of sarcopenia. Histone Kbhb enhances transcription of genes associated with mitochondrial pathways, including oxidative phosphorylation, ATP metabolic process and aerobic respiration. This ultimately leads to improve mitochondrial integrity and enhance mitochondrial respiration. The histone Kbhb are validated in mouse model with CR. Thus, we demonstrate that β-HB induces histone Kbhb, increases mitochondrial function, and reverses sarcopenia.

Creation of a novel zebrafish model with low DHA status to study the role of maternal nutrition during neurodevelopment

Docosahexaenoic acid (DHA), a dietary omega-3 fatty acid, is a major building block of brain cell membranes. Offspring rely on maternal DHA transfer to meet their neurodevelopmental needs, but DHA sources are lacking in the American diet. Low DHA status is linked to altered immune responses, white matter defects, impaired vision, and an increased risk of psychiatric disorders during development. However, the underlying cellular mechanisms involved are largely unknown, and advancements in the field have been limited by the existing tools and animal models. Zebrafish are an excellent model for studying neurodevelopmental mechanisms. Embryos undergo rapid external development and are optically transparent, enabling direct observation of individual cells and dynamic cell-cell interactions in a way that is not possible in rodents. Here, we create a novel DHA-deficient zebrafish model by 1) disrupting elovl2, a key gene in the DHA biosynthesis pathway, via CRISPR-Cas9 genome editing, and 2) feeding mothers a DHA-deficient diet. We show that low DHA status during development is associated with a small eye morphological phenotype and demonstrate that even the morphologically normal siblings exhibit dysregulated gene pathways related to vision and stress response. Future work using our zebrafish model could reveal the cellular and molecular mechanisms by which low DHA status leads to neurodevelopmental abnormalities and provide insight into maternal nutritional strategies that optimize infant brain health.

Fetal hypoplastic lungs have multilineage inflammation that is reversed by amniotic fluid stem cell extracellular vesicle treatment

Antenatal administration of extracellular vesicles from amniotic fluid stem cells (AFSC-EVs) reverses features of pulmonary hypoplasia in models of congenital diaphragmatic hernia (CDH). However, it remains unknown which lung cellular compartments and biological pathways are affected by AFSC-EV therapy. Herein, we conducted single-nucleus RNA sequencing (snRNA-seq) on rat fetal CDH lungs treated with vehicle or AFSC-EVs. We identified that intra-amniotically injected AFSC-EVs reach the fetal lung in rats with CDH, where they promote lung branching morphogenesis and epithelial cell differentiation. Moreover, snRNA-seq revealed that rat fetal CDH lungs have a multilineage inflammatory signature with macrophage enrichment, which is reversed by AFSC-EV treatment. Macrophage enrichment in CDH fetal rat lungs was confirmed by immunofluorescence, flow cytometry, and inhibition studies with GW2580. Moreover, we validated macrophage enrichment in human fetal CDH lung autopsy samples. Together, this study advances knowledge on the pathogenesis of pulmonary hypoplasia and further evidence on the value of an EV-based therapy for CDH fetuses.

Detecting novel cell type in single-cell chromatin accessibility data via open-set domain adaptation

Recent advances in single-cell technologies enable the rapid growth of multi-omics data. Cell type annotation is one common task in analyzing single-cell data. It is a challenge that some cell types in the testing set are not present in the training set (i.e. unknown cell types). Most scATAC-seq cell type annotation methods generally assign each cell in the testing set to one known type in the training set but neglect unknown cell types. Here, we present OVAAnno, an automatic cell types annotation method which utilizes open-set domain adaptation to detect unknown cell types in scATAC-seq data. Comprehensive experiments show that OVAAnno successfully identifies known and unknown cell types. Further experiments demonstrate that OVAAnno also performs well on scRNA-seq data. Our codes are available online at https://github.com/lisaber/OVAAnno/tree/master.

A comparison of scRNA-seq annotation methods based on experimentally labeled immune cell subtype dataset

Cell-type annotation is a critical step in single-cell data analysis. With the development of numerous cell annotation methods, it is necessary to evaluate these methods to help researchers use them effectively. Reference datasets are essential for evaluation, but currently, the cell labels of reference datasets mainly come from computational methods, which may have computational biases and may not reflect the actual cell-type outcomes. This study first constructed an experimentally labeled immune cell-subtype single-cell dataset of the same batch and systematically evaluated 18 cell annotation methods. We assessed those methods under five scenarios, including intra-dataset validation, immune cell-subtype validation, unsupervised clustering, inter-dataset annotation, and unknown cell-type prediction. Accuracy and ARI were evaluation metrics. The results showed that SVM, scBERT, and scDeepSort were the best-performing supervised methods. Seurat was the best-performing unsupervised clustering method, but it couldn't fully fit the actual cell-type distribution. Our results indicated that experimentally labeled immune cell-subtype datasets revealed the deficiencies of unsupervised clustering methods and provided new dataset support for supervised methods.

Molecular aspects of the exercise response and training adaptation in skeletal muscle

Skeletal muscle plasticity enables an enormous potential to adapt to various internal and external stimuli and perturbations. Most notably, changes in contractile activity evoke a massive remodeling of biochemical, metabolic and force-generating properties. In recent years, a large number of signals, sensors, regulators and effectors have been implicated in these adaptive processes. Nevertheless, our understanding of the molecular underpinnings of training adaptation remains rudimentary. Specifically, the mechanisms that underlie signal integration, output coordination, functional redundancy and other complex traits of muscle adaptation are unknown. In fact, it is even unclear how stimulus-dependent specification is brought about in endurance or resistance exercise. In this review, we will provide an overview on the events that describe the acute perturbations in single endurance and resistance exercise bouts. Furthermore, we will provide insights into the molecular principles of long-term training adaptation. Finally, current gaps in knowledge will be identified, and strategies for a multi-omic and -cellular analyses of the molecular mechanisms of skeletal muscle plasticity that are engaged in individual, acute exercise bouts and chronic training adaptation discussed.

PP2A-B55 phosphatase counteracts Ki-67-dependent chromosome individualization during mitosis

Cell cycle progression is regulated by the orderly balance between kinase and phosphatase activities. PP2A phosphatase holoenzymes containing the B55 family of regulatory B subunits function as major CDK1-counteracting phosphatases during mitotic exit in mammals. However, the identification of the specific mitotic roles of these PP2A-B55 complexes has been hindered by the existence of multiple B55 isoforms. Here, through the generation of loss-of-function genetic mouse models for the two ubiquitous B55 isoforms (B55α and B55δ), we report that PP2A-B55α and PP2A-B55δ complexes display overlapping roles in controlling the dynamics of proper chromosome individualization and clustering during mitosis. In the absence of PP2A-B55 activity, mitotic cells display increased chromosome individualization in the presence of enhanced phosphorylation and perichromosomal loading of Ki-67. These data provide experimental evidence for a regulatory mechanism by which the balance between kinase and PP2A-B55 phosphatase activity controls the Ki-67-mediated spatial organization of the mass of chromosomes during mitosis.

Insulin and leptin oscillations license food-entrained browning and metabolic flexibility

Timed feeding drives adipose browning, although the integrative mechanisms for the same remain unclear. Here, we show that twice-a-night (TAN) feeding generates biphasic oscillations of circulating insulin and leptin, representing their entrainment by timed feeding. Insulin and leptin surges lead to marked cellular, functional, and metabolic remodeling of subcutaneous white adipose tissue (sWAT), resulting in increased energy expenditure. Single-cell RNA-sequencing (scRNA-seq) analyses and flow cytometry demonstrate a role for insulin and leptin surges in innate lymphoid type 2 (ILC2) cell recruitment and sWAT browning, since sWAT depot denervation or loss of leptin or insulin receptor signaling or ILC2 recruitment each dampens TAN feeding-induced sWAT remodeling and energy expenditure. Consistently, recreating insulin and leptin oscillations via once-a-day timed co-injections is sufficient to favorably remodel innervated sWAT. Innervation is necessary for sWAT remodeling, since denervation of sWAT, but not brown adipose tissue (BAT), blocks TAN-induced sWAT remodeling and resolution of inflammation. In sum, reorganization of nutrient-sensitive pathways remodels sWAT and drives the metabolic benefits of timed feeding.

Cell cycle expression heterogeneity predicts degree of differentiation

Methods that predict fate potential or degree of differentiation from transcriptomic data have identified rare progenitor populations and uncovered developmental regulatory mechanisms. However, some state-of-the-art methods are too computationally burdensome for emerging large-scale data and all methods make inaccurate predictions in certain biological systems. We developed a method in R (stemFinder) that predicts single cell differentiation time based on heterogeneity in cell cycle gene expression. Our method is computationally tractable and is as good as or superior to competitors. As part of our benchmarking, we implemented four different performance metrics to assist potential users in selecting the tool that is most apt for their application. Finally, we explore the relationship between differentiation time and cell fate potential by analyzing a lineage tracing dataset with clonally labelled hematopoietic cells, revealing that metrics of differentiation time are correlated with the number of downstream lineages.

CAbiNet: joint clustering and visualization of cells and genes for single-cell transcriptomics

A fundamental analysis task for single-cell transcriptomics data is clustering with subsequent visualization of cell clusters. The genes responsible for the clustering are only inferred in a subsequent step. Clustering cells and genes together would be the remit of biclustering algorithms, which are often bogged down by the size of single-cell data. Here we present 'Correspondence Analysis based Biclustering on Networks' (CAbiNet) for joint clustering and visualization of single-cell RNA-sequencing data. CAbiNet performs efficient co-clustering of cells and their respective marker genes and jointly visualizes the biclusters in a non-linear embedding for easy and interactive visual exploration of the data.

MuSK-BMP signaling in adult muscle stem cells maintains quiescence and regulates myofiber size

A central question in adult stem cell biology is elucidating the signaling pathways regulating their dynamics and function in diverse physiological and age-related contexts. Muscle stem cells in adults (Satellite Cells; SCs) are generally quiescent but can activate and contribute to muscle repair and growth. Here we tested the role of the MuSK-BMP pathway in regulating adult SC quiescence by deletion of the BMP-binding MuSK Ig3 domain ('ΔIg3-MuSK'). At 3 months of age SC and myonuclei numbers and myofiber size were comparable to WT. However, at 5 months of age SC density was decreased while myofiber size, myonuclear number and grip strength were increased - indicating that SCs had activated and productively fused into the myofibers over this interval. Transcriptomic analysis showed that SCs from uninjured ΔIg3-MuSK mice exhibit signatures of activation. Regeneration experiments showed that ΔIg3-MuSK SCs maintain full stem cell function. Expression of ΔIg3-MuSK in adult SCs was sufficient to break quiescence and increase myofiber size. We conclude that the MuSK-BMP pathway regulates SC quiescence and myofiber size in a cell autonomous, age-dependent manner. Targeting MuSK-BMP signaling in muscle stem cells thus emerges a therapeutic strategy for promoting muscle growth and function in the settings of injury, disease, and aging.

Highlights

MuSK, in its role as a BMP co-receptor, regulates adult muscle stem cell quiescenceThe MuSK-BMP pathway acts cell autonomouslyIncreased muscle size and function with preservation of myonuclear density and stemness in mice with attenuated MuSK-BMP signaling.

Deciphering the impact of aging on splenic endothelial cell heterogeneity and immunosenescence through single-cell RNA sequencing analysis

BACKGROUND

Aging is associated with significant structural and functional changes in the spleen, leading to immunosenescence, yet the detailed effects on splenic vascular endothelial cells (ECs) and their immunomodulatory roles are not fully understood. In this study, a single-cell RNA (scRNA) atlas of EC transcriptomes from young and aged mouse spleens was constructed to reveal age-related molecular changes, including increased inflammation and reduced vascular development and also the potential interaction between splenic endothelial cells and immune cells.

RESULTS

Ten clusters of splenic endothelial cells were identified. DEGs analysis across different EC clusters revealed the molecular changes with aging, showing the increase in the overall inflammatory microenvironment and the loss in vascular development function of aged ECs. Notably, four EC clusters with immunological functions were identified, suggesting an Endothelial-to-Immune-like Cell Transition (EndICLT) potentially driven by aging. Pseudotime analysis of the Immunology4 cluster further indicated a possible aging-induced transitional state, potentially initiated by Ctss gene activation. Finally, the effects of aging on cell signaling communication between different EC clusters and immune cells were analyzed.

CONCLUSIONS

This comprehensive atlas elucidates the complex interplay between ECs and immune cells in the aging spleen, offering new insights into endothelial heterogeneity, reprogramming, and the mechanisms of immunosenescence.

Macrophage iron dyshomeostasis promotes aging-related renal fibrosis

Renal aging, marked by the accumulation of senescent cells and chronic low-grade inflammation, leads to renal interstitial fibrosis and impaired function. In this study, we investigate the role of macrophages, a key regulator of inflammation, in renal aging by analyzing kidney single-cell RNA sequencing data of C57BL/6J mice from 8 weeks to 24 months. Our findings elucidate the dynamic changes in the proportion of kidney cell types during renal aging and reveal that increased macrophage infiltration contributes to chronic low-grade inflammation, with these macrophages exhibiting senescence and activation of ferroptosis signaling. CellChat analysis indicates enhanced communications between macrophages and tubular cells during aging. Suppressing ferroptosis alleviates macrophage-mediated tubular partial epithelial-mesenchymal transition in vitro, thereby mitigating the expression of fibrosis-related genes. Using SCENIC analysis, we infer Stat1 as a key age-related transcription factor promoting iron dyshomeostasis and ferroptosis in macrophages by regulating the expression of Pcbp1, an iron chaperone protein that inhibits ferroptosis. Furthermore, through virtual screening and molecular docking from a library of anti-aging compounds, we construct a docking model targeting Pcbp1, which indicates that the natural small molecule compound Rutin can suppress macrophage senescence and ferroptosis by preserving Pcbp1. In summary, our study underscores the crucial role of macrophage iron dyshomeostasis and ferroptosis in renal aging. Our results also suggest Pcbp1 as an intervention target in aging-related renal fibrosis and highlight Rutin as a potential therapeutic agent in mitigating age-related renal chronic low-grade inflammation and fibrosis.

Key factor screening in mouse NASH model using single-cell sequencing combined with machine learning

Aims

To identify and analyze genes closely related to the progression of nonalcoholic steatohepatitis (NASH) by employing a combination of single-cell RNA sequencing and machine-learning algorithms.

Main methods

Single-cell RNA sequencing (scRNA-seq) analysis was performed to find the cell population with the most significant differences between the Chow and NASH groups. This approach was used to validate the developmental trajectory of this cell population and investigate changes in cellular communication and important signaling pathways among these cells. Subsequently, high dimensional Weighted Gene Co-expression Network Analysis (hdWGCNA) was used to find the key modules in NASH. Machine learning analyses were performed to further identify core genes. Deep learning techniques were applied to elucidate the correlation between core genes and immune cells. The accuracy of this correlation was further confirmed using deep learning techniques, specifically Convolutional Neural Networks.

Key findings

By comparing scRNA-seq data between the Chow and NASH groups, we have observed a notable distinction existing in the Kupffer cell population. Signaling interactions between hepatic macrophages and other cells were significantly heightened in the NASH group. Through subsequent analysis of macrophage subtypes and key modules, we identified 150 genes tightly associated with NASH. Finally, we highlighted the 16 most significant core genes using multiple iterations of machine learning. Furthermore, we pointed out the close relationship between core genes and immune cells.

Significances

Using scRNA-seq analysis and machine learning, we can distinguish NASH-related genes from large genetic datasets, providing theoretical support in finding potential targets for the development of novel therapies.

A versatile functional interaction between electrically silent KV subunits and KV7 potassium channels

Voltage-gated K+ (KV) channels govern K+ ion flux across cell membranes in response to changes in membrane potential. They are formed by the assembly of four subunits, typically from the same family. Electrically silent KV channels (KVS), however, are unable to conduct currents on their own. It has been assumed that these KVS must obligatorily assemble with subunits from the KV2 family into heterotetrameric channels, thereby giving rise to currents distinct from those of homomeric KV2 channels. Herein, we show that KVS subunits indeed also modulate the activity, biophysical properties and surface expression of recombinant KV7 isoforms in a subunit-specific manner. Employing co-immunoprecipitation, and proximity labelling, we unveil the spatial coexistence of KVS and KV7 within a single protein complex. Electrophysiological experiments further indicate functional interaction and probably heterotetramer formation. Finally, single-cell transcriptomic analyses identify native cell types in which this KVS and KV7 interaction may occur. Our findings demonstrate that KV cross-family interaction is much more versatile than previously thought-possibly serving nature to shape potassium conductance to the needs of individual cell types.

Astrocytic Neuroligin-3 influences gene expression and social behavior, but is dispensable for synapse number

Neuroligin-3 (Nlgn3) is an autism-associated cell-adhesion molecule that interacts with neurexins and is robustly expressed in both neurons and astrocytes. Neuronal Nlgn3 is an essential regulator of synaptic transmission but the function of astrocytic Nlgn3 is largely unknown. Given the high penetrance of Nlgn3 mutations in autism and the emerging role of astrocytes in neuropsychiatric disorders, we here asked whether astrocytic Nlgn3 might shape neural circuit properties in the cerebellum similar to neuronal Nlgn3. Imaging of tagged Nlgn3 protein produced by CRISPR/Cas9-mediated genome editing showed that Nlgn3 is enriched in the cell body but not the fine processes of cerebellar astrocytes (Bergmann glia). Astrocyte-specific knockout of Nlgn3 did not detectably alter the number of synapses, synaptic transmission, or astrocyte morphology in mouse cerebellum. However, spatial transcriptomic analyses revealed a significant shift in gene expression among multiple cerebellar cell types after the deletion of astrocytic Nlgn3. Hence, in contrast to neuronal Nlgn3, astrocytic Nlgn3 in the cerebellum is not involved in shaping synapses but may modulate gene expression in specific brain areas.

Growth Hormone Receptor Antagonist Markedly Improves Gemcitabine Response in a Mouse Xenograft Model of Human Pancreatic Cancer

Chemotherapy treatment against pancreatic ductal adenocarcinoma (PDAC) is thwarted by tumoral activation of multiple therapy resistance pathways. The growth hormone (GH)-GH receptor (GHR) pair is a covert driver of multimodal therapy resistance in cancer and is overexpressed in PDAC tumors, yet the therapeutic potential of targeting the same has not been explored. Here, we report that GHR expression is a negative prognostic factor in patients with PDAC. Combinations of gemcitabine with different GHR antagonists (GHRAs) markedly improve therapeutic outcomes in nude mice xenografts. Employing cultured cells, mouse xenografts, and analyses of the human PDAC transcriptome, we identified that attenuation of the multidrug transporter and epithelial-to-mesenchymal transition programs in the tumors underlie the observed augmentation of chemotherapy efficacy by GHRAs. Moreover, in human PDAC patients, GHR expression strongly correlates with a gene signature of tumor promotion and immune evasion, which corroborate with that in syngeneic tumors in wild-type vs. GH transgenic mice. Overall, we found that GH action in PDAC promoted a therapy-refractory gene signature in vivo, which can be effectively attenuated by GHR antagonism. Our results collectively present a proof of concept toward considering GHR antagonists to improve chemotherapeutic outcomes in the highly chemoresistant PDAC.

Site-specific genetic and functional signatures of aortic endothelial cells at aneurysm predilection sites in healthy and AngII ApoE-/- mice

Aortic aneurysm is characterized by a pathological dilation at specific predilection sites of the vessel and potentially results in life-threatening vascular rupture. Herein, we established a modified "Häutchen method" for the local isolation of endothelial cells (ECs) from mouse aorta to analyze their spatial heterogeneity and potential role in site-specific disease development. When we compared ECs from aneurysm predilection sites of healthy mice with adjacent control segments we found regulation of genes related to extracellular matrix remodeling, angiogenesis and inflammation, all pathways playing a critical role in aneurysm development. We also detected enhanced cortical stiffness of the endothelium at these sites. Gene expression of ECs from aneurysms of the AngII ApoE-/- model when compared to sham animals mimicked expression patterns from predilection sites of healthy animals. Thus, this work highlights a striking genetic and functional regional heterogeneity in aortic ECs of healthy mice, which defines the location of aortic aneurysm formation in disease.

DVsc: An Automated Framework for Efficiently Detecting Viral Infection from Single-cell Transcriptomics Data

Single-cell RNA sequencing (scRNA-seq) has emerged as a valuable tool for studying cellular heterogeneity in various fields, particularly in virological research. By studying the viral and cellular transcriptomes, the dynamics of viral infection can be investigated at a single-cell resolution. However, limited studies have been conducted to investigate whether RNA transcripts from clinical samples contain substantial amounts of viral RNAs, and a specific computational framework for efficiently detecting viral reads based on scRNA-seq data has not been developed. Hence, we introduce DVsc, an open-source framework for precise quantitative analysis of viral infection from single-cell transcriptomics data. When applied to approximately 200 diverse clinical samples that were infected by more than 10 different viruses, DVsc demonstrated high accuracy in systematically detecting viral infection across a wide array of cell types. This innovative bioinformatics pipeline could be crucial for addressing the potential effects of surreptitiously invading viruses on certain illnesses, as well as for designing novel medicines to target viruses in specific host cell subsets and evaluating the efficacy of treatment. DVsc supports the FASTQ format as an input and is compatible with multiple single-cell sequencing platforms. Moreover, it could also be applied to sequences from bulk RNA sequencing data. DVsc is available at http://62.234.32.33:5000/DVsc.

Endothelial Zmiz1 modulates physiological and pathophysiological angiogenesis during retinal development

Angiogenesis is a highly coordinated process involving the control of various endothelial cell behaviors. Mechanisms for transcription factor involvement in the regulation of endothelial cell dynamics and angiogenesis have become better understood, however much remains unknown, especially the role of non-DNA binding transcriptional cofactors. Here, we show that Zmiz1, a transcription cofactor, is enriched in the endothelium and critical for embryonic vascular development, postnatal retinal angiogenesis, and pathological angiogenesis in oxygen induced retinopathy (OIR). In mice, endothelial cell-specific deletion of Zmiz1 during embryogenesis led to lethality due to abnormal angiogenesis and vascular defects. Inducible endothelial cell-specific ablation of Zmiz1 postnatally resulted in impaired retinal vascular outgrowth, decreased vascular density, and increased vessel regression. In addition, angiogenic sprouting in the superficial and deep layers of the retina was markedly reduced. Correspondingly, vascular sprouting in fibrin bead assays was significantly reduced in the absence of Zmiz1, while further in vitro and in vivo evidence also suggested deficits in EC migration. In agreement with the defective sprouting angiogenesis phenotype, gene expression analysis of isolated retinal endothelial cells revealed downregulation of tip-cell enriched genes upon inactivation of Zmiz1. Lastly, our study suggested that endothelial Zmiz1 is critical for intraretinal revascularization following hypoxia exposure in the OIR model. Taken together, these findings begin to define the previously unspecified role of endothelial Zmiz1 in physiological and pathological angiogenesis.

Beaconet: A Reference-Free Method for Integrating Multiple Batches of Single-Cell Transcriptomic Data in Original Molecular Space

Integrating multiple single-cell datasets is essential for the comprehensive understanding of cell heterogeneity. Batch effect is the undesired systematic variations among technologies or experimental laboratories that distort biological signals and hinder the integration of single-cell datasets. However, existing methods typically rely on a selected dataset as a reference, leading to inconsistent integration performance using different references, or embed cells into uninterpretable low-dimensional feature space. To overcome these limitations, a reference-free method, Beaconet, for integrating multiple single-cell transcriptomic datasets in original molecular space by aligning the global distribution of each batch using an adversarial correction network is presented. Through extensive comparisons with 13 state-of-the-art methods, it is demonstrated that Beaconet can effectively remove batch effect while preserving biological variations and is superior to existing unsupervised methods using all possible references in overall performance. Furthermore, Beaconet performs integration in the original molecular feature space, enabling the characterization of cell types and downstream differential expression analysis directly using integrated data with gene-expression features. Additionally, when applying to large-scale atlas data integration, Beaconet shows notable advantages in both time- and space-efficiencies. In summary, Beaconet serves as an effective and efficient batch effect removal tool that can facilitate the integration of single-cell datasets in a reference-free and molecular feature-preserved mode.

A single-cell atlas of chromatin accessibility in mouse organogenesis

Organogenesis is a highly complex and precisely regulated process. Here we profiled the chromatin accessibility in >350,000 cells derived from 13 mouse embryos at four developmental stages from embryonic day (E) 10.5 to E13.5 by SPATAC-seq in a single experiment. The resulting atlas revealed the status of 830,873 candidate cis-regulatory elements in 43 major cell types. By integrating the chromatin accessibility atlas with the previous transcriptomic dataset, we characterized cis-regulatory sequences and transcription factors associated with cell fate commitment, such as Nr5a2 in the development of gastrointestinal tract, which was preliminarily supported by the in vivo experiment in zebrafish. Finally, we integrated this atlas with the previous single-cell chromatin accessibility dataset from 13 adult mouse tissues to delineate the developmental stage-specific gene regulatory programmes within and across different cell types and identify potential molecular switches throughout lineage development. This comprehensive dataset provides a foundation for exploring transcriptional regulation in organogenesis.

The marine-derived compound TAG alleviates Parkinson's disease by restoring RUBCN-mediated lipid metabolism homeostasis

Parkinson's disease (PD) is the second most common neurodegenerative disease, and its prevalence is increasing. Currently, no effective therapies for PD exist. Marine-derived natural compounds are considered important resources for the discovery of new drugs due to their distinctive structures and diverse activities. In this study, tetrahydroauroglaucin (TAG), a polyketide isolated from a marine sponge, was found to have notable neuroprotective effects on MPTP/MPP+-induced neurotoxicity. RNA sequencing analysis and metabolomics revealed that TAG significantly improved lipid metabolism disorder in PD models. Further investigation indicated that TAG markedly decreased the accumulation of lipid droplets (LDs), downregulated the expression of RUBCN, and promoted autophagic flux. Moreover, conditional knockdown of Rubcn notably attenuated PD-like symptoms and the accumulation of LDs, accompanied by blockade of the neuroprotective effect of TAG. Collectively, our results first indicated that TAG, a promising PD therapeutic candidate, could suppress the accumulation of LDs through the RUBCN-autophagy pathway, which highlighted a novel and effective strategy for PD treatment.