We present an integrated analysis of the clinical measurements, immune cells, and plasma multi-omics of 139 COVID-19 patients representing all levels of disease severity, from serial blood draws collected during the first week of infection following diagnosis. We identify a major shift between mild and moderate disease, at which point elevated inflammatory signaling is accompanied by the loss of specific classes of metabolites and metabolic processes. Within this stressed plasma environment at moderate disease, multiple unusual immune cell phenotypes emerge and amplify with increasing disease severity. We condensed over 120,000 immune features into a single axis to capture how different immune cell classes coordinate in response to SARS-CoV-2. This immune-response axis independently aligns with the major plasma composition changes, with clinical metrics of blood clotting, and with the sharp transition between mild and moderate disease. This study suggests that moderate disease may provide the most effective setting for therapeutic intervention.

We describe a microchip designed to quantify the levels of a dozen cytoplasmic and membrane proteins from single cells. We use the platform to assess protein–protein interactions associated with the EGF-receptor-mediated PI3K signaling pathway. Single-cell sensitivity is achieved by isolating a defined number of cells ( n = 0–5) in 2 nL volume chambers, each of which is patterned with two copies of a miniature antibody array. The cells are lysed on-chip, and the levels of released proteins are assayed using the antibody arrays. We investigate three isogenic cell lines representing the cancer glioblastoma multiforme, at the basal level, under EGF stimulation, and under erlotinib inhibition plus EGF stimulation. The measured protein abundances are consistent with previous work, and single-cell analysis uniquely reveals single-cell heterogeneity, and different types and strengths of protein–protein interactions. This platform helps provide a comprehensive picture of altered signal transduction networks in tumor cells and provides insight into the effect of targeted therapies on protein signaling networks.

Significance: This work provides biophysical insights into how BRAF mutant melanoma cells adapt to the stress of MAPK inhibition via a series of reversible phenotypic transitions toward drug-tolerant or drug-resistant cell states enriched for neural-crest factors and mesenchymal signatures. This adaptation is influenced by cell phenotype-specific drug selection and cell state interconversion, but not selection of genetically resistant clones. A panel of functional proteins, analyzed at the single-cell level, pointed to signaling network hubs that drive the initiation of the melanoma cell adaptive transition. Targeting those hubs halted the transition and arrested resistance development.

The determination of individual cell trajectories through a high-dimensional cell-state space is an outstanding challenge for understanding biological changes ranging from cellular differentiation to epigenetic responses of diseased cells upon drugging. We integrate experiments and theory to determine the trajectories that single BRAF^V600E mutant melanoma cancer cells take between drug-naive and drug-tolerant states. Although single-cell omics tools can yield snapshots of the cell-state landscape, the determination of individual cell trajectories through that space can be confounded by stochastic cell-state switching. We assayed for a panel of signaling, phenotypic, and metabolic regulators at points across 5 days of drug treatment to uncover a cell-state landscape with two paths connecting drug-naive and drug-tolerant states. The trajectory a given cell takes depends upon the drug-naive level of a lineage-restricted transcription factor. Each trajectory exhibits unique druggable susceptibilities, thus updating the paradigm of adaptive resistance development in an isogenic cell population.