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.

Continuous BRAF inhibition of BRAF mutant melanomas triggers a series of cell state changes that lead to therapy resistance and escape from immune control before establishing acquired resistance genetically. We used genome-wide transcriptomics and single-cell phenotyping to explore the response kinetics to BRAF inhibition for a panel of patient-derived BRAFV600-mutant melanoma cell lines. A subset of plastic cell lines, which followed a trajectory covering multiple known cell state transitions, provided models for more detailed biophysical investigations. Markov modeling revealed that the cell state transitions were reversible and mediated by both Lamarckian induction and nongenetic Darwinian selection of drug-tolerant states. Single-cell functional proteomics revealed activation of certain signaling networks shortly after BRAF inhibition, and before the appearance of drug-resistant phenotypes. Drug targeting those networks, in combination with BRAF inhibition, halted the adaptive transition and led to prolonged growth inhibition in multiple patient-derived cell lines.

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 drugtolerant 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.

Phenotypic plasticity is associated with non-genetic drug tolerance in several cancers. Such plasticity can arise from chromatin remodeling, transcriptomic reprogramming, and/or protein signaling rewiring, and is characterized as a cell state transition in response to molecular or physical perturbations. This, in turn, can confound interpretations of drug responses and resistance development. Using BRAF- mutant melanoma cell lines as the prototype, we report on a joint theoretical and experimental investigation of the cell-state transition dynamics associated with BRAF inhibitor drug tolerance. Thermodynamically motivated surprisal analysis of transcriptome data was used to treat the cell population as an entropy maximizing system under the influence of time-dependent constraints. This permits the extraction of an epigenetic potential landscape for drug-induced phenotypic evolution. Single-cell flow cytometry data of the same system were modeled with a modified Fokker-Planck-type kinetic model. The two approaches yield a consistent picture that accounts for the phenotypic heterogeneity observed over the course of drug tolerance development. The results reveal that, in certain plastic cancers, the population heterogeneity and evolution of cell phenotypes may be understood by accounting for the competing interactions of the epigenetic potential landscape and state-dependent cell proliferation. Accounting for such competition permits accurate, experimentally verifiable predictions that can potentially guide the design of effective treatment strategies.

Accurate prediction of chemo- or targeted therapy responses for patients with similar driver oncogenes through a simple and least-invasive assay represents an unmet need in the clinical diagnosis of non-small cell lung cancer. Using a single-cell on-chip metabolic cytometry and fluorescent metabolic probes, we show metabolic phenotyping on the rare disseminated tumor cells in pleural effusions across a panel of 32 lung adenocarcinoma patients. Our results reveal extensive metabolic heterogeneity of tumor cells that differentially engage in glycolysis and mitochondrial oxidation. The cell number ratio of the two metabolic phenotypes is found to be predictive for patient therapy response, physiological performance, and survival. Transcriptome analysis reveals that the glycolytic phenotype is associated with mesenchymal-like cell state with elevated expression of the resistant-leading receptor tyrosine kinase AXL and immune checkpoint ligands. Drug targeting AXL induces a significant cell killing in the glycolytic cells without affecting the cells with active mitochondrial oxidation.

We present a chemical approach to profile fatty acid uptake in single cells. We use azide-modified analogues to probe the fatty acid influx and surface-immobilized dendrimers with dibenzocyclooctyne (DBCO) groups for detection. A competition between the fatty acid probes and BHQ2-azide quencher molecules generates fluorescence signals in a concentration-dependent manner. By integrating this method onto a microfluidics-based multiplex protein analysis platform, we resolved the relationships between fatty acid influx, oncogenic signaling activities, and cell proliferation in single glioblastoma cells. We found that p70S6K and 4EBP1 differentially correlated with fatty acid uptake. We validated that cotargeting p70S6K and fatty acid metabolism synergistically inhibited cell proliferation. Our work provided the first example of studying fatty acid metabolism in the context of protein signaling at single-cell resolution and generated new insights into cancer biology.

The immunological picture of how different patients recover from COVID-19, and how those recovery trajectories are influenced by infection severity, remain unclear. We investigated 140 COVID-19 patients from diagnosis to convalescence using clinical data, viral load assessments, and multi-omic analyses of blood plasma and circulating immune cells. Immune-phenotype dynamics resolved four recovery trajectories. One trajectory signals a return to pre-infection healthy baseline, while the other three are characterized by differing fractions of persistent cytotoxic and proliferative T cells, distinct B cell maturation processes, and memory-like innate immunity. We resolve a small panel of plasma proteins that, when measured at diagnosis, can predict patient survival and recovery-trajectory commitment. Our study offers novel insights into post-acute immunological outcomes of COVID-19 that likely influence long-term adverse sequelae.

No abstract