Synthetic Lethal and Resistance Interactions with BET Bromodomain Inhibitors in Triple-Negative Breast Cancer.

BET bromodomain inhibitors (BBDIs) are candidate therapeutic agents for triple-negative breast cancer (TNBC) and other cancer types, but inherent and acquired resistance to BBDIs limits their potential clinical use. Using CRISPR and small-molecule inhibitor screens combined with comprehensive molecular profiling of BBDI response and resistance, we identified synthetic lethal interactions with BBDIs and genes that, when deleted, confer resistance. We observed synergy with regulators of cell cycle progression, YAP, AXL, and SRC signaling, and chemotherapeutic agents. We also uncovered functional similarities and differences among BRD2, BRD4, and BRD7. Although deletion of BRD2 enhances sensitivity to BBDIs, BRD7 loss leads to gain of TEAD-YAP chromatin binding and luminal features associated with BBDI resistance. Single-cell RNA-seq, ATAC-seq, and cellular barcoding analysis of BBDI responses in sensitive and resistant cell lines highlight significant heterogeneity among samples and demonstrate that BBDI resistance can be pre-existing or acquired.

Epigenetic silencing by SETDB1 suppresses tumour intrinsic immunogenicity.

Epigenetic dysregulation is a defining feature of tumorigenesis that is implicated in immune escape1,2. Here, to identify factors that modulate the immune sensitivity of cancer cells, we performed in vivo CRISPR-Cas9 screens targeting 936 chromatin regulators in mouse tumour models treated with immune checkpoint blockade. We identified the H3K9 methyltransferase SETDB1 and other members of the HUSH and KAP1 complexes as mediators of immune escape3-5. We also found that amplification of SETDB1 (1q21.3) in human tumours is associated with immune exclusion and resistance to immune checkpoint blockade. SETDB1 represses broad domains, primarily within the open genome compartment. These domains are enriched for transposable elements (TEs) and immune clusters associated with segmental duplication events, a central mechanism of genome evolution6. SETDB1 loss derepresses latent TE-derived regulatory elements, immunostimulatory genes, and TE-encoded retroviral antigens in these regions, and triggers TE-specific cytotoxic T cell responses in vivo. Our study establishes SETDB1 as an epigenetic checkpoint that suppresses tumour-intrinsic immunogenicity, and thus represents a candidate target for immunotherapy.

Avant-garde: an automated data-driven DIA data curation tool.

Several challenges remain in data-independent acquisition (DIA) data analysis, such as to confidently identify peptides, define integration boundaries, remove interferences, and control false discovery rates. In practice, a visual inspection of the signals is still required, which is impractical with large datasets. We present Avant-garde as a tool to refine DIA (and parallel reaction monitoring) data. Avant-garde uses a novel data-driven scoring strategy: signals are refined by learning from the dataset itself, using all measurements in all samples to achieve the best optimization. We evaluate the performance of Avant-garde using benchmark DIA datasets and show that it can determine the quantitative suitability of a peptide peak, and reach the same levels of selectivity, accuracy, and reproducibility as manual validation. Avant-garde is complementary to existing DIA analysis engines and aims to establish a strong foundation for subsequent analysis of quantitative mass spectrometry data.

Developmental control of polycomb subunit composition by GATA factors mediates a switch to non-canonical functions.

Polycomb repressive complex 2 (PRC2) plays crucial roles in transcriptional regulation and stem cell development. However, the context-specific functions associated with alternative subunits remain largely unexplored. Here we show that the related enzymatic subunits EZH1 and EZH2 undergo an expression switch during blood cell development. An erythroid-specific enhancer mediates transcriptional activation of EZH1, and a switch from GATA2 to GATA1 controls the developmental EZH1/2 switch by differential association with EZH1 enhancers. We further examine the in vivo stoichiometry of the PRC2 complexes by quantitative proteomics and reveal the existence of an EZH1-SUZ12 subcomplex lacking EED. EZH1 together with SUZ12 form a non-canonical PRC2 complex, occupy active chromatin, and positively regulate gene expression. Loss of EZH2 expression leads to repositioning of EZH1 to EZH2 targets. Thus, the lineage- and developmental stage-specific regulation of PRC2 subunit composition leads to a switch from canonical silencing to non-canonical functions during blood stem cell specification.

Study of mitotic chromatin supports a model of bookmarking by histone modifications and reveals nucleosome deposition patterns.

Mitosis encompasses key molecular changes including chromatin condensation, nuclear envelope breakdown, and reduced transcription levels. Immediately after mitosis, the interphase chromatin structure is reestablished and transcription resumes. The reestablishment of the interphase chromatin is probably achieved by "bookmarking," i.e., the retention of at least partial information during mitosis. To gain a deeper understanding of the contribution of histone modifications to the mitotic bookmarking process, we merged proteomics, immunofluorescence, and ChIP-seq approaches. We focused on key histone modifications and employed HeLa-S3 cells as a model system. Generally, in spite of the general hypoacetylation observed during mitosis, we observed a global concordance between the genomic organization of histone modifications in interphase and mitosis, suggesting that the epigenomic landscape may serve as a component of the mitotic bookmarking process. Next, we investigated the nucleosome that enters nucleosome depleted regions (NDRs) during mitosis. We observed that in ∼60% of the NDRs, the entering nucleosome is distinct from the surrounding highly acetylated nucleosomes and appears to have either low levels of acetylation or high levels of phosphorylation in adjacent residues (since adjacent phosphorylation may interfere with the ability to detect acetylation). Inhibition of histone deacetylases (HDACs) by the small molecule TSA reverts this pattern, suggesting that these nucleosomes are specifically deacetylated during mitosis. Altogether, by merging multiple approaches, our study provides evidence to support a model where histone modifications may play a role in mitotic bookmarking and uncovers new insights into the deposition of nucleosomes during mitosis.

Specter: linear deconvolution for targeted analysis of data-independent acquisition mass spectrometry proteomics.

Mass spectrometry with data-independent acquisition (DIA) is a promising method to improve the comprehensiveness and reproducibility of targeted and discovery proteomics, in theory by systematically measuring all peptide precursors in a biological sample. However, the analytical challenges involved in discriminating between peptides with similar sequences in convoluted spectra have limited its applicability in important cases, such as the detection of single-nucleotide polymorphisms (SNPs) and alternative site localizations in phosphoproteomics data. We report Specter (https://github.com/rpeckner-broad/Specter), an open-source software tool that uses linear algebra to deconvolute DIA mixture spectra directly through comparison to a spectral library, thus circumventing the problems associated with typical fragment-correlation-based approaches. We validate the sensitivity of Specter and its performance relative to that of other methods, and show that Specter is able to successfully analyze cases involving highly similar peptides that are typically challenging for DIA analysis methods.

Dual protease type XIII/pepsin digestion offers superior resolution and overlap for the analysis of histone tails by HX-MS.

The N-terminal regions of histone proteins (tails) are dynamic elements that protrude from the nucleosome and are involved in many aspects of chromatin organization. Their epigenetic role is well-established, and post-translational modifications (PTMs) present on these regions contribute to transcriptional regulation. While hydrogen/deuterium exchange mass spectrometry (HX-MS) is well-suited for the analysis of dynamic structures, it has seldom been employed to analyze histones due to the poor N-terminal coverage obtained using pepsin. Here, we test the applicability of a dual protease type XIII/pepsin digestion column to this class of proteins. We optimize online digestion conditions using the H4 monomer, and extend the method to the analysis of histones in monomeric states and nucleosome core particles (NCPs). We show that the dual protease column generates many short and overlapping N-terminal peptides. We evaluate our method by performing hydrogen exchange experiments of NCPs for different time points and present full coverage of the tails at excellent resolution. We further employ electron transfer dissociation and showcase an unprecedented degree of overlap across multiple peptides that is several fold higher than previously reported methods. The method we report here may be readily applied to the HX-MS investigation of histone dynamics and to the footprints of histone binding proteins on nucleosomes.

Chasing Tails: Cathepsin-L Improves Structural Analysis of Histones by HX-MS.

The N-terminal regions (tails) of histone proteins are dynamic elements that protrude from the nucleosome and are involved in many aspects of chromatin organization. Their epigenetic role is well-established, and post-translational modifications present on these regions contribute to transcriptional regulation. Considering their biological significance, relatively few structural details have been established for histone tails, mainly because of their inherently disordered nature. Although hydrogen/deuterium exchange mass spectrometry (HX-MS) is well-suited for the analysis of dynamic structures, it has seldom been employed in this context, presumably because of the poor N-terminal coverage provided by pepsin. Inspired from histone-clipping events, we profiled the activity of cathepsin-L under HX-MS quench conditions and characterized its specificity employing the four core histones (H2A, H2B, H3 and H4). Cathepsin-L demonstrated cleavage patterns that were substrate- and pH-dependent. Cathepsin-L generated overlapping N-terminal peptides about 20 amino acids long for H2A, H3, and H4 proving its suitability for the analysis of histone tails dynamics. We developed a comprehensive HX-MS method in combination with pepsin and obtained full sequence coverage for all histones. We employed our method to analyze histones H3 and H4. We observe rapid deuterium exchange of the N-terminal tails and cooperative unfolding (EX1 kinetics) in the histone-fold domains of histone monomers in-solution. Overall, this novel strategy opens new avenues for investigating the dynamic properties of histones that are not apparent from the crystal structures, providing insights into the structural basis of the histone code.

A Curated Resource for Phosphosite-specific Signature Analysis.

Signaling pathways are orchestrated by post-translational modifications (PTMs) such as phosphorylation. However, pathway analysis of PTM data sets generated by mass spectrometry (MS)-based proteomics is typically performed at a gene-centric level because of the lack of appropriately curated PTM signature databases and bioinformatic tools that leverage PTM site-specific information. Here we present the first version of PTMsigDB, a database of modification site-specific signatures of perturbations, kinase activities and signaling pathways curated from more than 2,500 publications. We adapted the widely used single sample Gene Set Enrichment Analysis approach to utilize PTMsigDB, enabling PTMSignature Enrichment Analysis (PTM-SEA) of quantitative MS data. We used a well-characterized data set of epidermal growth factor (EGF)-perturbed cancer cells to evaluate our approach and demonstrated better representation of signaling events compared with gene-centric methods. We then applied PTM-SEA to analyze the phosphoproteomes of cancer cells treated with cell-cycle inhibitors and detected mechanism-of-action specific signatures of cell cycle kinases. We also applied our methods to analyze the phosphoproteomes of PI3K-inhibited human breast cancer cells and detected signatures of compounds inhibiting PI3K as well as targets downstream of PI3K (AKT, MAPK/ERK) covering a substantial fraction of the PI3K pathway. PTMsigDB and PTM-SEA can be freely accessed at https://github.com/broadinstitute/ssGSEA2.0.

Activity of the PI3K-δ,γ inhibitor duvelisib in a phase 1 trial and preclinical models of T-cell lymphoma.

Duvelisib (IPI-145) is an oral inhibitor of phosphatidylinositol 3-kinase (PI3K)-δ/γ isoforms currently in clinical development. PI3K-δ/γ inhibition may directly inhibit malignant T-cell growth, making duvelisib a promising candidate for patients with peripheral (PTCL) or cutaneous (CTCL) T-cell lymphoma. Inhibition of either isoform may also contribute to clinical responses by modulating nonmalignant immune cells. We investigated these dual effects in a TCL cohort from a phase 1, open-label study of duvelisib in patients with relapsed or refractory PTCL (n = 16) and CTCL (n = 19), along with in vitro and in vivo models of TCL. The overall response rates in patients with PTCL and CTCL were 50.0% and 31.6%, respectively (P = .32). There were 3 complete responses, all among patients with PTCL. Activity was seen across a wide spectrum of subtypes. The most frequently observed grade 3 and 4 adverse events were transaminase increases (40% alanine aminotransferase, 17% aspartate aminotransferase), maculopapular rash (17%), and neutropenia (17%). Responders and nonresponders had markedly different changes in serum cytokine profiles induced by duvelisib. In vitro, duvelisib potently killed 3 of 4 TCL lines with constitutive phospho-AKT (pAKT) vs 0 of 7 lines lacking pAKT (P = .024) and exceeded cell killing by the PI3K-δ-specific inhibitor idelalisib. Administration of duvelisib to mice engrafted with a PTCL patient-derived xenograft resulted in a shift among tumor-associated macrophages from the immunosuppressive M2-like phenotype to the inflammatory M1-like phenotype. In summary, duvelisib demonstrated promising clinical activity and an acceptable safety profile in relapsed/refractory TCL, as well as preclinical evidence of both tumor cell-autonomous and immune-mediated effects. This trial was registered at www.clinicaltrials.gov as #NCT01476657.

Panorama Public: A Public Repository for Quantitative Data Sets Processed in Skyline.

To address the growing need for a centralized, community resource of published results processed with Skyline, and to provide reviewers and readers immediate visual access to the data behind published conclusions, we present Panorama Public (https://panoramaweb.org/public.url), a repository of Skyline documents supporting published results. Panorama Public is built on Panorama, an open source data management system for mass spectrometry data processed with the Skyline targeted mass spectrometry environment. The Panorama web application facilitates viewing, sharing, and disseminating results contained in Skyline documents via a web-browser. Skyline users can easily upload their documents to a Panorama server and allow other researchers to explore uploaded results in the Panorama web-interface through a variety of familiar summary graphs as well as annotated views of the chromatographic peaks processed with Skyline. This makes Panorama ideal for sharing targeted, quantitative results contained in Skyline documents with collaborators, reviewers, and the larger proteomics community. The Panorama Public repository employs the full data visualization capabilities of Panorama which facilitates sharing results with reviewers during manuscript review.

Reduced-representation Phosphosignatures Measured by Quantitative Targeted MS Capture Cellular States and Enable Large-scale Comparison of Drug-induced Phenotypes.

Profiling post-translational modifications represents an alternative dimension to gene expression data in characterizing cellular processes. Many cellular responses to drugs are mediated by changes in cellular phosphosignaling. We sought to develop a common platform on which phosphosignaling responses could be profiled across thousands of samples, and created a targeted MS assay that profiles a reduced-representation set of phosphopeptides that we show to be strong indicators of responses to chemical perturbagens.To develop the assay, we investigated the coordinate regulation of phosphosites in samples derived from three cell lines treated with 26 different bioactive small molecules. Phosphopeptide analytes were selected from these discovery studies by clustering and picking 1 to 2 proxy members from each cluster. A quantitative, targeted parallel reaction monitoring assay was developed to directly measure 96 reduced-representation probes. Sample processing for proteolytic digestion, protein quantification, peptide desalting, and phosphopeptide enrichment have been fully automated, making possible the simultaneous processing of 96 samples in only 3 days, with a plate phosphopeptide enrichment variance of 12%. This highly reproducible process allowed ∼95% of the reduced-representation phosphopeptide probes to be detected in ∼200 samples.The performance of the assay was evaluated by measuring the probes in new samples generated under treatment conditions from discovery experiments, recapitulating the observations of deeper experiments using a fraction of the analytical effort. We measured these probes in new experiments varying the treatments, cell types, and timepoints to demonstrate generalizability. We demonstrated that the assay is sensitive to disruptions in common signaling pathways (e.g. MAPK, PI3K/mTOR, and CDK). The high-throughput, reduced-representation phosphoproteomics assay provides a platform for the comparison of perturbations across a range of biological conditions, suitable for profiling thousands of samples. We believe the assay will prove highly useful for classification of known and novel drug and genetic mechanisms through comparison of phosphoproteomic signatures.

The genome of the green anole lizard and a comparative analysis with birds and mammals.

The evolution of the amniotic egg was one of the great evolutionary innovations in the history of life, freeing vertebrates from an obligatory connection to water and thus permitting the conquest of terrestrial environments. Among amniotes, genome sequences are available for mammals and birds, but not for non-avian reptiles. Here we report the genome sequence of the North American green anole lizard, Anolis carolinensis. We find that A. carolinensis microchromosomes are highly syntenic with chicken microchromosomes, yet do not exhibit the high GC and low repeat content that are characteristic of avian microchromosomes. Also, A. carolinensis mobile elements are very young and diverse-more so than in any other sequenced amniote genome. The GC content of this lizard genome is also unusual in its homogeneity, unlike the regionally variable GC content found in mammals and birds. We describe and assign sequence to the previously unknown A. carolinensis X chromosome. Comparative gene analysis shows that amniote egg proteins have evolved significantly more rapidly than other proteins. An anole phylogeny resolves basal branches to illuminate the history of their repeated adaptive radiations.

Deep, Quantitative Coverage of the Lysine Acetylome Using Novel Anti-acetyl-lysine Antibodies and an Optimized Proteomic Workflow.

Introduction of antibodies specific for acetylated lysine has significantly improved the detection of endogenous acetylation sites by mass spectrometry. Here, we describe a new, commercially available mixture of anti-lysine acetylation (Kac) antibodies and show its utility for in-depth profiling of the acetylome. Specifically, seven complementary monoclones with high specificity for Kac were combined into a final anti-Kac reagent which results in at least a twofold increase in identification of Kac peptides over a commonly used Kac antibody. We outline optimal antibody usage conditions, effective offline basic reversed phase separation, and use of state-of-the-art LC-MS technology for achieving unprecedented coverage of the acetylome. The methods were applied to quantify acetylation sites in suberoylanilide hydroxamic acid-treated Jurkat cells. Over 10,000 Kac peptides from over 3000 Kac proteins were quantified from a single stable isotope labeling by amino acids in cell culture labeled sample using 7.5 mg of peptide input per state. This constitutes the deepest coverage of acetylation sites in quantitative experiments obtained to-date. The approach was also applied to breast tumor xenograft samples using isobaric mass tag labeling of peptides (iTRAQ4, TMT6 and TMT10-plex reagents) for quantification. Greater than 6700 Kac peptides from over 2300 Kac proteins were quantified using 1 mg of tumor protein per iTRAQ 4-plex channel. The novel reagents and methods we describe here enable quantitative, global acetylome analyses with depth and sensitivity approaching that obtained for other well-studied post-translational modifications such as phosphorylation and ubiquitylation, and should have widespread application in biological and clinical studies employing mass spectrometry-based proteomics.

Functional Proteomic Analysis of Repressive Histone Methyltransferase Complexes Reveals ZNF518B as a G9A Regulator.

Cell-type specific gene silencing by histone H3 lysine 27 and lysine 9 methyltransferase complexes PRC2 and G9A-GLP is crucial both during development and to maintain cell identity. Although studying their interaction partners has yielded valuable insight into their functions, how these factors are regulated on a network level remains incompletely understood. Here, we present a new approach that combines quantitative interaction proteomics with global chromatin profiling to functionally characterize repressive chromatin modifying protein complexes in embryonic stem cells. We define binding stoichiometries of 9 new and 12 known interaction partners of PRC2 and 10 known and 29 new interaction partners of G9A-GLP, respectively. We demonstrate that PRC2 and G9A-GLP interact physically and share several interaction partners, including the zinc finger proteins ZNF518A and ZNF518B. Using global chromatin profiling by targeted mass spectrometry, we discover that even sub-stoichiometric binding partners such as ZNF518B can positively regulate global H3K9me2 levels. Biochemical analysis reveals that ZNF518B directly interacts with EZH2 and G9A. Our systematic analysis suggests that ZNF518B may mediate the structural association between PRC2 and G9A-GLP histone methyltransferases and additionally regulates the activity of G9A-GLP.

Atg16L1 T300A variant decreases selective autophagy resulting in altered cytokine signaling and decreased antibacterial defense.

A coding polymorphism (Thr300Ala) in the essential autophagy gene, autophagy related 16-like 1 (ATG16L1), confers increased risk for the development of Crohn disease, although the mechanisms by which single disease-associated polymorphisms contribute to pathogenesis have been difficult to dissect given that environmental factors likely influence disease initiation in these patients. Here we introduce a knock-in mouse model expressing the Atg16L1 T300A variant. Consistent with the human polymorphism, T300A knock-in mice do not develop spontaneous intestinal inflammation, but exhibit morphological defects in Paneth and goblet cells. Selective autophagy is reduced in multiple cell types from T300A knock-in mice compared with WT mice. The T300A polymorphism significantly increases caspase 3- and caspase 7-mediated cleavage of Atg16L1, resulting in lower levels of full-length Atg16Ll T300A protein. Moreover, Atg16L1 T300A is associated with decreased antibacterial autophagy and increased IL-1ß production in primary cells and in vivo. Quantitative proteomics for protein interactors of ATG16L1 identified previously unknown nonoverlapping sets of proteins involved in ATG16L1-dependent antibacterial autophagy or IL-1ß production. These findings demonstrate how the T300A polymorphism leads to cell type- and pathway-specific disruptions of selective autophagy and suggest a mechanism by which this polymorphism contributes to disease.

Automated Microchromatography Enables Multiplexing of Immunoaffinity Enrichment of Peptides to Greater than 150 for Targeted MS-Based Assays.

Immunoaffinity enrichment of peptides coupled with analysis by stable isotope dilution multiple reaction mass spectrometry has been shown to have analytical performance and detection limits suitable for many biomarker verification studies and biological applications. Prior studies have shown that antipeptide antibodies can be multiplexed up to 50 in a single assay without significant loss of performance. Achieving higher multiplex levels is relevant to all studies involving precious biological material as this minimizes the amount of sample that must be consumed to measure a given set of analytes and reduces the assay cost per analyte. Here we developed automated methods employing the Agilent AssayMAP Bravo microchromatography platform and used these methods to characterize the performance of immunoaffinity enrichment of peptides up to multiplex levels of 172. Median capture efficiency for the target peptides remained high (88%) even at levels of 150-plex and declined to 70% at 172-plex compared to antibody performance observed at standard lower multiplex levels (n = 25). Subsequently, we developed and analytically characterized a multiplexed immuno-multiple reaction monitoring-mass spectrometry (immuno-MRM-MS) assay (n = 110) and applied it to measure candidate protein biomarkers of cardiovascular disease in plasma of patients undergoing planned myocardial infarction. The median lower limit of detection of all peptides was 71.5 amol/µL (nM), and the coefficient of variation (CV) was less than 15% at the lower limit of quantification. The results demonstrate that high multiplexed immuno-MRM-MS assays are readily achievable using the optimized sample processing and peptide capture methods described here.

Integrated proteomic analysis of post-translational modifications by serial enrichment.

We report a mass spectrometry-based method for the integrated analysis of protein expression, phosphorylation, ubiquitination and acetylation by serial enrichments of different post-translational modifications (SEPTM) from the same biological sample. This technology enabled quantitative analysis of nearly 8,000 proteins and more than 20,000 phosphorylation, 15,000 ubiquitination and 3,000 acetylation sites per experiment, generating a holistic view of cellular signal transduction pathways as exemplified by analysis of bortezomib-treated human leukemia cells.

Building the Connectivity Map of epigenetics: chromatin profiling by quantitative targeted mass spectrometry.

Epigenetic control of genome function is an important regulatory mechanism in diverse processes such as lineage commitment and environmental sensing, and in disease etiologies ranging from neuropsychiatric disorders to cancer. Here we report a robust, high-throughput targeted, quantitative mass spectrometry (MS) method to rapidly profile modifications of the core histones of chromatin that compose the epigenetic landscape, enabling comparisons among cells with differing genetic backgrounds, genomic perturbations, and drug treatments.