Single-cell RNA-seq analysis identifies markers of resistance to targeted BRAF inhibitors in melanoma cell populations.

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 data sets also present additional challenges such as high transcript dropout rates, stochastic transcription events, and complex population substructures. Here, we present a single-cell RNA-seq analysis and klustering evaluation (SAKE), a robust method for scRNA-seq analysis that provides quantitative statistical metrics at each step of the analysis pipeline. Comparing SAKE to multiple single-cell analysis methods shows that most methods perform similarly across a wide range of 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 (BRAFi). Single-cell RNA-seq 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 analyses 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.

Single-cell profiling of environmental enteropathy reveals signatures of epithelial remodeling and immune activation.

Environmental enteropathy (EE) is a subclinical condition of the small intestine that is highly prevalent in low- and middle-income countries. It is thought to be a key contributing factor to childhood malnutrition, growth stunting, and diminished oral vaccine responses. Although EE has been shown to be the by-product of a recurrent enteric infection, its full pathophysiology remains unclear. Here, we mapped the cellular and molecular correlates of EE by performing high-throughput, single-cell RNA-sequencing on 33 small intestinal biopsies from 11 adults with EE in Lusaka, Zambia (eight HIV-negative and three HIV-positive), six adults without EE in Boston, United States, and two adults in Durban, South Africa, which we complemented with published data from three additional individuals from the same clinical site. We analyzed previously defined bulk-transcriptomic signatures of reduced villus height and decreased microbial translocation in EE and showed that these signatures may be driven by an increased abundance of surface mucosal cells-a gastric-like subset previously implicated in epithelial repair in the gastrointestinal tract. In addition, we determined cell subsets whose fractional abundances associate with EE severity, small intestinal region, and HIV infection. Furthermore, by comparing duodenal EE samples with those from three control cohorts, we identified dysregulated WNT and MAPK signaling in the EE epithelium and increased proinflammatory cytokine gene expression in a T cell subset highly expressing a transcriptional signature of tissue-resident memory cells in the EE cohort. Together, our work elucidates epithelial and immune correlates of EE and nominates cellular and molecular targets for intervention.

Integrated single-cell analysis of multicellular immune dynamics during hyperacute HIV-1 infection.

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.

Seq-Well: A Sample-Efficient, Portable Picowell Platform for Massively Parallel Single-Cell RNA Sequencing.

Seq-Well is a low-cost picowell platform that can be used to simultaneously profile the transcriptomes of thousands of cells from diverse, low input clinical samples. In Seq-Well, uniquely barcoded mRNA capture beads and cells are co-confined in picowells that are sealed using a semipermeable membrane, enabling efficient cell lysis and mRNA capture. The beads are subsequently removed and processed in parallel for sequencing, with each transcript's cell of origin determined via the unique barcodes. Due to its simplicity and portability, Seq-Well can be performed almost anywhere.

Induction of metabolic quiescence defines the transitional to follicular B cell switch.

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.

High-Frequency, Functional HIV-Specific T-Follicular Helper and Regulatory Cells Are Present Within Germinal Centers in Children but Not Adults.

Broadly neutralizing antibodies (bnAbs) against HIV-1 are an effective means of preventing transmission. To better understand the mechanisms by which HIV-specific bnAbs naturally develop, we investigated blood and lymphoid tissue in pediatric infection, since potent bnAbs develop with greater frequency in children than adults. As in adults, the frequency of circulating effector T-follicular helper cells (TFH) in HIV infected, treatment naïve children correlates with neutralization breadth. However, major differences between children and adults were also observed both in circulation, and in a small number of tonsil samples. In children, TFH cells are significantly more abundant, both in blood and in lymphoid tissue germinal centers, than in adults. Second, HIV-specific TFH cells are more frequent in pediatric than in adult lymphoid tissue and secrete the signature cytokine IL-21, which HIV-infected adults do not. Third, the enrichment of IL-21-secreting HIV-specific TFH in pediatric lymphoid tissue is accompanied by increased TFH regulation via more abundant regulatory follicular T-cells and HIV-specific CXCR5+ CD8 T-cells compared to adults. The relationship between regulation and neutralization breadth is also observed in the pediatric PBMC samples and correlates with neutralization breadth. Matching neutralization data from lymphoid tissue samples is not available. However, the distinction between infected children and adults in the magnitude, quality and regulation of HIV-specific TFH responses is consistent with the superior ability of children to develop high-frequency, potent bnAbs. These findings suggest the possibility that the optimal timing for next generation vaccine strategies designed to induce high-frequency, potent bnAbs to prevent HIV infection in adults would be in childhood.