Single-Cell Analysis Identifies LY6D as a Marker Linking Castration-Resistant Prostate Luminal Cells to Prostate Progenitors and Cancer

Enhanced prostatic Esr1+ luminal epithelial cells in the absence of SRD5A2

Steroid 5α reductase 2 (SRD5A2) converts testosterone to dihydrotestosterone and is crucial for prostatic development. 5α reductase inhibitors (5ARI) reduce prostate size in benign prostate hyperplasia (BPH) and ameliorate lower urinary tract symptoms secondary to BPH. However, the mechanisms of 5ARI functioning are still not fully understood. Here, we used a Srd5a2-/- mouse model and employed single-cell RNA sequencing to explore the impact of SRD5A2 absence on prostate cellular heterogeneity. Significant alterations in luminal epithelial cell (LE) populations were observed, alongside an increased proportion and proliferative phenotype of estrogen receptor 1 (ESR1)+ LE2 cells, following an SRD5A2-independent ESR1 differentiation trajectory. LE2 cells exhibited enhanced estrogen response gene signatures, suggesting an alternative pathway for prostate growth when SRD5A2 is absent. Human prostate biopsy analysis revealed an inverse correlation between the expressions of SRD5A2 and LE2 markers (ESR1/PKCα), and an inverse correlation between SRD5A2 and the clinical efficiency of 5ARI. These findings provide insights into 5ARI resistance mechanisms and potential alternative therapies for BPH-related lower urinary tract symptoms. © 2024 The Pathological Society of Great Britain and Ireland.

A Single-Cell Atlas of the Murine Pancreatic Ductal Tree Identifies Novel Cell Populations With Potential Implications in Pancreas Regeneration and Exocrine Pathogenesis

BACKGROUND & AIMS

Pancreatic ducts form an intricate network of tubules that secrete bicarbonate and drive acinar secretions into the duodenum. This network is formed by centroacinar cells; terminal, intercalated, intracalated ducts; and the main pancreatic duct. Ductal heterogeneity at the single-cell level has been poorly characterized; therefore, our understanding of the role of ductal cells in pancreas regeneration and exocrine pathogenesis has been hampered by the limited knowledge and unexplained diversity within the ductal network.

METHODS

We used small conditional RNA sequencing to comprehensively characterize mouse ductal heterogeneity at single-cell resolution of the entire ductal epithelium from centroacinar cells to the main duct. Moreover, we used organoid cultures, injury models, and pancreatic tumor samples to interrogate the role of novel ductal populations in pancreas regeneration and exocrine pathogenesis.

RESULTS

We have identified the coexistence of 15 ductal populations within the healthy pancreas and characterized their organoid formation capacity and endocrine differentiation potential. Cluster isolation and subsequent culturing let us identify ductal cell populations with high organoid formation capacity and endocrine and exocrine differentiation potential in vitro, including a Wnt-responsive population, a ciliated population, and FLRT3+ cells. Moreover, we have characterized the location of these novel ductal populations in healthy pancreas, chronic pancreatitis, and tumor samples. The expression of WNT-responsive, interferon-responsive, and epithelial-to-mesenchymal transition population markers increases in chronic pancreatitis and tumor samples.

CONCLUSIONS

In light of our discovery of previously unidentified ductal populations, we unmask potential roles of specific ductal populations in pancreas regeneration and exocrine pathogenesis. Thus, novel lineage-tracing models are needed to investigate ductal-specific populations in vivo.

Integrated single-cell analysis defines the epigenetic basis of castration-resistant prostate luminal cells

Understanding prostate response to castration and androgen receptor signaling inhibitors (ARSI) is critical to improving long-term prostate cancer (PCa) patient survival. Here, we use a multi-omics approach on 229,794 single cells to create a mouse single-cell reference atlas for interpreting mouse prostate biology and castration response. Our reference atlas refines single-cell annotations and provides a chromatin context, which, when coupled with mouse lineage tracing, demonstrates that castration-resistant luminal cells are distinct from the pre-existent urethra-proximal stem/progenitor cells. Molecular pathway analysis and therapeutic studies further implicate AP1 (JUN/FOS), WNT/β-catenin, FOXQ1, NF-κB, and JAK/STAT pathways as major drivers of castration-resistant luminal populations with relevance to human PCa. Our datasets, which can be explored through an interactive portal (https://visportal.roswellpark.org/data/tang/), can aid in developing combination treatments with ARSI for advanced PCa patients.

Obesogenic High-Fat Diet and MYC Cooperate to Promote Lactate Accumulation and Tumor Microenvironment Remodeling in Prostate Cancer

Cancer cells exhibit metabolic plasticity to meet oncogene-driven dependencies while coping with nutrient availability. A better understanding of how systemic metabolism impacts the accumulation of metabolites that reprogram the tumor microenvironment (TME) and drive cancer could facilitate development of precision nutrition approaches. Using the Hi-MYC prostate cancer mouse model, we demonstrated that an obesogenic high-fat diet (HFD) rich in saturated fats accelerates the development of c-MYC-driven invasive prostate cancer through metabolic rewiring. Although c-MYC modulated key metabolic pathways, interaction with an obesogenic HFD was necessary to induce glycolysis and lactate accumulation in tumors. These metabolic changes were associated with augmented infiltration of CD206+ and PD-L1+ tumor-associated macrophages (TAM) and FOXP3+ regulatory T cells, as well as with the activation of transcriptional programs linked to disease progression and therapy resistance. Lactate itself also stimulated neoangiogenesis and prostate cancer cell migration, which were significantly reduced following treatment with the lactate dehydrogenase inhibitor FX11. In patients with prostate cancer, high saturated fat intake and increased body mass index were associated with tumor glycolytic features that promote the infiltration of M2-like TAMs. Finally, upregulation of lactate dehydrogenase, indicative of a lactagenic phenotype, was associated with a shorter time to biochemical recurrence in independent clinical cohorts. This work identifies cooperation between genetic drivers and systemic metabolism to hijack the TME and promote prostate cancer progression through oncometabolite accumulation. This sets the stage for the assessment of lactate as a prognostic biomarker and supports strategies of dietary intervention and direct lactagenesis blockade in treating advanced prostate cancer.

SIGNIFICANCE

Lactate accumulation driven by high-fat diet and MYC reprograms the tumor microenvironment and promotes prostate cancer progression, supporting the potential of lactate as a biomarker and therapeutic target in prostate cancer. See related commentary by Frigo, p. 1742.

Molecular landscape for risk prediction and personalized therapeutics of castration-resistant prostate cancer: at a glance

Prostate cancer (PCa) is commonly occurred with high incidence in men worldwide, and many patients will be eventually suffered from the dilemma of castration-resistance with the time of disease progression. Castration-resistant PCa (CRPC) is an advanced subtype of PCa with heterogeneous carcinogenesis, resulting in poor prognosis and difficulties in therapy. Currently, disorders in androgen receptor (AR)-related signaling are widely acknowledged as the leading cause of CRPC development, and some non-AR-based strategies are also proposed for CRPC clinical analyses. The initiation of CRPC is a consequence of abnormal interaction and regulation among molecules and pathways at multi-biological levels. In this study, CRPC-associated genes, RNAs, proteins, and metabolites were manually collected and integrated by a comprehensive literature review, and they were functionally classified and compared based on the role during CRPC evolution, i.e., drivers, suppressors, and biomarkers, etc. Finally, translational perspectives for data-driven and artificial intelligence-powered CRPC systems biology analysis were discussed to highlight the significance of novel molecule-based approaches for CRPC precision medicine and holistic healthcare.

Integration Analysis of Single-Cell Multi-Omics Reveals Prostate Cancer Heterogeneity

Prostate cancer (PCa) is an extensive heterogeneous disease with a complex cellular ecosystem in the tumor microenvironment (TME). However, the manner in which heterogeneity is shaped by tumors and stromal cells, or vice versa, remains poorly understood. In this study, single-cell RNA sequencing, spatial transcriptomics, and bulk ATAC-sequence are integrated from a series of patients with PCa and healthy controls. A stemness subset of club cells marked with SOX9highARlow expression is identified, which is markedly enriched after neoadjuvant androgen-deprivation therapy (ADT). Furthermore, a subset of CD8+CXCR6+ T cells that function as effector T cells is markedly reduced in patients with malignant PCa. For spatial transcriptome analysis, machine learning and computational intelligence are comprehensively utilized to identify the cellular diversity of prostate cancer cells and cell-cell communication in situ. Macrophage and neutrophil state transitions along the trajectory of cancer progression are also examined. Finally, the immunosuppressive microenvironment in advanced PCa is found to be associated with the infiltration of regulatory T cells (Tregs), potentially induced by an FAP+ fibroblast subset. In summary, the cellular heterogeneity is delineated in the stage-specific PCa microenvironment at single-cell resolution, uncovering their reciprocal crosstalk with disease progression, which can be helpful in promoting PCa diagnosis and therapy.

Hyposecretion of cervical MUC5B is related to preterm birth in pregnant women after cervical excisional surgery

PROBLEM

Excisional surgery for cervical intraepithelial neoplasia is a risk factor for preterm birth in subsequent pregnancies. However, the underlying mechanisms of this association remain unclear. We previously showed that cervical MUC5B, a mucin protein, may be a barrier to ascending pathogens during pregnancy. We thus hypothesized that hyposecretion of cervical MUC5B is associated with preterm birth after cervical excisional surgery.

METHOD OF STUDY

This prospective nested case-control study (Study 1) included pregnant women who had previously undergone cervical excisional surgery across 11 hospitals. We used proteomics to compare cervicovaginal fluid at 18-22 weeks of gestation between the preterm and term birth groups. In another case-control analysis (Study 2), we compared MUC5B expression in nonpregnant uterine tissues between 15 women with a history of cervical excisional surgery and 26 women without a history of cervical surgery.

RESULTS

The abundance of MUC5B in cervicovaginal fluid was significantly decreased in the preterm birth group (fold change = 0.41, p = .035). Among the 480 quantified proteins, MUC5B had the second highest positive correlation with gestational age at delivery in the combined preterm and term groups. The cervicovaginal microbiome composition was not significantly different between the two groups. Cervical length was not correlated with gestational age at delivery (r = 0.18, p = .079). Histologically, the MUC5B-positive area in the nonpregnant cervix was significantly decreased in women with a history of cervical excisional surgery (0.85-fold, p = .048). The distribution of MUC5B-positive areas in the cervical tissues of 26 women without a history of cervical excisional surgery differed across individuals.

CONCLUSIONS

This study suggests that the primary mechanism by which cervical excisional surgery causes preterm birth is the hyposecretion of MUC5B due to loss of the cervical glands.

A single-cell atlas of the murine pancreatic ductal tree identifies novel cell populations with potential implications in pancreas regeneration and exocrine pathogenesis

Background and aims

Pancreatic ducts form an intricate network of tubules that secrete bicarbonate and drive acinar secretions into the duodenum. This network is formed by centroacinar cells, terminal, intercalated, intracalated ducts, and the main pancreatic duct. Ductal heterogeneity at the single-cell level has been poorly characterized; therefore, our understanding of the role of ductal cells in pancreas regeneration and exocrine pathogenesis has been hampered by the limited knowledge and unexplained diversity within the ductal network.

Methods

We used scRNA-seq to comprehensively characterize mouse ductal heterogeneity at single-cell resolution of the entire ductal epithelium from centroacinar cells to the main duct. Moreover, we used organoid cultures, injury models and pancreatic tumor samples to interrogate the role of novel ductal populations in pancreas regeneration and exocrine pathogenesis.

Results

We have identified the coexistence of 15 ductal populations within the healthy pancreas and characterized their organoid formation capacity and endocrine differentiation potential. Cluster isolation and subsequent culturing let us identify ductal cell populations with high organoid formation capacity and endocrine and exocrine differentiation potential in vitro , including Wnt-responsive-population, ciliated-population and FLRT3 + cells. Moreover, we have characterized the location of these novel ductal populations in healthy pancreas, chronic pancreatitis, and tumor samples, highlighting a putative role of WNT-responsive, IFN-responsive and EMT-populations in pancreatic exocrine pathogenesis as their expression increases in chronic pancreatitis and PanIN lesions.

Conclusions

In light of our discovery of previously unidentified ductal populations, we unmask the potential roles of specific ductal populations in pancreas regeneration and exocrine pathogenesis.

A nasal cell atlas reveals heterogeneity of tuft cells and their role in directing olfactory stem cell proliferation

The olfactory neuroepithelium serves as a sensory organ for odors and forms part of the nasal mucosal barrier. Olfactory sensory neurons are surrounded and supported by epithelial cells. Among them, microvillous cells (MVCs) are strategically positioned at the apical surface, but their specific functions are enigmatic, and their relationship to the other specialized epithelial cells is unclear. Here, we establish that the family of MVCs comprises tuft cells and ionocytes in both mice and humans. Integrating analysis of the respiratory and olfactory epithelia, we define the distinct receptor expression of TRPM5+ tuft-MVCs compared with Gɑ-gustducinhigh respiratory tuft cells and characterize a previously undescribed population of glandular DCLK1+ tuft cells. To establish how allergen sensing by tuft-MVCs might direct olfactory mucosal responses, we used an integrated single-cell transcriptional and protein analysis. Inhalation of Alternaria induced mucosal epithelial effector molecules including Chil4 and a distinct pathway leading to proliferation of the quiescent olfactory horizontal basal stem cell (HBC) pool, both triggered in the absence of olfactory apoptosis. Alternaria- and ATP-elicited HBC proliferation was dependent on TRPM5+ tuft-MVCs, identifying these specialized epithelial cells as regulators of olfactory stem cell responses. Together, our data provide high-resolution characterization of nasal tuft cell heterogeneity and identify a function of TRPM5+ tuft-MVCs in directing the olfactory mucosal response to allergens.

Silencing LY6D Expression Inhibits Colon Cancer in Xenograft Mice and Regulates Colon Cancer Stem Cells' Proliferation, Stemness, Invasion, and Apoptosis via the MAPK Pathway

This study explored the role of lymphocyte antigen 6 family member D (LY6D) in colon cancer stem cells' (CCSCs) proliferation and invasion. LY6D was knocked down using siRNA, and the down-regulation of LY6D was verified using Western blotting. After LY6D knockdown, CCSCs' proliferation, stemness, and invasion were suppressed, whereas apoptosis was increased. Gene Ontology (GO) enrichment analysis revealed that the differentially expressed genes (DEGs) between siLY6D and the negative control groups were significantly enriched in the cell-substrate adherens junction, focal adhesion, and cell-substrate junction terms. Meanwhile, the Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis revealed that the DEGs were significantly enriched in the MAPK pathway. In addition, Western blotting results showed that pBRAF and pERK1/2, cascade kinases of the MAPK pathway, were significantly down-regulated after LY6D knockdown. In addition, nude mice xenograft experiments showed that the siLY6D treatment decreased tumor sizes and weights and improved tumor-bearing mice survival rates compared with the control group. In conclusion, these findings indicate that LY6D, which is highly expressed in CCSCs, is a key factor involved in tumor growth and development and might be a potential cancer marker and therapeutic target for colon cancer.

Analysis of genetic biomarkers, polymorphisms in ADME-related genes and their impact on pharmacotherapy for prostate cancer

Prostate cancer (PCa) is a non-cutaneous malignancy in males with wide variation in incidence rates across the globe. It is the second most reported cause of cancer death. Its etiology may have been linked to genetic polymorphisms, which are not only dominating cause of malignancy casualties but also exerts significant effects on pharmacotherapy outcomes. Although many therapeutic options are available, but suitable candidates identified by useful biomarkers can exhibit maximum therapeutic efficacy. The single-nucleotide polymorphisms (SNPs) reported in androgen receptor signaling genes influence the effectiveness of androgen receptor pathway inhibitors and androgen deprivation therapy. Furthermore, SNPs located in genes involved in transport, drug metabolism, and efflux pumps also influence the efficacy of pharmacotherapy. Hence, SNPs biomarkers provide the basis for individualized pharmacotherapy. The pharmacotherapeutic options for PCa include hormonal therapy, chemotherapy (Docetaxel, Mitoxantrone, Cabazitaxel, and Estramustine, etc.), and radiotherapy. Here, we overview the impact of SNPs reported in various genes on the pharmacotherapy for PCa and evaluate current genetic biomarkers with an emphasis on early diagnosis and individualized treatment strategy in PCa.

LY6D is crucial for lipid accumulation and inflammation in nonalcoholic fatty liver disease

Nonalcoholic fatty liver disease (NAFLD) is a serious metabolic disorder characterized by excess fat accumulation in the liver. Over the past decade, NAFLD prevalence and incidence have risen globally. There are currently no effective licensed drugs for its treatment. Thus, further study is required to identify new targets for NAFLD prevention and treatment. In this study, we fed C57BL6/J mice one of three diets, a standard chow diet, high-sucrose diet, or high-fat diet, and then characterized them. The mice fed a high-sucrose diet had more severely compacted macrovesicular and microvesicular lipid droplets than those in the other groups. Mouse liver transcriptome analysis identified lymphocyte antigen 6 family member D (Ly6d) as a key regulator of hepatic steatosis and the inflammatory response. Data from the Genotype-Tissue Expression project database showed that individuals with high liver Ly6d expression had more severe NAFLD histology than those with low liver Ly6d expression. In AML12 mouse hepatocytes, Ly6d overexpression increased lipid accumulation, while Ly6d knockdown decreased lipid accumulation. Inhibition of Ly6d ameliorated hepatic steatosis in a diet-induced NAFLD mouse model. Western blot analysis showed that Ly6d phosphorylated and activated ATP citrate lyase, which is a key enzyme in de novo lipogenesis. In addition, RNA- and ATAC-sequencing analyses revealed that Ly6d drives NAFLD progression by causing genetic and epigenetic changes. In conclusion, Ly6d is responsible for the regulation of lipid metabolism, and inhibiting Ly6d can prevent diet-induced steatosis in the liver. These findings highlight Ly6d as a novel therapeutic target for NAFLD.

Autocrine amphiregulin signaling sustains castration-resistant Ly6d+ prostate cancer cells

Mechanisms underlying the development of castration-resistant prostate cancer (CRPC) remain incompletely understood. A recent study by Steiner et al. shows that Ly6d⁺ prostate tumor cells survive androgen deprivation through an autocrine amphiregulin signaling pathway.

Identification of a novel glycolysis-related prognosis risk signature in triple-negative breast cancer

Introduction

Triple-negative breast cancer (TNBC) is a particularly aggressive cluster of breast cancer characterized by significant molecular heterogeneity. Glycolysis is a metabolic pathway that is significantly associated with cancer progression, metastasis, recurrence and chemoresistance. However, the potential roles of glycolysis-related genes in TNBC remain unclear.

Methods

In the present study, we identified 108 glycolysis-related differentially expressed genes (DEGs) between breast cancer (BRCA) tumor tissues and normal tissues, and we divided patients into two different clusters with significantly distinct molecular characteristics, clinicopathological features, prognosis, immune cell infiltration and mutation burden. We then constructed a 10-gene signature that classified all TNBCs into low- and high-risk groups.

Results

The high-risk group had significantly lower survival than the low-risk group, which implied that the risk score was an independent prognostic indicator for TNBC patients. Consequently, we constructed and validated a prognostic nomogram, which accurately predicted individual overall survival (OS) of TNBC. Moreover, the risk score predicted the drug sensitivity of chemotherapeutic agents and immunotherapy for TNBC patients.

Discussion

The present comprehensive analysis of glycolysis-related DEGs in TNBC provides new methods for prognosis prediction and more effective treatment strategies.

Prostate organoids: emerging experimental tools for translational research

Organoid technology has provided new translational research opportunities in oncology, in part by enabling the development of patient-representative living biobanks. Prostate cancer research historically has been constrained to a small number of in vitro models, limiting the ability to translate experimental conclusions for contemporary, heterogeneous patient populations. The facility of organoid culture methods to maintain luminal prostate epithelia, the common lineage of prostate cancers, has greatly expanded the phenotypic and genotypic diversity of available tractable models, including luminal stem/progenitor cells and progressive patient-derived cancers. Biobanks of patient prostate cancer organoids enable increased accuracy in predicting therapeutic efficacy and informative clinical trial designs. Here, we discuss how prostate organoid technology is currently being used, the promising areas of future therapeutic applications, and the current obstacles to be overcome.

Autocrine activation of MAPK signaling mediates intrinsic tolerance to androgen deprivation in LY6D prostate cancer cells

The emergence of castration-resistant prostate cancer remains an area of unmet clinical need. We recently identified a subpopulation of normal prostate progenitor cells, characterized by an intrinsic resistance to androgen deprivation and expression of LY6D. We here demonstrate that conditional deletion of PTEN in the murine prostate epithelium causes an expansion of transformed LY6D+ progenitor cells without impairing stem cell properties. Transcriptomic analyses of LY6D+ luminal cells identified an autocrine positive feedback loop, based on the secretion of amphiregulin (AREG)-mediated activation of mitogen-activated protein kinase (MAPK) signaling, increasing cellular fitness and organoid formation. Pharmacological interference with this pathway overcomes the castration-resistant properties of LY6D+ cells with a suppression of organoid formation and loss of LY6D+ cells in vivo. Notably, LY6D+ tumor cells are enriched in high-grade and androgen-resistant prostate cancer, providing clinical evidence for their contribution to advanced disease. Our data indicate that early interference with MAPK inhibitors can prevent progression of castration-resistant prostate cancer.

A luminal intermediate cell state maintains long-term prostate homeostasis and contributes to tumorigenesis

Cellular heterogeneity poses tremendous challenges for developing cell-targeted therapies and biomarkers of clinically significant prostate cancer. The origins of this heterogeneity within normal adult and aging tissue remain unknown, leaving cellular states and transcriptional programs that allow expansions of malignant clones unidentified. To define cell states that contribute to early cancer development, we performed clonal analyses and single cell transcriptomics of normal prostate from genetically-engineered mouse models. We uncovered a luminal transcriptional state with a unique "basal-like" Wnt/p63 signaling ( luminal intermediate , LumI) which contributes to the maintenance of long-term prostate homeostasis. Moreover, LumI cells greatly expand during early stages of tumorigenesis in several mouse models of prostate cancer. Genetic ablation of p63 in vivo in luminal cells reduced the formation of aggressive clones in mouse prostate tumor models. Finally, the LumI cells and Wnt signaling appear to significantly increase in human aging prostate and prostate cancer samples, highlighting the importance of this hybrid cell state for human pathologies with potential translational impact.

Innate Immune Program in Formation of Tumor-Initiating Cells from Cells-of-Origin of Breast, Prostate, and Ovarian Cancers

Tumor-initiating cells (TICs), also known as cancer stem cells (CSCs), are cancer cells that can initiate a tumor, possess self-renewal capacity, and can contribute to tumor heterogeneity. TICs/CSCs are developed from their cells-of-origin. In breast, prostate, and ovarian cancers, progenitor cells for mammary alveolar cells, prostate luminal (secretory) cells, and fallopian tube secretory cells are the preferred cellular origins for their corresponding cancer types. These luminal progenitors (LPs) express common innate immune program (e.g., Toll-like receptor (TLR) signaling)-related genes. Microbes such as bacteria are now found in breast, prostate, and fallopian tube tissues and their corresponding cancer types, raising the possibility that their LPs may sense the presence of microbes and trigger their innate immune/TLR pathways, leading to an inflammatory microenvironment. Crosstalk between immune cells (e.g., macrophages) and affected epithelial cells (e.g., LPs) may eventually contribute to formation of TICs/CSCs from their corresponding LPs, in part via STAT3 and/or NFκB pathways. As such, TICs/CSCs can inherit expression of innate-immunity/TLR-pathway-related genes from their cells-of-origin; the innate immune program may also represent their unique vulnerability, which can be explored therapeutically (e.g., by enhancing immunotherapy via augmenting TLR signaling).

Prognostic biomarkers of pancreatic cancer identified based on a competing endogenous RNA regulatory network

BACKGROUND

Pancreatic cancer is an insidious and heterogeneous malignancy with poor prognosis that is often locally unresectable. Therefore, determining the underlying mechanisms and effective prognostic indicators of pancreatic cancer may help optimize clinical management. This study was conducted to develop a prognostic model for pancreatic cancer based on a competing endogenous RNA (ceRNA) network.

METHODS

We obtained transcriptomic data and corresponding clinicopathological information of pancreatic cancer samples from The Cancer Genome Atlas (TCGA) database (training set). Based on the ceRNA interaction network, we screened candidate genes to build prediction models. Univariate Cox regression analysis was performed to screen for genes associated with prognosis, and least absolute shrinkage and selection operator (LASSO) regression analysis was conducted to construct a predictive model. A receiver operating characteristic (ROC) curve was drawn, and the C-index was calculated to evaluate the accuracy of the prediction model. Furthermore, we downloaded transcriptomic data and related clinical information of pancreatic cancer samples from the Gene Expression Omnibus database (validation set) to evaluate the robustness of our prediction model.

RESULTS

Eight genes (ANLN, FHDC1, LY6D, SMAD6, ACKR4, RAB27B, AUNIP, and GPRIN3) were used to construct the prediction model, which was confirmed as an independent predictor for evaluating the prognosis of patients with pancreatic cancer through univariate and multivariate Cox regression analysis. By plotting the decision curve, we found that the risk score model is an independent predictor has the greatest impact on survival compared to pathological stage and targeted molecular therapy.

CONCLUSIONS

An eight-gene prediction model was constructed for effectively and independently predicting the prognosis of patients with pancreatic cancer. These eight genes identified show potential as diagnostic and therapeutic targets.

RNA-seq and ChIP-seq Identification of Unique and Overlapping Targets of GLI Transcription Factors in Melanoma Cell Lines

BACKGROUND

Despite significant progress in therapy, melanoma still has a rising incidence worldwide, and novel treatment strategies are needed. Recently, researchers have recognized the involvement of the Hedgehog-GLI (HH-GLI) signaling pathway in melanoma and its consistent crosstalk with the MAPK pathway. In order to further investigate the link between the two pathways and to find new target genes that could be considered for combination therapy, we set out to find transcriptional targets of all three GLI proteins in melanoma.

METHODS

We performed RNA sequencing on three melanoma cell lines (CHL-1, A375, and MEL224) with overexpressed GLI1, GLI2, and GLI3 and combined them with the results of ChIP-sequencing on endogenous GLI1, GLI2, and GLI3 proteins. After combining these results, 21 targets were selected for validation by qPCR.

RESULTS

RNA-seq revealed a total of 808 differentially expressed genes (DEGs) for GLI1, 941 DEGs for GLI2, and 58 DEGs for GLI3. ChIP-seq identified 527 genes that contained GLI1 binding sites in their promoters, 1103 for GLI2 and 553 for GLI3. A total of 15 of these targets were validated in the tested cell lines, 6 of which were detected by both RNA-seq and ChIP-seq.

CONCLUSIONS

Our study provides insight into the unique and overlapping transcriptional output of the GLI proteins in melanoma. We suggest that our findings could provide new potential targets to consider while designing melanoma-targeted therapy.

Lineage Tracing and Molecular Real-Time Imaging of Cancer Stem Cells

The cancer stem cells (CSC) are the roots of cancer. The CSC hypothesis may provide a model to explain the tumor cell heterogeneity. Understand the biological mechanism of CSC will help the early detection and cure of cancer. The discovery of the dynamic changes in CSC will be possible by the using of bio-engineering techniques-lineage tracing. However, it is difficult to obtain real-time, continuous, and dynamic live-imaging information using the traditional approaches that take snapshots of time points from different animals. The goal of molecular imaging is to monitor the in situ, continuous molecular changes of cells in vivo. Therefore, the most advanced bioengineering lineage tracing approach, while using a variety of molecular detection methods, will maximize the presentation of CSC. In this review, we first introduce the method of lineage tracing, and then introduce the various components of molecular images to dynamic detect the CSC. Finally, we analyze the current situation and look forward the future of CSC detection.

Expansion of mouse castration-resistant intermediate prostate stem cells in vitro

Background

Most castration-resistant prostate cancers (CRPCs) have a luminal phenotype with high androgen receptor (AR) and prostate-specific antigen (PSA) expression. Currently, it is difficult to culture castration-resistant luminal cells with AR and PSA expression.

Methods

We formulated a custom-made medium and isolated primary cells from the prostate of adult wild-type (WT) and TRAMP mice. The cells were characterized by immunofluorescence staining, transcriptomic analysis, and qRT-PCR verification. Their self-renewal and differentiation potential in vitro and in vivo were examined. We treated the cells with androgen deprivation and enzalutamide and performed immunofluorescence staining and western blotting to analyze their expression of AR and PSA.

Results

We isolated a novel type of castration-resistant intermediate prostate stem cells (CRIPSCs) from adult WT and TRAMP mice. The mouse CRIPSCs proliferated rapidly in two-dimensional (2D) culture dishes and can be cultured for more than six months. The mouse CRIPSCs expressed luminal markers (AR, PSA, and Dsg4), basal markers (CK5 and p63), Psca, and the intermediate cell marker (Ivl). Transcriptomic analysis showed that the mouse CRIPSCs had upregulated signaling pathways related to cancer development and drug resistance. In the long-term culture, TRAMP CRIPSCs had higher expression of the genes related to stem cells and cancers than WT mice. Both WT and TRAMP CRIPSCs formed organoids in Matrigel. WT CRIPSCs did not form prostate tissues when transplanted in vivo without urogenital sinus mesenchyme (UGM) cells. In contrast, TRAMP CRIPSCs formed prostate ducts in NOG mice without UGM  cells and differentiated into luminal, basal, and neuroendocrine cells. Androgens regulated AR translocation between the nucleus and cytoplasm in the mouse CRIPSCs. Treatment of androgen deprivation  (ADT) and enzalutamide reduced AR expression in WT and TRAMP CRIPSCs; however, this treatment promoted PSA expression in TRAMP, while not WT CRIPSCs, similar to the clinical observations of CRPC.

Conclusions

Our study established a method for isolating and expanding mouse CRIPSCs in 2D culture dishes. Mouse CRIPSCs had markers of basal and luminal cells, including AR and PSA, and can differentiate into prostate organoids and tissues. TRAMP CRIPSCs had elevated PSA expression upon ADT and enzalutamide treatment. Our method can be translated into clinical settings for CRPC precision medicine.

MUC1-C Dictates JUN and BAF-Mediated Chromatin Remodeling at Enhancer Signatures in Cancer Stem Cells

The oncogenic MUC1-C protein promotes dedifferentiation of castrate-resistant prostate cancer (CRPC) and triple-negative breast cancer (TNBC) cells. Chromatin remodeling is critical for the cancer stem cell (CSC) state; however, there is no definitive evidence that MUC1-C regulates chromatin accessibility and thereby expression of stemness-associated genes. We demonstrate that MUC1-C drives global changes in chromatin architecture in the dedifferentiation of CRPC and TNBC cells. Our results show that MUC1-C induces differentially accessible regions (DAR) across their genomes, which are significantly associated with differentially expressed genes (DEG). Motif and cistrome analysis further demonstrated MUC1-C-induced DARs align with genes regulated by the JUN/AP-1 family of transcription factors. MUC1-C activates the BAF chromatin remodeling complex, which is recruited by JUN in enhancer selection. In studies of the NOTCH1 gene, which is required for CRPC and TNBC cell self-renewal, we demonstrate that MUC1-C is necessary for (i) occupancy of JUN and ARID1A/BAF, (ii) increases in H3K27ac and H3K4me3 signals, and (iii) opening of chromatin accessibility on a proximal enhancer-like signature. Studies of the EGR1 and LY6E stemness-associated genes further demonstrate that MUC1-C-induced JUN/ARID1A complexes regulate chromatin accessibility on proximal and distal enhancer-like signatures. These findings uncover a role for MUC1-C in chromatin remodeling that is mediated at least in part by JUN/AP-1 and ARID1A/BAF in association with driving the CSC state.

IMPLICATIONS

These findings show that MUC1-C, which is necessary for the CRPC and TNBC CSC state, activates a novel pathway involving JUN/AP-1 and ARID1A/BAF that regulates chromatin accessibility of stemness-associated gene enhancers.

Intra-epithelial non-canonical Activin A signaling safeguards prostate progenitor quiescence

The healthy prostate is a relatively quiescent tissue. Yet, prostate epithelium overgrowth is a common condition during aging, associated with urinary dysfunction and tumorigenesis. For over thirty years, TGF-β ligands have been known to induce cytostasis in a variety of epithelia, but the intracellular pathway mediating this signal in the prostate, and its relevance for quiescence, have remained elusive. Here, using mouse prostate organoids to model epithelial progenitors, we find that intra-epithelial non-canonical Activin A signaling inhibits cell proliferation in a Smad-independent manner. Mechanistically, Activin A triggers Tak1 and p38 ΜAPK activity, leading to p16 and p21 nuclear import. Spontaneous evasion from this quiescent state occurs upon prolonged culture, due to reduced Activin A secretion, a condition associated with DNA replication stress and aneuploidy. Organoids capable to escape quiescence in vitro are also able to implant with increased frequency into immunocompetent mice. This study demonstrates that non-canonical Activin A signaling safeguards epithelial quiescence in the healthy prostate, with potential implications for the understanding of cancer initiation, and the development of therapies targeting quiescent tumor progenitors.

Identification of DDX31 as a Potential Oncogene of Invasive Metastasis and Proliferation in PDAC

Background: Pancreatic ductal adenocarcinoma (PDAC) is one of the deadliest malignant tumors worldwide and has poor prognosis. DEAD box proteins31 (DDX31) participate in cellular processes involving RNA secondary structure changes. However, the functions of DDX31 in PDAC remain to be elucidated.

Methods: The key gene DDX31 was identified using a combination of a risk model and weighted gene co-expression network analysis (WGCNA) with R software. The biological functions of DDX31 in PDAC were investigated through bioinformatics analysis and in vitro experiments.

Results: Combining with WGCNA and risk model, DDX31 was identified as a potential factor of the invasive metastasis properties of PDAC, and its expression was closely related to the malignant differentiation of PDAC. The results of gene set enrichment analysis (GSEA) showed that DDX31 was correlated with cell invasive metastasis and proliferation by activating MAPK signaling pathway. The inhibition of DDX31 inhibited the invasion and migration of PDAC cells. Survival analysis showed that DDX31 expression was negatively associated with the poor prognosis in patients with PDAC.

Interpretation:DDX31 may be a potential factor for PDAC. The inhibition of DDX31 may be a potential way to treat PDAC.

Prostate zones and cancer: lost in transition?

Localized prostate cancer shows great clinical, genetic and environmental heterogeneity; however, prostate cancer treatment is currently guided solely by clinical staging, serum PSA levels and histology. Increasingly, the roles of differential genomics, multifocality and spatial distribution in tumorigenesis are being considered to further personalize treatment. The human prostate is divided into three zones based on its histological features: the peripheral zone (PZ), the transition zone (TZ) and the central zone (CZ). Each zone has variable prostate cancer incidence, prognosis and outcomes, with TZ prostate tumours having better clinical outcomes than PZ and CZ tumours. Molecular and cell biological studies can improve understanding of the unique molecular, genomic and zonal cell type features that underlie the differences in tumour progression and aggression between the zones. The unique biology of each zonal tumour type could help to guide individualized treatment and patient risk stratification.

Heterogeneity and complexity of the prostate epithelium: New findings from single-cell RNA sequencing studies

The recent advent of single-cell RNA-sequencing technology has provided new fundamental insights into the heterogeneity of the prostate epithelium. Several independent studies have described extensive heterogeneity of the luminal epithelial compartment, including a major division between a novel population of luminal cells located in the proximal region of the prostate ducts versus luminal cells located more distally. Proximal luminal cells as well as novel periurethral cells display increased progenitor potential in organoid culture and tissue reconstitution assays, but not in lineage-tracing analyses during prostate homeostasis, suggesting context-dependent plasticity of these populations. Here we describe and synthesize recent findings regarding the epithelial cell populations in the mouse prostate, draw comparisons to the human prostate, and address the relevance of these findings to prostate diseases and cancer.

Single-cell sequencing and its applications in bladder cancer

Bladder cancer is the most common malignant tumour of the urinary system that is characterised by significant intra-tumoural heterogeneity. While large-scale sequencing projects have provided a preliminary understanding of tumour heterogeneity, these findings are based on the average signals obtained from the pooled populations of diverse cells. Recent advances in single-cell sequencing (SCS) technologies have been critical in this regard, opening up new ways of understanding the nuanced tumour biology by identifying distinct cellular subpopulations, dissecting the tumour microenvironment, and characterizing cellular genomic mutations. By integrating these novel insights, SCS technologies are expected to make powerful and meaningful changes to the current diagnosis and treatment of bladder cancer through the identification and usage of novel biomarkers as well as targeted therapeutics. SCS can discriminate complex heterogeneity in a large population of tumour cells and determine the key molecular properties that influence clinical outcomes. Here, we review the advances in single-cell technologies and discuss their applications in cancer research and clinical practice, with a specific focus on bladder cancer.

Experimental challenges to modeling prostate cancer heterogeneity

Tumor heterogeneity plays a key role in prostate cancer prognosis, therapy selection, relapse, and acquisition of treatment resistance. Prostate cancer presents a heterogeneous diversity at inter- and intra-tumor and inter-patient levels which are influenced by multiple intrinsic and/or extrinsic factors. Recent studies have started to characterize the complexity of prostate tumors and these different tiers of heterogeneity. In this review, we discuss the most common factors that contribute to tumoral diversity. Moreover, we focus on the description of the in vitro and in vivo approaches, as well as high-throughput technologies, that help to model intra-tumoral diversity. Further understanding tumor heterogeneities and the challenges they present will guide enhanced patient risk stratification, aid the design of more precise therapies, and ultimately help beat this chameleon-like disease.

Systematic Analysis of an Invasion-Related 3-Gene Signature and Its Validation as a Prognostic Model for Pancreatic Cancer

Background

Pancreatic adenocarcinoma (PAAD) is a malignant tumor of the digestive system that is associated with a poor prognosis in patients owing to its rapid progression and high invasiveness.

Methods

Ninety-seven invasive-related genes obtained from the CancerSEA database were clustered to obtain the molecular subtype of pancreatic cancer based on the RNA-sequencing (RNA-seq) data of The Cancer Genome Atlas (TCGA). The differentially expressed genes (DEGs) between subtypes were obtained using the limma package in R, and the multi-gene risk model based on DEGs was constructed by Lasso regression analysis. Independent datasets GSE57495 and GSE62452 were used to validate the prognostic value of the risk model. To further explore the expression of the hub genes, immunohistochemistry was performed on PAAD tissues obtained from a large cohort.

Results

The TCGA-PAAD samples were divided into two subtypes based on the expression of the invasion-related genes: C1 and C2. Most genes were overexpressed in the C1 subtype. The C1 subtype was mainly enriched in tumor-related signaling pathways, and the prognosis of patients with the C1 subtype was significantly worse than those with the C2 subtype. A 3-gene signature consisting of LY6D, BCAT1, and ITGB6 based on 538 DEGs between both subtypes serves as a stable prognostic marker in patients with pancreatic cancer across multiple cohorts. LY6D, BCAT1, and ITGB6 were over-expressed in 120 PAAD samples compared to normal samples.

Conclusions

The constructed 3-gene signature can be used as a molecular marker to assess the prognostic risk in patients with PAAD.

Recent advances in tissue stem cells

Stem cells are undifferentiated cells capable of self-renewal and differentiation, giving rise to specialized functional cells. Stem cells are of pivotal importance for organ and tissue development, homeostasis, and injury and disease repair. Tissue-specific stem cells are a rare population residing in specific tissues and present powerful potential for regeneration when required. They are usually named based on the resident tissue, such as hematopoietic stem cells and germline stem cells. This review discusses the recent advances in stem cells of various tissues, including neural stem cells, muscle stem cells, liver progenitors, pancreatic islet stem/progenitor cells, intestinal stem cells, and prostate stem cells, and the future perspectives for tissue stem cell research.

Construction of ceRNA network to identify the lncRNA and mRNA related to non-small cell lung cancer

BACKGROUND

Non-small cell lung cancer (NSCLC) harms human health, but its pathogenesis remains unclear. We wish to provide more molecular therapeutic targets for NSCLC.

METHODS

The NSCLC tissue and normal tissue samples were screened for genetic comparison in the TCGA database. The predicted lncRNA and mRNA in BEAS2B and A549 cells were detected.

RESULTS

Volcano plot displayed differentially expressed lncRNAs and mRNAs in adjacent tissues and NSCLC tissues. The survival curve showed that the lncRNA and mRNA had a significant impact on the patient's survival. The results of GO term enrichment analysis indicated that mRNA functions were enriched in cell cycle-related pathways. In the ceRNA interaction network, 13 lncRNAs and 20 miRNAs were found to have an interactive relationship. Finally, 3 significantly different lncRNAs (LINC00968, lnc-FAM92A-9 and lnc-PTGFR-1) and 6 mRNAs (CTCFL, KRT5, LY6D, TMEM, GBP6, and TMEM179) with potential therapeutic significance were screened out. And the cell experiment verified our results.

CONCLUSION

We screened out clinically significant 3 lncRNAs and 6 mRNAs involved in the ceRNA network, which were the key to our future research on the treatment of NSCLC.

Organoid technology for tissue engineering

For centuries, attempts have been continuously made to artificially reconstitute counterparts of in vivo organs from their tissues or cells. Only in the recent decade has organoid technology as a whole technological field systematically emerged and been shown to play important roles in tissue engineering. Based on their self-organizing capacities, stem cells of versatile organs, both harvested and induced, can form 3D structures that are structurally and functionally similar to their in vivo counterparts. These organoid models provide a powerful platform for elucidating the development mechanisms, modeling diseases, and screening drug candidates. In this review, we will summarize the advances of this technology for generating various organoids of tissues from the three germ layers and discuss their drawbacks and prospects for tissue engineering.

Inter- and intra-tumor heterogeneity of metastatic prostate cancer determined by digital spatial gene expression profiling

Metastatic prostate cancer (mPC) comprises a spectrum of diverse phenotypes. However, the extent of inter- and intra-tumor heterogeneity is not established. Here we use digital spatial profiling (DSP) technology to quantitate transcript and protein abundance in spatially-distinct regions of mPCs. By assessing multiple discrete areas across multiple metastases, we find a high level of intra-patient homogeneity with respect to tumor phenotype. However, there are notable exceptions including tumors comprised of regions with high and low androgen receptor (AR) and neuroendocrine activity. While the vast majority of metastases examined are devoid of significant inflammatory infiltrates and lack PD1, PD-L1 and CTLA4, the B7-H3/CD276 immune checkpoint protein is highly expressed, particularly in mPCs with high AR activity. Our results demonstrate the utility of DSP for accurately classifying tumor phenotype, assessing tumor heterogeneity, and identifying aspects of tumor biology involving the immunological composition of metastases.

The inter- and intra-tumor heterogeneity of metastatic prostate cancer (mPC) is underexplored. Here the authors use Digital Spatial Profiling to study gene and protein expression heterogeneity in 27 mPC patients, finding variation in associated pathways and potential immunotherapy targets.

Identification of UBE2C as hub gene in driving prostate cancer by integrated bioinformatics analysis

BACKGROUND

The aim of this study was to identify novel genes in promoting primary prostate cancer (PCa) progression and to explore its role in the prognosis of prostate cancer.

METHODS

Four microarray datasets containing primary prostate cancer samples and benign prostate samples were downloaded from Gene Expression Omnibus (GEO), then differentially expressed genes (DEGs) were identified by R software (version 3.6.2). Gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) were performed to identify the function of DEGs. Using STRING and Cytoscape (version 3.7.1), we constructed a protein-protein interaction (PPI) network and identified the hub gene of prostate cancer. Clinical data on GSE70770 and TCGA was collected to show the role of hub gene in prostate cancer progression. The correlations between hub gene and clinical parameters were also indicated by cox regression analysis. Gene Set Enrichment Analysis (GSEA) was performed to highlight the function of Ubiquitin-conjugating enzyme complex (UBE2C) in prostate cancer.

RESULTS

243 upregulated genes and 298 downregulated genes that changed in at least two microarrays have been identified. GO and KEGG analysis indicated significant changes in the oxidation-reduction process, angiogenesis, TGF-beta signaling pathway. UBE2C, PDZ-binding kinase (PBK), cyclin B1 (CCNB1), Cyclin-dependent kinase inhibitor 3 (CDKN3), topoisomerase II alpha (TOP2A), Aurora kinase A (AURKA) and MKI67 were identified as the candidate hub genes, which were all correlated with prostate cancer patient' disease-free survival in TCGA. In fact, only UBE2C was highly expressed in prostate cancer when compared with benign prostate tissue in TCGA and the expression of UBE2C was also in parallel with the Gleason score of prostate cancer. Cox regression analysis has indicated UBE2C could function as the independent prognostic factor of prostate cancer. GSEA showed UBE2C had played an important role in the pathway of prostate cancer, such as NOTCH signaling pathway, WNT-β-catenin signaling pathway.

CONCLUSIONS

UBE2C was pivotal for the progression of prostate cancer and the level of UBE2C was important to predict the prognosis of patients.

High mRNA expression of LY6 gene family is associated with overall survival outcome in pancreatic ductal adenocarcinoma

Pancreatic cancer ranks one of the worst in overall survival outcome with a 5 year survival rate being less than 10%. Pancreatic cancer faces unique challenges in its diagnosis and treatment, such as the lack of clinically validated biomarkers and the immensely immunosuppressive tumor microenvironment. Recently, the LY6 gene family has received increasing attention for its multi-faceted roles in cancer development, stem cell maintenance, immunomodulation, and association with more aggressive and hard-to-treat cancers. A detailed study of mRNA expression of LY6 gene family and its association with overall survival (OS) outcome in pancreatic cancers is lacking. We used publicly available clinical datasets to analyze the mRNA expression of a set of LY6 genes and its effect on OS outcome in the context of the tumor microenvironment and immunomodulation. We used web-based tools Kaplan-Meier Plotter, cBioPortal, Oncomine and R-programming to analyze copy number alterations, mRNA expression and its association with OS outcome in pancreatic cancer. These analyses demonstrated that high expression of LY6 genes is associated with OS and disease free survival (DFS) outcome. High expression of LY6 genes and their association with OS outcome is dependent on the composition of tumor microenvironment. Considering that LY6 proteins are anchored to the outer cell membrane or secreted, making them readily accessible, these findings highlight the potential of LY6 family members in the future of pancreatic cancer diagnosis and treatment.

Progenitors in prostate development and disease

The prostate develops by epithelial budding and branching processes that occur during fetal and postnatal stages. The adult prostate demonstrates remarkable regenerative capacity, with the ability to regrow to its original size over multiple cycles of castration and androgen administration. This capacity for controlled regeneration prompted the search for an androgen-independent epithelial progenitor in benign prostatic hyperplasia (BPH) and prostate cancer (PCa). BPH is hypothesized to be a reawakening of ductal branching, resulting in the formation of new proximal glands, all while androgen levels are decreasing in the aging male. Advanced prostate cancer can be slowed with androgen deprivation, but resistance eventually occurs, suggesting the existence of an androgen-independent progenitor. Recent studies indicate that there are multiple castration-insensitive epithelial cell types in the proximal area of the prostate, but not all act as progenitors during prostate development or regeneration. This review highlights how recent cellular and anatomical studies are changing our perspective on the identity of the prostate progenitor.

Vitamin D sufficiency enhances differentiation of patient-derived prostate epithelial organoids

Vitamin D is an essential steroid hormone that regulates systemic calcium homeostasis and cell fate decisions. The prostate gland is hormonally regulated, requiring steroids for proliferation and differentiation of secretory luminal cells. Vitamin D deficiency is associated with an increased risk of lethal prostate cancer, which exhibits a dedifferentiated pathology, linking vitamin D sufficiency to epithelial differentiation. To determine vitamin D regulation of prostatic epithelial differentiation, patient-derived benign prostate epithelial organoids were grown in vitamin D-deficient or -sufficient conditions. Organoids were assessed by phenotype and single-cell RNA sequencing. Mechanistic validation demonstrated that vitamin D sufficiency promoted organoid growth and accelerated differentiation by inhibiting canonical Wnt activity and suppressing Wnt family member DKK3. Wnt and DKK3 were also reduced by vitamin D in prostate tissue explants by spatial transcriptomics. Wnt dysregulation is a known contributor to aggressive prostate cancer, thus findings further link vitamin D deficiency to lethal disease.

The Transcriptomic Landscape of Prostate Cancer Development and Progression: An Integrative Analysis

Next-generation sequencing of primary tumors is now standard for transcriptomic studies, but microarray-based data still constitute the majority of available information on other clinically valuable samples, including archive material. Using prostate cancer (PC) as a model, we developed a robust analytical framework to integrate data across different technical platforms and disease subtypes to connect distinct disease stages and reveal potentially relevant genes not identifiable from single studies alone. We reconstructed the molecular profile of PC to yield the first comprehensive insight into its development, by tracking changes in mRNA levels from normal prostate to high-grade prostatic intraepithelial neoplasia, and metastatic disease. A total of nine previously unreported stage-specific candidate genes with prognostic significance were also found. Here, we integrate gene expression data from disparate sample types, disease stages and technical platforms into one coherent whole, to give a global view of the expression changes associated with the development and progression of PC from normal tissue through to metastatic disease. Summary and individual data are available online at the Prostate Integrative Expression Database (PIXdb), a user-friendly interface designed for clinicians and laboratory researchers to facilitate translational research.

Investigation of Candidate Genes and Pathways in Basal/TNBC Patients by Integrated Analysis

PURPOSE

This study aims to identify the key pathway and related genes and to further explore the potential molecular mechanisms of triple negative breast cancer (TNBC).

METHODS

The transcriptome data and clinical information of breast cancer patients were downloaded from the TCGA database, including 94 cases of paracancerous tissue, 225 cases of Basal like type, 151 cases of Her2 type, 318 cases of Luminal type A, 281 cases of Luminal type B, and 89 cases of Normal Like type. The differentially expressed genes (DEGs) were identified based on the criteria of |logFC|≥1.5 and adjust P < 0.001.Their functions were annotated by gene ontology (GO) analysis and Kyoto Encyclopedia of differentially expressed genes & Genomes (KEGG) pathway analysis. Cox regression univariate analysis and Kaplan-Meier survival curves (Log-rank method) were used for survival analysis. FOXD1, DLL3 and LY6D were silenced in breast cancer cell lines, and cell viability was assessed by CCK-8 assay. Further, the expression of FOXD1, DLL3 and LY6D were explored by immunohistochemistry on triple negative breast tumor tissue and normal breast tissue.

RESULTS

A total of 533 DEGs were identified. Functional annotation showed that DEGs were significantly enriched in intermediate filament cytoskeleton, DNA-binding transcription activator activity, epidermis development, and Neuroactive ligand-receptor interaction. Survival analysis found that FOXD1, DLL3, and LY6D showed significant correlation with the prognosis of patients with the Basal-like type (P < 0.05). CCK-8 assay showed that compared with Doxorubicin alone group, the cytotoxicity of Doxorubicin combined with siRNA-knockdown of FOXD1, DLL3, or LY6D was much significant.

CONCLUSION

The DEGs and their enriched functions and pathways identified in this study contribute to the understanding of the molecular mechanisms of TNBC. In addition, FOXD1, DLL3, and LY6D may be defined as the prognostic markers and potential therapeutic targets for TNBC patients.

Lung Adenocarcinoma Mouse Models Based on Orthotopic Transplantation of Syngeneic Tumor-Initiating Cells Expressing EpCAM, SCA-1, and Ly6d

Somatic mutations in EGFR and KRAS as well as chromosome rearrangements affecting ALK, ROS1, and RET have been identified in human lung adenocarcinoma (LUAD). We here developed organoid-based orthotopic and syngeneic mouse models for studies of the pathogenesis and treatment of LUAD. We isolated EpCAM-positive epithelial cells from mouse lungs and cultured them as organoids to maintain epithelial stem cell properties. These cells were transformed by KRAS(G12V) or EML4-ALK and then transplanted via the trachea into the lungs of the syngeneic mice, where they formed tumors that expressed the lung lineage marker TTF-1 and which closely recapitulated the pathology of human LUAD. Treatment with crizotinib suppressed the growth of tumors formed by the EML4-ALK-expressing lung epithelial cells in a subcutaneous transplantation model. Organoid culture of normal lung epithelial cells resulted in enrichment of EpCAM+SCA-1(Ly6a)+ cells as well as in that of cells expressing another member of the Ly6 protein family, Ly6d, which was found to be required for the growth of the LUAD-initiating cells expressing KRAS(G12V) or EML4-ALK. We also found that a high expression level of LY6D was associated with poor prognosis in human LUAD. Our results thus suggest that LY6D is a potential lung cancer stem cell marker.

Single-cell analysis supports a luminal-neuroendocrine transdifferentiation in human prostate cancer

Neuroendocrine prostate cancer is one of the most aggressive subtypes of prostate tumor. Although much progress has been made in understanding the development of neuroendocrine prostate cancer, the cellular architecture associated with neuroendocrine differentiation in human prostate cancer remain incompletely understood. Here, we use single-cell RNA sequencing to profile the transcriptomes of 21,292 cells from needle biopsies of 6 castration-resistant prostate cancers. Our analyses reveal that all neuroendocrine tumor cells display a luminal-like epithelial phenotype. In particular, lineage trajectory analysis suggests that focal neuroendocrine differentiation exclusively originate from luminal-like malignant cells rather than basal compartment. Further tissue microarray analysis validates the generality of the luminal phenotype of neuroendocrine cells. Moreover, we uncover neuroendocrine differentiation-associated gene signatures that may help us to further explore other intrinsic molecular mechanisms deriving neuroendocrine prostate cancer. In summary, our single-cell study provides direct evidence into the cellular states underlying neuroendocrine transdifferentiation in human prostate cancer.

Olfactomedin 4 mediation of prostate stem/progenitor-like cell proliferation and differentiation via MYC

Olfactomedin 4 (OLFM4) is expressed in normal prostate epithelial cells and immortalized normal human prostate epithelial cells (RWPE1), but the identity of OLFM4-expressing cells within these populations and OLFM4's physiological functions in these cells have not been elucidated. Using single-cell RNA sequencing analysis, we found here that OLFM4 was expressed in multiple stem/progenitor-like cell populations in both the normal prostate epithelium and RWPE1 cells and was frequently co-expressed with KRT13 and LY6D in RWPE1 cells. Functionally, OLFM4-knockout RWPE1 cells exhibited enhanced proliferation of the stem/progenitor-like cell population, shifts stem/progenitor-like cell division to favor symmetric division and differentiated into higher levels PSA expression cells in organoid assays compared with OLFM4-wild RWPE1 cells. Bulk-cell RNA sequencing analysis pinpointed that cMYC expression were enhanced in the OLFM4-knockout RWPE1 cells compared with OLFM4-wild cells. Molecular and signaling pathway studies revealed an increase in the WNT/APC/MYC signaling pathway gene signature, as well as that of MYC target genes that regulate multiple biological processes, in OLFM4-knockout RWPE1 cells. These findings indicated that OLFM4 is co-expressed with multiple stem/progenitor cell marker genes in prostate epithelial cells and acts as a novel mediator in prostate stem/progenitor cell proliferation and differentiation.

Fate and functional roles of Prominin 1+ cells in liver injury and cancer

Prominin 1 (PROM1) is one of a few clinically relevant progenitor markers in human alcoholic hepatitis (AH) and hepatocellular carcinoma (HCC), and mouse liver tumor initiating stem cell-like cells (TICs). However, the origin, fate and functions of PROM1⁺ cells in AH and HCC are unknown. Here we show by genetic lineage tracing that PROM1⁺ cells are derived in part from hepatocytes in AH and become tumor cells in mice with diethyl nitrosamine (DEN)-initiated, Western alcohol diet-promoted liver tumorigenesis. Our RNA sequencing analysis of mouse PROM1⁺ cells, reveals transcriptomic landscapes indicative of their identities as ductular reaction progenitors (DRPs) and TICs. Indeed, single-cell RNA sequencing reveals two subpopulations of Prom1⁺ Afp^– DRPs and Prom1⁺ Afp⁺ TICs in the DEN-WAD model. Integrated bioinformatic analysis identifies Discodin Domain Receptor 1 (DDR1) as a uniquely upregulated and patient-relevant gene in PROM1⁺ cells in AH and HCC. Translational relevance of DDR1 is supported by its marked elevation in HCC which is inversely associated with patient survival. Further, knockdown of Ddr1 suppresses the growth of TICs and TIC-derived tumor growth in mice. These results suggest the importance of PROM1⁺ cells in the evolution of liver cancer and DDR1 as a potential driver of this process.

Lymphocyte antigen 6 superfamily member D is a marker of urothelial and squamous differentiation: implications for risk stratification of bladder cancer

Background

Screening across a multitude of normal and malignant tissues revealed an enhanced expression of lymphocyte antigen 6 superfamily member D (LY6D) in squamous epithelium and urothelium, as well as in malignancies derived therefrom. The aim of this study was to further delineate the protein expression of LY6D in urothelial bladder cancer, with particular attention to its relationship with clinicopathological characteristics and patient outcome.

Methods

Immunohistochemical expression of LY6D was assessed in tissue microarrays with urothelial bladder cancer tumours from three independent patient cohorts; one with transurethral resection of the bladder (TURB) specimens of mixed tumour stages from 110 consecutive cases, one with tumours of mixed stages from 260 incident cases in a population-based cohort, and one with paired TURB specimens, resected tumours and a subset of lymph node metastases from 145 patients with muscle-invasive bladder cancer (MIBC). Chi-square and non-parametric tests were applied to examine associations of LY6D expression with clinicopathological characteristics. Kaplan-Meier and Cox regression analyses were applied to examine 5-year overall survival (OS) and recurrence free survival (RFS) in relation to LY6D expression.

Results

In the two cohorts with mixed stages, positive LY6D expression was denoted in 63 and 64% of the cases, respectively, and found to be significantly higher in low-grade and less invasive tumours. Negative LY6D expression was significantly associated with a reduced 5-year OS, although not independently of established prognostic factors. In the population-based cohort, LY6D expression was higher in tumours with squamous differentiation and lower in other variant histologies compared to pure urothelial tumours, and the association of LY6D expression with survival was somewhat enhanced after exclusion of the former. LY6D expression was generally lower in the MIBC cohort, and even more reduced in resected tumours compared to TURB specimens in patients who had not received neoadjuvant chemotherapy. There were no significant associations between LY6D expression and RFS, neither allover nor in relation to neoadjuvant chemotherapy.

Conclusion

LY6D is a marker of urothelial and squamous differentiation that may add useful diagnostic and prognostic information to better guide the clinical management of bladder cancer, given that the presence of variant histology is taken into account.

RUNX1 marks a luminal castration-resistant lineage established at the onset of prostate development

The characterization of prostate epithelial hierarchy and lineage heterogeneity is critical to understand its regenerative properties and malignancies. Here, we report that the transcription factor RUNX1 marks a specific subpopulation of proximal luminal cells (PLCs), enriched in the periurethral region of the developing and adult mouse prostate, and distinct from the previously identified NKX3.1+ luminal castration-resistant cells. Using scRNA-seq profiling and genetic lineage tracing, we show that RUNX1+ PLCs are unaffected by androgen deprivation, and do not contribute to the regeneration of the distal luminal compartments. Furthermore, we demonstrate that a transcriptionally similar RUNX1+ population emerges at the onset of embryonic prostate specification to populate the proximal region of the ducts. Collectively, our results reveal that RUNX1+ PLCs is an intrinsic castration-resistant and self-sustained lineage that emerges early during prostate development and provide new insights into the lineage relationships of the prostate epithelium.

A single-cell atlas of the mouse and human prostate reveals heterogeneity and conservation of epithelial progenitors

Understanding the cellular constituents of the prostate is essential for identifying the cell of origin for prostate adenocarcinoma. Here, we describe a comprehensive single-cell atlas of the adult mouse prostate epithelium, which displays extensive heterogeneity. We observe distal lobe-specific luminal epithelial populations (LumA, LumD, LumL, and LumV), a proximally enriched luminal population (LumP) that is not lobe-specific, and a periurethral population (PrU) that shares both basal and luminal features. Functional analyses suggest that LumP and PrU cells have multipotent progenitor activity in organoid formation and tissue reconstitution assays. Furthermore, we show that mouse distal and proximal luminal cells are most similar to human acinar and ductal populations, that a PrU-like population is conserved between species, and that the mouse lateral prostate is most similar to the human peripheral zone. Our findings elucidate new prostate epithelial progenitors, and help resolve long-standing questions about anatomical relationships between the mouse and human prostate.