Tianqi Lyu
Yuanbin Lin
Kerong Wu2- 1Cixi Institute of Biomedical Engineering, Ningbo Institute of Materials Technology and Engineering, Chinese Academy of Science (CAS), Ningbo, China
- 2Department of Urology, Ningbo First Hospital, School of Medicine Ningbo University, Zhejiang University Ningbo Hospital, Ningbo, China
Bladder cancer is among the most common malignant tumors with highly heterogeneous molecular characteristics. Despite advancements of the available therapeutic options, several bladder cancer patients exhibit unsatisfactory clinical outcomes. The lack of specific biomarkers for effective targeted therapy or immunotherapy remains a major obstacle in treating bladder cancer. The rapid development of single-cell techniques is transforming our understanding of the intra-tumoral heterogeneity, thereby providing us with a powerful high-throughput sequencing tool that can reveal tumorigenesis, progression, and invasion in bladder tumors. In this review, we summarise and discuss how single-cell sequencing technologies have been applied in bladder cancer research, to advance our collective knowledge on the heterogeneity of bladder tumor cells, as well as to provide new insights into the complex ecosystem of the tumor microenvironment. The application of single-cell approaches also uncovers the therapeutic resistance mechanism in bladder cancer and facilitates the detection of urinary-exfoliated tumor cells. Moreover, benefiting from the powerful technical advantages of single-cell techniques, several key therapeutic targets and prognostic models of bladder cancer have been identified. It is hoped that this paper can provide novel insights into the precision medicine of bladder cancer.
Introduction
Bladder cancer is the most common malignant tumor that affects the urinary system and is associated with high morbidity and mortality (Antoni et al., 2017). Several risk factors such as sex, age, lifestyle, pathogenic microbial infections, and drugs may augment the bladder cancer risk (Lin et al., 2006; Letašiová et al., 2012; Burger et al., 2013; Berdik, 2017). Despite advancements of the available therapeutic approaches, the clinical outcomes of bladder cancer patients remain unsatisfactory, and the corresponding 5-year survival rate remains almost unchanged. Bladder cancer is one of the highly heterogeneous human cancers, often involving frequent genomic alterations and molecular subtype diversity, and it may be the first cause of therapeutic failure; therefore, the elucidation of the molecular and cellular heterogeneity of bladder cancer is necessary (Genitsch et al., 2019). Probably, an in-depth understanding of the unique molecular characteristics and cellular heterogeneity involved in the initiation, development, and pathogenesis of bladder cancer can improve the clinical treatment and prognosis, as well as contribute to the promotion of its diagnosis and treatment towards precision medicine.
Over the past few decades, with the rapid development of next-generation sequencing and molecular biology technologies, research on the biological characteristics and genotyping of bladder cancer has made considerable strides (Pietzak et al., 2017; Roy et al., 2017; Wu et al., 2019; Li et al., 2021a; Tran et al., 2021). Bladder cancer can be categorised into three subtypes based on its histological features: urothelial carcinoma, squamous cell carcinoma, and adenocarcinoma. Of these subtypes, urothelial tumors account for approximately 90% of cases, while squamous and glandular-type tumors are rare (Kantor et al., 1988; Knowles and Hurst, 2015). The 2016 World Health Organisation grading system has stratified bladder cancer into non-muscle-invasive bladder cancer (NMIBC) and muscle-invasive bladder cancer (MIBC) (Compérat et al., 2019). At the first clinical diagnosis, approximately 85% of bladder cancer patients are diagnosed with NMIBC, and up to 80% of these patients have at least one recurrence, whereas disease progression into MIBC occurs in approximately 30% of these patients (Batista et al., 2020; van den Bosch and Alfred Witjes, 2011). It is difficult for doctors to predict the risk of recurrence and progression of NMIBC based on the conventional classification system.
As of today, enormous efforts are being undertaken to outline the diverse molecular subtypes and collect crucial biological information on bladder cancer. Several researchers have proposed that the molecular and genetic heterogeneities of bladder cancer cells are critical to diagnosis and treatment of the disease (Thomsen et al., 2017; Audenet et al., 2018; Sjödahl et al., 2019; Minoli et al., 2020). However, despite these advances, clinical therapies based on the molecular subtypes and genetic features of bladder cancer often show several limitations such as poor efficacy, drug resistance, and recurrence due to tumor heterogeneity (Casadei et al., 2019; Kamoun et al., 2019; Zhu et al., 2020). In general, the molecular mechanisms for the initiation and development of bladder cancer remain elusive, and the current molecular subtyping systems do not depict the cellular heterogeneity of bladder carcinoma. In the era of precision medicine, we should consider the vital role played by inter-tumoral and intra-tumoral heterogeneities. Decoding the heterogeneity of bladder cancer cells at the genetic and cellular levels, refinement of the molecular types, and discovering new ways to overcome resistance are the main issues that should be addressed in the studies on bladder cancer.
In 2009, Tang et al. reported the single-cell RNA sequencing (scRNA-seq) technology for the first time, which helped solve the problem of cell heterogeneity that was difficult to be solved by bulk RNA sequencing (Tang et al., 2009). The rapid developments in single-cell techniques allow us to interrogate the transcriptional, genomic, proteomic, epigenomic, metabolic, and multi-omic characteristics of thousands of individual cells (Kashima et al., 2020) (Table 1). This powerful high-resolution and high-throughput sequencing tool has revolutionised our understanding of cancer biology and tumor cytology (Stuart and Satija, 2019; Xing et al., 2020). These technologies have helped in successfully deciphering the compositional complexity and clonal heterogeneity in tumors (Lei et al., 2021). Indeed, single-cell techniques have multiple applications in cancer research, which include the identification of dynamic gene expression profiles during tumor progression, discovery of novel cell subpopulations, determination of cell status and phenotype switches, analysis of the regulatory pathways of key genes, and identification of potential therapeutic targets (Papalexi and Satija, 2018) (Figure 1). Owing to the excellent technical advantages, single-cell techniques have undoubtedly promoted the precision and accuracy of molecular cancer research to new levels. The past decade has seen unprecedented advances in bladder cancer studies, which can be ascribed mainly to single-cell techniques (Table 2). In this review, we summarise pioneering research in the field and discuss the latest progress in the aspects such as cellular heterogeneity, tumor microenvironment (TME), drug resistance, and prognosis of bladder cancer. All of these further inspire and foster research on clinically relevant diagnosis, prognosis prediction, targeted therapy, and personalised therapeutic strategies for bladder cancer.
TABLE 1. Single-cell omics technologies.
FIGURE 1. Schematic diagram of the single-cell sequencing workflow. Cells are dissolved from bladder tumors into single-cell suspension. Single cells are then processed into the reaction system followed by reverse transcription, amplification, library construction and sequencing. Sequencing data are processed via bioinformatic analysis such as cell clustering, cell trajectory, different gene expression, interactome or key marker genes analysis.
TABLE 2. Summary of the studies in bladder cancer using single-cell technologies.
Uncovering the heterogeneity of bladder urothelial cells through single-cell sequencing
The urothelium of the urinary bladder is a particular hierarchically organised tissue that serves as the strongest urine–blood barrier in the body. The adult urothelium can be histologically classified as the transitional epithelium consisting of three types of cells, namely basal, intermediate, and umbrella cells (Figure 2). Among these cells, the basal cells are small and polygonal and form a single layer that directly contacts the basement membrane. Intermediate cells are pyriform and form multiple cell layers. Umbrella cells are large binuclear or multinuclear cubic cells that form a single layer in direct contact with the urinary space (Apodaca, 2004; Khandelwal et al., 2009; Ho et al., 2012). Usually, the bladder urothelium regenerates rapidly upon injury, which depends on the proliferation of urothelial stem cells (Ho et al., 2012). Accumulating evidence suggests that minor subpopulations of basal cells, characterised by the expressions of keratin 5 (KRT5) and keratin 14 (KRT14), possess self-renewal capacity and give rise to all cell types of the urothelium during natural and injury-induced regeneration processes, which represents the cellular origins of urothelial carcinoma (Shin et al., 2011; Papafotiou et al., 2016). Reported research demonstrates that a vast majority of tumor lesions of the bladder arise from the urothelial cells (Castillo-Martin et al., 2010). Therefore, elucidation of the primitive cellular differentiation status, stem cell differentiation, and genetic alteration characteristics of bladder urothelial may provide critical biological insights for understanding the molecular mechanism of associated bladder diseases.
FIGURE 2. Schematic diagram of normal bladder urothelium. Stem cells are shown in blue, basal cells are shown in purple, intermediate cells are shown in brown, umbrella cells are shown in green, and stroma cells are shown in gray.
To gain an insight into the cell subtypes and urothelial differentiation characteristics, Li et al. performed scRNA-seq to detect urothelial cells digested and isolated from the mouse bladder urothelium (Li et al., 2021b). The authors discovered that the ASPM-labelled basal-like cells might be the progenitor cells of the other subpopulations that could participate in regeneration of the bladder urothelium. In particular, a novel superficial-like cell population (Plxna4+ urothelial cells) may play a priming role in the initiation of resisting inflammation and repairing of the responses of the bladder urothelium. Yu et al. compared the single-cell transcriptomic maps of human and mouse bladders to determine both the conservative and heterogeneous aspects of human and mouse bladder evolution (Yu et al., 2019). The authors identified two new types of human bladder cells: ADRA2A+ and HRH2+ interstitial cells and TNNT1+ epithelial cells. Of these, ADRA2A+ interstitial cells appear to play a special role in the human bladder, while TNNT1+ epithelial cells are also present in rat and mouse bladder tissues. Developmental trajectory analysis revealed that the bladder epithelial cells could be transformed from basal cells to intermediate cells and then into umbrella cells or TNNT1+ epithelial cells. Another study recapitulated the differentiation and function of urothelial cells by establishing mouse urothelial organoids, wherein the authors identified five cell clusters by utilising scRNA-seq and reported that the CD49f/CD44 high-labelled urothelial cell subpopulation exhibits the highest organoid-forming potential (Santos et al., 2019; Yu et al., 2019). In previous studies, CD49f and CD44 have been proposed as markers for aggressive basal/squamous bladder cancer subtypes (Lerner et al., 2016; Sjödahl et al., 2017). Notably, the investigation unveiled the essential role of Notch signalling in the differentiation of normal urothelial cells, which was elusive in the past bulk organoid transcriptomics studies. Furthermore, the study confirmed that gene mutations in the Notch pathway could induce urothelial tumors.
Recently, Cheng et al. explored the cell heterogeneity involved in the repair and regeneration processes of bladder urothelial cells and focused on the crosstalk between cell subsets (Cheng et al., 2021). When compared with the normal urothelial cells, two cycling cell populations (i.e., basal cells and luminal cells) were highly enriched under injured conditions. A study demonstrated that urothelium regeneration is mediated by distinct division patterns of the cycling basal and intermediate cells (Wang et al., 2018); in this study, the single-cell analysis confirmed that this discrepancy of cell cycle progression might be attributable to the difference in cell adhesion (Cheng et al., 2021). By analysing the lineage relationship among urothelial cell subpopulations, the authors found that cycling basal cells and cycling intermediate cells represent two different states of the urothelium. Notably, cycling intermediate cells could act as the direct progenitors of superficial cells, where the latter was directly repaired by the former, which verifies the previous hypothesis. According to cell–cell communication analyses, Acta2+Cd34+ myofibroblasts have the most significant influence on the proliferation of urothelial cells, and both Bmp and Fgfr signaling are involved in this process (Cheng et al., 2021). These results together supported the view that epithelial–myofibroblast crosstalk plays a key role in urothelial homeostasis and repair.
In summary, these results further reinforce the conclusion that bladder urothelial cells are highly heterogeneous and provide fundamental insights into the correlation between diverse cell phenotypes of the bladder urothelium and bladder diseases.
Single-cell sequencing in the cellular heterogeneity of bladder cancer
Discerning the heterogeneity of individual tumor cells carries clinical implications for predicting prognosis and treatment responses of bladder cancer. For instance, in the murine and human bladder cancer model systems, single-cell transcriptome analyses revealed the presence of multi-lineage gene expression phenotypes. In addition, among the epithelial cell clusters, cell subpopulations with concomitantly pronounced expressions of basal, luminal, and EMT claudin-related genes were characterised, and the cellular expression of these lineage markers was found to be dynamic either during tumor progression or in response to the treatment. Analysis of transplantable MIBC models demonstrated that these specific multi-lineage subtypes could define the cancer cells capable of tumor formation and lineage plasticity, which could be distinguished by the surface antigen CD49f in conjunction with the epithelial markers (Sfakianos et al., 2020). Finally, investigators affirmed the exact composition of the intra-tumoral cells derived from multi-lineage subpopulations as showing inherent heterogeneity and lineage plasticity of the epithelium. Data from a recent study revealed the heterogeneity of epithelium in bladder cancer patients with different pathological tumor stages (Lai et al., 2021). According to previous reports, bladder cancer is a urinary tract epithelium-derived cancer, and epithelial–mesenchymal transition (EMT) promotes the transition from NMIBC to MIBC (Cao et al., 2020). The malignant cell-specific gene expression suggested that basal-like cells were invasive precursors that displayed EMT characteristics in NMIBC. Recently, Wang et al. elucidated the spectrum of EMT associated with bladder cancer subtypes and identified that TCF7 promotes EMT in corroboration with single-cell ATAC sequencing (scATAC-seq) (Wang et al., 2021a). This evidence suggested that the epithelial cells represent the malignant cells of bladder cancer and that the upregulation of EMT is essential for driving local bladder tumor cells’ invasion and progression.
To further investigate the tumor heterogeneity involved in the transition from NMIBC to MIBC, Kerzeli et al. established a novel urothelial carcinoma model of Hgf-Cdk4R24C mice and identified eight heterogeneous cell clusters of urinary epithelium origin by utilising single-cell transcriptomic analyses (Kerzeli et al., 2021). Unlike the urothelial cells in healthy bladders, there were urothelial cell clusters with highly proliferative or cancer stemness gene expression characteristics in carcinogens-induced bladder tumors. Notably, these cells displaying cell cycle dysregulation gene expression profiles were considered asaggressive subclones. These results are in line with the reports of a study that observed intra-tumoral heterogeneity of human urothelial bladder cancer, identifying it as clonal mutations according to their genetic characteristics (Heide et al., 2019). Since the histopathological feature of Hgf-Cdk4R24C bladder tumors was urothelial carcinoma with squamous differentiation (UC/SCC), the proposed model would be a helpful tool for determining the tumor biological characteristics and the treatment response of UC/SCC.
Histologically, different from UC/SCC, tumors in bladder cancer patients with pure squamous cell components are defined as squamous cell carcinoma of the bladder (SCCB). Although SCCB accounts for <10% of primary bladder cancer, it shows frequent recurrence and metastasis when compared with urothelial carcinoma, which highlights the more complex intra-tumoral heterogeneity of SCCB (Minato et al., 2018; Lin et al., 2019). In 2016, a study conducted scRNA-seq on the tumor-normal paired tissues from an SCCB patient and revealed intra-tumoral heterogeneity and the potential mechanisms of SCCB (Zhang et al., 2016). The authors considered that genes with different expression patterns particularly enriched in the MAPK signaling pathway were involved in tumor evolution and heterogeneity formation. Especially, a keratinocyte-specific POU transcription factor, POU2F3, was associated with squamous epithelial stratification, indicating that POU2F3 may be a crucial biomarker of SCCB.
In general, these findings reveal the complex intra-tumoral and inter-tumor heterogeneities of bladder cancer, which can explain its diverse molecular and clinical phenotypes. These identified aggressive subclones and biomarkers may help enhance our understanding of the progression of bladder cancer, thereby facilitating individualised diagnosis and treatment.
Single-cell sequencing for investigating the genetic characteristics of bladder CSCs
Cancer stem cells (CSCs) possess the capability of self-renewal and are responsible for the initiation and development of tumors. Past studies have suggested that bladder CSCs possess high tumorigenicity, drug resistance, metastasis, and typical biomarkers (Ho et al., 2012). As early as 2012, Li et al. employed the single-cell exome sequencing technology to analyse the genetic landscape of bladder carcinoma and deciphered the mechanism of its occurrence and development (Li et al., 2012). Their study suggested that bladder cancer cells are derived from a single ancestral cell that formed two distinct tumor cell subgroups in the subsequent evolutionary process. Furthermore, they demonstrated that the ancestral bladder cancer cells had a monoclonal phenotype with multiple mutation driver gene candidates, among which ATM, COL6A3, and KIAA1958 were discovered as novel subclone-specific genes. These candidate cancer-related genes could drive the initiation of carcinogenesis and the development of subsequent cell lineages involved in cancer progression, which best matches the clonal evolution model. These findings explained that the genesis of bladder cancer was multi-gene mutation drive and multi-factorial, which in turn provided evidence for deciphering how tumors evolve into difficult-to-treat metastases. However, there remains an ambiguity related to whether ancestral cells inferred by the authors represent the initial tumor cells (Heng et al., 2006; Li et al., 2012; Ren et al., 2018; Zhang and Zhang, 2020).
To further investigate the genetic basis and origin of bladder CSCs, Yang et al. conducted single-cell sequencing on 59 cells from three human bladder cancer specimens and examined their phylogenetic status (Yang et al., 2017). Their research demonstrated that bladder CSCs originated from bladder epithelial stem cells or bladder cancer non-stem cells, thus providing genetic evidence to support the hypothesis of CSC origin. Moreover, 21 key gene mutations in bladder CSCs were discovered, which involved the genes related to cell cycle regulation, transcription regulation, chromatin remodeling, cell differentiation, and self-renewal. These specific mutations are critical to the acquisition of bladder CSC stemness. Moreover, the alterations of the three genes ARID1A, GPRC5A, and MLL2 play a crucial role in conferring stemness to the bladder cancer non-stem cells, as confirmed by co-mutation gene assays. Wang and others performed single-cell ATAC sequencing and found that bladder CSCs were enriched during tumor recurrence with elevated expression of EZH2. The specific molecular mechanism is that EZH2 maintains H3K27me3-mediated repression of the NCAM1 gene, thereby inactivating cellular invasiveness and stemness transcriptional programs (Wang et al., 2021a).
Overall, several potential key genes related to the characteristics of bladder CSCs have been uncovered so far. Nevertheless, further studies on the function and clinical applications are warranted to determine whether these candidate genes can act as bladder cancer biomarkers to guide targeted therapy.
Single-cell sequencing reveals the complexity of the bladder TME
The tumorigenesis, progression, and invasion of cancer cells are closely related to their surrounding microenvironment (Bryan, 2015; Song et al., 2019). Notably, bladder cancer is one of the least immune invasive cancers (Chen et al., 2017), and its complex TME may account for the poor response to immunotherapy. The composition of the bladder cancer TME was poorly understood. Recently, Chen et al. utilised the single-cell transcriptome sequencing technology to produce a cell atlas of the entire TME of bladder cancer (Chen et al., 2020). The authors found that bladder tumor cells express low levels of MHC-II molecules, suggesting that the downregulated immunogenicity may contribute to the avoidance of immune detection. Furthermore, the LAMP3+ dendritic cell subgroup expressed various cytokines, such as CCL17, CCL19, and CCL22, contributing to the formation of an immunosuppressive TME. Most importantly, two types of cancer-associated fibroblasts (CAFs) were identified. Among these fibroblasts, inflammatory CAFs (iCAFs) might promote the proliferation of bladder tumor cells and stromal cells and then recruit immune cells to the tumor region. Therefore, iCAFs were considered as the critical factor for the progression of bladder cancer and targeting CAFs might be an optimal choice for bladder cancer treatment in the future. Feng et al. integrated mass cytometry and imaging mass cytometry to investigate the MIBC TME at the single-cell proteomic level (Feng et al., 2021). They identified a specific CSC cluster (ALDH+ PD-L1+ ER-β−) that is associated with poor prognosis, improving the general understanding of the complexity of the TME in bladder cancer. Recently, Wang et al. also used single-cell mass cytometry to compare the characteristics of the TME between two groups MIBC patients (Wang et al., 2021b). They found that there was an immunosuppressive microenvironment in the immune-low group, with higher expression of PD-1 and Tim-3 on Tregs, and a higher proportion of PD-1+ Tregs.
Adaptive immune responses throughout the body are mainly coordinated by T lymphocytes. The phenotype of infiltrating T lymphocytes in the TME largely determines the response to immunotherapy (Hatogai and Sweis, 2020). Baras et al. for instance, indicated that the ratio of CD8 to Treg tumor-infiltrating lymphocyte density in the pre-treatment tissues can predict the response of bladder cancer to the platinum-based neoadjuvant chemotherapy (Baras et al., 2016). Notably, the current research on immunotherapy mainly focuses on the cytotoxic CD8+ T-cell-mediated response. Although the presence of cytotoxic CD4+ T cells has been validated in non-small cell lung cancer and hepatocellular carcinoma (Zheng et al., 2017; Guo et al., 2018), the extent of their heterogeneity and contribution to immunotherapy remains unclear. Through scRNA-seq and paired TCR-seq, Oh et al. analysed the tumor-infiltrating lymphocyte heterogeneity of bladder cancer (Oh et al., 2020). Notably, cytotoxic CD4+ T cells, rather than typical CD8+ T cells, were significantly enriched in bladder tumors. Cytotoxic CD4+ T cells recognised MHC-II antigens to kill the bladder tumor cells and lyse autologous tumor cells in a way that is inhibited by autologous Treg. Of note, signatures of these cytotoxic CD4+ T cells were significantly correlated with the clinical response to anti-PD-L1 therapy in bladder inflammation samples. Collectively, these findings have breakthrough research significance in the immunotherapy strategy of bladder cancer. In future treatments, the balance between cytotoxic CD4+ and regulatory T-cell status can be manipulated to provide therapeutic benefits to bladder cancer patients.
Aging can also lead to immune environment changes of the bladder, which induces corresponding bladder diseases (Maserejian et al., 2013). However, how and why bladder diseases become more prevalent with aging remains obscure. Single-cell transcriptomics analysis from the aged mouse bladders revealed the composition of the bladder immune cell repertoire, including novel subpopulations of macrophages and dendritic cells and notable changes in the immune repertoire. When compared with the young bladder tissues, T and B cells were highly enriched in aged bladders that constituted the organised bladder tertiary lymphoid tissues (bTLTs), which was supported by histological analyses (Ligon et al., 2020). Lymphoid aggregates have also been reported in some bladder tumor models, albeit they remain poorly characterised (Pitzalis et al., 2014). This report highlights that local immune environment dysfunction inside the bladder may be a potential mechanism that drives age-related bladder diseases. Moreover, the study provides an insight into the link among aging, inflammation, and bladder diseases.
In summary, the aforementioned single-cell studies have demonstrated the cell composition in the bladder TME and determined the molecular characteristics of immune cell subpopulations. Nevertheless, research in this field is still in its infancy.
Single-cell sequencing for studying urinary-exfoliated tumor cells of bladder cancer
Liquid biopsy for the detection of urinary-exfoliated tumor cells (UETCs) provides a hope for non-invasive screening and precision medicine for bladder cancer treatment. Before the development of bladder metastasis lesions, mutant tumor cells are released from the primary site of the bladder lesion or metastatic site to the patients’ urine. UETCs carry original information on each stage of the bladder cancer progression and have hence been utilised as potential biomarkers for urinary cytology detection (Critelli et al., 2016). However, UETCs are extremely rare and easily confounded by non-tumor cells, and the current related detection methods suffer from limited specificity or sensitivity (Yafi et al., 2015; Huang et al., 2021). To accurately and reliably identify UETCs, Chen et al. developed a novel microfluidic immunoassay method to separate and collect intact UETCs (Chen et al., 2018). This method combined the microfluidic immunoassay technology with two oncoproteins, namely CK20 and CD44v6, with frequent overexpression in the urothelial layers to achieve highly sensitive marking and separation of UETCs. Later, single-cell whole-genome sequencing and copy number variation (CNV) analysis were performed on 12 captured UETCs to achieve an accurate diagnosis of bladder cancer (with a specificity of 71.5%). To further improve the sensitivity of UETC detection and overcome the low-throughput limitation of single-cell DNA sequencing, Wang et al. performed single-cell CNV profile analysis combined with a cellular marker hexokinase 2 (HK2) to detect a total of 385 UETCs from eight urothelial carcinoma patients (Wang et al., 2020). Past studies have shown that, as a key enzyme in glucose metabolism, HK2 was significantly overexpressed in various types of cancers (Xu et al., 2018). Therefore, HK2 constitutes a promising biomarker for detecting urothelial cancer by identifying high glycolysis and metabolic abnormalities in UETCs. The study of the CNV profiles or oncogenic driver mutation signatures of bladder cancer, when combined with an HK2 threshold, achieved a specificity of >90% for UETC detection (Wang et al., 2020).
Overall, urine-based cellular or metabolic markers, when combined with single-cell sequencing, demonstrated superior sensitivity to the conventional urine cytology detection, which may enable high-sensitivity detection of UETCs in patients’ urine at the early and symptom-free stages.
Research on the therapeutic resistance of bladder cancer with single-cell sequencing
For non-metastatic bladder cancer patients, neoadjuvant cisplatin-based chemotherapy before radical cystectomy is the first-line of standard clinical treatment (Witjes et al., 2020). However, treatment failure due to drug resistance is common in bladder cancer, and approximately 40% of the patients show recurrence and metastasis after cystectomy (Kamoun et al., 2019). The second-line treatments for cisplatin-resistant and metastatic bladder cancer are the immune checkpoint inhibitors involving specific oncogene targeting agents, and approximately 20% of the patients exhibit a positive response to the treatment (Hargadon et al., 2018). Presently, the mechanisms of resistance to systemic therapy of bladder cancer are largely unknown, and molecular biomarkers for selecting effective second-line treatments are still lacking.
The emergence of drug resistance has been closely associated with intra-tumoral heterogeneity profiles. However, whether it results from pre-existing rare clones or new genome mutations remains controversial (Saunders et al., 2012). To investigate the mechanism of bladder cancer chemotherapy resistance, Tanaka et al. provided a single-cell atlas of bladder cancer treated with platinum-based chemotherapy and revealed the intra-tumoral heterogeneity states before and after treatment (Tanaka et al., 2018). By establishing an acquired platinum-resistant subline of urinary urothelial cancer cell line 5,637, the authors determined that 12 genes were considered as candidate platinum-resistance genes. Further evaluation of their clinical relevance demonstrated that only COX7B was associated with high mortality, while low levels of COX7B predicted a poor prognosis. Notably, the surrogate marker CD63 could distinguish low-COX7B subclones; thus, the authors speculated that COX7B and CD63 together contribute to platinum resistance in platinum-naïve bladder cancer.
The complex interactions between the TME and tumor cell-intrinsic mechanisms contribute to therapeutic failure and tumor evolution (Junttila and de Sauvage, 2013). To analyse tumor progression and TME changes during the clinical course of bladder cancer, Lee et al. depicted the tumor single-cell landscape of HRAS mutations in chemo-resistant metastatic MIBC (Lee et al., 2020). Clinically, the use of tipifarnib for advanced MIBC with HRAS mutations showed a significant therapeutic effect but was unable to achieve a complete response (Appels et al., 2005; Wang et al., 2017). Consequently, there is an urgent need to monitor the evolutionary trajectories of tumor cells and the surrounding TME during the clinical treatment course of MIBC. With the help of scRNA-seq and PDX models, a unique cell subpopulation with inherent resistance to tipifarnib was proven to exist at the time of tumor pre-treatment (Lee et al., 2020). Most importantly, the tumor subtypes with HRAS mutations show high expression levels of PD-L1, which indicates that alterations in tumor cells and TME by tipifarnib can be reversed by PD-1/PD-L1 immune checkpoint inhibitors. Therefore, the tipifarnib combined immunotherapy strategy may overcome the treatment failure and show favorable clinical responses.
Although immunotherapy offers great potential, not all MIBC patients can benefit from the current therapeutic regimens; in fact, improper treatment can increase mortality (Kamoun et al., 2019). For gaining in-depth understanding of the resistance mechanism of anti-PD-1/PD-L1 therapies, Wang et al. analysed the gene signatures of metastatic urothelial carcinoma resistance to PD-1/PD-L1 blockade through bulk transcriptome sequencing data (Wang et al., 2021c). Later, the cellular characteristics of the MIBC specimens were scrutinised through scRNA-seq. The results obtained suggested that the associated cells had diverse heterogeneity and that myeloid phagocyte could significantly upregulate the pro-tumorigenic inflammatory factors and downregulate the antigen-presenting genes. Hence, the authors anticipated that the PD-1/PD-L1 resistance of urothelial carcinoma was related to the balance of adaptive immunity and pro-tumorigenic inflammation in an individual TME.
In general, the occurrence of drug resistance in bladder cancer may be attributed to the pre-existing rare clones. Although the current findings are insufficient to explain the fundamental reasons for the treatment failure of bladder cancer, these studies can still provide novel insights into the drug resistance mechanism.
Single-cell sequencing in the prognosis of bladder cancer
A reliable prognostic prediction is essential to guide effective clinical treatment, which can safely avoid treatment failure and drug resistance (Wu et al., 2020). Zhou et al. collected the scRNA-seq profiles of 2075 cells from a bladder cancer patient using the GEO database and identified 14 different cell subgroups (Zhou et al., 2021). Then, with combined co-expression analyses and MATH value, 96 genes related to intra-tumoral heterogeneity, including several novel genes (such as KDELR3, GPSM3, RFC4, RPA3, IFI27L2, and APH1A) that may be correlated with bladder cancer progression, were screened. Moreover, a heterogeneity-related score model that could predict clinical treatment was established. This study revealed that the risk model based on scRNA-seq and intra-tumoral heterogeneity could help identify novel biomarkers and predict the clinical treatment outcomes, thereby providing new clues for understanding the intra-tumoral heterogeneity mechanisms of bladder cancer. Recently, Gouin III et al. using single-nucleus RNA-sequencing demonstrated that a novel epithelial cell phenotype marked by high expression of CDH12, which predicting poor prognosis for neoadjuvant chemotherapy (Gouin et al., 2021). In contrast, CDH12-enriched cells exhibit superior response to immune checkpoint therapy. CDH12-enriched cells highly express ligands for CD49a as well as PD-L1 and PD-L2, and co-localize with exhausted T-cells, providing one explanation for immune checkpoint blockade efficacy in these tumors.
The interaction between tumor cells and other cells in the TME has been widely studied (Wang et al., 2007); however, little attention has been paid to the communication among malignant cells. To analyse the guidance of cell communication between tumor cells on the prognosis, Lai et al. implemented single-cell analysis to study the ligand–receptor interaction in bladder cancer malignant cells (Lai et al., 2021). The authors discovered that six ligand genes and eight receptor genes were associated with the basal subtype tumors, while the high expression of the HBEGF-EGFR ligand–receptor pair predicted a poor prognosis. In addition, luminal subtype-related genes included three ligand genes and seven receptor genes, while the high expression of the FAM3B-LRP5 ligand–receptor pair indicated a more favorable outcome. These results highlight the relevance of cellular communication between tumor cells in the prognosis of bladder cancer for the first time.
Altogether, predictive models and biomarkers associated with the prognosis of bladder cancer have been identified so far. However, further prospective, large-scale, and multi-center clinical trials are warranted to validate their clinical utility.
Conclusion and perspectives
As a powerful high-throughput sequencing tool, single-cell techniques enable interrogating homeostatic and pathogenic cell populations with an exceptionally high resolution, thereby driving bladder cancer biomedical research to a new level of precision. These innovative methods have been extensively utilised to identify or interrogate the heterogeneity of bladder urothelial cells and bladder tumor cells, TME signatures, therapy resistance, UETCs, and CSCs (Figure 3). The application of scRNA-seq has provided a wide range of transcriptional expression profiles, which are essential for discovering rare cell subpopulations, discerning intra-tumor heterogeneity, monitoring the dynamic progression of bladder tumors, and revealing the drug resistance mechanisms. More importantly, these results of scRNA-seq could help identify novel carcinogenic driver biomarkers that may serve as therapeutic targets or prognostic factors ultimately guiding clinical therapies. In addition, single-cell DNA sequencing coupled with urine-based liquid biopsy has been found valuable in non-invasive screening for bladder cancer and monitoring the therapeutic efficacy. Overall, single-cell sequencing technologies have revolutionised our understanding of bladder tumor cytology, tumor biology, and tumor immunology, as they offer broader and deeper insights into prognosis prediction and treatment strategies to inspire future drug discovery.
FIGURE 3. Advances of single-cell sequencing in bladder cancer. Summary of the current application of single-cell sequencing technology in the cell heterogeneity, tumor microenvironment, drug resistance, cancer stem cells, urinary-exfoliated tumor cells, and prognosis of bladder cancer.
The past several years have witnessed unpredictable advances in single-cell sequencing technologies in the cancer field (Papalexi and Satija, 2018; Liu et al., 2021). However, the novel applications of single-cell technologies in characterising bladder cancer currently remain insufficient when compared with other human malignancies (Zang et al., 2021; Wei et al., 2022). First, existing single-cell sequencing studies have been limited to providing transcriptome or genomic information for a few specific types of bladder cancer patients, especially those lacking proteomic and cellular spatial information (Chen et al., 2018; Chen et al., 2020; Wang et al., 2020). Second, given the difficulty in obtaining the bladder urothelium of healthy individuals, single-cell sequencing studies on the link between the urothelium and bladder disease are limited to mice (Maserejian et al., 2013; Santos et al., 2019; Li et al., 2021b; Cheng et al., 2021). Third, drug resistance and metastasis often occur in bladder cancer patients, relevant single-cell research should not be limited to bladder cancer cell lines and PDX models (Tanaka et al., 2018; Lee et al., 2020). It remains indispensable to expand the study on the mechanism of bladder cancer drug resistance and metastasis through single-cell techniques in the future. As a result, our present understanding of bladder cancer cellular matrix components, TME, and drug resistance mechanisms are far from adequate. Therefore, it is a long and arduous process to comprehensively understand the underlying cellular or molecular mechanisms of bladder tumors, and plentiful studies need to be conducted.
Besides, due to the complex and expensive processes, single-cell sequencing could not be applied to profile large cohorts of tumor samples (Zhang et al., 2016; Chen et al., 2018; Chen et al., 2020; Oh et al., 2020; Wang et al., 2020). Even in real-world clinics, single-cell sequencing results cannot be applied to individualized therapy due to the rapidly deteriorating health status of patients (Lee et al., 2020). Notably, limited by the algorithm, it is difficult to identify accurate proportion of rare cell subsets in batch sequencing data (Li et al., 2012; Yang et al., 2017; Chen et al., 2020). This remains a big issue in bioinformatics analysis of single-cell sequencing data, and more proofs are needed to validate these single-cell analysis results. With the ongoing technological developments, single-cell sequencing is expected to become more high-throughput, with a low cost involved, simplified data analysis process, and reduced time expenditure, that can provide better accuracy and sensitivity to collectively promote the larger-scale single-cell sequencing efforts in bladder cancer and accelerate its clinical application. Furthermore, advanced techniques such as the application of single-cell multi-omics approaches and spatial transcriptomics will help us comprehensively reveal the pathogenesis and drug resistance mechanisms of bladder cancer. With the deepening of our understanding of the cellular dynamics of bladder tumor cells, the diagnosis and treatment of bladder cancer are expected to usher in a new era of precision medicine, improve the efficacy of personalised medicine, and ultimately save the lives of bladder cancer patients.
Author contributions
TL wrote the manuscript. JZ, YL, KW, and ZC took part in the manuscript modification. QZ helped with the preparation of the figure. All authors contributed to this paper and approved the submitted version.
Funding
This work was supported by the Postdoctoral Science Foundation of China (No. 2021M690158) and the “3315 Plan” Innovation Team Project of Ningbo City (No. 2019A-14-C).
Acknowledgments
Thanks to the contributions of all authors and the support fundings in this project.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
References
Antoni, S., Ferlay, J., Soerjomataram, I., Znaor, A., Jemal, A., and Bray, F. (2017). Bladder cancer incidence and mortality: A global overview and recent trends. Eur. Urol. 71, 96–108. doi:10.1016/j.eururo.2016.06.010
Apodaca, G. (2004). The uroepithelium: Not just a passive barrier. Traffic 5, 117–128. doi:10.1046/j.1600-0854.2003.00156.x
Appels, N. M., Beijnen, J. H., and Schellens, J. H. (2005). Development of farnesyl transferase inhibitors: A review. Oncologist 10, 565–578. doi:10.1634/theoncologist.10-8-565
Audenet, F., Attalla, K., and Sfakianos, J. P. (2018). The evolution of bladder cancer genomics: What have we learned and how can we use it? Urol. Oncol. 36, 313–320. doi:10.1016/j.urolonc.2018.02.017
Baras, A. S., Drake, C., Liu, J. J., Gandhi, N., Kates, M., Hoque, M. O., et al. (2016). The ratio of CD8 to Treg tumor-infiltrating lymphocytes is associated with response to cisplatin-based neoadjuvant chemotherapy in patients with muscle invasive urothelial carcinoma of the bladder. OncoImmunology 5, e1134412. doi:10.1080/2162402X.2015.1134412
Batista, R., Vinagre, N., Meireles, S., Vinagre, J., Prazeres, H., Leão, R., et al. (2020). Biomarkers for bladder cancer diagnosis and surveillance: A comprehensive review. Diagn. (Basel) 10, 39. doi:10.3390/diagnostics10010039
Bryan, R. T. (2015). Cell adhesion and urothelial bladder cancer: The role of cadherin switching and related phenomena. Philos. Trans. R. Soc. Lond. B Biol. Sci. 370, 20140042. doi:10.1098/rstb.2014.0042
Burger, M., Catto, J. W., Dalbagni, G., Grossman, H. B., Herr, H., Karakiewicz, P., et al. (2013). Epidemiology and risk factors of urothelial bladder cancer. Eur. Urol. 63, 234–241. doi:10.1016/j.eururo.2012.07.033
Cao, R., Yuan, L. S., Ma, B., Wang, G., Qiu, W., and Tian, Y. (2020). An EMT-related gene signature for the prognosis of human bladder cancer. J. Cell. Mol. Med. 24, 605–617. doi:10.1111/jcmm.14767
Casadei, C., Dizman, N., Schepisi, G., Cursano, M. C., Basso, U., Santini, D., et al. (2019). Targeted therapies for advanced bladder cancer: New strategies with FGFR inhibitors. Ther. Adv. Med. Oncol. 11, 1758835919890285. doi:10.1177/1758835919890285
Castillo-Martin, M., Domingo-Domenech, J., Karni-Schmidt, O., Matos, T., and Cordon-Cardo, C. (2010). Molecular pathways of urothelial development and bladder tumorigenesis. Urol. Oncol. 28, 401–408. doi:10.1016/j.urolonc.2009.04.019
Chen, A., Fu, G., Xu, Z., Sun, Y., Chen, X., Cheng, K., et al. (2018). Detection of urothelial bladder carcinoma via microfluidic immunoassay and single-cell DNA copy-number alteration analysis of captured urinary-exfoliated tumor cells. Cancer Res. 78, 4073–4085. doi:10.1158/0008-5472.CAN-17-2615
Chen, Y. P., Zhang, Y., Lv, J. W., Li, Y. Q., Wang, Y. Q., He, Q. M., et al. (2017). Genomic analysis of tumor microenvironment immune types across 14 solid cancer types: Immunotherapeutic implications. Theranostics 7, 3585–3594. doi:10.7150/thno.21471
Chen, Z., Zhou, L., Liu, L., Hou, Y., Xiong, M., Yang, Y., et al. (2020). Single-cell RNA sequencing highlights the role of inflammatory cancer-associated fibroblasts in bladder urothelial carcinoma. Nat. Commun. 11 (1), 5077. doi:10.1038/s41467-020-18916-5
Cheng, X., Lai, H., Luo, W., Zhang, M., Miao, J., Song, W., et al. (2021). Single-cell analysis reveals urothelial cell heterogeneity and regenerative cues following cyclophosphamide-induced bladder injury. Cell Death Dis. 12, 446. doi:10.1038/s41419-021-03740-6
Compérat, E. M., Burger, M., Gontero, P., Mostafid, A. H., Palou, J., Rouprêt, M., et al. (2019). Grading of urothelial carcinoma and the new “World Health Organisation classification of tumours of the urinary system and male genital organs 2016”. Eur. Urol. Focus 5, 457–466. doi:10.1016/j.euf.2018.01.003
Critelli, R., Fasanelli, F., Oderda, M., Polidoro, S., Assumma, M. B., Viberti, C., et al. (2016). Detection of multiple mutations in urinary exfoliated cells from male bladder cancer patients at diagnosis and during follow-up. Oncotarget 7, 67435–67448. doi:10.18632/oncotarget.11883
Feng, C., Wang, X., Tao, Y., Xie, Y., Lai, Z., Li, Z., et al. (2021). Single-cell proteomic analysis dissects the complexity of tumor microenvironment in muscle invasive bladder cancer. Cancers (Basel) 13 (21), 5440. doi:10.3390/cancers13215440
Gouin, K. H., Ing, N., Plummer, J. T., Rosser, C., Cheikh, B., Oh, C., et al. (2021). An N-Cadherin 2 expressing epithelial cell subpopulation predicts response to surgery, chemotherapy and immunotherapy in bladder cancer. Nat. Commun. 12 (1), 4906. doi:10.1038/s41467-021-25103-7
Guo, X., Zhang, Y., Zheng, L., Zheng, C., Song, J., Zhang, Q., et al. (2018). Global characterization of T cells in non-small-cell lung cancer by single-cell sequencing. Nat. Med. 24, 978–985. doi:10.1038/s41591-018-0045-3
Hargadon, K. M., Johnson, C. E., and Williams, C. J. (2018). Immune checkpoint blockade therapy for cancer: An overview of FDA-approved immune checkpoint inhibitors. Int. Immunopharmacol. 62, 29–39. doi:10.1016/j.intimp.2018.06.001
Hatogai, K., and Sweis, R. F. (2020). The tumor microenvironment of bladder cancer. Adv. Exp. Med. Biol. 1296, 275–290. doi:10.1007/978-3-030-59038-3_17
Heide, T., Maurer, A., Eipel, M., Knoll, K., Geelvink, M., Veeck, J., et al. (2019). Multiregion human bladder cancer sequencing reveals tumour evolution, bladder cancer phenotypes and implications for targeted therapy. J. Pathol. 248, 230–242. doi:10.1002/path.5250
Heng, H. H., Stevens, J. B., Liu, G., Bremer, S. W., Ye, K. J., Reddy, P. V., et al. (2006). Stochastic cancer progression driven by non-clonal chromosome aberrations. J. Cell. Physiol. 208 (2), 461–472. doi:10.1002/jcp.20685
Ho, P. L., Kurtova, A., and Chan, K. S. (2012). Normal and neoplastic urothelial stem cells: Getting to the root of the problem. Nat. Rev. Urol. 9, 583–594. doi:10.1038/nrurol.2012.142
Huang, R., Wang, L., He, J., and Gao, W. (2021). Application and prospects of single cell sequencing in tumors. Biomark. Res. 9 (1), 88. doi:10.1186/s40364-021-00336-2
Junttila, M. R., and de Sauvage, F. J. (2013). Influence of tumour micro-environment heterogeneity on therapeutic response. Nature 501, 346–354. doi:10.1038/nature12626
Kamoun, A., de Reyniès, A., Allory, Y., Sjödahl, G., Robertson, A. G., Seiler, R., et al. (2019). A consensus molecular classification of muscle-invasive bladder cancer. Eur. Urol. 77, 420–433. doi:10.1016/j.eururo.2019.09.006
Kantor, A. F., Hartge, P., Hoover, R. N., and Fraumeni, J. F. (1988). Epidemiological characteristics of squamous cell carcinoma and adenocarcinoma of the bladder. Cancer Res. 48, 3853–3855.
Kashima, Y., Sakamoto, Y., Kaneko, K., Seki, M., Suzuki, Y., and Suzuki, A. (2020). Single-cell sequencing techniques from individual to multiomics analyses. Exp. Mol. Med. 52 (9), 1419–1427. doi:10.1038/s12276-020-00499-2
Kerzeli, I. K., Lord, M., Doroszko, M., Elgendy, R., Chourlia, A., Stepanek, I., et al. (2021). Single-cell RNAseq and longitudinal proteomic analysis of a novel semi-spontaneous urothelial cancer model reveals tumor cell heterogeneity and pretumoral urine protein alterations. PLoS One 16, 0253178. doi:10.1371/journal.pone.0253178
Khandelwal, P., Abraham, S. N., and Apodaca, G. (2009). Cell biology and physiology of the uroepithelium. Am. J. Physiol. Ren. Physiol. 297, 1477–1501. doi:10.1152/ajprenal.00327.2009
Knowles, M., and Hurst, C. (2015). Molecular biology of bladder cancer: New insights into pathogenesis and clinical diversity. Nat. Rev. Cancer 15, 25–41. doi:10.1038/nrc3817
Lai, H., Cheng, X., Liu, Q., Luo, W., Liu, M., Zhang, M., et al. (2021). Single-cell RNA sequencing reveals the epithelial cell heterogeneity and invasive subpopulation in human bladder cancer. Int. J. Cancer 149, 2099–2115. doi:10.1002/ijc.33794
Lee, H. W., Chung, W., Lee, H. O., Jeong, D. E., Jo, A., Lim, J. E., et al. (2020). Single-cell RNA sequencing reveals the tumor microenvironment and facilitates strategic choices to circumvent treatment failure in a chemorefractory bladder cancer patient. Genome Med. 12, 47. doi:10.1186/s13073-020-00741-6
Lei, Y., Tang, R., Xu, J., Wang, W., Zhang, B., Liu, J., et al. (2021). Applications of single-cell sequencing in cancer research: Progress and perspectives. J. Hematol. Oncol. 14, 91. doi:10.1186/s13045-021-01105-2
Lerner, S. P., McConkey, D. J., Hoadley, K. A., Chan, K. S., Kim, W. Y., Radvanyi, F., et al. (2016). Bladder cancer molecular taxonomy: Summary from a consensus meeting. Bladder Cancer 2, 37–47. doi:10.3233/BLC-150037
Letašiová, S., Medve'ová, A., Šovčíková, A., Dušinská, M., Volkovová, K., Mosoiu, C., et al. (2012). Bladder cancer, a review of the environmental risk factors. Environ. Health 11, 11. doi:10.1186/1476-069X-11-S1-S11
Li, D. Y., Yang, F., Liao, W. Q., Zhou, X. F., Li, W. B., Cai, J. R., et al. (2021). Deep genomic sequencing of bladder urothelial carcinoma in southern Chinese patients: A single-center study. Front. Oncol. 11, 538927. doi:10.3389/fonc.2021.538927
Li, Y., Liu, Y., Gao, Z., Zhang, L., Chen, L., Wu, Z., et al. (2021). Single-cell transcriptomes of mouse bladder urothelium uncover novel cell type markers and urothelial differentiation characteristics. Cell Prolif. 54, 13007. doi:10.1111/cpr.13007
Li, Y., Xu, X., Song, L., Hou, Y., Li, Z., Tsang, S., et al. (2012). Single-cell sequencing analysis characterizes common and cell-lineage-specific mutations in a muscle-invasive bladder cancer. Gigascience 1, 12. doi:10.1186/2047-217X-1-12
Ligon, M. M., Wang, C., DeJong, E. N., Schulz, C., Bowdish, D., and Mysorekar, I. (2020). Single cell and tissue-transcriptomic analysis of murine bladders reveals age- and TNFα-dependent but microbiota-independent tertiary lymphoid tissue formation. Mucosal Immunol. 13, 908–918. doi:10.1038/s41385-020-0290-x
Lin, J., Spitz, M. R., Dinney, C. P., Etzel, C. J., Grossman, H. B., and Wu, X. (2006). Bladder cancer risk as modified by family history and smoking. Cancer 107, 705–711. doi:10.1002/cncr.22071
Lin, X., Deng, T., Wu, S., Lin, S., Wang, D., Wu, C-L., et al. (2019). The clinicopathological characteristics and prognostic value of squamous differentiation in patients with bladder urothelial carcinoma: A meta-analysis. World J. Urol. 38 (2), 323–333. doi:10.1007/s00345-019-02771-1
Liu, J., Qu, S., Zhang, T., Gao, Y., Shi, H., Song, K., et al. (2021). Applications of single-cell omics in tumor immunology. Front. Immunol. 12, 697412. doi:10.3389/fimmu.2021.697412
Maserejian, N. N., Chen, S., Chiu, G. R., Wager, C. G., Kupelian, V., Araujo, A. B., et al. (2013). Incidence of lower urinary tract symptoms in a population-based study of men and women. Urology 82 (3), 560–564. doi:10.1016/j.urology.2013.05.009
Minato, A., Noguchi, H., Tomisaki, I., Fukuda, A., Kubo, T., Nakayama, T., et al. (2018). Clinical significance of squamous differentiation in urothelial carcinoma of the bladder. Cancer control. 25, 1073274818800269. doi:10.1177/1073274818800269
Minoli, M., Kiener, M., Thalmann, G. N., Kruithof-de Julio, M., and Seiler, R. (2020). Evolution of urothelial bladder cancer in the context of molecular classifications. Int. J. Mol. Sci. 21, 5670. doi:10.3390/ijms21165670
Oh, D. Y., Kwek, S. S., Raju, S. S., Li, T., McCarthy, A., Chow, E., et al. (2020). Intratumoral CD4+ T cells mediate anti-tumor cytotoxicity in human bladder cancer. Cell 181, 1612–1625. e13. doi:10.1016/j.cell.2020.05.017
Papafotiou, G., Paraskevopoulou, V., Vasilaki, E., Kanaki, Z., Paschalidis, N., and Klinakis, A. (2016). KRT14 marks a subpopulation of bladder basal cells with pivotal role in regeneration and tumorigenesis. Nat. Commun. 7, 11914. doi:10.1038/ncomms11914
Papalexi, E., and Satija, R. (2018). Single-cell RNA sequencing to explore immune cell heterogeneity. Nat. Rev. Immunol. 18, 35–45. doi:10.1038/nri.2017.76
Pietzak, E. J., Bagrodia, A., Cha, E. K., Drill, E. N., Iyer, G., Isharwal, S., et al. (2017). Next-generation sequencing of nonmuscle invasive bladder cancer reveals potential biomarkers and rational therapeutic targets. Eur. Urol. 72, 952–959. doi:10.1016/j.eururo.2017.05.032
Pitzalis, C., Jones, G. W., Bombardieri, M., and Jones, S. A. (2014). Ectopic lymphoid-like structures in infection, cancer and autoimmunity. Nat. Rev. Immunol. 14, 447–462. doi:10.1038/nri3700
Ren, X., Kang, B., and Zhang, Z. (2018). Understanding tumor ecosystems by single-cell sequencing: Promises and limitations. Genome Biol. 19, 211. doi:10.1186/s13059-018-1593-z
Roy, S., Pradhan, D., Ernst, W. L., Mercurio, S., Najjar, Y., Parikh, R., et al. (2017). Next-generation sequencing-based molecular characterization of primary urinary bladder adenocarcinoma. Mod. Pathol. 30, 1133–1143. doi:10.1038/modpathol.2017.33
Santos, C. P., Lapi, E., Martínez de Villarreal, J., Álvaro-Espinosa, L., Fernández-Barral, A., Barbáchano, A., et al. (2019). Urothelial organoids originating from Cd49fhigh mouse stem cells display Notch-dependent differentiation capacity. Nat. Commun. 10, 4407. doi:10.1038/s41467-019-12307-1
Saunders, N. A., Simpson, F., Thompson, E. W., Hill, M. M., Endo-Munoz, L., Leggatt, G., et al. (2012). Role of intratumoural heterogeneity in cancer drug resistance: Molecular and clinical perspectives. EMBO Mol. Med. 4, 675–684. doi:10.1002/emmm.201101131
Sfakianos, J. P., Daza, J., Yang, H., Anastos, H., Bryant, G., Bareja, R., et al. (2020). Epithelial plasticity can generate multi-lineage phenotypes in human and murine bladder cancers. Nat. Commun. 11, 2540. doi:10.1038/s41467-020-16162-3
Shin, K., Lee, J., Guo, N., Kim, J., Lim, A., Qu, L., et al. (2011). Hedgehog/Wnt feedback supports regenerative proliferation of epithelial stem cells in bladder. Nature 472, 110–114. doi:10.1038/nature09851
Sjödahl, G., Eriksson, P., Liedberg, F., and Höglund, M. (2017). Molecular classification of urothelial carcinoma: Global mRNA classification versus tumour-cell phenotype classification. J. Pathol. 242, 113–125. doi:10.1002/path.4886
Sjödahl, G., Jackson, C. L., Bartlett, J. M., Siemens, D. R., and Berman, D. M. (2019). Molecular profiling in muscle-invasive bladder cancer: More than the sum of its parts. J. Pathol. 247, 563–573. doi:10.1002/path.5230
Song, D., Powles, T., Shi, L., Zhang, L., Ingersoll, , and Lu, Y. (2019). Bladder cancer, a unique model to understand cancer immunity and develop immunotherapy approaches. J. Pathol. 249 (2), 151–165. doi:10.1002/path.5306
Stuart, T., and Satija, R. (2019). Integrative single-cell analysis. Nat. Rev. Genet. 20, 257–272. doi:10.1038/s41576-019-0093-7
Tanaka, N., Katayama, S., Reddy, A., Nishimura, K., Niwa, N., Hongo, H., et al. (2018). Single‐cell RNA‐seq analysis reveals the platinum resistance gene COX7B and the surrogate marker CD63. Cancer Med. 7, 6193–6204. doi:10.1002/cam4.1828
Tang, F., Barbacioru, C., Wang, Y., Nordman, E., Lee, C., Xu, N., et al. (2009). mRNA-seq whole-transcriptome analysis of a single cell. Nat. Methods 6, 377–382. doi:10.1038/nmeth.1315
Thomsen, M. B. H., Nordentoft, I., Lamy, P., Vang, S., Reinert, L., Mapendano, C. K., et al. (2017). Comprehensive multiregional analysis of molecular heterogeneity in bladder cancer. Sci. Rep. 7, 11702. doi:10.1038/s41598-017-11291-0
Tran, L., Xiao, J. F., Agarwal, N., Duex, J. E., and Theodorescu, D. (2021). Advances in bladder cancer biology and therapy. Nat. Rev. Cancer 21, 104–121. doi:10.1038/s41568-020-00313-1
van den Bosch, S., and Alfred Witjes, J. (2011). Long-term cancer-specific survival in patients with high-risk, non-muscle-invasive bladder cancer and tumour progression: A systematic review. Eur. Urol. 60, 493–500. doi:10.1016/j.eururo.2011.05.045
Wang, H., Mei, Y., Luo, C., Huang, Q., Wang, Z., Lu, G., et al. (2021). Single-cell analyses reveal mechanisms of cancer stem cell maintenance and epithelial-mesenchymal transition in recurrent bladder cancer. Clin. Cancer Res. 27 (22), 6265–6278. doi:10.1158/1078-0432.CCR-20-4796
Wang, J., Batourina, E., Schneider, K., Souza, S., Swayne, T., Liu, C., et al. (2018). Polyploid superficial cells that maintain the urothelial barrier are produced via incomplete cytokinesis and endoreplication. Cell Rep. 25, 464–477. e4. doi:10.1016/j.celrep.2018.09.042
Wang, J., Yao, X., and Huang, J. (2017). New tricks for human farnesyltransferase inhibitor: Cancer and beyond. Medchemcomm 8, 841–854. doi:10.1039/c7md00030h
Wang, L., Sfakianos, J. P., Beaumont, K. G., Akturk, G., Horowitz, A., Sebra, R. P., et al. (2021). Myeloid cell-associated resistance to PD-1/PD-L1 blockade in urothelial cancer revealed through bulk and single-cell RNA sequencing. Clin. Cancer Res. 27 (15), 4287–4300. doi:10.1158/1078-0432.CCR-20-4574
Wang, P. F., Nishitani, M. A., Tanimoto, S., Kishimoto, T., Fukumori, T., Takahashi, M., et al. (2007). Bladder cancer cell invasion is enhanced by cross-talk with fibroblasts through hepatocyte growth factor. Urology 69, 780–784. doi:10.1016/j.urology.2007.01.063
Wang, X., Pan, L., Lu, Q., Huang, H., Feng, C., Tao, Y., et al. (2021). A combination of ssGSEA and mass cytometry identifies immune microenvironment in muscle-invasive bladder cancer. J. Clin. Lab. Anal. 35 (5), e23754. doi:10.1002/jcla.23754
Wang, Z., Chen, J., Yang, L., Cao, M., Yu, Y., Zhang, R., et al. (2020). Single-cell sequencing-enabled Hexokinase 2 assay for noninvasive bladder cancer diagnosis and screening by detecting rare malignant cells in urine. Anal. Chem. 92, 16284–16292. doi:10.1021/acs.analchem.0c04282
Wei, W., Rong, Y., Sanhe, L., Chunxiu, Y., Thokerunga, E., Cui, D., et al. (2022). Single-cell sequencing and its applications in bladder cancer. Expert Rev. Mol. Med. 24, e6. doi:10.1017/erm.2021.23
Witjes, J. A., Bruins, H. M., Cathomas, R., Compérat, E. M., Cowan, N. C., Gakis, G., et al. (2020). European association of urology guidelines on muscle-invasive and metastatic bladder cancer: Summary of the 2020 guidelines. Eur. Urol. 79, 82–104. doi:10.1016/j.eururo.2020.03.055
Wu, S., Ou, T., Xing, N., Lu, J., Wan, S., Wang, C., et al. (2019). Whole-genome sequencing identifies ADGRG6 enhancer mutations and FRS2 duplications as angiogenesis-related drivers in bladder cancer. Nat. Commun. 10, 720. doi:10.1038/s41467-019-08576-5
Wu, X., Lv, D., Cai, C., Zhao, Z., Wang, M., Chen, W., et al. (2020). A TP53-associated immune prognostic signature for the prediction of overall survival and therapeutic responses in muscle-invasive bladder cancer. Front. Immunol. 11, 590618. doi:10.3389/fimmu.2020.590618
Xing, Q. R., Cipta, N. O., Hamashima, K., Liou, Y. C., Koh, C. G., and Loh, Y. H. (2020). Unraveling heterogeneity in transcriptome and its regulation through single-cell multi-omics technologies. Front. Genet. 11, 662. doi:10.3389/fgene.2020.00662
Xu, S., Catapang, A., Braas, D., Stiles, L., Doh, H., Lee, J., et al. (2018). A precision therapeutic strategy for hexokinase 1-null, hexokinase 2-positive cancers. Cancer Metab. 6, 7. doi:10.1186/s40170-018-0181-8
Yafi, F. A., Brimo, F., Steinberg, J., Aprikian, A. G., Tanguay, S., and Kassouf, W. (2015). Prospective analysis of sensitivity and specificity of urinary cytology and other urinary biomarkers for bladder cancer. Urol. Oncol. 33, 6625–e31. doi:10.1016/j.urolonc.2014.06.008
Yang, Z., Li, C., Fan, Z., Liu, H., Zhang, X., Cai, Z., et al. (2017). Single-cell sequencing reveals variants in ARID1A, GPRC5A and MLL2 driving self-renewal of human bladder cancer stem cells. Eur. Urol. 71 (1), 8–12. doi:10.1016/j.eururo.2016.06.025
Yu, Z., Liao, J., Chen, Y., Zhang, L., Zou, C., Zhang, H., et al. (2019). Single-cell transcriptomic map of the human and mouse bladders. J. Am. Soc. Nephrol. 30, 2159–2176. doi:10.1681/ASN.2019040335
Zang, J., Ye, K., Fei, Y., Zhang, R., Chen, H., and Zhuang, G. (2021). Immunotherapy in the treatment of urothelial bladder cancer: Insights from single-cell analysis. Front. Oncol. 11, 696716. doi:10.3389/fonc.2021.696716
Zhang, X., Meng, Z., Yong, H., Xu, L., Li, W., Zou, Z., et al. (2016). Single-cell analyses of transcriptional heterogeneity in squamous cell carcinoma of urinary bladder. Oncotarget 7 (40), 66069–66076. doi:10.18632/oncotarget.11803
Zhang, Y., and Zhang, Z. (2020). The history and advances in cancer immunotherapy: Understanding the characteristics of tumor-infiltrating immune cells and their therapeutic implications. Cell. Mol. Immunol. 17, 807–821. doi:10.1038/s41423-020-0488-6
Zheng, C., Zheng, L., Yoo, J. K., Guo, H., Zhang, Y., Guo, X., et al. (2017). Landscape of infiltrating T cells in liver cancer revealed by single-cell sequencing. Cell 169, 1342–1356. e16. doi:10.1016/j.cell.2017.05.035
Zhou, R., Liang, J., Chen, Q., Tian, H., Yang, C., and Liu, C. (2021). Development and validation of an intra-tumor heterogeneity-related signature to predict prognosis of bladder cancer: A study based on single-cell RNA-seq. Aging 13, 19415–19441. doi:10.18632/aging.203353
Zhu, S., Yu, W., Yang, X., Wu, C., and Cheng, F. (2020). Traditional classification and novel subtyping systems for bladder cancer. Front. Oncol. 10, 102. doi:10.3389/fonc.2020.00102
Keywords: single-cell sequencing, bladder cancer, heterogeneity, biomarkers, precision medicine
Citation: Lyu T, Lin Y, Wu K, Cao Z, Zhang Q and Zheng J (2022) Single-cell sequencing technologies in bladder cancer research: Applications and challenges. Front. Genet. 13:1027909. doi: 10.3389/fgene.2022.1027909
Received: 25 August 2022; Accepted: 06 October 2022;
Published: 19 October 2022.
Edited by:
Sathiya Pandi Narayanan, Sidra Medicine, QatarReviewed by:
Liangkuan Bi, Second Affiliated Hospital of Anhui Medical University, ChinaHung N. Luu, University of Pittsburgh, United States
Copyright © 2022 Lyu, Lin, Wu, Cao, Zhang and Zheng. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Jianping Zheng, emhlbmdqaWFucGluZ0BuaW10ZS5hYy5jbg==