Suchit Jhunjhunwala

Predicting immunogenic tumour mutations by combining mass spectrometry and exome sequencing

Human tumours typically harbour a remarkable number of somatic mutations. If presented on major histocompatibility complex class I molecules (MHCI), peptides containing these mutations could potentially be immunogenic as they should be recognized as ‘non-self’ neo-antigens by the adaptive immune system. Recent work has confirmed that mutant peptides can serve as T-cell epitopes. However, few mutant epitopes have been described because their discovery required the laborious screening of patient tumour-infiltrating lymphocytes for their ability to recognize antigen libraries constructed following tumour exome sequencing. We sought to simplify the discovery of immunogenic mutant peptides by characterizing their general properties. We developed an approach that combines whole-exome and transcriptome sequencing analysis with mass spectrometry to identify neo-epitopes in two widely used murine tumour models. Of the >1,300 amino acid changes identified, ∼13% were predicted to bind MHCI, a small fraction of which were confirmed by mass spectrometry. The peptides were then structurally modelled bound to MHCI. Mutations that were solvent-exposed and therefore accessible to T-cell antigen receptors were predicted to be immunogenic. Vaccination of mice confirmed the approach, with each predicted immunogenic peptide yielding therapeutically active T-cell responses. The predictions also enabled the generation of peptide–MHCI dextramers that could be used to monitor the kinetics and distribution of the anti-tumour T-cell response before and after vaccination. These findings indicate that a suitable prediction algorithm may provide an approach for the pharmacodynamic monitoring of T-cell responses as well as for the development of personalized vaccines in cancer patients.

TRPS1 Targeting by miR-221/222 Promotes the Epithelial-to-Mesenchymal Transition in Breast Cancer

The microRNAs miR-221 and miR-222 promote a phenotype associated with metastasis and are found in a clinically aggressive form of breast cancer. Parsing Breast Cancer Subtype with MicroRNAs MicroRNAs (miRNAs), short noncoding RNAs that bind to and silence target mRNAs, have emerged as playing crucial regulatory roles not only in normal cellular processes but also in pathological conditions, such as cancer. Stinson et al. analyzed miRNA expression in different types of human breast cancer and found that miR-221 and miR-222 (miR-221/222) abundance was increased in the clinically aggressive basal-like subtype compared to the less aggressive luminal subtype. They determined that signaling through the epidermal growth factor receptor (EGFR)–RAS–extracellular signal–regulated kinase (ERK) pathway increased miR-221/222 transcription, and they defined a transcriptional regulatory pathway through which miR-221/222 promoted a phenotype associated with cancer cell invasion and metastasis. Their data suggest that combining inhibition of the EGFR-RAS-ERK pathway with standard chemotherapy could, by limiting miR-221/222 production, provide a strategy to combat metastasis in the basal-like subtype of breast cancer. The basal-like subtype of breast cancer has an aggressive clinical behavior compared to that of the luminal subtype. We identified the microRNAs (miRNAs) miR-221 and miR-222 (miR-221/222) as basal-like subtype–specific miRNAs and showed that expression of miR-221/222 decreased expression of epithelial-specific genes and increased expression of mesenchymal-specific genes, and increased cell migration and invasion in a manner characteristic of the epithelial-to-mesenchymal transition (EMT). The transcription factor FOSL1 (also known as Fra-1), which is found in basal-like breast cancers but not in the luminal subtype, stimulated the transcription of miR-221/222, and the abundance of these miRNAs decreased with inhibition of the epidermal growth factor receptor (EGFR) or MEK (mitogen-activated or extracellular signal–regulated protein kinase kinase), placing miR-221/222 downstream of the RAS pathway. Furthermore, miR-221/222–mediated reduction in E-cadherin abundance depended on their targeting the 3′ untranslated region of the GATA family transcriptional repressor TRPS1 (tricho-rhino-phalangeal syndrome type 1), which inhibited EMT by decreasing ZEB2 (zinc finger E-box–binding homeobox2) expression. We conclude that by promoting EMT, miR-221/222 may contribute to the more aggressive clinical behavior of basal-like breast cancers.

The effects of hepatitis B virus integration into the genomes of hepatocellular carcinoma patients

Hepatitis B virus (HBV) infection is a leading risk factor for hepatocellular carcinoma (HCC). HBV integration into the host genome has been reported, but its scale, impact and contribution to HCC development is not clear. Here, we sequenced the tumor and nontumor genomes (>80× coverage) and transcriptomes of four HCC patients and identified 255 HBV integration sites. Increased sequencing to 240× coverage revealed a proportionally higher number of integration sites. Clonal expansion of HBV-integrated hepatocytes was found specifically in tumor samples. We observe a diverse collection of genomic perturbations near viral integration sites, including direct gene disruption, viral promoter-driven human transcription, viral-human transcript fusion, and DNA copy number alteration. Thus, we report the most comprehensive characterization of HBV integration in hepatocellular carcinoma patients. Such widespread random viral integration will likely increase carcinogenic opportunities in HBV-infected individuals.

Genome and transcriptome sequencing of lung cancers reveal diverse mutational and splicing events

Lung cancer is a highly heterogeneous disease in terms of both underlying genetic lesions and response to therapeutic treatments. We performed deep whole-genome sequencing and transcriptome sequencing on 19 lung cancer cell lines and three lung tumor/normal pairs. Overall, our data show that cell line models exhibit similar mutation spectra to human tumor samples. Smoker and never-smoker cancer samples exhibit distinguishable patterns of mutations. A number of epigenetic regulators, including KDM6A, ASH1L, SMARCA4, and ATAD2, are frequently altered by mutations or copy number changes. A systematic survey of splice-site mutations identified 106 splice site mutations associated with cancer specific aberrant splicing, including mutations in several known cancer-related genes. RAC1b, an isoform of the RAC1 GTPase that includes one additional exon, was found to be preferentially up-regulated in lung cancer. We further show that its expression is significantly associated with sensitivity to a MAP2K (MEK) inhibitor PD-0325901. Taken together, these data present a comprehensive genomic landscape of a large number of lung cancer samples and further demonstrate that cancer-specific alternative splicing is a widespread phenomenon that has potential utility as therapeutic biomarkers. The detailed characterizations of the lung cancer cell lines also provide genomic context to the vast amount of experimental data gathered for these lines over the decades, and represent highly valuable resources for cancer biology.

miR-221/222 Targeting of Trichorhinophalangeal 1 (TRPS1) Promotes Epithelial-to-Mesenchymal Transition in Breast Cancer

MicroRNAs miR-221 and miR-222 are associated with a clinically aggressive form of breast cancer and promote epithelial-to-mesenchymal transition. Compared with the luminal subtype, the basal-like subtype of breast cancer has an aggressive clinical behavior, but the reasons for this difference between the two subtypes are poorly understood. We identified microRNAs (miRNAs) miR-221 and miR-222 (miR-221/222) as basal-like subtype-specific miRNAs that decrease expression of epithelial-specific genes and increase expression of mesenchymal-specific genes. In addition, expression of these miRNAs increased cell migration and invasion, which collectively are characteristics of the epithelial-to-mesenchymal transition (EMT). The basal-like transcription factor FOSL1 (also known as Fra-1) directly stimulated the transcription of miR-221/222, and the abundance of these miRNAs decreased with inhibition of MEK (mitogen-activated or extracellular signal–regulated protein kinase kinase), placing miR-221/222 downstream of the RAS pathway. The miR-221/222–mediated reduction in E-cadherin abundance depended on their targeting of the 3′ untranslated region (3′UTR) of TRPS1 (trichorhinophalangeal syndrome type 1), which is a member of the GATA family of transcriptional repressors. TRPS1 inhibited EMT by directly repressing expression of ZEB2 (Zinc finger E-box–binding homeobox 2). Therefore, miR-221/222 may contribute to the aggressive clinical behavior of basal-like breast cancers.

TGFβ attenuates tumour response to PD-L1 blockade by contributing to exclusion of T cells

Therapeutic antibodies that block the programmed death-1 (PD-1)–programmed death-ligand 1 (PD-L1) pathway can induce robust and durable responses in patients with various cancers, including metastatic urothelial cancer. However, these responses only occur in a subset of patients. Elucidating the determinants of response and resistance is key to improving outcomes and developing new treatment strategies. Here we examined tumours from a large cohort of patients with metastatic urothelial cancer who were treated with an anti-PD-L1 agent (atezolizumab) and identified major determinants of clinical outcome. Response to treatment was associated with CD8+ T-effector cell phenotype and, to an even greater extent, high neoantigen or tumour mutation burden. Lack of response was associated with a signature of transforming growth factor β (TGFβ) signalling in fibroblasts. This occurred particularly in patients with tumours, which showed exclusion of CD8+ T cells from the tumour parenchyma that were instead found in the fibroblast- and collagen-rich peritumoural stroma; a common phenotype among patients with metastatic urothelial cancer. Using a mouse model that recapitulates this immune-excluded phenotype, we found that therapeutic co-administration of TGFβ-blocking and anti-PD-L1 antibodies reduced TGFβ signalling in stromal cells, facilitated T-cell penetration into the centre of tumours, and provoked vigorous anti-tumour immunity and tumour regression. Integration of these three independent biological features provides the best basis for understanding patient outcome in this setting and suggests that TGFβ shapes the tumour microenvironment to restrain anti-tumour immunity by restricting T-cell infiltration.

A novel anti-CD22 anthracycline-based antibody-drug conjugate (ADC) that overcomes resistance to auristatin based ADCs

Purpose: We are interested in identifying mechanisms of resistance to the current generation of antibody–drug conjugates (ADC) and developing ADCs that can overcome this resistance. Experimental Design: Pinatuzumab vedotin (anti-CD22-vc-MMAE) and polatuzumab vedotin (anti-CD79b-vc-MMAE) are ADCs that contain the microtubule inhibitor monomethyl auristatin E (MMAE) attached to the antibody by the protease-cleavable linker maleimidocaproyl-valine-citrulline-p-aminobenzoyloxycarbonyl (MC-vc-PAB). Early clinical trial data suggest that these ADCs have promising efficacy for the treatment of non-Hodgkin lymphoma (NHL); however, some patients do not respond or become resistant to the ADCs. Anthracyclines are very effective in NHL, but ADCs containing the anthracycline doxorubicin were not clinically efficacious probably due to the low drug potency and inadequate linker technology. The anthracycline analogue PNU-159682 is thousands of times more cytotoxic than doxorubicin, so we used it to develop a new class of ADCs. We used the same MC-vc-PAB linker and antibody in pinatuzumab vedotin but replaced the MMAE with a derivative of PNU-159682 to make anti-CD22-NMS249 and tested it for in vivo efficacy in xenograft tumors resistant to MMAE-based ADCs. Results: We derived cell lines from in vivo xenograft tumors that were made resistant to anti-CD22-vc-MMAE and anti-CD79b-vc-MMAE. We identified P-gp (ABCB1/MDR1) as the major driver of resistance to the vc-MMAE–based conjugates. Anti-CD22-NMS249 was at least as effective as anti-CD22-vc-MMAE in xenograft models of the parental cell lines and maintained its efficacy in the resistant cell lines. Conclusions: These studies provide proof of concept for an anthracycline-based ADC that could be used to treat B-cell malignancies that are resistant to vc-MMAE conjugates. Clin Cancer Res; 21(14); 3298–306. ©2015 AACR.

Integrated exome and transcriptome sequencing reveals ZAK isoform usage in gastric cancer

Gastric cancer is the second leading cause of worldwide cancer mortality, yet the underlying genomic alterations remain poorly understood. Here we perform exome and transcriptome sequencing and SNP array assays to characterize 51 primary gastric tumours and 32 cell lines. Meta-analysis of exome data and previously published data sets reveals 24 significantly mutated genes in microsatellite stable (MSS) tumours and 16 in microsatellite instable (MSI) tumours. Over half the patients in our collection could potentially benefit from targeted therapies. We identify 55 splice site mutations accompanied by aberrant splicing products, in addition to mutation-independent differential isoform usage in tumours. ZAK kinase isoform TV1 is preferentially upregulated in gastric tumours and cell lines relative to normal samples. This pattern is also observed in colorectal, bladder and breast cancers. Overexpression of this particular isoform activates multiple cancer-related transcription factor reporters, while depletion of ZAK in gastric cell lines inhibits proliferation. These results reveal the spectrum of genomic and transcriptomic alterations in gastric cancer, and identify isoform-specific oncogenic properties of ZAK.

Diverse modes of genomic alteration in hepatocellular carcinoma

BackgroundHepatocellular carcinoma (HCC) is a heterogeneous disease with high mortality rate. Recent genomic studies have identified TP53, AXIN1, and CTNNB1 as the most frequently mutated genes. Lower frequency mutations have been reported in ARID1A, ARID2 and JAK1. In addition, hepatitis B virus (HBV) integrations into the human genome have been associated with HCC.ResultsHere, we deep-sequence 42 HCC patients with a combination of whole genome, exome and transcriptome sequencing to identify the mutational landscape of HCC using a reasonably large discovery cohort. We find frequent mutations in TP53, CTNNB1 and AXIN1, and rare but likely functional mutations in BAP1 and IDH1. Besides frequent hepatitis B virus integrations at TERT, we identify translocations at the boundaries of TERT. A novel deletion is identified in CTNNB1 in a region that is heavily mutated in multiple cancers. We also find multiple high-allelic frequency mutations in the extracellular matrix protein LAMA2. Lower expression levels of LAMA2 correlate with a proliferative signature, and predict poor survival and higher chance of cancer recurrence in HCC patients, suggesting an important role of the extracellular matrix and cell adhesion in tumor progression of a subgroup of HCC patients.ConclusionsThe heterogeneous disease of HCC features diverse modes of genomic alteration. In addition to common point mutations, structural variations and methylation changes, there are several virus-associated changes, including gene disruption or activation, formation of chimeric viral-human transcripts, and DNA copy number changes. Such a multitude of genomic events likely contributes to the heterogeneous nature of HCC.

MBASED: allele-specific expression detection in cancer tissues and cell lines

Allele-specific gene expression, ASE, is an important aspect of gene regulation. We developed a novel method MBASED, meta-analysis based allele-specific expression detection for ASE detection using RNA-seq data that aggregates information across multiple single nucleotide variation loci to obtain a gene-level measure of ASE, even when prior phasing information is unavailable. MBASED is capable of one-sample and two-sample analyses and performs well in simulations. We applied MBASED to a panel of cancer cell lines and paired tumor-normal tissue samples, and observed extensive ASE in cancer, but not normal, samples, mainly driven by genomic copy number alterations.