Mapping Cellular Interactions from Spatially Resolved Transcriptomics Data.

Cell-cell communication (CCC) is essential to how life forms and functions. However, accurate, high-throughput mapping of how expression of all genes in one cell affects expression of all genes in another cell is made possible only recently, through the introduction of spatially resolved transcriptomics technologies (SRTs), especially those that achieve single cell resolution. However, significant challenges remain to analyze such highly complex data properly. Here, we introduce a Bayesian multi-instance learning framework, spacia, to detect CCCs from data generated by SRTs, by uniquely exploiting their spatial modality. We highlight spacia's power to overcome fundamental limitations of popular analytical tools for inference of CCCs, including losing single-cell resolution, limited to ligand-receptor relationships and prior interaction databases, high false positive rates, and most importantly the lack of consideration of the multiple-sender-to-one-receiver paradigm. We evaluated the fitness of spacia for all three commercialized single cell resolution ST technologies: MERSCOPE/Vizgen, CosMx/Nanostring, and Xenium/10X. Spacia unveiled how endothelial cells, fibroblasts and B cells in the tumor microenvironment contribute to Epithelial-Mesenchymal Transition and lineage plasticity in prostate cancer cells. We deployed spacia in a set of pan-cancer datasets and showed that B cells also participate in PDL1/PD1 signaling in tumors. We demonstrated that a CD8+ T cell/PDL1 effectiveness signature derived from spacia analyses is associated with patient survival and response to immune checkpoint inhibitor treatments in 3,354 patients. We revealed differential spatial interaction patterns between γδ T cells and liver hepatocytes in healthy and cancerous contexts. Overall, spacia represents a notable step in advancing quantitative theories of cellular communications.

5-Alpha reductase inhibitors induce a prostate luminal to club cell transition in human benign prostatic hyperplasia.

Benign prostatic hyperplasia (BPH) is a progressive expansion of peri-urethral prostate tissue common in aging men. Patients with enlarged prostates are treated with 5-alpha reductase inhibitors (5ARIs) to shrink prostate volume by blocking the conversion of testosterone to dihydrotestosterone (DHT). A reduction in DHT levels can elicit atrophy and apoptosis of prostate secretory luminal cells, which results in a favorable clinical response characterized by improved lower urinary tract symptoms. However, the histologic response to 5ARI treatment is often heterogeneous across prostate acini and lower urinary tract symptoms can persist to require surgical intervention. We used two spatial profiling approaches to characterize gene expression changes across histologically normal and atrophied regions in prostates from 5ARI-treated men. Objective transcriptomic profiling using the Visium spatial gene expression platform showed that 5ARI-induced atrophy of prostate luminal cells correlated with reduced androgen receptor signaling and increased expression of urethral club cell genes including LTF, PIGR, OLFM4, SCGB1A1, and SCGB3A1. Prostate luminal cells within atrophied acini adapted to decreased DHT conditions by increasing NF-κB signaling and anti-apoptotic BCL2 expression, which may explain their survival. Using GeoMx digital spatial profiling with a probe set to assess ~18 000 RNA targets, we confirmed that atrophied acini expressing SCGB3A1 displayed higher levels of club cell markers compared with histologically normal acini with NKX3-1 expression. In addition, club-like cells within regions of 5ARI-induced atrophy closely resembled true club cells from the prostatic urethra. A comparison of histologically normal regions from 5ARI-treated men and histologically normal regions from untreated men revealed few transcriptional differences. Taken together, our results describe a heterogeneous response to 5ARI treatment where cells in atrophied acini undergo an adaptation from a prostate secretory luminal to a club cell-like state in response to 5ARI treatment. © 2021 The Pathological Society of Great Britain and Ireland.

MP44-09 UNDERSTANDING PRUNE BELLY SYNDROME AT SINGLE CELL RESOLUTION.

INTRODUCTION AND OBJECTIVE: Prune Belly Syndrome (PBS) is characterized by bladder dysmyogenesis, yielding a dysfunctional compliant thick wall with excess collagen deposition. To dissect the cellular heterogeneity and gene expression networks altered in PBS, we report the cell type composition and transcriptional activity of PBS human bladder by using single cell RNA sequencing (scRNA-seq). METHODS: Using IRB-approved methods, bladder dome from 2 PBS and 6 non-PBS control (CO) males underwent fresh single-cell digestion. scRNA-seq was performed and 5277 and 31828 bladder cells from PBS and CO patients was detected, respectively. Cell type clusters were graphically displayed by Uniform Manifold Approximation and Projection (UMap) plot and differentially expressed genes (DEGs) were generated to assign each cluster identity. KEGG (Kyoto Encyclopedia of Genes and Genomes) pathway analysis was performed for PBS affected genes. RESULTS: We identified 17 distinct bladder cell clusters, including 6 fibroblast (1, 2, 3, and 4, immunofibroblast, myofibroblast), 1 smooth muscle (SM), and 2 urothelial (umbrella and basal+intermediate) clusters (Fig 1A-B). Counts of individual cell types were expressed as relative proportions, identifying significant PBS fibroblast enrichment, (67% PBS vs 40% CO). Five of 6 PBS fibroblast sub-types are proportionately fewer in number than in CO. The exception is a dominant fibroblast sub-type we label as fibroblast 4, (61% of all PBS fibroblasts vs <10% CO fibroblast subtypes). SM and urothelial cell populations are dramatically reduced in PBS (SM: 5% PBS vs 11% CO and urothelial: <1% PBS vs 7% CO) (Fig 1C-E). PBS fibroblast DEGs, but not SM cells, are enriched in collagen genes. Fibroblast markers (DCN and PLA2G2A) and SM genes (DES, TPM2 and TAGLN) are reduced (by 4, 13, 2, 4, and 2x respectively) in PBS cells (Fig 1G). KEGG pathways analysis for fibroblasts and SM showed a highly similar enrichment for neurodegenerative disease pathways (Fig 1H-I). CONCLUSIONS: Using scRNA-seq, we identified and characterized the disarrayed cell type populations in PBS bladders, generating their unbiased transcriptomic signatures which highlight commonality with neurodegenerative diseases. This PBS transcriptomic map is a step toward potential markers for diagnosis and therapeutic intervention.[Figure: see text]Source of Funding:NIH DK100483, DK127589 PI: Baker, L.

Single-cell analysis of mouse and human prostate reveals novel fibroblasts with specialized distribution and microenvironment interactions.

Stromal-epithelial interactions are critical to the morphogenesis, differentiation, and homeostasis of the prostate, but the molecular identity and anatomy of discrete stromal cell types is poorly understood. Using single-cell RNA sequencing, we identified and validated the in situ localization of three smooth muscle subtypes (prostate smooth muscle, pericytes, and vascular smooth muscle) and two novel fibroblast subtypes in human prostate. Peri-epithelial fibroblasts (APOD+) wrap around epithelial structures, whereas interstitial fibroblasts (C7+) are interspersed in extracellular matrix. In contrast, the mouse displayed three fibroblast subtypes with distinct proximal-distal and lobe-specific distribution patterns. Statistical analysis of mouse and human fibroblasts showed transcriptional correlation between mouse prostate (C3+) and urethral (Lgr5+) fibroblasts and the human interstitial fibroblast subtype. Both urethral fibroblasts (Lgr5+) and ductal fibroblasts (Wnt2+) in the mouse contribute to a proximal Wnt/Tgfb signaling niche that is absent in human prostate. Instead, human peri-epithelial fibroblasts express secreted WNT inhibitors SFRPs and DKK1, which could serve as a buffer against stromal WNT ligands by creating a localized signaling niche around individual prostate glands. We also identified proximal-distal fibroblast density differences in human prostate that could amplify stromal signaling around proximal prostate ducts. In human benign prostatic hyperplasia, fibroblast subtypes upregulate critical immunoregulatory pathways and show distinct distributions in stromal and glandular phenotypes. A detailed taxonomy of leukocytes in benign prostatic hyperplasia reveals an influx of myeloid dendritic cells, T cells and B cells, resembling a mucosal inflammatory disorder. A receptor-ligand interaction analysis of all cell types revealed a central role for fibroblasts in growth factor, morphogen, and chemokine signaling to endothelia, epithelia, and leukocytes. These data are foundational to the development of new therapeutic targets in benign prostatic hyperplasia. © 2021 The Pathological Society of Great Britain and Ireland. Published by John Wiley & Sons, Ltd.

Identification and characterization of cellular heterogeneity within the developing renal interstitium.

Kidney formation requires the coordinated growth of multiple cell types including the collecting ducts, nephrons, vasculature and interstitium. There is a long-held belief that interactions between progenitors of the collecting ducts and nephrons are primarily responsible for kidney development. However, over the last several years, it has become increasingly clear that multiple aspects of kidney development require signaling from the interstitium. How the interstitium orchestrates these various roles is poorly understood. Here, we show that during development the interstitium is a highly heterogeneous patterned population of cells that occupies distinct positions correlated to the adjacent parenchyma. Our analysis indicates that the heterogeneity is not a mere reflection of different stages in a linear developmental trajectory but instead represents several novel differentiated cell states. Further, we find that ß-catenin has a cell autonomous role in the development of a medullary subset of the interstitium and that this non-autonomously affects the development of the adjacent epithelia. These findings suggest the intriguing possibility that the different interstitial subtypes may create microenvironments that play unique roles in development of the adjacent epithelia and endothelia.

Urethral luminal epithelia are castration-insensitive cells of the proximal prostate.

BACKGROUND: Castration-insensitive epithelial progenitors capable of regenerating the prostate have been proposed to be concentrated in the proximal region based on facultative assays. Functional characterization of prostate epithelial populations isolated with individual cell surface markers has failed to provide a consensus on the anatomical and transcriptional identity of proximal prostate progenitors. METHODS: Here, we use single-cell RNA sequencing to obtain a complete transcriptomic profile of all epithelial cells in the mouse prostate and urethra to objectively identify cellular subtypes. Pan-transcriptomic comparison to human prostate cell types identified a mouse equivalent of human urethral luminal cells, which highly expressed putative prostate progenitor markers. Validation of the urethral luminal cell cluster was performed using immunostaining and flow cytometry. RESULTS: Our data reveal that previously identified facultative progenitors marked by Trop2, Sca-1, KRT4, and PSCA are actually luminal epithelial cells of the urethra that extend into the proximal region of the prostate, and are resistant to castration-induced androgen deprivation. Mouse urethral luminal cells were identified to be the equivalent of previously identified human club and hillock cells that similarly extend into proximal prostate ducts. Benign prostatic hyperplasia (BPH) has long been considered an "embryonic reawakening," but the cellular origin of the hyperplastic growth concentrated in the periurethral region is unclear. We demonstrate an increase in urethral luminal cells within glandular nodules from BPH patients. Urethral luminal cells are further increased in patients treated with a 5-α reductase inhibitor. CONCLUSIONS: Our data demonstrate that cells of the proximal prostate that express putative progenitor markers, and are enriched by castration in the proximal prostate, are urethral luminal cells and that these cells may play an important role in the etiology of human BPH.

A Cellular Anatomy of the Normal Adult Human Prostate and Prostatic Urethra.

A comprehensive cellular anatomy of normal human prostate is essential for solving the cellular origins of benign prostatic hyperplasia and prostate cancer. The tools used to analyze the contribution of individual cell types are not robust. We provide a cellular atlas of the young adult human prostate and prostatic urethra using an iterative process of single-cell RNA sequencing (scRNA-seq) and flow cytometry on ∼98,000 cells taken from different anatomical regions. Immunohistochemistry with newly derived cell type-specific markers revealed the distribution of each epithelial and stromal cell type on whole mounts, revising our understanding of zonal anatomy. Based on discovered cell surface markers, flow cytometry antibody panels were designed to improve the purification of each cell type, with each gate confirmed by scRNA-seq. The molecular classification, anatomical distribution, and purification tools for each cell type in the human prostate create a powerful resource for experimental design in human prostate disease.

Molecular Evolution of Early-Onset Prostate Cancer Identifies Molecular Risk Markers and Clinical Trajectories.

Identifying the earliest somatic changes in prostate cancer can give important insights into tumor evolution and aids in stratifying high- from low-risk disease. We integrated whole genome, transcriptome and methylome analysis of early-onset prostate cancers (diagnosis ≤55 years). Characterization across 292 prostate cancer genomes revealed age-related genomic alterations and a clock-like enzymatic-driven mutational process contributing to the earliest mutations in prostate cancer patients. Our integrative analysis identified four molecular subgroups, including a particularly aggressive subgroup with recurrent duplications associated with increased expression of ESRP1, which we validate in 12,000 tissue microarray tumors. Finally, we combined the patterns of molecular co-occurrence and risk-based subgroup information to deconvolve the molecular and clinical trajectories of prostate cancer from single patient samples.

Molecular pathogenesis of human prostate basal cell hyperplasia.

BACKGROUND: Understanding the molecular pathogenesis of distinct phenotypes in human benign prostatic hyperplasia (BPH) is essential to improving therapeutic intervention. Current therapies target smooth muscle and luminal epithelia for relief of lower urinary tract symptoms (LUTS) due to BPH, but basal cell hyperplasia (BCH) remains untargeted. The incidence of has been reported at 8-10%, but a molecular and cellular characterization has not been performed on this phenotype. METHODS: Using freshly digested tissue from surgical specimens, we performed RNA-seq analysis of flow cytometry-purified basal epithelia from 3 patients with and 4 patients without a majority BCH phenotype. qPCR was performed on 28 genes identified as significant from 13 non-BCH and 7 BCH specimens to confirm transcriptomic analysis. IHC was performed on several non-BCH and BCH specimens for 3 proteins identified as significant by transcriptomic analysis. RESULTS: A total of 141 human BPH specimens were analyzed for the presence of BCH. Clinical characteristics of non-BCH and BCH cohorts revealed no significant differences in age, PSA, prostate volume, medical treatment, or comorbidities. Quantitation of cellular subsets by flow cytometry in 11 BCH patients vs. 11 non-BCH patients demonstrated a significant increase in the ratio of basal to luminal epithelia in patients with BCH (P <0.05), but no significant differences in the total number of leukocytes. RNA-seq data from flow cytometry isolated basal epithelia from patients with and without BCH were subjected to gene set enrichment analysis of differentially expressed genes, which revealed increased expression of members of the epidermal differentiation complex. Transcriptomic data were complemented by immunohistochemistry for members of the epidermal differentiation complex, revealing a morphological similarity to other stratified squamous epithelial layers. CONCLUSIONS: Increased expression of epidermal differentiation complex members and altered epithelial stratification resembles the progression of other metaplastic diseases. These data provide insight into the plasticity of the human prostate epithelium and suggest a classification of basal cell hyperplasia as a metaplasia.

Non-Cell-Autonomous Regulation of Prostate Epithelial Homeostasis by Androgen Receptor.

Prostate inflammation has been suggested as an etiology for benign prostatic hyperplasia (BPH). We show that decreased expression of the androgen receptor (AR) in luminal cells of human BPH specimens correlates with a higher degree of regional prostatic inflammation. However, the cause-and-effect relationship between the two events remains unclear. We investigated specifically whether attenuating AR activity in prostate luminal cells induces inflammation. Disrupting luminal cell AR signaling in mouse models promotes cytokine production cell-autonomously, impairs epithelial barrier function, and induces immune cell infiltration, which further augments local production of cytokines and chemokines including Il-1 and Ccl2. This inflammatory microenvironment promotes AR-independent prostatic epithelial proliferation, which can be abolished by ablating IL-1 signaling or depleting its major cellular source, the macrophages. This study demonstrates that disrupting luminal AR signaling promotes prostate inflammation, which may serve as a mechanism for resistance to androgen-targeted therapy for prostate-related diseases.