Mapping transcriptomic vector fields of single cells.

Single-cell (sc)RNA-seq, together with RNA velocity and metabolic labeling, reveals cellular states and transitions at unprecedented resolution. Fully exploiting these data, however, requires kinetic models capable of unveiling governing regulatory functions. Here, we introduce an analytical framework dynamo (https://github.com/aristoteleo/dynamo-release), which infers absolute RNA velocity, reconstructs continuous vector fields that predict cell fates, employs differential geometry to extract underlying regulations, and ultimately predicts optimal reprogramming paths and perturbation outcomes. We highlight dynamo's power to overcome fundamental limitations of conventional splicing-based RNA velocity analyses to enable accurate velocity estimations on a metabolically labeled human hematopoiesis scRNA-seq dataset. Furthermore, differential geometry analyses reveal mechanisms driving early megakaryocyte appearance and elucidate asymmetrical regulation within the PU.1-GATA1 circuit. Leveraging the least-action-path method, dynamo accurately predicts drivers of numerous hematopoietic transitions. Finally, in silico perturbations predict cell-fate diversions induced by gene perturbations. Dynamo, thus, represents an important step in advancing quantitative and predictive theories of cell-state transitions.

Differential encoding in prefrontal cortex projection neuron classes across cognitive tasks.

Single-cell transcriptomics has been widely applied to classify neurons in the mammalian brain, while systems neuroscience has historically analyzed the encoding properties of cortical neurons without considering cell types. Here we examine how specific transcriptomic types of mouse prefrontal cortex (PFC) projection neurons relate to axonal projections and encoding properties across multiple cognitive tasks. We found that most types projected to multiple targets, and most targets received projections from multiple types, except PFC→PAG (periaqueductal gray). By comparing Ca2+ activity of the molecularly homogeneous PFC→PAG type against two heterogeneous classes in several two-alternative choice tasks in freely moving mice, we found that all task-related signals assayed were qualitatively present in all examined classes. However, PAG-projecting neurons most potently encoded choice in cued tasks, whereas contralateral PFC-projecting neurons most potently encoded reward context in an uncued task. Thus, task signals are organized redundantly, but with clear quantitative biases across cells of specific molecular-anatomical characteristics.

Therapy-Induced Evolution of Human Lung Cancer Revealed by Single-Cell RNA Sequencing.

Lung cancer, the leading cause of cancer mortality, exhibits heterogeneity that enables adaptability, limits therapeutic success, and remains incompletely understood. Single-cell RNA sequencing (scRNA-seq) of metastatic lung cancer was performed using 49 clinical biopsies obtained from 30 patients before and during targeted therapy. Over 20,000 cancer and tumor microenvironment (TME) single-cell profiles exposed a rich and dynamic tumor ecosystem. scRNA-seq of cancer cells illuminated targetable oncogenes beyond those detected clinically. Cancer cells surviving therapy as residual disease (RD) expressed an alveolar-regenerative cell signature suggesting a therapy-induced primitive cell-state transition, whereas those present at on-therapy progressive disease (PD) upregulated kynurenine, plasminogen, and gap-junction pathways. Active T-lymphocytes and decreased macrophages were present at RD and immunosuppressive cell states characterized PD. Biological features revealed by scRNA-seq were biomarkers of clinical outcomes in independent cohorts. This study highlights how therapy-induced adaptation of the multi-cellular ecosystem of metastatic cancer shapes clinical outcomes.

In situ tumour arrays reveal early environmental control of cancer immunity.

The immune phenotype of a tumour is a key predictor of its response to immunotherapy1-4. Patients who respond to checkpoint blockade generally present with immune-inflamed5-7 tumours that are highly infiltrated by T cells. However, not all inflamed tumours respond to therapy, and even lower response rates occur among tumours that lack T cells (immune desert) or that spatially exclude T cells to the periphery of the tumour lesion (immune excluded)8. Despite the importance of these tumour immune phenotypes in patients, little is known about their development, heterogeneity or dynamics owing to the technical difficulty of tracking these features in situ. Here we introduce skin tumour array by microporation (STAMP)-a preclinical approach that combines high-throughput time-lapse imaging with next-generation sequencing of tumour arrays. Using STAMP, we followed the development of thousands of arrayed tumours in vivo to show that tumour immune phenotypes and outcomes vary between adjacent tumours and are controlled by local factors within the tumour microenvironment. Particularly, the recruitment of T cells by fibroblasts and monocytes into the tumour core was supportive of T cell cytotoxic activity and tumour rejection. Tumour immune phenotypes were dynamic over time and an early conversion to an immune-inflamed phenotype was predictive of spontaneous or therapy-induced tumour rejection. Thus, STAMP captures the dynamic relationships of the spatial, cellular and molecular components of tumour rejection and has the potential to translate therapeutic concepts into successful clinical strategies.

Molecular hallmarks of heterochronic parabiosis at single-cell resolution.

The ability to slow or reverse biological ageing would have major implications for mitigating disease risk and maintaining vitality1. Although an increasing number of interventions show promise for rejuvenation2, their effectiveness on disparate cell types across the body and the molecular pathways susceptible to rejuvenation remain largely unexplored. Here we performed single-cell RNA sequencing on 20 organs to reveal cell-type-specific responses to young and aged blood in heterochronic parabiosis. Adipose mesenchymal stromal cells, haematopoietic stem cells and hepatocytes are among those cell types that are especially responsive. On the pathway level, young blood invokes new gene sets in addition to reversing established ageing patterns, with the global rescue of genes encoding electron transport chain subunits pinpointing a prominent role of mitochondrial function in parabiosis-mediated rejuvenation. We observed an almost universal loss of gene expression with age that is largely mimicked by parabiosis: aged blood reduces global gene expression, and young blood restores it in select cell types. Together, these data lay the groundwork for a systemic understanding of the interplay between blood-borne factors and cellular integrity.

Single-cell proteomics enabled by next-generation sequencing or mass spectrometry.

In the last decade, single-cell RNA sequencing routinely performed on large numbers of single cells has greatly advanced our understanding of the underlying heterogeneity of complex biological systems. Technological advances have also enabled protein measurements, further contributing to the elucidation of cell types and states present in complex tissues. Recently, there have been independent advances in mass spectrometric techniques bringing us one step closer to characterizing single-cell proteomes. Here we discuss the challenges of detecting proteins in single cells by both mass spectrometry and sequencing-based methods. We review the state of the art for these techniques and propose that there is a space for technological advancements and complementary approaches that maximize the advantages of both classes of technologies.

Ageing hallmarks exhibit organ-specific temporal signatures.

Ageing is the single greatest cause of disease and death worldwide, and understanding the associated processes could vastly improve quality of life. Although major categories of ageing damage have been identified-such as altered intercellular communication, loss of proteostasis and eroded mitochondrial function1-these deleterious processes interact with extraordinary complexity within and between organs, and a comprehensive, whole-organism analysis of ageing dynamics has been lacking. Here we performed bulk RNA sequencing of 17 organs and plasma proteomics at 10 ages across the lifespan of Mus musculus, and integrated these findings with data from the accompanying Tabula Muris Senis2-or 'Mouse Ageing Cell Atlas'-which follows on from the original Tabula Muris3. We reveal linear and nonlinear shifts in gene expression during ageing, with the associated genes clustered in consistent trajectory groups with coherent biological functions-including extracellular matrix regulation, unfolded protein binding, mitochondrial function, and inflammatory and immune response. Notably, these gene sets show similar expression across tissues, differing only in the amplitude and the age of onset of expression. Widespread activation of immune cells is especially pronounced, and is first detectable in white adipose depots during middle age. Single-cell RNA sequencing confirms the accumulation of T cells and B cells in adipose tissue-including plasma cells that express immunoglobulin J-which also accrue concurrently across diverse organs. Finally, we show how gene expression shifts in distinct tissues are highly correlated with corresponding protein levels in plasma, thus potentially contributing to the ageing of the systemic circulation. Together, these data demonstrate a similar yet asynchronous inter- and intra-organ progression of ageing, providing a foundation from which to track systemic sources of declining health at old age.

Tuning MPL signaling to influence hematopoietic stem cell differentiation and inhibit essential thrombocythemia progenitors.

Thrombopoietin (TPO) and the TPO-receptor (TPO-R, or c-MPL) are essential for hematopoietic stem cell (HSC) maintenance and megakaryocyte differentiation. Agents that can modulate TPO-R signaling are highly desirable for both basic research and clinical utility. We developed a series of surrogate protein ligands for TPO-R, in the form of diabodies (DBs), that homodimerize TPO-R on the cell surface in geometries that are dictated by the DB receptor binding epitope, in effect "tuning" downstream signaling responses. These surrogate ligands exhibit diverse pharmacological properties, inducing graded signaling outputs, from full to partial TPO agonism, thus decoupling the dual functions of TPO/TPO-R. Using single-cell RNA sequencing and HSC self-renewal assays we find that partial agonistic diabodies preserved the stem-like properties of cultured HSCs, but also blocked oncogenic colony formation in essential thrombocythemia (ET) through inverse agonism. Our data suggest that dampening downstream TPO signaling is a powerful approach not only for HSC preservation in culture, but also for inhibiting oncogenic signaling through the TPO-R.

MARS: discovering novel cell types across heterogeneous single-cell experiments.

Although tremendous effort has been put into cell-type annotation, identification of previously uncharacterized cell types in heterogeneous single-cell RNA-seq data remains a challenge. Here we present MARS, a meta-learning approach for identifying and annotating known as well as new cell types. MARS overcomes the heterogeneity of cell types by transferring latent cell representations across multiple datasets. MARS uses deep learning to learn a cell embedding function as well as a set of landmarks in the cell embedding space. The method has a unique ability to discover cell types that have never been seen before and annotate experiments that are as yet unannotated. We apply MARS to a large mouse cell atlas and show its ability to accurately identify cell types, even when it has never seen them before. Further, MARS automatically generates interpretable names for new cell types by probabilistically defining a cell type in the embedding space.

PS1 FAD mutants decrease ephrinB2-regulated angiogenic functions, ischemia-induced brain neovascularization and neuronal survival.

Microvascular pathology and ischemic lesions contribute substantially to neuronal dysfunction and loss that lead to Alzheimer disease (AD). To facilitate recovery, the brain stimulates neovascularization of damaged tissue via sprouting angiogenesis, a process regulated by endothelial cell (EC) sprouting and the EphB4/ephrinB2 system. Here, we show that in cultures of brain ECs, EphB4 stimulates the VE-cadherin/Rok-α angiogenic complexes known to mediate sprouting angiogenesis. Importantly, brain EC cultures expressing PS1 FAD mutants decrease the EphB4-stimulated γ-secretase cleavage of ephrinB2 and reduce production of the angiogenic peptide ephrinB2/CTF2, the VE-cadherin angiogenic complexes and EC sprouting and tube formation. These data suggest that FAD mutants may attenuate ischemia-induced brain angiogenesis. Supporting this hypothesis, ischemia-induced VE-cadherin angiogenic complexes, levels of neoangiogenesis marker Endoglin, vascular density, and cerebral blood flow recovery, are all decreased in brains of mouse models expressing PS1 FAD mutants. Ischemia-induced brain neuronal death and cognitive deficits also increase in these mice. Furthermore, a small peptide comprising the C-terminal sequence of peptide ephrinB2/CTF2 rescues angiogenic functions of brain ECs expressing PS1 FAD mutants. Together, our data show that PS1 FAD mutations impede the EphB4/ephrinB2-mediated angiogenic functions of ECs and impair brain neovascularization, neuronal survival and cognitive recovery following ischemia. Furthermore, our data reveal a novel brain angiogenic mechanism targeted by PS1 FAD mutants and a potential therapeutic target for ischemia-induced neurodegeneration. Importantly, FAD mutant effects occur in absence of neuropathological hallmarks of AD, supporting that such hallmarks may form downstream of mutant effects on neoangiogenesis and neuronal survival.

Single-cell RNAseq reveals cell adhesion molecule profiles in electrophysiologically defined neurons.

In brain, signaling mediated by cell adhesion molecules defines the identity and functional properties of synapses. The specificity of presynaptic and postsynaptic interactions that is presumably mediated by cell adhesion molecules suggests that there exists a logic that could explain neuronal connectivity at the molecular level. Despite its importance, however, the nature of such logic is poorly understood, and even basic parameters, such as the number, identity, and single-cell expression profiles of candidate synaptic cell adhesion molecules, are not known. Here, we devised a comprehensive list of genes involved in cell adhesion, and used single-cell RNA sequencing (RNAseq) to analyze their expression in electrophysiologically defined interneurons and projection neurons. We compared the cell type-specific expression of these genes with that of genes involved in transmembrane ion conductances (i.e., channels), exocytosis, and rho/rac signaling, which regulates the actin cytoskeleton. Using these data, we identified two independent, developmentally regulated networks of interacting genes encoding molecules involved in cell adhesion, exocytosis, and signal transduction. Our approach provides a framework for a presumed cell adhesion and signaling code in neurons, enables correlating electrophysiological with molecular properties of neurons, and suggests avenues toward understanding synaptic specificity.

A survey of human brain transcriptome diversity at the single cell level.

The human brain is a tissue of vast complexity in terms of the cell types it comprises. Conventional approaches to classifying cell types in the human brain at single cell resolution have been limited to exploring relatively few markers and therefore have provided a limited molecular characterization of any given cell type. We used single cell RNA sequencing on 466 cells to capture the cellular complexity of the adult and fetal human brain at a whole transcriptome level. Healthy adult temporal lobe tissue was obtained during surgical procedures where otherwise normal tissue was removed to gain access to deeper hippocampal pathology in patients with medical refractory seizures. We were able to classify individual cells into all of the major neuronal, glial, and vascular cell types in the brain. We were able to divide neurons into individual communities and show that these communities preserve the categorization of interneuron subtypes that is typically observed with the use of classic interneuron markers. We then used single cell RNA sequencing on fetal human cortical neurons to identify genes that are differentially expressed between fetal and adult neurons and those genes that display an expression gradient that reflects the transition between replicating and quiescent fetal neuronal populations. Finally, we observed the expression of major histocompatibility complex type I genes in a subset of adult neurons, but not fetal neurons. The work presented here demonstrates the applicability of single cell RNA sequencing on the study of the adult human brain and constitutes a first step toward a comprehensive cellular atlas of the human brain.

Multiple recognition assay reveals prostasomes as promising plasma biomarkers for prostate cancer.

Prostasomes are microvesicles (mean diameter, 150 nm) that are produced and secreted by normal and malignant prostate acinar cells. It has been hypothesized that invasive growth of malignant prostate cells may cause these microvesicles, normally released into seminal fluid, to appear in interstitial space and therewith into peripheral circulation. The suitability of prostasomes as blood biomarkers in patients with prostate cancer was tested by using an expanded variant of the proximity ligation assay (PLA). We developed an extremely sensitive and specific assay (4PLA) for detection of complex target structures such as microvesicles in which the target is first captured via an immobilized antibody and subsequently detected by using four other antibodies with attached DNA strands. The requirement for coincident binding by five antibodies to generate an amplifiable reporter results in both increased specificity and sensitivity. The assay successfully detected significantly elevated levels of prostasomes in blood samples from patients with prostate cancer before radical prostatectomy, compared with controls and men with benign biopsy results. The medians for prostasome levels in blood plasma of patients with prostate cancer were 2.5 to sevenfold higher compared with control samples in two independent studies, and the assay also distinguished patients with high and medium prostatectomy Gleason scores (8/9 and 7, respectively) from those with low score (≤ 6), thus reflecting disease aggressiveness. This approach that enables detection of prostasomes in peripheral blood may be useful for early diagnosis and assessment of prognosis in organ-confined prostate cancer.

Coordinated wound responses in a regenerative animal-algal holobiont.

Animal regeneration involves coordinated responses across cell types throughout the animal body. In endosymbiotic animals, whether and how symbionts react to host injury and how cellular responses are integrated across species remain unexplored. Here, we study the acoel Convolutriloba longifissura, which hosts symbiotic Tetraselmis sp. green algae and can regenerate entire bodies from tissue fragments. We show that animal injury causes a decline in the photosynthetic efficiency of the symbiotic algae, alongside two distinct, sequential waves of transcriptional responses in acoel and algal cells. The initial algal response is characterized by the upregulation of a cohort of photosynthesis-related genes, though photosynthesis is not necessary for regeneration. A conserved animal transcription factor, runt, is induced after injury and required for acoel regeneration. Knockdown of Cl-runt dampens transcriptional responses in both species and further reduces algal photosynthetic efficiency post-injury. Our results suggest that the holobiont functions as an integrated unit of biological organization by coordinating molecular networks across species through the runt-dependent animal regeneration program.

IL-4-induced SOX9 confers lineage plasticity to aged adult lung stem cells.

Wound healing in response to acute injury is mediated by the coordinated and transient activation of parenchymal, stromal, and immune cells that resolves to homeostasis. Environmental, genetic, and epigenetic factors associated with inflammation and aging can lead to persistent activation of the microenvironment and fibrosis. Here, we identify opposing roles of interleukin-4 (IL-4) cytokine signaling in interstitial macrophages and type II alveolar epithelial cells (ATIIs). We show that IL4Ra signaling in macrophages promotes regeneration of the alveolar epithelium after bleomycin-induced lung injury. Using organoids and mouse models, we show that IL-4 directly acts on a subset of ATIIs to induce the expression of the transcription factor SOX9 and reprograms them toward a progenitor-like state with both airway and alveolar lineage potential. In the contexts of aging and bleomycin-induced lung injury, this leads to aberrant epithelial cell differentiation and bronchiolization, consistent with cellular and histological changes observed in interstitial lung disease.

Single-cell analysis of breast cancer metastasis reveals epithelial-mesenchymal plasticity signatures associated with poor outcomes.

Metastasis is the leading cause of cancer-related deaths. It is unclear how intratumor heterogeneity (ITH) contributes to metastasis and how metastatic cells adapt to distant tissue environments. The study of these adaptations is challenged by the limited access to patient material and a lack of experimental models that appropriately recapitulate ITH. To investigate metastatic cell adaptations and the contribution of ITH to metastasis, we analyzed single-cell transcriptomes of matched primary tumors and metastases from patient-derived xenograft models of breast cancer. We found profound transcriptional differences between the primary tumor and metastatic cells. Primary tumors upregulated several metabolic genes, whereas motility pathway genes were upregulated in micrometastases, and stress response signaling was upregulated during progression. Additionally, we identified primary tumor gene signatures that were associated with increased metastatic potential and correlated with patient outcomes. Immune-regulatory control pathways were enriched in poorly metastatic primary tumors, whereas genes involved in epithelial-mesenchymal transition were upregulated in highly metastatic tumors. We found that ITH was dominated by epithelial-mesenchymal plasticity (EMP), which presented as a dynamic continuum with intermediate EMP cell states characterized by specific genes such as CRYAB and S100A2. Elevated expression of an intermediate EMP signature correlated with worse patient outcomes. Our findings identified inhibition of the intermediate EMP cell state as a potential therapeutic target to block metastasis.

Transcriptional dynamics of delaminating neuroblasts in the mouse otic vesicle.

An abundance of research has recently highlighted the susceptibility of cochleovestibular ganglion (CVG) neurons to noise damage and aging in the adult cochlea, resulting in hearing deficits. Furthering our understanding of the transcriptional cascades that contribute to CVG development may provide insight into how these cells can be regenerated to treat inner ear dysfunction. Here we perform a high-depth single-cell RNA sequencing analysis of the E10.5 otic vesicle and its surrounding tissues, including CVG precursor neuroblasts and emerging CVG neurons. Clustering and trajectory analysis of otic-lineage cells reveals otic markers and the changes in gene expression that occur from neuroblast delamination toward the development of the CVG. This dataset provides a valuable resource for further identifying the mechanisms associated with CVG development from neurosensory competent cells within the otic vesicle.

Generation of Schwann cell derived melanocytes from hPSCs identifies pro-metastatic factors in melanoma.

The neural crest (NC) is highly multipotent and generates diverse lineages in the developing embryo. However, spatiotemporally distinct NC populations display differences in fate potential, such as increased gliogenic and parasympathetic potential from later migrating, nerve-associated Schwann cell precursors (SCPs). Interestingly, while melanogenic potential is shared by both early migrating NC and SCPs, differences in melanocyte identity resulting from differentiation through these temporally distinct progenitors have not been determined. Here, we leverage a human pluripotent stem cell (hPSC) model of NC temporal patterning to comprehensively characterize human NC heterogeneity, fate bias, and lineage development. We captured the transition of NC differentiation between temporally and transcriptionally distinct melanogenic progenitors and identified modules of candidate transcription factor and signaling activity associated with this transition. For the first time, we established a protocol for the directed differentiation of melanocytes from hPSCs through a SCP intermediate, termed trajectory 2 (T2) melanocytes. Leveraging an existing protocol for differentiating early NC-derived melanocytes, termed trajectory 1 (T1), we performed the first comprehensive comparison of transcriptional and functional differences between these distinct melanocyte populations, revealing differences in pigmentation and unique expression of transcription factors, ligands, receptors and surface markers. We found a significant link between the T2 melanocyte transcriptional signature and decreased survival in melanoma patients in the cancer genome atlas (TCGA). We performed an in vivo CRISPRi screen of T1 and T2 melanocyte signature genes in a human melanoma cell line and discovered several T2-specific markers that promote lung metastasis in mice. We further demonstrated that one of these factors, SNRPB, regulates the splicing of transcripts involved in metastasis relevant functions such as migration, cell adhesion and proliferation. Overall, this study identifies distinct developmental trajectories as a source of diversity in melanocytes and implicates the unique molecular signature of SCP-derived melanocytes in metastatic melanoma.

DNA-assisted protein detection technologies.

Improved protein assays promise to offer new insights into biological processes as well as the identification of new, clinically important biomarkers. In recent years, a number of approaches have been developed where protein-binding reagents, typically antibodies, are equipped with DNA strands to enable protein analyses via powerful nucleic acid detection reactions for improved performance. In this review, we provide a background to this emerging field, and we describe several different ways in which these reagents can improve protein analyses by lowering detection thresholds, improving multiplexing and extending the range of biomolecules available for analysis, both in research settings and in clinical routine.