DoubletFinder: Doublet Detection in Single-Cell RNA Sequencing Data Using Artificial Nearest Neighbors

Cerebral hypoperfusion exacerbates traumatic brain injury in male but not female mice

Mild-moderate traumatic brain injuries (TBIs) are prevalent, and while many individuals recover, there is evidence that a significant number experience long-term health impacts, including increased vulnerability to neurodegenerative diseases. These effects are influenced by other risk factors, such as cardiovascular disease. Our study tested the hypothesis that a pre-injury reduction in cerebral blood flow (CBF), mimicking cardiovascular disease, worsens TBI recovery. We induced bilateral carotid artery stenosis (BCAS) and a mild-moderate closed-head TBI in male and female mice, either alone or in combination, and analyzed CBF, spatial learning, memory, axonal damage, and gene expression. Findings showed that BCAS and TBI independently caused a ~10% decrease in CBF. Mice subjected to both BCAS and TBI experienced more significant CBF reductions, notably affecting spatial learning and memory, particularly in males. Additionally, male mice showed increased axonal damage with both BCAS and TBI compared to either condition alone. Females exhibited spatial memory deficits due to BCAS, but these were not worsened by subsequent TBI. Gene expression analysis in male mice highlighted that TBI and BCAS individually altered neuronal and glial profiles. However, the combination of BCAS and TBI resulted in markedly different transcriptional patterns. Our results suggest that mild cerebrovascular impairments, serving as a stand-in for preexisting cardiovascular conditions, can significantly worsen TBI outcomes in males. This highlights the potential for mild comorbidities to modify TBI outcomes and increase the risk of secondary diseases.

Overexpression of COX7A1 Promotes the Resistance of Gastric Cancer to Oxaliplatin and Weakens the Efficacy of Immunotherapy

Gastric cancer (GC) is one of the most common clinical malignant tumors worldwide, with high morbidity and mortality. Presently, the overall response rate to immunotherapy is low, and current methods for predicting the prognosis of GC are not optimal. Therefore, novel biomarkers with accuracy, efficiency, stability, performance ratio, and wide clinical application are needed. Based on public data sets, the chemotherapy cohort and immunotherapy cohort from Sun Yat-sen University Cancer Center, a series of bioinformatics analyses, such as differential expression analysis, survival analysis, drug sensitivity prediction, enrichment analysis, tumor immune dysfunction and exclusion analysis, single-sample gene set enrichment analysis, stemness index calculation, and immune cell infiltration analysis, were performed for screening and preliminary exploration. Immunohistochemical staining and in vitro experiments were performed for further verification. Overexpression of COX7A1 promoted the resistance of GC cells to Oxaliplatin. COX7A1 may induce immune escape by regulating the number of fibroblasts and their cellular communication with immune cells. In summary, measuring the expression levels of COX7A1 in the clinic may be useful in predicting the prognosis of GC patients, the degree of chemotherapy resistance, and the efficacy of immunotherapy.

Cranioencephalic functional lymphoid units in glioblastoma

The ecosystem of brain tumors is considered immunosuppressed, but our current knowledge may be incomplete. Here we analyzed clinical cell and tissue specimens derived from patients presenting with glioblastoma or nonmalignant intracranial disease to report that the cranial bone (CB) marrow, in juxtaposition to treatment-naive glioblastoma tumors, harbors active lymphoid populations at the time of initial diagnosis. Clinical and anatomical imaging, single-cell molecular and immune cell profiling and quantification of tumor reactivity identified CD8+ T cell clonotypes in the CB that were also found in the tumor. These were characterized by acute and durable antitumor response rooted in the entire T cell developmental spectrum. In contrast to distal bone marrow, the CB niche proximal to the tumor showed increased frequencies of tumor-reactive CD8+ effector types expressing the lymphoid egress marker S1PR1. In line with this, cranial enhancement of CXCR4 radiolabel may serve as a surrogate marker indicating focal association with improved progression-free survival. The data of this study advocate preservation and further exploitation of these cranioencephalic units for the clinical care of glioblastoma.

Macrophage-fibroblast crosstalk drives Arg1-dependent lung fibrosis via ornithine loading

Monocyte-derived macrophages recruited to injured tissues induce a maladaptive fibrotic response characterized by excessive production of collagen by local fibroblasts. Macrophages initiate this programming via paracrine factors, but it is unknown whether reciprocal responses from fibroblasts enhance profibrotic polarization of macrophages. We identify macrophage-fibroblast crosstalk necessary for injury-associated fibrosis, in which macrophages induced interleukin 6 ( IL-6 ) expression in fibroblasts via purinergic receptor P2rx4 signaling, and IL-6, in turn, induced arginase 1 ( Arg1 ) expression in macrophages. Arg1 contributed to fibrotic responses by metabolizing arginine to ornithine, which fibroblasts used as a substrate to synthesize proline, a uniquely abundant constituent of collagen. Imaging of idiopathic pulmonary fibrosis (IPF) lung samples confirmed expression of ARG1 in myeloid cells, and arginase inhibition suppressed collagen expression in cultured precision-cut IPF lung slices. Taken together, we define a circuit between macrophages and fibroblasts that facilitates cross-feeding metabolism necessary for injury-associated fibrosis.

Single-cell profiling identifies a CD8bright CD244bright Natural Killer cell subset that reflects disease activity in HLA-A29-positive birdshot chorioretinopathy

Birdshot chorioretinopathy is an inflammatory eye condition strongly associated with MHC-I allele HLA-A29. The striking association with MHC-I suggests involvement of T cells, whereas natural killer (NK) cell involvement remains largely unstudied. Here we show that HLA-A29-positive birdshot chorioretinopathy patients have a skewed NK cell pool containing expanded CD16 positive NK cells which produce more proinflammatory cytokines. These NK cells contain populations that express CD8A which is involved in MHC-I recognition on target cells, display gene signatures indicative of high cytotoxic activity (GZMB, PRF1 and ISG15), and signaling through NK cell receptor CD244 (SH2D1B). Long-term monitoring of a cohort of birdshot chorioretinopathy patients with active disease identifies a population of CD8bright CD244bright NK cells, which rapidly declines to normal levels upon clinical remission following successful treatment. Collectively, these studies implicate CD8bright CD244bright NK cells in birdshot chorioretinopathy.

Vaccine Therapy for Heart Failure Targeting the Inflammatory Cytokine Igfbp7

BACKGROUND

The heart comprises many types of cells such as cardiomyocytes, endothelial cells (ECs), fibroblasts, smooth muscle cells, pericytes, and blood cells. Every cell type responds to various stressors (eg, hemodynamic overload and ischemia) and changes its properties and interrelationships among cells. To date, heart failure research has focused mainly on cardiomyocytes; however, other types of cells and their cell-to-cell interactions might also be important in the pathogenesis of heart failure.

METHODS

Pressure overload was imposed on mice by transverse aortic constriction and the vascular structure of the heart was examined using a tissue transparency technique. Functional and molecular analyses including single-cell RNA sequencing were performed on the hearts of wild-type mice and EC-specific gene knockout mice. Metabolites in heart tissue were measured by capillary electrophoresis-time of flight-mass spectrometry system. The vaccine was prepared by conjugating the synthesized epitope peptides with keyhole limpet hemocyanin and administered to mice with aluminum hydroxide as an adjuvant. Tissue samples from heart failure patients were used for single-nucleus RNA sequencing to examine gene expression in ECs and perform pathway analysis in cardiomyocytes.

RESULTS

Pressure overload induced the development of intricately entwined blood vessels in murine hearts, leading to the accumulation of replication stress and DNA damage in cardiac ECs. Inhibition of cell proliferation by a cyclin-dependent kinase inhibitor reduced DNA damage in ECs and ameliorated transverse aortic constriction-induced cardiac dysfunction. Single-cell RNA sequencing analysis revealed upregulation of Igfbp7 (insulin-like growth factor-binding protein 7) expression in the senescent ECs and downregulation of insulin signaling and oxidative phosphorylation in cardiomyocytes of murine and human failing hearts. Overexpression of Igfbp7 in the murine heart using AAV9 (adeno-associated virus serotype 9) exacerbated cardiac dysfunction, while EC-specific deletion of Igfbp7 and the vaccine targeting Igfbp7 ameliorated cardiac dysfunction with increased oxidative phosphorylation in cardiomyocytes under pressure overload.

CONCLUSIONS

Igfbp7 produced by senescent ECs causes cardiac dysfunction and vaccine therapy targeting Igfbp7 may be useful to prevent the development of heart failure.

Single-cell RNA sequencing reveals that MYBL2 in malignant epithelial cells is involved in the development and progression of ovarian cancer

Background

Ovarian carcinoma (OC) is a prevalent gynecological malignancy associated with high recurrence rates and mortality, often diagnosed at advanced stages. Despite advances in immunotherapy, immune exhaustion remains a significant challenge in achieving optimal tumor control. However, the exploration of intratumoral heterogeneity of malignant epithelial cells and the ovarian cancer tumor microenvironment is still limited, hindering our comprehensive understanding of the disease.

Materials and methods

Utilizing single-cell RNA sequencing (scRNA-seq), we comprehensively investigated the cellular composition across six ovarian cancer patients with omental metastasis. Our focus centered on analysis of the malignant epithelial cells. Employing CytoTRACE and slingshot pseudotime analyses, we identified critical subpopulations and explored associated transcription factors (TFs) influencing ovarian cancer progression. Furthermore, by integrating clinical factors from a large cohort of bulk RNA sequencing data, we have established a novel prognostic model to investigate the impact of the tumor immune microenvironment on ovarian cancer patients. Furthermore, we have investigated the condition of immunological exhaustion.

Results

Our study identified a distinct and highly proliferative subgroup of malignant epithelial cells, known as C2 TOP2A+ TCs. This subgroup primarily consisted of patients who hadn't received neoadjuvant chemotherapy. Ovarian cancer patients with elevated TOP2A expression exhibited heightened sensitivity to neoadjuvant chemotherapy (NACT). Moreover, the transcription factor MYBL2 in this subgroup played a critical role in ovarian cancer development. Additionally, we developed an independent prognostic indicator, the TOP2A TCs Risk Score (TTRS), which revealed a correlation between the High TTRS Group and unfavorable outcomes. Furthermore, immune infiltration and drug sensitivity analyses demonstrated increased responsiveness to Paclitaxel, Cisplatin, and Gemcitabine in the Low TTRS Group.

Conclusion

This research deepens our understanding of malignant epithelial cells in ovarian cancer and enhances our knowledge of the ovarian cancer immune microenvironment and immune exhaustion. We have revealed the heightened susceptibility of the C2 TOP2A+ TCs subgroup to neoadjuvant chemotherapy and emphasized the role of MYBL2 within the C2 subgroup in promoting the occurrence and progression of ovarian cancer. These insights provide valuable guidance for the management of ovarian cancer treatment.

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

BACKGROUND

Salmonid species have followed markedly divergent evolutionary trajectories in their interactions with sea lice. While sea lice parasitism poses significant economic, environmental, and animal welfare challenges for Atlantic salmon (Salmo salar) aquaculture, coho salmon (Oncorhynchus kisutch) exhibit near-complete resistance to sea lice, achieved through a potent epithelial hyperplasia response leading to rapid louse detachment. The molecular mechanisms underlying these divergent responses to sea lice are unknown.

RESULTS

We characterized the cellular and molecular responses of Atlantic salmon and coho salmon to sea lice using single-nuclei RNA sequencing. Juvenile fish were exposed to copepodid sea lice (Lepeophtheirus salmonis), and lice-attached pelvic fin and skin samples were collected 12 h, 24 h, 36 h, 48 h, and 60 h after exposure, along with control samples. Comparative analysis of control and treatment samples revealed an immune and wound-healing response that was common to both species, but attenuated in Atlantic salmon, potentially reflecting greater sea louse immunomodulation. Our results revealed unique but complementary roles of three layers of keratinocytes in the epithelial hyperplasia response leading to rapid sea lice rejection in coho salmon. Our results suggest that basal keratinocytes direct the expansion and mobility of intermediate and, especially, superficial keratinocytes, which eventually encapsulate the parasite.

CONCLUSIONS

Our results highlight the key role of keratinocytes in coho salmon's sea lice resistance and the diverged biological response of the two salmonid host species when interacting with this parasite. This study has identified key pathways and candidate genes that could be manipulated using various biotechnological solutions to improve Atlantic salmon sea lice resistance.

Cytotoxic Lymphocyte-Monocyte Complex Reflects the Dynamics of Coronavirus Disease 2019 Systemic Immune Response

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infection causes a variety of clinical manifestations, many of which originate from altered immune responses, either locally or systemically. Immune cell cross-talk occurs mainly in lymphoid organs. However, systemic cell interaction specific to coronavirus disease 2019 has not been well characterized. Here, by employing single-cell RNA sequencing and imaging flow cytometry analysis, we unraveled, in peripheral blood, a heterogeneous group of cell complexes formed by the adherence of CD14+ monocytes to different cytotoxic lymphocytes, including SARS-CoV-2-specific CD8+ T cells, γδ T cells, and natural killer T cells. These lymphocytes attached to CD14+ monocytes that showed enhanced inflammasome activation and pyroptosis-induced cell death in progression stage; in contrast, in the convalescent phase, CD14+ monocytes with elevated antigen presentation potential were targeted by cytotoxic lymphocytes, thereby restricting the excessive immune activation. Collectively, our study reports previously unrecognized cell-cell interplay in the SARS-CoV-2-specific immune response, providing new insight into the intricacy of dynamic immune cell interaction representing antiviral defense.

Ag-driven CD8+ T cell clonal expansion is a prominent feature of MASH in humans and mice

BACKGROUND AND AIMS

Chronic liver disease due to metabolic dysfunction-associated steatohepatitis (MASH) is a rapidly increasing global epidemic. MASH progression is a consequence of the complex interplay between inflammatory insults and dysregulated hepatic immune responses. T lymphocytes have been shown to accumulate in the liver during MASH, but the cause and consequence of T cell accumulation in the liver remain unclear. Our study aimed to define the phenotype and T cell receptor diversity of T cells from human cirrhotic livers and an animal model of MASH to begin resolving their function in disease.

APPROACH AND RESULTS

In these studies, we evaluated differences in T cell phenotype in the context of liver disease. Accordingly, we isolated liver resident T cell populations from humans with cirrhosis and from mice with diet-induced MASH. Using both 5' single-cell sequencing and flow cytometry, we defined the phenotype and T cell receptor repertoire of liver resident T cells during health and disease.

CONCLUSIONS

MASH-induced human cirrhosis and diet-induced MASH in mice resulted in the accumulation of activated and clonally expanded T cells in the liver. The clonally expanded T cells in the liver expressed markers of chronic antigenic stimulation, including PD1, TIGIT, and TOX. Overall, this study establishes for the first time that T cells undergo Ag-dependent clonal expansion and functional differentiation during the progression of MASH. These studies could lead to the identification of antigenic targets that drive T cell activation, clonal expansion, and recruitment to the liver during MASH.

Single-nucleus RNA sequencing reveals glial cell type-specific responses to ischemic stroke in male rodents

Neuroglia critically shape the brain´s response to ischemic stroke. However, their phenotypic heterogeneity impedes a holistic understanding of the cellular composition of the early ischemic lesion. Here we present a single cell resolution transcriptomics dataset of the brain´s acute response to infarction. Oligodendrocyte lineage cells and astrocytes range among the most transcriptionally perturbed populations and exhibit infarction- and subtype-specific molecular signatures. Specifically, we find infarction restricted proliferating oligodendrocyte precursor cells (OPCs), mature oligodendrocytes and reactive astrocytes, exhibiting transcriptional commonalities in response to ischemic injury. OPCs and reactive astrocytes are involved in a shared immuno-glial cross talk with stroke-specific myeloid cells. Within the perilesional zone, osteopontin positive myeloid cells accumulate in close proximity to CD44+ proliferating OPCs and reactive astrocytes. In vitro, osteopontin increases the migratory capacity of OPCs. Collectively, our study highlights molecular cross talk events which might govern the cellular composition of acutely infarcted brain tissue.

Prognostic value and microenvironmental crosstalk of exosome-related signatures in human epidermal growth factor receptor 2 positive breast cancer

This study aimed to determine the prognostic value and microenvironmental crosstalk of exosome-related signatures in human epidermal growth factor receptor 2 positive breast cancer (HER2+ BC). Transcriptome sequencing and clinicopathological data were downloaded from the Cancer Genome Atlas. The 10X single cell sequencing dataset was downloaded from the National Center for Biotechnology Information Gene Expression Omnibus. Exosomes-Related Genes were extracted from the ExoCarta and Gene Set Enrichment Analysis databases. FGF9, SF3B4, and EPCAM were found and deemed the most accurate predictive signatures. Patients with HER2+ BC were subtyped into three groupings by exosome prognostic gene (EPGs). The expression of SF3B4 was positively linked with the infiltration of macrophages, neutrophils, and CD4+ T cells. The expression characteristics of EPGs were associated with the biological process of "response to xenobiotic stimuli." Interactions were relatively high between malignant epithelial cells and fibroblasts, endothelial cells, monocytes, and macrophages. Malignant epithelial cells interact more with fibroblasts and endothelial cells. The migration inhibitory factor pathway was the primary outgoing signaling pattern, while the C-C motif chemokine ligand pathway was the primary incoming signaling pattern for communication between malignant epithelial cells and macrophages. This study described the role of exosome signatures in the prognosis and microenvironment of HER2+ BC and provided a basis for future research.

Glioblastoma Neurovascular Progenitor Orchestrates Tumor Cell Type Diversity

Glioblastoma (GBM) is the deadliest form of primary brain tumor with limited treatment options. Recent studies have profiled GBM tumor heterogeneity, revealing numerous axes of variation that explain the molecular and spatial features of the tumor. Here, we seek to bridge descriptive characterization of GBM cell type heterogeneity with the functional role of individual populations within the tumor. Our lens leverages a gene program-centric meta-atlas of published transcriptomic studies to identify commonalities between diverse tumors and cell types in order to decipher the mechanisms that drive them. This approach led to the discovery of a tumor-derived stem cell population with mixed vascular and neural stem cell features, termed a neurovascular progenitor (NVP). Following in situ validation and molecular characterization of NVP cells in GBM patient samples, we characterized their function in vivo. Genetic depletion of NVP cells resulted in altered tumor cell composition, fewer cycling cells, and extended survival, underscoring their critical functional role. Clonal analysis of primary patient tumors in a human organoid tumor transplantation system demonstrated that the NVP has dual potency, generating both neuronal and vascular tumor cells. Although NVP cells comprise a small fraction of the tumor, these clonal analyses demonstrated that they strongly contribute to the total number of cycling cells in the tumor and generate a defined subset of the whole tumor. This study represents a paradigm by which cell type-specific interrogation of tumor populations can be used to study functional heterogeneity and therapeutically targetable vulnerabilities of GBM.

Single-cell omics: experimental workflow, data analyses and applications

Cells are the fundamental units of biological systems and exhibit unique development trajectories and molecular features. Our exploration of how the genomes orchestrate the formation and maintenance of each cell, and control the cellular phenotypes of various organismsis, is both captivating and intricate. Since the inception of the first single-cell RNA technology, technologies related to single-cell sequencing have experienced rapid advancements in recent years. These technologies have expanded horizontally to include single-cell genome, epigenome, proteome, and metabolome, while vertically, they have progressed to integrate multiple omics data and incorporate additional information such as spatial scRNA-seq and CRISPR screening. Single-cell omics represent a groundbreaking advancement in the biomedical field, offering profound insights into the understanding of complex diseases, including cancers. Here, we comprehensively summarize recent advances in single-cell omics technologies, with a specific focus on the methodology section. This overview aims to guide researchers in selecting appropriate methods for single-cell sequencing and related data analysis.

Intercellular interaction between FAP+ fibroblasts and CD150+ inflammatory monocytes mediates fibrostenosis in Crohn's disease

Crohn's disease (CD) is marked by recurring intestinal inflammation and tissue injury, often resulting in fibrostenosis and bowel obstruction, necessitating surgical intervention with high recurrence rates. To elucidate the mechanisms underlying fibrostenosis in CD, we analyzed the transcriptome of cells isolated from the transmural ileum of patients with CD, including a trio of lesions from each patient: non-affected, inflamed, and stenotic ileum samples, and compared them with samples from patients without CD. Our computational analysis revealed that profibrotic signals from a subset of monocyte-derived cells expressing CD150 induced a disease-specific fibroblast population, resulting in chronic inflammation and tissue fibrosis. The transcription factor TWIST1 was identified as a key modulator of fibroblast activation and extracellular matrix (ECM) deposition. Genetic and pharmacological inhibition of TWIST1 prevents fibroblast activation, reducing ECM production and collagen deposition. Our findings suggest that the myeloid-stromal axis may offer a promising therapeutic target to prevent fibrostenosis in CD.

Cell-specific regulation of the circadian clock by BMAL1 in the paraventricular nucleus: Implications for regulation of systemic biological rhythms

Circadian rhythms are internal biological rhythms driving temporal tissue-specific, metabolic programs. Loss of the circadian transcription factor BMAL1 in the paraventricular nucleus (PVN) of the hypothalamus reveals its importance in metabolic rhythms, but its functions in individual PVN cells are poorly understood. Here, loss of BMAL1 in the PVN results in arrhythmicity of processes controlling energy balance and alters peripheral diurnal gene expression. BMAL1 chromatin immunoprecipitation sequencing (ChIP-seq) and single-nucleus RNA sequencing (snRNA-seq) reveal its temporal regulation of target genes, including oxytocin (OXT), and restoring circulating OXT peaks in BMAL1-PVN knockout (KO) mice rescues absent activity rhythms. While glutamatergic neurons undergo day/night changes in expression of genes involved in cell morphogenesis, astrocytes and oligodendrocytes show gene expression changes in cytoskeletal organization and oxidative phosphorylation. Collectively, our findings show diurnal gene regulation in neuronal and non-neuronal PVN cells and that BMAL1 contributes to diurnal OXT secretion, which is important for systemic diurnal rhythms.

ONECUT2 acts as a lineage plasticity driver in adenocarcinoma as well as neuroendocrine variants of prostate cancer

Androgen receptor- (AR-) indifference is a mechanism of resistance to hormonal therapy in prostate cancer (PC). Here we demonstrate that ONECUT2 (OC2) activates resistance through multiple drivers associated with adenocarcinoma, stem-like and neuroendocrine (NE) variants. Direct OC2 gene targets include the glucocorticoid receptor (GR; NR3C1) and the NE splicing factor SRRM4, which are key drivers of lineage plasticity. Thus, OC2, despite its previously described NEPC driver function, can indirectly activate a portion of the AR cistrome through epigenetic activation of GR. Mechanisms by which OC2 regulates gene expression include promoter binding, enhancement of genome-wide chromatin accessibility, and super-enhancer reprogramming. Pharmacologic inhibition of OC2 suppresses lineage plasticity reprogramming induced by the AR signaling inhibitor enzalutamide. These results demonstrate that OC2 activation promotes a range of drug resistance mechanisms associated with treatment-emergent lineage variation in PC and support enhanced efforts to therapeutically target OC2 as a means of suppressing treatment-resistant disease.

A cross-species transcriptomic analysis reveals a novel 2-dimensional classification system explaining the invasiveness heterogeneity of pancreatic neuroendocrine tumor

Pancreatic neuroendocrine tumors (PanNETs), the second most common type of primary pancreatic tumors, display notable heterogeneity in invasiveness. Current knowledge regarding genomic alterations, including DAXX/ATRX, MEN1 mutations, and copy number variations (CNVs), provides some insights into tumor invasiveness. However, the underlying reasons for the significant variation in invasiveness between insulinoma and other types of PanNETs remain unclear. To construct a comprehensive model for the stratification of prognosis, we employed analysis of both the well-established Rip1-Tag2 (RT2) mouse model of PanNETs and human PanNETs with various functional types. Firstly, by applying single-cell and bulk RNA sequencing in PanNETs from different ages and strains of RT2 mice and human PanNETs, we introduced a 2-dimensional (2D) classification system. Based on the 2-D classification system, human PanNETs were mainly classified as benign insulinomas or non-insulinomas subclusters. Non-insulinomas subtypes mainly included gastrinomas, glucagonomas, VIPomas, and NF-PanNETs, which all exhibited potential invasiveness. In addition, we discovered an enrichment of specific CNV patterns and mutations in corresponding human PanNET subclusters. Then we denoted somatic DAXX/ATRX as the 'second hit' and confounding factors for invasiveness. Finally, by combining the 2D system, DAXX/ATRX mutation status, and tumor diameter, a group of indolent PanNETs with minimal recurrence risk was identified. In conclusion, our current work constructed a comprehensive model to elucidate the heterogeneity of invasiveness in PanNETs and improve prognostic stratification.

Inflammatory Alterations to Renal Lymphatic Endothelial Cell Gene Expression in Mouse Models of Hypertension

INTRODUCTION

Hypertension (HTN) is a major cardiovascular disease that can cause and be worsened by renal damage and inflammation. We previously reported that renal lymphatic endothelial cells (LECs) increase in response to HTN and that augmenting lymphangiogenesis in the kidneys reduces blood pressure and renal pro-inflammatory immune cells in mice with various forms of HTN. Our aim was to evaluate the specific changes that renal LECs undergo in HTN.

METHODS

We performed single-cell RNA sequencing. Using the angiotensin II-induced and salt-sensitive mouse models of HTN, we isolated renal CD31+ and podoplanin+ cells.

RESULTS

Sequencing of these cells revealed three distinct cell types with unique expression profiles, including LECs. The number and transcriptional diversity of LECs increased in samples from mice with HTN, as demonstrated by 597 differentially expressed genes (p < 0.01), 274 significantly enriched pathways (p < 0.01), and 331 regulons with specific enrichment in HTN LECs. These changes demonstrate a profound inflammatory response in renal LECs in HTN, leading to an increase in genes and pathways associated with inflammation-driven growth and immune checkpoint activity in LECs.

CONCLUSION

These results reinforce and help to further explain the benefits of renal LECs and lymphangiogenesis in HTN.

A single-cell pan-cancer analysis to show the variability of tumor-infiltrating myeloid cells in immune checkpoint blockade

Myeloid cells are vital components of the immune system and have pivotal functions in orchestrating immune responses. Understanding their functions within the tumor microenvironment and their interactions with tumor-infiltrating lymphocytes presents formidable challenges across diverse cancer types, particularly with regards to cancer immunotherapies. Here, we explore tumor-infiltrating myeloid cells (TIMs) by conducting a pan-cancer analysis using single-cell transcriptomics across eight distinct cancer types, encompassing a total of 192 tumor samples from 129 patients. By examining gene expression patterns and transcriptional activities of TIMs in different cancer types, we discern notable alterations in abundance of TIMs and kinetic behaviors prior to and following immunotherapy. We also identify specific cell-cell interaction targets in immunotherapy; unique and shared regulatory profiles critical for treatment response; and TIMs associated with survival outcomes. Overall, our study illuminates the heterogeneity of TIMs and improves our understanding of tissue-specific and cancer-specific myeloid subsets within the context of tumor immunotherapies.

Detection of neoplastic-immune hybrid cells with metastatic properties in uveal melanoma

BACKGROUND

Uveal melanoma is the most common non-cutaneous melanoma and is an intraocular malignancy affecting nearly 7,000 individuals per year worldwide. Of these, approximately 50% will progress to metastatic disease for which there are currently no effective curative therapies. Despite advances in molecular profiling and metastatic stratification of uveal melanoma tumors, little is known regarding their underlying biology of metastasis. Our group has identified a disseminated neoplastic cell population characterized by co-expression of immune and melanoma proteins, circulating hybrid cells (hybrids), in patients with uveal melanoma. Compared to circulating tumor cells, which lack expression of immune proteins, hybrids are detected at an increased prevalence in peripheral blood and can be used as a non-invasive biomarker to predict metastatic progression.

METHODS

To ascertain mechanisms underlying enhanced hybrid cell dissemination we identified hybrid cells within primary uveal melanoma tumors using single cell RNA sequencing (n = 8) and evaluated their gene expression and predicted ligand-receptor interactions in relation to other melanoma and immune cells within the primary tumor. We then verified expression of upregulated hybrid pathways within patient-matched tumor and peripheral blood hybrids (n = 4) using cyclic immunofluorescence and quantified their protein expression relative to other non-hybrid tumor and disseminated tumor cells.

RESULTS

Among the top upregulated genes and pathways in hybrid cells were those involved in enhanced cell motility and cytoskeletal rearrangement, immune evasion, and altered cellular metabolism. In patient-matched tumor and peripheral blood, we verified gene expression by examining concordant protein expression for each pathway category: TMSB10 (cell motility), CD74 (immune evasion) and GPX1 (metabolism). Both TMSB10 and GPX1 were expressed on significantly higher numbers of disseminated hybrid cells compared to circulating tumor cells, and CD74 and GPX1 were expressed on more disseminated hybrids than tumor-resident hybrids. Lastly, we identified that hybrid cells express ligand-receptor signaling pathways implicated in promoting metastasis including GAS6-AXL, CXCL12-CXCR4, LGALS9-P4HB and IGF1-IGFR1.

CONCLUSION

These findings highlight the importance of TMSB10, GPX1 and CD74 for successful hybrid cell dissemination and survival in circulation. Our results contribute to the understanding of uveal melanoma tumor progression and interactions between tumor cells and immune cells in the tumor microenvironment that may promote metastasis.

Single-cell transcriptomic insights into chemotherapy-induced remodeling of the osteosarcoma tumor microenvironment

PURPOSE

Neoadjuvant chemotherapy serves as an effective strategy for treating osteosarcoma (OS) not only by targeting cancerous cells but also by influencing the tumor's immune and stromal elements. Gaining insights into how chemotherapy reshapes the tumor's local environment is crucial for advancing OS treatment protocols.

METHODS

Using single-cell RNA sequencing, this study analyzed tumor samples from patients with advanced osteosarcoma collected both before and after chemotherapy.

RESULTS

The results revealed that chemotherapy caused the remaining OS cells to express higher levels of genes associated with stemness. Additionally, this process enhances the presence of cancer-associated fibroblasts, increasing their ability to modify the extracellular matrix (ECM). Chemotherapy also increases the number of endothelial cells, albeit with compromised differentiation capabilities. Importantly, the treatment reduced the immune cell population, including myeloid and T/NK cells, particularly impacting the subpopulations with tumor-fighting capabilities.

CONCLUSION

These findings highlight the complex reaction of the tumor environment to chemotherapy, providing valuable insights into how chemotherapy influences OS cells and the tumor microenvironment (TME). This knowledge is essential for understanding OS resistance mechanisms to treatments, potentially guiding the development of novel therapies for managing advanced OS.

Single-nucleus transcriptomics uncovers a geroprotective role of YAP in primate gingival aging

Aging has a profound impact on the gingiva and significantly increases its susceptibility to periodontitis, a worldwide prevalent inflammatory disease. However, a systematic characterization and comprehensive understanding of the regulatory mechanism underlying gingival aging is still lacking. Here, we systematically dissected the phenotypic characteristics of gingiva during aging in primates and constructed the first single-nucleus transcriptomic landscape of gingival aging, by which a panel of cell type-specific signatures were elucidated. Epithelial cells were identified as the most affected cell types by aging in the gingiva. Further analyses pinpointed the crucial role of YAP in epithelial self-renew and homeostasis, which declined during aging in epithelial cells, especially in basal cells. The decline of YAP activity during aging was confirmed in the human gingival tissues, and downregulation of YAP in human primary gingival keratinocytes recapitulated the major phenotypic defects observed in the aged primate gingiva while overexpression of YAP showed rejuvenation effects. Our work provides an in-depth understanding of gingival aging and serves as a rich resource for developing novel strategies to combat aging-associated gingival diseases, with the ultimate goal of advancing periodontal health and promoting healthy aging.

Inhibition of CD226 Co-Stimulation Suppresses Diabetes Development in the NOD Mouse by Augmenting Tregs and Diminishing Effector T Cell Function

Aims/hypothesis

Immunotherapeutics targeting T cells are crucial for inhibiting autoimmune disease progression proximal to disease onset in type 1 diabetes. A growing number of T cell-directed therapeutics have demonstrated partial therapeutic efficacy, with anti-CD3 (α-CD3) representing the only regulatory agency-approved drug capable of slowing disease progression through a mechanism involving the induction of partial T cell exhaustion. There is an outstanding need to augment the durability and effectiveness of T cell targeting by directly restraining proinflammatory T helper type 1 (Th1) and type 1 cytotoxic CD8 + T cell (Tc1) subsets, while simultaneously augmenting regulatory T cell (Treg) activity. Here, we present a novel strategy for reducing diabetes incidence in the NOD mouse model using a blocking monoclonal antibody targeting the type 1 diabetes-risk associated T cell co-stimulatory receptor, CD226.

Methods

Female NOD mice were treated with anti-CD226 between 7-8 weeks of age and then monitored for diabetes incidence and therapeutic mechanism of action.

Results

Compared to isotype-treated controls, anti-CD226 treated NOD mice showed reduced insulitis severity at 12 weeks and decreased disease incidence at 30 weeks. Flow cytometric analysis performed five weeks post-treatment demonstrated reduced proliferation of CD4 + and CD8 + effector memory T cells in spleens of anti-CD226 treated mice. Phenotyping of pancreatic Tregs revealed increased CD25 expression and STAT5 phosphorylation following anti-CD226, with splenic Tregs displaying augmented suppression of CD4 + T cell responders in vitro. Anti-CD226 treated mice exhibited reduced frequencies of islet-specific glucose-6-phosphatase catalytic subunit related protein (IGRP)-reactive CD8 + T cells in the pancreas, using both ex vivo tetramer staining and single-cell T cell receptor sequencing (scTCR-seq) approaches. 51 Cr-release assays demonstrated reduced cell-mediated lysis of beta-cells by anti-CD226-treated autoreactive cytotoxic T lymphocytes.

Conclusions/interpretation

CD226 blockade reduces T cell cytotoxicity and improves Treg function, representing a targeted and rational approach for restoring immune regulation in type 1 diabetes.

Research in Context

What is already known about this subject?: The co-stimulatory receptor CD226 is upregulated upon activation and is highly expressed on NK cell subsets, myeloid cells, and effector T cells. A single nucleotide polymorphism in CD226 ( rs763361 ; C>T) results in a Gly307Ser missense mutation linked to genetic susceptibility for type 1 diabetes. Global knockout of Cd226 and conditional Cd226 knockout in FoxP3 + Tregs reduced insulitis severity and diabetes incidence in NOD mice, indicating a crucial role for CD226 in disease pathogenesis. What is the key question?: Can CD226 blockade reduce T cell cytotoxicity and improve Treg function to diminish diabetes incidence in NOD mice?What are the new findings?: Anti-CD226 treatment reduced insulitis, decreased disease incidence, and inhibited splenic CD4 + and CD8 + effector memory T cell proliferation. Pancreatic Tregs from anti-CD226 treated mice exhibited increased CD25 expression; splenic Tregs displayed augmented STAT5 phosphorylation and suppressive capacity in vitro . Anti-CD226 treatment reduced IGRP-specific pancreatic CD8 + T cell frequencies, and reduced autoreactive CD8 + T cell-mediated lysis of beta-cells in vitro . How might this impact on clinical practice in the foreseeable future?: CD226 blockade could reduce autoreactive T cell cytotoxicity, enhance Treg function, and slow disease progression in high-risk or recent-onset type 1 diabetes cases.

Integrating multi-omics data to analyze the potential pathogenic mechanism of CTSH gene involved in type 1 diabetes in the exocrine pancreas

Type 1 diabetes (T1D) is an autoimmune disease caused by the destruction of insulin-producing pancreatic islet beta cells. Despite significant advancements, the precise pathogenesis of the disease remains unknown. This work integrated data from expression quantitative trait locus (eQTL) studies with Genome wide association study (GWAS) summary data of T1D and single-cell transcriptome data to investigate the potential pathogenic mechanisms of the CTSH gene involved in T1D in exocrine pancreas. Using the summary data-based Mendelian randomization (SMR) approach, we obtained four potential causative genes associated with T1D: BTN3A2, PGAP3, SMARCE1 and CTSH. To further investigate these genes'roles in T1D development, we validated them using a scRNA-seq dataset from pancreatic tissues of both T1D patients and healthy controls. The analysis showed a significantly high expression of the CTSH gene in T1D acinar cells, whereas the other three genes showed no significant changes in the scRNA-seq data. Moreover, single-cell WGCNA analysis revealed the strongest positive correlation between the module containing CTSH and T1D. In addition, we found cellular ligand-receptor interactions between the acinar cells and different cell types, especially ductal cells. Finally, based on functional enrichment analysis, we hypothesized that the CTSH gene in the exocrine pancreas enhances the antiviral response, leading to the overexpression of pro-inflammatory cytokines and the development of an inflammatory microenvironment. This process promotes β cells injury and ultimately the development of T1D. Our findings offer insights into the underlying pathogenic mechanisms of T1D.

CCN1 Is a Therapeutic Target for Reperfused Ischemic Brain Injury

Ischemic stroke can lead to systemic inflammation, which can activate peripheral immune cells, causing neuroinflammation and brain injury. Meningeal lymphatics play a crucial role in transporting solutes and immune cells out of the brain and draining them into cervical lymph nodes (CLNs). However, the role of meningeal lymphatics in regulating systemic inflammation during the reperfusion stage after ischemia is not well understood. In this study, we demonstrated that brain infarct size, neuronal loss, and the effector function of inflammatory macrophage subsets were reduced after ischemia-reperfusion and disruption of meningeal lymphatics. Spatial memory function was improved in the late stage of ischemic stroke following meningeal lymphatic disruption. Brain-infiltrating immune cells, including neutrophils, monocytes, and T and natural killer cells, were reduced after cerebral ischemia-reperfusion and meningeal lymphatic disruption. Single-cell RNA sequencing analysis revealed that meningeal lymphatic disruption reprogrammed the transcriptome profile related to chemotaxis and leukocyte migration in CLN lymphatic endothelial cells (LECs), and it also decreased chemotactic CCN1 expression in floor LECs. Replenishment of CCN1 through intraventricular injection increased brain infarct size and neuronal loss, while restoring numbers of macrophages/microglia in the brains of meningeal lymphatic-disrupted mice after ischemic stroke. Blocking CCN1 in cerebrospinal fluid reduced brain infarcts and improves spatial memory function after ischemia-reperfusion injury. In summary, this study indicates that CCN1-mediated detrimental inflammation was alleviated after cerebral ischemia-reperfusion injury and meningeal lymphatic disruption. CCN1 represents a novel therapeutic target for inhibiting systemic inflammation in the brain-CLN axis after ischemia-reperfusion injury.

Mucosal signatures of pathogenic T cells in HLA-B*27+ anterior uveitis and axial spondyloarthritis

HLA-B*27 was one of the first HLA alleles associated with an autoimmune disease, i.e., axial spondyloarthritis (axSpA) and acute anterior uveitis (B27AAU), which cause joint and eye inflammation, respectively. Gastrointestinal inflammation has been suggested as a trigger of axSpA. We recently identified a bacterial peptide (YeiH) that can be presented by HLA-B*27 to expanded public T cell receptors in the joint in axSpA and the eye in B27AAU. While YeiH is present in enteric microbiota and pathogens, additional evidence that pathogenic T cells in HLA-B*27-associated autoimmunity may have had a prior antigenic encounter within the gastrointestinal tract remains lacking. Here, we analyzed ocular, synovial, and blood T cells in B27AAU and axSpA, showing that YeiH-specific CD8+ T cells express a mucosal gene set and surface proteins consistent with intestinal differentiation, including CD161, integrin α4β7, and CCR6. In addition, we found an expansion of YeiH-specific CD8+ T cells in axSpA and B27AAU blood compared with that from individuals acting as healthy controls, whereas influenza-specific CD8+ T cells were equivalent across groups. Finally, we demonstrated the dispensability of TRBV9 for antigen recognition. Collectively, our data suggest that, in HLA-B27-associated autoimmunity, early antigen exposure and differentiation of pathogenic CD8+ T cells may occur in enteric organs.

Integration of single-cell and bulk RNA sequencing revealed immune heterogeneity and its association with disease activity in rheumatoid arthritis patients

Rheumatoid arthritis (RA) is a chronic, inflammatory, systemic autoimmune disease characterized by cartilage, bone damage, synovial inflammation, hyperplasia, autoantibody production, and systemic features. To obtain an overall profile of the immune environment in RA patients and its association with clinical features, we performed single-cell transcriptome and T-cell receptor sequencing of mononuclear cells from peripheral blood (PBMC) and synovial fluid (SF) from RA patients, integrated with two large cohorts with bulk RNA sequencing for further validation and investigation. Dendritic cells (DCs) exhibited relatively high functional heterogeneity and tissue specificity in relation to both antigen presentation and proinflammatory functions. Peripheral helper T cells (TPHs) are likely to originate from synovial tissue, undergo activation and exhaustion, and are subsequently released into the peripheral blood. Notably, among all immune cell types, TPHs were found to have the most intense associations with disease activity. In addition, CD8 effector T cells could be clustered into two groups with different cytokine expressions and play distinct roles in RA development. By integrating single-cell data with bulk sequencing from two large cohorts, we identified interactions among TPHs, CD8 cells, CD16 monocytes, and DCs that strongly contribute to the proinflammatory local environment in RA joints. Of note, the swollen 28-joint counts exhibited a more pronounced association with this immune environment compared to other disease activity indexes. The immune environment alternated significantly from PBMCs to SF, which indicated that a series of immune cells was involved in proinflammatory responses in the local joints of RA patients. By integrating single-cell data with two large cohorts, we have uncovered associations between specific immune cell populations and clinical features. This integration provides a rapid and precise methodology for assessing local immune activation, offering valuable insights into the pathophysiological mechanisms at play in RA.

Mapping Single-Cell Transcriptomes of Endometrium Reveals Potential Biomarkers in Endometrial Cancer

Background

The heterogeneity and dynamic changes of endometrial cells have a significant impact on health as they determine the normal function of the endometrium during the menstrual cycle. Dysfunction of the endometrium can lead to the occurrence of various gynecological diseases. Therefore, deconvolution of immune microenvironment that drives transcriptional programs throughout the menstrual cycle is key to understand regulatory biology of endometrium.

Methods

Herein, we comprehensively analyzed single-cell transcriptome of 59,397 cells across ten human endometrium samples and revealed the dynamic cellular heterogeneity throughout the menstrual cycle.

Results

We identified two perivascular cell subtypes, four epithelial subtypes and four fibroblast cell types in endometrium. Moreover, we inferred the cell type-specific transcription factor (TF) activities and linked critical TFs to transcriptional output of diverse immune cell types, highlighting the importance of transcriptional regulation in endometrium. Dynamic interactions between various types of cells in endometrium contribute to a range of biological pathways regulating differentiation of secretory. Integration of the molecular biomarkers identified in endometrium and bulk transcriptome of 535 endometrial cancers (EC), we revealed five RNA-based molecular subtypes of EC with highly intratumoral heterogeneity and different clinical manifestations. Mechanism analysis uncovered clinically relevant pathways for pathogenesis of EC.

Conclusion

In summary, our results revealed the dynamic immune microenvironment of endometrium and provided novel insights into future development of RNA-based treatments for endometriosis and endometrial carcinoma.

Genomic and single-cell analyses reveal genetic signatures of swimming pattern and diapause strategy in jellyfish

Jellyfish exhibit innovative swimming patterns that contribute to exploring the origins of animal locomotion. However, the genetic and cellular basis of these patterns remains unclear. Herein, we generated chromosome-level genome assemblies of two jellyfish species, Turritopsis rubra and Aurelia coerulea, which exhibit straight and free-swimming patterns, respectively. We observe positive selection of numerous genes involved in statolith formation, hair cell ciliogenesis, ciliary motility, and motor neuron function. The lineage-specific absence of otolith morphogenesis- and ciliary movement-related genes in T. rubra may be associated with homeostatic structural statocyst loss and straight swimming pattern. Notably, single-cell transcriptomic analyses covering key developmental stages reveal the enrichment of diapause-related genes in the cyst during reverse development, suggesting that the sustained diapause state favours the development of new polyps under favourable conditions. This study highlights the complex relationship between genetics, locomotion patterns and survival strategies in jellyfish, thereby providing valuable insights into the evolutionary lineages of movement and adaptation in the animal kingdom.

Integrated analysis of single-cell and bulk RNA-sequencing identifies a signature based on drug response genes to predict prognosis and therapeutic response in ovarian cancer

Ovarian cancer represents a severe gynecological malignancy with a dire prognosis, underscoring the imperative need for dependable biomarkers that can accurately predict drug response and guide therapeutic choices. In this study, we harnessed online single-cell RNA sequencing (scRNAseq) and bulk RNA sequencing (RNAseq) datasets, applying the Scissor algorithm to identify cells responsive to paclitaxel. From these cells, we derived a gene signature, subsequently used to construct a prognostic model that demonstrated high sensitivity and specificity in predicting patient outcomes. Moreover, we conducted pathway and functional enrichment analyses to uncover potential molecular mechanisms driving the prognostic gene signature. This study illustrates the critical role of scRNAseq and bulk RNAseq in developing precise prognostic models for ovarian cancer, potentially transforming clinical decision-making.

Msh3 and Pms1 Set Neuronal CAG-repeat Migration Rate to Drive Selective Striatal and Cortical Pathogenesis in HD Mice

Modifiers of Huntington's disease (HD) include mismatch repair (MMR) genes; however, their underlying disease-altering mechanisms remain unresolved. Knockout (KO) alleles for 9 HD GWAS modifiers/MMR genes were crossed to the Q140 Huntingtin (mHtt) knock-in mice to probe such mechanisms. Four KO mice strongly ( Msh3 and Pms1 ) or moderately ( Msh2 and Mlh1 ) rescue a triad of adult-onset, striatal medium-spiny-neuron (MSN)-selective phenotypes: somatic Htt DNA CAG-repeat expansion, transcriptionopathy, and mHtt protein aggregation. Comparatively, Q140 cortex also exhibits an analogous, but later-onset, pathogenic triad that is Msh3 -dependent. Remarkably, Q140/homozygous Msh3-KO lacks visible mHtt aggregates in the brain, even at advanced ages (20-months). Moreover, Msh3 -deficiency prevents striatal synaptic marker loss, astrogliosis, and locomotor impairment in HD mice. Purified Q140 MSN nuclei exhibit highly linear age-dependent mHtt DNA repeat expansion (i.e. repeat migration), with modal-CAG increasing at +8.8 repeats/month (R 2 =0.98). This linear rate is reduced to 2.3 and 0.3 repeats/month in Q140 with Msh3 heterozygous and homozygous alleles, respectively. Our study defines somatic Htt CAG-repeat thresholds below which there are no detectable mHtt nuclear or neuropil aggregates. Mild transcriptionopathy can still occur in Q140 mice with stabilized Htt 140-CAG repeats, but the majority of transcriptomic changes are due to somatic repeat expansion. Our analysis reveals 479 genes with expression levels highly correlated with modal-CAG length in MSNs. Thus, our study mechanistically connects HD GWAS genes to selective neuronal vulnerability in HD, in which Msh3 and Pms1 set the linear rate of neuronal mHtt CAG-repeat migration to drive repeat-length dependent pathogenesis; and provides a preclinical platform for targeting these genes for HD suppression across brain regions.

One Sentence Summary

Msh3 and Pms1 are genetic drivers of sequential striatal and cortical pathogenesis in Q140 mice by mediating selective CAG-repeat migration in HD vulnerable neurons.

Unraveling the molecular architecture of autoimmune thyroid diseases at spatial resolution

Autoimmune thyroid diseases (AITD) such as Graves' disease (GD) or Hashimoto's thyroiditis (HT) are organ-specific diseases that involve complex interactions between distinct components of thyroid tissue. Here, we use spatial transcriptomics to explore the molecular architecture, heterogeneity and location of different cells present in the thyroid tissue, including thyroid follicular cells (TFCs), stromal cells such as fibroblasts, endothelial cells, and thyroid infiltrating lymphocytes. We identify damaged antigen-presenting TFCs with upregulated CD74 and MIF expression in thyroid samples from AITD patients. Furthermore, we discern two main fibroblast subpopulations in the connective tissue including ADIRF+ myofibroblasts, mainly enriched in GD, and inflammatory fibroblasts, enriched in HT patients. We also demonstrate an increase of fenestrated PLVAP+ vessels in AITD, especially in GD. Our data unveil stromal and thyroid epithelial cell subpopulations that could play a role in the pathogenesis of AITD.

Synthetic DNA barcodes identify singlets in scRNA-seq datasets and evaluate doublet algorithms

Single-cell RNA sequencing (scRNA-seq) datasets contain true single cells, or singlets, in addition to cells that coalesce during the protocol, or doublets. Identifying singlets with high fidelity in scRNA-seq is necessary to avoid false negative and false positive discoveries. Although several methodologies have been proposed, they are typically tested on highly heterogeneous datasets and lack a priori knowledge of true singlets. Here, we leveraged datasets with synthetically introduced DNA barcodes for a hitherto unexplored application: to extract ground-truth singlets. We demonstrated the feasibility of our framework, "singletCode," to evaluate existing doublet detection methods across a range of contexts. We also leveraged our ground-truth singlets to train a proof-of-concept machine learning classifier, which outperformed other doublet detection algorithms. Our integrative framework can identify ground-truth singlets and enable robust doublet detection in non-barcoded datasets.

Single cell transcriptomes and multiscale networks from persons with and without Alzheimer's disease

The emergence of single nucleus RNA sequencing (snRNA-seq) offers to revolutionize the study of Alzheimer's disease (AD). Integration with complementary multiomics data such as genetics, proteomics and clinical data provides powerful opportunities to link cell subpopulations and molecular networks with a broader disease-relevant context. We report snRNA-seq profiles from superior frontal gyrus samples from 101 well characterized subjects from the Banner Brain and Body Donation Program in combination with whole genome sequences. We report findings that link common AD risk variants with CR1 expression in oligodendrocytes as well as alterations in hematological parameters. We observed an AD-associated CD83(+) microglial subtype with unique molecular networks and which is associated with immunoglobulin IgG4 production in the transverse colon. Our major observations were replicated in two additional, independent snRNA-seq data sets. These findings illustrate the power of multi-tissue molecular profiling to contextualize snRNA-seq brain transcriptomics and reveal disease biology.

Genome-scale metabolic network of human carotid plaque reveals the pivotal role of glutamine/glutamate metabolism in macrophage modulating plaque inflammation and vulnerability

BACKGROUND

Metabolism is increasingly recognized as a key regulator of the function and phenotype of the primary cellular constituents of the atherosclerotic vascular wall, including endothelial cells, smooth muscle cells, and inflammatory cells. However, a comprehensive analysis of metabolic changes associated with the transition of plaque from a stable to a hemorrhaged phenotype is lacking.

METHODS

In this study, we integrated two large mRNA expression and protein abundance datasets (BIKE, n = 126; MaasHPS, n = 43) from human atherosclerotic carotid artery plaque to reconstruct a genome-scale metabolic network (GEM). Next, the GEM findings were linked to metabolomics data from MaasHPS, providing a comprehensive overview of metabolic changes in human plaque.

RESULTS

Our study identified significant changes in lipid, cholesterol, and inositol metabolism, along with altered lysosomal lytic activity and increased inflammatory activity, in unstable plaques with intraplaque hemorrhage (IPH+) compared to non-hemorrhaged (IPH-) plaques. Moreover, topological analysis of this network model revealed that the conversion of glutamine to glutamate and their flux between the cytoplasm and mitochondria were notably compromised in hemorrhaged plaques, with a significant reduction in overall glutamate levels in IPH+ plaques. Additionally, reduced glutamate availability was associated with an increased presence of macrophages and a pro-inflammatory phenotype in IPH+ plaques, suggesting an inflammation-prone microenvironment.

CONCLUSIONS

This study is the first to establish a robust and comprehensive GEM for atherosclerotic plaque, providing a valuable resource for understanding plaque metabolism. The utility of this GEM was illustrated by its ability to reliably predict dysregulation in the cholesterol hydroxylation, inositol metabolism, and the glutamine/glutamate pathway in rupture-prone hemorrhaged plaques, a finding that may pave the way to new diagnostic or therapeutic measures.

Pathophysiology of hypereosinophilia-associated heart disease

Background

Cardiac complications in patients with hypereosinophilia cause significant morbidity and mortality. However, mechanisms of how eosinophilic inflammation causes heart damage are poorly understood.

Methods

We developed a model of hypereosinophilia-associated heart disease by challenging hypereosinophilic mice with peptide from the cardiac myosin heavy chain. Disease outcomes were measured by histology, immunohistochemistry, flow cytometry, and measurement of cells and biomarkers in peripheral blood. Eosinophil dependence was determined by using eosinophil-deficient mice (ΔdblGATA). Single cells from heart were subjected to single cell RNA sequencing to assess cell composition, subtypes and expression profiles.

Results

Mice challenged with myocarditic and control peptide had peripheral blood leukocytosis, but only those challenged with myocarditic peptide had heart inflammation. Heart tissue was infiltrated by eosinophil-rich inflammatory infiltrates associated with cardiomyocyte damage. Disease penetrance and severity were dependent on the presence of eosinophils. Single cell RNA sequencing showed enrichment of myeloid cells, T-cells and granulocytes (neutrophils and eosinophils) in the myocarditic mice. Macrophages were M2 skewed, and eosinophils had an activated phenotype. Gene enrichment analysis identified several pathways potentially involved in pathophysiology of disease.

Conclusion

Eosinophils are required for heart damage in hypereosinophilia-associated heart disease. Additionally, myeloid cells, granulocytes and T-cell cooperatively or independently participate in the pathogenesis of hypereosinophilia-associated heart disease.

Human yolk sac-derived innate lymphoid-biased multipotent progenitors emerge prior to hematopoietic stem cell formation

Hematopoietic stem cell (HSC)-independent lymphopoiesis has been elucidated in murine embryos. However, our understanding regarding human embryonic counterparts remains limited. Here, we demonstrated the presence of human yolk sac-derived lymphoid-biased progenitors (YSLPs) expressing CD34, IL7R, LTB, and IRF8 at Carnegie stage 10, much earlier than the first HSC emergence. The number and lymphopoietic potential of these progenitors were both significantly higher in the yolk sac than the embryo proper at this early stage. Importantly, single-cell/bulk culture and CITE-seq have elucidated the tendency of YSLP to differentiate into innate lymphoid cells and dendritic cells. Notably, lymphoid progenitors in fetal liver before and after HSC seeding displayed distinct transcriptional features, with the former closely resembling those of YSLPs. Overall, our data identified the origin, potential, and migratory dynamics of innate lymphoid-biased multipotent progenitors in human yolk sac before HSC emergence, providing insights for understanding the stepwise establishment of innate immune system in humans.

Control of innate olfactory valence by segregated cortical amygdala circuits

Animals perform innate behaviors that are stereotyped responses to specific evolutionarily relevant stimuli in the absence of prior learning or experience. These behaviors can be reduced to an axis of valence, whereby specific odors evoke approach or avoidance. The cortical amygdala (plCoA) mediates innate attraction and aversion to odor. However, little is known about how this brain area gives rise to behaviors of opposing motivational valence. Here, we sought to define the circuit features of plCoA that give rise to innate olfactory behaviors of valence. We characterized the physiology, gene expression, and projections of this structure, identifying a divergent, topographic organization that selectively controls innate attraction and avoidance to odor. First, we examined odor-evoked responses in these areas and found sparse encoding of odor identity, but not valence. We next considered a topographic organization and found that optogenetic stimulation of the anterior and posterior domains of plCoA elicits attraction and avoidance, respectively, suggesting a functional axis for valence. Using single cell and spatial RNA sequencing, we identified the molecular cell types in plCoA, revealing an anteroposterior gradient in cell types, whereby anterior glutamatergic neurons preferentially express Slc17a6 and posterior neurons express Slc17a7. Activation of these respective cell types recapitulates appetitive and aversive valence behaviors, and chemogenetic inhibition reveals partial necessity for valence responses to innate appetitive or aversive odors. Finally, we identified topographically organized circuits defined by projections, whereby anterior neurons preferentially project to medial amygdala, and posterior neurons preferentially project to nucleus accumbens, which are respectively sufficient and necessary for innate negative and positive olfactory valence. Together, these data advance our understanding of how the olfactory system generates stereotypic, hardwired attraction and avoidance, and supports a model whereby distinct, topographically distributed plCoA populations direct innate olfactory valence responses by signaling to divergent valence-specific targets, linking upstream olfactory identity to downstream valence behaviors, through a population code. This represents a novel circuit motif in which valence encoding is represented not by the firing properties of individual neurons, but by population level identity encoding that is routed through divergent targets to mediate distinct valence.

Hypoxia inhibits the iMo/cDC2/CD8+ TRMs immune axis in the tumor microenvironment of human esophageal cancer

BACKGROUND

Esophageal cancer (ESCA) is a form of malignant tumor associated with chronic inflammation and immune dysregulation. However, the specific immune status and key mechanisms of immune regulation in this disease require further exploration.

METHODS

To investigate the features of the human ESCA tumor immune microenvironment and its possible regulation, we performed mass cytometry by time of flight, single-cell RNA sequencing, multicolor fluorescence staining of tissue, and flow cytometry analyses on tumor and paracancerous tissue from treatment-naïve patients.

RESULTS

We depicted the immune landscape of the ESCA and revealed that CD8+ (tissue-resident memory CD8+ T cells (CD8+ TRMs) were closely related to disease progression. We also revealed the heterogeneity of CD8+ TRMs in the ESCA tumor microenvironment (TME), which was associated with their differentiation and function. Moreover, the subset of CD8+ TRMs in tumor (called tTRMs) that expressed high levels of granzyme B and immune checkpoints was markedly decreased in the TME of advanced ESCA. We showed that tTRMs are tumor effector cells preactivated in the TME. We then demonstrated that conventional dendritic cells (cDC2s) derived from intermediate monocytes (iMos) are essential for maintaining the proliferation of CD8+ TRMs in the TME. Our preliminary study showed that hypoxia can promote the apoptosis of iMos and impede the maturation of cDC2s, which in turn reduces the proliferative capacity of CD8+ TRMs, thereby contributing to the progression of cancer.

CONCLUSIONS

Our study revealed the essential antitumor roles of CD8+ TRMs and preliminarily explored the regulation of the iMo/cDC2/CD8+ TRM immune axis in the human ESCA TME.

Integrating single-nucleus RNA sequencing and spatial transcriptomics to elucidate a specialized subpopulation of astrocytes, microglia and vascular cells in brains of mouse model of lipopolysaccharide-induced sepsis-associated encephalopathy

BACKGROUND

Understanding the mechanism behind sepsis-associated encephalopathy (SAE) remains a formidable task. This study endeavors to shed light on the complex cellular and molecular alterations that occur in the brains of a mouse model with SAE, ultimately unraveling the underlying mechanisms of this condition.

METHODS

We established a murine model using intraperitoneal injection of lipopolysaccharide (LPS) in wild type and Anxa1-/- mice and collected brain tissues for analysis at 0-hour, 12-hour, 24-hour, and 72-hour post-injection. Utilizing advanced techniques such as single-nucleus RNA sequencing (snRNA-seq) and Stereo-seq, we conducted a comprehensive characterization of the cellular responses and molecular patterns within the brain.

RESULTS

Our study uncovered notable temporal differences in the response to LPS challenge between Anxa1-/- (annexin A1 knockout) and wild type mice, specifically at the 12-hour and 24-hour time points following injection. We observed a significant increase in the proportion of Astro-2 and Micro-2 cells in these mice. These cells exhibited a colocalization pattern with the vascular subtype Vas-1, forming a distinct region known as V1A2M2, where Astro-2 and Micro-2 cells surrounded Vas-1. Moreover, through further analysis, we discovered significant upregulation of ligands and receptors such as Timp1-Cd63, Timp1-Itgb1, Timp1-Lrp1, as well as Ccl2-Ackr1 and Cxcl2-Ackr1 within this region. In addition, we observed a notable increase in the expression of Cd14-Itgb1, Cd14-Tlr2, and Cd14-C3ar1 in regions enriched with Micro-2 cells. Additionally, Cxcl10-Sdc4 showed broad upregulation in brain regions containing both Micro-2 and Astro-2 cells. Notably, upon LPS challenge, there was an observed increase in Anxa1 expression in the mouse brain. Furthermore, our study revealed a noteworthy increase in mortality rates following Anxa1 knockdown. However, we did not observe substantial differences in the types, numbers, or distribution of other brain cells between Anxa1-/- and wildtype mice over time. Nevertheless, when comparing the 24-hour post LPS injection time point, we observed a significant decrease in the proportion and distribution of Micro-2 and Astro-2 cells in the vicinity of blood vessels in Anxa1-/- mice. Additionally, we noted reduced expression levels of several ligand-receptor pairs including Cd14-Tlr2, Cd14-C3ar1, Cd14-Itgb1, Cxcl10-Sdc4, Ccl2-Ackr1, and Cxcl2-Ackr1.

CONCLUSIONS

By combining snRNA-seq and Stereo-seq techniques, our study successfully identified a distinctive cellular colocalization, referred to as a special pathological niche, comprising Astro-2, Micro-2, and Vas-1 cells. Furthermore, we observed an upregulation of ligand-receptor pairs within this niche. These findings suggest a potential association between this cellular arrangement and the underlying mechanisms contributing to SAE or the increased mortality observed in Anxa1 knockdown mice.

Single cell regulatory architecture of human pancreatic islets suggests sex differences in β cell function and the pathogenesis of type 2 diabetes

Type 2 and type 1 diabetes (T2D, T1D) exhibit sex differences in insulin secretion, the mechanisms of which are unknown. We examined sex differences in human pancreatic islets from 52 donors with and without T2D combining single cell RNA-seq (scRNA-seq), single nucleus assay for transposase-accessible chromatin sequencing (snATAC-seq), hormone secretion, and bioenergetics. In nondiabetic (ND) donors, sex differences in islet cells gene accessibility and expression predominantly involved sex chromosomes. Islets from T2D donors exhibited similar sex differences in sex chromosomes differentially expressed genes (DEGs), but also exhibited sex differences in autosomal genes. Comparing β cells from T2D vs. ND donors, gene enrichment of female β cells showed suppression in mitochondrial respiration, while male β cells exhibited suppressed insulin secretion. Thus, although sex differences in gene accessibility and expression of ND β cells predominantly affect sex chromosomes, the transition to T2D reveals sex differences in autosomes highlighting mitochondrial failure in females.

Metabolic Reprogramming of Tumor-Associated Macrophages Using Glutamine Antagonist JHU083 Drives Tumor Immunity in Myeloid-Rich Prostate and Bladder Cancers

Glutamine metabolism in tumor microenvironments critically regulates antitumor immunity. Using the glutamine-antagonist prodrug JHU083, we report potent tumor growth inhibition in urologic tumors by JHU083-reprogrammed tumor-associated macrophages (TAMs) and tumor-infiltrating monocytes. We show JHU083-mediated glutamine antagonism in tumor microenvironments induced by TNF, proinflammatory, and mTORC1 signaling in intratumoral TAM clusters. JHU083-reprogrammed TAMs also exhibited increased tumor cell phagocytosis and diminished proangiogenic capacities. In vivo inhibition of TAM glutamine consumption resulted in increased glycolysis, a broken tricarboxylic acid (TCA) cycle, and purine metabolism disruption. Although the antitumor effect of glutamine antagonism on tumor-infiltrating T cells was moderate, JHU083 promoted a stem cell-like phenotype in CD8+ T cells and decreased the abundance of regulatory T cells. Finally, JHU083 caused a global shutdown in glutamine-utilizing metabolic pathways in tumor cells, leading to reduced HIF-1α, c-MYC phosphorylation, and induction of tumor cell apoptosis, all key antitumor features. Altogether, our findings demonstrate that targeting glutamine with JHU083 led to suppressed tumor growth as well as reprogramming of immunosuppressive TAMs within prostate and bladder tumors that promoted antitumor immune responses. JHU083 can offer an effective therapeutic benefit for tumor types that are enriched in immunosuppressive TAMs.

A unified model-based framework for doublet or multiplet detection in single-cell multiomics data

Droplet-based single-cell sequencing techniques rely on the fundamental assumption that each droplet encapsulates a single cell, enabling individual cell omics profiling. However, the inevitable issue of multiplets, where two or more cells are encapsulated within a single droplet, can lead to spurious cell type annotations and obscure true biological findings. The issue of multiplets is exacerbated in single-cell multiomics settings, where integrating cross-modality information for clustering can inadvertently promote the aggregation of multiplet clusters and increase the risk of erroneous cell type annotations. Here, we propose a compound Poisson model-based framework for multiplet detection in single-cell multiomics data. Leveraging experimental cell hashing results as the ground truth for multiplet status, we conducted trimodal DOGMA-seq experiments and generated 17 benchmarking datasets from two tissues, involving a total of 280,123 droplets. We demonstrated that the proposed method is an essential tool for integrating cross-modality multiplet signals, effectively eliminating multiplet clusters in single-cell multiomics data-a task at which the benchmarked single-omics methods proved inadequate.

Investigating the characteristics of mild intervertebral disc degeneration at various age stages using single-cell genomics

Introduction: Intervertebral disc degeneration often occurs in the elderly population, but in recent years, there has been an increasing incidence of disc degeneration in younger individuals, primarily with mild degeneration. Methods: In order to explore the underlying mechanisms of disc degeneration in both young and aging individuals, we collected four types of nucleus pulposus (NP) single-cell sequencing samples for analysis based on Pfirrmann grading: normal-young (NY) (Grade I), normal-old (NO) (Grade I), mild degenerative-young (MY) (Grade II-III), and mild degenerative-old (MO) (Grade II-III). Results: We found that most NP cells in NO and MY samples exhibited oxidative stress, which may be important pathogenic factors in NO and MY groups. On the other hand, NP cells in MO group exhibited endoplasmic reticulum stress. In terms of inflammation, myeloid cells were mainly present in the degenerative group, with the MY group showing a stronger immune response compared to the MO group. Interestingly, dendritic cells in the myeloid lineage played a critical role in the process of mild degeneration. Discussion: Our study investigated the molecular mechanisms of intervertebral disc degeneration from an age perspective, providing insights for improving treatment strategies for patients with disc degeneration at different age groups.

Complexin-1 enhances ultrasound neurotransmission in the mammalian auditory pathway

Unlike megabats, which rely on well-developed vision, microbats use ultrasonic echolocation to navigate and locate prey. To study ultrasound perception, here we compared the auditory cortices of microbats and megabats by constructing reference genomes and single-nucleus atlases for four species. We found that parvalbumin (PV)+ neurons exhibited evident cross-species differences and could respond to ultrasound signals, whereas their silencing severely affected ultrasound perception in the mouse auditory cortex. Moreover, megabat PV+ neurons expressed low levels of complexins (CPLX1-CPLX4), which can facilitate neurotransmitter release, while microbat PV+ neurons highly expressed CPLX1, which improves neurotransmission efficiency. Further perturbation of Cplx1 in PV+ neurons impaired ultrasound perception in the mouse auditory cortex. In addition, CPLX1 functioned in other parts of the auditory pathway in microbats but not megabats and exhibited convergent evolution between echolocating microbats and whales. Altogether, we conclude that CPLX1 expression throughout the entire auditory pathway can enhance mammalian ultrasound neurotransmission.

Reduced adipocyte glutaminase activity promotes energy expenditure and metabolic health

Glutamine and glutamate are interconverted by several enzymes and alterations in this metabolic cycle are linked to cardiometabolic traits. Herein, we show that obesity-associated insulin resistance is characterized by decreased plasma and white adipose tissue glutamine-to-glutamate ratios. We couple these stoichiometric changes to perturbed fat cell glutaminase and glutamine synthase messenger RNA and protein abundance, which together promote glutaminolysis. In human white adipocytes, reductions in glutaminase activity promote aerobic glycolysis and mitochondrial oxidative capacity via increases in hypoxia-inducible factor 1α abundance, lactate levels and p38 mitogen-activated protein kinase signalling. Systemic glutaminase inhibition in male and female mice, or genetically in adipocytes of male mice, triggers the activation of thermogenic gene programs in inguinal adipocytes. Consequently, the knockout mice display higher energy expenditure and improved glucose tolerance compared to control littermates, even under high-fat diet conditions. Altogether, our findings highlight white adipocyte glutamine turnover as an important determinant of energy expenditure and metabolic health.

Influence of METTL3 knockdown on PDLSC osteogenesis in E. coli LPS-induced inflammation

OBJECTIVE

This study aimed to investigate the effect of METTL3 knockdown on osteogenic differentiation of human periodontal ligament stem cells (PDLSCs) in the weak inflammation microenvironments, as well as the underlying mechanisms.

MATERIALS AND METHODS

PDLSCs were stimulated by lipopolysaccharide from Escherichia coli (E. coli LPS), followed by quantification of METTL3. METTL3 expression was assessed using RT-qPCR and Western blot analysis in periodontitis. METTL3 knockdown PDLSCs were stimulated with or without E. coli LPS. The evaluation included proinflammatory cytokines, osteogenic markers, ALP activity, and mineralized nodules. Bioinformatics analysis and Western blot determined the association between METTL3 and the PI3K/Akt pathway.

RESULTS

METTL3 was overexpressed in periodontitis. METTL3 knockdown in PDLSCs reduced proinflammatory cytokines, osteogenic markers, ALP activity, and mineralized nodules in both environments. Bioinformatics analysis suggested a link between METTL3 and the PI3K/Akt pathway. METTL3 knockdown inhibited PI3K/Akt signaling pathway activation.

CONCLUSION

METTL3 knockdown might inhibit osteogenesis in PDLSCs through the inactivation of PI3K/Akt signaling pathway. Concomitant findings might shed novel light on the roles and potential mechanisms of METTL3 in the LPS-stimulated inflammatory microenvironments of PDLSCs.

Single-cell and spatial atlases of spinal cord injury in the Tabulae Paralytica

Here, we introduce the Tabulae Paralytica-a compilation of four atlases of spinal cord injury (SCI) comprising a single-nucleus transcriptome atlas of half a million cells, a multiome atlas pairing transcriptomic and epigenomic measurements within the same nuclei, and two spatial transcriptomic atlases of the injured spinal cord spanning four spatial and temporal dimensions. We integrated these atlases into a common framework to dissect the molecular logic that governs the responses to injury within the spinal cord1. The Tabulae Paralytica uncovered new biological principles that dictate the consequences of SCI, including conserved and divergent neuronal responses to injury; the priming of specific neuronal subpopulations to upregulate circuit-reorganizing programs after injury; an inverse relationship between neuronal stress responses and the activation of circuit reorganization programs; the necessity of re-establishing a tripartite neuroprotective barrier between immune-privileged and extra-neural environments after SCI and a failure to form this barrier in old mice. We leveraged the Tabulae Paralytica to develop a rejuvenative gene therapy that re-established this tripartite barrier, and restored the natural recovery of walking after paralysis in old mice. The Tabulae Paralytica provides a window into the pathobiology of SCI, while establishing a framework for integrating multimodal, genome-scale measurements in four dimensions to study biology and medicine.

Revealing the molecular landscape of human placenta: a systematic review and meta-analysis of single-cell RNA sequencing studies

BACKGROUND

With increasing significance of developmental programming effects associated with placental dysfunction, more investigations are devoted to improving the characterization and understanding of placental signatures in health and disease. The placenta is a transitory but dynamic organ adapting to the shifting demands of fetal development and available resources of the maternal supply throughout pregnancy. Trophoblasts (cytotrophoblasts, syncytiotrophoblasts, and extravillous trophoblasts) are placental-specific cell types responsible for the main placental exchanges and adaptations. Transcriptomic studies with single-cell resolution have led to advances in understanding the placenta's role in health and disease. These studies, however, often show discrepancies in characterization of the different placental cell types.

OBJECTIVE AND RATIONALE

We aim to review the knowledge regarding placental structure and function gained from the use of single-cell RNA sequencing (scRNAseq), followed by comparing cell-type-specific genes, highlighting their similarities and differences. Moreover, we intend to identify consensus marker genes for the various trophoblast cell types across studies. Finally, we will discuss the contributions and potential applications of scRNAseq in studying pregnancy-related diseases.

SEARCH METHODS

We conducted a comprehensive systematic literature review to identify different cell types and their functions at the human maternal-fetal interface, focusing on all original scRNAseq studies on placentas published before March 2023 and published reviews (total of 28 studies identified) using PubMed search. Our approach involved curating cell types and subtypes that had previously been defined using scRNAseq and comparing the genes used as markers or identified as potential new markers. Next, we reanalyzed expression matrices from the six available scRNAseq raw datasets with cell annotations (four from first trimester and two at term), using Wilcoxon rank-sum tests to compare gene expression among studies and annotate trophoblast cell markers in both first trimester and term placentas. Furthermore, we integrated scRNAseq raw data available from 18 healthy first trimester and nine term placentas, and performed clustering and differential gene expression analysis. We further compared markers obtained with the analysis of annotated and raw datasets with the literature to obtain a common signature gene list for major placental cell types.

OUTCOMES

Variations in the sampling site, gestational age, fetal sex, and subsequent sequencing and analysis methods were observed between the studies. Although their proportions varied, the three trophoblast types were consistently identified across all scRNAseq studies, unlike other non-trophoblast cell types. Notably, no marker genes were shared by all studies for any of the investigated cell types. Moreover, most of the newly defined markers in one study were not observed in other studies. These discrepancies were confirmed by our analysis on trophoblast cell types, where hundreds of potential marker genes were identified in each study but with little overlap across studies. From 35 461 and 23 378 cells of high quality in the first trimester and term placentas, respectively, we obtained major placental cell types, including perivascular cells that previously had not been identified in the first trimester. Importantly, our meta-analysis provides marker genes for major placental cell types based on our extensive curation.

WIDER IMPLICATIONS

This review and meta-analysis emphasizes the need for establishing a consensus for annotating placental cell types from scRNAseq data. The marker genes identified here can be deployed for defining human placental cell types, thereby facilitating and improving the reproducibility of trophoblast cell annotation.