Urine scRNAseq reveals new insights into the bladder tumor immune microenvironment.

Due to bladder tumors' contact with urine, urine-derived cells (UDCs) may serve as a surrogate for monitoring the tumor microenvironment (TME) in bladder cancer (BC). However, the composition of UDCs and the extent to which they mirror the tumor remain poorly characterized. We generated the first single-cell RNA-sequencing of BC patient UDCs with matched tumor and peripheral blood mononuclear cells (PBMC). BC urine was more cellular than healthy donor (HD) urine, containing multiple immune populations including myeloid cells, CD4+ and CD8+ T cells, natural killer (NK) cells, B cells, and dendritic cells (DCs) in addition to tumor and stromal cells. Immune UDCs were transcriptionally more similar to tumor than blood. UDCs encompassed cytotoxic and activated CD4+ T cells, exhausted and tissue-resident memory CD8+ T cells, macrophages, germinal-center-like B cells, tissue-resident and adaptive NK cells, and regulatory DCs found in tumor but lacking or absent in blood. Our findings suggest BC UDCs may be surrogates for the TME and serve as therapeutic biomarkers.

A master regulator of opioid reward in the ventral prefrontal cortex.

In addition to their intrinsic rewarding properties, opioids can also evoke aversive reactions that protect against misuse. Cellular mechanisms that govern the interplay between opioid reward and aversion are poorly understood. We used whole-brain activity mapping in mice to show that neurons in the dorsal peduncular nucleus (DPn) are highly responsive to the opioid oxycodone. Connectomic profiling revealed that DPn neurons innervate the parabrachial nucleus (PBn). Spatial and single-nuclei transcriptomics resolved a population of PBn-projecting pyramidal neurons in the DPn that express µ-opioid receptors (µORs). Disrupting µOR signaling in the DPn switched oxycodone from rewarding to aversive and exacerbated the severity of opioid withdrawal. These findings identify the DPn as a key substrate for the abuse liability of opioids.

Single-cell transcriptomic-informed deconvolution of bulk data identifies immune checkpoint blockade resistance in urothelial cancer.

Interactions within the tumor microenvironment (TME) significantly influence tumor progression and treatment responses. While single-cell RNA sequencing (scRNA-seq) and spatial genomics facilitate TME exploration, many clinical cohorts are assessed at the bulk tissue level. Integrating scRNA-seq and bulk tissue RNA-seq data through computational deconvolution is essential for obtaining clinically relevant insights. Our method, ProM, enables the examination of major and minor cell types. Through evaluation against existing methods using paired single-cell and bulk RNA sequencing of human urothelial cancer (UC) samples, ProM demonstrates superiority. Application to UC cohorts treated with immune checkpoint inhibitors reveals pre-treatment cellular features associated with poor outcomes, such as elevated SPP1 expression in macrophage/monocytes (MM). Our deconvolution method and paired single-cell and bulk tissue RNA-seq dataset contribute novel insights into TME heterogeneity and resistance to immune checkpoint blockade.

High-Resolution Profiling of Human Vocal Fold Cellular Landscapes With Single-Nuclei RNA Sequencing.

INTRODUCTION: The function of the vocal folds (VFs) is determined by the phenotype, abundance, and distribution of differentiated cells within specific microenvironments. Identifying this histologic framework is crucial in understanding laryngeal disease. A paucity of studies investigating VF cellular heterogeneity has been undertaken. Here, we examined the cellular landscape of human VFs by utilizing single-nuclei RNA-sequencing. METHODS: Normal true VF tissue was excised from five patients undergoing pitch elevation surgery. Tissue was snap frozen in liquid nitrogen and subjected to cellular digestion and nuclear extraction. Nuclei were processed for single-nucleus sequencing using the 10X Genomics Chromium platform. Sequencing reads were assembled using cellranger and analyzed with the scanpy package in python. RESULTS: RNA sequencing revealed 18 global cell clusters. While many were of epithelial origin, expected cell types, such as fibroblasts, immune cells, muscle cells, and endothelial cells were present. Subcluster analysis defined unique epithelial, immune, and fibroblast subpopulations. CONCLUSION: This study evaluated the cellular heterogeneity of normal human VFs by utilizing single-nuclei RNA-sequencing. With further confirmation through additional spatial sequencing and microscopic imaging, a novel cellular map of the VFs may provide insight into new cellular targets for VF disease. LEVEL OF EVIDENCE: NA Laryngoscope, 134:3193-3200, 2024.

Molecular profiling of human substantia nigra identifies diverse neuron types associated with vulnerability in Parkinson's disease.

Parkinson's disease (PD) is characterized pathologically by the loss of dopaminergic (DA) neurons in the substantia nigra (SN). Whether cell types beyond DA neurons in the SN show vulnerability in PD remains unclear. Through transcriptomic profiling of 315,867 high-quality single nuclei in the SN from individuals with and without PD, we identified cell clusters representing various neuron types, glia, endothelial cells, pericytes, fibroblasts, and T cells and investigated cell type-dependent alterations in gene expression in PD. Notably, a unique neuron cluster marked by the expression of RIT2, a PD risk gene, also displayed vulnerability in PD. We validated RIT2-enriched neurons in midbrain organoids and the mouse SN. Our results demonstrated distinct transcriptomic signatures of the RIT2-enriched neurons in the human SN and implicated reduced RIT2 expression in the pathogenesis of PD. Our study sheds light on the diversity of cell types, including DA neurons, in the SN and the complexity of molecular and cellular changes associated with PD pathogenesis.

Nuclear RNA catabolism controls endogenous retroviruses, gene expression asymmetry, and dedifferentiation.

Endogenous retroviruses (ERVs) are remnants of ancient parasitic infections and comprise sizable portions of most genomes. Although epigenetic mechanisms silence most ERVs by generating a repressive environment that prevents their expression (heterochromatin), little is known about mechanisms silencing ERVs residing in open regions of the genome (euchromatin). This is particularly important during embryonic development, where induction and repression of distinct classes of ERVs occur in short temporal windows. Here, we demonstrate that transcription-associated RNA degradation by the nuclear RNA exosome and Integrator is a regulatory mechanism that controls the productive transcription of most genes and many ERVs involved in preimplantation development. Disrupting nuclear RNA catabolism promotes dedifferentiation to a totipotent-like state characterized by defects in RNAPII elongation and decreased expression of long genes (gene-length asymmetry). Our results indicate that RNA catabolism is a core regulatory module of gene networks that safeguards RNAPII activity, ERV expression, cell identity, and developmental potency.

Sexually dimorphic role of diet and stress on behavior, energy metabolism, and the ventromedial hypothalamus.

Scientific evidence underscores the influence of biological sex on the interplay between stress and metabolic dysfunctions. However, there is limited understanding of how diet and stress jointly contribute to metabolic dysregulation in both males and females. To address this gap, our study aimed to investigate the combined effects of a high-fat diet (HFD) and repeated footshock stress on fear-related behaviors and metabolic outcomes in male and female mice. Using a robust rodent model that recapitulates key aspects of post-traumatic stress disorder (PTSD), we subjected mice to footshock stressor followed by weekly reminder footshock stressor or no stressor for 14 weeks while on either an HFD or chow diet. Our findings revealed that HFD impaired fear memory extinction in male mice that received initial stressor but not in female mice. Blood glucose levels were influenced by both diet and sex, with HFD-fed female mice displaying elevated levels that returned to baseline in the absence of stress, a pattern not observed in male mice. Male mice on HFD exhibited higher energy expenditure, while HFD-fed female mice showed a decreased respiratory exchange ratio (RER). Sex-specific alterations in pro-inflammatory markers and abundance of hematopoietic stem cells were observed in chronically stressed mice on an HFD in different peripheral tissues, indicating the manifestation of distinct comorbid disorders. Single-nuclei RNA sequencing of the ventromedial hypothalamus from stressed mice on an HFD provided insights into sex-specific glial cell activation and cell-type-specific transcriptomic changes. In conclusion, our study offers a comprehensive understanding of the intricate interactions between stress, diet, sex, and various physiological and behavioral outcomes, shedding light on a potential brain region coordinating these interactions.

Spinal cord injury regulates circular RNA expression in axons.

Introduction: Neurons transport mRNA and translational machinery to axons for local translation. After spinal cord injury (SCI), de novo translation is assumed to enable neurorepair. Knowledge of the identity of axonal mRNAs that participate in neurorepair after SCI is limited. We sought to identify and understand how axonal RNAs play a role in axonal regeneration. Methods: We obtained preparations enriched in axonal mRNAs from control and SCI rats by digesting spinal cord tissue with cold-active protease (CAP). The digested samples were then centrifuged to obtain a supernatant that was used to identify mRNA expression. We identified differentially expressed genes (DEGS) after SCI and mapped them to various biological processes. We validated the DEGs by RT-qPCR and RNA-scope. Results: The supernatant fraction was highly enriched for mRNA from axons. Using Gene Ontology, the second most significant pathway for all DEGs was axonogenesis. Among the DEGs was Rims2, which is predominately a circular RNA (circRNA) in the CNS. We show that Rims2 RNA within spinal cord axons is circular. We found an additional 200 putative circRNAs in the axonal-enriched fraction. Knockdown in primary rat cortical neurons of the RNA editing enzyme ADAR1, which inhibits formation of circRNAs, significantly increased axonal outgrowth and increased the expression of circRims2. Using Rims2 as a prototype we used Circular RNA Interactome to predict miRNAs that bind to circRims2 also bind to the 3'UTR of GAP-43, PTEN or CREB1, all known regulators of axonal outgrowth. Axonally-translated GAP-43 supports axonal elongation and we detect GAP-43 mRNA in the rat axons by RNAscope. Discussion: By enriching for axonal RNA, we detect SCI induced DEGs, including circRNA such as Rims2. Ablation of ADAR1, the enzyme that regulates circRNA formation, promotes axonal outgrowth of cortical neurons. We developed a pathway model using Circular RNA Interactome that indicates that Rims2 through miRNAs can regulate the axonal translation GAP-43 to regulate axonal regeneration. We conclude that axonal regulatory pathways will play a role in neurorepair.

MacroH2A restricts inflammatory gene expression in melanoma cancer-associated fibroblasts by coordinating chromatin looping.

MacroH2A has established tumour suppressive functions in melanoma and other cancers, but an unappreciated role in the tumour microenvironment. Using an autochthonous, immunocompetent mouse model of melanoma, we demonstrate that mice devoid of macroH2A variants exhibit increased tumour burden compared with wild-type counterparts. MacroH2A-deficient tumours accumulate immunosuppressive monocytes and are depleted of functional cytotoxic T cells, characteristics consistent with a compromised anti-tumour response. Single cell and spatial transcriptomics identify increased dedifferentiation along the neural crest lineage of the tumour compartment and increased frequency and activation of cancer-associated fibroblasts following macroH2A loss. Mechanistically, macroH2A-deficient cancer-associated fibroblasts display increased myeloid chemoattractant activity as a consequence of hyperinducible expression of inflammatory genes, which is enforced by increased chromatin looping of their promoters to enhancers that gain H3K27ac. In summary, we reveal a tumour suppressive role for macroH2A variants through the regulation of chromatin architecture in the tumour stroma with potential implications for human melanoma.

HIV-1 activates oxidative phosphorylation in infected CD4 T cells in a human tonsil explant model.

Introduction: Human immunodeficiency virus type 1 (HIV-1) causes a chronic, incurable infection leading to immune activation and chronic inflammation in people with HIV-1 (PWH), even with virologic suppression on antiretroviral therapy (ART). The role of lymphoid structures as reservoirs for viral latency and immune activation has been implicated in chronic inflammation mechanisms. Still, the specific transcriptomic changes induced by HIV-1 infection in different cell types within lymphoid tissue remain unexplored. Methods: In this study, we utilized human tonsil explants from healthy human donors and infected them with HIV-1 ex vivo. We performed single-cell RNA sequencing (scRNA-seq) to analyze the cell types represented in the tissue and to investigate the impact of infection on gene expression profiles and inflammatory signaling pathways. Results: Our analysis revealed that infected CD4+ T cells exhibited upregulation of genes associated with oxidative phosphorylation. Furthermore, macrophages exposed to the virus but uninfected showed increased expression of genes associated with the NLRP3 inflammasome pathway. Discussion: These findings provide valuable insights into the specific transcriptomic changes induced by HIV-1 infection in different cell types within lymphoid tissue. The activation of oxidative phosphorylation in infected CD4+ T cells and the proinflammatory response in macrophages may contribute to the chronic inflammation observed in PWH despite ART. Understanding these mechanisms is crucial for developing targeted therapeutic strategies to eradicate HIV-1 infection in PWH.

A distal lung organoid model to study interstitial lung disease, viral infection and human lung development.

Organoids have been an exciting advancement in stem cell research. Here we describe a strategy for directed differentiation of human pluripotent stem cells into distal lung organoids. This protocol recapitulates lung development by sequentially specifying human pluripotent stem cells to definitive endoderm, anterior foregut endoderm, ventral anterior foregut endoderm, lung bud organoids and finally lung organoids. The organoids take ~40 d to generate and can be maintained more than 180 d, while progressively maturing up to a stage consistent with the second trimester of human gestation. They are unique because of their branching morphology, the near absence of non-lung endodermal lineages, presence of mesenchyme and capacity to recapitulate interstitial lung diseases. This protocol can be performed by anyone familiar with cell culture techniques, is conducted in serum-free conditions and does not require lineage-specific reporters or enrichment steps. We also provide a protocol for the generation of single-cell suspensions for single-cell RNA sequencing.

Human pluripotent stem cell-derived respiratory airway progenitors generate alveolar epithelial cells and recapitulate features of idiopathic pulmonary fibrosis.

Human lungs contain unique cell populations in distal respiratory airways (RAs). These populations accumulate in patients with lung injury, chronic obstructive pulmonary disease (COPD) and idiopathic pulmonary fibrosis (IPF). Their lineage potentials and roles are unknown, however. As they are absent in rodents, deeper understanding of these cells requires a human in vitro model. Here we report the generation from human pluripotent stem cells (hPSCs) of expandable spheres (induced respiratory airway progenitors (iRAPs)) consisting of all RA-associated cell types. iRAPs could differentiate into type 1 (AT1) and type 2 alveolar (AT2) epithelial cells in defined conditions, showing that alveolar cells can be derived from RAs. iRAPs with deletion of HPS1, which causes pulmonary fibrosis in humans, display defects that are hallmarks of IPF, indicating involvement of intrinsic dysfunction of RA-associated cells in IPF. iRAPs thus provide a model to gain insight into human lung regeneration and into pathogenesis of IPF.

The Breast Cancer Single-Cell Atlas: Defining cellular heterogeneity within model cell lines and primary tumors to inform disease subtype, stemness, and treatment options.

PURPOSE: Breast Cancer (BC) is the most diagnosed cancer in women; however, through significant research, relative survival rates have significantly improved. Despite progress, there remains a gap in our understanding of BC subtypes and personalized treatments. This manuscript characterized cellular heterogeneity in BC cell lines through scRNAseq to resolve variability in subtyping, disease modeling potential, and therapeutic targeting predictions. METHODS: We generated a Breast Cancer Single-Cell Cell Line Atlas (BSCLA) to help inform future BC research. We sequenced over 36,195 cells composed of 13 cell lines spanning the spectrum of clinical BC subtypes and leveraged publicly available data comprising 39,214 cells from 26 primary tumors. RESULTS: Unsupervised clustering identified 49 subpopulations within the cell line dataset. We resolve ambiguity in subtype annotation comparing expression of Estrogen Receptor, Progesterone Receptor, and Human Epidermal Growth Factor Receptor 2 genes. Gene correlations with disease subtype highlighted S100A7 and MUCL1 overexpression in HER2 + cells as possible cell motility and localization drivers. We also present genes driving populational drifts to generate novel gene vectors characterizing each subpopulation. A global Cancer Stem Cell (CSC) scoring vector was used to identify stemness potential for subpopulations and model multi-potency. Finally, we overlay the BSCLA dataset with FDA-approved targets to identify to predict the efficacy of subpopulation-specific therapies. CONCLUSION: The BSCLA defines the heterogeneity within BC cell lines, enhancing our overall understanding of BC cellular diversity to guide future BC research, including model cell line selection, unintended sample source effects, stemness factors between cell lines, and cell type-specific treatment response.

An atlas of late prenatal human neurodevelopment resolved by single-nucleus transcriptomics.

Late prenatal development of the human neocortex encompasses a critical period of gliogenesis and cortical expansion. However, systematic single-cell analyses to resolve cellular diversity and gliogenic lineages of the third trimester are lacking. Here, we present a comprehensive single-nucleus RNA sequencing atlas of over 200,000 nuclei derived from the proliferative germinal matrix and laminating cortical plate of 15 prenatal, non-pathological postmortem samples from 17 to 41 gestational weeks, and 3 adult controls. This dataset captures prenatal gliogenesis with high temporal resolution and is provided as a resource for further interrogation. Our computational analysis resolves greater complexity of glial progenitors, including transient glial intermediate progenitor cell (gIPC) and nascent astrocyte populations in the third trimester of human gestation. We use lineage trajectory and RNA velocity inference to further characterize specific gIPC subpopulations preceding both oligodendrocyte (gIPC-O) and astrocyte (gIPC-A) lineage differentiation. We infer unique transcriptional drivers and biological pathways associated with each developmental state, validate gIPC-A and gIPC-O presence within the human germinal matrix and cortical plate in situ, and demonstrate gIPC states being recapitulated across adult and pediatric glioblastoma tumors.

Hedgehog-interacting protein acts in the habenula to regulate nicotine intake.

Hedgehog-interacting protein (HHIP) sequesters Hedgehog ligands to repress Smoothened (SMO)-mediated recruitment of the GLI family of transcription factors. Allelic variation in HHIP confers risk of chronic obstructive pulmonary disease and other smoking-related lung diseases, but underlying mechanisms are unclear. Using single-cell and cell-type-specific translational profiling, we show that HHIP expression is highly enriched in medial habenula (MHb) neurons, particularly MHb cholinergic neurons that regulate aversive behavioral responses to nicotine. HHIP deficiency dysregulated the expression of genes involved in cholinergic signaling in the MHb and disrupted the function of nicotinic acetylcholine receptors (nAChRs) through a PTCH-1/cholesterol-dependent mechanism. Further, CRISPR/Cas9-mediated genomic cleavage of the Hhip gene in MHb neurons enhanced the motivational properties of nicotine in mice. These findings suggest that HHIP influences vulnerability to smoking-related lung diseases in part by regulating the actions of nicotine on habenular aversion circuits.

Genomic signature of Fanconi anaemia DNA repair pathway deficiency in cancer.

Fanconi anaemia (FA), a model syndrome of genome instability, is caused by a deficiency in DNA interstrand crosslink repair resulting in chromosome breakage1-3. The FA repair pathway protects against endogenous and exogenous carcinogenic aldehydes4-7. Individuals with FA are hundreds to thousands fold more likely to develop head and neck (HNSCC), oesophageal and anogenital squamous cell carcinomas8 (SCCs). Molecular studies of SCCs from individuals with FA (FA SCCs) are limited, and it is unclear how FA SCCs relate to sporadic HNSCCs primarily driven by tobacco and alcohol exposure or infection with human papillomavirus9 (HPV). Here, by sequencing genomes and exomes of FA SCCs, we demonstrate that the primary genomic signature of FA repair deficiency is the presence of high numbers of structural variants. Structural variants are enriched for small deletions, unbalanced translocations and fold-back inversions, and are often connected, thereby forming complex rearrangements. They arise in the context of TP53 loss, but not in the context of HPV infection, and lead to somatic copy-number alterations of HNSCC driver genes. We further show that FA pathway deficiency may lead to epithelial-to-mesenchymal transition and enhanced keratinocyte-intrinsic inflammatory signalling, which would contribute to the aggressive nature of FA SCCs. We propose that the genomic instability in sporadic HPV-negative HNSCC may arise as a result of the FA repair pathway being overwhelmed by DNA interstrand crosslink damage caused by alcohol and tobacco-derived aldehydes, making FA SCC a powerful model to study tumorigenesis resulting from DNA-crosslinking damage.

High-resolution transcriptomics informs glial pathology in human temporal lobe epilepsy.

The pathophysiology of epilepsy underlies a complex network dysfunction between neurons and glia, the molecular cell type-specific contributions of which remain poorly defined in the human disease. In this study, we validated a method that simultaneously isolates neuronal (NEUN +), astrocyte (PAX6 + NEUN-), and oligodendroglial progenitor (OPC) (OLIG2 + NEUN-) enriched nuclei populations from non-diseased, fresh-frozen human neocortex and then applied it to characterize the distinct transcriptomes of such populations isolated from electrode-mapped temporal lobe epilepsy (TLE) surgical samples. Nuclear RNA-seq confirmed cell type specificity and informed both common and distinct pathways associated with TLE in astrocytes, OPCs, and neurons. Compared to postmortem control, the transcriptome of epilepsy astrocytes showed downregulation of mature astrocyte functions and upregulation of development-related genes. To gain further insight into glial heterogeneity in TLE, we performed single cell transcriptomics (scRNA-seq) on four additional human TLE samples. Analysis of the integrated TLE dataset uncovered a prominent subpopulation of glia that express a hybrid signature of both reactive astrocyte and OPC markers, including many cells with a mixed GFAP + OLIG2 + phenotype. A further integrated analysis of this TLE scRNA-seq dataset and a previously published normal human temporal lobe scRNA-seq dataset confirmed the unique presence of hybrid glia only in TLE. Pseudotime analysis revealed cell transition trajectories stemming from this hybrid population towards both OPCs and reactive astrocytes. Immunofluorescence studies in human TLE samples confirmed the rare presence of GFAP + OLIG2 + glia, including some cells with proliferative activity, and functional analysis of cells isolated directly from these samples disclosed abnormal neurosphere formation in vitro. Overall, cell type-specific isolation of glia from surgical epilepsy samples combined with transcriptomic analyses uncovered abnormal glial subpopulations with de-differentiated phenotype, motivating further studies into the dysfunctional role of reactive glia in temporal lobe epilepsy.

NKG2A and HLA-E define an alternative immune checkpoint axis in bladder cancer.

Programmed cell death protein 1 (PD-1)/programmed death-ligand 1 (PD-L1)-blockade immunotherapies have limited efficacy in the treatment of bladder cancer. Here, we show that NKG2A associates with improved survival and responsiveness to PD-L1 blockade immunotherapy in bladder tumors that have high abundance of CD8+ T cells. In bladder tumors, NKG2A is acquired on CD8+ T cells later than PD-1 as well as other well-established immune checkpoints. NKG2A+ PD-1+ CD8+ T cells diverge from classically defined exhausted T cells through their ability to react to human leukocyte antigen (HLA) class I-deficient tumors using T cell receptor (TCR)-independent innate-like mechanisms. HLA-ABC expression by bladder tumors is progressively diminished as disease progresses, framing the importance of targeting TCR-independent anti-tumor functions. Notably, NKG2A+ CD8+ T cells are inhibited when HLA-E is expressed by tumors and partly restored upon NKG2A blockade in an HLA-E-dependent manner. Overall, our study provides a framework for subsequent clinical trials combining NKG2A blockade with other T cell-targeted immunotherapies, where tumors express higher levels of HLA-E.

Dissecting mechanisms of chamber-specific cardiac differentiation and its perturbation following retinoic acid exposure.

The specification of distinct cardiac lineages occurs before chamber formation and acquisition of bona fide atrial or ventricular identity. However, the mechanisms underlying these early specification events remain poorly understood. Here, we performed single cell analysis at the murine cardiac crescent, primitive heart tube and heart tube stages to uncover the transcriptional mechanisms underlying formation of atrial and ventricular cells. We find that progression towards differentiated cardiomyocytes occurs primarily based on heart field progenitor identity, and that progenitors contribute to ventricular or atrial identity through distinct differentiation mechanisms. We identify new candidate markers that define such differentiation processes and examine their expression dynamics using computational lineage trajectory methods. We further show that exposure to exogenous retinoic acid causes defects in ventricular chamber size, dysregulation in FGF signaling and a shunt in differentiation towards orthogonal lineages. Retinoic acid also causes defects in cell-cycle exit resulting in formation of hypomorphic ventricles. Collectively, our data identify, at a single cell level, distinct lineage trajectories during cardiac specification and differentiation, and the precise effects of manipulating cardiac progenitor patterning via retinoic acid signaling.