Atypical acute myeloid leukemia-specific transcripts generate shared and immunogenic MHC class-I-associated epitopes.

Acute myeloid leukemia (AML) has not benefited from innovative immunotherapies, mainly because of the lack of actionable immune targets. Using an original proteogenomic approach, we analyzed the major histocompatibility complex class I (MHC class I)-associated immunopeptidome of 19 primary AML samples and identified 58 tumor-specific antigens (TSAs). These TSAs bore no mutations and derived mainly (86%) from supposedly non-coding genomic regions. Two AML-specific aberrations were instrumental in the biogenesis of TSAs, intron retention, and epigenetic changes. Indeed, 48% of TSAs resulted from intron retention and translation, and their RNA expression correlated with mutations of epigenetic modifiers (e.g., DNMT3A). AML TSA-coding transcripts were highly shared among patients and were expressed in both blasts and leukemic stem cells. In AML patients, the predicted number of TSAs correlated with spontaneous expansion of cognate T cell receptor clonotypes, accumulation of activated cytotoxic T cells, immunoediting, and improved survival. These TSAs represent attractive targets for AML immunotherapy.

A functional screen to identify novel effectors of hematopoietic stem cell activity.

Despite tremendous progress made toward the identification of the molecular circuitry that governs cell fate in embryonic stem cells, genes controlling this process in the adult hematopoietic stem cell have proven to be more difficult to unmask. We now report the results of a novel gain-of-function screening approach, which identified a series of 18 nuclear factors that affect hematopoietic stem cell activity. Overexpression of ten of these factors resulted in an increased repopulating activity compared to unmanipulated cells. Interestingly, at least four of the 18 factors, Fos, Tcfec, Hmgb1, and Sfpi1, show non-cell-autonomous functions. The utilization of this screening method together with the creation of a database enriched for potential determinants of hematopoietic stem cell self-renewal will serve as a resource to uncover regulatory networks in these cells.

Immunopeptidomic analyses of colorectal cancers with and without microsatellite instability.

Colorectal cancer is the second leading cause of cancer death worldwide, and the incidence of this disease is expected to increase as global socioeconomic changes occur. Immune checkpoint inhibition therapy is effective in treating a minority of colorectal cancer tumors; however, microsatellite stable tumors do not respond well to this treatment. Emerging cancer immunotherapeutic strategies aim to activate a cytotoxic T cell response against tumor-specific antigens, presented exclusively at the cell surface of cancer cells. These antigens are rare and are most effectively identified with a mass spectrometry-based approach, which allows the direct sampling and sequencing of these peptides. Although the few tumor-specific antigens identified to date are derived from coding regions of the genome, recent findings indicate that a large proportion of tumor-specific antigens originate from allegedly noncoding regions. Here, we employed a novel proteogenomic approach to identify tumor antigens in a collection of colorectal cancer-derived cell lines and biopsy samples consisting of matched tumor and normal adjacent tissue. The generation of personalized cancer databases paired with mass spectrometry analyses permitted the identification of more than 30,000 unique MHC I-associated peptides. We identified 19 tumor-specific antigens in both microsatellite stable and unstable tumors, over two-thirds of which were derived from noncoding regions. Many of these peptides were derived from source genes known to be involved in colorectal cancer progression, suggesting that antigens from these genes could have therapeutic potential in a wide range of tumors. These findings could benefit the development of T cell-based vaccines, in which T cells are primed against these antigens to target and eradicate tumors. Such a vaccine could be used in tandem with existing immune checkpoint inhibition therapies, to bridge the gap in treatment efficacy across subtypes of colorectal cancer with varying prognoses. Data are available via ProteomeXchange with identifier PXD028309.

MAPDP: A Cloud-Based Computational Platform for Immunopeptidomics Analyses.

The immunopeptidome corresponds to the repertoire of peptides presented at the cell surface by the major histocompatibility complex (MHC) molecules. Cytotoxic T cells scan this repertoire to identify nonself antigens that can arise from tumors or infected cells. The identification of actionable antigenic targets is key to the development of therapeutic cancer vaccines, T-cell therapy, and other T-cell receptor-based biologics. The growing clinical interest for immunopeptidomics has accelerated the development of high throughput proteogenomic platforms that provide a system-level analysis of MHC-associated peptides. Improvement in sensitivity and throughput of mass spectrometers now allows the detection of a few thousands of peptides from less than 100 million cells. To manage the amount of data generated by these instruments, we have developed the MHC-associated peptide discovery platform (MAPDP), a novel open-source cloud-based computational platform for immunopeptidomic analyses. It provides convenient access from a web portal to immunopeptidomes stored in the database, filtering tools, various visualizations, annotations (e.g., IEDB, dbSNP, gnomAD), peptide-binding affinity prediction (mhcflurry, NetMHC), HLA genotyping, and the generation of personalized proteome databases. MAPDP functionalities are demonstrated here by the discovery of MHC peptides featuring new genetic variants identified in two previously published ovarian carcinoma data sets.

Major multilevel molecular divergence between THP-1 cells from different biorepositories.

The THP-1 cell line is broadly used as a model for acute myeloid leukemia (AML) with MLL fusion and to study monocyte differentiation and function. We studied THP-1 cells obtained from two major biorepositories. The two cell lines were closely related with a percentage match of short tandem repeat (STR) profiles ranging from 93.75% to 100%, depending on the algorithm used. Nevertheless, we found that the two cell lines presented discordant HLA type, cytogenetic aberrations and AML-related gene expression (including critical targets of MLL fusion). These discrepancies resulted mainly from loss of heterozygosity (LOH) involving five chromosomal regions. In view of their aberrant expression of key "leukemia" genes (e.g., LIN28B, MEIS1 and SPARC), we argue that one of the THP-1 cell lines may not be a reliable model for studying leukemia. Their defective expression of HLA molecules and abnormal adhesion properties is also a caveat for studies of antigen presentation. In a more general perspective, our findings show that seemingly minor discrepancies in STR profiles among cell lines may be the sign of major genetic drift, of sufficient magnitude to affect the reliability of cell line-based research.

The expanding roles of Gßγ subunits in G protein-coupled receptor signaling and drug action.

Gßγ subunits from heterotrimeric G proteins perform a vast array of functions in cells with respect to signaling, often independently as well as in concert with Gα subunits. However, the eponymous term "Gßγ" does not do justice to the fact that 5 Gß and 12 Gγ isoforms have evolved in mammals to serve much broader roles beyond their canonical roles in cellular signaling. We explore the phylogenetic diversity of Gßγ subunits with a view toward understanding these expanded roles in different cellular organelles. We suggest that the particular content of distinct Gßγ subunits regulates cellular activity, and that the granularity of individual Gß and Gγ action is only beginning to be understood. Given the therapeutic potential of targeting Gßγ action, this larger view serves as a prelude to more specific development of drugs aimed at individual isoforms.

MHC I-associated peptides preferentially derive from transcripts bearing miRNA response elements.

MHC I-associated peptides (MIPs) play an essential role in normal homeostasis and diverse pathologic conditions. MIPs derive mainly from defective ribosomal products (DRiPs), a subset of nascent proteins that fail to achieve a proper conformation and the physical nature of which remains elusive. In the present study, we used high-throughput proteomic and transcriptomic methods to unravel the structure and biogenesis of MIPs presented by HLA-A and HLA-B molecules on human EBV-infected B lymphocytes from 4 patients. We found that although HLA-different subjects present distinctive MIPs derived from different proteins, these MIPs originate from proteins that are functionally interconnected and implicated in similar biologic pathways. Secondly, the MIP repertoire of human B cells showed no bias toward conserved versus polymorphic genomic sequences, were derived preferentially from abundant transcripts, and conveyed to the cell surface a cell-type-specific signature. Finally, we discovered that MIPs derive preferentially from transcripts bearing miRNA response elements. Furthermore, whereas MIPs of HLA-disparate subjects are coded by different sets of transcripts, these transcripts are regulated by mostly similar miRNAs. Our data support an emerging model in which the generation of MIPs by a transcript depends on its abundance and DRiP rate, which is regulated to a large extent by miRNAs.

Asymmetric segregation and self-renewal of hematopoietic stem and progenitor cells with endocytic Ap2a2.

The stem cell-intrinsic model of self-renewal via asymmetric cell division (ACD) posits that fate determinants be partitioned unequally between daughter cells to either activate or suppress the stemness state. ACD is a purported mechanism by which hematopoietic stem cells (HSCs) self-renew, but definitive evidence for this cellular process remains open to conjecture. To address this issue, we chose 73 candidate genes that function within the cell polarity network to identify potential determinants that may concomitantly alter HSC fate while also exhibiting asymmetric segregation at cell division. Initial gene-expression profiles of polarity candidates showed high and differential expression in both HSCs and leukemia stem cells. Altered HSC fate was assessed by our established in vitro to in vivo screen on a subcohort of candidate polarity genes, which revealed 6 novel positive regulators of HSC function: Ap2a2, Gpsm2, Tmod1, Kif3a, Racgap1, and Ccnb1. Interestingly, live-cell videomicroscopy of the endocytic protein AP2A2 shows instances of asymmetric segregation during HSC/progenitor cell cytokinesis. These results contribute further evidence that ACD is functional in HSC self-renewal, suggest a role for Ap2a2 in HSC activity, and provide a unique opportunity to prospectively analyze progeny from HSC asymmetric divisions.

Breast cancer immunopeptidomes contain numerous shared tumor antigens.

Hormone receptor-positive breast cancer (HR+) is immunologically cold and has not benefited from advances in immunotherapy. In contrast, subsets of triple-negative breast cancer (TNBC) display high leukocytic infiltration and respond to checkpoint blockade. CD8+ T cells, the main effectors of anticancer responses, recognize MHC I-associated peptides (MAPs). Our work aimed to characterize the repertoire of MAPs presented by HR+ and TNBC tumors. Using mass spectrometry, we identified 57,094 unique MAPs in 26 primary breast cancer samples. MAP source genes highly overlapped between both subtypes. We identified 25 tumor-specific antigens (TSAs) mainly deriving from aberrantly expressed regions. TSAs were most frequently identified in TNBC samples and were more shared among The Cancer Genome Atlas (TCGA) database TNBC than HR+ samples. In the TNBC cohort, the predicted number of TSAs positively correlated with leukocytic infiltration and overall survival, supporting their immunogenicity in vivo. We detected 49 tumor-associated antigens (TAAs), some of which derived from cancer-associated fibroblasts. Functional expansion of specific T cell assays confirmed the in vitro immunogenicity of several TSAs and TAAs. Our study identified attractive targets for cancer immunotherapy in both breast cancer subtypes. The higher prevalence of TSAs in TNBC tumors provides a rationale for their responsiveness to checkpoint blockade.

Autophagy degrades immunogenic endogenous retroelements induced by 5-azacytidine in acute myeloid leukemia.

The hypomethylating agent 5-azacytidine (AZA) is the first-line treatment for AML patients unfit for intensive chemotherapy. The effect of AZA results in part from T-cell cytotoxic responses against MHC-I-associated peptides (MAPs) deriving from hypermethylated genomic regions such as cancer-testis antigens (CTAs), or endogenous retroelements (EREs). However, evidence supporting higher ERE MAPs presentation after AZA treatment is lacking. Therefore, using proteogenomics, we examined the impact of AZA on the repertoire of MAPs and their source transcripts. AZA-treated AML upregulated both CTA and ERE transcripts, but only CTA MAPs were presented at greater levels. Upregulated ERE transcripts triggered innate immune responses against double-stranded RNAs but were degraded by autophagy, and not processed into MAPs. Autophagy resulted from the formation of protein aggregates caused by AZA-dependent inhibition of DNMT2. Autophagy inhibition had an additive effect with AZA on AML cell proliferation and survival, increased ERE levels, increased pro-inflammatory responses, and generated immunogenic tumor-specific ERE-derived MAPs. Finally, autophagy was associated with a lower abundance of CD8+ T-cell markers in AML patients expressing high levels of EREs. This work demonstrates that AZA-induced EREs are degraded by autophagy and shows that inhibiting autophagy can improve the immune recognition of AML blasts in treated patients.

Transposable elements regulate thymus development and function.

Transposable elements (TEs) are repetitive sequences representing ~45% of the human and mouse genomes and are highly expressed by medullary thymic epithelial cells (mTECs). In this study, we investigated the role of TEs on T-cell development in the thymus. We performed multiomic analyses of TEs in human and mouse thymic cells to elucidate their role in T-cell development. We report that TE expression in the human thymus is high and shows extensive age- and cell lineage-related variations. TE expression correlates with multiple transcription factors in all cell types of the human thymus. Two cell types express particularly broad TE repertoires: mTECs and plasmacytoid dendritic cells (pDCs). In mTECs, transcriptomic data suggest that TEs interact with transcription factors essential for mTEC development and function (e.g., PAX1 and REL), and immunopeptidomic data showed that TEs generate MHC-I-associated peptides implicated in thymocyte education. Notably, AIRE, FEZF2, and CHD4 regulate small yet non-redundant sets of TEs in murine mTECs. Human thymic pDCs homogenously express large numbers of TEs that likely form dsRNA, which can activate innate immune receptors, potentially explaining why thymic pDCs constitutively secrete IFN ɑ/ß. This study highlights the diversity of interactions between TEs and the adaptive immune system. TEs are genetic parasites, and the two thymic cell types most affected by TEs (mTEcs and pDCs) are essential to establishing central T-cell tolerance. Therefore, we propose that orchestrating TE expression in thymic cells is critical to prevent autoimmunity in vertebrates.

Genome-wide interrogation of Mammalian stem cell fate determinants by nested chromosome deletions.

Understanding the function of important DNA elements in mammalian stem cell genomes would be enhanced by the availability of deletion collections in which segmental haploidies are precisely characterized. Using a modified Cre-loxP-based system, we now report the creation and characterization of a collection of ∼1,300 independent embryonic stem cell (ESC) clones enriched for nested chromosomal deletions. Mapping experiments indicate that this collection spans over 25% of the mouse genome with good representative coverage of protein-coding genes, regulatory RNAs, and other non-coding sequences. This collection of clones was screened for in vitro defects in differentiation of ESC into embryoid bodies (EB). Several putative novel haploinsufficient regions, critical for EB development, were identified. Functional characterization of one of these regions, through BAC complementation, identified the ribosomal gene Rps14 as a novel haploinsufficient determinant of embryoid body formation. This new library of chromosomal deletions in ESC (DelES: http://bioinfo.iric.ca/deles) will serve as a unique resource for elucidation of novel protein-coding and non-coding regulators of ESC activity.

BamQuery: a proteogenomic tool to explore the immunopeptidome and prioritize actionable tumor antigens.

MHC-I-associated peptides deriving from non-coding genomic regions and mutations can generate tumor-specific antigens, including neoantigens. Quantifying tumor-specific antigens' RNA expression in malignant and benign tissues is critical for discriminating actionable targets. We present BamQuery, a tool attributing an exhaustive RNA expression to MHC-I-associated peptides of any origin from bulk and single-cell RNA-sequencing data. We show that many cryptic and mutated tumor-specific antigens can derive from multiple discrete genomic regions, abundantly expressed in normal tissues. BamQuery can also be used to predict MHC-I-associated peptides immunogenicity and identify actionable tumor-specific antigens de novo.

The MHC I immunopeptidome conveys to the cell surface an integrative view of cellular regulation.

Self/non-self discrimination is a fundamental requirement of life. Endogenous peptides presented by major histocompatibility complex class I (MHC I) molecules represent the essence of self for CD8 T lymphocytes. These MHC I peptides (MIPs) are collectively referred to as the immunopeptidome. From a systems-level perspective, very little is known about the origin, composition and plasticity of the immunopeptidome. Here, we show that the immunopeptidome, and therefore the nature of the immune self, is plastic and moulded by cellular metabolic activity. By using a quantitative high-throughput mass spectrometry-based approach, we found that altering cellular metabolism via the inhibition of the mammalian target of rapamycin results in dynamic changes in the cell surface MIPs landscape. Moreover, we provide systems-level evidence that the immunopeptidome projects at the cell surface a representation of biochemical networks and metabolic events regulated at multiple levels inside the cell. Our findings open up new perspectives in systems immunology and predictive biology. Indeed, predicting variations in the immunopeptidome in response to cell-intrinsic and -extrinsic factors could be relevant to the rational design of immunotherapeutic interventions.

Two types of human TCR differentially regulate reactivity to self and non-self antigens.

Based on analyses of TCR sequences from over 1,000 individuals, we report that the TCR repertoire is composed of two ontogenically and functionally distinct types of TCRs. Their production is regulated by variations in thymic output and terminal deoxynucleotidyl transferase (TDT) activity. Neonatal TCRs derived from TDT-negative progenitors persist throughout life, are highly shared among subjects, and are reported as disease-associated. Thus, 10%-30% of most frequent cord blood TCRs are associated with common pathogens and autoantigens. TDT-dependent TCRs present distinct structural features and are less shared among subjects. TDT-dependent TCRs are produced in maximal numbers during infancy when thymic output and TDT activity reach a summit, are more abundant in subjects with AIRE mutations, and seem to play a dominant role in graft-versus-host disease. Factors decreasing thymic output (age, male sex) negatively impact TCR diversity. Males compensate for their lower repertoire diversity via hyperexpansion of selected TCR clonotypes.

Induced pluripotent stem cells display a distinct set of MHC I-associated peptides shared by human cancers.

Previous reports showed that mouse vaccination with pluripotent stem cells (PSCs) induces durable anti-tumor immune responses via T cell recognition of some elusive oncofetal epitopes. We characterize the MHC I-associated peptide (MAP) repertoire of human induced PSCs (iPSCs) using proteogenomics. Our analyses reveal a set of 46 pluripotency-associated MAPs (paMAPs) absent from the transcriptome of normal tissues and adult stem cells but expressed in PSCs and multiple adult cancers. These paMAPs derive from coding and allegedly non-coding (48%) transcripts involved in pluripotency maintenance, and their expression in The Cancer Genome Atlas samples correlates with source gene hypomethylation and genomic aberrations common across cancer types. We find that several of these paMAPs were immunogenic. However, paMAP expression in tumors coincides with activation of pathways instrumental in immune evasion (WNT, TGF-ß, and CDK4/6). We propose that currently available inhibitors of these pathways could synergize with immune targeting of paMAPs for the treatment of poorly differentiated cancers.

PSMB11 regulates gene expression in cortical thymic epithelial cells.

The PSMB11 proteasomal subunit, expressed only in cortical thymic epithelial cells (cTECs), is essential for the development of functional CD8+ T cells. An attractive yet unproven theory holds that PSMB11 generates unique major histocompatibility complex class I (MHC I)-associated peptides required for positive selection. We recently reported that PSMB11 regulates the expression of hundreds of genes in cTECs, mainly by differential proteolysis of transcription factors. Thereby, PSMB11 maintains the distinctness of cTECs relative to medullary TECs (mTECs) and promotes cortex-to-medulla migration of developing thymocytes. These conclusions have been challenged by Ohigashi and colleagues, who suggest that their data show that PSMB11 uniquely controls antigen presentation without affecting cTEC biology. Here, we perform a comprehensive reanalysis of transcriptomic and proteomic data from the Ohigashi lab and confirm our original conclusions. This Matters Arising paper is in response to Ohigashi et al. (2019), published in Cell Reports. See also the response by Ohigashi and Takahama (2021), published in this issue of Cell Reports.

Widespread and tissue-specific expression of endogenous retroelements in human somatic tissues.

BACKGROUND: Endogenous retroelements (EREs) constitute about 42% of the human genome and have been implicated in common human diseases such as autoimmunity and cancer. The dominant paradigm holds that EREs are expressed in embryonic stem cells (ESCs) and germline cells but are repressed in differentiated somatic cells. Despite evidence that some EREs can be expressed at the RNA and protein levels in specific contexts, a system-level evaluation of their expression in human tissues is lacking. METHODS: Using RNA sequencing data, we analyzed ERE expression in 32 human tissues and cell types, including medullary thymic epithelial cells (mTECs). A tissue specificity index was computed to identify tissue-restricted ERE families. We also analyzed the transcriptome of mTECs in wild-type and autoimmune regulator (AIRE)-deficient mice. Finally, we developed a proteogenomic workflow combining RNA sequencing and mass spectrometry (MS) in order to evaluate whether EREs might be translated and generate MHC I-associated peptides (MAP) in B-lymphoblastoid cell lines (B-LCL) from 16 individuals. RESULTS: We report that all human tissues express EREs, but the breadth and magnitude of ERE expression are very heterogeneous from one tissue to another. ERE expression was particularly high in two MHC I-deficient tissues (ESCs and testis) and one MHC I-expressing tissue, mTECs. In mutant mice, we report that the exceptional expression of EREs in mTECs was AIRE-independent. MS analyses identified 103 non-redundant ERE-derived MAPs (ereMAPs) in B-LCLs. These ereMAPs preferentially derived from sense translation of intronic EREs. Notably, detailed analyses of their amino acid composition revealed that ERE-derived MAPs presented homology to viral MAPs. CONCLUSIONS: This study shows that ERE expression in somatic tissues is more pervasive and heterogeneous than anticipated. The high and diversified expression of EREs in mTECs and their ability to generate MAPs suggest that EREs may play an important role in the establishment of self-tolerance. The viral-like properties of ERE-derived MAPs suggest that those not expressed in mTECs can be highly immunogenic.

Proteogenomics Uncovers a Vast Repertoire of Shared Tumor-Specific Antigens in Ovarian Cancer.

High-grade serous ovarian cancer (HGSC), the principal cause of death from gynecologic malignancies in the world, has not significantly benefited from advances in cancer immunotherapy. Although HGSC infiltration by lymphocytes correlates with superior survival, the nature of antigens that can elicit anti-HGSC immune responses is unknown. The goal of this study was to establish the global landscape of HGSC tumor-specific antigens (TSA) using a mass spectrometry pipeline that interrogated all reading frames of all genomic regions. In 23 HGSC tumors, we identified 103 TSAs. Classic TSA discovery approaches focusing only on mutated exonic sequences would have uncovered only three of these TSAs. Other mutated TSAs resulted from out-of-frame exonic translation (n = 2) or from noncoding sequences (n = 7). One group of TSAs (n = 91) derived from aberrantly expressed unmutated genomic sequences, which were not expressed in normal tissues. These aberrantly expressed TSAs (aeTSA) originated primarily from nonexonic sequences, in particular intronic (29%) and intergenic (22%) sequences. Their expression was regulated at the transcriptional level by variations in gene copy number and DNA methylation. Although mutated TSAs were unique to individual tumors, aeTSAs were shared by a large proportion of HGSCs. Taking into account the frequency of aeTSA expression and HLA allele frequencies, we calculated that, in Caucasians, the median number of aeTSAs per tumor would be five. We conclude that, in view of their number and the fact that they are shared by many tumors, aeTSAs may be the most attractive targets for HGSC immunotherapy.

Noncoding regions are the main source of targetable tumor-specific antigens.

Tumor-specific antigens (TSAs) represent ideal targets for cancer immunotherapy, but few have been identified thus far. We therefore developed a proteogenomic approach to enable the high-throughput discovery of TSAs coded by potentially all genomic regions. In two murine cancer cell lines and seven human primary tumors, we identified a total of 40 TSAs, about 90% of which derived from allegedly noncoding regions and would have been missed by standard exome-based approaches. Moreover, most of these TSAs derived from nonmutated yet aberrantly expressed transcripts (such as endogenous retroelements) that could be shared by multiple tumor types. Last, we demonstrated that, in mice, the strength of antitumor responses after TSA vaccination was influenced by two parameters that can be estimated in humans and could serve for TSA prioritization in clinical studies: TSA expression and the frequency of TSA-responsive T cells in the preimmune repertoire. In conclusion, the strategy reported herein could considerably facilitate the identification and prioritization of actionable human TSAs.