Dictionary learning for integrative, multimodal and scalable single-cell analysis

Characterization of a Human Respiratory Mucosa Model to Study Odorant Metabolism

Nasal xenobiotic metabolizing enzymes (XMEs) are important for the sense of smell because they influence odorant availability and quality. Since the major part of the human nasal cavity is lined by a respiratory mucosa, we hypothesized that this tissue contributed to nasal odorant metabolism through XME activity. Thus, we built human respiratory tissue models and characterized the XME profiles using single-cell RNA sequencing. We focused on the XMEs dicarbonyl and l-xylulose reductase, aldehyde dehydrogenase (ALDH) 1A1, and ALDH3A1, which play a role in food odorant metabolism. We demonstrated protein abundance and localization in the tissue models and showed the metabolic activity of the corresponding enzyme families by exposing the models to the odorants 3,4-hexandione and benzaldehyde. Using gas chromatography coupled with mass spectrometry, we observed, for example, a significantly higher formation of the corresponding metabolites 4-hydroxy-3-hexanone (39.03 ± 1.5%, p = 0.0022), benzyl alcohol (10.05 ± 0.88%, p = 0.0008), and benzoic acid (8.49 ± 0.57%, p = 0.0004) in odorant-treated tissue models compared to untreated controls (0 ± 0, 0.12 ± 0.12, and 0.18 ± 0.18%, respectively). This is the first study that reveals the XME profile of tissue-engineered human respiratory mucosa models and demonstrates their suitability to study nasal odorant metabolism.

Single-cell and spatial transcriptomics reveal a high glycolysis B cell and tumor-associated macrophages cluster correlated with poor prognosis and exhausted immune microenvironment in diffuse large B-cell lymphoma

BACKGROUND

Diffuse large B-cell lymphoma (DLBCL) is a heterogeneous malignancy characterized by varied responses to treatment and prognoses. Understanding the metabolic characteristics driving DLBCL progression is crucial for developing personalized therapies.

METHODS

This study utilized multiple omics technologies including single-cell transcriptomics (n = 5), bulk transcriptomics (n = 966), spatial transcriptomics (n = 10), immunohistochemistry (n = 34), multiple immunofluorescence (n = 20) and to elucidate the metabolic features of highly malignant DLBCL cells and tumor-associated macrophages (TAMs), along with their associated tumor microenvironment. Metabolic pathway analysis facilitated by scMetabolism, and integrated analysis via hdWGCNA, identified glycolysis genes correlating with malignancy, and the prognostic value of glycolysis genes (STMN1, ENO1, PKM, and CDK1) and TAMs were verified.

RESULTS

High-glycolysis malignant DLBCL tissues exhibited an immunosuppressive microenvironment characterized by abundant IFN_TAMs (CD68+CXCL10+PD-L1+) and diminished CD8+ T cell infiltration. Glycolysis genes were positively correlated with malignancy degree. IFN_TAMs exhibited high glycolysis activity and closely communicating with high-malignancy DLBCL cells identified within datasets. The glycolysis score, evaluated by seven genes, emerged as an independent prognostic factor (HR = 1.796, 95% CI: 1.077-2.995, p = 0.025 and HR = 2.631, 95% CI: 1.207-5.735, p = 0.015) along with IFN_TAMs were positively correlated with poor survival (p < 0.05) in DLBCL. Immunohistochemical validation of glycolysis markers (STMN1, ENO1, PKM, and CDK1) and multiple immunofluorescence validation of IFN_TAMs underscored their prognostic value (p < 0.05) in DLBCL.

CONCLUSIONS

This study underscores the significance of glycolysis in tumor progression and modulation of the immune microenvironment. The identified glycolysis genes and IFN_TAMs represent potential prognostic markers and therapeutic targets in DLBCL.

Intrinsic Gata4 expression sensitizes the aortic root to dilation in a Loeys-Dietz syndrome mouse model

Loeys-Dietz syndrome (LDS) is an aneurysm disorder caused by mutations that decrease transforming growth factor-β (TGF-β) signaling. Although aneurysms develop throughout the arterial tree, the aortic root is a site of heightened risk. To identify molecular determinants of this vulnerability, we investigated the heterogeneity of vascular smooth muscle cells (VSMCs) in the aorta of Tgfbr1 M318R/+ LDS mice by single cell and spatial transcriptomics. Reduced expression of components of the extracellular matrix-receptor apparatus and upregulation of stress and inflammatory pathways were observed in all LDS VSMCs. However, regardless of genotype, a subset of Gata4-expressing VSMCs predominantly located in the aortic root intrinsically displayed a less differentiated, proinflammatory profile. A similar population was also identified among aortic VSMCs in a human scRNAseq dataset. Postnatal VSMC-specific Gata4 deletion reduced aortic root dilation in LDS mice, suggesting that this factor sensitizes the aortic root to the effects of impaired TGF-β signaling.

Single-cell and spatiotemporal profile of ovulation in the mouse ovary

Ovulation is a spatiotemporally coordinated process that involves several tightly controlled events, including oocyte meiotic maturation, cumulus expansion, follicle wall rupture and repair, and ovarian stroma remodeling. To date, no studies have detailed the precise window of ovulation at single-cell resolution. Here, we performed parallel single-cell RNA-seq and spatial transcriptomics on paired mouse ovaries across an ovulation time course to map the spatiotemporal profile of ovarian cell types. We show that major ovarian cell types exhibit time-dependent transcriptional states enriched for distinct functions and have specific localization profiles within the ovary. We also identified gene markers for ovulation-dependent cell states and validated these using orthogonal methods. Finally, we performed cell-cell interaction analyses to identify ligand-receptor pairs that may drive ovulation, revealing previously unappreciated interactions. Taken together, our data provides a rich and comprehensive resource of murine ovulation that can be mined for discovery by the scientific community.

Recent advances in exploring transcriptional regulatory landscape of crops

Crop breeding entails developing and selecting plant varieties with improved agronomic traits. Modern molecular techniques, such as genome editing, enable more efficient manipulation of plant phenotype by altering the expression of particular regulatory or functional genes. Hence, it is essential to thoroughly comprehend the transcriptional regulatory mechanisms that underpin these traits. In the multi-omics era, a large amount of omics data has been generated for diverse crop species, including genomics, epigenomics, transcriptomics, proteomics, and single-cell omics. The abundant data resources and the emergence of advanced computational tools offer unprecedented opportunities for obtaining a holistic view and profound understanding of the regulatory processes linked to desirable traits. This review focuses on integrated network approaches that utilize multi-omics data to investigate gene expression regulation. Various types of regulatory networks and their inference methods are discussed, focusing on recent advancements in crop plants. The integration of multi-omics data has been proven to be crucial for the construction of high-confidence regulatory networks. With the refinement of these methodologies, they will significantly enhance crop breeding efforts and contribute to global food security.

APOE Christchurch enhances a disease-associated microglial response to plaque but suppresses response to tau pathology

Background

Apolipoprotein E ε4 (APOE4) is the strongest genetic risk factor for late-onset Alzheimer's disease (LOAD). A recent case report identified a rare variant in APOE, APOE3-R136S (Christchurch), proposed to confer resistance to autosomal dominant Alzheimer's Disease (AD). However, it remains unclear whether and how this variant exerts its protective effects.

Methods

We introduced the R136S variant into mouse Apoe (ApoeCh) and investigated its effect on the development of AD-related pathology using the 5xFAD model of amyloidosis and the PS19 model of tauopathy. We used immunohistochemical and biochemical analysis along with single-cell spatial transcriptomics and proteomics to explore the impact of the ApoeCh variant on AD pathological development and the brain's response to plaques and tau.

Results

In 5xFAD mice, ApoeCh enhances a Disease-Associated Microglia (DAM) phenotype in microglia surrounding plaques, and reduces plaque load, dystrophic neurites, and plasma neurofilament light chain. By contrast, in PS19 mice, ApoeCh suppresses the microglial and astrocytic responses to tau-laden neurons and does not reduce tau accumulation or phosphorylation, but partially rescues tau-induced synaptic and myelin loss. We compared how microglia responses differ between the two mouse models to elucidate the distinct DAM signatures induced by ApoeCh. We identified upregulation of antigen presentation-related genes in the DAM response in a PS19 compared to a 5xFAD background, suggesting a differential response to amyloid versus tau pathology that is modulated by the presence of ApoeCh.

Conclusions

These findings highlight the ability of the ApoeCh variant to modulate microglial responses based on the type of pathology, enhancing DAM reactivity in amyloid models and dampening neuroinflammation to promote protection in tau models. This suggests that the Christchurch variant's protective effects likely involve multiple mechanisms, including changes in receptor binding and microglial programming.

Changes in the cellular makeup of motor patterning circuits drive courtship song evolution in Drosophila

How evolutionary changes in genes and neurons encode species variation in complex motor behaviors is largely unknown. Here, we develop genetic tools that permit a neural circuit comparison between the model species Drosophila melanogaster and the closely related species D. yakuba, which has undergone a lineage-specific loss of sine song, one of the two major types of male courtship song in Drosophila. Neuroanatomical comparison of song-patterning neurons called TN1 across the phylogeny demonstrates a link between the loss of sine song and a reduction both in the number of TN1 neurons and the neurites supporting the sine circuit connectivity. Optogenetic activation confirms that TN1 neurons in D. yakuba have lost the ability to drive sine song, although they have maintained the ability to drive the singing wing posture. Single-cell transcriptomic comparison shows that D. yakuba specifically lacks a cell type corresponding to TN1A neurons, the TN1 subtype that is essential for sine song. Genetic and developmental manipulation reveals a functional divergence of the sex determination gene doublesex in D. yakuba to reduce TN1 number by promoting apoptosis. Our work illustrates the contribution of motor patterning circuits and cell type changes in behavioral evolution and uncovers the evolutionary lability of sex determination genes to reconfigure the cellular makeup of neural circuits.

Spatial Deconvolution of Cell Types and Cell States at Scale Utilizing TACIT

Identifying cell types and states remains a time-consuming and error-prone challenge for spatial biology. While deep learning is increasingly used, it is difficult to generalize due to variability at the level of cells, neighborhoods, and niches in health and disease. To address this, we developed TACIT, an unsupervised algorithm for cell annotation using predefined signatures that operates without training data, using unbiased thresholding to distinguish positive cells from background, focusing on relevant markers to identify ambiguous cells in multiomic assays. Using five datasets (5,000,000-cells; 51-cell types) from three niches (brain, intestine, gland), TACIT outperformed existing unsupervised methods in accuracy and scalability. Integration of TACIT-identified cell with a novel Shiny app revealed new phenotypes in two inflammatory gland diseases. Finally, using combined spatial transcriptomics and proteomics, we discover under- and overrepresented immune cell types and states in regions of interest, suggesting multimodality is essential for translating spatial biology to clinical applications.

Pan-cancer analysis of telomere maintenance mechanisms

Telomeres, protective caps at chromosome ends, maintain genomic stability and control cell lifespan. Dysregulated telomere maintenance mechanisms (TMMs) are cancer hallmarks, enabling unchecked cell proliferation. We conducted a pan-cancer evaluation of TMM using RNA sequencing data from The Cancer Genome Atlas for 33 different cancer types and analyzed the activities of telomerase-dependent (TEL) and alternative lengthening of telomeres (ALT) TMM pathways in detail. To further characterize the TMM profiles, we categorized the tumors based on their ALT and TEL TMM pathway activities into five major phenotypes: ALT high TEL low, ALT low TEL low, ALT middle TEL middle, ALT high TEL high, and ALT low TEL high. These phenotypes refer to variations in telomere maintenance strategies, shedding light on the heterogeneous nature of telomere regulation in cancer. Moreover, we investigated the clinical implications of TMM phenotypes by examining their associations with clinical characteristics and patient outcomes. Specific TMM profiles were linked to specific survival patterns, emphasizing the potential of TMM profiling as a prognostic indicator and aiding in personalized cancer treatment strategies. Gene ontology analysis of the TMM phenotypes unveiled enriched biological processes associated with cell cycle regulation (both TEL and ALT), DNA replication (TEL), and chromosome dynamics (ALT) showing that telomere maintenance is tightly intertwined with cellular processes governing proliferation and genomic stability. Overall, our study provides an overview of the complexity of transcriptional regulation of telomere maintenance mechanisms in cancer.

Neuropilin-1high monocytes protect against neonatal inflammation

Neonates are susceptible to inflammatory disorders such as necrotizing enterocolitis (NEC) due to their immature immune system. The timely appearance of regulatory immune cells in early life contributes to the control of inflammation in neonates, yet the underlying mechanisms of which remain poorly understood. In this study, we identified a subset of neonatal monocytes characterized by high levels of neuropilin-1 (Nrp1), termed Nrp1high monocytes. Compared with their Nrp1low counterparts, Nrp1high monocytes displayed potent immunosuppressive activity. Nrp1 deficiency in myeloid cells aggravated the severity of NEC, whereas adoptive transfer of Nrp1high monocytes led to remission of NEC. Mechanistic studies showed that Nrp1, by binding to its ligand Sema4a, induced intracellular p38-MAPK/mTOR signaling and activated the transcription factor KLF4. KLF4 transactivated Nos2 and enhanced the production of nitric oxide (NO), a key mediator of immunosuppression in monocytes. These findings reveal an important immunosuppressive axis in neonatal monocytes and provide a potential therapeutic strategy for treating inflammatory disorders in neonates.

Single-cell mapping of cell-type specific chromatin architecture in the central nervous system

Determining how chromatin is structured in the nucleus is critical to studying its role in gene regulation. Recent advances in the analysis of single-cell chromatin architecture have considerably improved our understanding of cell-type-specific chromosome conformation and nuclear architecture. In this review, we discuss the methods used for analysis of 3D chromatin conformation, including sequencing-based methods, imaging-based techniques, and computational approaches. We further review the application of these methods in the study of the role of chromatin topology in neural development and disorders.

Evidences of neurological injury caused by COVID-19 from glioma tissues and glioma organoids

INTRODUCTION

Despite the extensive neurological symptoms induced by COVID-19 and the identification of SARS-CoV-2 in post-mortem brain samples from COVID-19 patients months after death, the precise mechanisms of SARS-CoV-2 invasion into the central nervous system remain unclear due to the lack of research models.

METHODS

We collected glioma tissue samples from glioma patients who had a recent history of COVID-19 and examined the presence of the SARS-CoV-2 spike protein. Subsequently, spatial transcriptomic analyses were conducted on normal brain tissues, glioma tissues, and glioma tissues from glioma patients with recent COVID-19 history. Additionally, single-cell sequencing data from both glioma tissues and glioma organoids were collected and analyzed. Glioma organoids were utilized to evaluate the efficacy of potential COVID-19 blocking agents.

RESULTS

Glioma tissues from glioma patients with recent COVID-19 history exhibited the presence of the SARS-CoV-2 spike protein. Differences between glioma tissues from glioma patients who had a recent history of COVID-19 and healthy brain tissues primarily manifested in neuronal cells. Notably, neuronal cells within glioma tissues of COVID-19 history demonstrated heightened susceptibility to Alzheimer's disease, depression, and synaptic dysfunction, indicative of neuronal aberrations. Expressions of SARS-CoV-2 entry factors were confirmed in both glioma tissues and glioma organoids. Moreover, glioma organoids were susceptible to pseudo-SARS-CoV-2 infection and the infections could be partly blocked by the potential COVID-19 drugs.

CONCLUSIONS

Gliomas had inherent traits that render them susceptible to SARS-CoV-2 infection, leading to their representability of COVID-19 neurological symptoms. This established a biological foundation for the rationality and feasibility of utilization of glioma organoids as research and blocking drug testing model in SARS-CoV-2 infection within the central nervous system.

Nanoscale Methods for Longitudinal Extraction of Intracellular Contents

Longitudinal analysis of intracellular contents including gene and protein expression is crucial for deciphering the fundamentally dynamic nature of cells. This offers invaluable insights into complex tissue composition and behavior, and drives progress in disease diagnosis, biomarker discovery, and drug development. Traditional longitudinal analysis workflows, involving the destruction of cells at various timepoints, limit insights to singular moments and fail to account for cellular heterogeneity. Current non-destructive approaches, like temporal modeling with single-cell ribonucleic acid sequencing (RNA-seq) and live-cell fluorescence imaging, either rely on biological assumptions or possess the risk of cellular perturbation. Recent advances in nanoscale technologies for non-destructive intracellular content extraction offer a promising solution to these challenges. These novel methods work at the nanoscale to non-destructively access cellular membranes and can be broadly classified into three mechanisms: tip-facilitated aspiration, membrane-based, and probe-based methods. This perspective focuses on these emerging nanotechnologies for repeated intracellular content extraction. Their potential in longitudinal analysis is discussed, the critical requirements for effective repeated sampling are addressed, and the suitability of each technique for various applications is explored. Furthermore, unresolved challenges in repeated sampling are highlighted to encourage further research in this growing field.

Aggressive high-grade NF2 mutant meningiomas downregulate oncogenic YAP signaling via the upregulation of VGLL4 and FAT3/4

Meningiomas are the most common primary brain tumors in adults. Although generally benign, a subset of meningiomas is of higher grade, shows aggressive growth behavior and recurs even after multiple surgeries. Around half of all meningiomas harbor inactivating mutations in NF2. While benign low-grade NF2 mutant meningiomas exhibit few genetic events in addition to NF2 inactivation, aggressive high-grade NF2 mutant meningiomas frequently harbor a highly aberrant genome. We and others have previously shown that NF2 inactivation leads to YAP1 activation and that YAP1 acts as the pivotal oncogenic driver in benign NF2 mutant meningiomas. Using bulk and single-cell RNA-Seq data from a large cohort of human meningiomas, we show that aggressive NF2 mutant meningiomas harbor decreased levels YAP1 activity compared to their benign counterparts. Decreased expression levels of YAP target genes are significantly associated with an increased risk of recurrence. We then identify the increased expression of the YAP1 competitor VGLL4 as well as the YAP1 upstream regulators FAT3/4 as a potential mechanism for the downregulation of YAP activity in aggressive NF2 mutant meningiomas. High expression of these genes is significantly associated with an increased risk of recurrence. In vitro, overexpression of VGLL4 resulted in the downregulation of YAP activity in benign NF2 mutant meningioma cells, confirming the direct link between VGLL4 expression and decreased levels of YAP activity observed in aggressive NF2 mutant meningiomas. Our results shed new insight on the biology of benign and aggressive NF2 mutant meningiomas and may have important implications for the efficacy of therapies targeting oncogenic YAP1 activity in NF2 mutant meningiomas.

Uncover spatially informed variations for single-cell spatial transcriptomics with STew

Motivation

The recent spatial transcriptomics (ST) technologies have enabled characterization of gene expression patterns and spatial information, advancing our understanding of cell lineages within diseased tissues. Several analytical approaches have been proposed for ST data, but effectively utilizing spatial information to unveil the shared variation with gene expression remains a challenge.

Results

We introduce STew, a Spatial Transcriptomic multi-viEW representation learning method, to jointly analyze spatial information and gene expression in a scalable manner, followed by a data-driven statistical framework to measure the goodness of model fit. Through benchmarking using human dorsolateral prefrontal cortex and mouse main olfactory bulb data with true manual annotations, STew achieved superior performance in both clustering accuracy and continuity of identified spatial domains compared with other methods. STew is also robust to generate consistent results insensitive to model parameters, including sparsity constraints. We next applied STew to various ST data acquired from 10× Visium, Slide-seqV2, and 10× Xenium, encompassing single-cell and multi-cellular resolution ST technologies, which revealed spatially informed cell type clusters and biologically meaningful axes. In particular, we identified a proinflammatory fibroblast spatial niche using ST data from psoriatic skins. Moreover, STew scales almost linearly with the number of spatial locations, guaranteeing its applicability to datasets with thousands of spatial locations to capture disease-relevant niches in complex tissues.

Availability and implementation

Source code and the R software tool STew are available from github.com/fanzhanglab/STew.

SCIMAP: A Python Toolkit for Integrated Spatial Analysis of Multiplexed Imaging Data

Multiplexed imaging data are revolutionizing our understanding of the composition and organization of tissues and tumors ("Catching up with Multiplexed Tissue Imaging," 2022). A critical aspect of such "tissue profiling" is quantifying the spatial relationships among cells at different scales from the interaction of neighboring cells to recurrent communities of cells of multiple types. This often involves statistical analysis of 107 or more cells in which up to 100 biomolecules (commonly proteins) have been measured. While software tools currently cater to the analysis of spatial transcriptomics data (Liu et al., 2022), there remains a need for toolkits explicitly tailored to the complexities of multiplexed imaging data including the need to seamlessly integrate image visualization with data analysis and exploration. We introduce SCIMAP, a Python package specifically crafted to address these challenges. With SCIMAP, users can efficiently preprocess, analyze, and visualize large datasets, facilitating the exploration of spatial relationships and their statistical significance. SCIMAP's modular design enables the integration of new algorithms, enhancing its capabilities for spatial analysis.

TIGIT+Tfh show poor B-helper function and negatively correlate with SARS-CoV-2 antibody titre

Circulating follicular helper T cells (cTfh) can show phenotypic alterations in disease settings, including in the context of tissue-damaging autoimmune or anti-viral responses. Using severe COVID-19 as a paradigm of immune dysregulation, we have explored how cTfh phenotype relates to the titre and quality of antibody responses. Severe disease was associated with higher titres of neutralising S1 IgG and evidence of increased T cell activation. ICOS, CD38 and HLA-DR expressing cTfh correlated with serum S1 IgG titres and neutralising strength, and interestingly expression of TIGIT by cTfh showed a negative correlation. TIGIT+cTfh expressed increased IFNγ and decreased IL-17 compared to their TIGIT-cTfh counterparts, and showed reduced capacity to help B cells in vitro. Additionally, TIGIT+cTfh expressed lower levels of CD40L than TIGIT-cTfh, providing a potential explanation for their poor B-helper function. These data identify phenotypic changes in polyclonal cTfh that correlate with specific antibody responses and reveal TIGIT as a marker of cTfh with altered function.

Trans-omic profiling uncovers molecular controls of early human cerebral organoid formation

Defining the molecular networks orchestrating human brain formation is crucial for understanding neurodevelopment and neurological disorders. Challenges in acquiring early brain tissue have incentivized the use of three-dimensional human pluripotent stem cell (hPSC)-derived neural organoids to recapitulate neurodevelopment. To elucidate the molecular programs that drive this highly dynamic process, here, we generate a comprehensive trans-omic map of the phosphoproteome, proteome, and transcriptome of the exit of pluripotency and neural differentiation toward human cerebral organoids (hCOs). These data reveal key phospho-signaling events and their convergence on transcriptional factors to regulate hCO formation. Comparative analysis with developing human and mouse embryos demonstrates the fidelity of our hCOs in modeling embryonic brain development. Finally, we demonstrate that biochemical modulation of AKT signaling can control hCO differentiation. Together, our data provide a comprehensive resource to study molecular controls in human embryonic brain development and provide a guide for the future development of hCO differentiation protocols.

A multi-platform analysis of human gingival crevicular fluid reveals ferroptosis as a relevant regulated cell death mechanism during the clinical progression of periodontitis

Ferroptosis is implicated in the pathogenesis of numerous chronic-inflammatory diseases, yet its association with progressive periodontitis remains unexplored. To investigate the involvement and significance of ferroptosis in periodontitis progression, we assessed sixteen periodontitis-diagnosed patients. Disease progression was clinically monitored over twelve weeks via weekly clinical evaluations and gingival crevicular fluid (GCF) collection was performed for further analyses. Clinical metrics, proteomic data, in silico methods, and bioinformatics tools were combined to identify protein profiles linked to periodontitis progression and to explore their potential connection with ferroptosis. Subsequent western blot analyses validated key findings. Finally, a single-cell RNA sequencing (scRNA-seq) dataset (GSE164241) for gingival tissues was analyzed to elucidate cellular dynamics during periodontitis progression. Periodontitis progression was identified as occurring at a faster rate than traditionally thought. GCF samples from progressing and non-progressing periodontal sites showed quantitative and qualitatively distinct proteomic profiles. In addition, specific biological processes and molecular functions during progressive periodontitis were revealed and a set of hub proteins, including SNCA, CA1, HBB, SLC4A1, and ANK1 was strongly associated with the clinical progression status of periodontitis. Moreover, we found specific proteins - drivers or suppressors - associated with ferroptosis (SNCA, FTH1, HSPB1, CD44, and GCLC), revealing the co-occurrence of this specific type of regulated cell death during the clinical progression of periodontitis. Additionally, the integration of quantitative proteomic data with scRNA-seq analysis suggested the susceptibility of fibroblasts to ferroptosis. Our analyses reveal proteins and processes linked to ferroptosis for the first time in periodontal patients, which offer new insights into the molecular mechanisms of progressive periodontal disease. These findings may lead to novel diagnostic and therapeutic strategies.

Spatial analysis of recurrent glioblastoma reveals perivascular niche organization

Tumor evolution is driven by genetic variation; however, it is the tumor microenvironment (TME) that provides the selective pressure contributing to evolution in cancer. Despite high histopathological heterogeneity within glioblastoma (GBM), the most aggressive brain tumor, the interactions between the genetically distinct GBM cells and the surrounding TME are not fully understood. To address this, we analyzed matched primary and recurrent GBM archival tumor tissues with imaging-based techniques aimed to simultaneously evaluate tumor tissues for the presence of hypoxic, angiogenic, and inflammatory niches, extracellular matrix (ECM) organization, TERT promoter mutational status, and several oncogenic amplifications on the same slide and location. We found that the relationships between genetic and TME diversity are different in primary and matched recurrent tumors. Interestingly, the texture of the ECM, identified by label-free reflectance imaging, was predictive of single-cell genetic traits present in the tissue. Moreover, reflectance of ECM revealed structured organization of the perivascular niche in recurrent GBM, enriched in immunosuppressive macrophages. Single-cell spatial transcriptomics further confirmed the presence of the niche-specific macrophage populations and identified interactions between endothelial cells, perivascular fibroblasts, and immunosuppressive macrophages. Our results underscore the importance of GBM tissue organization in tumor evolution and highlight genetic and spatial dependencies.

scPLAN: a hierarchical computational framework for single transcriptomics data annotation, integration and cell-type label refinement

MOTIVATION

In the past decade, single-cell RNA sequencing (scRNA-seq) has emerged as a pivotal method for transcriptomic profiling in biomedical research. Precise cell-type identification is crucial for subsequent analysis of single-cell data. And the integration and refinement of annotated data are essential for building comprehensive databases. However, prevailing annotation techniques often overlook the hierarchical organization of cell types, resulting in inconsistent annotations. Meanwhile, most existing integration approaches fail to integrate datasets with different annotation depths and none of them can enhance the labels of outdated data with lower annotation resolutions using more intricately annotated datasets or novel biological findings.

RESULTS

Here, we introduce scPLAN, a hierarchical computational framework designed for scRNA-seq data analysis. scPLAN excels in annotating unlabeled scRNA-seq data using a reference dataset structured along a hierarchical cell-type tree. It identifies potential novel cell types in a systematic, layer-by-layer manner. Additionally, scPLAN effectively integrates annotated scRNA-seq datasets with varying levels of annotation depth, ensuring consistent refinement of cell-type labels across datasets with lower resolutions. Through extensive annotation and novel cell detection experiments, scPLAN has demonstrated its efficacy. Two case studies have been conducted to showcase how scPLAN integrates datasets with diverse cell-type label resolutions and refine their cell-type labels.

AVAILABILITY

https://github.com/michaelGuo1204/scPLAN.

Single-cell RNA analysis of chromodomain-encoding genes in mesenchymal stromal cells of the mouse dental pulp

The chromodomain helicase DNA-binding (CHD) and chromobox (CBX) families of proteins play crucial roles in cell fate decisions, differentiation, and cell proliferation in a broad variety of tissues and cell types. CHD proteins are ATP-dependent epigenetic enzymes actively engaged in transcriptional regulation, DNA replication, and DNA damage repair, whereas CBX proteins are transcriptional repressors mainly involved in the formation of heterochromatin. The pleiotropic effects of CHD and CBX proteins are largely dependent on their versatility to interact with other key components of the epigenetic and transcriptional machinery. Although the function and regulatory modes of CHD and CBX factors are well established in many cell types, little is known about their roles during osteogenic differentiation. A single-cell RNA-sequencing analysis of the mouse incisor dental pulp revealed distinct spatiotemporal expression patterns of CHD- and CBX-encoding genes within different clusters of mesenchymal stromal cells (MSCs) representing various stages of osteogenic differentiation. Additionally, genes encoding interaction partners of CHD and CBX proteins, such as subunits of the trithorax-COMPASS and polycomb chromatin remodeling complexes, exhibited differential co-expression behaviors within MSC subpopulations. Thus, CHD- and CBX-encoding genes show partially overlapping but distinct expression patterns in MSCs, suggesting their differential roles in osteogenic cell fate decisions.

Multi-omic profiling reveals the endogenous and neoplastic responses to immunotherapies in cutaneous T cell lymphoma

Cutaneous T cell lymphomas (CTCLs) are skin cancers with poor survival rates and limited treatments. While immunotherapies have shown some efficacy, the immunological consequences of administering immune-activating agents to CTCL patients have not been systematically characterized. We apply a suite of high-dimensional technologies to investigate the local, cellular, and systemic responses in CTCL patients receiving either mono- or combination anti-PD-1 plus interferon-gamma (IFN-γ) therapy. Neoplastic T cells display no evidence of activation after immunotherapy. IFN-γ induces muted endogenous immunological responses, while anti-PD-1 elicits broader changes, including increased abundance of CLA+CD39+ T cells. We develop an unbiased multi-omic profiling approach enabling discovery of immune modules stratifying patients. We identify an enrichment of activated regulatory CLA+CD39+ T cells in non-responders and activated cytotoxic CLA+CD39+ T cells in leukemic patients. Our results provide insights into the effects of immunotherapy in CTCL patients and a generalizable framework for multi-omic analysis of clinical trials.

Transcriptomic profiling and regulatory pathways of cardiac resident macrophages in aging

Cardiovascular diseases are an array of age-related disorders, and accumulating evidence suggests a link between cardiac resident macrophages (CRMs) and the age-related disorders. However, how does CRMs alter with aging remains elusive. In the present study, aged mice (20 months old) have been employed to check for their cardiac structural and functional alterations, and the changes in the proportion of CRM subsets as well, followed by sorting of CRMs, including C-C Motif Chemokine Receptor 2 (CCR2)+ and CCR2- CRMs, which were subjected to Smart-Seq. Integrated analysis of the Smart-Seq data with three publicly available single-cell RNA-seq datasets revealed that inflammatory genes were drastic upregulated for both CCR2+ and CCR2- CRMs with aging, but genes germane to wound healing were downregulated for CCR2- CRMs, suggesting the differential functions of these two subsets. More importantly, inflammatory genes involved in damage sensing, complement cascades, and phagocytosis were largely upregulated in CCR2- CRMs, implying the imbalance of inflammatory response upon aging. Our work provides a comprehensive framework and transcriptional resource for assessing the impact of aging on CRMs with a potential for further understanding cardiac aging.

Applications of single-cell RNA sequencing in spermatogenesis and molecular evolution

Spermatogenic cell heterogeneity is determined by the complex process of spermatogenesis differentiation. However, effectively revealing the regulatory mechanisms underlying mammalian spermatogenic cell development and differentiation via traditional methods is difficult. Advances in technology have led to the emergence of many single-cell transcriptome sequencing protocols, which have partially addressed these challenges. In this review, we detail the principles of 10x Genomics technology and summarize the methods for downstream analysis of single-cell transcriptome sequencing data. Furthermore, we explore the role of single-cell transcriptome sequencing in revealing the heterogeneity of testicular ecological niche cells, delineating the establishment and disruption of testicular immune homeostasis during human spermatogenesis, investigating abnormal spermatogenesis in humans, and, ultimately, elucidating the molecular evolution of mammalian spermatogenesis.

Early infiltrating NKT lymphocytes attenuate bone regeneration through secretion of CXCL2

Trauma rapidly mobilizes the immune response of surrounding tissues and activates regeneration program. Manipulating immune response to promote tissue regeneration shows a broad application prospect. However, the understanding of bone healing dynamics at cellular level remains limited. Here, we characterize the landscape of immune cells after alveolar bone injury and reveal a pivotal role of infiltrating natural killer T (NKT) cells. We observe a rapid increase in NKT cells after injury, which inhibit osteogenic differentiation of mesenchymal stem cells (MSCs) and impair alveolar bone healing. Cxcl2 is up-regulated in NKT cells after injury. Systemic administration of CXCL2-neutralizing antibody or genetic deletion of Cxcl2 improves the bone healing process. In addition, we fabricate a gelatin-based porous hydrogel to deliver NK1.1 depletion antibody, which successfully promotes alveolar bone healing. In summary, our study highlights the importance of NKT cells in the early stage of bone healing and provides a potential therapeutic strategy for accelerating bone regeneration.

Age-related and species-specific methylation changes in the protein-coding marmoset sperm epigenome

The sperm epigenome is thought to affect the developmental programming of the resulting embryo, influencing health and disease in later life. Age-related methylation changes in the sperm of old fathers may mediate the increased risks for reproductive and offspring medical problems. The impact of paternal age on sperm methylation has been extensively studied in humans and, to a lesser extent, in rodents and cattle. Here, we performed a comparative analysis of paternal age effects on protein-coding genes in the human and marmoset sperm methylomes. The marmoset has gained growing importance as a non-human primate model of aging and age-related diseases. Using reduced representation bisulfite sequencing, we identified age-related differentially methylated transcription start site (ageTSS) regions in 204 marmoset and 27 human genes. The direction of methylation changes was the opposite, increasing with age in marmosets and decreasing in humans. None of the identified ageTSS was differentially methylated in both species. Although the average methylation levels of all TSS regions were highly correlated between marmosets and humans, with the majority of TSS being hypomethylated in sperm, more than 300 protein-coding genes were endowed with species-specifically (hypo)methylated TSS. Several genes of the glycosphingolipid (GSL) biosynthesis pathway, which plays a role in embryonic stem cell differentiation and regulation of development, were hypomethylated (<5%) in human and fully methylated (>95%) in marmoset sperm. The expression levels and patterns of defined sets of GSL genes differed considerably between human and marmoset pre-implantation embryo stages and blastocyst tissues, respectively.

Tendon-associated gene expression precedes osteogenesis in mid-palatal suture establishment

Orthodontic maxillary expansion relies on intrinsic mid-palatal suture mechanobiology to induce guided osteogenesis, yet establishment of the mid-palatal suture within the continuous secondary palate and causes of maxillary insufficiency remain poorly understood. In contrast, advances in cranial suture research hold promise to improve surgical repair of prematurely fused cranial sutures in craniosynostosis to potentially restore the obliterated signaling environment and ensure continual success of the intervention. We hypothesized that mid-palatal suture establishment is governed by shared principles with calvarial sutures and involves functional linkage between expanding primary ossification centres with the midline mesenchyme. We characterized establishment of the mid-palatal suture from late embryonic to early postnatal timepoints. Suture establishment was visualized using histological techniques and multimodal transcriptomics. We identified that mid-palatal suture formation depends on a spatiotemporally controlled signalling milieu in which tendon-associated genes play a significant role. We mapped relationships between extracellular matrix-encoding gene expression, tenocyte markers, and novel suture patency candidate genes. We identified similar expression patterns in FaceBase-deposited scRNA-seq datasets from cranial sutures. These findings demonstrate shared biological principles for suture establishment, providing further avenues for future development and understanding of maxillofacial interventions.

Enhanced Complement Expression in the Tumor Microenvironment Following Neoadjuvant Therapy: Implications for Immunomodulation and Survival in Pancreatic Ductal Adenocarcinoma

Background

Neoadjuvant therapy (NAT) is increasingly being used for pancreatic ductal adenocarcinoma (PDAC) treatment. However, its specific effects on carcinoma cells and the tumor microenvironment (TME) are not fully understood. This study aims to investigate how NAT differentially impacts PDAC's carcinoma cells and TME.

Methods

Spatial transcriptomics was used to compare gene expression profiles in carcinoma cells and the TME between 23 NAT-treated and 13 NAT-naïve PDAC patients, correlating with their clinicopathologic features. Analysis of an online single-nucleus RNA sequencing (snRNA-seq) dataset was performed for validation of the specific cell types responsible for NAT-induced gene expression alterations.

Results

NAT not only induces apoptosis and inhibits proliferation in carcinoma cells but also significantly remodels the TME. Notably, NAT induces a coordinated upregulation of multiple key complement genes (C3, C1S, C1R, C4B and C7) in the TME, making the complement pathway one of the most significantly affected pathways by NAT. Patients with higher TME complement expression following NAT exhibit improved overall survival. These patients also exhibit increased immunomodulatory and neurotrophic cancer-associated fibroblasts (CAFs); more CD4+ T cells, monocytes, and mast cells; and reduced immune exhaustion gene expression. snRNA-seq analysis demonstrates C3 complement was specifically upregulated in CAFs but not in other stroma cell types.

Conclusions

NAT can enhance complement production and signaling within the TME, which is associated with reduced immunosuppression in PDAC. These findings suggest that local complement dynamics could serve as a novel biomarker for prognosis, evaluating treatment response and resistance, and guiding therapeutic strategies in NAT-treated PDAC patients.

Phenome-wide identification of therapeutic genetic targets, leveraging knowledge graphs, graph neural networks, and UK Biobank data

The ongoing expansion of human genomic datasets propels therapeutic target identification; however, extracting gene-disease associations from gene annotations remains challenging. Here, we introduce Mantis-ML 2.0, a framework integrating AstraZeneca's Biological Insights Knowledge Graph and numerous tabular datasets, to assess gene-disease probabilities throughout the phenome. We use graph neural networks, capturing the graph's holistic structure, and train them on hundreds of balanced datasets via a robust semi-supervised learning framework to provide gene-disease probabilities across the human exome. Mantis-ML 2.0 incorporates natural language processing to automate disease-relevant feature selection for thousands of diseases. The enhanced models demonstrate a 6.9% average classification power boost, achieving a median receiver operating characteristic (ROC) area under curve (AUC) score of 0.90 across 5220 diseases from Human Phenotype Ontology, OpenTargets, and Genomics England. Notably, Mantis-ML 2.0 prioritizes associations from an independent UK Biobank phenome-wide association study (PheWAS), providing a stronger form of triaging and mitigating against underpowered PheWAS associations. Results are exposed through an interactive web resource.

Unbiased metastatic niche-labeling identifies estrogen receptor-positive macrophages as a barrier of T cell infiltration during bone colonization

Microenvironment niches determine cellular fates of metastatic cancer cells. However, robust and unbiased approaches to identify niche components and their molecular profiles are lacking. We established Sortase A-Based Microenvironment Niche Tagging (SAMENT), which selectively labels cells encountered by cancer cells during metastatic colonization. SAMENT was applied to multiple cancer models colonizing the same organ and the same cancer to different organs. Common metastatic niche features include macrophage enrichment and T cell depletion. Macrophage niches are phenotypically diverse between different organs. In bone, macrophages express the estrogen receptor alpha (ERα) and exhibit active ERα signaling in male and female hosts. Conditional knockout of Esr1 in macrophages significantly retarded bone colonization by allowing T cell infiltration. ERα expression was also discovered in human bone metastases of both genders. Collectively, we identified a unique population of ERα+ macrophages in the metastatic niche and functionally tied ERα signaling in macrophages to T cell exclusion during metastatic colonization.

HIGHLIGHTS

SAMENT is a robust metastatic niche-labeling approach amenable to single-cell omics.Metastatic niches are typically enriched with macrophages and depleted of T cells.Direct interaction with cancer cells induces ERα expression in niche macrophages. Knockout of Esr1 in macrophages allows T cell infiltration and retards bone colonization.

The future of rapid and automated single-cell data analysis using reference mapping

As the number of single-cell datasets continues to grow rapidly, workflows that map new data to well-curated reference atlases offer enormous promise for the biological community. In this perspective, we discuss key computational challenges and opportunities for single-cell reference-mapping algorithms. We discuss how mapping algorithms will enable the integration of diverse datasets across disease states, molecular modalities, genetic perturbations, and diverse species and will eventually replace manual and laborious unsupervised clustering pipelines.

Endothelial deletion of p53 generates transitional endothelial cells and improves lung development during neonatal hyperoxia

Bronchopulmonary dysplasia (BPD), a prevalent and chronic lung disease affecting premature newborns, results in vascular rarefaction and alveolar simplification. Although the vasculature has been recognized as a main player in this disease, the recently found capillary heterogeneity and cellular dynamics of endothelial subpopulations in BPD remain unclear. Here, we show Cap2 cells are damaged during neonatal hyperoxic injury, leading to their replacement by Cap1 cells which, in turn, significantly decline. Single-cell RNA-seq identifies the activation of numerous p53 target genes in endothelial cells, including Cdkn1a (p21). While global deletion of p53 results in worsened vasculature, endothelial-specific deletion of p53 reverses the vascular phenotype and improves alveolar simplification during hyperoxia. This recovery is associated with the emergence of a transitional EC state, enriched for oxidative stress response genes and growth factors. These findings implicate the p53 pathway in EC type transition during injury-repair and highlights the endothelial contributions to BPD.

A metabolic atlas of mouse aging

Humans are living longer, but this is accompanied by an increased incidence of age-related chronic diseases. Many of these diseases are influenced by age-associated metabolic dysregulation, but how metabolism changes in multiple organs during aging in males and females is not known. Answering this could reveal new mechanisms of aging and age-targeted therapeutics. In this study, we describe how metabolism changes in 12 organs in male and female mice at 5 different ages. Organs show distinct patterns of metabolic aging that are affected by sex differently. Hydroxyproline shows the most consistent change across the dataset, decreasing with age in 11 out of 12 organs investigated. We also developed a metabolic aging clock that predicts biological age and identified alpha-ketoglutarate, previously shown to extend lifespan in mice, as a key predictor of age. Our results reveal fundamental insights into the aging process and identify new therapeutic targets to maintain organ health.

IL-4Rα signaling promotes barrier-altering oncostatin M and IL-6 production in aspirin-exacerbated respiratory disease

BACKGROUND

Aspirin-exacerbated respiratory disease (AERD) is a severe disease involving dysregulated type 2 inflammation. However, the role other inflammatory pathways play in AERD is poorly understood.

OBJECTIVE

We sought to broadly define the inflammatory milieu of the upper respiratory tract in AERD and to determine the effects of IL-4Rα inhibition on mediators of nasal inflammation.

METHODS

Twenty-two AERD patients treated with dupilumab for 3 months were followed over 3 visits and compared to 10 healthy controls. Nasal fluid was assessed for 45 cytokines and chemokines using Olink Target 48. Blood neutrophils and cultured human mast cells, monocytes/macrophages, and nasal fibroblasts were assessed for response to IL-4/13 stimulation in vitro.

RESULTS

Of the nasal fluid cytokines measured, nearly one third were higher in AERD patients compared to healthy controls, including IL-6 and the IL-6 family-related cytokine oncostatin M (OSM), both of which correlated with nasal albumin levels, a marker of epithelial barrier dysregulation. Dupilumab significantly decreased many nasal mediators, including OSM and IL-6. IL-4 stimulation induced OSM production from mast cells and macrophages but not from neutrophils, and OSM and IL-13 stimulation induced IL-6 production from nasal fibroblasts.

CONCLUSION

In addition to type 2 inflammation, innate and IL-6-related cytokines are also elevated in the respiratory tract in AERD. Both OSM and IL-6 are locally produced in nasal polyps and likely promote pathology by negatively affecting epithelial barrier function. IL-4Rα blockade, although seemingly directed at type 2 inflammation, also decreases mediators of innate inflammation and epithelial dysregulation, which may contribute to dupilumab's therapeutic efficacy in AERD.

Scalable integration of multiomic single-cell data using generative adversarial networks

MOTIVATION

Single-cell profiling has become a common practice to investigate the complexity of tissues, organs, and organisms. Recent technological advances are expanding our capabilities to profile various molecular layers beyond the transcriptome such as, but not limited to, the genome, the epigenome, and the proteome. Depending on the experimental procedure, these data can be obtained from separate assays or the very same cells. Yet, integration of more than two assays is currently not supported by the majority of the computational frameworks avaiable.

RESULTS

We here propose a Multi-Omic data integration framework based on Wasserstein Generative Adversarial Networks suitable for the analysis of paired or unpaired data with a high number of modalities (>2). At the core of our strategy is a single network trained on all modalities together, limiting the computational burden when many molecular layers are evaluated.

AVAILABILITY AND IMPLEMENTATION

Source code of our framework is available at https://github.com/vgiansanti/MOWGAN.

Sox9 marks limbal stem cells and is required for asymmetric cell fate switch in the corneal epithelium

Adult tissues with high cellular turnover require a balance between stem cell renewal and differentiation, yet the mechanisms underlying this equilibrium are unclear. The cornea exhibits a polarized lateral flow of progenitors from the peripheral stem cell niche to the center; attributed to differences in cellular fate. To identify genes that are critical for regulating the asymmetric fates of limbal stem cells and their transient amplified progeny in the central cornea, we utilized an in vivo cell cycle reporter to isolate proliferating basal cells across the anterior ocular surface epithelium and performed single-cell transcriptional analysis. This strategy greatly increased the resolution and revealed distinct basal cell identities with unique expression profiles of structural genes and transcription factors. We focused on Sox9; a transcription factor implicated in stem cell regulation across various organs. Sox9 was found to be differentially expressed between limbal stem cells and their progeny in the central corneal. Lineage tracing analysis confirmed that Sox9 marks long-lived limbal stem cells and conditional deletion led to abnormal differentiation and squamous metaplasia in the central cornea. These data suggest a requirement for Sox9 for the switch to asymmetric fate and commitment toward differentiation, as transient cells exit the limbal niche. By inhibiting terminal differentiation of corneal progenitors and forcing them into perpetual symmetric divisions, we replicated the Sox9 loss-of-function phenotype. Our findings reveal an essential role for Sox9 for the spatial regulation of asymmetric fate in the corneal epithelium that is required to sustain tissue homeostasis.

Spatial Biology of Breast Cancer

Spatial findings have shaped on our understanding of breast cancer. In this review, we discuss how spatial methods, including spatial transcriptomics and proteomics and the resultant understanding of spatial relationships, have contributed to concepts regarding cancer progression and treatment. In addition to discussing traditional approaches, we examine how emerging multiplex imaging technologies have contributed to the field and how they might influence future research.

CD4+ T Cells From Individuals With Type 1 Diabetes Respond to a Novel Class of Deamidated Peptides Formed in Pancreatic Islets

The β-cell plays a crucial role in the pathogenesis of type 1 diabetes, in part through the posttranslational modification of self-proteins by biochemical processes such as deamidation. These neoantigens are potential triggers for breaking immune tolerance. We report the detection by LC-MS/MS of 16 novel Gln and 27 novel Asn deamidations in 14 disease-related proteins within inflammatory cytokine-stressed human islets of Langerhans. T-cell clones responsive against one Gln- and three Asn-deamidated peptides could be isolated from peripheral blood of individuals with type 1 diabetes. Ex vivo HLA class II tetramer staining detected higher T-cell frequencies in individuals with the disease compared with control individuals. Furthermore, there was a positive correlation between the frequencies of T cells specific for deamidated peptides, insulin antibody levels at diagnosis, and duration of disease. These results highlight that stressed human islets are prone to enzymatic and biochemical deamidation and suggest that both Gln- and Asn-deamidated peptides can promote the activation and expansion of autoreactive CD4+ T cells. These findings add to the growing evidence that posttranslational modifications undermine tolerance and may open the road for the development of new diagnostic and therapeutic applications for individuals living with type 1 diabetes.

ARTICLE HIGHLIGHTS

Emx2 underlies the development and evolution of marsupial gliding membranes

Phenotypic variation among species is a product of evolutionary changes to developmental programs1,2. However, how these changes generate novel morphological traits remains largely unclear. Here we studied the genomic and developmental basis of the mammalian gliding membrane, or patagium-an adaptative trait that has repeatedly evolved in different lineages, including in closely related marsupial species. Through comparative genomic analysis of 15 marsupial genomes, both from gliding and non-gliding species, we find that the Emx2 locus experienced lineage-specific patterns of accelerated cis-regulatory evolution in gliding species. By combining epigenomics, transcriptomics and in-pouch marsupial transgenics, we show that Emx2 is a critical upstream regulator of patagium development. Moreover, we identify different cis-regulatory elements that may be responsible for driving increased Emx2 expression levels in gliding species. Lastly, using mouse functional experiments, we find evidence that Emx2 expression patterns in gliders may have been modified from a pre-existing program found in all mammals. Together, our results suggest that patagia repeatedly originated through a process of convergent genomic evolution, whereby regulation of Emx2 was altered by distinct cis-regulatory elements in independently evolved species. Thus, different regulatory elements targeting the same key developmental gene may constitute an effective strategy by which natural selection has harnessed regulatory evolution in marsupial genomes to generate phenotypic novelty.

IL-1α facilitates GSH synthesis to counteract oxidative stress in oral squamous cell carcinoma under glucose-deprivation

Understanding the intrinsic mechanisms underpinning cancer metabolism and therapeutic resistance is of central importance for effective nutrition-starvation therapies. Here, we report that Interleukin 1A (IL1A) mRNA and IL-1α protein facilitate glutathione (GSH) synthesis to counteract oxidative stress and resistance against nutrition-starvation therapy in oral squamous cell carcinoma (OSCC). The expression of IL1A mRNA was elevated in the case of OSCC associated with unfavorable clinical outcomes. Both IL1A mRNA and IL-1α protein expression were increased under glucose-deprivation in vitro and in vivo. The transcription of IL1A mRNA was regulated in an NRF2-dependent manner in OSCC cell lines under glucose-deprivation. Moreover, the IL-1α conferred resistance to oxidative stress via GSH synthesis in OSCC cell lines. The intratumoral administration of siRNAs against IL1A mRNA markedly reversed GSH production and sensitized OSCC cells to Anlotinib in HN6 xenograft models. Overall, the current study demonstrates novel evidence that the autocrine IL-1α favors endogenous anti-oxidative process and confers therapeutic resistance to nutrition-starvation in OSCCs.

Evolution of a novel adrenal cell type that promotes parental care

Cell types with specialized functions fundamentally regulate animal behaviour, and yet the genetic mechanisms that underlie the emergence of novel cell types and their consequences for behaviour are not well understood1. Here we show that the monogamous oldfield mouse (Peromyscus polionotus) has recently evolved a novel cell type in the adrenal gland that expresses the enzyme AKR1C18, which converts progesterone into 20α-hydroxyprogesterone. We then demonstrate that 20α-hydroxyprogesterone is more abundant in oldfield mice, where it induces monogamous-typical parental behaviours, than in the closely related promiscuous deer mice (Peromyscus maniculatus). Using quantitative trait locus mapping in a cross between these species, we ultimately find interspecific genetic variation that drives expression of the nuclear protein GADD45A and the glycoprotein tenascin N, which contribute to the emergence and function of this cell type in oldfield mice. Our results provide an example by which the recent evolution of a new cell type in a gland outside the brain contributes to the evolution of social behaviour.

Emerging Roles of Spatial Transcriptomics in Liver Research

Spatial transcriptomics, leveraging sequencing- and imaging-based techniques, has emerged as a groundbreaking technology for mapping gene expression within the complex architectures of tissues. This approach provides an in-depth understanding of cellular and molecular dynamics across various states of healthy and diseased livers. Through the integration of sophisticated bioinformatics strategies, it enables detailed exploration of cellular heterogeneity, transitions in cell states, and intricate cell-cell interactions with remarkable precision. In liver research, spatial transcriptomics has been particularly revelatory, identifying distinct zonated functions of hepatocytes that are crucial for understanding the metabolic and detoxification processes of the liver. Moreover, this technology has unveiled new insights into the pathogenesis of liver diseases, such as the role of lipid-associated macrophages in steatosis and endothelial cell signals in liver regeneration and repair. In the domain of liver cancer, spatial transcriptomics has proven instrumental in delineating intratumor heterogeneity, identifying supportive microenvironmental niches and revealing the complex interplay between tumor cells and the immune system as well as susceptibility to immune checkpoint inhibitors. In conclusion, spatial transcriptomics represents a significant advance in hepatology, promising to enhance our understanding and treatment of liver diseases.

Heterogeneity of immune cells and their communications unveiled by transcriptome profiling in acute inflammatory lung injury

Background

Acute Respiratory Distress Syndrome (ARDS) or its earlier stage Acute lung injury (ALI), is a worldwide health concern that jeopardizes human well-being. Currently, the treatment strategies to mitigate the incidence and mortality of ARDS are severely restricted. This limitation can be attributed, at least in part, to the substantial variations in immunity observed in individuals with this syndrome.

Methods

Bulk and single cell RNA sequencing from ALI mice and single cell RNA sequencing from ARDS patients were analyzed. We utilized the Seurat program package in R and cellmarker 2.0 to cluster and annotate the data. The differential, enrichment, protein interaction, and cell-cell communication analysis were conducted.

Results

The mice with ALI caused by pulmonary and extrapulmonary factors demonstrated differential expression including Clec4e, Retnlg, S100a9, Coro1a, and Lars2. We have determined that inflammatory factors have a greater significance in extrapulmonary ALI, while multiple pathways collaborate in the development of pulmonary ALI. Clustering analysis revealed significant heterogeneity in the relative abundance of immune cells in different ALI models. The autocrine action of neutrophils plays a crucial role in pulmonary ALI. Additionally, there was a significant increase in signaling intensity between B cells and M1 macrophages, NKT cells and M1 macrophages in extrapulmonary ALI. The CXCL, CSF3 and MIF, TGFβ signaling pathways play a vital role in pulmonary and extrapulmonary ALI, respectively. Moreover, the analysis of human single-cell revealed DCs signaling to monocytes and neutrophils in COVID-19-associated ARDS is stronger compared to sepsis-related ARDS. In sepsis-related ARDS, CD8+ T and Th cells exhibit more prominent signaling to B-cell nucleated DCs. Meanwhile, both MIF and CXCL signaling pathways are specific to sepsis-related ARDS.

Conclusion

This study has identified specific gene signatures and signaling pathways in animal models and human samples that facilitate the interaction between immune cells, which could be targeted therapeutically in ARDS patients of various etiologies.

Flotillin-2 dampens T cell antigen-sensitivity and functionality

T cell receptor (TCR) engagement triggers T cell responses, yet how TCR-mediated activation is regulated at the plasma membrane remains unclear. Here, we report that deleting the membrane scaffolding protein Flotillin-2 (Flot2) increases T cell antigen sensitivity, resulting in enhanced TCR signaling and effector function to weak TCR stimulation. T cell-specific Flot2-deficient mice exhibited reduced tumor growth and enhanced immunity to infection. Flot2-null CD4 + T cells exhibited increased T helper 1 polarization, proliferation, Nur77 induction, and phosphorylation of ZAP70 and LCK upon weak TCR stimulation, indicating a sensitized TCR-triggering threshold. Single cell-RNA sequencing suggested that Flot2 - null CD4 + T cells follow a similar route of activation as wild-type CD4 + T cells but exhibit higher occupancy of a discrete activation state under weak TCR stimulation. Given prior reports that TCR clustering influences sensitivity of T cells to stimuli, we evaluated TCR distribution with super-resolution microscopy. Flot2 ablation increased the number of surface TCR nanoclusters on naïve CD4 + T cells. Collectively, we posit that Flot2 modulates T cell functionality to weak TCR stimulation, at least in part, by regulating surface TCR clustering. Our findings have implications for improving T cell reactivity in diseases with poor antigenicity, such as cancer and chronic infections.

scCompressSA: dual-channel self-attention based deep autoencoder model for single-cell clustering by compressing gene-gene interactions

BACKGROUND

Single-cell clustering has played an important role in exploring the molecular mechanisms about cell differentiation and human diseases. Due to highly-stochastic transcriptomics data, accurate detection of cell types is still challenged, especially for RNA-sequencing data from human beings. In this case, deep neural networks have been increasingly employed to mine cell type specific patterns and have outperformed statistic approaches in cell clustering.

RESULTS

Using cross-correlation to capture gene-gene interactions, this study proposes the scCompressSA method to integrate topological patterns from scRNA-seq data, with support of self-attention (SA) based coefficient compression (CC) block. This SA-based CC block is able to extract and employ static gene-gene interactions from scRNA-seq data. This proposed scCompressSA method has enhanced clustering accuracy in multiple benchmark scRNA-seq datasets by integrating topological and temporal features.

CONCLUSION

Static gene-gene interactions have been extracted as temporal features to boost clustering performance in single-cell clustering For the scCompressSA method, dual-channel SA based CC block is able to integrate topological features and has exhibited extraordinary detection accuracy compared with previous clustering approaches that only employ temporal patterns.

PROSTATE CELL HETEROGENEITY AND CXCL17 UPREGULATION IN MOUSE STEROID HORMONE IMBALANCE

Benign prostatic hyperplasia (BPH) is a prevalent age-related condition often characterized by debilitating urinary symptoms. Its etiology is believed to stem from hormonal imbalance, particularly an elevated estradiol-to-testosterone ratio and chronic inflammation. Our previous studies using a mouse steroid hormone imbalance model identified a specific increase in macrophages that migrate and accumulate in the prostate lumen where they differentiate into lipid-laden foam cells in mice implanted with testosterone and estradiol pellets, but not in sham animals. The current study focused on further characterizing the cellular heterogeneity of the prostate in this model as well as identifying the specific transcriptomic signature of the recruited foam cells. Moreover, we aimed to identify the epithelia-derived signals that drive macrophage infiltration and luminal translocation. Male C57BL/6J mice were implanted with slow-release testosterone and estradiol pellets (T+E2) and harvested the ventral prostates two weeks later for scRNA-seq analysis, or performed sham surgery. We identified Ear2+ and Cd72+ macrophages that were elevated in response to steroid hormone imbalance, whereas a Mrc1+ resident macrophage population did not change. In addition, an Spp1+ foam cell cluster was almost exclusively found in T+E2 mice. Further markers of foam cells were also identified, including Gpnmb and Trem2, and GPNMB was confirmed as a novel histological marker with immunohistochemistry. Foam cells were also shown to express known pathological factors Vegf, Tgfb1, Ccl6, Cxcl16 and Mmp12. Intriguingly, a screen for chemokines identified the upregulation of epithelial-derived Cxcl17, a known monocyte attractant, in T+E2 prostates suggesting that it might be responsible for the elevated macrophage number as well as their translocation to the lumen. Our study identified macrophage subsets that respond to steroid hormone imbalance as well as further confirmed a potential pathological role of luminal foam cells in the prostate. These results underscore a pathological role of the identified prostate foam cells and suggests CXCL17-mediated macrophage migration as a critical initiating event.

A novel pyroptosis-related gene signature exhibits distinct immune cells infiltration landscape in Wilms' tumor

BACKGROUND

Wilms' tumor (WT) is the most common renal tumor in childhood. Pyroptosis, a type of inflammation-characterized and immune-related programmed cell death, has been extensively studied in multiple tumors. In the current study, we aim to construct a pyroptosis-related gene signature for predicting the prognosis of Wilms' tumor.

METHODS

We acquired RNA-seq data from TARGET kidney tumor projects for constructing a gene signature, and snRNA-seq data from GEO database for validating signature-constructing genes. Pyroptosis-related genes (PRGs) were collected from three online databases. We constructed the gene signature by Lasso Cox regression and then established a nomogram. Underlying mechanisms by which gene signature is related to overall survival states of patients were explored by immune cell infiltration analysis, differential expression analysis, and functional enrichment analysis.

RESULTS

A pyroptosis-related gene signature was constructed with 14 PRGs, which has a moderate to high predicting capacity with 1-, 3-, and 5-year area under the curve (AUC) values of 0.78, 0.80, and 0.83, respectively. A prognosis-predicting nomogram was established by gender, stage, and risk score. Tumor-infiltrating immune cells were quantified by seven algorithms, and the expression of CD8( +) T cells, B cells, Th2 cells, dendritic cells, and type 2 macrophages are positively or negatively correlated with risk score. Two single nuclear RNA-seq samples of different histology were harnessed for validation. The distribution of signature genes was identified in various cell types.

CONCLUSIONS

We have established a pyroptosis-related 14-gene signature in WT. Moreover, the inherent roles of immune cells (CD8( +) T cells, B cells, Th2 cells, dendritic cells, and type 2 macrophages), functions of differentially expressed genes (tissue/organ development and intercellular communication), and status of signaling pathways (proteoglycans in cancer, signaling pathways regulating pluripotent of stem cells, and Wnt signaling pathway) have been elucidated, which might be employed as therapeutic targets in the future.

Cancer SLC6A6-mediated taurine uptake transactivates immune checkpoint genes and induces exhaustion in CD8+ T cells

Taurine is used to bolster immunity, but its effects on antitumor immunity are unclear. Here, we report that cancer-related taurine consumption causes T cell exhaustion and tumor progression. The taurine transporter SLC6A6 is correlated with aggressiveness and poor outcomes in multiple cancers. SLC6A6-mediated taurine uptake promotes the malignant behaviors of tumor cells but also increases the survival and effector function of CD8+ T cells. Tumor cells outcompete CD8+ T cells for taurine by overexpressing SLC6A6, which induces T cell death and malfunction, thereby fueling tumor progression. Mechanistically, taurine deficiency in CD8+ T cells increases ER stress, promoting ATF4 transcription in a PERK-JAK1-STAT3 signaling-dependent manner. Increased ATF4 transactivates multiple immune checkpoint genes and induces T cell exhaustion. In gastric cancer, we identify a chemotherapy-induced SP1-SLC6A6 regulatory axis. Our findings suggest that tumoral-SLC6A6-mediated taurine deficiency promotes immune evasion and that taurine supplementation reinvigorates exhausted CD8+ T cells and increases the efficacy of cancer therapies.