Highlights d We build the genomic and transcriptomic landscape of 465 primary TNBCs d Chinese TNBC cases demonstrate more PIK3CA mutations and LAR subtype d Transcriptomic data classify TNBCs into four subtypes d Multi-omics profiling identifies potential targets within specific TNBC subtypes

The identification of MHC class II restricted peptide epitopes is an important goal in immunological research. A number of computational tools have been developed for this purpose, but there is a lack of large-scale systematic evaluation of their performance. Herein, we used a comprehensive dataset consisting of more than 10,000 previously unpublished MHC-peptide binding affinities, 29 peptide/MHC crystal structures, and 664 peptides experimentally tested for CD4+ T cell responses to systematically evaluate the performances of publicly available MHC class II binding prediction tools. While in selected instances the best tools were associated with AUC values up to 0.86, in general, class II predictions did not perform as well as historically noted for class I predictions. It appears that the ability of MHC class II molecules to bind variable length peptides, which requires the correct assignment of peptide binding cores, is a critical factor limiting the performance of existing prediction tools. To improve performance, we implemented a consensus prediction approach that combines methods with top performances. We show that this consensus approach achieved best overall performance. Finally, we make the large datasets used publicly available as a benchmark to facilitate further development of MHC class II binding peptide prediction methods.

BackgroundMHC class II binding predictions are widely used to identify epitope candidates in infectious agents, allergens, cancer and autoantigens. The vast majority of prediction algorithms for human MHC class II to date have targeted HLA molecules encoded in the DR locus. This reflects a significant gap in knowledge as HLA DP and DQ molecules are presumably equally important, and have only been studied less because they are more difficult to handle experimentally.ResultsIn this study, we aimed to narrow this gap by providing a large scale dataset of over 17,000 HLA-peptide binding affinities for a set of 11 HLA DP and DQ alleles. We also expanded our dataset for HLA DR alleles resulting in a total of 40,000 MHC class II binding affinities covering 26 allelic variants. Utilizing this dataset, we generated prediction tools utilizing several machine learning algorithms and evaluated their performance.ConclusionWe found that 1) prediction methodologies developed for HLA DR molecules perform equally well for DP or DQ molecules. 2) Prediction performances were significantly increased compared to previous reports due to the larger amounts of training data available. 3) The presence of homologous peptides between training and testing datasets should be avoided to give real-world estimates of prediction performance metrics, but the relative ranking of different predictors is largely unaffected by the presence of homologous peptides, and predictors intended for end-user applications should include all training data for maximum performance. 4) The recently developed NN-align prediction method significantly outperformed all other algorithms, including a naïve consensus based on all prediction methods. A new consensus method dropping the comparably weak ARB prediction method could outperform the NN-align method, but further research into how to best combine MHC class II binding predictions is required.

TCRαβ thymocytes differentiate to either CD8αβ cytotoxic T lymphocytes or CD4+ T helper cells. This functional dichotomy is controlled by key transcription factors, including the T helper master regulator, ThPOK, which suppresses the cytolytic program in MHC class II-restricted CD4+ thymocytes. ThPOK continues to repress CD8-lineage genes in mature CD4+ T cells, even as they differentiate to T helper effector subsets. Here we show that the T helper-fate was not fixed and that mature antigen-stimulated CD4+ T cells could terminate Thpok expression and reactivate CD8-lineage genes. This unexpected plasticity resulted in the post-thymic termination of the T helper-program and the functional differentiation of distinct MHC class II-restricted CD4+ cytotoxic T lymphocytes.

Epithelial-to-mesenchymal transition (EMT) is a complex multistep process in which phenotype switches are mediated by a network of transcription factors (TFs). Systematic characterization of all dynamic TFs controlling EMT state transitions, especially for the intermediate partial-EMT state, represents a highly relevant yet largely unexplored task. Here, we performed a computational analysis that integrated time-course EMT transcriptomic data with public cistromic data and identified three synergistic master TFs (ETS2, HNF4A and JUNB) that regulate the transition through the partial-EMT state. Overexpression of these regulators predicted a poor clinical outcome, and their elimination readily abolished TGF-β-induced EMT. Importantly, these factors utilized a clique motif, physically interact and their cumulative binding generally characterized EMT-associated genes. Furthermore, analyses of H3K27ac ChIP-seq data revealed that ETS2, HNF4A and JUNB are associated with super-enhancers and the administration of BRD4 inhibitor readily abolished TGF-β-induced EMT. These findings have implications for systematic discovery of master EMT regulators and super-enhancers as novel targets for controlling metastasis.