Search

E3 ligases regulate key processes, but many of their roles remain unknown. Using Perturb-seq, we interrogated the function of 1,130 E3 ligases, partners and substrates in the inflammatory response in primary dendritic cells (DCs). Dozens impacted the balance of DC1, DC2, migratory DC and macrophage states and a gradient of DC maturation. Family members grouped into co-functional modules that were enriched for physical interactions and impacted specific programs through substrate transcription factors. E3s and their adaptors co-regulated the same processes, but partnered with different substrate recognition adaptors to impact distinct aspects of the DC life cycle. Genetic interactions were more prevalent within than between modules, and a deep learning model, comβVAE, predicts the outcome of new combinations by leveraging modularity. The E3 regulatory network was associated with heritable variation and aberrant gene expression in immune cells in human inflammatory diseases. Our study provides a general approach to dissect gene function.

Pooled CRISPR screens with single-cell RNA-seq readout (Perturb-seq) have emerged as a key technique in functional genomics, but are limited in scale by cost and combinatorial complexity. Here, we reimagine Perturb-seq’s design through the lens of algorithms applied to random, low-dimensional observations. We present compressed Perturb-seq, which measures multiple random perturbations per cell or multiple cells per droplet and computationally decompresses these measurements by leveraging the sparse structure of regulatory circuits. Applied to 598 genes in the immune response to bacterial lipopolysaccharide, compressed Perturb-seq achieves the same accuracy as conventional Perturb-seq at 4 to 20-fold reduced cost, with greater power to learn genetic interactions. We identify known and novel regulators of immune responses and uncover evolutionarily constrained genes with downstream targets enriched for immune disease heritability, including many missed by existing GWAS or trans-eQTL studies. Our framework enables new scales of interrogation for a foundational method in functional genomics.

The balance between T helper type 1 (Th1) cells and other Th cells is critical for antiviral and anti-tumor responses, but how this fine balance is achieved remains poorly understood. Here, we dissected the dynamic regulation of Th1 cell differentiation during in vitro polarization, as well as in vivo differentiation upon acute viral infection, using scRNA-seq and multiple in vitro and in vivo CRISPR screens. We confirmed the role of known regulators and identified 5 novel regulators for Th1 differentiation. Among the novel regulators the neuropeptide receptor RAMP3, which is a component of the receptor for calcitonin gene-related peptide (CGRP), plays a cell-intrinsic role in Th1 cell-fate determination. Using a unique Th1/Th2 dichotomous culture system, we identified that RAMP3-CGRP interaction directly restricted the differentiation of Th2 cells but promoted Th1 differentiation through activation of downstream cyclic AMP (cAMP) signaling in T cells. Mechanistically, RAMP3 and cAMP signaling resulted in global changes in chromatin accessibility, blocking Th2 genes and specific induction of Th1 programs through activation of IFNγ-STAT1 pathway. Furthermore, both CGRP and RAMP3 were required for inducing effective anti-viral T cell responses to control acute viral infection. Our work reveals a novel neuro-immune circuit in which tissue itself participates in T cell fate determination by producing a neuropeptide during acute viral infection, which acts on RAMP3 expressing T cells to induce an effective anti-viral Th1 response.

Here we introduce a mostly natural sequencing-by-synthesis (mnSBS) method for single-cell RNA sequencing (scRNA-seq), adapted to the Ultima genomics platform, and systematically benchmark it against current scRNA-seq technology. mnSBS uses mostly natural, unmodified nucleotides and only a low fraction of fluorescently labeled nucleotides, which allows for high polymerase processivity and lower costs. We demonstrate successful application in four scRNA-seq case studies of different technical and biological types, including 5′ and 3′ scRNA-seq, human peripheral blood mononuclear cells from a single individual and in multiplex, as well as Perturb-Seq. Benchmarking shows that results from mnSBS-based scRNA-seq are very similar to those using Illumina sequencing, with minor differences in results related to the position of reads relative to annotated gene boundaries, owing to single-end reads of Ultima being closer to gene ends than reads from Illumina. The method is thus compatible with state-of-the-art scRNA-seq libraries independent of the sequencing technology. We expect mnSBS to be of particular utility for cost-effective large-scale scRNA-seq projects.

A collection of programs is presented to analyze the thermodynamics of folding of linear repeat proteins using a 1D Ising model to determine intrinsic folding and interfacial coupling free energies. Expressions for folding transitions are generated for a series of constructs with different repeat numbers and are globally fitted to transitions for these constructs. These programs are designed to analyze Ising parameters for capped homopolymeric consensus repeat constructs as well as heteropolymeric constructs that contain point substitutions, providing a rigorous framework for analysis of the effects of mutation on intrinsic and directional (i.e., N- versus C-terminal) interfacial coupling free-energies. A bootstrap analysis is provided to estimate parameter uncertainty as well as correlations among fitted parameters. Rigorous statistical analysis is essential for interpreting fits using the complex models required for Ising analysis of repeat proteins, especially heteropolymeric repeat proteins. Programs described here are available athttps://github.com/barricklab-at-jhu/Ising_programs.

A collection of programs is presented to analyze the thermodynamics of folding of linear repeat proteins using a 1D Ising model to determine intrinsic folding and interfacial coupling free energies. Expressions for folding transitions are generated for a series of constructs with different repeat numbers and are globally fitted to transitions for these constructs. These programs are designed to analyze Ising parameters for capped homopolymeric consensus repeat constructs as well as heteropolymeric constructs that contain point substitutions, providing a rigorous framework for analysis of the effects of mutation on intrinsic and directional (i.e., N- vs. C-terminal) interfacial coupling free-energies. A bootstrap analysis is provided to estimate parameter uncertainty as well as correlations among fitted parameters. Rigorous statistical analysis is essential for interpreting fits using the complex models required for Ising analysis of repeat proteins, especially heteropolymeric repeat proteins. Programs described here are available at https://github.com/barricklab-at-jhu/Ising_programs.

Assays to study cancer cell responses to pharmacologic or genetic perturbations are typically restricted to using simple phenotypic readouts such as proliferation rate or the expression of a marker gene. Information-rich assays, such as gene-expression profiling, are generally not amenable to efficient profiling of a given perturbation across multiple cellular contexts. Here, we developed MIX-Seq, a method for multiplexed transcriptional profiling of post-perturbation responses across a mixture of samples with single-cell resolution, using SNP-based computational demultiplexing of single-cell RNA-sequencing data. We show that MIX-Seq can be used to profile responses to chemical or genetic perturbations across pools of 100 or more cancer cell lines, and combine it with Cell Hashing to further multiplex additional experimental conditions, such as multiple post-treatment time points or drug doses. Analyzing the high-content readout of scRNA-seq reveals both shared and context-specific transcriptional response components that can identify drug mechanism of action and can be used to predict long-term cell viability from short-term transcriptional responses to treatment.

Transcription activator-like effector proteins (TALEs) contain large numbers of repeats that bind double-stranded DNA, wrapping around DNA to form a continuous superhelix. Since unbound TALEs retain superhelical structure, it seems likely that DNA binding requires a significant structural distortion or partial unfolding. In this study, we use nearest-neighbor “Ising” analysis of consensus TALE (cTALE) repeat unfolding to quantify intrinsic folding free energies, coupling energies between repeats, and the free energy distribution of partly unfolded states, and to determine how those energies depend on the sequence that determines DNA-specificity (called the “RVD”). We find a moderate level of cooperativity for both the HD and NS RVD sequences (stabilizing interfaces combined with unstable repeats), as has been seen in other linear repeat proteins. Surprisingly, RVD sequence identity influences both the overall stability and the balance of intrinsic repeat stability and interfacial coupling energy. Using parameters from the Ising analysis, we have analyzed the distribution of partly folded states as a function of cTALE length and RVD sequence. We find partly unfolded states where one or more repeats are unfolded to be energetically accessible. Mixing repeats with different RVD sequences increases the population of partially folded states. Local folding free energies plateau for central repeats, suggesting that TALEs access partially folded states where a single internal repeat is unfolded while adjacent repeats remain folded. This breakage should allow TALEs to access superhelically-broken states, and may facilitate DNA binding.

Designed helical repeats (DHRs) are modular helix–loop–helix–loop protein structures that are tandemly repeated to form a superhelical array. Structures combining tandem DHRs demonstrate a wide range of molecular geometries, many of which are not observed in nature. Understanding cooperativity of DHR proteins provides insight into the molecular origins of Rosetta-based protein design hyperstability and facilitates comparison of energy distributions in artificial and naturally occurring protein folds. Here, we use a nearest-neighbor Ising model to quantify the intrinsic and interfacial free energies of four different DHRs. We measure the folding free energies of constructs with varying numbers of internal and terminal capping repeats for four different DHR folds, using guanidine-HCl and glycerol as destabilizing and solubilizing cosolvents. One-dimensional Ising analysis of these series reveals that, although interrepeat coupling energies are within the range seen for naturally occurring repeat proteins, the individual repeats of DHR proteins are intrinsically stable. This favorable intrinsic stability, which has not been observed for naturally occurring repeat proteins, adds to stabilizing interfaces, resulting in extraordinarily high stability. Stable repeats also impart a downhill shape to the energy landscape for DHR folding. These intrinsic stability differences suggest that part of the success of Rosetta-based design results from capturing favorable local interactions.

Rad6 is a yeast E2 ubiquitin conjugating enzyme that monoubiquitinates histone H2B in conjunction with the E3, Bre1, but can non-specifically modify histones on its own. We determined the crystal structure of a Rad6∼Ub thioester mimic, which revealed a network of interactions in the crystal in which the ubiquitin in one conjugate contacts Rad6 in another. The region of Rad6 contacted is located on the distal face of Rad6 opposite the active site, but differs from the canonical E2 backside that mediates free ubiquitin binding and polyubiquitination activity in other E2 enzymes. We find that free ubiquitin interacts weakly with both non-canonical and canonical backside residues of Rad6 and that mutations of non-canonical residues have deleterious effects on Rad6 activity comparable to those observed to mutations in the canonical E2 backside. The effect of non-canonical backside mutations is similar in the presence and absence of Bre1, indicating that contacts with non-canonical backside residues govern the intrinsic activity of Rad6. Our findings shed light on the determinants of intrinsic Rad6 activity and reveal new ways in which contacts with an E2 backside can regulate ubiquitin conjugating activity.