Niche-DE: niche-differential gene expression analysis in spatial transcriptomics data identifies context-dependent cell-cell interactions.

Existing methods for analysis of spatial transcriptomic data focus on delineating the global gene expression variations of cell types across the tissue, rather than local gene expression changes driven by cell-cell interactions. We propose a new statistical procedure called niche-differential expression (niche-DE) analysis that identifies cell-type-specific niche-associated genes, which are differentially expressed within a specific cell type in the context of specific spatial niches. We further develop niche-LR, a method to reveal ligand-receptor signaling mechanisms that underlie niche-differential gene expression patterns. Niche-DE and niche-LR are applicable to low-resolution spot-based spatial transcriptomics data and data that is single-cell or subcellular in resolution.

A clinical trial of therapeutic vaccination in lymphoma with serial tumor sampling and single-cell analysis.

ABSTRACT: In situ vaccination (ISV) triggers an immune response to tumor-associated antigens at 1 tumor site, which can then tackle the disease throughout the body. Here, we report clinical and biological results of a phase 1/2 ISV trial in patients with low-grade lymphoma, combining an intratumoral toll-like receptor 9 (TLR9) agonist with local low-dose radiation and ibrutinib (an inhibitor of B- and T-cell kinases). Adverse events were predominately low grade. The overall response rate was 50%, including 1 complete response. All patients experienced tumor reduction at distant sites. Single-cell analyses of serial fine needle aspirates from injected and uninjected tumors revealed correlates of clinical response, such as lower CD47 and higher major histocompatibility complex class II expression on tumor cells, enhanced T-cell and natural killer cell effector function, and reduced immune suppression from transforming growth factor ß and inhibitory T regulatory 1 cells. Although changes at the local injected site were more pronounced, changes at distant uninjected sites were more often associated with clinical responses. Functional immune response assays and tracking of T-cell receptor sequences provided evidence of treatment-induced tumor-specific T-cell responses. Induction of immune effectors and reversal of negative regulators were both important in producing clinically meaningful tumor responses. The trial was registered at www.clinicaltrials.gov as #NCT02927964.

GITR and TIGIT immunotherapy provokes divergent multicellular responses in the tumor microenvironment of gastrointestinal cancers.

BACKGROUND: Understanding the mechanistic effects of novel immunotherapy agents is critical to improving their successful clinical translation. These effects need to be studied in preclinical models that maintain the heterogenous tumor microenvironment (TME) and dysfunctional cell states found in a patient's tumor. We investigated immunotherapy perturbations targeting co-stimulatory molecule GITR and co-inhibitory immune checkpoint TIGIT in a patient-derived ex vivo system that maintains the TME in its near-native state. Leveraging single-cell genomics, we identified cell type-specific transcriptional reprogramming in response to immunotherapy perturbations. METHODS: We generated ex vivo tumor slice cultures from fresh surgical resections of gastric and colon cancer and treated them with GITR agonist or TIGIT antagonist antibodies. We applied paired single-cell RNA and TCR sequencing to the original surgical resections, control, and treated ex vivo tumor slice cultures. We additionally confirmed target expression using multiplex immunofluorescence and validated our findings with RNA in situ hybridization. RESULTS: We confirmed that tumor slice cultures maintained the cell types, transcriptional cell states and proportions of the original surgical resection. The GITR agonist was limited to increasing effector gene expression only in cytotoxic CD8 T cells. Dysfunctional exhausted CD8 T cells did not respond to GITR agonist. In contrast, the TIGIT antagonist increased TCR signaling and activated both cytotoxic and dysfunctional CD8 T cells. This included cells corresponding to TCR clonotypes with features indicative of potential tumor antigen reactivity. The TIGIT antagonist also activated T follicular helper-like cells and dendritic cells, and reduced markers of immunosuppression in regulatory T cells. CONCLUSIONS: We identified novel cellular mechanisms of action of GITR and TIGIT immunotherapy in the patients' TME. Unlike the GITR agonist that generated a limited transcriptional response, TIGIT antagonist orchestrated a multicellular response involving CD8 T cells, T follicular helper-like cells, dendritic cells, and regulatory T cells. Our experimental strategy combining single-cell genomics with preclinical models can successfully identify mechanisms of action of novel immunotherapy agents. Understanding the cellular and transcriptional mechanisms of response or resistance will aid in prioritization of targets and their clinical translation.

A spatially mapped gene expression signature for intestinal stem-like cells identifies high-risk precursors of gastric cancer.

Objective: Gastric intestinal metaplasia (GIM) is a precancerous lesion that increases gastric cancer (GC) risk. The Operative Link on GIM (OLGIM) is a combined clinical-histopathologic system to risk-stratify patients with GIM. The identification of molecular biomarkers that are indicators for advanced OLGIM lesions may improve cancer prevention efforts. Methods: This study was based on clinical and genomic data from four cohorts: 1) GAPS, a GIM cohort with detailed OLGIM severity scoring (N=303 samples); 2) the Cancer Genome Atlas (N=198); 3) a collation of in-house and publicly available scRNA-seq data (N=40), and 4) a spatial validation cohort (N=5) consisting of annotated histology slides of patients with either GC or advanced GIM. We used a multi-omics pipeline to identify, validate and sequentially parse a highly-refined signature of 26 genes which characterize high-risk GIM. Results: Using standard RNA-seq, we analyzed two separate, non-overlapping discovery (N=88) and validation (N=215) sets of GIM. In the discovery phase, we identified 105 upregulated genes specific for high-risk GIM (defined as OLGIM III-IV), of which 100 genes were independently confirmed in the validation set. Spatial transcriptomic profiling revealed 36 of these 100 genes to be expressed in metaplastic foci in GIM. Comparison with bulk GC sequencing data revealed 26 of these genes to be expressed in intestinal-type GC. Single-cell profiling resolved the 26-gene signature to both mature intestinal lineages (goblet cells, enterocytes) and immature intestinal lineages (stem-like cells). A subset of these genes was further validated using single-molecule multiplex fluorescence in situ hybridization. We found certain genes (TFF3 and ANPEP) to mark differentiated intestinal lineages, whereas others (OLFM4 and CPS1) localized to immature cells in the isthmic/crypt region of metaplastic glands, consistent with the findings from scRNAseq analysis. Conclusions: using an integrated multi-omics approach, we identified a novel 26-gene expression signature for high-OLGIM precursors at increased risk for GC. We found this signature localizes to aberrant intestinal stem-like cells within the metaplastic microenvironment. These findings hold important translational significance for future prevention and early detection efforts.

Follicular lymphoma evolves with a surmountable dependency on acquired glycosylation motifs in the B-cell receptor.

ABSTRACT: An early event in the genesis of follicular lymphoma (FL) is the acquisition of new glycosylation motifs in the B-cell receptor (BCR) due to gene rearrangement and/or somatic hypermutation. These N-linked glycosylation motifs (N-motifs) contain mannose-terminated glycans and can interact with lectins in the tumor microenvironment, activating the tumor BCR pathway. N-motifs are stable during FL evolution, suggesting that FL tumor cells are dependent on them for their survival. Here, we investigated the dynamics and potential impact of N-motif prevalence in FL at the single-cell level across distinct tumor sites and over time in 17 patients. Although most patients had acquired at least 1 N-motif as an early event, we also found (1) cases without N-motifs in the heavy or light chains at any tumor site or time point and (2) cases with discordant N-motif patterns across different tumor sites. Inferring phylogenetic trees of the patients with discordant patterns, we observed that both N-motif-positive and N-motif-negative tumor subclones could be selected and expanded during tumor evolution. Comparing N-motif-positive with N-motif-negative tumor cells within a patient revealed higher expression of genes involved in the BCR pathway and inflammatory response, whereas tumor cells without N-motifs had higher activity of pathways involved in energy metabolism. In conclusion, although acquired N-motifs likely support FL pathogenesis through antigen-independent BCR signaling in most patients with FL, N-motif-negative tumor cells can also be selected and expanded and may depend more heavily on altered metabolism for competitive survival.

Direct measurement of engineered cancer mutations and their transcriptional phenotypes in single cells.

Genome sequencing studies have identified numerous cancer mutations across a wide spectrum of tumor types, but determining the phenotypic consequence of these mutations remains a challenge. Here, we developed a high-throughput, multiplexed single-cell technology called TISCC-seq to engineer predesignated mutations in cells using CRISPR base editors, directly delineate their genotype among individual cells and determine each mutation's transcriptional phenotype. Long-read sequencing of the target gene's transcript identifies the engineered mutations, and the transcriptome profile from the same set of cells is simultaneously analyzed by short-read sequencing. Through integration, we determine the mutations' genotype and expression phenotype at single-cell resolution. Using cell lines, we engineer and evaluate the impact of >100 TP53 mutations on gene expression. Based on the single-cell gene expression, we classify the mutations as having a functionally significant phenotype.

Single-cell multi-gene identification of somatic mutations and gene rearrangements in cancer.

In this proof-of-concept study, we developed a single-cell method that provides genotypes of somatic alterations found in coding regions of messenger RNAs and integrates these transcript-based variants with their matching cell transcriptomes. We used nanopore adaptive sampling on single-cell complementary DNA libraries to validate coding variants in target gene transcripts, and short-read sequencing to characterize cell types harboring the mutations. CRISPR edits for 16 targets were identified using a cancer cell line, and known variants in the cell line were validated using a 352-gene panel. Variants in primary cancer samples were validated using target gene panels ranging from 161 to 529 genes. A gene rearrangement was also identified in one patient, with the rearrangement occurring in two distinct tumor sites.

Single cell and spatial alternative splicing analysis with long read sequencing.

Long-read sequencing has become a powerful tool for alternative splicing analysis. However, technical and computational challenges have limited our ability to explore alternative splicing at single cell and spatial resolution. The higher sequencing error of long reads, especially high indel rates, have limited the accuracy of cell barcode and unique molecular identifier (UMI) recovery. Read truncation and mapping errors, the latter exacerbated by the higher sequencing error rates, can cause the false detection of spurious new isoforms. Downstream, there is yet no rigorous statistical framework to quantify splicing variation within and between cells/spots. In light of these challenges, we developed Longcell, a statistical framework and computational pipeline for accurate isoform quantification for single cell and spatial spot barcoded long read sequencing data. Longcell performs computationally efficient cell/spot barcode extraction, UMI recovery, and UMI-based truncation- and mapping-error correction. Through a statistical model that accounts for varying read coverage across cells/spots, Longcell rigorously quantifies the level of inter-cell/spot versus intra-cell/ spot diversity in exon-usage and detects changes in splicing distributions between cell populations. Applying Longcell to single cell long-read data from multiple contexts, we found that intra-cell splicing heterogeneity, where multiple isoforms co-exist within the same cell, is ubiquitous for highly expressed genes. On matched single cell and Visium long read sequencing for a tissue of colorectal cancer metastasis to the liver, Longcell found concordant signals between the two data modalities. Finally, on a perturbation experiment for 9 splicing factors, Longcell identified regulatory targets that are validated by targeted sequencing.

GITR and TIGIT immunotherapy provokes divergent multi-cellular responses in the tumor microenvironment of gastrointestinal cancers.

Understanding the cellular mechanisms of novel immunotherapy agents in the human tumor microenvironment (TME) is critical to their clinical success. We examined GITR and TIGIT immunotherapy in gastric and colon cancer patients using ex vivo slice tumor slice cultures derived from cancer surgical resections. This primary culture system maintains the original TME in a near-native state. We applied paired single-cell RNA and TCR sequencing to identify cell type specific transcriptional reprogramming. The GITR agonist was limited to increasing effector gene expression only in cytotoxic CD8 T cells. The TIGIT antagonist increased TCR signaling and activated both cytotoxic and dysfunctional CD8 T cells, including clonotypes indicative of potential tumor antigen reactivity. The TIGIT antagonist also activated T follicular helper-like cells and dendritic cells, and reduced markers of immunosuppression in regulatory T cells. Overall, we identified cellular mechanisms of action of these two immunotherapy targets in the patients' TME.

Colorectal Cancer Metastases in the Liver Establish Immunosuppressive Spatial Networking between Tumor-Associated SPP1+ Macrophages and Fibroblasts.

PURPOSE: The liver is the most frequent metastatic site for colorectal cancer. Its microenvironment is modified to provide a niche that is conducive for colorectal cancer cell growth. This study focused on characterizing the cellular changes in the metastatic colorectal cancer (mCRC) liver tumor microenvironment (TME). EXPERIMENTAL DESIGN: We analyzed a series of microsatellite stable (MSS) mCRCs to the liver, paired normal liver tissue, and peripheral blood mononuclear cells using single-cell RNA sequencing (scRNA-seq). We validated our findings using multiplexed spatial imaging and bulk gene expression with cell deconvolution. RESULTS: We identified TME-specific SPP1-expressing macrophages with altered metabolism features, foam cell characteristics, and increased activity in extracellular matrix (ECM) organization. SPP1+ macrophages and fibroblasts expressed complementary ligand-receptor pairs with the potential to mutually influence their gene-expression programs. TME lacked dysfunctional CD8 T cells and contained regulatory T cells, indicative of immunosuppression. Spatial imaging validated these cell states in the TME. Moreover, TME macrophages and fibroblasts had close spatial proximity, which is a requirement for intercellular communication and networking. In an independent cohort of mCRCs in the liver, we confirmed the presence of SPP1+ macrophages and fibroblasts using gene-expression data. An increased proportion of TME fibroblasts was associated with the worst prognosis in these patients. CONCLUSIONS: We demonstrated that mCRC in the liver is characterized by transcriptional alterations of macrophages in the TME. Intercellular networking between macrophages and fibroblasts supports colorectal cancer growth in the immunosuppressed metastatic niche in the liver. These features can be used to target immune-checkpoint-resistant MSS tumors.

Integrative single-cell analysis of allele-specific copy number alterations and chromatin accessibility in cancer.

Cancer progression is driven by both somatic copy number aberrations (CNAs) and chromatin remodeling, yet little is known about the interplay between these two classes of events in shaping the clonal diversity of cancers. We present Alleloscope, a method for allele-specific copy number estimation that can be applied to single-cell DNA- and/or transposase-accessible chromatin-sequencing (scDNA-seq, ATAC-seq) data, enabling combined analysis of allele-specific copy number and chromatin accessibility. On scDNA-seq data from gastric, colorectal and breast cancer samples, with validation using matched linked-read sequencing, Alleloscope finds pervasive occurrence of highly complex, multiallelic CNAs, in which cells that carry varying allelic configurations adding to the same total copy number coevolve within a tumor. On scATAC-seq from two basal cell carcinoma samples and a gastric cancer cell line, Alleloscope detected multiallelic copy number events and copy-neutral loss-of-heterozygosity, enabling dissection of the contributions of chromosomal instability and chromatin remodeling to tumor evolution.

Single-cell analysis can define distinct evolution of tumor sites in follicular lymphoma.

Tumor heterogeneity complicates biomarker development and fosters drug resistance in solid malignancies. In lymphoma, our knowledge of site-to-site heterogeneity and its clinical implications is still limited. Here, we profiled 2 nodal, synchronously acquired tumor samples from 10 patients with follicular lymphoma (FL) using single-cell RNA, B-cell receptor (BCR) and T-cell receptor sequencing, and flow cytometry. By following the rapidly mutating tumor immunoglobulin genes, we discovered that BCR subclones were shared between the 2 tumor sites in some patients, but in many patients, the disease had evolved separately with limited tumor cell migration between the sites. Patients exhibiting divergent BCR evolution also exhibited divergent tumor gene-expression and cell-surface protein profiles. While the overall composition of the tumor microenvironment did not differ significantly between sites, we did detect a specific correlation between site-to-site tumor heterogeneity and T follicular helper (Tfh) cell abundance. We further observed enrichment of particular ligand-receptor pairs between tumor and Tfh cells, including CD40 and CD40LG, and a significant correlation between tumor CD40 expression and Tfh proliferation. Our study may explain discordant responses to systemic therapies, underscores the difficulty of capturing a patient's disease with a single biopsy, and furthers our understanding of tumor-immune networks in FL.

Joint single cell DNA-seq and RNA-seq of gastric cancer cell lines reveals rules of in vitro evolution.

Cancer cell lines are not homogeneous nor are they static in their genetic state and biological properties. Genetic, transcriptional and phenotypic diversity within cell lines contributes to the lack of experimental reproducibility frequently observed in tissue-culture-based studies. While cancer cell line heterogeneity has been generally recognized, there are no studies which quantify the number of clones that coexist within cell lines and their distinguishing characteristics. We used a single-cell DNA sequencing approach to characterize the cellular diversity within nine gastric cancer cell lines and integrated this information with single-cell RNA sequencing. Overall, we sequenced the genomes of 8824 cells, identifying between 2 and 12 clones per cell line. Using the transcriptomes of more than 28 000 single cells from the same cell lines, we independently corroborated 88% of the clonal structure determined from single cell DNA analysis. For one of these cell lines, we identified cell surface markers that distinguished two subpopulations and used flow cytometry to sort these two clones. We identified substantial proportions of replicating cells in each cell line, assigned these cells to subclones detected among the G0/G1 population and used the proportion of replicating cells per subclone as a surrogate of each subclone's growth rate.

Single-Cell Genomic Characterization Reveals the Cellular Reprogramming of the Gastric Tumor Microenvironment.

PURPOSE: The tumor microenvironment (TME) consists of a heterogenous cellular milieu that can influence cancer cell behavior. Its characteristics have an impact on treatments such as immunotherapy. These features can be revealed with single-cell RNA sequencing (scRNA-seq). We hypothesized that scRNA-seq analysis of gastric cancer together with paired normal tissue and peripheral blood mononuclear cells (PBMC) would identify critical elements of cellular deregulation not apparent with other approaches. EXPERIMENTAL DESIGN: scRNA-seq was conducted on seven patients with gastric cancer and one patient with intestinal metaplasia. We sequenced 56,167 cells comprising gastric cancer (32,407 cells), paired normal tissue (18,657 cells), and PBMCs (5,103 cells). Protein expression was validated by multiplex immunofluorescence. RESULTS: Tumor epithelium had copy number alterations, a distinct gene expression program from normal, with intratumor heterogeneity. Gastric cancer TME was significantly enriched for stromal cells, macrophages, dendritic cells (DC), and Tregs. TME-exclusive stromal cells expressed distinct extracellular matrix components than normal. Macrophages were transcriptionally heterogenous and did not conform to a binary M1/M2 paradigm. Tumor DCs had a unique gene expression program compared to PBMC DCs. TME-specific cytotoxic T cells were exhausted with two heterogenous subsets. Helper, cytotoxic T, Treg, and NK cells expressed multiple immune checkpoint or co-stimulatory molecules. Receptor-ligand analysis revealed TME-exclusive intercellular communication. CONCLUSIONS: Single-cell gene expression studies revealed widespread reprogramming across multiple cellular elements in the gastric cancer TME. Cellular remodeling was delineated by changes in cell numbers, transcriptional states, and intercellular interactions. This characterization facilitates understanding of tumor biology and enables identification of novel targets including for immunotherapy.

scPred: accurate supervised method for cell-type classification from single-cell RNA-seq data.

Single-cell RNA sequencing has enabled the characterization of highly specific cell types in many tissues, as well as both primary and stem cell-derived cell lines. An important facet of these studies is the ability to identify the transcriptional signatures that define a cell type or state. In theory, this information can be used to classify an individual cell based on its transcriptional profile. Here, we present scPred, a new generalizable method that is able to provide highly accurate classification of single cells, using a combination of unbiased feature selection from a reduced-dimension space, and machine-learning probability-based prediction method. We apply scPred to scRNA-seq data from pancreatic tissue, mononuclear cells, colorectal tumor biopsies, and circulating dendritic cells and show that scPred is able to classify individual cells with high accuracy. The generalized method is available at https://github.com/powellgenomicslab/scPred/.

Functional genomics identifies predictive markers and clinically actionable resistance mechanisms to CDK4/6 inhibition in bladder cancer.

BACKGROUND: CDK4/6 inhibitors are a promising treatment strategy in tumor therapy but are hampered by resistance mechanisms. This study was performed to reveal predictive markers, mechanisms of resistance and to develop rational combination therapies for a personalized therapy approach in bladder cancer. METHODS: A genome-scale CRISPR-dCas9 activation screen for resistance to the CDK4/6 inhibitor Palbociclib was performed in the bladder cancer derived cell line T24. sgRNA counts were analyzed using next generation sequencing and MAGeCK-VISPR. Significantly enriched sgRNAs were cloned and validated on a molecular and functional level for mediating resistance to Palbociclib treatment. Analysis was done in vitro and in vivo in the chorioallantois membrane model of the chicken embryo. Comparison of screen hits to signaling pathways and clinically relevant molecular alterations was performed using DAVID, Reactome, DGIdb and cBioPortal. RESULTS: In the screen, 1024 sgRNAs encoding for 995 genes were significantly enriched indicative of mediators of resistance. 8 random sgRNAs were validated, revealing partial rescue to Palbociclib treatment. Within this gene panel, members of Receptor-Tyrosine Kinases, PI3K-Akt, Ras/MAPK, JAK/STAT or Wnt signaling pathways were identified. Combination of Palbociclib with inhibitors against these signaling pathways revealed beneficial effects in vitro and in in vivo xenografts. CONCLUSIONS: Identification of potential predictive markers, resistance mechanisms and rational combination therapies could be achieved by applying a CRISPR-dCas9 screening approach in bladder cancer.

Parallel PI3K, AKT and mTOR inhibition is required to control feedback loops that limit tumor therapy.

Targeting the PI3K pathway has achieved limited success in cancer therapy. One reason for the disappointing activity of drugs that interfere with molecules that are important player in this pathway is the induction of multiple feedback loops that have been only partially understood. To understand these limitations and develop improved treatment strategies, we comprehensively characterized molecular mechanisms of PI3K pathway signaling in bladder cancer cell lines upon using small molecule inhibitors and RNAi technologies against all key molecules and protein complexes within the pathway and analyzed functional and molecular consequences. When targeting either mTORC1, mTOR, AKT or PI3K, only S6K1 phosphorylation was affected in most cell lines examined. Dephosphorylation of 4E-BP1 required combined inhibition of PI3K and mTORC1, independent from AKT, and resulted in a robust reduction in cell viability. Long-term inhibition of PI3K however resulted in a PDK1-dependent, PIP3 and mTORC2 independent rephosphorylation of AKT. AKT rephosphorylation could also be induced by mTOR or PDK1 inhibition. Combining PI3K/mTOR inhibitors with AKT or PDK1 inhibitors suppressed this rephosphorylation, induced apoptosis, decreased colony formation, cell viability and growth of tumor xenografts. Our findings reveal novel molecular mechanisms that explain the requirement for simultaneous targeting of PI3K, AKT and mTORC1 to achieve effective tumor growth inhibition.

Targeting the PI3K/AKT/mTOR Pathway in Bladder Cancer.

The PI3K/AKT/mTOR signaling pathway shows frequent molecular alterations and increased activity in cancer. Given its role in the regulation of cell growth, survival and metastasis, molecules within this pathway are promising targets for pharmacologic intervention. Metastatic bladder cancer (BLCA) continues to have few treatment options. Although various molecular alterations in PI3K/AKT/mTOR signaling have been described in BLCA, clinical trials with small molecule inhibitors have not met their endpoints. In this article, we summarize results from preclinical studies and clinical trials that examined PI3K pathway inhibitors in BLCA focusing on technical challenges that might result in contradictory findings in preclinical studies. Based on published data from our group, we also address challenges that need to be overcome to optimize PI3K inhibition in BLCA and enable its successful translation into the clinic.

Non-invasive imaging of engineered human tumors in the living chicken embryo.

The growing interest in engineered tumor models prompted us to devise a method for the non-invasive assessment of such models. Here, we report on bioluminescence imaging (BLI) for the assessment of engineered tumor models in the fertilized chicken egg, i.e, chick chorioallantoic membrane (CAM) assay. One prostate cancer (PC-3) and two osteosarcoma (MG63 and HOS) cell lines were modified with luciferase reporter genes. To create engineered tumors, these cell lines were seeded either onto basement membrane extract (BME) or gelfoam scaffolds, and subsequently grafted in vivo onto the CAM. BLI enabled non-invasive, specific detection of the engineered tumors on the CAM in the living chicken embryo. Further, BLI permitted daily, quantitative monitoring of the engineered tumors over the course of up to 7 days. Data showed that an extracellular matrix (ECM) composed of BME supported growth of reporter gene marked PC-3 tumors but did not support MG63 or HOS tumor growth. However, MG63 tumors engineered on the collagen-based gelfoam ECM showed a temporal proliferation burst in MG63 tumors. Together, the data demonstrated imaging of engineered human cancer models in living chicken embryos. The combination of CAM assay and BLI holds significant potential for the examination of a broad range of engineered tumor models.

CDK4/6 Inhibitors in Cancer Therapy: A Novel Treatement Strategy for Bladder Cancer.

Patients with metastatic bladder cancer (mBC) treated with cisplatin-based chemotherapy have a limited median survival of only around 14 months [1]. Despite over 30 years of basic and clinical research, until recently no therapeutic options beyond cisplatin-based therapy had entered clinical routine and, at least in the US, none of the tested agents had been approved for second-line treatment. This has changed with the advent of immune checkpoint blockade, including especially PD-1/PD-L1 inhibitors. The high response rates of 24% over a 14.4 month follow up led to the first US Food and Drug Administration (FDA) approval for a second line therapy for these patients, and it is likely that this marks the beginning of a new era in the systemic treatment of muscle-invasive bladder cancer [2-4]. The strong clinical need to improve the medical management of this disease for those patients, not responding to current therapy has led to an increased molecular understanding of bladder cancer and has forstered the development of many potential molecular manipulations and targeted strategies beyond the new immune-oncologic approaches. Among the molecular alterations indentified in bladder cancer, cell cycle deregulation appears to be a key driver of disease progression. Target-directed therapy against CDK4/6 is an emerging strategy to regain control of cell cycle deregulation. Here, we provide an overview of the current status of CDK4/6 inhibitors in cancer therapy, their potential use in mBC and the challenges for their clinical use.