Single-cell mRNA quantification and differential analysis with Census

Employing single-cell RNA sequencing coupled with an array of bioinformatics approaches to ascertain the shared genetic characteristics between osteoporosis and obesity

Aims

This study examined the relationship between obesity (OB) and osteoporosis (OP), aiming to identify shared genetic markers and molecular mechanisms to facilitate the development of therapies that target both conditions simultaneously.

Methods

Using weighted gene co-expression network analysis (WGCNA), we analyzed datasets from the Gene Expression Omnibus (GEO) database to identify co-expressed gene modules in OB and OP. These modules underwent Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment and protein-protein interaction analysis to discover Hub genes. Machine learning refined the gene selection, with further validation using additional datasets. Single-cell analysis emphasized specific cell subpopulations, and enzyme-linked immunosorbent assay (ELISA), protein blotting, and cellular staining were used to investigate key genes.

Results

WGCNA revealed critical gene modules for OB and OP, identifying the Toll-like receptor (TLR) signalling pathway as a common factor. TLR2 was the most significant gene, with a pronounced expression in macrophages. Elevated TLR2 expression correlated with increased adipose accumulation, inflammation, and osteoclast differentiation, linking it to OP development.

Conclusion

Our study underscores the pivotal role of TLR2 in connecting OP and OB. It highlights the influence of TLR2 in macrophages, driving both diseases through a pro-inflammatory mechanism. These insights propose TLR2 as a potential dual therapeutic target for treating OP and OB.

Single-cell spatiotemporal analysis of the lungs reveals Slamf9+ macrophages involved in viral clearance and inflammation resolution

How the lung achieves immune homeostasis after a pulmonary infection is not fully understood. Here, we analyzed the spatiotemporal changes in the lungs over a 2-week natural recovery from severe pneumonia in a Syrian hamster model of SARS-CoV-2 infection. We find that SARS-CoV-2 infects multiple cell types and causes massive cell death at the early stage, including alveolar macrophages. We identify a group of monocyte-derived Slamf9+ macrophages, which are induced after SARS-CoV-2 infection and resistant to impairment caused by SARS-CoV-2. Slamf9+ macrophages contain SARS-CoV-2, recruit and interact with Isg12+Cst7+ neutrophils to clear the viruses. After viral clearance, Slamf9+ macrophages differentiate into Trem2+ and Fbp1+ macrophages, contributing to inflammation resolution at the late stage, and finally replenish alveolar macrophages. These findings are validated in a SARS-CoV-2-infected hACE2 mouse model and confirmed with publicly available human autopsy single-cell RNA-seq data, demonstrating the potential role of Slamf9+ macrophages and their coordination with neutrophils in post-injury tissue repair and inflammation resolution.

A prospective study of neoadjuvant pembrolizumab plus chemotherapy for resectable esophageal squamous cell carcinoma: The Keystone-001 trial

In this phase II study, 47 patients with locally advanced, resectable esophageal squamous cell carcinoma (ESCC) received three cycles of pembrolizumab plus chemotherapy, followed by Da Vinci robot-assisted surgery. The primary endpoints were safety and major pathological response (MPR). Key secondary endpoints included complete pathological response (pCR) and survival. No grade ≥3 adverse events or surgical delays occurred during neoadjuvant therapy. Among 46 patients studied for efficacy, the MPR and pCR rates were 72% and 41%, respectively. After a median follow-up of 27.2 months, the 2-year overall survival (OS) and disease-free survival (DFS) rates were 91% and 89%, respectively. Expansion of TRGC2+ NKT cells in peripheral blood correlated with neoadjuvant treatment effectiveness, which was validated by in vitro organoid experiments and external cancer datasets, and its functional classification and mechanism of action were further explored. These findings show preoperative pembrolizumab plus chemotherapy is a promising therapeutic strategy for resectable ESCC.

Resting natural killer cells promote the progress of colon cancer liver metastasis by elevating tumor-derived stem cell factor

The abundance and biological contribution of natural killer (NK) cells in cancer are controversial. Here, we aim to uncover clinical relevance and cellular roles of NK cells in colon cancer liver metastasis (CCLM). Here, we integrated single-cell RNA-sequencing, spatial transcriptomics (ST), and bulk RNA-sequencing datasets to investigate NK cells' biological properties and functions in the microenvironment of primary and liver metastatic tumors. Results were validated through an in vitro co-culture experiment based on bioinformatics analysis. Useing single-cell RNA-sequencing and ST, we mapped the immune cellular landscape of colon cancer and well-matched liver metastatic cancer. We discovered that GZMK+ resting NK cells increased significantly in tumor tissues and were enriched in the tumor regions of both diseases. After combining bulk RNA and clinical data, we observed that these NK cell subsets contributed to a worse prognosis. Meanwhile, KIR2DL4+ activated NK cells exhibited the opposite position and relevance. Pseudotime cell trajectory analysis revealed the evolution of activated to resting NK cells. In vitro experiments further confirmed that tumor-cell-co-cultured NK cells exhibited a decidual-like status, as evidenced by remarkable increasing CD9 expression. Functional experiments finally revealed that NK cells exhibited tumor-activating characteristics by promoting the dissociation of SCF (stem cell factor) on the tumor cells membrane depending on cell-to-cell interaction, as the supernatant of the co-culture system enhanced tumor progression. In summary, our findings revealed resting NK cells exhibited a clinical relevance with CCLM, which may be exploited for novel strategies to improve therapeutic outcomes for patients with CCLM.

Prognostic risk model of LIHC T-cells based on scRNA-seq and RNA-seq and the regulation of the tumor immune microenvironment

BACKGROUND

T-cell-related genes play a crucial role in LIHC development. However, a reliable prognostic profile based on risk models of these genes has yet to be identified.

METHODS

Single-cell datasets from both tumor and normal tissue samples were obtained from the GEO database. We identified T-cell marker genes and developed a genetic risk model using the TCGA-LIHC dataset, which was subsequently validated with an independent GEO dataset. We also explored the relationship between risk model predictions and immune responses.

RESULTS

We constructed a prognostic risk model using eight gene features identified through screening 860 T-cell marker genes via scRNA-seq and RNA-seq, which was subsequently integrated with the TCGA dataset. Its validity was independently confirmed using GEO and ICGC datasets. The TCGA dataset was stratified into high-risk and low-risk groups based on the risk model. Multivariate Cox regression analysis confirmed the risk score as an independent prognostic factor. GSEA indicated ribosomal transporter metabolism enrichment in the high-risk group and significant transcriptional activation in the low-risk group. ESTIMATE analysis showed higher ESTIMATE, immune, and stromal scores in the low-risk group, which also exhibited lower tumor purity than the high-risk group. Immunophenotyping revealed distinct patterns of immune cell infiltration and an immunosuppressive environment in the high-risk group.

CONCLUSIONS

This study introduces a T-cell marker-based prognostic risk model for LIHC patients. This model effectively predicted survival outcomes and immunotherapy effectiveness in LIHC patients, aligning with diverse immune responses and the distinct immunological profiles observed in the high-risk group.

Single-cell RNA sequencing reveals vascularization-associated cell subpopulations in dental pulp: PDGFRβ+ DPSCs with activated PI3K/AKT pathway

BACKGROUND

This study aims to address challenges in dental pulp regeneration therapy. The heterogeneity of DPSCs poses challenges, especially in stem cell transplantation for clinical use, particularly when sourced from donors of different ages and conditions.

METHODS

Pseudotime analysis was employed to analyze single-cell sequencing data, and immunohistochemical studies were conducted to investigate the expression of fibronectin 1 (FN1). We performed in vitro sorting of PDGFRβ+ DPSCs using flow cytometry. A series of functional assays, including cell proliferation, scratch, and tube formation assays, were performed to experimentally validate the vasculogenic capabilities of the identified PDGFRβ+ DPSC subset. Furthermore, gene-edited mouse models were utilized to demonstrate the importance of PDGFRβ+ DPSCs. Transcriptomic sequencing was conducted to compare the differences between PDGFRβ+ DPSCs and P1-DPSCs.

RESULTS

Single-cell sequencing analysis unveiled a distinct subset, PDGFRβ+ DPSCs, characterized by significantly elevated FN1 expression during dental pulp development. Subsequent cell experiments demonstrated that this subset possesses remarkable abilities to promote HUVEC proliferation, migration, and tube formation. Gene-edited mouse models confirmed the vital role of PDGFRβ+ DPSCs in dental pulp development. Transcriptomic sequencing and in vitro experiments demonstrated that the PDGFR/PI3K/AKT signaling pathway is a crucial factor mediating the proliferation rate and pro-angiogenic properties of PDGFRβ+ DPSCs.

CONCLUSION

We defined a new subset, PDGFRβ+ DPSCs, characterized by strong proliferative activity and pro-angiogenic capabilities, demonstrating significant clinical translational potential.

Pre-ciliated tubal epithelial cells are prone to initiation of high-grade serous ovarian carcinoma

The distal region of the uterine (Fallopian) tube is commonly associated with high-grade serous carcinoma (HGSC), the predominant and most aggressive form of ovarian or extra-uterine cancer. Specific cell states and lineage dynamics of the adult tubal epithelium (TE) remain insufficiently understood, hindering efforts to determine the cell of origin for HGSC. Here, we report a comprehensive census of cell types and states of the mouse uterine tube. We show that distal TE cells expressing the stem/progenitor cell marker Slc1a3 can differentiate into both secretory (Ovgp1+) and ciliated (Fam183b+) cells. Inactivation of Trp53 and Rb1, whose pathways are commonly altered in HGSC, leads to elimination of targeted Slc1a3+ cells by apoptosis, thereby preventing their malignant transformation. In contrast, pre-ciliated cells (Krt5+, Prom1+, Trp73+) remain cancer-prone and give rise to serous tubal intraepithelial carcinomas and overt HGSC. These findings identify transitional pre-ciliated cells as a cancer-prone cell state and point to pre-ciliation mechanisms as diagnostic and therapeutic targets.

Serial single-cell RNA sequencing unveils drug resistance and metastatic traits in stage IV breast cancer

Metastasis is a complex process that remains poorly understood at the molecular levels. We profiled single-cell transcriptomic, genomic, and epigenomic changes associated with cancer cell progression, chemotherapy resistance, and metastasis from a Stage IV breast cancer patient. Pretreatment- and posttreatment-specimens from the primary tumor and distant metastases were collected for single-cell RNA sequencing and subsequent cell clustering, copy number variation (CNV) estimation, transcriptomic factor estimation, and pseudotime analyses. CNV analysis revealed that a small population of pretreatment cancer cells resisted chemotherapy and expanded. New clones including Metastatic Precursor Cells (MPCs), emerged in the posttreatment primary tumors in CNV similar to metastatic cells. MPCs exhibited expression profiles indicative of epithelial-mesenchymal transition. Comparison of MPCs with metastatic cancer cells also revealed dynamic changes in transcription factors and calcitonin pathway gene expression. These findings demonstrate the utility of single-patient clinical sample analysis for understanding tumor drug resistance, regrowth, and metastasis.

Single-cell analysis reveals hypoxia-induced immunosuppressive microenvironment in intrahepatic cholangiocarcinoma

The role of hypoxia in the tumor microenvironment of intrahepatic cholangiocarcinoma (iCCA) remains unclear. Here, we generated a comprehensive atlas of the entire tumor microenvironment and delineated the multifaceted cell-cell interactions to decipher hypoxia-induced pro-tumor immune suppression. We discovered hypoxia is significantly associated with iCCA progression via the activation of HIF1A expression. Moreover, hypoxia-dependent PPARγ-mediated fatty acid oxidation in APOE+ TAMs promoted M2 macrophage polarization by activating the HIF1A-PPARG-CD36 axis. These polarized APOE+ TAMs recruited Treg cell infiltration via the CCL3-CCR5 pair to form an immunosuppressive microenvironment. APOE+ TAMs tended to co-localize spatially with Treg cells in the malignant tissue based on spatial transcriptome data and immunofluorescence analysis results. We identified tumor-reactive CXCL13+ CD8-PreTex with specific high expression of ENTPD1 and ITGAE, which acted as precursors of CD8-Tex and had higher cytotoxicity, lower exhaustion, and more vigorous proliferation. Consequently, CXCL13+ CD8-PreTex functioned as a positive regulator of antitumor immunity by expressing the pro-inflammatory cytokines IFNG and TNF, associated with a better survival outcome. Our study reveals the mechanisms involved in hypoxia-induced immunosuppression and suggests that targeting precursor-exhausted CXCL13+CD8+ T cells might provide a pratical immunotherapeutic approach.

Single-cell transcriptome and chromatin accessibility mapping of upper lip and primary palate fusion

Cleft lip and/or primary palate (CL/P) represent a prevalent congenital malformation, the aetiology of which is highly intricate. Although it is generally accepted that the condition arises from failed fusion between the upper lip and primary palate, the precise mechanism underlying this fusion process remains enigmatic. In this study, we utilized transposase-accessible chromatin sequencing (scATAC-seq) and single-cell RNA sequencing (scRNA-seq) to interrogate lambdoidal junction tissue derived from C57BL/6J mouse embryos at critical stages of embryogenesis (10.5, 11.5 and 12.5 embryonic days). We successfully identified distinct subgroups of mesenchymal and ectodermal cells involved in the fusion process and characterized their unique transcriptional profiles. Furthermore, we conducted cell differentiation trajectory analysis, revealing a dynamic repertoire of genes that are sequentially activated or repressed during pseudotime, facilitating the transition of relevant cell types. Additionally, we employed scATAC data to identify key genes associated with the fusion process and demonstrated differential chromatin accessibility across major cell types. Finally, we constructed a dynamic intercellular communication network and predicted upstream transcriptional regulators of critical genes involved in important signalling pathways. Our findings provide a valuable resource for future studies on upper lip and primary palate development, as well as congenital defects.

Single-cell transcriptomic analysis of flowering regulation and vernalization in Chinese cabbage shoot apex

In Chinese cabbage development the interplay between shoot apex activity and vernalization is pivotal for flowering timing. The intricate relationship between various cell types in the shoot apex meristem and their roles in regulating flowering gene expression in Chinese cabbage is not yet fully understood. A thorough analysis of single-cell types in the Chinese cabbage shoot apex and their influence on flowering genes and vernalization is essential for deeper insight. Our study first established a single-cell transcriptomic atlas of Chinese cabbage after 25 days of non-vernalization. Analyzing 19 602 single cells, we differentiated them into 15 distinct cell clusters using established marker genes. We found that key genes in shoot apex development and flowering were primarily present in shoot meristematic cells (SMCs), companion cells (CCs), and mesophyll cells (MCs). MADS-box protein FLOWERING LOCUS C 2 (BrFLC2), a gene suppressing flowering, was observed in CCs, mirroring patterns found in Arabidopsis. By mapping developmental trajectories of SMCs, CCs, and MCs, we elucidated the evolutionary pathways of crucial genes in shoot apex development and flowering. The creation of a single-cell transcriptional atlas of the Chinese cabbage shoot apex under vernalization revealed distinct alterations in the expression of known flowering genes, such as VERNALIZATION INSENSITIVE 3 (VIN3), VERNALIZATION 1 (VRN1), VERNALIZATION 2 (VRN2), BrFLC, and FLOWERING LOCUS T (FT), which varied by cell type. Our study underscores the transformative impact of single-cell RNA sequencing (scRNA-seq) for unraveling the complex differentiation and vernalization processes in the Chinese cabbage shoot apex. These insights are pivotal for enhancing breeding strategies and cultivation management of this vital vegetable.

The effects of age and dysfunction on meibomian gland population dynamics

PURPOSE

While meibomian gland dysfunction (MGD) is widely recognized as a major cause of evaporative dry eye disease, little is known about normal gland differentiation and lipid synthesis or the mechanism underlying gland atrophy and abnormal lipid secretion. The purpose of this study was to use single-cell and spatial transcriptomics to probe changes in cell composition, differentiation, and gene expression associated with two murine models of MGD: age-related gland atrophy in wild-type mice and altered meibum quality in acyl-CoA wax alcohol acyltransferase 2 (Awat2) knockout (KO) mice.

METHODS

Young (6 month) and old (22 month) wild type, C57Bl/6 mice and young (3 month) and old (13 month) Awat2 KO mice were used in these studies. For single-cell analysis, the tarsal plate was dissected from the upper and lower eyelids, and single cells isolated and submitted to the UCI Genomic Core, while for the spatial analysis frozen tissue sections were shipped to Resolve Biosciences on dry ice and sections probed in duplicate using a meibomian gland specific, 100 gene Molecular Chartography panel.

RESULTS

Analysis of gene expression patterns identified the stratified expression of lipogenic genes during meibocyte differentiation, which may control the progressive synthesis of meibum lipids; an age-related decrease in meibocytes; and increased immune cell infiltration. Additionally, we detected unique immune cell populations in the Awat2 KO mouse suggesting activation of psoriasis-like, inflammatory pathways perhaps caused by ductal dilation and hyperplasia.

CONCLUSION

Together these findings support novel mechanism controlling gland function and dysfunction.

Single-cell RNA Sequencing Identifies a Novel Subtype of Microglia with High Cd74 Expression that Facilitates White Matter Inflammation During Chronic Cerebral Hypoperfusion

Vascular dementia (VaD) causes progressive cognitive decline in the elderly population, but there is short of available therapeutic measures. Microglia-mediated neuroinflammation is vigorously involved in the pathogenesis of VaD, but the traditional classification of microglial M1/M2 phenotypes remains restrictive and controversial. This study aims to investigate whether microglia transform into novel subtypes in VaD. Chronic cerebral hypoperfusion (CCH) rat model was constructed to mimic VaD. Microglia were isolated via magnetic-activated cell sorting and analyzed by single-cell RNA sequencing (scRNA-seq) and bioinformatics. The findings inferred from scRNA-seq and bioinformatics were further validated through in vivo experiments. In this study, microglia were divided into eight clusters. The proportion of MG5 cluster was significantly increased in the white matter of the CCH group compared with the Sham group and was named chronic ischemia-associated microglia (CIAM). Immunity- and inflammation-related genes, including RT1-Db1, RT1-Da, RT1-Ba, Cd74, Spp1, C3, and Cd68, were markedly upregulated in CIAM. Enrichment analysis illustrated that CIAM possessed the function of evoking neuroinflammation. Further studies unveiled that Cd74 is associated with the most abundant GO terms involved in inflammation as well as cell proliferation and differentiation. In addition, microglia-specific Cd74 knockdown mediated by adeno-associated virus decreased the abundance of CIAM in the white matter, thereby mitigating inflammatory cytokine levels, alleviating white matter lesions, and improving cognitive impairment for CCH rats. These findings indicate that Cd74 is the core molecule of CIAM to trigger neuroinflammation and induce microglial differentiation to CIAM, suggesting that Cd74 may be a potential therapeutic target for VaD.

Prostate Cancer Progression Modeling Provides Insight into Dynamic Molecular Changes Associated with Progressive Disease States

Prostate cancer is a significant health concern and the most commonly diagnosed cancer in men worldwide. Understanding the complex process of prostate tumor evolution and progression is crucial for improved diagnosis, treatments, and patient outcomes. Previous studies have focused on unraveling the dynamics of prostate cancer evolution using phylogenetic or lineage analysis approaches. However, those approaches have limitations in capturing the complete disease process or incorporating genomic and transcriptomic variations comprehensively. In this study, we applied a novel computational approach to derive a prostate cancer progression model using multidimensional data from 497 prostate tumor samples and 52 tumor-adjacent normal samples obtained from The Cancer Genome Atlas study. The model was validated using data from an independent cohort of 545 primary tumor samples. By integrating transcriptomic and genomic data, our model provides a comprehensive view of prostate tumor progression, identifies crucial signaling pathways and genetic events, and uncovers distinct transcription signatures associated with disease progression. Our findings have significant implications for cancer research and hold promise for guiding personalized treatment strategies in prostate cancer.

SIGNIFICANCE

We developed and validated a progression model of prostate cancer using >1,000 tumor and normal tissue samples. The model provided a comprehensive view of prostate tumor evolution and progression.

Single-cell RNA sequencing reveals the MIF/ACKR3 receptor-ligand interaction between neutrophils and nucleus pulposus cells in intervertebral disc degeneration

OBJECTIVES

To unravel the heterogeneity and function of microenvironmental neutrophils during intervertebral disc degeneration (IDD).

METHODS

Single-cell RNA sequencing (scRNA-seq) was utilized to dissect the cellular landscape of neutrophils in intervertebral disc (IVD) tissues and their crosstalk with nucleus pulposus cells (NPCs). The expression levels of macrophage migration inhibitory factor (MIF) and ACKR3 in IVD tissues were detected. The MIF/ACKR3 axis was identified and its effects on IDD were investigated in vitro and in vivo.

RESULTS

We sequenced here 71520 single cells from 5 control and 9 degenerated IVD samples using scRNA-seq. We identified a unique cluster of neutrophils abundant in degenerated IVD tissues that highly expressed MIF and was functionally enriched in extracellular matrix organization (ECMO). Cell-to-cell communication analyses showed that this ECMO-neutrophil subpopulation was closely interacted with an effector NPCs subtype, which displayed high expression of ACKR3. Further analyses revealed that MIF was positively correlated with ACKR3 and functioned via directly binding to ACKR3 on effector NPCs. MIF inhibition attenuated degenerative changes of NPCs and extracellular matrix, which could be partially reversed by ACKR3 overexpression. Clinically, a significant correlation of high MIF/ACKR3 expression with advanced IDD grade was observed. Furthermore, we also found a positive association between MIF+ ECMO-neutrophil counts and ACKR3+ effector NPCs density as well as higher expression of the MIF/ACKR3 signaling in areas where these two cell types were neighbors.

CONCLUSIONS

These data suggest that ECMO-neutrophil promotes IDD progression by their communication with NPCs via the MIF/ACKR3 axis, which may shed light on therapeutic strategies.

Development of a CD39 nanobody and its enhancement to chimeric antigen receptor T cells efficacy against ovarian cancer in preclinical studies

Rationale: CD39, a key ectonucleotidase that drives adenosine production, acts as a critical immunosuppressive checkpoint in cancer. Although it has shown promise as a therapeutic target, clinical trials are demonstrating the need for more potent targeting approaches. This need is driving innovation towards the development of novel antibodies and the exploration of strategic combinations with a range of immunotherapies. Methods: An anti-CD39 nanobody was screened and tested for its affinity and binding ability using biolayer interferometry, ELISA and flow cytometry. Blocking ability against soluble and membrane-bound CD39 was measured after CD39 blockade. Internalization was detected using immunofluorescence. The reversal of T-cell function by the anti-CD39 antibody was assessed by CFSE-based T-cell proliferation, CD25 expression and IFN-γ secretion. The in vivo function of tumor growth inhibition was further tested in a mouse model and we also tested the phenotype of immune cells after CD39 antibody administration from tumor tissue, draining lymph nodes and peripheral blood. We inserted the antibody sequence into the chimeric antigen receptor (CAR) construct to induce MSLN CAR-T cells to secret the CD39 antibody, and the efficacy was measured in xenograft models of ovarian cancer. Results: We screened human CD39 antibodies using a VHH library and developed a single-epitope anti-CD39 nanobody, named huCD39 mAb, with high affinity and potent binding and blocking ability. The huCD39 mAb was internalized in a time-dependent manner. The in vitro study revealed that the huCD39 mAb was highly effective in enhancing T-cell proliferation and functionality. In vivo, the huCD39 mAb showed significant anti-tumor efficacy in an immunocompetent mouse model. Flow cytometry analysis demonstrated downregulated CD39 expression in immune cells after antibody administration. We also observed increased CD39 expression in ovarian cancer tissue and in activated CAR T cells. Subsequently, we developed a type of MSLN CAR-T cells secreting huCD39 mAb which showed effective eradication or inhibition in ovarian tumor xenografts. Conclusions: A novel huCD39 mAb with strong blocking ability against human CD39 and potent inhibition of tumor growth has been developed. Furthermore, a modified huCD39 mAb-secreting CAR-T cell has been generated, exhibiting superior efficacy against ovarian cancer. This provides a promising strategy for optimizing immunotherapies in ovarian cancer and potentially other malignancies.

Leveraging cell type-specificity for gene set analysis of single cell transcriptomics

Although single cell RNA-sequencing (scRNA-seq) provides unprecedented insights into the biology of complex tissues, analyzing such data on a gene-by-gene basis is challenging due to the large number of tested hypotheses and consequent low statistical power and difficult interpretation. These issues are magnified by the increased noise, significant sparsity and multi-modal distributions characteristic of single cell data. One promising approach for addressing these challenges is gene set testing, or pathway analysis. Unfortunately, statistical and biological differences between single cell and bulk transcriptomic data make it challenging to use existing gene set collections, which were developed for bulk tissue analysis, on scRNA-seq data. In this paper, we describe a procedure for customizing gene set collections originally created for bulk tissue analysis to reflect the structure of gene activity within specific cell types. Our approach leverages information about mean gene expression in the 81 human cell types profiled via scRNA-seq by the Human Protein Atlas (HPA) Single Cell Type Atlas. This HPA information is used to compute cell type-specific gene and gene set weights that can be used to filter or weight gene set collections. As demonstrated through the analysis of immune cell scRNA-seq data using gene sets from the Molecular Signatures Database (MSigDB), accounting for cell type-specificity can significantly improve gene set testing power and interpretability. An example vignette along with gene and gene set weights for the 81 HPA SCTA cell types and the MSigDB collections are available at https://hrfrost.host.dartmouth.edu/SCGeneSetOpt/.

Functional heterogeneity of cancer-associated fibroblasts with distinct neoadjuvant immunotherapy plus chemotherapy response in esophageal squamous cell carcinoma

Novel neoadjuvant immunotherapy combined with chemotherapy (neoICT) has improved outcomes for patients with esophageal squamous-cell carcinoma (ESCC), but challenges persist in low response rates and therapy resistance. Little is known about the intra-tumoral heterogeneity in the ESCC tumor microenvironment (TME) that underlies differential responses to neoadjuvant therapy. We applied single-cell RNA sequencing (scRNA-seq) profiling and multiplexed immunofluorescence staining to thoroughly decipher the TME in ESCC specimens from a neoadjuvant anti-PD1 combination therapy clinical trial. The cancer-associated fibroblasts (CAFs) population showed the significant alteration in abundance following neoadjuvant therapy. Specifically, IL6 + CCL2 + immunomodulatory CAFs and a novel CD248 + mechanoresponsive CAFs subset exhibited increasing infiltration. Mechanistically, CD248 + mechanoresponsive CAFs approached and lined the tumor nest to physically block the infiltration of CD8 + T cells and drug delivery, while IL6 + CCL2 + immunomodulatory CAFs induced therapeutic resistance with distinct IL-6 expression. Among patients treated with neoICT, we observed prominent CAF-T cell interactions. In particular, the NECTIN2-TIGIT ligand-receptor pair was enriched in treated samples, and TIGIT was identified as the major inhibitory checkpoint of T cells. Our findings demonstrate distinct alterations in TME constituent responses to neoadjuvant immunotherapy and identify functional phenotypes of CAFs associated with unfavorable therapeutic responses in patients. This provides potential targets to enhance responses to neoadjuvant therapy in ESCC.

Deciphering the role of sphingolipid metabolism in the immune microenvironment and prognosis of esophageal cancer via single-cell sequencing and bulk data analysis

BACKGROUND

Esophageal squamous cell carcinoma (ESCC) stands as a significant global health challenge, distinguished by its aggressive progression from the esophageal epithelium. Central to this malignancy is sphingolipid metabolism, a critical pathway that governs key cellular processes, including apoptosis and immune regulation, thereby influencing tumor behavior. The advent of single-cell and transcriptome sequencing technologies has catalyzed significant advancements in oncology research, offering unprecedented insights into the molecular underpinnings of cancer.

METHODS

We explored sphingolipid metabolism-related genes in ESCC using scRNA-seq data from GEO and transcriptome data from TCGA. We assessed 97 genes in epithelial cells with AUCell, UCell, and singscore algorithms, followed by bulk RNA-seq and differential analysis to identify prognosis-related genes. Immune infiltration and potential immunotherapeutic strategies were also investigated, and tumor gene mutations and drug treatment strategies were analyzed.

RESULT

Our study identified distinct gene expression patterns, highlighting ARSD, CTSA, DEGS1, and PPTQ's roles in later cellular stages. We identified seven independent prognostic genes and created a precise nomogram for prognosis.

CONCLUSION

This study integrates single-cell and transcriptomic data to provide a reliable prognostic model associated with sphingolipid metabolism and to inform immunotherapy and pharmacotherapy for ESCC at the genetic level. The findings have significant implications for precision therapy in esophageal cancer.

Single Cell Expression Analysis of Ductal Carcinoma in Situ Identifies Complex Genotypic-Phenotypic Relationships Altering Epithelial Composition

To identify mechanisms underlying the growth of ductal carcinoma in situ (DCIS) and properties that lead to progression to invasive cancer, we performed single-cell RNA-sequencing (scRNA-seq) on DCIS lesions and matched synchronous normal breast tissue. Using inferred copy number variations (CNV), we identified neoplastic epithelial cells from the clinical specimens which contained a mixture of DCIS and normal ducts. Phylogenetic analysis based on the CNVs demonstrated intratumoral clonal heterogeneity was associated with significant gene expression differences. We also classified epithelial cells into mammary cell states and found that individual genetic clones contained a mixture of cell states suggesting an ongoing pattern of differentiation after neoplastic transformation. Cell state proportions were significantly different based on estrogen receptor (ER) expression with ER-DCIS more closely resembling the distribution in the normal breast, particularly with respect to cells with basal characteristics. Using deconvolution from bulk RNA-seq in archival DCIS specimens, we show that specific alterations in cell state proportions are associated with progression to invasive cancer. Loss of an intact basement membrane (BM) is the functional definition of invasive breast cancer (IBC) and scRNA-seq data demonstrated that ongoing transcription of key BM genes occurs in specific subsets of epithelial cell states. Examining BM in archival microinvasive breast cancers and an in vitro model of invasion, we found that passive loss of BM gene expression due to cell state proportion alterations is associated with loss of the structural integrity of the duct leading to an invasive phenotype. Our analyses provide detailed insight into DCIS biology.

SIGNIFICANCE

Single cell analysis reveals that preinvasive breast cancer is comprised of multiple genetic clones and there is substantial phenotypic diversity both within and between these clones. Ductal carcinoma in situ (DCIS) of the breast is a non-invasive condition commonly identified through mammographic screening. A primary diagnosis of DCIS carries little mortality risk on its own, but its presence is a risk factor for subsequent clonally related invasive breast cancer (IBC) (1-5).

Identification of the tumor metastasis-related tumor subgroups overexpressed NENF in triple-negative breast cancer by single-cell transcriptomics

Tumor metastasis is a continuous and dynamic process and is a major cause of tumor-related death in triple-negative breast cancer. However, this biological process remains largely unknown in triple-negative breast cancer. The emergence of single-cell sequencing enables a deeper understanding of the tumor microenvironment and provides a new strategy for discovering the potential mechanism of tumor metastasis. Herein, we integrated the single-cell expression profiling of primary and metastatic triple-negative breast cancer by Seurat package. Nine tumor cell subgroups were identified. Enrichment analysis suggested tumor subgroups (C0, C4) were associated with tumor metastasis with poor prognosis in TNBC. Weighted gene co-expression network was constructed and identified NENF was a metastasis-related gene. Subsequently, RT-qPCR, Immunohistochemistry, and western blot confirmed NENF is highly expressed in TNBC tissues. And cell function assays indicated NENF promote cell invasion and migration through regulating EMT in TNBC. Finally, TIDE and Connectivity Map database suggest the candidate drugs for targeting NENF. In conclusion, our findings provide a new insight into the progression and metastasis of TNBC and uncover NENF may be a prognostic biomarker and potential therapy targets.

Analyses of single-cell and bulk RNA sequencing combined with machine learning reveal the expression patterns of disrupted mitophagy in schizophrenia

Background

Mitochondrial dysfunction is an important factor in the pathogenesis of schizophrenia. However, the relationship between mitophagy and schizophrenia remains to be elucidated.

Methods

Single-cell RNA sequencing datasets of peripheral blood and brain organoids from SCZ patients and healthy controls were retrieved. Mitophagy-related genes that were differentially expressed between the two groups were screened. The diagnostic model based on key mitophagy genes was constructed using two machine learning methods, and the relationship between mitophagy and immune cells was analyzed. Single-cell RNA sequencing data of brain organoids was used to calculate the mitophagy score (Mitoscore).

Results

We found 7 key mitophagy genes to construct a diagnostic model. The mitophagy genes were related to the infiltration of neutrophils, activated dendritic cells, resting NK cells, regulatory T cells, resting memory T cells, and CD8 T cells. In addition, we identified 12 cell clusters based on the Mitoscore, and the most abundant neurons were further divided into three subgroups. Results at the single-cell level showed that Mitohigh_Neuron established a novel interaction with endothelial cells via SPP1 signaling pathway, suggesting their distinct roles in SCZ pathogenesis.

Conclusion

We identified a mitophagy signature for schizophrenia that provides new insights into disease pathogenesis and new possibilities for its diagnosis and treatment.

Quantifiable blood TCR repertoire components associate with immune aging

T cell senescence alters the homeostasis of distinct T cell populations and results in decayed adaptive immune protection in older individuals, but a link between aging and dynamic T cell clone changes has not been made. Here, using a newly developed computational framework, Repertoire Functional Units (RFU), we investigate over 6500 publicly available TCR repertoire sequencing samples from multiple human cohorts and identify age-associated RFUs consistently across different cohorts. Quantification of RFU reduction with aging reveals accelerated loss under immunosuppressive conditions. Systematic analysis of age-associated RFUs in clinical samples manifests a potential link between these RFUs and improved clinical outcomes, such as lower ICU admission and reduced risk of complications, during acute viral infections. Finally, patients receiving bone marrow transplantation show a secondary expansion of the age-associated clones upon stem cell transfer from younger donors. Together, our results suggest the existence of a 'TCR clock' that could reflect the immune functions in aging populations.

Impact of Hashimoto's thyroiditis on the tumor microenvironment in papillary thyroid cancer: insights from single-cell analysis

This study investigates the impact of Hashimoto's thyroiditis (HT), an autoimmune disorder, on the papillary thyroid cancer (PTC) microenvironment using a dataset of 140,456 cells from 11 patients. By comparing PTC cases with and without HT, we identify HT-specific cell populations (HASCs) and their role in creating a TSH-suppressive environment via mTE3, nTE0, and nTE2 thyroid cells. These cells facilitate intricate immune-stromal communication through the MIF-(CD74+CXCR4) axis, emphasizing immune regulation in the TSH context. In the realm of personalized medicine, our HASC-focused analysis within the TCGA-THCA dataset validates the utility of HASC profiling for guiding tailored therapies. Moreover, we introduce a novel, objective method to determine K-means clustering coefficients in copy number variation inference from bulk RNA-seq data, mitigating the arbitrariness in conventional coefficient selection. Collectively, our research presents a detailed single-cell atlas illustrating HT-PTC interactions, deepening our understanding of HT's modulatory effects on PTC microenvironments. It contributes to our understanding of autoimmunity-carcinogenesis dynamics and charts a course for discovering new therapeutic targets in PTC, advancing cancer genomics and immunotherapy research.

Podocyte-specific KLF6 primes proximal tubule CaMK1D signaling to attenuate diabetic kidney disease

Diabetic kidney disease (DKD) is the main cause of chronic kidney disease worldwide. While injury to the podocytes, visceral epithelial cells that comprise the glomerular filtration barrier, drives albuminuria, proximal tubule (PT) dysfunction is the critical mediator of DKD progression. Here, we report that the podocyte-specific induction of human KLF6, a zinc-finger binding transcription factor, attenuates podocyte loss, PT dysfunction, and eventual interstitial fibrosis in a male murine model of DKD. Utilizing combination of snRNA-seq, snATAC-seq, and tandem mass spectrometry, we demonstrate that podocyte-specific KLF6 triggers the release of secretory ApoJ to activate calcium/calmodulin dependent protein kinase 1D (CaMK1D) signaling in neighboring PT cells. CaMK1D is enriched in the first segment of the PT, proximal to the podocytes, and is critical to attenuating mitochondrial fission and restoring mitochondrial function under diabetic conditions. Targeting podocyte-PT signaling by enhancing ApoJ-CaMK1D might be a key therapeutic strategy in attenuating the progression of DKD.

The use of ectopic volar fibroblasts to modify skin identity

Skin identity is controlled by intrinsic features of the epidermis and dermis and their interactions. Modifying skin identity has clinical potential, such as the conversion of residual limb and stump (nonvolar) skin of amputees to pressure-responsive palmoplantar (volar) skin to enhance prosthesis use and minimize skin breakdown. Greater keratin 9 (KRT9) expression, higher epidermal thickness, keratinocyte cytoplasmic size, collagen length, and elastin are markers of volar skin and likely contribute to volar skin resiliency. Given fibroblasts' capacity to modify keratinocyte differentiation, we hypothesized that volar fibroblasts influence these features. Bioprinted skin constructs confirmed the capacity of volar fibroblasts to induce volar keratinocyte features. A clinical trial of healthy volunteers demonstrated that injecting volar fibroblasts into nonvolar skin increased volar features that lasted up to 5 months, highlighting a potential cellular therapy.

Scoping Review: Methods and Applications of Spatial Transcriptomics in Tumor Research

Spatial transcriptomics (ST) examines gene expression within its spatial context on tissue, linking morphology and function. Advances in ST resolution and throughput have led to an increase in scientific interest, notably in cancer research. This scoping study reviews the challenges and practical applications of ST, summarizing current methods, trends, and data analysis techniques for ST in neoplasm research. We analyzed 41 articles published by the end of 2023 alongside public data repositories. The findings indicate cancer biology is an important focus of ST research, with a rising number of studies each year. Visium (10x Genomics, Pleasanton, CA, USA) is the leading ST platform, and SCTransform from Seurat R library is the preferred method for data normalization and integration. Many studies incorporate additional data types like single-cell sequencing and immunohistochemistry. Common ST applications include discovering the composition and function of tumor tissues in the context of their heterogeneity, characterizing the tumor microenvironment, or identifying interactions between cells, including spatial patterns of expression and co-occurrence. However, nearly half of the studies lacked comprehensive data processing protocols, hindering their reproducibility. By recommending greater transparency in sharing analysis methods and adapting single-cell analysis techniques with caution, this review aims to improve the reproducibility and reliability of future studies in cancer research.

A primate-specific endogenous retroviral envelope protein sequesters SFRP2 to regulate human cardiomyocyte development

Endogenous retroviruses (ERVs) occupy a significant part of the human genome, with some encoding proteins that influence the immune system or regulate cell-cell fusion in early extra-embryonic development. However, whether ERV-derived proteins regulate somatic development is unknown. Here, we report a somatic developmental function for the primate-specific ERVH48-1 (SUPYN/Suppressyn). ERVH48-1 encodes a fragment of a viral envelope that is expressed during early embryonic development. Loss of ERVH48-1 led to impaired mesoderm and cardiomyocyte commitment and diverted cells to an ectoderm-like fate. Mechanistically, ERVH48-1 is localized to sub-cellular membrane compartments through a functional N-terminal signal peptide and binds to the WNT antagonist SFRP2 to promote its polyubiquitination and degradation, thus limiting SFRP2 secretion and blocking repression of WNT/β-catenin signaling. Knockdown of SFRP2 or expression of a chimeric SFRP2 with the ERVH48-1 signal peptide rescued cardiomyocyte differentiation. This study demonstrates how ERVH48-1 modulates WNT/β-catenin signaling and cell type commitment in somatic development.

Loss of Annexin A1 in macrophages restrains efferocytosis and remodels immune microenvironment in pancreatic cancer by activating the cGAS/STING pathway

OBJECTIVE

Pancreatic cancer is an incurable malignant disease with extremely poor prognosis and a complex tumor microenvironment. We sought to characterize the role of Annexin A1 (ANXA1) in pancreatic cancer, including its ability to promote efferocytosis and antitumor immune responses.

METHODS

The tumor expression of ANXA1 and cleaved Caspase-3 (c-Casp3) and numbers of tumor-infiltrating CD68+ macrophages in 151 cases of pancreatic cancer were examined by immunohistochemistry and immunofluorescence. The role of ANXA1 in pancreatic cancer was investigated using myeloid-specific ANXA1-knockout mice. The changes in tumor-infiltrating immune cell populations induced by ANXA1 deficiency in macrophages were assessed by single-cell RNA sequencing and flow cytometry.

RESULTS

ANXA1 expression in pancreatic cancer patient samples correlated with the number of CD68+ macrophages. The percentage of ANXA1+ tumor-infiltrating macrophages negatively correlated with c-Casp3 expression and was significantly associated with worse survival. In mice, myeloid-specific ANXA1 deficiency inhibited tumor growth and was accompanied by the accumulation of apoptotic cells in pancreatic tumor tissue caused by inhibition of macrophage efferocytosis, which was dependent on cGAS-STING pathway-induced type I interferon signaling. ANXA1 deficiency significantly remodeled the intratumoral lymphocyte and macrophage compartments in tumor-bearing mice by increasing the number of effector T cells and pro-inflammatory macrophages. Furthermore, combination therapy of ANXA1 knockdown with gemcitabine and anti-programmed cell death protein-1 antibody resulted in synergistic inhibition of pancreatic tumor growth.

CONCLUSION

This research uncovers a novel role of macrophage ANXA1 in pancreatic cancer. ANXA1-mediated regulation of efferocytosis by tumor-associated macrophages promotes antitumor immune response via STING signaling, suggesting potential treatment strategies for pancreatic cancer.

Discovery of an unconventional lamprey lymphocyte lineage highlights divergent features in vertebrate adaptive immune system evolution

Lymphocyte receptors independently evolved in both jawed and jawless vertebrates with similar adaptive immune responses. However, the diversity of functional subtypes and molecular architecture in jawless vertebrate lymphocytes, comparable to jawed species, is not well defined. Here, we profile the gills, intestines, and blood of the lamprey, Lampetra morii, with single-cell RNA sequencing, using a full-length transcriptome as a reference. Our findings reveal higher tissue-specific heterogeneity among T-like cells in contrast to B-like cells. Notably, we identify a unique T-like cell subtype expressing a homolog of the nonlymphoid hematopoietic growth factor receptor, MPL-like (MPL-L). These MPL-L+ T-like cells exhibit features distinct from T cells of jawed vertebrates, particularly in their elevated expression of hematopoietic genes. We further discovered that MPL-L+ VLRA+ T-like cells are widely present in the typhlosole, gill, liver, kidney, and skin of lamprey and they proliferate in response to both a T cell mitogen and recombinant human thrombopoietin. These findings provide new insights into the adaptive immune response in jawless vertebrates, shedding new light on the evolution of adaptive immunity.

Multi-omics Analysis Revealed that the CCN Family Regulates Cell Crosstalk, Extracellular Matrix, and Immune Escape, Leading to a Poor Prognosis of Glioma

The CCN family is a group of matricellular proteins associated with the extracellular matrix. This study aims to explore the role of the CCN family in glioma development and its implications in the tumor microenvironment. Through analysis of bulk RNA-seq cohorts, correlations between CCN family expression and glioma subtypes, patient survival, and bioactive pathway enrichment were investigated. Additionally, single-cell datasets were employed to identify novel cell subgroups, followed by analyses of cell communication and transcription factors. Spatial transcriptomic analysis was utilized to validate the CCN family's involvement in glioma. Results indicate overexpression of CYR61,CTGF, and WISP1 in glioma, associated with unfavorable subtypes and reduced survival. Enrichment analyses revealed associations with oncogenic pathways, while CTGF and WISP1 expression correlated with increased infiltration of regulatory T cells and M2 macrophages. Single-cell analysis identified MES-like cells as the highest CCN expression. Moreover, intercellular signal transduction analysis demonstrated active pathways, including SPP1-CD44, in cell subgroups with elevated CYR61 and CTGF expression. Spatial transcriptomic analysis confirmed co-localization of CYR61,CTGF and SPP1-CD44 with high oncogenic pathway activity. These findings suggest that CCN family members may serve as potential prognostic biomarkers and therapeutic targets for glioma.

Tumor immune microenvironment of NSCLC with EGFR exon 20 insertions may predict efficacy of first-line ICI-combined regimen

OBJECTIVES

Non-small cell lung cancer (NSCLC) patients with exon 20 insertion mutations (ex20ins) of the epidermal growth factor receptor (EGFR) were resistant to monotherapy of immune checkpoint inhibitor (ICI). However, recent reports have shown that the combination of ICI and chemotherapy (ICI-combined regimen) exhibited certain efficacy for NSCLC with EGFR ex20ins. The mechanisms behind this phenomenon have not been thoroughly clarified. Hence, we conducted this study tofind correlations between the tumor immune microenvironment of EGFR ex20ins and the efficacy of ICI-combined regimen.

METHODS

We performed single-cell transcriptome sequencing and multiplex immunofluorescence staining (mIF) to investigate the immune microenvironment of NSCLC patients with EGFR ex20ins, L858R, and EGFR wild-type. We analyzed 15 treatment-naïve NSCLC samples utilizing single-cell RNA sequencing (scRNA-seq). Another 30 cases of EGFR L858R and 4 cases of wild-type were recruited to compare the immune microenvironment with that of EGFR ex20ins (28 cases) by mIF.

RESULTS

We observed that cell components, function and interactions varied between EGFR ex20ins, L858R, and wild-type NSCLC.We discovered similar T cell and CD8+ T cell distributions among groups but found noninferior or even better T cell activation in ex20ins patients. Infiltrating CD8+ FOXP3- T cells were significantly lower in the tumor region of EGFR ex20ins compared to wild-type. T cells from the ex20ins group had a greater tendency to promote cancer cell inflammation and epithelial-mesenchymal transition (EMT) compared to wild-type group. For macrophages, there were more M2-like macrophages in ex20ins patients. M1-like macrophages in ex20ins group produced fewer antitumor cytokines than in other groups.

CONCLUSIONS

The immune microenvironment of EGFR ex20ins is more suppressive than that of L858R and wild-type, suggesting that ICI monotherapy may not be sufficient for these patients. ICI-combined regimen might be a treatment option for EGFR ex20ins due to tumor-promoting inflammation and noninferior T cell functions in the immune microenvironment.

Advancing laser ablation assessment in hyperspectral imaging through machine learning

Hyperspectral imaging (HSI) is gaining increasing relevance in medicine, with an innovative application being the intraoperative assessment of the outcome of laser ablation treatment used for minimally invasive tumor removal. However, the high dimensionality and complexity of HSI data create a need for end-to-end image processing workflows specifically tailored to handle these data. This study addresses this challenge by proposing a multi-stage workflow for the analysis of hyperspectral data and allows investigating the performance of different components and modalities for ablation detection and segmentation. To address dimensionality reduction, we integrated principal component analysis (PCA) and t-distributed stochastic neighbor embedding (t-SNE) to capture dominant variations and reveal intricate structures, respectively. Additionally, we employed the Faster Region-based Convolutional Neural Network (Faster R-CNN) to accurately localize ablation areas. The two-stage detection process of Faster R-CNN, along with the choice of dimensionality reduction technique and data modality, significantly influenced the performance in detecting ablation areas. The evaluation of the ablation detection on an independent test set demonstrated a mean average precision of approximately 0.74, which validates the generalization ability of the models. In the segmentation component, the Mean Shift algorithm showed high quality segmentation without manual cluster definition. Our results prove that the integration of PCA, t-SNE, and Faster R-CNN enables improved interpretation of hyperspectral data, leading to the development of reliable ablation detection and segmentation systems.

Single-cell RNA sequencing reveals the heterogeneity of MYH11+ tumour-associated fibroblasts between left-sided and right-sided colorectal cancer

Colorectal cancer (CRC) exhibits considerable heterogeneity on tumour location. However, there is still a lack of comprehensive annotation regarding the characteristics and differences between the left-sided (L-CRC) and right-sided (R-CRC) CRC. Here, we performed single-cell RNA sequencing (scRNA-seq) on immune and stromal cells from 12 L-CRC and 10 R-CRC patients. We found that L-CRC exhibited stronger tumour invasion and poor prognosis compared with R-CRC. In addition, functional enrichment analysis of a normal cohort showed that fibroblasts of left colon are associated with tumour-related pathways. This suggested that the heterogeneity observed in both L-CRC and R-CRC may be influenced by the specific location within the colon itself. Further, we identified a potentially novel MYH11+ cancer-associated fibroblast (CAF) subset predominantly enriched in L-CRC. Moreover, we found that MYH11+ CAFs may promote tumour migration via interacting with macrophages, and was associated with poor prognosis in CRC. In summary, our study revealed the crucial role of MYH11+ CAFs in predicting a poor prognosis, thereby contributing valuable insights to the exploration of heterogeneity in L-CRC and R-CRC.

Deciphering cellular heterogeneity in Spodoptera frugiperda midgut cell line through single cell RNA sequencing

Using the 10x Genomics Chromium single-cell RNA sequencing (scRNA-seq) platform, we discovered unexpected heterogeneity in an established cell line developed from the midgut of the Fall armyworm, Spodoptera frugiperda, a major global pest. We analyzed the sequences of 18,794 cells and identified ten unique cellular clusters, including stem cells, enteroblasts, enterocytes and enteroendocrine cells, characterized by the expression of specific marker genes. Additionally, these studies addressed an important knowledge gap by investigating the expression of genes coding for respiratory and midgut membrane insecticide targets classified by the Insecticide Resistance Action Committee. Dual-fluorescence tagging method, fluorescence microscopy and fluorescence-activated cell sorting confirmed the expression of midgut cell type-specific genes. Stem cells were isolated from the heterogeneous population of SfMG-0617 cells. Our results, validated by KEGG and Gene Ontology analyses and supported by Monocle 3.0, advance the fields of midgut cellular biology and establish standards for scRNA-seq studies in non-model organisms.

Influenza vaccination stimulates maturation of the human T follicular helper cell response

The differentiation and specificity of human CD4+ T follicular helper cells (TFH cells) after influenza vaccination have been poorly defined. Here we profiled blood and draining lymph node (LN) samples from human volunteers for over 2 years after two influenza vaccines were administered 1 year apart to define the evolution of the CD4+ TFH cell response. The first vaccination induced an increase in the frequency of circulating TFH (cTFH) and LN TFH cells at week 1 postvaccination. This increase was transient for cTFH cells, whereas the LN TFH cells further expanded during week 2 and remained elevated in frequency for at least 3 months. We observed several distinct subsets of TFH cells in the LN, including pre-TFH cells, memory TFH cells, germinal center (GC) TFH cells and interleukin-10+ TFH cell subsets beginning at baseline and at all time points postvaccination. The shift toward a GC TFH cell phenotype occurred with faster kinetics after the second vaccine compared to the first vaccine. We identified several influenza-specific TFH cell clonal lineages, including multiple responses targeting internal influenza virus proteins, and found that each TFH cell state was attainable within a clonal lineage. Thus, human TFH cells form a durable and dynamic multitissue network.

Prostate cancer cancer-associated fibroblasts with stable markers post-androgen deprivation therapy associated with tumor progression and castration resistant prostate cancer

The specificity and clinical relevance of cancer-associated fibroblasts (CAFs) in prostate cancer (PCa), as well as the effect of androgen deprivation therapy (ADT) on CAFs, remain to be fully elucidated. Using cell lineage diversity and weighted gene co-expression network analysis (WGCNA), we pinpointed a unique CAF signature exclusive to PCa. The specificity of this CAF signature was validated through single-cell RNA sequencing (scRNA-seq), cell line RNA sequencing, and immunohistochemistry. This signature associates CAFs with tumor progression, elevated Gleason scores, and the emergence of castration resistant prostate cancer (CRPC). Using scRNA-seq on collected samples, we demonstrated that the CAF-specific signature is not altered by ADT, maintaining its peak signal output. Identifying a PCa-specific CAF signature and observing signaling changes in CAFs after ADT lay essential groundwork for further PCa studies.

Multi-omics analysis reveals the landscape of tumor microenvironments in left-sided and right-sided colon cancer

Background

Distinct clinical features and molecular characteristics of left-sided colon cancer (LCC) and right-sided colon cancer (RCC) suggest significant variations in their tumor microenvironments (TME). These differences can impact the efficacy of immunotherapy, making it essential to investigate and understand these disparities.

Methods

We conducted a multi-omics analysis, including bulk RNA sequencing (bulk RNA-seq), single-cell RNA sequencing (scRNA-seq), and whole-exome sequencing (WES), to investigate the constituents and characteristic differences of the tumor microenvironment (TME) in left-sided colon cancer (LCC) and right-sided colon cancer (RCC).

Result

Deconvolution algorithms revealed significant differences in infiltrated immune cells between left-sided colon cancer (LCC) and right-sided colon cancer (RCC), including dendritic cells, neutrophils, natural killer (NK) cells, CD4 and CD8 T cells, and M1 macrophages (P < 0.05). Notably, whole-exome sequencing (WES) data analysis showed a significantly higher mutation frequency in RCC compared to LCC (82,187/162 versus 18,726/115, P < 0.01). Single-cell analysis identified predominant tumor cell subclusters in RCC characterized by heightened proliferative potential and increased expression of major histocompatibility complex class I molecules. However, the main CD8 + T cell subpopulations in RCC exhibited a highly differentiated state, marked by T cell exhaustion and recent activation, defined as tumor-specific cytotoxic T lymphocytes (CTLs). Immunofluorescence and flow cytometry results confirmed this trend. Additionally, intercellular communication analysis demonstrated a greater quantity and intensity of interactions between tumor-specific CTLs and tumor cells in RCC.

Conclusion

RCC patients with an abundance of tumor-specific cytotoxic T lymphocytes (CTLs) and increased immunogenicity of tumor cells in the TME may be better candidates for immune checkpoint inhibitor therapy.

Single-cell RNA sequencing reveals dynamics of gene expression for 2D elongation and 3D growth in Physcomitrium patens

The transition from two-dimensional (2D) to 3D growth likely facilitated plants to colonize land, but its heterogeneity is not well understood. In this study, we utilized single-cell RNA sequencing to analyze the moss Physcomitrium patens, whose morphogenesis involves a transition from 2D to 3D growth. We profiled over 17,000 single cells covering all major vegetative tissues, including 2D filaments (chloronema and caulonema) and 3D structures (bud and gametophore). Pseudotime analyses revealed larger numbers of candidate genes that determine cell fates for 2D tip elongation or 3D bud differentiation. Using weighted gene co-expression network analysis, we identified a module that connects β-type carbonic anhydrases (βCAs) with auxin. We further validated the cellular expression patterns of βCAs and demonstrated their roles in 3D gametophore development. Overall, our study provides insights into cellular heterogeneity in a moss and identifies molecular signatures that underpin the 2D-to-3D growth transition at single-cell resolution.

Lipid-associated macrophages for osimertinib resistance and leptomeningeal metastases in NSCLC

Leptomeningeal metastases (LMs) remain a devastating complication of non-small cell lung cancer (NSCLC), particularly following osimertinib resistance. We conducted single-cell RNA sequencing on cerebrospinal fluid (CSF) from EGFR-mutant NSCLC with central nervous system metastases. We found that macrophages of LMs displayed functional and phenotypic heterogeneity and enhanced immunosuppressive properties. A population of lipid-associated macrophages, namely RNASE1_M, were linked to osimertinib resistance and LM development, which was regulated by Midkine (MDK) from malignant epithelial cells. MDK exhibited significant elevation in both CSF and plasma among patients with LMs, with higher MDK levels correlating to poorer outcomes in an independent cohort. Moreover, MDK could promote macrophage M2 polarization with lipid metabolism and phagocytic function. Furthermore, malignant epithelial cells in CSF, particularly after resistance to osimertinib, potentially achieved immune evasion through CD47-SIRPA interactions with RNASE1_M. In conclusion, we revealed a specific subtype of macrophages linked to osimertinib resistance and LM development, providing a potential target to overcome LMs.

Synaptic plasticity in human thalamocortical assembloids

Synaptic plasticities, such as long-term potentiation (LTP) and depression (LTD), tune synaptic efficacy and are essential for learning and memory. Current studies of synaptic plasticity in humans are limited by a lack of adequate human models. Here, we modeled the thalamocortical system by fusing human induced pluripotent stem cell-derived thalamic and cortical organoids. Single-nucleus RNA sequencing revealed that >80% of cells in thalamic organoids were glutamatergic neurons. When fused to form thalamocortical assembloids, thalamic and cortical organoids formed reciprocal long-range axonal projections and reciprocal synapses detectable by light and electron microscopy, respectively. Using whole-cell patch-clamp electrophysiology and two-photon imaging, we characterized glutamatergic synaptic transmission. Thalamocortical and corticothalamic synapses displayed short-term plasticity analogous to that in animal models. LTP and LTD were reliably induced at both synapses; however, their mechanisms differed from those previously described in rodents. Thus, thalamocortical assembloids provide a model system for exploring synaptic plasticity in human circuits.

HMGB2 drives tumor progression and shapes the immunosuppressive microenvironment in hepatocellular carcinoma: insights from multi-omics analysis

Background

Hepatocellular carcinoma (HCC) poses a significant health burden globally, with high mortality rates despite various treatment options. Immunotherapy, particularly immune-checkpoint inhibitors (ICIs), has shown promise, but resistance and metastasis remain major challenges. Understanding the intricacies of the tumor microenvironment (TME) is imperative for optimizing HCC management strategies and enhancing patient prognosis.

Methods

This study employed a comprehensive approach integrating multi-omics approaches, including single-cell RNA sequencing (scRNA-seq), bulk RNA sequencing (Bulk RNA-seq), and validation in clinical samples using spatial transcriptomics (ST) and multiplex immunohistochemistry (mIHC). The analysis aimed to identify key factors influencing the immunosuppressive microenvironment associated with HCC metastasis and immunotherapy resistance.

Results

HMGB2 is significantly upregulated in HCCTrans, a transitional subgroup associated with aggressive metastasis. Furthermore, HMGB2 expression positively correlates with an immunosuppressive microenvironment, particularly evident in exhausted T cells. Notably, HMGB2 expression correlated positively with immunosuppressive markers and poor prognosis in HCC patients across multiple cohorts. ST combined with mIHC validated the spatial expression patterns of HMGB2 within the TME, providing additional evidence of its role in HCC progression and immune evasion.

Conclusion

HMGB2 emerges as a critical player of HCC progression, metastasis, and immunosuppression. Its elevated expression correlates with aggressive tumor behavior and poor patient outcomes, suggesting its potential as both a therapeutic target and a prognostic indicator in HCC management.

TP63 transcriptionally regulates SLC7A5 to suppress ferroptosis in head and neck squamous cell carcinoma

Background

Most head and neck squamous cell carcinoma (HNSCC) patients are diagnosed at an advanced local stage. While immunotherapy has improved survival rates, only a minority of patients respond durably to targeted immunotherapies, posing substantial clinical challenges. We investigated the heterogeneity of the tumor microenvironment in HNSCC cohorts before and after immunotherapy by analyzing single-cell RNA sequencing (scRNA-seq) data and bulk RNA sequencing datasets retrieved from public databases.

Methods

We constructed a single-cell transcriptome landscape of HNSCC patients before and after immunotherapy and analyzed the cellular composition, developmental trajectories, gene regulatory networks, and communication patterns of different cell type subpopulations. Additionally, we assessed the expression levels of relevant indicators in HNSCC cells via western blot, ELISA, and fluorescent probe techniques.

Results

At the single-cell level, we identified a subpopulation of TP63+ SLC7A5+ HNSCC that exhibited a ferroptosis-resistant phenotype. This subpopulation suppresses ferroptosis in malignant cells through the transcriptional upregulation of SLC7A5 mediated by high TP63 expression, thereby promoting tumor growth and resistance to immunotherapy. The experimental results demonstrated that the overexpression of TP63 upregulated the expression of SLC7A5 and suppressed the concentrations of Fe2+ and ROS in HNSCC cells. By integrating bulk transcriptome data, we developed a clinical scoring model based on TP63 and SLC7A5, which are closely associated with tumor stage, revealing the significant prognostic efficacy of the TP63+ SLC7A5+ HNSCC-mediated ferroptosis mechanism in HNSCC patients.

Conclusion

Our research elucidates the TME in HNSCC before and after immunotherapy, revealing a novel mechanism by which TP63+ SLC7A5+ HNSCC inhibits ferroptosis and enhances tumor resistance via TP63-induced SLC7A5 upregulation. These insights lay the foundation for the development of more effective treatments for HNSCC.

Single-cell transcriptomics of blood identified IFIT1+ neutrophil subcluster expansion in NTM-PD patients

OBJECTIVE

Non-tuberculous mycobacterial pulmonary disease (NTM-PD) is caused by an imbalance between pathogens and impaired host immune responses. Mycobacterium avium complex (MAC) and Mycobacterium abscessus (MAB) are the two major pathogens that cause NTM-PD. In this study, we sought to dissect the transcriptomes of peripheral blood immune cells at the single-cell resolution in NTM-PD patients and explore potential clinical markers for NTM-PD diagnosis and treatment.

METHODS

Peripheral blood samples were collected from six NTM-PD patients, including three MAB-PD patients, three MAC-PD patients, and two healthy controls. We employed single-cell RNA sequencing (scRNA-seq) to define the transcriptomic landscape at a single-cell resolution. A comprehensive scRNA-seq analysis was performed, and flow cytometry was conducted to validate the results of scRNA-seq.

RESULTS

A total of 27,898 cells were analyzed. Nine T-cells, six mononuclear phagocytes (MPs), and four neutrophil subclusters were defined. During NTM infection, naïve T-cells were reduced, and effector T-cells increased. High cytotoxic activities were shown in T-cells of NTM-PD patients. The proportion of inflammatory and activated MPs subclusters was enriched in NTM-PD patients. Among neutrophil subclusters, an IFIT1+ neutrophil subcluster was expanded in NTM-PD compared to healthy controls. This suggests that IFIT1+ neutrophil subcluster might play an important role in host defense against NTM. Functional enrichment analysis of this subcluster suggested that it is related to interferon response. Cell-cell interaction analysis revealed enhanced CXCL8-CXCR1/2 interactions between the IFIT1+ neutrophil subcluster and NK cells, NKT cells, classical mononuclear phagocytes subcluster 1 (classical Mo1), classical mononuclear phagocytes subcluster 2 (classical Mo2) in NTM-PD patients compared to healthy controls.

CONCLUSIONS

Our data revealed disease-specific immune cell subclusters and provided potential new targets of NTM-PD. Specific expansion of IFIT1+ neutrophil subclusters and the CXCL8-CXCR1/2 axis may be involved in the pathogenesis of NTM-PD. These insights may have implications for the diagnosis and treatment of NTM-PD.

Discovery of an embryonically derived bipotent population of endothelial-macrophage progenitor cells in postnatal aorta

Converging evidence indicates that extra-embryonic yolk sac is the source of both macrophages and endothelial cells in adult mouse tissues. Prevailing views are that these embryonically derived cells are maintained after birth by proliferative self-renewal in their differentiated states. Here we identify clonogenic endothelial-macrophage (EndoMac) progenitor cells in the adventitia of embryonic and postnatal mouse aorta, that are independent of Flt3-mediated bone marrow hematopoiesis and derive from an early embryonic CX3CR1+ and CSF1R+ source. These bipotent progenitors are proliferative and vasculogenic, contributing to adventitial neovascularization and formation of perfused blood vessels after transfer into ischemic tissue. We establish a regulatory role for angiotensin II, which enhances their clonogenic and differentiation properties and rapidly stimulates their proliferative expansion in vivo. Our findings demonstrate that embryonically derived EndoMac progenitors participate in local vasculogenic responses in the aortic wall by contributing to the expansion of endothelial cells and macrophages postnatally.

Single-cell RNA Sequencing Identifies Natural Kill Cell-Related Transcription Factors Associated With Age-Related Macular Degeneration

Background

Age-related Macular Degeneration (AMD) poses a growing global health concern as the leading cause of central vision loss in elderly people.

Objection

This study focuses on unraveling the intricate involvement of Natural Killer (NK) cells in AMD, shedding light on their immune responses and cytokine regulatory roles.

Methods

Transcriptomic data from the Gene Expression Omnibus database were utilized, employing single-cell RNA-seq analysis. High-dimensional weighted gene co-expression network analysis (hdWGCNA) and single-cell regulatory network inference and clustering (SCENIC) analysis were applied to reveal the regulatory mechanisms of NK cells in early-stage AMD patients. Machine learning models, such as random forests and decision trees, were employed to screen hub genes and key transcription factors (TFs) associated with AMD.

Results

Distinct cell clusters were identified in the present study, especially the T/NK cluster, with a notable increase in NK cell abundance observed in AMD. Cell-cell communication analyses revealed altered interactions, particularly in NK cells, indicating their potential role in AMD pathogenesis. HdWGCNA highlighted the turquoise module, enriched in inflammation-related pathways, as significantly associated with AMD in NK cells. The SCENIC analysis identified key TFs in NK cell regulatory networks. The integration of hub genes and TFs identified CREM, FOXP1, IRF1, NFKB2, and USF2 as potential predictors for AMD through machine learning.

Conclusion

This comprehensive approach enhances our understanding of NK cell dynamics, signaling alterations, and potential predictive models for AMD. The identified TFs provide new avenues for molecular interventions and highlight the intricate relationship between NK cells and AMD pathogenesis. Overall, this study contributes valuable insights for advancing our understanding and management of AMD.

Single-Cell RNA Sequencing Reveals a Unique Fibroblastic Subset and Immune Disorder in Lichen Sclerosus Urethral Stricture

Purpose

Lichen sclerosus urethral stricture disease (LS USD) is a refractory and progressive disease primarily affecting the anterior urethra in males. Various potential etiological factors, such as genetics, autoimmunity, infection, and exposure to infectious urine, have been suggested. However, the accurate etiology of LS in the male urethra remains unclear.

Patients and Methods

In this study, we conducted single-cell RNA sequencing to identify the transcriptional profiles of three patients with LS USD and three patients with non-LS USD. Immunofluorescence was used to confirm the single-cell sequence results.

Results

Our study revealed distinct subsets of vein endothelial cells (ECs), smooth muscle cells (SMCs), and fibroblasts (FBs) with high proportions in LS USD, contributing to the tissue microenvironment primarily involved in proinflammatory and immune responses. In particular, FBs displayed a unique subset, Fib7, which is exclusively present in LS USD, and exhibited high expression levels of SAA1 and SAA2. The accumulation of macrophages, along with the dysregulated ratios of M1/M2-like phenotype macrophages, may be engaged in the pathogenesis of LS USD. Through cell-cell communication analysis, we identified significant interactions involving CXCL8/ACKR1 and CCR7/CCL19 in LS USD. Remarkably, Fib7 exhibited exclusive communication with IL-1B macrophages through the SAA1/FPR2 receptor-ligand pair.

Conclusion

Our study provides a profound understanding of the tissue microenvironment in LS USD, which may be valuable for understanding the pathogenesis of LS USD.

The significance of CD16+ monocytes in the occurrence and development of chronic thromboembolic pulmonary hypertension: insights from single-cell RNA sequencing

Background

Chronic thromboembolic pulmonary hypertension (CTEPH) is a serious pulmonary vascular disease characterized by residual thrombi in the pulmonary arteries and distal pulmonary microvascular remodeling. The pathogenesis of CTEPH remains unclear, but many factors such as inflammation, immunity, coagulation and angiogenesis may be involved. Monocytes are important immune cells that can differentiate into macrophages and dendritic cells and play an important role in thrombus formation. However, the distribution, gene expression profile and differentiation trajectory of monocyte subsets in CTEPH patients have not been systematically studied. This study aims to reveal the characteristics and functions of monocytes in CTEPH patients using single-cell sequencing technology, and to provide new insights for the diagnosis and treatment of CTEPH.

Methods

Single-cell RNA sequencing (scRNA-seq) were performed to analyze the transcriptomic features of peripheral blood mononuclear cells (PBMCs) from healthy controls, CTEPH patients and the tissues from CTEPH patients after the pulmonary endarterectomy (PEA). We established a CTEPH rat model with chronic pulmonary embolism caused by repeated injection of autologous thrombi through a central venous catheter, and used flow cytometry to detect the proportion changes of monocyte subsets in CTEPH patients and CTEPH rat model. We also observed the infiltration degree of macrophage subsets in thrombus tissue and their differentiation relationship with peripheral blood monocyte subsets by immunofluorescence staining.

Results

The results showed that the monocyte subsets in peripheral blood of CTEPH patients changed significantly, especially the proportion of CD16+ monocyte subset increased. This monocyte subset had unique functional features at the transcriptomic level, involving processes such as cell adhesion, T cell activation, coagulation response and platelet activation, which may play an important role in pulmonary artery thrombus formation and pulmonary artery intimal remodeling. In addition, we also found that the macrophage subsets in pulmonary endarterectomy tissue of CTEPH patients showed pro-inflammatory and lipid metabolism reprogramming features, which may be related to the persistence and insolubility of pulmonary artery thrombi and the development of pulmonary hypertension. Finally, we also observed that CD16+ monocyte subset in peripheral blood of CTEPH patients may be recruited to pulmonary artery intimal tissue and differentiate into macrophage subset with high expression of IL-1β, participating in disease progression.

Conclusion

CD16+ monocytes subset had significant gene expression changes in CTEPH patients, related to platelet activation, coagulation response and inflammatory response. And we also found that these cells could migrate to the thrombus and differentiate into macrophages with high expression of IL-1β involved in CTEPH disease progression. We believe that CD16+ monocytes are important participants in CTEPH and potential therapeutic targets.

PTBP1-mediated repression of neuron-specific CDC42 splicing constitutes a genomic alteration-independent, developmentally conserved vulnerability in IDH-wildtype glioblastoma

Gene co-expression networks may encode hitherto inadequately recognized vulnerabilities for adult gliomas. By identifying evolutionally conserved gene co-expression modules around EGFR (EM) or PDGFRA (PM), we recently proposed an EM/PM classification scheme, which assigns IDH-wildtype glioblastomas (GBM) into the EM subtype committed in neural stem cell compartment, IDH-mutant astrocytomas and oligodendrogliomas into the PM subtype committed in early oligodendrocyte lineage. Here, we report the identification of EM/PM subtype-specific gene co-expression networks and the characterization of hub gene polypyrimidine tract-binding protein 1 (PTBP1) as a genomic alteration-independent vulnerability in IDH-wildtype GBM. Supervised by the EM/PM classification scheme, we applied weighted gene co-expression network analysis to identify subtype-specific global gene co-expression modules. These gene co-expression modules were characterized for their clinical relevance, cellular origin and conserved expression pattern during brain development. Using lentiviral vector-mediated constitutive or inducible knockdown, we characterized the effects of PTBP1 on the survival of IDH-wildtype GBM cells, which was complemented with the analysis of PTBP1-depedent splicing pattern and overexpression of splicing target neuron-specific CDC42 (CDC42-N) isoform.  Transcriptomes of adult gliomas can be robustly assigned into 4 large gene co-expression modules that are prognostically relevant and are derived from either malignant cells of the EM/PM subtypes or tumor microenvironment. The EM subtype is associated with a malignant cell-intrinsic gene module involved in pre-mRNA splicing, DNA replication and damage response, and chromosome segregation, and a microenvironment-derived gene module predominantly involved in extracellular matrix organization and infiltrating immune cells. The PM subtype is associated with two malignant cell-intrinsic gene modules predominantly involved in transcriptional regulation and mRNA translation, respectively. Expression levels of these gene modules are independent prognostic factors and malignant cell-intrinsic gene modules are conserved during brain development. Focusing on the EM subtype, we identified PTBP1 as the most significant hub for the malignant cell-intrinsic gene module. PTBP1 is not altered in most glioma genomes. PTBP1 represses the conserved splicing of CDC42-N. PTBP1 knockdown or CDC42-N overexpression disrupts actin cytoskeleton dynamics, causing accumulation of reactive oxygen species and cell apoptosis. PTBP1-mediated repression of CDC42-N splicing represents a potential genomic alteration-independent, developmentally conserved vulnerability in IDH-wildtype GBM.

RankCompV3: a differential expression analysis algorithm based on relative expression orderings and applications in single-cell RNA transcriptomics

BACKGROUND

Effective identification of differentially expressed genes (DEGs) has been challenging for single-cell RNA sequencing (scRNA-seq) profiles. Many existing algorithms have high false positive rates (FPRs) and often fail to identify weak biological signals.

RESULTS

We present a novel method for identifying DEGs in scRNA-seq data called RankCompV3. It is based on the comparison of relative expression orderings (REOs) of gene pairs which are determined by comparing the expression levels of a pair of genes in a set of single-cell profiles. The numbers of genes with consistently higher or lower expression levels than the gene of interest are counted in two groups in comparison, respectively, and the result is tabulated in a 3 × 3 contingency table which is tested by McCullagh's method to determine if the gene is dysregulated. In both simulated and real scRNA-seq data, RankCompV3 tightly controlled the FPR and demonstrated high accuracy, outperforming 11 other common single-cell DEG detection algorithms. Analysis with either regular single-cell or synthetic pseudo-bulk profiles produced highly concordant DEGs with the ground-truth. In addition, RankCompV3 demonstrates higher sensitivity to weak biological signals than other methods. The algorithm was implemented using Julia and can be called in R. The source code is available at https://github.com/pathint/RankCompV3.jl .

CONCLUSIONS

The REOs-based algorithm is a valuable tool for analyzing single-cell RNA profiles and identifying DEGs with high accuracy and sensitivity.