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.

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.

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.

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.

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.

Neuroendocrinology of the lung revealed by single-cell RNA sequencing.

Pulmonary neuroendocrine cells (PNECs) are sensory epithelial cells that transmit airway status to the brain via sensory neurons and locally via calcitonin gene-related peptide (CGRP) and γ- aminobutyric acid (GABA). Several other neuropeptides and neurotransmitters have been detected in various species, but the number, targets, functions, and conservation of PNEC signals are largely unknown. We used scRNAseq to profile hundreds of the rare mouse and human PNECs. This revealed over 40 PNEC neuropeptide and peptide hormone genes, most cells expressing unique combinations of 5-18 genes. Peptides are packaged in separate vesicles, their release presumably regulated by the distinct, multimodal combinations of sensors we show are expressed by each PNEC. Expression of the peptide receptors predicts an array of local cell targets, and we show the new PNEC signal angiotensin directly activates one subtype of innervating sensory neuron. Many signals lack lung targets so may have endocrine activity like those of PNEC-derived carcinoid tumors. PNECs are an extraordinarily rich and diverse signaling hub rivaling the enteroendocrine system.

Population-wide gene disruption in the murine lung epithelium via AAV-mediated delivery of CRISPR-Cas9 components.

With the aim of expediting drug target discovery and validation for respiratory diseases, we developed an optimized method for in situ somatic gene disruption in murine lung epithelial cells via AAV6-mediated CRISPR-Cas9 delivery. Efficient gene editing was observed in lung type II alveolar epithelial cells and distal airway cells following assessment of single- or dual-guide AAV vector formats, Cas9 variants, and a sequential dosing strategy with combinatorial guide RNA expression cassettes. In particular, we were able to demonstrate population-wide gene disruption within distinct epithelial cell types for separate targets in Cas9 transgenic animals, with minimal to no associated inflammation. We also observed and characterized AAV vector integration events that occurred within directed double-stranded DNA break sites in lung cells, highlighting a complicating factor with AAV-mediated delivery of DNA nucleases. Taken together, we demonstrate a uniquely effective approach for somatic engineering of the murine lung, which will greatly facilitate the modeling of disease and therapeutic intervention.

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.

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.

Leveraging the Cell Ontology to classify unseen cell types.

Single cell technologies are rapidly generating large amounts of data that enables us to understand biological systems at single-cell resolution. However, joint analysis of datasets generated by independent labs remains challenging due to a lack of consistent terminology to describe cell types. Here, we present OnClass, an algorithm and accompanying software for automatically classifying cells into cell types that are part of the controlled vocabulary that forms the Cell Ontology. A key advantage of OnClass is its capability to classify cells into cell types not present in the training data because it uses the Cell Ontology graph to infer cell type relationships. Furthermore, OnClass can be used to identify marker genes for all the cell ontology categories, regardless of whether the cell types are present or absent in the training data, suggesting that OnClass goes beyond a simple annotation tool for single cell datasets, being the first algorithm capable to identify marker genes specific to all terms of the Cell Ontology and offering the possibility of refining the Cell Ontology using a data-centric approach.

Human melanocyte development and melanoma dedifferentiation at single-cell resolution.

In humans, epidermal melanocytes are responsible for skin pigmentation, defence against ultraviolet radiation and the deadliest common skin cancer, melanoma. Although there is substantial overlap in melanocyte development pathways between different model organisms, species-dependent differences are frequent and the conservation of these processes in human skin remains unresolved. Here, we used a single-cell enrichment and RNA-sequencing pipeline to study human epidermal melanocytes directly from the skin, capturing transcriptomes across different anatomical sites, developmental age, sexes and multiple skin tones. We uncovered subpopulations of melanocytes that exhibit anatomical site-specific enrichment that occurs during gestation and persists through adulthood. The transcriptional signature of the volar-enriched subpopulation is retained in acral melanomas. Furthermore, we identified human melanocyte differentiation transcriptional programs that are distinct from gene signatures generated from model systems. Finally, we used these programs to define patterns of dedifferentiation that are predictive of melanoma prognosis and response to immune checkpoint inhibitor therapy.

Impaired immune signaling and changes in the lung microbiome precede secondary bacterial pneumonia in COVID-19.

Secondary bacterial infections, including ventilator-associated pneumonia (VAP), lead to worse clinical outcomes and increased mortality following viral respiratory infections including in patients with coronavirus disease 2019 (COVID-19). Using a combination of tracheal aspirate bulk and single-cell RNA sequencing (scRNA-seq) we assessed lower respiratory tract immune responses and microbiome dynamics in 28 COVID-19 patients, 15 of whom developed VAP, and eight critically ill uninfected controls. Two days before VAP onset we observed a transcriptional signature of bacterial infection. Two weeks prior to VAP onset, following intubation, we observed a striking impairment in immune signaling in COVID-19 patients who developed VAP. Longitudinal metatranscriptomic analysis revealed disruption of lung microbiome community composition in patients with VAP, providing a connection between dysregulated immune signaling and outgrowth of opportunistic pathogens. These findings suggest that COVID-19 patients who develop VAP have impaired antibacterial immune defense detectable weeks before secondary infection onset.

Impaired immune signaling and changes in the lung microbiome precede secondary bacterial pneumonia in COVID-19.

Secondary bacterial infections, including ventilator-associated pneumonia (VAP), lead to worse clinical outcomes and increased mortality following viral respiratory infections including in patients with coronavirus disease 2019 (COVID-19). Using a combination of tracheal aspirate bulk and single-cell RNA sequencing we assessed lower respiratory tract immune responses and microbiome dynamics in 23 COVID-19 patients, 10 of whom developed VAP, and eight critically ill uninfected controls. At a median of three days (range: 2-4 days) before VAP onset we observed a transcriptional signature of bacterial infection. At a median of 15 days prior to VAP onset (range: 8-38 days), we observed a striking impairment in immune signaling in COVID-19 patients who developed VAP. Longitudinal metatranscriptomic analysis revealed disruption of lung microbiome community composition in patients with VAP, providing a connection between dysregulated immune signaling and outgrowth of opportunistic pathogens. These findings suggest that COVID-19 patients who develop VAP have impaired antibacterial immune defense detectable weeks before secondary infection onset.

Mouse aging cell atlas analysis reveals global and cell type-specific aging signatures.

Aging is associated with complex molecular and cellular processes that are poorly understood. Here we leveraged the Tabula Muris Senis single-cell RNA-seq data set to systematically characterize gene expression changes during aging across diverse cell types in the mouse. We identified aging-dependent genes in 76 tissue-cell types from 23 tissues and characterized both shared and tissue-cell-specific aging behaviors. We found that the aging-related genes shared by multiple tissue-cell types also change their expression congruently in the same direction during aging in most tissue-cell types, suggesting a coordinated global aging behavior at the organismal level. Scoring cells based on these shared aging genes allowed us to contrast the aging status of different tissues and cell types from a transcriptomic perspective. In addition, we identified genes that exhibit age-related expression changes specific to each functional category of tissue-cell types. Altogether, our analyses provide one of the most comprehensive and systematic characterizations of the molecular signatures of aging across diverse tissue-cell types in a mammalian system.

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.

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.

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.

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.

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.