Cellular immunity is critical for controlling intracellular pathogens, but individual cellular dynamics and cell–cell cooperativity in evolving human immune responses remain poorly understood. Single-cell RNA-sequencing (scRNA-seq) represents a powerful tool for dissecting complex multicellular behaviors in health and disease1,2 and nominating testable therapeutic targets3. Its application to longitudinal samples could afford an opportunity to uncover cellular factors associated with the evolution of disease progression without potentially confounding inter-individual variability4. Here, we present an experimental and computational methodology that uses scRNA-seq to characterize dynamic cellular programs and their molecular drivers, and apply it to HIV infection. By performing scRNA-seq on peripheral blood mononuclear cells from four untreated individuals before and longitudinally during acute infection5, we were powered within each to discover gene response modules that vary by time and cell subset. Beyond previously unappreciated individual- and cell-type-specific interferon-stimulated gene upregulation, we describe temporally aligned gene expression responses obscured in bulk analyses, including those involved in proinflammatory T cell differentiation, prolonged monocyte major histocompatibility complex II upregulation and persistent natural killer (NK) cell cytolytic killing. We further identify response features arising in the first weeks of infection, for example proliferating natural killer cells, which potentially may associate with future viral control. Overall, our approach provides a unified framework for characterizing multiple dynamic cellular responses and their coordination.

Single-cell RNA-seq's (scRNA-seq) unprecedented cellular resolution at a genome wide scale enables us to address questions about cellular heterogeneity that are inaccessible using methods that average over bulk tissue extracts.However, scRNA-seq datasets also present additional challenges such as high transcript dropout rates, stochastic transcription events, and complex population substructures. Here, we present SAKE (Single-cell RNA-seq Analysis and Klustering Evaluation): a robust method for scRNA-seq analysis that provides quantitative statistical metrics at each step of the scRNA-seq analysis pipeline.Comparing SAKE to multiple single-cell analysis methods shows that most methods perform similarly across a wide range cellular contexts, with SAKE outperforming these methods in the case of large complex populations. We next applied the SAKE algorithms to identify drug-resistant cellular populations as human melanoma cells respond to targeted BRAF inhibitors. Single-cell RNA-seq Cold Spring Harbor Laboratory Press on June 7, 2019 -Published by genome.cshlp.org Downloaded from data from both the Fluidigm C1 and 10x Genomics platforms were analyzed with SAKE to dissect this problem at multiple scales. Data from both platforms indicate that BRAF inhibitor resistant cells can emerge from rare populations already present before drug application, with SAKE identifying both novel and known markers of resistance. These experimentally validated markers of BRAFi resistance share overlap with previous analysis in different melanoma cell lines, demonstrating the generality of these findings and highlighting the utility of single-cell analysis to elucidate mechanisms of BRAFi resistance.

Transitional B cells must actively undergo selection for self-tolerance before maturing into their resting follicular B cell successors. We found that metabolic quiescence was acquired at the follicular B cell stage in both humans and mice. In follicular B cells, the expression of genes involved in ribosome biogenesis, aerobic respiration, and mammalian target of rapamycin complex 1 (mTORC1) signaling was reduced when compared to that in transitional B cells. Functional metabolism studies, profiling of whole-cell metabolites, and analysis of cell surface proteins in human B cells suggested that this transition was also associated with increased extracellular adenosine salvage. Follicular B cells increased the abundance of the cell surface ectonucleotidase CD73, which coincided with adenosine 5′-monophosphate–activated protein kinase (AMPK) activation. Differentiation to the follicular B cell stage in vitro correlated with surface acquisition of CD73 on human transitional B cells and was augmented with the AMPK agonist, AICAR. Last, individuals with gain-of-function PIK3CD (PI3Kδ) mutations and increased pS6 activation exhibited a near absence of circulating follicular B cells. Together, our data suggest that mTORC1 attenuation may be necessary for human follicular B cell development. These data identify a distinct metabolic switch during human B cell development at the transitional to follicular stages, which is characterized by an induction of extracellular adenosine salvage, AMPK activation, and the acquisition of metabolic quiescence.