Cellular circuits sense the environment, process signals, and compute decisions using networks of interacting proteins. To model such a system, the abundance of each activated protein species can be described as a stochastic function of the abundance of other proteins. High-dimensional single-cell technologies, such as mass cytometry, offer an opportunity to characterize signaling circuit-wide. However, the challenge of developing and applying computational approaches to interpret such complex data remains. Here, we developed computational methods, based on established statistical concepts, to characterize signaling network relationships by quantifying the strengths of network edges and deriving signaling response functions. In comparing signaling between naïve and antigen-exposed CD4(+) T lymphocytes, we find that although these two cell subtypes had similarly wired networks, naïve cells transmitted more information along a key signaling cascade than did antigen-exposed cells. We validated our characterization on mice lacking the extracellular-regulated mitogen-activated protein kinase (MAPK) ERK2, which showed stronger influence of pERK on pS6 (phosphorylated-ribosomal protein S6), in naïve cells as compared with antigen-exposed cells, as predicted. We demonstrate that by using cell-to-cell variation inherent in single-cell data, we can derive response functions underlying molecular circuits and drive the understanding of how cells process signals.

Signaling from the T-cell receptor (TCR) conditions T-cell differentiation and activation, requiring exquisite sensitivity and discrimination. Using mass cytometry, a high-dimensional technique that can probe multiple signaling nodes at the single-cell level, we interrogate TCR signaling dynamics in control C57BL/6 and autoimmunity-prone nonobese diabetic (NOD) mice, which show ineffective ERK activation after TCR triggering. By quantitating signals at multiple steps along the signaling cascade and parsing the phosphorylation level of each node as a function of its predecessors, we show that a small impairment in initial pCD3ζ activation resonates farther down the signaling cascade and results in larger defects in activation of the ERK1/2-S6 and IκBα modules. This nonlinear property of TCR signaling networks, which magnifies small initial differences during signal propagation, also applies in cells from B6 mice activated at different levels of intensity. Impairment in pCD3ζ and pSLP76 is not a feedback consequence of a primary deficiency in ERK activation because no proximal signaling defect was observed in Erk2 KO T cells. These defects, which were manifest at all stages of T-cell differentiation from early thymic pre-T cells to memory T cells, may condition the imbalanced immunoregulation and tolerance in NOD T cells. More generally, this amplification of small initial differences in signal intensity may explain how T cells discriminate between closely related ligands and adopt strongly delineated cell fates.

Tumor responses to programmed cell death protein 1 (PD-1) blockade therapy are mediated by T cells, which we characterized in 102 tumor biopsies obtained from 53 patients treated with pembrolizumab, an antibody to PD-1. Biopsies were dissociated, and single-cell infiltrates were analyzed by multicolor flow cytometry using two computational approaches to resolve the leukocyte phenotypes at the single-cell level. There was a statistically significant increase in the frequency of T cells in patients who responded to therapy. The frequency of intratumoral B cells and monocytic myeloid-derived suppressor cells significantly increased in patients' biopsies taken on treatment. The percentage of cells with a regulatory T-cell phenotype, monocytes, and natural killer cells did not change while on PD-1 blockade therapy. CD8(+) memory T cells were the most prominent phenotype that expanded intratumorally on therapy. However, the frequency of CD4(+) effector memory T cells significantly decreased on treatment, whereas CD4(+) effector T cells significantly increased in nonresponding tumors on therapy. In peripheral blood, an unusual population of blood cells expressing CD56 was detected in two patients with regressing melanoma. In conclusion, PD-1 blockade increases the frequency of T cells, B cells, and myeloid-derived suppressor cells in tumors, with the CD8(+) effector memory T-cell subset being the major T-cell phenotype expanded in patients with a response to therapy.

Chromosomal instability is a hallmark of cancer that results from ongoing errors in chromosome segregation during mitosis. Although chromosomal instability is a major driver of tumour evolution, its role in metastasis has not been established. Here we show that chromosomal instability promotes metastasis by sustaining a tumour cell-autonomous response to cytosolic DNA. Errors in chromosome segregation create a preponderance of micronuclei whose rupture spills genomic DNA into the cytosol. This leads to the activation of the cGAS-STING (cyclic GMP-AMP synthase-stimulator of interferon genes) cytosolic DNA-sensing pathway and downstream noncanonical NF-κB signalling. Genetic suppression of chromosomal instability markedly delays metastasis even in highly aneuploid tumour models, whereas continuous chromosome segregation errors promote cellular invasion and metastasis in a STING-dependent manner. By subverting lethal epithelial responses to cytosolic DNA, chromosomally unstable tumour cells co-opt chronic activation of innate immune pathways to spread to distant organs.

The Human Developmental Cell Atlas (HDCA) initiative, which is part of the Human Cell Atlas, aims to create a comprehensive reference map of cells during development. This will be critical to understanding normal organogenesis, the effect of mutations, environmental factors and infectious agents on human development, congenital and childhood disorders, and the cellular basis of ageing, cancer and regenerative medicine. Here we outline the HDCA initiative and the challenges of mapping and modelling human development using state-of-the-art technologies to create a reference atlas across gestation. Similar to the Human Genome Project, the HDCA will integrate the output from a growing community of scientists who are mapping human development into a unified atlas. We describe the early milestones that have been achieved and the use of human stem-cell-derived cultures, organoids and animal models to inform the HDCA, especially for prenatal tissues that are hard to acquire. Finally, we provide a roadmap towards a complete atlas of human development.

Crucial transitions in cancer-including tumor initiation, local expansion, metastasis, and therapeutic resistance-involve complex interactions between cells within the dynamic tumor ecosystem. Transformative single-cell genomics technologies and spatial multiplex in situ methods now provide an opportunity to interrogate this complexity at unprecedented resolution. The Human Tumor Atlas Network (HTAN), part of the National Cancer Institute (NCI) Cancer Moonshot Initiative, will establish a clinical, experimental, computational, and organizational framework to generate informative and accessible three-dimensional atlases of cancer transitions for a diverse set of tumor types. This effort complements both ongoing efforts to map healthy organs and previous large-scale cancer genomics approaches focused on bulk sequencing at a single point in time. Generating single-cell, multiparametric, longitudinal atlases and integrating them with clinical outcomes should help identify novel predictive biomarkers and features as well as therapeutically relevant cell types, cell states, and cellular interactions across transitions. The resulting tumor atlases should have a profound impact on our understanding of cancer biology and have the potential to improve cancer detection, prevention, and therapeutic discovery for better precision-medicine treatments of cancer patients and those at risk for cancer.

The recent advent of methods for high-throughput single-cell molecular profiling has catalyzed a growing sense in the scientific community that the time is ripe to complete the 150-year-old effort to identify all cell types in the human body. The Human Cell Atlas Project is an international collaborative effort that aims to define all human cell types in terms of distinctive molecular profiles (such as gene expression profiles) and to connect this information with classical cellular descriptions (such as location and morphology). An open comprehensive reference map of the molecular state of cells in healthy human tissues would propel the systematic study of physiological states, developmental trajectories, regulatory circuitry and interactions of cells, and also provide a framework for understanding cellular dysregulation in human disease. Here we describe the idea, its potential utility, early proofs-of-concept, and some design considerations for the Human Cell Atlas, including a commitment to open data, code, and community.

Single-cell technologies have transformed our understanding of human tissues. Yet, studies typically capture only a limited number of donors and disagree on cell type definitions. Integrating many single-cell datasets can address these limitations of individual studies and capture the variability present in the population. Here we present the integrated Human Lung Cell Atlas (HLCA), combining 49 datasets of the human respiratory system into a single atlas spanning over 2.4 million cells from 486 individuals. The HLCA presents a consensus cell type re-annotation with matching marker genes, including annotations of rare and previously undescribed cell types. Leveraging the number and diversity of individuals in the HLCA, we identify gene modules that are associated with demographic covariates such as age, sex and body mass index, as well as gene modules changing expression along the proximal-to-distal axis of the bronchial tree. Mapping new data to the HLCA enables rapid data annotation and interpretation. Using the HLCA as a reference for the study of disease, we identify shared cell states across multiple lung diseases, including SPP1⁺ profibrotic monocyte-derived macrophages in COVID-19, pulmonary fibrosis and lung carcinoma. Overall, the HLCA serves as an example for the development and use of large-scale, cross-dataset organ atlases within the Human Cell Atlas.