Single Cell Analysis of Treatment-Resistant Prostate Cancer: Implications of Cell State Changes for Cell Surface Antigen Targeted Therapies.

Targeting cell surface molecules using radioligand and antibody-based therapies has yielded considerable success across cancers. However, it remains unclear how the expression of putative lineage markers, particularly cell surface molecules, varies in the process of lineage plasticity, wherein tumor cells alter their identity and acquire new oncogenic properties. A notable example of lineage plasticity is the transformation of prostate adenocarcinoma (PRAD) to neuroendocrine prostate cancer (NEPC)--a growing resistance mechanism that results in the loss of responsiveness to androgen blockade and portends dismal patient survival. To understand how lineage markers vary across the evolution of lineage plasticity in prostate cancer, we applied single cell analyses to 21 human prostate tumor biopsies and two genetically engineered mouse models, together with tissue microarray analysis (TMA) on 131 tumor samples. Not only did we observe a higher degree of phenotypic heterogeneity in castrate-resistant PRAD and NEPC than previously anticipated, but also found that the expression of molecules targeted therapeutically, namely PSMA, STEAP1, STEAP2, TROP2, CEACAM5, and DLL3, varied within a subset of gene-regulatory networks (GRNs). We also noted that NEPC and small cell lung cancer (SCLC) subtypes shared a set of GRNs, indicative of conserved biologic pathways that may be exploited therapeutically across tumor types. While this extreme level of transcriptional heterogeneity, particularly in cell surface marker expression, may mitigate the durability of clinical responses to novel antigen-directed therapies, its delineation may yield signatures for patient selection in clinical trials, potentially across distinct cancer types.

Increasing Fluid Viscosity Ensures Consistent Single-Cell Encapsulation.

High-throughput single-cell analysis typically relies on the isolation of cells of interest in separate compartments for subsequent phenotypic or genotypic characterization. Using microfluidics, this is achieved by isolating individual cells in microdroplets or microwells. However, due to cell-to-cell variability in size, shape, and density, the cell capture efficiencies may vary significantly. This variability can negatively impact the measurements and introduce undesirable artifacts when trying to isolate and characterize heterogeneous cell populations. In this study, we show that single-cell isolation biases in microfluidics can be circumvented by increasing the viscosity of fluids in which cells are dispersed. At a viscosity of 40-50 cP (cP), the cell sedimentation is effectively reduced, resulting in a steady cell flow inside the microfluidics chip and consistent encapsulation in water-in-oil droplets over extended periods of time. This approach allows nearly all cells in a sample to be isolated with the same efficiency, irrespective of their type. Our results show that increased fluid viscosity, rather than cell-adjusted density, provides a more reliable approach to mitigate single-cell isolation biases.

Deep generative model deciphers derailed trajectories in acute myeloid leukemia.

Single-cell genomics has the potential to map cell states and their dynamics in an unbiased way in response to perturbations like disease. However, elucidating the cell-state transitions from healthy to disease requires analyzing data from perturbed samples jointly with unperturbed reference samples. Existing methods for integrating and jointly visualizing single-cell datasets from distinct contexts tend to remove key biological differences or do not correctly harmonize shared mechanisms. We present Decipher, a model that combines variational autoencoders with deep exponential families to reconstruct derailed trajectories (https://github.com/azizilab/decipher). Decipher jointly represents normal and perturbed single-cell RNA-seq datasets, revealing shared and disrupted dynamics. It further introduces a novel approach to visualize data, without the need for methods such as UMAP or TSNE. We demonstrate Decipher on data from acute myeloid leukemia patient bone marrow specimens, showing that it successfully characterizes the divergence from normal hematopoiesis and identifies transcriptional programs that become disrupted in each patient when they acquire NPM1 driver mutations.

Epigenetic plasticity cooperates with cell-cell interactions to direct pancreatic tumorigenesis.

The response to tumor-initiating inflammatory and genetic insults can vary among morphologically indistinguishable cells, suggesting as yet uncharacterized roles for epigenetic plasticity during early neoplasia. To investigate the origins and impact of such plasticity, we performed single-cell analyses on normal, inflamed, premalignant, and malignant tissues in autochthonous models of pancreatic cancer. We reproducibly identified heterogeneous cell states that are primed for diverse, late-emerging neoplastic fates and linked these to chromatin remodeling at cell-cell communication loci. Using an inference approach, we revealed signaling gene modules and tissue-level cross-talk, including a neoplasia-driving feedback loop between discrete epithelial and immune cell populations that was functionally validated in mice. Our results uncover a neoplasia-specific tissue-remodeling program that may be exploited for pancreatic cancer interception.

Conserved transcriptional connectivity of regulatory T cells in the tumor microenvironment informs new combination cancer therapy strategies.

While regulatory T (Treg) cells are traditionally viewed as professional suppressors of antigen presenting cells and effector T cells in both autoimmunity and cancer, recent findings of distinct Treg cell functions in tissue maintenance suggest that their regulatory purview extends to a wider range of cells and is broader than previously assumed. To elucidate tumoral Treg cell 'connectivity' to diverse tumor-supporting accessory cell types, we explored immediate early changes in their single-cell transcriptomes upon punctual Treg cell depletion in experimental lung cancer and injury-induced inflammation. Before any notable T cell activation and inflammation, fibroblasts, endothelial and myeloid cells exhibited pronounced changes in their gene expression in both cancer and injury settings. Factor analysis revealed shared Treg cell-dependent gene programs, foremost, prominent upregulation of VEGF and CCR2 signaling-related genes upon Treg cell deprivation in either setting, as well as in Treg cell-poor versus Treg cell-rich human lung adenocarcinomas. Accordingly, punctual Treg cell depletion combined with short-term VEGF blockade showed markedly improved control of PD-1 blockade-resistant lung adenocarcinoma progression in mice compared to the corresponding monotherapies, highlighting a promising factor-based querying approach to elucidating new rational combination treatments of solid organ cancers.

RNA cytometry of single-cells using semi-permeable microcapsules.

Analytical tools for gene expression profiling of individual cells are critical for studying complex biological systems. However, the techniques enabling rapid measurements of gene expression on thousands of single-cells are lacking. Here, we report a high-throughput RNA cytometry for digital profiling of single-cells isolated in liquid droplets enveloped by a thin semi-permeable membrane (microcapsules). Due to the selective permeability of the membrane, the desirable enzymes and reagents can be loaded, or replaced, in the microcapsule at any given step by simply changing the reaction buffer in which the microcapsules are dispersed. Therefore, complex molecular biology workflows can be readily adapted to conduct nucleic acid analysis on encapsulated mammalian cells, or other biological species. The microcapsules support sequential multi-step enzymatic reactions and remain intact under different biochemical conditions, freezing, thawing, and thermocycling. Combining microcapsules with conventional FACS provides a high-throughput approach for conducting RNA cytometry of individual cells based on their digital gene expression signature.

Protocol to dissociate, process, and analyze the human lung tissue using single-cell RNA-seq.

We report a protocol for obtaining high-quality single-cell transcriptomics data from human lung biospecimens acquired from core needle biopsies, fine-needle aspirates, surgical resection, and pleural effusions. The protocol relies upon the brief mechanical and enzymatic disruption of tissue, enrichment of live cells by fluorescence-activated cell sorting (FACS), and droplet-based single-cell RNA sequencing (scRNA-seq). The protocol also details a procedure for analyzing the scRNA-seq data. For complete details on the use and execution of this protocol, please refer to Chan et al. (2021).

Lineage plasticity in prostate cancer depends on JAK/STAT inflammatory signaling.

Drug resistance in cancer is often linked to changes in tumor cell state or lineage, but the molecular mechanisms driving this plasticity remain unclear. Using murine organoid and genetically engineered mouse models, we investigated the causes of lineage plasticity in prostate cancer and its relationship to antiandrogen resistance. We found that plasticity initiates in an epithelial population defined by mixed luminal-basal phenotype and that it depends on increased Janus kinase (JAK) and fibroblast growth factor receptor (FGFR) activity. Organoid cultures from patients with castration-resistant disease harboring mixed-lineage cells reproduce the dependency observed in mice by up-regulating luminal gene expression upon JAK and FGFR inhibitor treatment. Single-cell analysis confirms the presence of mixed-lineage cells with increased JAK/STAT (signal transducer and activator of transcription) and FGFR signaling in a subset of patients with metastatic disease, with implications for stratifying patients for clinical trials.

High-throughput single-cell antibody secretion quantification and enrichment using droplet microfluidics-based FRET assay.

High-throughput screening and enrichment of antibody-producing cells have many important applications. Herein, we present a droplet microfluidic approach for high-throughput screening and sorting of antibody-secreting cells using a Förster resonance electron transfer (FRET)-based assay. The FRET signal is mediated by the specific binding of the secreted antibody to two fluorescently labeled probes supplied within a droplet. Functional hybridoma cells expressing either membrane-bound or secreted monoclonal antibodies (mAbs), or both, were efficiently differentiated in less than 30 min. The antibody secretion rate by individual hybridoma cells was recorded in the range of 14,000 Abs/min, while the density of membrane-bound fraction was approximately 100 Abs/µm2. Combining the FRET assay with droplet-based single-cell sorting, an 800-fold enrichment of antigen-specific cells was achieved after one round of sorting. The presented system overcomes several key limitations observed in conventional FACS-based screening methods and should be applicable to assaying various other secreted proteins.

Microfluidics for Cancer Biomarker Discovery, Research, and Clinical Application.

Currently, cancer is the leading cause of death and its incidence and mortality is growing rapidly all over the world. One of the confounding factors contributing to the failure of conventional cancer diagnostics and treatment strategies is a high degree of intratumoral and intertumoral heterogeneity at the single-cell and molecular levels. Recent innovations in microfluidic techniques have revolutionized single-cell and single-molecule research and challenged the conventional definition of a "biomarker." Alongside classic cancer biomarkers such as circulating tumor DNA or circulating tumor cells (CTC), tumor cell heterogeneity, transcriptional and epigenetic cell states and their abundance in the tumor microenvironment have been demonstrated to impact disease progression and treatment response. Utilizing high-throughput, robust microfluidic techniques for the detection, isolation, and analysis of various cancer biomarkers, valuable information about the tumor can be obtained for clinical decision-making. This chapter presents clinically relevant advances of cancer biomarker research using microfluidics technology and identifies the emerging applications for disease diagnosis, monitoring, and personalized treatment.

Single-Cell Transcriptional Profiling Reveals Signatures of Helper, Effector, and Regulatory MAIT Cells during Homeostasis and Activation.

Mucosal-associated invariant T (MAIT) cells are innate-like lymphocytes that recognize microbial vitamin B metabolites and have emerging roles in infectious disease, autoimmunity, and cancer. Although MAIT cells are identified by a semi-invariant TCR, their phenotypic and functional heterogeneity is not well understood. Here we present an integrated single cell transcriptomic analysis of over 76,000 human MAIT cells during early and prolonged Ag-specific activation with the MR1 ligand 5-OP-RU and nonspecific TCR stimulation. We show that MAIT cells span a broad range of homeostatic, effector, helper, tissue-infiltrating, regulatory, and exhausted phenotypes, with distinct gene expression programs associated with CD4+ or CD8+ coexpression. During early activation, MAIT cells rapidly adopt a cytotoxic phenotype characterized by high expression of GZMB, IFNG and TNF In contrast, prolonged stimulation induces heterogeneous states defined by proliferation, cytotoxicity, immune modulation, and exhaustion. We further demonstrate a FOXP3 expressing MAIT cell subset that phenotypically resembles conventional regulatory T cells. Moreover, scRNAseq-defined MAIT cell subpopulations were also detected in individuals recently exposed to Mycobacterium tuberculosis, confirming their presence during human infection. To our knowledge, our study provides the first comprehensive atlas of human MAIT cells in activation conditions and defines substantial functional heterogeneity, suggesting complex roles in health and disease.

AD-linked R47H-TREM2 mutation induces disease-enhancing microglial states via AKT hyperactivation.

The hemizygous R47H variant of triggering receptor expressed on myeloid cells 2 (TREM2), a microglia-specific gene in the brain, increases risk for late-onset Alzheimer's disease (AD). Using transcriptomic analysis of single nuclei from brain tissues of patients with AD carrying the R47H mutation or the common variant (CV)­TREM2, we found that R47H-associated microglial subpopulations had enhanced inflammatory signatures reminiscent of previously identified disease-associated microglia (DAM) and hyperactivation of AKT, one of the signaling pathways downstream of TREM2. We established a tauopathy mouse model with heterozygous knock-in of the human TREM2 with the R47H mutation or CV and found that R47H induced and exacerbated TAU-mediated spatial memory deficits in female mice. Single-cell transcriptomic analysis of microglia from these mice also revealed transcriptomic changes induced by R47H that had substantial overlaps with R47H microglia in human AD brains, including robust increases in proinflammatory cytokines, activation of AKT signaling, and elevation of a subset of DAM signatures. Pharmacological AKT inhibition with MK-2206 largely reversed the enhanced inflammatory signatures in primary R47H microglia treated with TAU fibrils. In R47H heterozygous tauopathy mice, MK-2206 treatment abolished a tauopathy-dependent microglial subcluster and rescued tauopathy-induced synapse loss. By uncovering disease-enhancing mechanisms of the R47H mutation conserved in human and mouse, our study supports inhibitors of AKT signaling as a microglial modulating strategy to treat AD.

Inhibition of Carbonic Anhydrase IX Suppresses Breast Cancer Cell Motility at the Single-Cell Level.

Protein Carbonic Anhydrase IX (CA IX), which is expressed in various hypoxic solid tumors in order to maintain proper pH, is also related to cancer cell adhesion, invasion, and metastasis processes. Here, we investigated whether CA IX inhibition by a highly CA IX selective agent benzenesulfonamide VD11-4-2 triggers changes in individual cell motility. We seeded breast cancer cells on an extracellular matrix-coated glass-bottomed dish and in a microfluidic device with a gradient flow of epidermal growth factor (EGF), tracked individual cell movement, calculated their migration speeds, and/or followed movement direction. Our results showed that the inhibitor VD11-4-2 decreased the speed of CA IX positive breast cancer cells by 20-26% while not affecting non-cancerous cell migration. The inhibitor suppressed the cell migration velocity increment and hindered cells from reaching their maximum speed. VD11-4-2 also reduced CA IX, expressing cell movement towards the growth factor as a chemoattractant. Such a single cell-based migration assay enabled the comprehensive investigation of the cell motility and revealed that VD11-4-2 shows the ability to suppress breast cancer cell migration at a lower concentration than previously tested CA IX inhibitors.

Signatures of plasticity, metastasis, and immunosuppression in an atlas of human small cell lung cancer.

Small cell lung cancer (SCLC) is an aggressive malignancy that includes subtypes defined by differential expression of ASCL1, NEUROD1, and POU2F3 (SCLC-A, -N, and -P, respectively). To define the heterogeneity of tumors and their associated microenvironments across subtypes, we sequenced 155,098 transcriptomes from 21 human biospecimens, including 54,523 SCLC transcriptomes. We observe greater tumor diversity in SCLC than lung adenocarcinoma, driven by canonical, intermediate, and admixed subtypes. We discover a PLCG2-high SCLC phenotype with stem-like, pro-metastatic features that recurs across subtypes and predicts worse overall survival. SCLC exhibits greater immune sequestration and less immune infiltration than lung adenocarcinoma, and SCLC-N shows less immune infiltrate and greater T cell dysfunction than SCLC-A. We identify a profibrotic, immunosuppressive monocyte/macrophage population in SCLC tumors that is particularly associated with the recurrent, PLCG2-high subpopulation.

A unified atlas of CD8 T cell dysfunctional states in cancer and infection.

CD8 T cells play an essential role in defense against viral and bacterial infections and in tumor immunity. Deciphering T cell loss of functionality is complicated by the conspicuous heterogeneity of CD8 T cell states described across experimental and clinical settings. By carrying out a unified analysis of over 300 assay for transposase-accessible chromatin sequencing (ATAC-seq) and RNA sequencing (RNA-seq) experiments from 12 studies of CD8 T cells in cancer and infection, we defined a shared differentiation trajectory toward dysfunction and its underlying transcriptional drivers and revealed a universal early bifurcation of functional and dysfunctional T cell states across models. Experimental dissection of acute and chronic viral infection using single-cell ATAC (scATAC)-seq and allele-specific single-cell RNA (scRNA)-seq identified state-specific drivers and captured the emergence of similar TCF1+ progenitor-like populations at an early branch point, at which functional and dysfunctional T cells diverge. Our atlas of CD8 T cell states will facilitate mechanistic studies of T cell immunity and translational efforts.

A gene-environment-induced epigenetic program initiates tumorigenesis.

Tissue damage increases the risk of cancer through poorly understood mechanisms1. In mouse models of pancreatic cancer, pancreatitis associated with tissue injury collaborates with activating mutations in the Kras oncogene to markedly accelerate the formation of early neoplastic lesions and, ultimately, adenocarcinoma2,3. Here, by integrating genomics, single-cell chromatin assays and spatiotemporally controlled functional perturbations in autochthonous mouse models, we show that the combination of Kras mutation and tissue damage promotes a unique chromatin state in the pancreatic epithelium that distinguishes neoplastic transformation from normal regeneration and is selected for throughout malignant evolution. This cancer-associated epigenetic state emerges within 48 hours of pancreatic injury, and involves an 'acinar-to-neoplasia' chromatin switch that contributes to the early dysregulation of genes that define human pancreatic cancer. Among the factors that are most rapidly activated after tissue damage in the pre-malignant pancreatic epithelium is the alarmin cytokine interleukin 33, which recapitulates the effects of injury in cooperating with mutant Kras to unleash the epigenetic remodelling program of early neoplasia and neoplastic transformation. Collectively, our study demonstrates how gene-environment interactions can rapidly produce gene-regulatory programs that dictate early neoplastic commitment, and provides a molecular framework for understanding the interplay between genetic and environmental cues in the initiation of cancer.

Fully defined human pluripotent stem cell-derived microglia and tri-culture system model C3 production in Alzheimer's disease.

Aberrant inflammation in the CNS has been implicated as a major player in the pathogenesis of human neurodegenerative disease. We developed a new approach to derive microglia from human pluripotent stem cells (hPSCs) and built a defined hPSC-derived tri-culture system containing pure populations of hPSC-derived microglia, astrocytes, and neurons to dissect cellular cross-talk along the neuroinflammatory axis in vitro. We used the tri-culture system to model neuroinflammation in Alzheimer's disease with hPSCs harboring the APPSWE+/+ mutation and their isogenic control. We found that complement C3, a protein that is increased under inflammatory conditions and implicated in synaptic loss, is potentiated in tri-culture and further enhanced in APPSWE+/+ tri-cultures due to microglia initiating reciprocal signaling with astrocytes to produce excess C3. Our study defines the major cellular players contributing to increased C3 in Alzheimer's disease and presents a broadly applicable platform to study neuroinflammation in human disease.

Multi-step processing of single cells using semi-permeable capsules.

Droplet microfluidics technology provides a powerful approach to isolate and process millions of single cells simultaneously. Despite many exciting applications that have emerged based on this technology, workflows based on multi-step operations, including molecular biology and cell-based phenotypic screening assays, cannot be easily adapted to droplet format. Here, we present a microfluidics-based technique to isolate single cells, or biological samples, into semi-permeable hydrogel capsules and perform multi-step biological workflows on thousands to millions of individual cells simultaneously. The biochemical reactions are performed by changing the aqueous buffer surrounding the capsules, without needing sophisticated equipment. The semi-permeable nature of the capsules' shell retains large encapsulated biomolecules (such as genome) while allowing smaller molecules (such as proteins) to passively diffuse. In contrast to conventional hydrogel bead assays, the approach presented here improves bacterial cell retention during multi-step procedures as well as the efficiency of biochemical reactions. We showcase two examples of capsule use for single genome amplification of bacteria, and expansion of individual clones into isogenic microcolonies for later screening for biodegradable plastic production.

Molecular Fingerprint and Developmental Regulation of the Tegmental GABAergic and Glutamatergic Neurons Derived from the Anterior Hindbrain.

Tegmental nuclei in the ventral midbrain and anterior hindbrain control motivated behavior, mood, memory, and movement. These nuclei contain inhibitory GABAergic and excitatory glutamatergic neurons, whose molecular diversity and development remain largely unraveled. Many tegmental neurons originate in the embryonic ventral rhombomere 1 (r1), where GABAergic fate is regulated by the transcription factor (TF) Tal1. We used single-cell mRNA sequencing of the mouse ventral r1 to characterize the Tal1-dependent and independent neuronal precursors. We describe gene expression dynamics during bifurcation of the GABAergic and glutamatergic lineages and show how active Notch signaling promotes GABAergic fate selection in post-mitotic precursors. We identify GABAergic precursor subtypes that give rise to distinct tegmental nuclei and demonstrate that Sox14 and Zfpm2, two TFs downstream of Tal1, are necessary for the differentiation of specific tegmental GABAergic neurons. Our results provide a framework for understanding the development of cellular diversity in the tegmental nuclei.