Systematic discovery and validation of T cell targets directed against oncogenic KRAS mutations.

Oncogenic mutations in KRAS can be recognized by T cells on specific class I human leukocyte antigen (HLA-I) molecules, leading to tumor control. To date, the discovery of T cell targets from KRAS mutations has relied on occasional T cell responses in patient samples or the use of transgenic mice. To overcome these limitations, we have developed a systematic target discovery and validation pipeline. We evaluate the presentation of mutant KRAS peptides on individual HLA-I molecules using targeted mass spectrometry and identify 13 unpublished KRASG12C/D/R/V mutation/HLA-I pairs and nine previously described pairs. We assess immunogenicity, generating T cell responses to nearly all targets. Using cytotoxicity assays, we demonstrate that KRAS-specific T cells and T cell receptors specifically recognize endogenous KRAS mutations. The discovery and validation of T cell targets from KRAS mutations demonstrate the potential for this pipeline to aid the development of immunotherapies for important cancer targets.

Sequence-based prediction of SARS-CoV-2 vaccine targets using a mass spectrometry-based bioinformatics predictor identifies immunogenic T cell epitopes.

BACKGROUND: The ongoing COVID-19 pandemic has created an urgency to identify novel vaccine targets for protective immunity against SARS-CoV-2. Early reports identify protective roles for both humoral and cell-mediated immunity for SARS-CoV-2. METHODS: We leveraged our bioinformatics binding prediction tools for human leukocyte antigen (HLA)-I and HLA-II alleles that were developed using mass spectrometry-based profiling of individual HLA-I and HLA-II alleles to predict peptide binding to diverse allele sets. We applied these binding predictors to viral genomes from the Coronaviridae family and specifically focused on T cell epitopes from SARS-CoV-2 proteins. We assayed a subset of these epitopes in a T cell induction assay for their ability to elicit CD8+ T cell responses. RESULTS: We first validated HLA-I and HLA-II predictions on Coronaviridae family epitopes deposited in the Virus Pathogen Database and Analysis Resource (ViPR) database. We then utilized our HLA-I and HLA-II predictors to identify 11,897 HLA-I and 8046 HLA-II candidate peptides which were highly ranked for binding across 13 open reading frames (ORFs) of SARS-CoV-2. These peptides are predicted to provide over 99% allele coverage for the US, European, and Asian populations. From our SARS-CoV-2-predicted peptide-HLA-I allele pairs, 374 pairs identically matched what was previously reported in the ViPR database, originating from other coronaviruses with identical sequences. Of these pairs, 333 (89%) had a positive HLA binding assay result, reinforcing the validity of our predictions. We then demonstrated that a subset of these highly predicted epitopes were immunogenic based on their recognition by specific CD8+ T cells in healthy human donor peripheral blood mononuclear cells (PBMCs). Finally, we characterized the expression of SARS-CoV-2 proteins in virally infected cells to prioritize those which could be potential targets for T cell immunity. CONCLUSIONS: Using our bioinformatics platform, we identify multiple putative epitopes that are potential targets for CD4+ and CD8+ T cells, whose HLA binding properties cover nearly the entire population. We also confirm that our binding predictors can predict epitopes eliciting CD8+ T cell responses from multiple SARS-CoV-2 proteins. Protein expression and population HLA allele coverage, combined with the ability to identify T cell epitopes, should be considered in SARS-CoV-2 vaccine design strategies and immune monitoring.

Utility of patient-derived lymphoblastoid cell lines as an ex vivo capecitabine sensitivity prediction model for breast cancer patients.

Capecitabine is commonly used in treating breast cancer; however, therapeutic response varies among patients and there is no clinically validated model to predict individual outcomes. Here, we investigated whether drug sensitivity quantified in ex vivo patients' blood-derived cell lines can predict response to capecitabine in vivo. Lymphoblastoid cell lines (LCLs) were established from a cohort of metastatic breast cancer patients (n = 53) who were prospectively monitored during treatment with single agent capecitabine at 2000 mg/m2/day. LCLs were treated with increasing concentrations of 5'-DFUR, a major capecitabine metabolite, to assess patients' ex vivo sensitivity to this drug. Subsequently, ex vivo phenotype was compared to observed patient disease response and drug induced-toxicities. We acquired an independent cohort of breast cancer cell lines and LCLs derived from the same donors from ATCC, compared their sensitivity to 5'-DFUR. As seen in the patient population, we observed large inter-individual variability in response to 5'-DFUR treatment in patient-derived LCLs. Patients whose LCLs were more sensitive to 5'-DFUR had a significantly longer median progression free survival (9-month vs 6-month, log rank p-value = 0.017). In addition, this significant positive correlation for 5'-DFUR sensitivity was replicated in an independent cohort of 8 breast cancer cell lines and LCLs derived from the same donor. Our data suggests that at least a portion of the individual sensitivity to capecitabine is shared between germline tissue and tumor tissue. It also supports the utility of patient-derived LCLs as a predictive model for capecitabine treatment efficacy in breast cancer patients.

Predicting Response to Histone Deacetylase Inhibitors Using High-Throughput Genomics.

BACKGROUND: Many disparate biomarkers have been proposed as predictors of response to histone deacetylase inhibitors (HDI); however, all have failed when applied clinically. Rather than this being entirely an issue of reproducibility, response to the HDI vorinostat may be determined by the additive effect of multiple molecular factors, many of which have previously been demonstrated. METHODS: We conducted a large-scale gene expression analysis using the Cancer Genome Project for discovery and generated another large independent cancer cell line dataset across different cancers for validation. We compared different approaches in terms of how accurately vorinostat response can be predicted on an independent out-of-batch set of samples and applied the polygenic marker prediction principles in a clinical trial. RESULTS: Using machine learning, the small effects that aggregate, resulting in sensitivity or resistance, can be recovered from gene expression data in a large panel of cancer cell lines.This approach can predict vorinostat response accurately, whereas single gene or pathway markers cannot. Our analyses recapitulated and contextualized many previous findings and suggest an important role for processes such as chromatin remodeling, autophagy, and apoptosis. As a proof of concept, we also discovered a novel causative role for CHD4, a helicase involved in the histone deacetylase complex that is associated with poor clinical outcome. As a clinical validation, we demonstrated that a common dose-limiting toxicity of vorinostat, thrombocytopenia, can be predicted (r = 0.55, P = .004) several days before it is detected clinically. CONCLUSION: Our work suggests a paradigm shift from single-gene/pathway evaluation to simultaneously evaluating multiple independent high-throughput gene expression datasets, which can be easily extended to other investigational compounds where similar issues are hampering clinical adoption.

MicroRNA biogenesis and cellular proliferation.

Given the fundamental roles of microRNAs (miRNAs) in physiological, developmental, and pathologic processes, we hypothesized that genes involved in miRNA biogenesis contribute to human complex traits. For 13 such genes, we evaluated the relationship between transcription and 2 classes of complex traits, namely cellular growth and sensitivity to various chemotherapeutic agents in a set of lymphoblastoid cell lines. We found a highly significant correlation between argonaute RNA-induced silencing complex catalytic component 2 (AGO2) expression and cellular growth rate (Bonferroni-adjusted P < 0.05), and report additional miRNA biogenesis genes with suggestive associations with either cellular growth rate or chemotherapeutic sensitivity. AGO2 expression was found to be correlated with multiple drug sensitivity phenotypes. Furthermore, small interfering RNA knockdown of AGO2 resulted in cellular growth inhibition in an ovarian cancer cell line (OVCAR-3), supporting the role of this miRNA biogenesis gene in cell proliferation in cancer cells. Expression quantitative trait loci mapping indicated that genetic variation (in the form of both single-nucleotide polymorphisms and copy number variations) that may regulate the expression of AGO2 can have downstream effects on cellular growth-dependent complex phenotypes.

Peripheral blood IRF1 expression as a marker for glucocorticoid sensitivity.

OBJECTIVE: Despite of the common usage of glucocorticoids (GCs), a significant portion of asthma patients exhibit GC insensitivity. This could be mediated by diverse mechanisms, including genomics. Recent work has suggested that measuring changes in gene expression may provide more predictive information about GC insensitivity than baseline gene expression alone, and that expression changes in peripheral blood may be reflective of those in the airway. METHODS: We performed in silico discovery using gene expression omnibus (GEO) data that evaluated GC effect on gene expression in multiple tissue types. Subsequently, candidate genes whose expression levels are affected by GC were examined in cell lines and in primary cells derived from human airway and blood. RESULTS: Through gene expression omnibus analysis, we identified interferon regulator factor 1 (IRF1), whose expression is affected by GC treatment in airway smooth muscle cells, normal human bronchial epithelial (NHBE) cells, and lymphoblastoid cell lines (LCLs). Significant IRF1 downregulation post GC exposure was confirmed in two cultured airway epithelial cell lines and primary NHBE cells (P<0.05). We observed large interindividual variation in GC-induced IRF1 expression changes among primary NHBE cells tested. Significant downregulation of IRF1 was also observed in six randomly selected LCLs (P<0.05), with variable degrees of downregulation among different samples. In peripheral blood mononuclear cells obtained from healthy volunteers, variable downregulation of IRF1 by GC was also shown. NFKB1, a gene whose expression is known to be downregulated by GC and the degree of downregulation of which is reflective of GC response, was used as a control in our study. IRF1 shows more consistent downregulation across tissue types when compared with NFKB1. CONCLUSION: Our results suggest that GC-induced IRF1 gene expression changes in peripheral blood could be used as a marker to reflect GC response in the airway.

Identifying and validating a combined mRNA and microRNA signature in response to imatinib treatment in a chronic myeloid leukemia cell line.

Imatinib, a targeted tyrosine kinase inhibitor, is the gold standard for managing chronic myeloid leukemia (CML). Despite its wide application, imatinib resistance occurs in 20-30% of individuals with CML. Multiple potential biomarkers have been identified to predict imatinib response; however, the majority of them remain externally uncorroborated. In this study, we set out to systematically identify gene/microRNA (miRNA) whose expression changes are related to imatinib response. Through a Gene Expression Omnibus search, we identified two genome-wide expression datasets that contain expression changes in response to imatinib treatment in a CML cell line (K562): one for mRNA and the other for miRNA. Significantly differentially expressed transcripts/miRNAs post imatinib treatment were identified from both datasets. Three additional filtering criteria were applied 1) miRbase/miRanda predictive algorithm; 2) opposite direction of imatinib effect for genes and miRNAs; and 3) literature support. These criteria narrowed our candidate gene-miRNA to a single pair: IL8 and miR-493-5p. Using PCR we confirmed the significant up-regulation and down-regulation of miR-493-5p and IL8 by imatinib treatment, respectively in K562 cells. In addition, IL8 expression was significantly down-regulated in K562 cells 24 hours after miR-493-5p mimic transfection (p = 0.002). Furthermore, we demonstrated significant cellular growth inhibition after IL8 inhibition through either gene silencing or by over-expression of miR-493-5p (p = 0.0005 and p = 0.001 respectively). The IL8 inhibition also further sensitized K562 cells to imatinib cytotoxicity (p < 0.0001). Our study combined expression changes in transcriptome and miRNA after imatinib exposure to identify a potential gene-miRNA pair that is a critical target in imatinib response. Experimental validation supports the relationships between IL8 and miR-493-5p and between this gene-miRNA pair and imatinib sensitivity in a CML cell line. Our data suggests integrative analysis of multiple omic level data may provide new insight into biomarker discovery as well as mechanisms of imatinib resistance.

Integrative analyses of genetic variation, epigenetic regulation, and the transcriptome to elucidate the biology of platinum sensitivity.

BACKGROUND: Using genome-wide genetic, gene expression, and microRNA expression (miRNA) data, we developed an integrative approach to investigate the genetic and epigenetic basis of chemotherapeutic sensitivity. RESULTS: Through a sequential multi-stage framework, we identified genes and miRNAs whose expression correlated with platinum sensitivity, mapped these to genomic loci as quantitative trait loci (QTLs), and evaluated the associations between these QTLs and platinum sensitivity. A permutation analysis showed that top findings from our approach have a much lower false discovery rate compared to those from a traditional GWAS of drug sensitivity. Our approach identified five SNPs associated with 10 miRNAs and the expression level of 15 genes, all of which were associated with carboplatin sensitivity. Of particular interest was one SNP (rs11138019), which was associated with the expression of both miR-30d and the gene ABCD2, which were themselves correlated with both carboplatin and cisplatin drug-specific phenotype in the HapMap samples. Functional study found that knocking down ABCD2 in vitro led to increased apoptosis in ovarian cancer cell line SKOV3 after cisplatin treatment. Over-expression of miR-30d in vitro caused a decrease in ABCD2 expression, suggesting a functional relationship between the two. CONCLUSIONS: We developed an integrative approach to the investigation of the genetic and epigenetic basis of human complex traits. Our approach outperformed standard GWAS and provided hints at potential biological function. The relationships between ABCD2 and miR-30d, and ABCD2 and platin sensitivity were experimentally validated, suggesting a functional role of ABCD2 and miR-30d in sensitivity to platinating agents.

The impact of microRNA expression on cellular proliferation.

As an important class of non-coding regulatory RNAs, microRNAs (miRNAs) play a key role in a range of biological processes. These molecules serve as post-transcriptional regulators of gene expression and their regulatory activity has been implicated in disease pathophysiology and pharmacological traits. We sought to investigate the impact of miRNAs on cellular proliferation to gain insight into the molecular basis of complex traits that depend on cellular growth, including, most prominently, cancer. We examined the relationship between miRNA expression and intrinsic cellular growth (iGrowth) in the HapMap lymphoblastoid cell lines derived from individuals of different ethnic backgrounds. We found a substantial enrichment for miRNAs (53 miRNAs, FDR < 0.05) correlated with cellular proliferation in pooled CEU (Caucasian of northern and western European descent) and YRI (individuals from Ibadan, Nigeria) samples. Specifically, 119 miRNAs (59 %) were significantly correlated with iGrowth in YRI; of these miRNAs, 18 were correlated with iGrowth in CEU. To gain further insight into the effect of miRNAs on cellular proliferation in cancer, we showed that over-expression of miR-22, one of the top iGrowth-associated miRNAs, leads to growth inhibition in an ovarian cancer cell line (SKOV3). Furthermore, over-expression of miR-22 down-regulates the expression of its target genes (MXI1 and SLC25A37) in this ovarian cancer cell line, highlighting an miRNA-mediated regulatory network potentially important for cellular proliferation. Importantly, our study identified miRNAs that can be used as molecular targets in cancer therapy.

Genetic variation is the major determinant of individual differences in leukocyte endothelial adhesion.

OBJECTIVE: To determine the genetic contribution to leukocyte endothelial adhesion. METHODS: Leukocyte endothelial adhesion was assessed through a novel cell-based assay using human lymphoblastoid cell lines. A high-throughput screening method was developed to evaluate the inter-individual variability in leukocyte endothelial adhesion using lymphoblastoid cell lines derived from different donors. To assess heritability, ninety-two lymphoblastoid cell lines derived from twenty-three monozygotic twin pairs and twenty-three sibling pairs were compared. These lymphoblastoid cell lines were plated with the endothelial cell line EA.hy926 and labeled with Calcein AM dye. Fluorescence was assessed to determine endothelial cell adhesion to each lymphoblastoid cell line. Intra-pair similarity was determined for monozygotic twins and siblings using Pearson pairwise correlation coefficients. RESULTS: A leukocyte endothelial adhesion assay for lymphoblastoid cell lines was developed and optimized (CV = 8.68, Z'-factor = 0.67, SNR = 18.41). A higher adhesion correlation was found between the twins than that between the siblings. Intra-pair similarity for leukocyte endothelial adhesion in monozygotic twins was 0.60 compared to 0.25 in the siblings. The extent to which these differences are attributable to underlying genetic factors was quantified and the heritability of leukocyte endothelial adhesion was calculated to be 69.66% (p-value<0.0001). CONCLUSIONS: There is a heritable component to leukocyte endothelial adhesion. Underlying genetic predisposition plays a significant role in inter-individual variability of leukocyte endothelial adhesion.