Single-Cell RNA-seq Reveals a Subpopulation of Prostate Cancer Cells with Enhanced Cell-Cycle-Related Transcription and Attenuated Androgen Response

Integrated single-cell transcriptomic analyses identify a novel lineage plasticity-related cancer cell type involved in prostate cancer progression

BACKGROUND

Cancer cell plasticity is the ability of neoplastic cells to alter their identity and acquire new biological properties under microenvironmental pressures. In prostate cancer (PCa), lineage plasticity often results in therapy resistance and trans-differentiation to neuroendocrine (NE) lineage. However, identifying the cancer cells harboring lineage plasticity-related status remains challenging.

METHODS

Based on 13 multi-center human PCa bulk transcriptomic cohorts (samples = 3314) and 9 bulk transcriptomic datasets derived from PCa experimental models, we established an integrated lineage plasticity-related gene signature, termed LPSig. Leveraging this gene signature, AUCell enrichment analysis was applied to identify the cell population with high lineage plasticity from a comprehensive single-cell RNA-sequencing (scRNA-seq) meta-atlas assembled by us, which consisted of 10 public human PCa scRNA-seq datasets (samples = 93, cells = 222,529). Moreover, additional scRNA-seq dataset of human PCa, multiplex immunohistochemistry staining for human PCa tissues, in vitro and in vivo functional experiments, as well as qPCR and Western blot analyses were employed to validate our findings.

FINDINGS

We found that LPSig could finely capture the dynamics of tumor lineage plasticity throughout the progression of PCa, accurately estimating the status of lineage plasticity. Based on LPSig, we identified a previously undefined minority population of lineage plasticity-related PCa cells (LPCs) from the human PCa scRNA-seq meta-atlas assembled by this study. Furthermore, in-depth dissection revealed pivotal roles of LPCs in trans-differentiation, tumor recurrence, and poor patient survival during PCa progression. Furthermore, we identified HMMR as a representative cell surface marker for LPCs, which was validated using additional scRNA-seq datasets and multiplexed immunohistochemistry. Moreover, HMMR was transcriptionally inhibited by androgen receptor (AR), and was required for the aggressive adenocarcinoma features and NE phenotype.

INTERPRETATION

Our study uncovers a novel population of lineage plasticity-related cells with low AR activity, stemness-like traits, and elevated HMMR expression, that may facilitate poor prognosis in PCa.

FUNDING

This work was supported by National Key R&D Program of China (2022YFA0807000), National Natural Science Foundation of China (82160584), Advanced Prostate Cancer Diagnosis and Treatment Technology Innovation Team of Kunming Medical University (CXTD202216), and Reserve Talents of Young and Middle-aged Academic Leaders in Yunnan Province (202105AC160013).

Changes in immunophenotypes after neoadjuvant endocrine therapy for prostate cancer and their clinical significance

Background

To investigate changes in the immunophenotypes of androgen receptor (AR), prostate-specific antigen (PSA), synaptophysin (Syn), chromogranin A (CgA), p53 and Ki-67 after neoadjuvant endocrine therapy (NET) for prostate cancer (PCa) and to analyze their clinical significance.

Methods

Paired paraffin samples were collected from 40 PCa patients before and after NET, and immunohistochemistry were used to detect AR, PSA, Syn, CgA, p53 and Ki-67 expression. Based on The Cancer Genome Atlas (TCGA), Kaplan‒Meier survival curves were plotted for analysis of PSA and Ki-67 expression in relation to progression-free survival (PFS).

Results

After NET, the mean scores for PSA and Ki-67 expression in PCa patients were lower than those before NET (P < 0.05), while the mean scores for Syn and CgA expression were higher than those before NET (P < 0.05). The mean Gleason score and WHO/ISUP (World Health Organization/International Society of Urological Pathology) grade after NET were lower than those before NET (P < 0.05). In PCa patients who had not yet received NET, PSA expression correlated positively with Gleason score and WHO/ISUP grade and negatively with Ki-67 expression (P < 0.05); p53 expression correlated negatively with Gleason score and WHO/ISUP grade (P < 0.05). TCGA showed that PFS was lower in PCa patients with high PSA and Ki-67 expression (P < 0.05).

Conclusions

PSA and Ki-67 protein expressions decreased significantly in PCa patients after NET and can be used as biological markers for prognostic assessment of PCa patients. NETs may induce a neuroendocrine (NE) phenotype in PCa. Monitoring the immunophenotypes of PCa patients after NET may inform assessment of efficacy and prognosis.

Correlation between abdominal computed tomography signs and postoperative prognosis for patients with colorectal cancer

BACKGROUND

Patients with different stages of colorectal cancer (CRC) exhibit different abdominal computed tomography (CT) signs. Therefore, the influence of CT signs on CRC prognosis must be determined.

AIM

To observe abdominal CT signs in patients with CRC and analyze the correlation between the CT signs and postoperative prognosis.

METHODS

The clinical history and CT imaging results of 88 patients with CRC who underwent radical surgery at Xingtan Hospital Affiliated to Shunde Hospital of Southern Medical University were retrospectively analyzed. Univariate and multivariate Cox regression analyses were used to explore the independent risk factors for postoperative death in patients with CRC. The three-year survival rate was analyzed using the Kaplan-Meier curve, and the correlation between postoperative survival time and abdominal CT signs in patients with CRC was analyzed using Spearman correlation analysis.

RESULTS

For patients with CRC, the three-year survival rate was 73.86%. The death group exhibited more severe characteristics than the survival group. A multivariate Cox regression model analysis showed that body mass index (BMI), degree of periintestinal infiltration, tumor size, and lymph node CT value were independent factors influencing postoperative death (P < 0.05 for all). Patients with characteristics typical to the death group had a low three-year survival rate (log-rank χ 2 = 66.487, 11.346, 12.500, and 27.672, respectively, P < 0.05 for all). The survival time of CRC patients was negatively correlated with BMI, degree of periintestinal infiltration, tumor size, lymph node CT value, mean tumor long-axis diameter, and mean tumor short-axis diameter (r = -0.559, 0.679, -0.430, -0.585, -0.425, and -0.385, respectively, P < 0.05 for all). BMI was positively correlated with the degree of periintestinal invasion, lymph node CT value, and mean tumor short-axis diameter (r = 0.303, 0.431, and 0.437, respectively, P < 0.05 for all).

CONCLUSION

The degree of periintestinal infiltration, tumor size, and lymph node CT value are crucial for evaluating the prognosis of patients with CRC.

Cancer stem cells: advances in knowledge and implications for cancer therapy

Cancer stem cells (CSCs), a small subset of cells in tumors that are characterized by self-renewal and continuous proliferation, lead to tumorigenesis, metastasis, and maintain tumor heterogeneity. Cancer continues to be a significant global disease burden. In the past, surgery, radiotherapy, and chemotherapy were the main cancer treatments. The technology of cancer treatments continues to develop and advance, and the emergence of targeted therapy, and immunotherapy provides more options for patients to a certain extent. However, the limitations of efficacy and treatment resistance are still inevitable. Our review begins with a brief introduction of the historical discoveries, original hypotheses, and pathways that regulate CSCs, such as WNT/β-Catenin, hedgehog, Notch, NF-κB, JAK/STAT, TGF-β, PI3K/AKT, PPAR pathway, and their crosstalk. We focus on the role of CSCs in various therapeutic outcomes and resistance, including how the treatments affect the content of CSCs and the alteration of related molecules, CSCs-mediated therapeutic resistance, and the clinical value of targeting CSCs in patients with refractory, progressed or advanced tumors. In summary, CSCs affect therapeutic efficacy, and the treatment method of targeting CSCs is still difficult to determine. Clarifying regulatory mechanisms and targeting biomarkers of CSCs is currently the mainstream idea.

Analysis of the heterogeneity and complexity of murine extraorbital lacrimal gland via single-cell RNA sequencing

PURPOSE

The lacrimal gland is essential for maintaining ocular surface health and avoiding external damage by secreting an aqueous layer of the tear film. However, a healthy lacrimal gland's inventory of cell types and heterogeneity remains understudied.

METHODS

Here, 10X Genome-based single-cell RNA sequencing was used to generate an unbiased classification of cellular diversity in the extraorbital lacrimal gland (ELG) of C57BL/6J mice. From 43,850 high-quality cells, we produced an atlas of cell heterogeneity and defined cell types using classic marker genes. The possible functions of these cells were analyzed through bioinformatics analysis. Additionally, the CellChat was employed for a preliminary analysis of the cell-cell communication network in the ELG.

RESULTS

Over 37 subclasses of cells were identified, including seven types of glandular epithelial cells, three types of fibroblasts, ten types of myeloid-derived immune cells, at least eleven types of lymphoid-derived immune cells, and five types of vascular-associated cell subsets. The cell-cell communication network analysis revealed that fibroblasts and immune cells play a pivotal role in the dense intercellular communication network within the mouse ELG.

CONCLUSIONS

This study provides a comprehensive transcriptome atlas and related database of the mouse ELG.

Comprehensive data mining reveals RTK/RAS signaling pathway as a promoter of prostate cancer lineage plasticity through transcription factors and CNV

Prostate cancer lineage plasticity is a key driver in the transition to neuroendocrine prostate cancer (NEPC), and the RTK/RAS signaling pathway is a well-established cancer pathway. Nevertheless, the comprehensive link between the RTK/RAS signaling pathway and lineage plasticity has received limited investigation. In particular, the intricate regulatory network governing the interplay between RTK/RAS and lineage plasticity remains largely unexplored. The multi-omics data were clustered with the coefficient of argument and neighbor joining algorithm. Subsequently, the clustered results were analyzed utilizing the GSEA, gene sets related to stemness, multi-lineage state datasets, and canonical cancer pathway gene sets. Finally, a comprehensive exploration of the data based on the ssGSEA, WGCNA, GSEA, VIPER, prostate cancer scRNA-seq data, and the GPSAdb database was conducted. Among the six modules in the clustering results, there are 300 overlapping genes, including 3 previously unreported prostate cancer genes that were validated to be upregulated in prostate cancer through RT-qPCR. Function Module 6 shows a positive correlation with prostate cancer cell stemness, multi-lineage states, and the RTK/RAS signaling pathway. Additionally, the 19 leading-edge genes of the RTK/RAS signaling pathway promote prostate cancer lineage plasticity through a complex network of transcriptional regulation and copy number variations. In the transcriptional regulation network, TP63 and FOXO1 act as suppressors of prostate cancer lineage plasticity, whereas RORC exerts a promoting effect. This study provides a comprehensive perspective on the role of the RTK/RAS pathway in prostate cancer lineage plasticity and offers new clues for the treatment of NEPC.

scNovel: a scalable deep learning-based network for novel rare cell discovery in single-cell transcriptomics

Single-cell RNA sequencing has achieved massive success in biological research fields. Discovering novel cell types from single-cell transcriptomics has been demonstrated to be essential in the field of biomedicine, yet is time-consuming and needs prior knowledge. With the unprecedented boom in cell atlases, auto-annotation tools have become more prevalent due to their speed, accuracy and user-friendly features. However, existing tools have mostly focused on general cell-type annotation and have not adequately addressed the challenge of discovering novel rare cell types. In this work, we introduce scNovel, a powerful deep learning-based neural network that specifically focuses on novel rare cell discovery. By testing our model on diverse datasets with different scales, protocols and degrees of imbalance, we demonstrate that scNovel significantly outperforms previous state-of-the-art novel cell detection models, reaching the most AUROC performance(the only one method whose averaged AUROC results are above 94%, up to 16.26% more comparing to the second-best method). We validate scNovel's performance on a million-scale dataset to illustrate the scalability of scNovel further. Applying scNovel on a clinical COVID-19 dataset, three potential novel subtypes of Macrophages are identified, where the COVID-related differential genes are also detected to have consistent expression patterns through deeper analysis. We believe that our proposed pipeline will be an important tool for high-throughput clinical data in a wide range of applications.

SARIFA as a new histopathological biomarker is associated with adverse clinicopathological characteristics, tumor-promoting fatty-acid metabolism, and might predict a metastatic pattern in pT3a prostate cancer

BACKGROUND

Recently, we introduced Stroma-AReactive-Invasion-Front-Areas (SARIFA) as a novel hematoxylin-eosin (H&E)-based histopathologic prognostic biomarker for various gastrointestinal cancers, closely related to lipid metabolism. To date, no studies on SARIFA, which is defined as direct tumor-adipocyte-interaction, beyond the alimentary tract exist. Hence, the objective of our current investigation was to study the significance of SARIFA in pT3a prostate cancer (PCa) and explore its association with lipid metabolism in PCa as lipid metabolism plays a key role in PCa development and progression.

METHODS

To this end, we evaluated SARIFA-status in 301 radical prostatectomy specimens and examined the relationship between SARIFA-status, clinicopathological characteristics, overall survival, and immunohistochemical expression of FABP4 and CD36 (proteins closely involved in fatty-acid metabolism). Additionally, we investigated the correlation between SARIFA and biochemical recurrence-free survival (BRFS) and PSMA-positive recurrences in PET/CT imaging in a patient subgroup. Moreover, a quantitative SARIFA cut-off was established to further understand the underlying tumor biology.

RESULTS

SARIFA positivity occurred in 59.1% (n = 178) of pT3a PCas. Our analysis demonstrated that SARIFA positivity is strongly associated with established high-risk features, such as R1 status, extraprostatic extension, and higher initial PSA values. Additionally, we observed an upregulation of immunohistochemical CD36 expression specifically at SARIFAs (p = 0.00014). Kaplan-Meier analyses revealed a trend toward poorer outcomes, particularly in terms of BRFS (p = 0.1). More extensive tumor-adipocyte interaction, assessed as quantity-dependent SARIFA-status on H&E slides, is also significantly associated with high-risk features, such as lymph node metastasis, and seems to be associated with worse survival outcomes (p = 0.16). Moreover, SARIFA positivity appeared to be linked to more distant lymph node and bone metastasis, although statistical significance was slightly not achieved (both p > 0.05).

CONCLUSIONS

This is the first study to introduce SARIFA as easy-and-fast-to-assess H&E-based biomarker in locally advanced PCa. SARIFA as the histopathologic correlate of a distinct tumor biology, closely related to lipid metabolism, could pave the way to a more detailed patient stratification and to the development of novel drugs targeting lipid metabolism in pT3a PCa. On the basis of this biomarker discovery study, further research efforts on the prognostic and predictive role of SARIFA in PCa can be designed.

De novo identification of expressed cancer somatic mutations from single-cell RNA sequencing data

Identifying expressed somatic mutations from single-cell RNA sequencing data de novo is challenging but highly valuable. We propose RESA - Recurrently Expressed SNV Analysis, a computational framework to identify expressed somatic mutations from scRNA-seq data. RESA achieves an average precision of 0.77 on three in silico spike-in datasets. In extensive benchmarking against existing methods using 19 datasets, RESA consistently outperforms them. Furthermore, we applied RESA to analyze intratumor mutational heterogeneity in a melanoma drug resistance dataset. By enabling high precision detection of expressed somatic mutations, RESA substantially enhances the reliability of mutational analysis in scRNA-seq. RESA is available at https://github.com/ShenLab-Genomics/RESA .

CamoTSS: analysis of alternative transcription start sites for cellular phenotypes and regulatory patterns from 5' scRNA-seq data

Five-prime single-cell RNA-seq (scRNA-seq) has been widely employed to profile cellular transcriptomes, however, its power of analysing transcription start sites (TSS) has not been fully utilised. Here, we present a computational method suite, CamoTSS, to precisely identify TSS and quantify its expression by leveraging the cDNA on read 1, which enables effective detection of alternative TSS usage. With various experimental data sets, we have demonstrated that CamoTSS can accurately identify TSS and the detected alternative TSS usages showed strong specificity in different biological processes, including cell types across human organs, the development of human thymus, and cancer conditions. As evidenced in nasopharyngeal cancer, alternative TSS usage can also reveal regulatory patterns including systematic TSS dysregulations.

Integrative multi-omics analysis unveils stemness-associated molecular subtypes in prostate cancer and pan-cancer: prognostic and therapeutic significance

BACKGROUND

Prostate cancer (PCA) is the fifth leading cause of cancer-related deaths worldwide, with limited treatment options in the advanced stages. The immunosuppressive tumor microenvironment (TME) of PCA results in lower sensitivity to immunotherapy. Although molecular subtyping is expected to offer important clues for precision treatment of PCA, there is currently a shortage of dependable and effective molecular typing methods available for clinical practice. Therefore, we aim to propose a novel stemness-based classification approach to guide personalized clinical treatments, including immunotherapy.

METHODS

An integrative multi-omics analysis of PCA was performed to evaluate stemness-level heterogeneities. Unsupervised hierarchical clustering was used to classify PCAs based on stemness signature genes. To make stemness-based patient classification more clinically applicable, a stemness subtype predictor was jointly developed by using four PCA datasets and 76 machine learning algorithms.

RESULTS

We identified stemness signatures of PCA comprising 18 signaling pathways, by which we classified PCA samples into three stemness subtypes via unsupervised hierarchical clustering: low stemness (LS), medium stemness (MS), and high stemness (HS) subtypes. HS patients are sensitive to androgen deprivation therapy, taxanes, and immunotherapy and have the highest stemness, malignancy, tumor mutation load (TMB) levels, worst prognosis, and immunosuppression. LS patients are sensitive to platinum-based chemotherapy but resistant to immunotherapy and have the lowest stemness, malignancy, and TMB levels, best prognosis, and the highest immune infiltration. MS patients represent an intermediate status of stemness, malignancy, and TMB levels with a moderate prognosis. We further demonstrated that these three stemness subtypes are conserved across pan-tumor. Additionally, the 9-gene stemness subtype predictor we developed has a comparable capability to 18 signaling pathways to make tumor diagnosis and to predict tumor recurrence, metastasis, progression, prognosis, and efficacy of different treatments.

CONCLUSIONS

The three stemness subtypes we identified have the potential to be a powerful tool for clinical tumor molecular classification in PCA and pan-cancer, and to guide the selection of immunotherapy or other sensitive treatments for tumor patients.

Identification of Niche-Specific Gene Signatures between Malignant Tumor Microenvironments by Integrating Single Cell and Spatial Transcriptomics Data

The tumor microenvironment significantly affects the transcriptomic states of tumor cells. Single-cell RNA sequencing (scRNA-seq) helps elucidate the transcriptomes of individual cancer cells and their neighboring cells. However, cell dissociation results in the loss of information on neighboring cells. To address this challenge and comprehensively assess the gene activity in tissue samples, it is imperative to integrate scRNA-seq with spatial transcriptomics. In our previous study on physically interacting cell sequencing (PIC-seq), we demonstrated that gene expression in single cells is affected by neighboring cell information. In the present study, we proposed a strategy to identify niche-specific gene signatures by harmonizing scRNA-seq and spatial transcriptomic data. This approach was applied to the paired or matched scRNA-seq and Visium platform data of five cancer types: breast cancer, gastrointestinal stromal tumor, liver hepatocellular carcinoma, uterine corpus endometrial carcinoma, and ovarian cancer. We observed distinct gene signatures specific to cellular niches and their neighboring counterparts. Intriguingly, these niche-specific genes display considerable dissimilarity to cell type markers and exhibit unique functional attributes independent of the cancer types. Collectively, these results demonstrate the potential of this integrative approach for identifying novel marker genes and their spatial relationships.

DLGAP5 Regulates the Proliferation, Migration, Invasion, and Cell Cycle of Breast Cancer Cells via the JAK2/STAT3 Signaling Axis

The aim of this study was to discover new biomarkers to detect breast cancer (BC), which is an aggressive cancer with a high mortality rate. In this study, bioinformatic analyses (differential analysis, weighted gene co-expression network analysis, and machine learning) were performed to identify potential candidate genes for BC to study their molecular mechanisms. Furthermore, Quantitative Real-time PCR and immunohistochemistry assays were used to examine the protein and mRNA expression levels of a particular candidate gene (DLGAP5). And the effects of DLGAP5 on cell proliferation, migration, invasion, and cell cycle were further assessed using the Cell Counting Kit-8 assay, colony formation, Transwell, wound healing, and flow cytometry assays. Moreover, the changes in the JAK2/STAT3 signaling-pathway-related proteins were detected by Western Blot. A total of 44 overlapping genes were obtained by differential analysis and weighted gene co-expression network analysis, of which 25 genes were found in the most tightly connected cluster. Finally, NEK2, CKS2, UHRF1, DLGAP5, and FAM83D were considered as potential biomarkers of BC. Moreover, DLGAP5 was highly expressed in BC. The down-regulation of DLGAP5 may inhibit the proliferation, migration, invasion, and cell cycle of BC cells, and the opposite was true for DLGAP5 overexpression. Correspondingly, silencing or overexpression of the DLGAP5 gene inhibited or activated the JAK2/STAT3 signaling pathway, respectively. DLGAP5, as a potential biomarker of BC, may impact the cell proliferation, migration, invasion, cell cycle, and BC development by modulating the JAK2/STAT3 signaling pathway.

A scalable sparse neural network framework for rare cell type annotation of single-cell transcriptome data

Automatic cell type annotation methods are increasingly used in single-cell RNA sequencing (scRNA-seq) analysis due to their fast and precise advantages. However, current methods often fail to account for the imbalance of scRNA-seq datasets and ignore information from smaller populations, leading to significant biological analysis errors. Here, we introduce scBalance, an integrated sparse neural network framework that incorporates adaptive weight sampling and dropout techniques for auto-annotation tasks. Using 20 scRNA-seq datasets with varying scales and degrees of imbalance, we demonstrate that scBalance outperforms current methods in both intra- and inter-dataset annotation tasks. Additionally, scBalance displays impressive scalability in identifying rare cell types in million-level datasets, as shown in the bronchoalveolar cell landscape. scBalance is also significantly faster than commonly used tools and comes in a user-friendly format, making it a superior tool for scRNA-seq analysis on the Python-based platform.

HER2 as a potential therapeutic target on quiescent prostate cancer cells

Quiescent prostate cancer (PCa) cells are common in tumors but are often resistant to chemotherapy. Quiescent PCa cells are also enriched for a stem-like tumor initiating population, and can lead to recurrence after dormancy. Unfortunately, quiescent PCa cells are difficult to identify and / or target with treatment in part because the relevant markers are intracellular and regulated by protein stability. We addressed this problem by utilizing PCa cells expressing fluorescent markers for CDKN1B (p27) and CDT1, which can separate viable PCa cells into G0, G1, or combined S/G2/M populations. We used FACS to collect G1 and G0 PC3 PCa cells, isolated membrane proteins, and analyzed protein abundance in G0 vs G1 cells by gas chromatography mass spectrometry. Enrichment analysis identified nucleocytoplasmic transport as the most significantly different pathway. To identify cell surface proteins potentially identifying quiescent PCa cells for future patient samples or for antibody based therapeutic research, we focused on differentially abundant plasma membrane proteins, and identified ERBB2 (HER2) as a cell surface protein enriched on G0 PCa cells. High HER2 on the cell membrane is associated with quiescence in PCa cells and likely induced by the bone microenvironment. Using a drug conjugated anti-HER2 antibody (trastuzumab emtansine) in a mouse PCa xenograft model delayed metastatic tumor growth, suggesting approaches that target HER2-high cells may be beneficial in treating PCa. We propose that HER2 is deserving of further study in PCa as a target on quiescent cells to prevent recurrence, decrease chemotherapy resistance, or eradicate minimal residual disease.

Small-Molecule Inhibition of Androgen Receptor Dimerization as a Strategy against Prostate Cancer

The clinically used androgen receptor (AR) antagonists for the treatment of prostate cancer (PCa) are all targeting the AR ligand binding pocket (LBP), resulting in various drug-resistant problems. Therefore, a new strategy to combat PCa is urgently needed. Enlightened by the gain-of-function mutations of androgen insensitivity syndrome, we discovered for the first time small-molecule antagonists toward a prospective pocket on the AR dimer interface named the dimer interface pocket (DIP) via molecular dynamics (MD) simulation, structure-based virtual screening, structure-activity relationship exploration, and bioassays. The first-in-class antagonist M17-B15 targeting the DIP is capable of effectively disrupting AR self-association, thereby suppressing AR signaling. Furthermore, M17-B15 exhibits extraordinary anti-PCa efficacy in vitro and also in mouse xenograft tumor models, demonstrating that AR dimerization disruption by small molecules targeting the DIP is a novel and valid strategy against PCa.

The Molecular Biology of Prostate Cancer Stem Cells: From the Past to the Future

Prostate cancer (PCa) continues to rank as the second leading cause of cancer-related mortality in western countries, despite the golden treatment using androgen deprivation therapy (ADT) or anti-androgen therapy. With decades of research, scientists have gradually realized that the existence of prostate cancer stem cells (PCSCs) successfully explains tumor recurrence, metastasis and therapeutic failure of PCa. Theoretically, eradication of this small population may improve the efficacy of current therapeutic approaches and prolong PCa survival. However, several characteristics of PCSCs make their diminishment extremely challenging: inherent resistance to anti-androgen and chemotherapy treatment, over-activation of the survival pathway, adaptation to tumor micro-environments, escape from immune attack and being easier to metastasize. For this end, a better understanding of PCSC biology at the molecular level will definitely inspire us to develop PCSC targeted approaches. In this review, we comprehensively summarize signaling pathways responsible for homeostatic regulation of PCSCs and discuss how to eliminate these fractional cells in clinical practice. Overall, this study deeply pinpoints PCSC biology at the molecular level and provides us some research perspectives.

SMURF: embedding single-cell RNA-seq data with matrix factorization preserving self-consistency

The advance in single-cell RNA-sequencing (scRNA-seq) sheds light on cell-specific transcriptomic studies of cell developments, complex diseases and cancers. Nevertheless, scRNA-seq techniques suffer from 'dropout' events, and imputation tools are proposed to address the sparsity. Here, rather than imputation, we propose a tool, SMURF, to extract the low-dimensional embeddings from cells and genes utilizing matrix factorization with a mixture of Poisson-Gamma divergent as objective while preserving self-consistency. SMURF exhibits feasible cell subpopulation discovery efficacy with obtained cell embeddings on replicated in silico and eight web lab scRNA datasets with ground truth cell types. Furthermore, SMURF can reduce the cell embedding to a 1D-oval space to recover the time course of cell cycle. SMURF can also serve as an imputation tool; the in silico data assessment shows that SMURF parades the most robust gene expression recovery power with low root mean square error and high Pearson correlation. Moreover, SMURF recovers the gene distribution for the WM989 Drop-seq data. SMURF is available at https://github.com/deepomicslab/SMURF.

Vacuolar protein sorting 35 (VPS35) acts as a tumor promoter via facilitating cell cycle progression in pancreatic ductal adenocarcinoma

Pancreatic ductal adenocarcinoma (PDAC) is insidious and highly malignant with extremely poor prognosis and drug resistance to current chemotherapies. Therefore, there is a critical need to investigate the molecular mechanism underlying PDAC progression to develop promising diagnostic and therapeutic interventions. In parallel, vacuolar protein sorting (VPS) proteins, involved in the sorting, transportation, and localization of membrane proteins, have gradually attracted the attention of researchers in the development of cancers. Although VPS35 has been reported to promote carcinoma progression, the specific molecular mechanism is still unclear. Here, we determined the impact of VPS35 on the tumorigenesis of PDAC and explored the underlying molecular mechanism. We performed a pan-cancer analysis of 46 VPS genes using RNAseq data from GTEx (control) and TCGA (tumor) and predicted potential functions of VPS35 in PDAC by enrichment analysis. Furthermore, cell cloning experiments, gene knockout, cell cycle analysis, immunohistochemistry, and other molecular and biochemical experiments were used to validate the function of VPS35. Consequently, VPS35 was found overexpressed in multiple cancers and correlated with the poor prognosis of PDAC. Meanwhile, we verified that VPS35 could modulate the cell cycle and promote tumor cell growth in PDAC. Collectively, we provide solid evidence that VPS35 facilitates the cell cycle progression as a critical novel target in PDAC clinical therapy.

Cancer stem cells: Recent insights and therapies

Tumors are intricate ecosystems containing malignant components that generate adaptive and evolutionarily driven abnormal tissues. Through self-renewal and differentiation, cancers are reconstructed by a dynamic subset of stem-like cells that enforce tumor heterogeneity and remodel the tumor microenvironment (TME). Through recent technology advances, we are now better equipped to investigate the fundamental role of cancer stem cells (CSCs) in cancer biology. In this review, we discuss the latest insights into characteristics, markers and mechanism of CSCs and describe the crosstalk between CSCs and other cells in TME. Additionally, we explore the performance of single-cell sequencing and spatial transcriptome analysis in CSCs studies and summarize the therapeutic strategies to eliminate CSCs, which could broaden the understanding of CSCs and exploit for therapeutic benefit.

Detection of genes with differential expression dispersion unravels the role of autophagy in cancer progression

The majority of gene expression studies focus on the search for genes whose mean expression is different between two or more populations of samples in the so-called "differential expression analysis" approach. However, a difference in variance in gene expression may also be biologically and physiologically relevant. In the classical statistical model used to analyze RNA-sequencing (RNA-seq) data, the dispersion, which defines the variance, is only considered as a parameter to be estimated prior to identifying a difference in mean expression between conditions of interest. Here, we propose to evaluate four recently published methods, which detect differences in both the mean and dispersion in RNA-seq data. We thoroughly investigated the performance of these methods on simulated datasets and characterized parameter settings to reliably detect genes with a differential expression dispersion. We applied these methods to The Cancer Genome Atlas datasets. Interestingly, among the genes with an increased expression dispersion in tumors and without a change in mean expression, we identified some key cellular functions, most of which were related to catabolism and were overrepresented in most of the analyzed cancers. In particular, our results highlight autophagy, whose role in cancerogenesis is context-dependent, illustrating the potential of the differential dispersion approach to gain new insights into biological processes and to discover new biomarkers.

SR9009 inhibits lethal prostate cancer subtype 1 by regulating the LXRα/FOXM1 pathway independently of REV-ERBs

Perturbations of the circadian clock are linked to multiple diseases, including cancers. Pharmacological activation of REV-ERB nuclear receptors, the core components of the circadian clock, has antitumor effects on various malignancies, while the impact of SR9009 on prostate cancer (PCa) remains unknown. Here, we found that SR9009 was specifically lethal to PCa cell lines but had no cytotoxic effect on prostate cells. SR9009 significantly inhibited colony formation, the cell cycle, and cell migration and promoted apoptosis in PCa cells. SR9009 treatment markedly inhibited prostate cancer subtype 1 (PCS1), the most lethal and aggressive PCa subtype, through FOXM1 pathway blockade, while it had no impacts on PCS2 and PCS3. Seven representative genes, including FOXM1, CENPA, CENPF, CDK1, CCNB1, CCNB2, and BIRC5, were identified as the shared genes involved in the FOXM1 pathway and PCS1. All of these genes were upregulated in PCa tissues, associated with worse clinicopathological outcomes and downregulated after SR9009 treatment. Nevertheless, knockdown or knockout of REV-ERB could not rescue the anticancer effect of SR9009 in PCa. Further analysis confirmed that it was LXRα rather than REV-ERBs which has been activated by SR9009. The expression levels of these seven genes were changed correspondingly after LXRα knockdown and SR9009 treatment. An in vivo study validated that SR9009 restrained tumor growth in 22RV1 xenograft models and inhibited FOXM1 and its targeted gene expression. In summary, SR9009 can serve as an effective treatment option for highly aggressive and lethal PCS1 tumors through mediating the LXRα/FOXM1 pathway independently of REV-ERBs.

Deep autoencoder for interpretable tissue-adaptive deconvolution and cell-type-specific gene analysis

Single-cell RNA-sequencing has become a powerful tool to study biologically significant characteristics at explicitly high resolution. However, its application on emerging data is currently limited by its intrinsic techniques. Here, we introduce Tissue-AdaPtive autoEncoder (TAPE), a deep learning method connecting bulk RNA-seq and single-cell RNA-seq to achieve precise deconvolution in a short time. By constructing an interpretable decoder and training under a unique scheme, TAPE can predict cell-type fractions and cell-type-specific gene expression tissue-adaptively. Compared with popular methods on several datasets, TAPE has a better overall performance and comparable accuracy at cell type level. Additionally, it is more robust among different cell types, faster, and sensitive to provide biologically meaningful predictions. Moreover, through the analysis of clinical data, TAPE shows its ability to predict cell-type-specific gene expression profiles with biological significance. We believe that TAPE will enable and accelerate the precise analysis of high-throughput clinical data in a wide range.

PRR11 promotes cell proliferation by regulating PTTG1 through interacting with E2F1 transcription factor in pan-cancer

The upregulated proline rich 11 (PRR11) plays a critical role in cancer progression. The relevant biological functions of PRR11 in pan-cancer development are not well understood. In the current study, we found that PRR11 was upregulated in 19 cancer types compared with that of normal tissues and high-expressed PRR11 was a predictor of poor prognosis in 10 cancer types by bioinformatics. Then we showed that interfering PRR11 on three cancer cell lines could greatly inhibit cell proliferation and migration and arrest cells to S phase in vivo. Based on RNA-seq, downregulation of PRR11 expression could extremely suppress the expression of PTTG1 and the cell cycle pathway identified by a differentially expressed gene analysis and an enrichment analysis. The expression of PRR11 and PTTG1 was positively correlated in TCGA and independent GEO data sets. Importantly, we revealed that the PRR11 could express itself in the nucleus and interact with E2F1 on the PTTG1 promoter region to increase the expression of PTTG1. Further results indicated that the expression of PTTG1 was also associated with poor prognosis in 10 cancer types, while downregulation of PTTG1 expression could inhibit cancer cell proliferation and migration. Therefore, we found that PRR11 served as an oncogene in pan-cancer and could influence the cell cycle progression through regulating the expression of PTTG1 by interacting with the transcription factor E2F1.

The Prognostic Assessment of CDC20 in Patients with Renal Clear Cell Carcinoma and Its Relationship with Body Immunity

This article analyzes the relationship between cell division cycle (CDC20) molecules and oncology outcomes in patients with renal clear cell carcinoma (KIRC). CDC20 appears to act as a regulatory protein interacting with many other proteins at multiple points in the cycle. The RNA sequencing data and corresponding clinical information of CDC20 molecules were obtained from The Cancer Genome Atlas (TCGA) database. The expression of CDC20 in kidney renal clear cell carcinoma tissue and adjacent normal tissue was detected by immunohistochemical methods. Logistic analysis was performed to analyze the role of CDC20 in the clinicopathological characteristics and prognosis of KIRC. Gene Set Enrichment Analysis (GSEA) was used to identify the signal pathways which were related to CDC20. Independent prognostic factors were evaluated using univariate and multivariate Cox regression analysis. A nomogram involved in CDC20 expression and clinicopathological variables was conducted to predict overall survival (OS) in KIRC patients at 1, 3, and 5 years. Furthermore, the relation between CDC20 and immunity was also studied. Our results showed that CDC20 was upregulated in kidney renal clear cell carcinoma tissues, accompanying shorter OS (all P < 0.05). According to the results obtained by immunohistochemistry and TCGA database, CDC20 was significantly upregulated in kidney renal clear cell carcinoma tissues compared with neighboring normal kidney tissues. Univariate and multivariate Cox regression analysis showed that high expression of CDC20 was an independent prognostic factor of poor prognosis in kidney renal clear cell carcinoma patients (all P < 0.05). GSEA analysis suggested that the high expression of CDC20 was related to eight multiple signaling pathways. In addition, CDC20 was linked to tumour mutation burden (TMB), immune checkpoint molecules, tumour microenvironment, and immunological infiltration.

NEAT1/MALAT1/XIST/PKD--Hsa-Mir-101-3p--DLGAP5 Axis as a Novel Diagnostic and Prognostic Biomarker Associated With Immune Cell Infiltration in Bladder Cancer

Background: The clinical value of the biomarkers of bladder cancer (BC) is limited due to their low sensitivity or specificity. As a biomarker, DLG associated protein 5 (DLGAP5) is a potential cell cycle regulator in cancer cell carcinogenesis. However, its functional part in BC remains unclear. Therefore, this study aims to identify DLGAP5 expression in BC and its potential diagnostic and prognostic values. Eventually, it predicts the possible RNA regulatory pathways of BC.

Methods: Data on DLGAP5 expression levels in BC and normal bladder tissues were obtained from The Cancer Genome Atlas (TCGA) and the Gene Expression Omnibus (GEO) databases. The receiver operating characteristic (ROC), Kaplan–Meier survival curves, and the univariate and multivariate Cox regression analysis determined the diagnostic and prognostic values of DLGAP5 in BC patients. Finally, the StarBase predicted the target RNAs and constructed networks using Cytoscape.

Results: DLGAP5 expression was significantly upregulated in BC tissue, verified by the TCGA (p &lt; 0.001), GSE3167, GSE7476, and GSE65635 datasets (p &lt; 0.01). BC patients with increased DLGAP5 had poor overall survival (OS) (p = 0.01), disease specific survival (DSS) (p = 0.006) and progress free interval (DFI) (p = 0.007). The area under the ROC curve (AUC) was 0.913. The multivariate Cox analysis identified that lymphovascular invasion (p = 0.007) and DLGAP5 (p = 0.002) were independent prognostic factors.

Conclusion: Increased DLGAP5 expression was closely associated with a poor prognosis in BC patients. In this case, DLGAP5 might be a diagnostic and prognostic biomarker for BC. DLGAP5 expression might be regulated by NEAT1/MALAT1/XIST/PKD--Hsa-mir-101-3p pathways.

Understanding and targeting prostate cancer cell heterogeneity and plasticity

Prostate cancer (PCa) is a prevalent malignancy that occurs primarily in old males. Prostate tumors in different patients manifest significant inter-patient heterogeneity with respect to histo-morphological presentations and molecular architecture. An individual patient tumor also harbors genetically distinct clones in which PCa cells display intra-tumor heterogeneity in molecular features and phenotypic marker expression. This inherent PCa cell heterogeneity, e.g., in the expression of androgen receptor (AR), constitutes a barrier to the long-term therapeutic efficacy of AR-targeting therapies. Furthermore, tumor progression as well as therapeutic treatments induce PCa cell plasticity such that AR-positive PCa cells may turn into AR-negative cells and prostate tumors may switch lineage identity from adenocarcinomas to neuroendocrine-like tumors. This induced PCa cell plasticity similarly confers resistance to AR-targeting and other therapies. In this review, I first discuss PCa from the perspective of an abnormal organ development and deregulated cellular differentiation, and discuss the luminal progenitor cells as the likely cells of origin for PCa. I then focus on intrinsic PCa cell heterogeneity in treatment-naïve tumors with the presence of prostate cancer stem cells (PCSCs). I further elaborate on PCa cell plasticity induced by genetic alterations and therapeutic interventions, and present potential strategies to therapeutically tackle PCa cell heterogeneity and plasticity. My discussions will make it clear that, to achieve enduring clinical efficacy, both intrinsic PCa cell heterogeneity and induced PCa cell plasticity need to be targeted with novel combinatorial approaches.

Characterization of gene cluster heterogeneity in single-cell transcriptomic data within and across cancer types

Despite the remarkable progress in probing tumor transcriptomic heterogeneity by single-cell RNA sequencing (sc-RNAseq) data, several gaps exist in prior studies. Tumor heterogeneity is frequently mentioned but not quantified. Clustering analyses typically target cells rather than genes, and differential levels of transcriptomic heterogeneity of gene clusters are not characterized. Relations between gene clusters inferred from multiple datasets remain less explored. We provided a series of quantitative methods to analyze cancer sc-RNAseq data. First, we proposed two quantitative measures to assess intra-tumoral heterogeneity/homogeneity. Second, we established a hierarchy of gene clusters from sc-RNAseq data, devised an algorithm to reduce the gene cluster hierarchy to a compact structure, and characterized the gene clusters with functional enrichment and heterogeneity. Third, we developed an algorithm to align the gene cluster hierarchies from multiple datasets to a small number of meta gene clusters. By applying these methods to nine cancer sc-RNAseq datasets, we discovered that cancer cell transcriptomes were more homogeneous within tumors than the accompanying normal cells. Furthermore, many gene clusters from the nine datasets were aligned to two large meta gene clusters, which had high and low heterogeneity and were enriched with distinct functions. Finally, we found the homogeneous meta gene cluster retained stronger expression coherence and associations with survival times in bulk level RNAseq data than the heterogeneous meta gene cluster, yet the combinatorial expression patterns of breast cancer subtypes in bulk level data were not preserved in single-cell data. The inference outcomes derived from nine cancer sc-RNAseq datasets provide insights about the contributing factors for transcriptomic heterogeneity of cancer cells and complex relations between bulk level and single-cell RNAseq data. They demonstrate the utility of our methods to enable a comprehensive characterization of co-expressed gene clusters in a wide range of sc-RNAseq data in cancers and beyond.

Summary: We propose quantitative methods to analyze cancer sc-RNAseq data: measures of intra-tumoral heterogeneity, characterization of a hierarchy of gene clusters, and alignment of gene cluster hierarchies from multiple datasets.

Human liver single nucleus and single cell RNA sequencing identify a hepatocellular carcinoma-associated cell-type affecting survival

BACKGROUND

Hepatocellular carcinoma (HCC) is a common primary liver cancer with poor overall survival. We hypothesized that there are HCC-associated cell-types that impact patient survival.

METHODS

We combined liver single nucleus (snRNA-seq), single cell (scRNA-seq), and bulk RNA-sequencing (RNA-seq) data to search for cell-type differences in HCC. To first identify cell-types in HCC, adjacent non-tumor tissue, and normal liver, we integrated single-cell level data from a healthy liver cohort (n = 9 non-HCC samples) collected in the Strasbourg University Hospital; an HCC cohort (n = 1 non-HCC, n = 14 HCC-tumor, and n = 14 adjacent non-tumor samples) collected in the Singapore General Hospital and National University; and another HCC cohort (n = 3 HCC-tumor and n = 3 adjacent non-tumor samples) collected in the Dumont-UCLA Liver Cancer Center. We then leveraged these single cell level data to decompose the cell-types in liver bulk RNA-seq data from HCC patients' tumor (n = 361) and adjacent non-tumor tissue (n = 49) from the Cancer Genome Atlas (TCGA) multi-center cohort. For replication, we decomposed 221 HCC and 209 adjacent non-tumor liver microarray samples from the Liver Cancer Institute (LCI) cohort collected by the Liver Cancer Institute and Zhongshan Hospital of Fudan University.

RESULTS

We discovered a tumor-associated proliferative cell-type, Prol (80.4% tumor cells), enriched for cell cycle and mitosis genes. In the liver bulk tissue from the TCGA cohort, the proportion of the Prol cell-type is significantly increased in HCC and associates with a worse overall survival. Independently from our decomposition analysis, we reciprocally show that Prol nuclei/cells significantly over-express both tumor-elevated and survival-decreasing genes obtained from the bulk tissue. Our replication analysis in the LCI cohort confirmed that an increased estimated proportion of the Prol cell-type in HCC is a significant marker for a shorter overall survival. Finally, we show that somatic mutations in the tumor suppressor genes TP53 and RB1 are linked to an increase of the Prol cell-type in HCC.

CONCLUSIONS

By integrating liver single cell, single nucleus, and bulk expression data from multiple cohorts we identified a proliferating cell-type (Prol) enriched in HCC tumors, associated with a decreased overall survival, and linked to TP53 and RB1 somatic mutations.

Construction of Bone Metastasis-Specific Regulation Network Based on Prognostic Stemness-Related Signatures in Prostate Cancer

Background

We planned to uncover the cancer stemness-related genes (SRGs) in prostate cancer (PCa) and its underlying mechanism in PCa metastasis.

Methods

We acquired the RNA-seq data of 406 patients with PCa from the TCGA database. Based on the mRNA stemness index (mRNAsi) calculated by one-class logistic regression (OCLR) algorithm, SRGs in PCa were extracted by WGCNA. Univariate and multivariate regression analyses were applied to uncover OS-associated SRGs. Gene Set Variation Analysis (GSVA), Gene Set Enrichment Analysis (GSEA), and Pearson's correlation analysis were performed to discover the possible mechanism of PCa metastasis. The significantly correlated transcription factors of OS-associated SRGs were also identified by Pearson's correlation analysis. ChIP-seq was applied to validate the binding relationship of TFs and OS-associated SRGs and spatial transcriptome and single-cell sequencing were performed to uncover the location of key biomarkers expression. Lastly, we explored the specific inhibitors for SRGs using CMap algorithm.

Results

We identified 538 differentially expressed genes (DEGs) between non-metastatic and metastatic PCa. Furthermore, OS-associated SRGs were identified. The Pearson correlation analysis revealed that FOXM1 was significantly correlated with NEIL3 (correlation efficient =0.89, p < 0.001) and identified hallmark_E2F_targets as the potential pathway mechanism of NEIL3 promoting PCa metastasis (correlation efficient =0.58, p < 0.001). Single-cell sequencing results indicated that FOXM1 regulating NEIL3 may get involved in the antiandrogen resistance of PCa. Rottlerin was discovered to be a potential target drug for PCa.

Conclusion

We constructed a regulatory network based on SRGs associated with PCa metastasis and explored possible mechanism.

miR‑374c‑5p regulates PTTG1 and inhibits cell growth and metastasis in hepatocellular carcinoma by regulating epithelial‑mesenchymal transition

A large number of studies have reported that microRNA (miR)-374c-5p plays an important role in the occurrence and development of malignant tumors, but there is no research on the role of miR-374c-5p in hepatocellular carcinoma (HCC). The aim of the present study was to investigate the role of miR-374c-5p in HCC and the underlying molecular mechanism. The expression of miR-374c-5p in HCC tissues and HCC cell lines was analyzed via reverse transcription-quantitative PCR. The association between miR-374c-5p and clinical pathology was also analyzed in patients with HCC. Kaplan-Meier analysis and Cox multivariate analysis were used to evaluate the prognostic significance of miR-374c-5p in HCC. The biological functions of miR-374c-5p, including cell proliferation, migration and invasion and its potential molecular mechanism were analyzed in vivo and in vitro. In addition, the molecular mechanism of miR-374c-5p in HCC was further explored. The results demonstrated that miR-374c-5p expression was lower in HCC than in matched adjacent tissue samples. Patients with low expression of miR-374c-5p had poor prognosis and short survival time. Overexpression of miR-374c-5p inhibited HCC cell proliferation, migration and invasion in vitro. In vivo, it was found that overexpression of miR-374c-5p significantly inhibited the growth and proliferation of HCC cells. Dual-luciferase reporter assays verified that miR-374c-5p directly targets the 3′-untranslated region of pituitary tumor-transforming 1 (PTTG1) and regulates PTTG1 expression. In general, it was revealed that miR-374c-5p regulates the malignant biological behavior of HCC through PTTG1, thereby affecting epithelial-mesenchymal transition. Thus, miR-374c-5p is a potential biological indicator to predict poor prognosis in patients with HCC.

Single-cell RNA sequencing technologies and applications: A brief overview

Single-cell RNA sequencing (scRNA-seq) technology has become the state-of-the-art approach for unravelling the heterogeneity and complexity of RNA transcripts within individual cells, as well as revealing the composition of different cell types and functions within highly organized tissues/organs/organisms. Since its first discovery in 2009, studies based on scRNA-seq provide massive information across different fields making exciting new discoveries in better understanding the composition and interaction of cells within humans, model animals and plants. In this review, we provide a concise overview about the scRNA-seq technology, experimental and computational procedures for transforming the biological and molecular processes into computational and statistical data. We also provide an explanation of the key technological steps in implementing the technology. We highlight a few examples on how scRNA-seq can provide unique information for better understanding health and diseases. One important application of the scRNA-seq technology is to build a better and high-resolution catalogue of cells in all living organism, commonly known as atlas, which is key resource to better understand and provide a solution in treating diseases. While great promises have been demonstrated with the technology in all areas, we further highlight a few remaining challenges to be overcome and its great potentials in transforming current protocols in disease diagnosis and treatment.

The Therapeutic Effect of iMSC-Derived Small Extracellular Vesicles on Tendinopathy Related Pain Through Alleviating Inflammation: An in vivo and in vitro Study

Background

Tendinopathy is a common cause of tendon pain. However, there is a lack of effective therapies for managing tendinopathy pain, despite the pain being the most common complaint of patients. This study aimed to evaluate the therapeutic effect of small extracellular vesicles released from induced pluripotent stem cell-derived mesenchymal stem cells (iMSC-sEVs) on tendinopathy pain and explore the underlying mechanisms.

Methods

Rat tendinopathy model was established and underwent the injection of iMSC-sEVs to the quadriceps tendon one week after modeling. Pain-related behaviors were measured for the following four weeks. Tendon histology was assessed four weeks after the injection. To further investigate the potential mechanism, tenocytes were stimulated with IL-1β to mimic tendinopathy in vitro. The effect of iMSC-sEVs on tenocyte proliferation and the expression of proinflammatory cytokines were measured by CCK-8, RT-qPCR, and ELISA. RNA-seq was further performed to systematically analyze the related global changes and underlying mechanisms.

Results

Local injection of iMSC-sEVs was effective in alleviating pain in the tendinopathy rats compared with the vehicle group. Tendon histology showed ameliorated tendinopathy characteristics. Upon iMSC-sEVs treatment, significantly increased tenocyte proliferation and less expression of proinflammatory cytokines were observed. Transcriptome analysis revealed that iMSC-sEVs treatment upregulated the expression of genes involved in cell proliferation and downregulated the expression of genes involved in inflammation and collagen degeneration.

Conclusion

Collectively, this study demonstrated iMSC-sEVs protect tenocytes from inflammatory stimulation and promote cell proliferation as well as collagen synthesis, thereby relieving pain derived from tendinopathy. As a cell-free regenerative treatment, iMSC-sEVs might be a promising therapeutic candidate for tendinopathy.

Single-cell gene fusion detection by scFusion

Gene fusions can play important roles in tumor initiation and progression. While fusion detection so far has been from bulk samples, full-length single-cell RNA sequencing (scRNA-seq) offers the possibility of detecting gene fusions at the single-cell level. However, scRNA-seq data have a high noise level and contain various technical artifacts that can lead to spurious fusion discoveries. Here, we present a computational tool, scFusion, for gene fusion detection based on scRNA-seq. We evaluate the performance of scFusion using simulated and five real scRNA-seq datasets and find that scFusion can efficiently and sensitively detect fusions with a low false discovery rate. In a T cell dataset, scFusion detects the invariant TCR gene recombinations in mucosal-associated invariant T cells that many methods developed for bulk data fail to detect; in a multiple myeloma dataset, scFusion detects the known recurrent fusion IgH-WHSC1, which is associated with overexpression of the WHSC1 oncogene. Our results demonstrate that scFusion can be used to investigate cellular heterogeneity of gene fusions and their transcriptional impact at the single-cell level.

Androgen receptor negatively regulates mitotic checkpoint signaling to induce docetaxel resistance in castration-resistant prostate cancer

BACKGROUND

Despite multiple treatment advances for castration-resistant prostate cancer (CRPC), there are currently no curative therapies and patients ultimately to succumb to the disease. Docetaxel (DTX) is the standard first-line chemotherapy for patients with metastatic CRPC; however, drug resistance is inevitable and often develops rapidly, leading to disease progression in nearly all patients. In contrast, when DTX is deployed with androgen deprivation therapy in castration-sensitive disease, more durable responses and improved outcomes are observed, suggesting that aberrant androgen receptor (AR) signaling accelerates DTX resistance in CRPC. In this study, we demonstrate that AR dysregulates the mitotic checkpoint, a critical pathway involved in the anticancer action of DTX.

METHODS

Androgen-dependent and independent cell lines were used to evaluate the role of AR in DTX resistance. Impact of drug treatment on cell viability, survival, and cell-cycle distribution were determined by plate-based viability assay, clonogenic assay, and cell-cycle analysis by flow cytometry, respectively. Mitotic checkpoint kinase signal transduction and apoptosis activation was evaluated by Western blotting. Pathway gene expression analysis was evaluated by RT-PCR. A Bliss independence model was used to calculate synergy scores for drug combination studies.

RESULTS

Activation of AR in hormone-sensitive cells induces a rescue phenotype by increasing cell viability and survival and attenuating G2/M arrest in response to DTX. Analysis of mitotic checkpoint signaling shows that AR negatively regulates spindle checkpoint signaling, resulting in premature mitotic progression and evasion of apoptosis. This phenotype is characteristic of mitotic slippage and is also observed in CRPC cell lines where we demonstrate involvement of AR splice variant AR-v7 in dysregulation of checkpoint signaling. Our findings suggest that DTX resistance is mediated through mechanisms that drive premature mitotic exit. Using pharmacologic inhibitors of anaphase-promoting complex/cyclosome and polo-like kinase 1, we show that blocking mitotic exit induces mitotic arrest, apoptosis, and synergistically inhibits cell survival in combination with DTX.

CONCLUSION

Our results suggest that targeting the mechanisms of dysregulated mitotic checkpoint signaling in AR-reactivated tumors has significant clinical potential to extend treatment benefit with DTX and improve outcomes in patients with lethal prostate 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.

Capturing the third dimension in drug discovery: Spatially-resolved tools for interrogation of complex 3D cell models

Advanced three-dimensional (3D) cell models have proven to be capable of depicting architectural and microenvironmental features of several tissues. By providing data of higher physiological and pathophysiological relevance, 3D cell models have been contributing to a better understanding of human development, pathology onset and progression mechanisms, as well as for 3D cell-based assays for drug discovery. Nonetheless, the characterization and interrogation of these tissue-like structures pose major challenges on the conventional analytical methods, pushing the development of spatially-resolved technologies. Herein, we review recent advances and pioneering technologies suitable for the interrogation of multicellular 3D models, while capable of retaining biological spatial information. We focused on imaging technologies and omics tools, namely transcriptomics, proteomics and metabolomics. The advantages and shortcomings of these novel methodologies are discussed, alongside the opportunities to intertwine data from the different tools.

Alternative RNA splicing in stem cells and cancer stem cells: Importance of transcript-based expression analysis

Alternative ribonucleic acid (RNA) splicing can lead to the assembly of different protein isoforms with distinctive functions. The outcome of alternative splicing (AS) can result in a complete loss of function or the acquisition of new functions. There is a gap in knowledge of abnormal RNA splice variants promoting cancer stem cells (CSCs), and their prospective contribution in cancer progression. AS directly regulates the self-renewal features of stem cells (SCs) and stem-like cancer cells. Notably, octamer-binding transcription factor 4A spliced variant of octamer-binding transcription factor 4 contributes to maintaining stemness properties in both SCs and CSCs. The epithelial to mesenchymal transition pathway regulates the AS events in CSCs to maintain stemness. The alternative spliced variants of CSCs markers, including cluster of differentiation 44, aldehyde dehydrogenase, and doublecortin-like kinase, α6β1 integrin, have pivotal roles in increasing self-renewal properties and maintaining the pluripotency of CSCs. Various splicing analysis tools are considered in this study. LeafCutter software can be considered as the best tool for differential splicing analysis and identification of the type of splicing events. Additionally, LeafCutter can be used for efficient mapping splicing quantitative trait loci. Altogether, the accumulating evidence re-enforces the fact that gene and protein expression need to be investigated in parallel with alternative splice variants.

A New Stemness-Related Prognostic Model for Predicting the Prognosis in Pancreatic Ductal Adenocarcinoma

Objective

This study is aimed at identifying stemness-related genes in pancreatic ductal adenocarcinoma (PDAC).

Methods

The RNA-seq data of PADC patients were downloaded from The Cancer Genome Atlas (TCGA) and Genotype-Tissue Expression (GTEx) databases. The mRNA expression-based stemness index (mRNAsi) and epigenetically regulated mRNAsi (EREG-mRNAsi) of PADC patients were evaluated. The mRNAsi-related gene sets in PADC were identified by weighted gene coexpression network analysis (WGCNA). The key genes were further analyzed using functional enrichment analysis. The Kaplan-Meier survival analysis and the Cox proportional hazards model were used to evaluate the prognostic value of the key genes. Prognostic hub genes were used to establish nomograms. The receiver operating characteristic (ROC) curves, concordance index (C-index), and calibration curves were used to assess the discrimination and accuracy of the nomogram. Finally, these results were validated in the Gene Expression Omnibus (GEO) database.

Results

A total of 36 key genes related to mRNAsi were identified by WGCNA. A prognostic gene signature compromising seven genes (TPX2, ZWINT, UBE2C, CCNB2, CDK1, BUB1, and BIRC5) was established to predict the overall survival (OS) of PADC patients. The Cox regression analysis revealed that the risk score was an independent prognostic factor for PADC. Patients were then divided into the high-risk and low-risk groups. The ROC curves, C-index, and calibration curves indicated good performance of the prognostic signature in the TCGA and GEO datasets. Moreover, the nomogram incorporating clinical parameters showed better sensitivity and specificity for predicting the OS of PADC patients.

Conclusion

The stemness-related prognostic model successfully predicted the OS of PADC patients and could be used for the treatment of PADC.

Screening and predicted value of potential biomarkers for breast cancer using bioinformatics analysis

Breast cancer is the most common cancer and the leading cause of cancer-related deaths in women. Increasing molecular targets have been discovered for breast cancer prognosis and therapy. However, there is still an urgent need to identify new biomarkers. Therefore, we evaluated biomarkers that may aid the diagnosis and treatment of breast cancer. We searched three mRNA microarray datasets (GSE134359, GSE31448 and GSE42568) and identified differentially expressed genes (DEGs) by comparing tumor and non-tumor tissues using GEO2R. Functional and pathway enrichment analyses of the DEGs were performed using the DAVID database. The protein-protein interaction (PPI) network was plotted with STRING and visualized using Cytoscape. Module analysis of the PPI network was done using MCODE. The associations between the identified genes and overall survival (OS) were analyzed using an online Kaplan-Meier tool. The redundancy analysis was conducted by DepMap. Finally, we verified the screened HUB gene at the protein level. A total of 268 DEGs were identified, which were mostly enriched in cell division, cell proliferation, and signal transduction. The PPI network comprised 236 nodes and 2132 edges. Two significant modules were identified in the PPI network. Elevated expression of the genes Discs large-associated protein 5 (DLGAP5), aurora kinase A (AURKA), ubiquitin-conjugating enzyme E2 C (UBE2C), ribonucleotide reductase regulatory subunit M2(RRM2), kinesin family member 23(KIF23), kinesin family member 11(KIF11), non-structural maintenance of chromosome condensin 1 complex subunit G (NCAPG), ZW10 interactor (ZWINT), and denticleless E3 ubiquitin protein ligase homolog(DTL) are associated with poor OS of breast cancer patients. The enriched functions and pathways included cell cycle, oocyte meiosis and the p53 signaling pathway. The DEGs in breast cancer have the potential to become useful targets for the diagnosis and treatment of breast cancer.

Single-cell ATAC and RNA sequencing reveal pre-existing and persistent cells associated with prostate cancer relapse

Prostate cancer is heterogeneous and patients would benefit from methods that stratify those who are likely to respond to systemic therapy. Here, we employ single-cell assays for transposase-accessible chromatin (ATAC) and RNA sequencing in models of early treatment response and resistance to enzalutamide. In doing so, we identify pre-existing and treatment-persistent cell subpopulations that possess regenerative potential when subjected to treatment. We find distinct chromatin landscapes associated with enzalutamide treatment and resistance that are linked to alternative transcriptional programs. Transcriptional profiles characteristic of persistent cells are able to stratify the treatment response of patients. Ultimately, we show that defining changes in chromatin and gene expression in single-cell populations from pre-clinical models can reveal as yet unrecognized molecular predictors of treatment response. This suggests that the application of single-cell methods with high analytical resolution in pre-clinical models may powerfully inform clinical decision-making.

Cancer biology deciphered by single-cell transcriptomic sequencing

Tumors are complex ecosystems in which heterogeneous cancer cells interact with their microenvironment composed of diverse immune, endothelial, and stromal cells. Cancer biology had been studied using bulk genomic and gene expression profiling, which however mask the cellular diversity and average the variability among individual molecular programs. Recent advances in single-cell transcriptomic sequencing have enabled a detailed dissection of tumor ecosystems and promoted our understanding of tumorigenesis at single-cell resolution. In the present review, we discuss the main topics of recent cancer studies that have implemented single-cell RNA sequencing (scRNA-seq). To study cancer cells, scRNA-seq has provided novel insights into the cancer stem-cell model, treatment resistance, and cancer metastasis. To study the tumor microenvironment, scRNA-seq has portrayed the diverse cell types and complex cellular states of both immune and non-immune cells interacting with cancer cells, with the promise to discover novel targets for future immunotherapy.

Regulation of cancer stem cells in triple negative breast cancer

Triple Negative Breast Cancer (TNBC) is the most lethal subtype of breast cancer. Despite the successes of emerging targeted therapies, relapse, recurrence, and therapy failure rates in TNBC significantly outpace other subtypes of breast cancer. Mounting evidence suggests accumulation of therapy resistant Cancer Stem Cell (CSC) populations within TNBCs contributes to poor clinical outcomes. These CSCs are enriched in TNBC compared to non-TNBC breast cancers. The mechanisms underlying CSC accumulation have been well-characterized and discussed in other reviews. In this review, we focus on TNBC-specific mechanisms that allow the expansion and activity of self-renewing CSCs. We highlight cellular signaling pathways and transcription factors, specifically enriched in TNBC over non-TNBC breast cancer, contributing to stemness. We also analyze publicly available single-cell RNA-seq data from basal breast cancer tumors to highlight the potential of emerging bioinformatic approaches in identifying novel drivers of stemness in TNBC and other cancers.

CDC20 and PTTG1 are Important Biomarkers and Potential Therapeutic Targets for Metastatic Prostate Cancer

INTRODUCTION

Metastatic prostate cancer (mPCa) is responsible for most prostate cancer (PCa) deaths worldwide. The present study aims to explore the molecular differences between mPCa and PCa.

METHODS

The authors downloaded GSE6752, GSE6919, and GSE32269 from the Gene Expression Omnibus and employed integrated analysis to identify differentially expressed genes (DEGs) between mPCa and PCa. Functional and pathway-enrichment analyses were performed, and a protein-protein interaction (PPI) network and modules were constructed. Clinical mPCa specimens were collected to verify the results by performing RT-qPCR. The Cancer Genome Atlas database was used to conduct a survival analysis, and an immunohistochemical assay was performed. The invasion ability of PCa cells was verified by Transwell assay.

RESULTS

One-hundred six consistently DEGs were found in mPCa compared with PCa. DEGs significantly enriched the positive regulation of cell proliferation, cell division, and cell adhesion in small cell lung cancer and PCa. Cell division, nucleoplasm, and cell cycle were selected from the PPI network, and the top 10 hub genes were selected. CDC20 and PTTG1 with genetic alterations were significantly associated with poorer disease-free survival. Immunohistochemical assay results showed that the expression levels of CDC20 and PTTG1 in mPCa were higher than those in PCa. The results of the migration assay indicated that CDC20 and PTTG1 could enhance the migration ability of PCa cells.

CONCLUSION

The present study revealed that CDC20 and PTTG1 contribute more to migration, progression, and poorer prognoses in mPCa compared with PCa. CDC20 and PTTG1 could represent therapeutic targets in mPCa medical research and clinical studies.

High copy number variations, particular transcription factors, and low immunity contribute to the stemness of prostate cancer cells

Background

Tumor metastasis is the main cause of death of cancer patients, and cancer stem cells (CSCs) is the basis of tumor metastasis. However, systematic analysis of the stemness of prostate cancer cells is still not abundant. In this study, we explore the effective factors related to the stemness of prostate cancer cells by comprehensively mining the multi-omics data from TCGA database.

Methods

Based on the prostate cancer transcriptome data in TCGA, gene expression modules that strongly relate to the stemness of prostate cancer cells are obtained with WGCNA and stemness scores. Copy number variation of stemness genes of prostate cancer is calculated and the difference of transcription factors between prostate cancer and normal tissues is evaluated by using CNV (copy number variation) data and ATAC-seq data. The protein interaction network of stemness genes in prostate cancer is constructed using the STRING database. Meanwhile, the correlation between stemness genes of prostate cancer and immune cells is analyzed.

Results

Prostate cancer with higher Gleason grade possesses higher cell stemness. The gene set highly related to prostate cancer stemness has higher CNV in prostate cancer samples than that in normal samples. Although the transcription factors of stemness genes have similar expressions, they have different contributions between normal and prostate cancer tissues; and particular transcription factors enhance the stemness of prostate cancer, such as PUM1, CLOCK, SP1, TCF12, and so on. In addition, the lower tumor immune microenvironment is conducive to the stemness of prostate cancer. CD8 + T cells and M1 macrophages may play more important role in the stemness of prostate cancer than other immune cells in the tumor microenvironment. Finally, EZH2 is found to play a central role in stemness genes and is negatively correlated with resting mast cells and positively correlated with activated memory CD4 + T cells.

Conclusions

Based on the systematic and combined analysis of multi-omics data, we find that high copy number variation, specific transcription factors, and low immune microenvironment jointly contribute to the stemness of prostate cancer cells. These findings may provide us new clues and directions for the future research on stemness of prostate cancer.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12967-021-02870-x.

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.

The current landscape of single-cell transcriptomics for cancer immunotherapy

Immunotherapies such as immune checkpoint blockade and adoptive cell transfer have revolutionized cancer treatment, but further progress is hindered by our limited understanding of tumor resistance mechanisms. Emerging technologies now enable the study of tumors at the single-cell level, providing unprecedented high-resolution insights into the genetic makeup of the tumor microenvironment and immune system that bulk genomics cannot fully capture. Here, we highlight the recent key findings of the use of single-cell RNA sequencing to deconvolute heterogeneous tumors and immune populations during immunotherapy. Single-cell RNA sequencing has identified new crucial factors and cellular subpopulations that either promote tumor progression or leave tumors vulnerable to immunotherapy. We anticipate that the strategic use of single-cell analytics will promote the development of the next generation of successful, rationally designed immunotherapeutics.

Toward uncharted territory of cellular heterogeneity: advances and applications of single-cell RNA-seq

Among single-cell analysis technologies, single-cell RNA-seq (scRNA-seq) has been one of the front runners in technical inventions. Since its induction, scRNA-seq has been well received and undergone many fast-paced technical improvements in cDNA synthesis and amplification, processing and alignment of next generation sequencing reads, differentially expressed gene calling, cell clustering, subpopulation identification, and developmental trajectory prediction. scRNA-seq has been exponentially applied to study global transcriptional profiles in all cell types in humans and animal models, healthy or with diseases, including cancer. Accumulative novel subtypes and rare subpopulations have been discovered as potential underlying mechanisms of stochasticity, differentiation, proliferation, tumorigenesis, and aging. scRNA-seq has gradually revealed the uncharted territory of cellular heterogeneity in transcriptomes and developed novel therapeutic approaches for biomedical applications. This review of the advancement of scRNA-seq methods provides an exploratory guide of the quickly evolving technical landscape and insights of focused features and strengths in each prominent area of progress.

I-Impute: a self-consistent method to impute single cell RNA sequencing data

Background

Single-cell RNA-sequencing (scRNA-seq) is becoming indispensable in the study of cell-specific transcriptomes. However, in scRNA-seq techniques, only a small fraction of the genes are captured due to “dropout” events. These dropout events require intensive treatment when analyzing scRNA-seq data. For example, imputation tools have been proposed to estimate dropout events and de-noise data. The performance of these imputation tools are often evaluated, or fine-tuned, using various clustering criteria based on ground-truth cell subgroup labels. This limits their effectiveness in the cases where we lack cell subgroup knowledge. We consider an alternative strategy which requires the imputation to follow a “self-consistency” principle; that is, the imputation process is to refine its results until there is no internal inconsistency or dropouts from the data.

Results

We propose the use of “self-consistency” as a main criteria in performing imputation. To demonstrate this principle we devised I-Impute, a “self-consistent” method, to impute scRNA-seq data. I-Impute optimizes continuous similarities and dropout probabilities, in iterative refinements until a self-consistent imputation is reached. On the in silico data sets, I-Impute exhibited the highest Pearson correlations for different dropout rates consistently compared with the state-of-art methods SAVER and scImpute. Furthermore, we collected three wetlab datasets, mouse bladder cells dataset, embryonic stem cells dataset, and aortic leukocyte cells dataset, to evaluate the tools. I-Impute exhibited feasible cell subpopulation discovery efficacy on all the three datasets. It achieves the highest clustering accuracy compared with SAVER and scImpute.

Conclusions

A strategy based on “self-consistency”, captured through our method, I-Impute, gave imputation results better than the state-of-the-art tools. Source code of I-Impute can be accessed at https://github.com/xikanfeng2/I-Impute.