Understanding the effects of genetic perturbations on the cellular state has been challenging using traditional pooled screens, which typically rely on the delivery of a single perturbation per cell and unidimensional phenotypic readouts. Here, we use barcoded open reading frame overexpression libraries coupled with single-cell RNA sequencing to assay cell state and fitness, a technique we call SEUSS (scalable functional screening by sequencing). Using SEUSS, we perturbed hPSCs with a library of developmentally critical transcription factors (TFs) and assayed the impact of TF overexpression on fitness and transcriptomic states. We further leveraged the versatility of the ORF library approach to assay mutant genes and whole gene families. From the transcriptomic responses, we built genetic co-regulatory networks to identify altered gene modules and found that KLF4 and SNAI2 drive opposing effects along the epithelial-mesenchymal transition axis. From the fitness responses, we identified ETV2 as a driver of reprogramming toward an endothelial-like state.

The metabolic pathways fueling tumor growth have been well characterized, but the specific impact of transforming events on network topology and enzyme essentiality remains poorly understood. To this end, we performed combinatorial CRISPR-Cas9 screens on a set of 51 carbohydrate metabolism genes that represent glycolysis and the pentose phosphate pathway (PPP). This high-throughput methodology enabled systems-level interrogation of metabolic gene dispensability, interactions, and compensation across multiple cell types. The metabolic impact of specific combinatorial knockouts was validated using ¹³C and ²H isotope tracing, and these assays together revealed key nodes controlling redox homeostasis along the KEAP-NRF2 signaling axis. Specifically, targeting KEAP1 in combination with oxidative PPP genes mitigated the deleterious effects of these knockouts on growth rates. These results demonstrate how our integrated framework, combining genetic, transcriptomic, and flux measurements, can improve elucidation of metabolic network alterations and guide precision targeting of metabolic vulnerabilities based on tumor genetics.

We developed a systematic approach to map human genetic networks by combinatorial CRISPR-Cas9 perturbations coupled to robust analysis of growth kinetics. We targeted all pairs of 73 cancer genes with dual guide RNAs in three cell lines, comprising 141,912 tests of interaction. Numerous therapeutically relevant interactions were identified, and these patterns replicated with combinatorial drugs at 75% precision. From these results, we anticipate that cellular context will be critical to synthetic-lethal therapies.