Immune features are associated with response to neoadjuvant chemo-immunotherapy for muscle-invasive bladder cancer.

Neoadjuvant cisplatin-based chemotherapy is standard of care for muscle-invasive bladder cancer (MIBC). Immune checkpoint inhibition (ICI) alone, and ICI in combination with chemotherapy, have demonstrated promising pathologic response (

LENS: Landscape of Effective Neoantigens Software.

MOTIVATION: Elimination of cancer cells by T cells is a critical mechanism of anti-tumor immunity and cancer immunotherapy response. T cells recognize cancer cells by engagement of T cell receptors with peptide epitopes presented by major histocompatibility complex molecules on the cancer cell surface. Peptide epitopes can be derived from antigen proteins coded for by multiple genomic sources. Bioinformatics tools used to identify tumor-specific epitopes via analysis of DNA and RNA-sequencing data have largely focused on epitopes derived from somatic variants, though a smaller number have evaluated potential antigens from other genomic sources. RESULTS: We report here an open-source workflow utilizing the Nextflow DSL2 workflow manager, Landscape of Effective Neoantigens Software (LENS), which predicts tumor-specific and tumor-associated antigens from single nucleotide variants, insertions and deletions, fusion events, splice variants, cancer-testis antigens, overexpressed self-antigens, viruses, and endogenous retroviruses. The primary advantage of LENS is that it expands the breadth of genomic sources of discoverable tumor antigens using genomics data. Other advantages include modularity, extensibility, ease of use, and harmonization of relative expression level and immunogenicity prediction across multiple genomic sources. We present an analysis of 115 acute myeloid leukemia samples to demonstrate the utility of LENS. We expect LENS will be a valuable platform and resource for T cell epitope discovery bioinformatics, especially in cancers with few somatic variants where tumor-specific epitopes from alternative genomic sources are an elevated priority. AVAILABILITY AND IMPLEMENTATION: More information about LENS, including code, workflow documentation, and instructions, can be found at (https://gitlab.com/landscape-of-effective-neoantigens-software).

Associations of minor histocompatibility antigens with outcomes following allogeneic hematopoietic cell transplantation.

The role of minor histocompatibility antigens (mHAs) in mediating graft versus leukemia and graft versus host disease (GvHD) following allogeneic hematopoietic cell transplantation (alloHCT) is recognized but not well-characterized. By implementing improved methods for mHA prediction in two large patient cohorts, this study aimed to comprehensively explore the role of mHAs in alloHCT by analyzing whether (1) the number of predicted mHAs, or (2) individual mHAs are associated with clinical outcomes. The study population consisted of 2249 donor-recipient pairs treated for acute myeloid leukemia and myelodysplastic syndrome with alloHCT. A Cox proportional hazard model showed that patients with a class I mHA count greater than the population median had an increased hazard of GvHD mortality (hazard ratio [HR] = 1.39, 95% confidence interval [CI] = 1.01, 1.77, p = .046). Competing risk analyses identified the class I mHAs DLRCKYISL (GSTP), WEHGPTSLL (CRISPLD2), and STSPTTNVL (SERPINF2) were associated with increased GVHD mortality (HR = 2.84, 95% CI = 1.52, 5.31, p = .01), decreased leukemia-free survival (LFS) (HR = 1.94, 95% CI = 1.27, 2.95, p = .044), and increased disease-related mortality (DRM) (HR = 2.32, 95% CI = 1.5, 3.6, p = .008), respectively. One class II mHA YQEIAAIPSAGRERQ (TACC2) was associated with increased risk of treatment-related mortality (TRM) (HR = 3.05, 95% CI = 1.75, 5.31, p = .02). WEHGPTSLL and STSPTTNVL were both present within HLA haplotype B*40:01-C*03:04 and showed a positive dose-response relationship with increased all-cause mortality and DRM and decreased LFS, indicating these two mHAs contribute to the risk of mortality in an additive manner. Our study reports the first large-scale investigation of the associations of predicted mHA peptides with clinical outcomes following alloHCT.

Evaluation of tumor antigen-specific antibody responses in patients with metastatic triple negative breast cancer treated with cyclophosphamide and pembrolizumab.

The role of B cells in antitumor immunity is becoming increasingly appreciated, as B cell populations have been associated with response to immune checkpoint blockade (ICB) in patients with breast cancer and murine models of breast cancer. Deeper understanding of antibody responses to tumor antigens is needed to clarify the function of B cells in determining response to immunotherapy. We evaluated tumor antigen-specific antibody responses in patients with metastatic triple negative breast cancer treated with pembrolizumab following low-dose cyclophosphamide therapy using computational linear epitope prediction and custom peptide microarrays. We found that a minority of predicted linear epitopes were associated with antibody signal, and signal was associated with both neoepitopes and self-peptides. No association was observed between signal presence and subcellular localization or RNA expression of parent proteins. Patient-specific patterns of antibody signal boostability were observed that were independent of clinical response. Intriguingly, measures of cumulative antibody signal intensity relative to immunotherapy treatment showed that the one complete responder in the trial had the greatest increase in total antibody signal, which supports a potential association between ICB-dependent antibody boosting and clinical response. The antibody boost in the complete responder was largely driven by increased levels of IgG specific to a sequence of N-terminal residues in native Epidermal Growth Factor Receptor Pathway Substrate 8 (EPS8) protein, a known oncogene in several cancer types including breast cancer. Structural protein prediction showed that the targeted epitope of EPS8 was in a region of the protein with mixed linear/helical structure, and that this region was solvent-exposed and not predicted to bind to interacting macromolecules. This study highlights the potential importance of the humoral immune response targeting neoepitopes as well as self epitopes in shaping clinical response to immunotherapy.

Shared graft-versus-leukemia minor histocompatibility antigens in DISCOVeRY-BMT.

T-cell responses to minor histocompatibility antigens (mHAs) mediate graft-versus-leukemia (GVL) effects and graft-versus-host disease (GVHD) in allogeneic hematopoietic cell transplantation. Therapies that boost T-cell responses improve allogeneic hematopoietic cell transplant (alloHCT) efficacy but are limited by concurrent increases in the incidence and severity of GVHD. mHAs with expression restricted to hematopoietic tissue (GVL mHAs) are attractive targets for driving GVL without causing GVHD. Prior work to identify mHAs has focused on a small set of mHAs or population-level single-nucleotide polymorphism-association studies. We report the discovery of a large set of novel GVL mHAs based on predicted immunogenicity, tissue expression, and degree of sharing among donor-recipient pairs (DRPs) in the DISCOVeRY-BMT data set of 3231 alloHCT DRPs. The total number of predicted mHAs varied by HLA allele, and the total number and number of each class of mHA significantly differed by recipient genomic ancestry group. From the pool of predicted mHAs, we identified the smallest sets of GVL mHAs needed to cover 100% of DRPs with a given HLA allele. We used mass spectrometry to search for high-population frequency mHAs for 3 common HLA alleles. We validated 24 predicted novel GVL mHAs that are found cumulatively within 98.8%, 60.7%, and 78.9% of DRPs within DISCOVeRY-BMT that express HLA-A∗02:01, HLA-B∗35:01, and HLA-C∗07:02, respectively. We confirmed the immunogenicity of an example novel mHA via T-cell coculture with peptide-pulsed dendritic cells. This work demonstrates that the identification of shared mHAs is a feasible and promising technique for expanding mHA-targeting immunotherapeutics.

Comprehensive Analysis of the Immunogenomics of Triple-Negative Breast Cancer Brain Metastases From LCCC1419.

Background: Triple negative breast cancer (TNBC) is an aggressive variant of breast cancer that lacks the expression of estrogen and progesterone receptors (ER and PR) and HER2. Nearly 50% of patients with advanced TNBC will develop brain metastases (BrM), commonly with progressive extracranial disease. Immunotherapy has shown promise in the treatment of advanced TNBC; however, the immune contexture of BrM remains largely unknown. We conducted a comprehensive analysis of TNBC BrM and matched primary tumors to characterize the genomic and immune landscape of TNBC BrM to inform the development of immunotherapy strategies in this aggressive disease. Methods: Whole-exome sequencing (WES) and RNA sequencing were conducted on formalin-fixed, paraffin-embedded samples of BrM and primary tumors of patients with clinical TNBC (n = 25, n = 9 matched pairs) from the LCCC1419 biobank at UNC-Chapel Hill. Matched blood was analyzed by DNA sequencing as a comparison for tumor WES for the identification of somatic variants. A comprehensive genomics assessment, including mutational and copy number alteration analyses, neoantigen prediction, and transcriptomic analysis of the tumor immune microenvironment were performed. Results: Primary and BrM tissues were confirmed as TNBC (23/25 primaries, 16/17 BrM) by immunohistochemistry and of the basal intrinsic subtype (13/15 primaries and 16/19 BrM) by PAM50. Compared to primary tumors, BrM demonstrated a higher tumor mutational burden. TP53 was the most frequently mutated gene and was altered in 50% of the samples. Neoantigen prediction showed elevated cancer testis antigen- and endogenous retrovirus-derived MHC class I-binding peptides in both primary tumors and BrM and predicted that single-nucleotide variant (SNV)-derived peptides were significantly higher in BrM. BrM demonstrated a reduced immune gene signature expression, although a signature associated with fibroblast-associated wound healing was elevated in BrM. Metrics of T and B cell receptor diversity were also reduced in BrM. Conclusions: BrM harbored higher mutational burden and SNV-derived neoantigen expression along with reduced immune gene signature expression relative to primary TNBC. Immune signatures correlated with improved survival, including T cell signatures. Further research will expand these findings to other breast cancer subtypes in the same biobank. Exploration of immunomodulatory approaches including vaccine applications and immune checkpoint inhibition to enhance anti-tumor immunity in TNBC BrM is warranted.

NeoSplice: a bioinformatics method for prediction of splice variant neoantigens.

Motivation: Splice variant neoantigens are a potential source of tumor-specific antigen (TSA) that are shared between patients in a variety of cancers, including acute myeloid leukemia. Current tools for genomic prediction of splice variant neoantigens demonstrate promise. However, many tools have not been well validated with simulated and/or wet lab approaches, with no studies published that have presented a targeted immunopeptidome mass spectrometry approach designed specifically for identification of predicted splice variant neoantigens. Results: In this study, we describe NeoSplice, a novel computational method for splice variant neoantigen prediction based on (i) prediction of tumor-specific k-mers from RNA-seq data, (ii) alignment of differentially expressed k-mers to the splice graph and (iii) inference of the variant transcript with MHC binding prediction. NeoSplice demonstrates high sensitivity and precision (>80% on average across all splice variant classes) through in silico simulated RNA-seq data. Through mass spectrometry analysis of the immunopeptidome of the K562.A2 cell line compared against a synthetic peptide reference of predicted splice variant neoantigens, we validated 4 of 37 predicted antigens corresponding to 3 of 17 unique splice junctions. Lastly, we provide a comparison of NeoSplice against other splice variant prediction tools described in the literature. NeoSplice provides a well-validated platform for prediction of TSA vaccine targets for future cancer antigen vaccine studies to evaluate the clinical efficacy of splice variant neoantigens. Availability and implementation: https://github.com/Benjamin-Vincent-Lab/NeoSplice. Supplementary information: Supplementary data are available at Bioinformatics Advances online.

Landscape and selection of vaccine epitopes in SARS-CoV-2.

BACKGROUND: Early in the pandemic, we designed a SARS-CoV-2 peptide vaccine containing epitope regions optimized for concurrent B cell, CD4+ T cell, and CD8+ T cell stimulation. The rationale for this design was to drive both humoral and cellular immunity with high specificity while avoiding undesired effects such as antibody-dependent enhancement (ADE). METHODS: We explored the set of computationally predicted SARS-CoV-2 HLA-I and HLA-II ligands, examining protein source, concurrent human/murine coverage, and population coverage. Beyond MHC affinity, T cell vaccine candidates were further refined by predicted immunogenicity, sequence conservation, source protein abundance, and coverage of high frequency HLA alleles. B cell epitope regions were chosen from linear epitope mapping studies of convalescent patient serum, followed by filtering for surface accessibility, sequence conservation, spatial localization near functional domains of the spike glycoprotein, and avoidance of glycosylation sites. RESULTS: From 58 initial candidates, three B cell epitope regions were identified. From 3730 (MHC-I) and 5045 (MHC-II) candidate ligands, 292 CD8+ and 284 CD4+ T cell epitopes were identified. By combining these B cell and T cell analyses, as well as a manufacturability heuristic, we proposed a set of 22 SARS-CoV-2 vaccine peptides for use in subsequent murine studies. We curated a dataset of ~ 1000 observed T cell epitopes from convalescent COVID-19 patients across eight studies, showing 8/15 recurrent epitope regions to overlap with at least one of our candidate peptides. Of the 22 candidate vaccine peptides, 16 (n = 10 T cell epitope optimized; n = 6 B cell epitope optimized) were manually selected to decrease their degree of sequence overlap and then synthesized. The immunogenicity of the synthesized vaccine peptides was validated using ELISpot and ELISA following murine vaccination. Strong T cell responses were observed in 7/10 T cell epitope optimized peptides following vaccination. Humoral responses were deficient, likely due to the unrestricted conformational space inhabited by linear vaccine peptides. CONCLUSIONS: Overall, we find our selection process and vaccine formulation to be appropriate for identifying T cell epitopes and eliciting T cell responses against those epitopes. Further studies are needed to optimize prediction and induction of B cell responses, as well as study the protective capacity of predicted T and B cell epitopes.

Landscape and Selection of Vaccine Epitopes in SARS-CoV-2.

There is an urgent need for a vaccine with efficacy against SARS-CoV-2. We hypothesize that peptide vaccines containing epitope regions optimized for concurrent B cell, CD4+ T cell, and CD8+ T cell stimulation would drive both humoral and cellular immunity with high specificity, potentially avoiding undesired effects such as antibody-dependent enhancement (ADE). Additionally, such vaccines can be rapidly manufactured in a distributed manner. In this study, we combine computational prediction of T cell epitopes, recently published B cell epitope mapping studies, and epitope accessibility to select candidate peptide vaccines for SARS-CoV-2. We begin with an exploration of the space of possible T cell epitopes in SARS-CoV-2 with interrogation of predicted HLA-I and HLA-II ligands, overlap between predicted ligands, protein source, as well as concurrent human/murine coverage. Beyond MHC affinity, T cell vaccine candidates were further refined by predicted immunogenicity, viral source protein abundance, sequence conservation, coverage of high frequency HLA alleles and co-localization of CD4+ and CD8+ T cell epitopes. B cell epitope regions were chosen from linear epitope mapping studies of convalescent patient serum, followed by filtering to select regions with surface accessibility, high sequence conservation, spatial localization near functional domains of the spike glycoprotein, and avoidance of glycosylation sites. From 58 initial candidates, three B cell epitope regions were identified. By combining these B cell and T cell analyses, as well as a manufacturability heuristic, we propose a set of SARS-CoV-2 vaccine peptides for use in subsequent murine studies and clinical trials.

Recent progress and open questions in Drosophila dosage compensation.

Sexual dimorphism is observed in many traits across diverse taxa, and often it is quite extreme. Within a species, individuals of opposing sex can appear strikingly different, reflecting differences at the molecular level that may be similarly striking. Among the most extreme cases of such molecular sexual dimorphism is the quantity of sex chromosomes that each sex possesses. Hemizygous sex chromosomes are common to many species, and various mechanisms have evolved to regulate transcriptional activity to ensure appropriate sex chromosome-to-autosome gene expression stoichiometry. Among the most thoroughly investigated of these mechanisms is Drosophila melanogaster's male-specific lethal (MSL) complex-mediated dosage compensation. In Drosophila, the male X chromosome transcription is upregulated approximately two-fold in somatic tissues to counterbalance the effects of sex chromosome hemizygosity on transcript abundance. Despite dramatic advances in our understanding of the Drosophila dosage compensation, many questions remain unanswered, and our understanding of its molecular underpinnings remains incomplete. In this review, we synthesize recent progress in the field as a means to highlight open questions, including how the MSL complex targets the X chromosome, how dosage compensation has shaped evolution of X-linked genes, and the degree to which MSL complex-mediated dosage compensation varies in activity across somatic tissues.

No evidence for a global male-specific lethal complex-mediated dosage compensation contribution to the demasculinization of the Drosophila melanogaster X chromosome.

In Drosophila melanogaster males, the expression of X-linked genes is regulated by mechanisms that operate on a chromosomal scale. One such mechanism, male-specific lethal complex-dependent X-linked dosage compensation, is thought to broadly enhance the expression of male X-linked genes through two-fold transcriptional upregulation. The evolutionary consequences of this form of dosage compensation are not well understood, particularly with regard to genes more highly expressed in males. It has been observed the X chromosome arrangement of these male-biased genes is non-random, consistent with what one might expect if there is a selective advantage for male-biased genes to avoid dosage compensation. Separately, it has been noted that the male-specific lethal complex and its dosage compensation mechanism appear absent in some male tissues, thus providing a control for the selection hypothesis. Here we utilized publicly available datasets to reassess the arrangement of X-linked male-biased expressed genes after accounting for expression in tissues not dosage compensated by the male-specific lethal complex. Our results do not corroborate previous observations supporting organismal-wide detrimental effects by dosage compensation on X-linked male-biased expressed genes. We instead find no evidence that dosage compensation has played a role in the arrangement of dosage compensated male-biased genes on the X chromosome.

Conservation and sex-specific splicing of the transformer gene in the calliphorids Cochliomyia hominivorax, Cochliomyia macellaria and Lucilia sericata.

Transformer (TRA) promotes female development in several dipteran species including the Australian sheep blowfly Lucilia cuprina, the Mediterranean fruit fly, housefly and Drosophila melanogaster. tra transcripts are sex-specifically spliced such that only the female form encodes full length functional protein. The presence of six predicted TRA/TRA2 binding sites in the sex-specific female intron of the L. cuprina gene suggested that tra splicing is auto-regulated as in medfly and housefly. With the aim of identifying conserved motifs that may play a role in tra sex-specific splicing, here we have isolated and characterized the tra gene from three additional blowfly species, L. sericata, Cochliomyia hominivorax and C. macellaria. The blowfly adult male and female transcripts differ in the choice of splice donor site in the first intron, with males using a site downstream of the site used in females. The tra genes all contain a single TRA/TRA2 site in the male exon and a cluster of four to five sites in the male intron. However, overall the sex-specific intron sequences are poorly conserved in closely related blowflies. The most conserved regions are around the exon/intron junctions, the 3' end of the intron and near the cluster of TRA/TRA2 sites. We propose a model for sex specific regulation of tra splicing that incorporates the conserved features identified in this study. In L. sericata embryos, the male tra transcript was first detected at around the time of cellular blastoderm formation. RNAi experiments showed that tra is required for female development in L. sericata and C. macellaria. The isolation of the tra gene from the New World screwworm fly C. hominivorax, a major livestock pest, will facilitate the development of a "male-only" strain for genetic control programs.