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 to 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. This article is protected by copyright. All rights reserved.

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.

Using single-cell RNA sequencing, immunofluorescence, and flow cytometry, Henry et al. create a cellular anatomy of the normal human prostate and provide the tools to identify, isolate, and localize every cell type. They identify two additional epithelial cell types enriched in the prostatic urethra and proximal prostatic ducts.

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.

Accumulation and localisation of maternal RNA in oocytes of Japanese quail

Accumulation of total RNA and poly(A+)RNA was determined in the oocytes of Japanese quail (Coturnix coturnix japonica) during oogenesis, by a standard spectrophotometric method, after RNA extraction. Intensive RNA accumulation was observed in the oocytes 0.25-2.0 mm in diameter (correlated with the presence of lampbrush chromosomes), followed by a plateau in 2.0-3.0 mm oocytes (correlated with the condensation of lampbrush chromosomes). Resumption of the RNA accumulation in the 3.0-5.0 mm oocytes is interpreted as a result of external RNA delivery by the granulosa cells. Most of the RNA (approximately 3/4, 3-4 micrograms) contained in the ovulated oocyte was found under the vitelline membrane surrounding the yolk. This RNA has been designated as 'extraembryonic RNA', as it is located outside the germinal disc region where the embryo is formed. The extraembryonic RNA is rapidly degraded within 24 h, from ovulation until oviposition, during egg passage through the oviduct, while the RNA present in a germinal disc (approximately 1.1 micrograms) is stable. In bird oocytes the presence of two, territorially separated pools of RNA has been postulated: one contained in a germinal disc and needed for early embryo development, and the second present in the cytoplasmic layer around the yolk supporting oocyte growth and development during oogenesis.

Maturation and ovulation of Japanese quail oocytes under in vitro conditions

1. An in vitro system for ovulation and maturation of Japanese quail oocytes is described. 2. Ovarian follicles removed from the ovary at 2, 4 or 6 h before the estimated time of ovulation may ovulate under in vitro conditions. 3. The presence of progesterone in the medium had a stimulatory effect on the process of maturation, as has been shown for Xenopus oocytes.

RNase A activity in Japanese quail oocytes

The presence of RNase A activity was studied in vitro in homogenates of quail oocytes and early embryos using [3H]poly(U) as a substrate. The activity was measured by adsorption of the undegraded substrate onto DE-81 filter paper discs and by chromatographic separation on a Sephadex G-50 column. RNase A activity examined by these methods was almost undetectable in quail previtellogenic, vitellogenic and ovulated oocytes as well as in the embryos from laid eggs. It is estimated to be about 1.1 x 10(-5) Kunitz units per ovulated oocyte. Higher activity starts to appear in gastrulating embryos. These findings are discussed in relation to other studies demonstrating the high stability of maternal RNA during early development, especially in growing oocytes.