Epigenetic modulators of B cell fate identified through coupled phenotype-transcriptome analysis

High-throughput methodologies are the cornerstone of screening approaches to identify novel compounds that regulate immune cell function. To identify novel targeted therapeutics to treat immune disorders and haematological malignancies, there is a need to integrate functional cellular information with the molecular mechanisms that regulate changes in immune cell phenotype. We facilitate this goal by combining quantitative methods for dissecting complex simultaneous cell phenotypic effects with genomic analysis. This combination strategy we term Multiplexed Analysis of Cells sequencing (MAC-seq), a modified version of Digital RNA with perturbation of Genes (DRUGseq). We applied MAC-seq to screen compounds that target the epigenetic machinery of B cells and assess altered humoral immunity by measuring changes in proliferation, survival, differentiation and transcription. This approach revealed that polycomb repressive complex 2 (PRC2) inhibitors promote antibody secreting cell (ASC) differentiation in both murine and human B cells in vitro. This is further validated using T cell-dependent immunization in mice. Functional dissection of downstream effectors of PRC2 using arrayed CRISPR screening uncovered novel regulators of B cell differentiation, including Mybl1, Myof, Gas7 and Atoh8. Together, our findings demonstrate that integrated phenotype-transcriptome analyses can be effectively combined with drug screening approaches to uncover the molecular circuitry that drives lymphocyte fate decisions. Overall design: To examine the effect of EMCs on B cells, DRUGseq was performed on purified naive murine B cells stimulated with LPS (15ug/mL) with 1µM of various EMCs. To examine the changes in transcriptome and epigenome as a result of GSK126 treatment, purified B cells were stimulated with LPS (15ug/mL) and cultured with various concentrations of GSK126. The samples are harvested at indicated timepoints for RNA-seq and ChIP-seq (H3K27me3).