Impact of predicted HLA class I immunopeptidome on viral reservoir in a cohort of people living with HIV in Italy.

The class I HLA genotype has been widely recognized as a factor influencing HIV disease progression in treatment-naïve subjects. However, little is known regarding its role in HIV disease course and how it influences the size of the viral reservoir once anti-retroviral therapy (ART) is started. Here, leveraging on cutting-edge bioinformatic tools, we explored the relationship between HLA class I and the HIV reservoir in a cohort of 90 people living with HIV (PLWH) undergoing ART and who achieved viral suppression. Analysis of HLA allele distribution among patients with high and low HIV reservoir allowed us to document a predominant role of HLA-B and -C genes in regulating the size of HIV reservoir. We then focused on the analysis of HIV antigen (Ag) repertoire, by investigating immunogenetic parameters such as the degree of homozygosity, HLA evolutionary distance and Ag load. In particular, we used two different bioinformatic algorithms, NetMHCpan and MixMHCpred, to predict HLA presentation of immunogenic HIV-derived peptides and identified HLA-B*57:01 and HLA-B*58:01 among the highest ranking HLAs in terms of total load, suggesting that their previously reported protective role against HIV disease progression might be linked to a more effective viral recognition and presentation to Cytotoxic T lymphocytes (CTLs). Further, we speculated that some peptide-HLA complexes, including those produced by the interaction between HLA-B*27 and the HIV Gag protein, might be particularly relevant for the efficient regulation of HIV replication and containment of the HIV reservoir. Last, we provide evidence of a possible synergistic effect between the CCR5 ∆32 mutation and Ag load in controlling HIV reservoir.

Harnessing T cell exhaustion and trogocytosis to isolate patient-derived tumor-specific TCR.

To study and then harness the tumor-specific T cell dynamics after allogeneic hematopoietic stem cell transplant, we typed the frequency, phenotype, and function of lymphocytes directed against tumor-associated antigens (TAAs) in 39 consecutive transplanted patients, for 1 year after transplant. We showed that TAA-specific T cells circulated in 90% of patients but display a limited effector function associated to an exhaustion phenotype, particularly in the subgroup of patients deemed to relapse, where exhausted stem cell memory T cells accumulated. Accordingly, cancer-specific cytolytic functions were relevant only when the TAA-specific T cell receptors (TCRs) were transferred into healthy, genome-edited T cells. We then exploited trogocytosis and ligandome-on-chip technology to unveil the specificities of tumor-specific TCRs retrieved from the exhausted T cell pool. Overall, we showed that harnessing circulating TAA-specific and exhausted T cells allow to isolate TCRs against TAAs and previously not described acute myeloid leukemia antigens, potentially relevant for T cell-based cancer immunotherapy.

Driver mutations in GNAQ and GNA11 genes as potential targets for precision immunotherapy in uveal melanoma patients.

Uveal melanoma (UM) is the most common ocular malignancy in adults. Nearly 95% of UM patients carry the mutually exclusive mutations in the homologous genes GNAQ (amino acid change Q209L/Q209P) and GNA11 (aminoacid change Q209L). UM is located in an immunosuppressed organ and does not suffer immunoediting. Therefore, we hypothesize that driver mutations in GNAQ/11 genes could be recognized by the immune system. Genomic and transcriptomic data from primary uveal tumors were collected from the TCGA-UM dataset (n = 80) and used to assess the immunogenic potential for GNAQ/GNA11 Q209L/Q209P mutations using a variety of tools and HLA type information. All prediction tools showed stronger GNAQ/11 Q209L binding to HLA than GNAQ/11 Q209P. The immunogenicity analysis revealed that Q209L is likely to be presented by more than 73% of individuals in 1000 G databases whereas Q209P is only predicted to be presented in 24% of individuals. GNAQ/11 Q209L showed a higher likelihood to be presented by HLA-I molecules than almost all driver mutations analyzed. Finally, samples carrying Q209L had a higher immune-reactive phenotype. Regarding cancer risk, seven HLA genotypes with low Q209L affinity show higher frequency in uveal melanoma patients than in the general population. However, no clear association was found between any HLA genotype and survival. Results suggest a high potential immunogenicity of the GNAQ/11 Q209L variant that could allow the generation of novel therapeutic tools to treat UM like neoantigen vaccinations.

Phase II trial of hypomethylating agent combined with nivolumab for acute myeloid leukaemia relapse after allogeneic haematopoietic cell transplantation-Immune signature correlates with response.

Acute myeloid leukaemia (AML) relapse after allogeneic haematopoietic cell transplantation (allo-HCT) is often driven by immune-related mechanisms and associated with poor prognosis. Immune checkpoint inhibitors combined with hypomethylating agents (HMA) may restore or enhance the graft-versus-leukaemia effect. Still, data about using this combination regimen after allo-HCT are limited. We conducted a prospective, phase II, open-label, single-arm study in which we treated patients with haematological AML relapse after allo-HCT with HMA plus the anti-PD-1 antibody nivolumab. The response was correlated with DNA-, RNA- and protein-based single-cell technology assessments to identify biomarkers associated with therapeutic efficacy. Sixteen patients received a median number of 2 (range 1-7) nivolumab applications. The overall response rate (CR/PR) at day 42 was 25%, and another 25% of the patients achieved stable disease. The median overall survival was 15.6 months. High-parametric cytometry documented a higher frequency of activated (ICOS+ , HLA-DR+ ), low senescence (KLRG1- , CD57- ) CD8+ effector T cells in responders. We confirmed these findings in a preclinical model. Single-cell transcriptomics revealed a pro-inflammatory rewiring of the expression profile of T and myeloid cells in responders. In summary, the study indicates that the post-allo-HCT HMA/nivolumab combination induces anti-AML immune responses in selected patients and could be considered as a bridging approach to a second allo-HCT. Trial-registration: EudraCT-No. 2017-002194-18.

Preexisting Immunity Drives the Response to Neoadjuvant Chemotherapy in Esophageal Adenocarcinoma.

Current treatment for patients with locally advanced esophageal adenocarcinoma (EAC) is neoadjuvant chemotherapy (nCT), alone or combined with radiotherapy, before surgery. However, fewer than 30% of treated patients show a pathologic complete response to nCT, which correlates with increased 5-year survival compared with nonresponders. Understanding the mechanisms of response to nCT is pivotal to better stratify patients and inform more efficacious therapies. Here, we investigated the immune mechanisms involved in nCT response by multidimensional profiling of pretreatment tumor biopsies and blood from 68 patients with EAC (34 prospectively and 34 retrospectively collected), comparing complete responders versus nonresponders to nCT. At the tumor level, complete response to nCT was associated with molecular signatures of immune response and proliferation, increased putative antitumor tissue-resident memory CD39+ CD103+ CD8+ T cells, and reduced immunosuppressive T regulatory cells (Treg) and M2-like macrophages. Systemically, complete responders showed higher frequencies of immunostimulatory CD14+ CD11c+ HLA-DRhigh cells, and reduced programmed cell death ligand 1-positive (PD-L1+) monocytic myeloid-derived suppressor cells, along with high plasma GM-CSF (proinflammatory) and low IL4, CXCL10, C3a, and C5a (suppressive). Plasma proinflammatory and suppressive cytokines correlated directly and inversely, respectively, with the frequency of tumor-infiltrating CD39+ CD103+ CD8+ T cells. These results suggest that preexisting immunity in baseline tumor drives the clinical activity of nCT in locally advanced EAC. Furthermore, it may be possible to stratify patients based on predictive immune signatures, enabling tailored neoadjuvant and/or adjuvant regimens. SIGNIFICANCE: Multidimensional profiling of pretreatment esophageal adenocarcinoma shows patient response to nCT is correlated with active preexisting immunity and indicates molecular pathways of resistance that may be targeted to improve clinical outcomes.

T-cell receptor determinants of response to chemoradiation in locally-advanced HPV16-driven malignancies.

Background: The effect of chemoradiation on the anti-cancer immune response is being increasingly acknowledged; however, its clinical implications in treatment responses are yet to be fully understood. Human papillomavirus (HPV)-driven malignancies express viral oncogenic proteins which may serve as tumor-specific antigens and represent ideal candidates for monitoring the peripheral T-cell receptor (TCR) changes secondary to chemoradiotherapy (CRT). Methods: We performed intra-tumoral and pre- and post-treatment peripheral TCR sequencing in a cohort of patients with locally-advanced HPV16-positive cancers treated with CRT. An in silico computational pipeline was used to cluster TCR repertoire based on epitope-specificity and to predict affinity between these clusters and HPV16-derived epitopes. Results: Intra-tumoral repertoire diversity, intra-tumoral and post-treatment peripheral CDR3ß similarity clustering were predictive of response. In responders, CRT triggered an increase peripheral TCR clonality and clonal relatedness. Post-treatment expansion of baseline peripheral dominant TCRs was associated with response. Responders showed more baseline clustered structures of TCRs maintained post-treatment and displayed significantly more maintained clustered structures. When applying clustering by TCR-specificity methods, responders displayed a higher proportion of intra-tumoral TCRs predicted to recognise HPV16 peptides. Conclusions: Baseline TCR characteristics and changes in the peripheral T-cell clones triggered by CRT are associated with treatment outcome. Maintenance and boosting of pre-existing clonotypes are key elements of an effective anti-cancer immune response driven by CRT, supporting a paradigm in which the immune system plays a central role in the success of CRT in current standard-of-care protocols.

Targeting MDM2 enhances antileukemia immunity after allogeneic transplantation via MHC-II and TRAIL-R1/2 upregulation.

Patients with acute myeloid leukemia (AML) often achieve remission after allogeneic hematopoietic cell transplantation (allo-HCT) but subsequently die of relapse driven by leukemia cells resistant to elimination by allogeneic T cells based on decreased major histocompatibility complex II (MHC-II) expression and apoptosis resistance. Here we demonstrate that mouse-double-minute-2 (MDM2) inhibition can counteract immune evasion of AML. MDM2 inhibition induced MHC class I and II expression in murine and human AML cells. Using xenografts of human AML and syngeneic mouse models of leukemia, we show that MDM2 inhibition enhanced cytotoxicity against leukemia cells and improved survival. MDM2 inhibition also led to increases in tumor necrosis factor-related apoptosis-inducing ligand receptor-1 and -2 (TRAIL-R1/2) on leukemia cells and higher frequencies of CD8+CD27lowPD-1lowTIM-3low T cells, with features of cytotoxicity (perforin+CD107a+TRAIL+) and longevity (bcl-2+IL-7R+). CD8+ T cells isolated from leukemia-bearing MDM2 inhibitor-treated allo-HCT recipients exhibited higher glycolytic activity and enrichment for nucleotides and their precursors compared with vehicle control subjects. T cells isolated from MDM2 inhibitor-treated AML-bearing mice eradicated leukemia in secondary AML-bearing recipients. Mechanistically, the MDM2 inhibitor-mediated effects were p53-dependent because p53 knockdown abolished TRAIL-R1/2 and MHC-II upregulation, whereas p53 binding to TRAILR1/2 promotors increased upon MDM2 inhibition. The observations in the mouse models were complemented by data from human individuals. Patient-derived AML cells exhibited increased TRAIL-R1/2 and MHC-II expression on MDM2 inhibition. In summary, we identified a targetable vulnerability of AML cells to allogeneic T-cell-mediated cytotoxicity through the restoration of p53-dependent TRAIL-R1/2 and MHC-II production via MDM2 inhibition.

The Immunogenic Potential of Recurrent Cancer Drug Resistance Mutations: An In Silico Study.

Cancer somatic mutations have been identified as a source of antigens that can be targeted by cancer immunotherapy. In this work, expanding on previous studies, we analyze the HLA-presentation properties of mutations that are known to drive resistance to cancer targeted-therapies. We survey a large dataset of mutations that confer resistance to different drugs and occur in numerous genes and tumor types. We show that a significant number of them are predicted in silico to be potentially immunogenic across a large proportion of the human population. Further, by analyzing a cohort of patients carrying a small subset of these resistance mutations, we provide evidence that what is observed in the general population may be indicative of the mutations' immunogenic potential in resistant patients. Two of the mutations in our dataset had previously been experimentally validated by others and it was confirmed that some of their associated neopeptides elicit T-cell responses in vitro. The identification of potent cancer-specific antigens can be instrumental for developing more effective immunotherapies. In this work, we propose a novel list of drug-resistance mutations, several of which are recurrent, that could be of particular interest in the context of off-the-shelf precision immunotherapies such as therapeutic cancer vaccines.

Clinical BRCA1/2 Reversion Analysis Identifies Hotspot Mutations and Predicted Neoantigens Associated with Therapy Resistance.

Reversion mutations in BRCA1 or BRCA2 are associated with resistance to PARP inhibitors and platinum. To better understand the nature of these mutations, we collated, codified, and analyzed more than 300 reversions. This identified reversion "hotspots" and "deserts" in regions encoding the N and C terminus, respectively, of BRCA2, suggesting that pathogenic mutations in these regions may be at higher or lower risk of reversion. Missense and splice-site pathogenic mutations in BRCA1/2 also appeared less likely to revert than truncating mutations. Most reversions were <100 bp deletions. Although many deletions exhibited microhomology, this was not universal, suggesting that multiple DNA-repair processes cause reversion. Finally, we found that many reversions were predicted to encode immunogenic neopeptides, suggesting a route to the treatment of reverted disease. As well as providing a freely available database for the collation of future reversion cases, these observations have implications for how drug resistance might be managed in BRCA-mutant cancers. SIGNIFICANCE: Reversion mutations in BRCA genes are a major cause of clinical platinum and PARP inhibitor resistance. This analysis of all reported clinical reversions suggests that the position of BRCA2 mutations affects the risk of reversion. Many reversions are also predicted to encode tumor neoantigens, providing a potential route to targeting resistance.This article is highlighted in the In This Issue feature, p. 1426.

Author Correction: Extreme intratumour heterogeneity and driver evolution in mismatch repair deficient gastro-oesophageal cancer.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

Extreme intratumour heterogeneity and driver evolution in mismatch repair deficient gastro-oesophageal cancer.

Mismatch repair deficient (dMMR) gastro-oesophageal adenocarcinomas (GOAs) show better outcomes than their MMR-proficient counterparts and high immunotherapy sensitivity. The hypermutator-phenotype of dMMR tumours theoretically enables high evolvability but their evolution has not been investigated. Here we apply multi-region exome sequencing (MSeq) to four treatment-naive dMMR GOAs. This reveals extreme intratumour heterogeneity (ITH), exceeding ITH in other cancer types >20-fold, but also long phylogenetic trunks which may explain the exquisite immunotherapy sensitivity of dMMR tumours. Subclonal driver mutations are common and parallel evolution occurs in RAS, PIK3CA, SWI/SNF-complex genes and in immune evasion regulators. MSeq data and evolution analysis of single region-data from 64 MSI GOAs show that chromosome 8 gains are early genetic events and that the hypermutator-phenotype remains active during progression. MSeq may be necessary for biomarker development in these heterogeneous cancers. Comparison with other MSeq-analysed tumour types reveals mutation rates and their timing to determine phylogenetic tree morphologies.

Identification of single nucleotide variants using position-specific error estimation in deep sequencing data.

BACKGROUND: Targeted deep sequencing is a highly effective technology to identify known and novel single nucleotide variants (SNVs) with many applications in translational medicine, disease monitoring and cancer profiling. However, identification of SNVs using deep sequencing data is a challenging computational problem as different sequencing artifacts limit the analytical sensitivity of SNV detection, especially at low variant allele frequencies (VAFs). METHODS: To address the problem of relatively high noise levels in amplicon-based deep sequencing data (e.g. with the Ion AmpliSeq technology) in the context of SNV calling, we have developed a new bioinformatics tool called AmpliSolve. AmpliSolve uses a set of normal samples to model position-specific, strand-specific and nucleotide-specific background artifacts (noise), and deploys a Poisson model-based statistical framework for SNV detection. RESULTS: Our tests on both synthetic and real data indicate that AmpliSolve achieves a good trade-off between precision and sensitivity, even at VAF below 5% and as low as 1%. We further validate AmpliSolve by applying it to the detection of SNVs in 96 circulating tumor DNA samples at three clinically relevant genomic positions and compare the results to digital droplet PCR experiments. CONCLUSIONS: AmpliSolve is a new tool for in-silico estimation of background noise and for detection of low frequency SNVs in targeted deep sequencing data. Although AmpliSolve has been specifically designed for and tested on amplicon-based libraries sequenced with the Ion Torrent platform it can, in principle, be applied to other sequencing platforms as well. AmpliSolve is freely available at https://github.com/dkleftogi/AmpliSolve .

Genomic and Transcriptomic Determinants of Therapy Resistance and Immune Landscape Evolution during Anti-EGFR Treatment in Colorectal Cancer.

Despite biomarker stratification, the anti-EGFR antibody cetuximab is only effective against a subgroup of colorectal cancers (CRCs). This genomic and transcriptomic analysis of the cetuximab resistance landscape in 35 RAS wild-type CRCs identified associations of NF1 and non-canonical RAS/RAF aberrations with primary resistance and validated transcriptomic CRC subtypes as non-genetic predictors of benefit. Sixty-four percent of biopsies with acquired resistance harbored no genetic resistance drivers. Most of these had switched from a cetuximab-sensitive transcriptomic subtype at baseline to a fibroblast- and growth factor-rich subtype at progression. Fibroblast-supernatant conferred cetuximab resistance in vitro, confirming a major role for non-genetic resistance through stromal remodeling. Cetuximab treatment increased cytotoxic immune infiltrates and PD-L1 and LAG3 immune checkpoint expression, potentially providing opportunities to treat cetuximab-resistant CRCs with immunotherapy.

Microenvironmental niche divergence shapes BRCA1-dysregulated ovarian cancer morphological plasticity.

How tumor microenvironmental forces shape plasticity of cancer cell morphology is poorly understood. Here, we conduct automated histology image and spatial statistical analyses in 514 high grade serous ovarian samples to define cancer morphological diversification within the spatial context of the microenvironment. Tumor spatial zones, where cancer cell nuclei diversify in shape, are mapped in each tumor. Integration of this spatially explicit analysis with omics and clinical data reveals a relationship between morphological diversification and the dysregulation of DNA repair, loss of nuclear integrity, and increased disease mortality. Within the Immunoreactive subtype, spatial analysis further reveals significantly lower lymphocytic infiltration within diversified zones compared with other tumor zones, suggesting that even immune-hot tumors contain cells capable of immune escape. Our findings support a model whereby a subpopulation of morphologically plastic cancer cells with dysregulated DNA repair promotes ovarian cancer progression through positive selection by immune evasion.

MetalDetector: a web server for predicting metal-binding sites and disulfide bridges in proteins from sequence.

UNLABELLED: The web server MetalDetector classifies histidine residues in proteins into one of two states (free or metal bound) and cysteines into one of three states (free, metal bound or disulfide bridged). A decision tree integrates predictions from two previously developed methods (DISULFIND and Metal Ligand Predictor). Cross-validated performance assessment indicates that our server predicts disulfide bonding state at 88.6% precision and 85.1% recall, while it identifies cysteines and histidines in transition metal-binding sites at 79.9% precision and 76.8% recall, and at 60.8% precision and 40.7% recall, respectively. AVAILABILITY: Freely available at http://metaldetector.dsi.unifi.it. SUPPLEMENTARY INFORMATION: Details and data can be found at http://metaldetector.dsi.unifi.it/help.php.

Identifying cysteines and histidines in transition-metal-binding sites using support vector machines and neural networks.

Accurate predictions of metal-binding sites in proteins by using sequence as the only source of information can significantly help in the prediction of protein structure and function, genome annotation, and in the experimental determination of protein structure. Here, we introduce a method for identifying histidines and cysteines that participate in binding of several transition metals and iron complexes. The method predicts histidines as being in either of two states (free or metal bound) and cysteines in either of three states (free, metal bound, or in disulfide bridges). The method uses only sequence information by utilizing position-specific evolutionary profiles as well as more global descriptors such as protein length and amino acid composition. Our solution is based on a two-stage machine-learning approach. The first stage consists of a support vector machine trained to locally classify the binding state of single histidines and cysteines. The second stage consists of a bidirectional recurrent neural network trained to refine local predictions by taking into account dependencies among residues within the same protein. A simple finite state automaton is employed as a postprocessing in the second stage in order to enforce an even number of disulfide-bonded cysteines. We predict histidines and cysteines in transition-metal-binding sites at 73% precision and 61% recall. We observe significant differences in performance depending on the ligand (histidine or cysteine) and on the metal bound. We also predict cysteines participating in disulfide bridges at 86% precision and 87% recall. Results are compared to those that would be obtained by using expert information as represented by PROSITE motifs and, for disulfide bonds, to state-of-the-art methods.

Molecular modeling studies on CNG channel from bovine retinal rod: a structural model of the cyclic nucleotide-binding domain.

A dimeric model of the cyclic nucleotide-binding domain of the all-alpha homomeric cyclic nucleotide-gated channel from bovine retinal rod is constructed. The model, based on the structure of the fairly homologous catabolite gene activator protein (Weber and Steitz, J Mol Biol 1987;198:311-326), is obtained by use of comparative modeling and molecular dynamics simulations. Our model provides a structural basis for the experimentally measured difference in activity between cAMP and cGMP, as well as the different solvent accessibilities of GLY597 in the complex with cGMP, with cAMP and in the protein in free state. In addition, it provides support for the rearrangement of the domain C helix on ligand binding and releasing proposed by Matulef et al. (Neuron 1999;24:443-452).

Movement of the C-helix during the gating of cyclic nucleotide-gated channels.

Movements within the cyclic nucleotide-binding domain of cyclic nucleotide-gated channels are thought to underlie the initial phase of channel gating (Tibbs, G. R., D. T. Liu, B. G. Leypold, and S. A. Siegelbaum. 1998. J. Biol. Chem. 273:4497-4505; Zong, X., H. Zucker, F. Hofmann, and M. Biel. 1998. EMBO J. 17:353-362; Matulef, K., G. E. Flynn, and W. N. Zagotta. 1999. Neuron. 24:443-452; Paoletti, P., E. C. Young, and S. A. Siegelbaum. 1999. J. Gen. Physiol. 113:17-33; Johnson, J. P., and W. N. Zagotta. 2001. Nature. 412:917-921). To investigate these movements, cysteine mutation was performed on each of the 28 residues (Leu-583 to Asn-610), which span the agonist-binding domain of the alpha-subunit of the bovine rod cyclic nucleotide-gated channel. The effects of Cd(2+) ions, 2-trimethylammonioethylmethane thiosulfonate (MTSET) and copper phenanthroline (CuP) on channel activity were examined, in excised inside-out patches in the presence and in the absence of a saturating concentration of cGMP. The application of 100 microM Cd(2+) in the presence of saturating concentration of cGMP caused an irreversible and almost complete reduction of the current in mutant channels E594C, I600C, and L601C. In the absence of cGMP, the presence of 100 microM Cd(2+) caused a strong current reduction in all cysteine mutants from Asp-588 to Leu-607, with the exception of mutant channels A589C, M592C, M602C, K603C, and L606C. The selective effect of Cd(2+) ions was very similar to that observed when adding the oxidizing agent CuP to the bath medium, except for mutant channel G597C, where CuP caused a stronger current decrease (67 +/- 7%) than Cd(2+) (23 +/- 4%). In the absence of cGMP, MTSET caused a reduction of the current by >40% in mutant channels L607C, L601C, I600C, G597C, and E594C, whereas in the presence of cGMP only mutant channel L601C was affected. The application of MTSET protected many mutant channels from the effects of Cd(2+) and CuP. These results suggest that, when CNG channels are in the open state, residues from Asp-588 to Leu-607 are in an alpha-helical structure, homologous to the C-helix of the catabolite gene activator protein (Weber, I. T., and T. A. Steitz. 1987. J. Mol. Biol. 198:311-326). Furthermore, residues Glu-594, Gly-597, Ile-600, and Leu-601 of these helices belonging to two different subunits must be in close proximity. In the closed state the C-helices are in a different configuration and undergo significant fluctuations.