Immunogenomic pan-cancer landscape reveals immune escape mechanisms and immunoediting histories.

Immune reactions in the tumor microenvironment are an important hallmark of cancer, and emerging immune therapies have been proven effective against several types of cancers. To investigate cancer genome-immune interactions and the role of immunoediting or immune escape mechanisms in cancer development, we analyzed 2834 whole genome and RNA sequencing datasets across 31 distinct tumor types with respect to key immunogenomic aspects and provided comprehensive immunogenomic profiles of pan-cancers. We found that selective copy number changes in immune-related genes may contribute to immune escape. Furthermore, we developed an index of the immunoediting history of each tumor sample based on the information of mutations in exonic regions and pseudogenes and evaluated the immunoediting history of each tumor. Our immuno-genomic analyses of pan-cancers have the potential to identify a subset of tumors with immunogenicity and diverse backgrounds or intrinsic pathways associated with their immune status and immunoediting history.

Development of an MSI-positive colon tumor with aberrant DNA methylation in a PPAP patient.

Polymerase proofreading-associated polyposis (PPAP) is a disease caused by germline variations in the POLE and POLD1 genes that encode catalytic subunits of DNA polymerases. Studies of cancer genomes have identified somatic mutations in these genes, suggesting the importance of polymerase proofreading of DNA replication in suppressing tumorigenesis. Here, we identified a germline frameshift variation in the POLE gene (c.4191_4192delCT, p.Tyr1398*) in a case with multiple adenomatous polyps and three synchronous colon cancers. Interestingly, one of the colon cancers showed microsatellite instability-high (MSI-H) and another microsatellite stable. Immunohistochemical staining revealed that the MSI-H tumor cells lost the expression of MLH1 protein. Whole genome sequencing of the MSI-H tumor did not find pathogenic somatic mutations in mismatch repair genes but found frameshift mutations in the TET genes that catalyze 5-methylcytosine hydroxylation. Bisulfite sequencing of the tumor corroborated an increase in the number of hypermethylated regions including the MLH1 promoter. These data indicate that PPAP patients might develop MSI-positive tumors through epigenetic silencing of MLH1. These findings will contribute to comprehensive understanding of the molecular basis of tumors that involve deficiency of proofreading activity of DNA polymerases.

ALPHLARD-NT: Bayesian Method for Human Leukocyte Antigen Genotyping and Mutation Calling through Simultaneous Analysis of Normal and Tumor Whole-Genome Sequence Data.

Human leukocyte antigen (HLA) genes provide useful information on the relationship between cancer and the immune system. Despite the ease of obtaining these data through next-generation sequencing methods, interpretation of these relationships remains challenging owing to the complexity of HLA genes. To resolve this issue, we developed a Bayesian method, ALPHLARD-NT, to identify HLA germline and somatic mutations as well as HLA genotypes from whole-exome sequencing (WES) and whole-genome sequencing (WGS) data. ALPHLARD-NT showed 99.2% accuracy for WGS-based HLA genotyping and detected five HLA somatic mutations in 25 colon cancer cases. In addition, ALPHLARD-NT identified 88 HLA somatic mutations, including recurrent mutations and a novel HLA-B type, from WES data of 343 colon adenocarcinoma cases. These results demonstrate the potential of ALPHLARD-NT for conducting an accurate analysis of HLA genes even from low-coverage data sets. This method can become an essential tool for comprehensive analyses of HLA genes from WES and WGS data, helping to advance understanding of immune regulation in cancer as well as providing guidance for novel immunotherapy strategies.

Antigen delivery targeted to tumor-associated macrophages overcomes tumor immune resistance.

Immune checkpoint inhibitors and adoptive transfer of gene-engineered T cells have emerged as novel therapeutic modalities for hard-to-treat solid tumors; however, many patients are refractory to these immunotherapies, and the mechanisms underlying tumor immune resistance have not been fully elucidated. By comparing the tumor microenvironment of checkpoint inhibition-sensitive and -resistant murine solid tumors, we observed that the resistant tumors had low immunogenicity. We identified antigen presentation by CD11b+F4/80+ tumor-associated macrophages (TAMs) as a key factor correlated with immune resistance. In the resistant tumors, TAMs remained inactive and did not exert antigen-presenting activity. Targeted delivery of a long peptide antigen to TAMs by using a nano-sized hydrogel (nanogel) in the presence of a TLR agonist activated TAMs, induced their antigen-presenting activity, and thereby transformed the resistant tumors into tumors sensitive to adaptive immune responses such as adoptive transfer of tumor-specific T cell receptor-engineered T cells. These results indicate that the status and function of TAMs have a significant impact on tumor immune sensitivity and that manipulation of TAM functions would be an effective approach for improving the efficacy of immunotherapies.

ALPHLARD: a Bayesian method for analyzing HLA genes from whole genome sequence data.

BACKGROUND: Although human leukocyte antigen (HLA) genotyping based on amplicon, whole exome sequence (WES), and RNA sequence data has been achieved in recent years, accurate genotyping from whole genome sequence (WGS) data remains a challenge due to the low depth. Furthermore, there is no method to identify the sequences of unknown HLA types not registered in HLA databases. RESULTS: We developed a Bayesian model, called ALPHLARD, that collects reads potentially generated from HLA genes and accurately determines a pair of HLA types for each of HLA-A, -B, -C, -DPA1, -DPB1, -DQA1, -DQB1, and -DRB1 genes at 3rd field resolution. Furthermore, ALPHLARD can detect rare germline variants not stored in HLA databases and call somatic mutations from paired normal and tumor sequence data. We illustrate the capability of ALPHLARD using 253 WES data and 25 WGS data from Illumina platforms. By comparing the results of HLA genotyping from SBT and amplicon sequencing methods, ALPHLARD achieved 98.8% for WES data and 98.5% for WGS data at 2nd field resolution. We also detected three somatic point mutations and one case of loss of heterozygosity in the HLA genes from the WGS data. CONCLUSIONS: ALPHLARD showed good performance for HLA genotyping even from low-coverage data. It also has a potential to detect rare germline variants and somatic mutations in HLA genes. It would help to fill in the current gaps in HLA reference databases and unveil the immunological significance of somatic mutations identified in HLA genes.

Identification of an immunogenic neo-epitope encoded by mouse sarcoma using CXCR3 ligand mRNAs as sensors.

The CXCR3 ligands CXCL9, 10, and 11 play critical roles in the amplification of immune responses by recruiting CXCR3+ immune effector cells to the tumor site. Taking advantage of this property of CXCR3 ligands, we aimed to establish a novel approach to identify immunogenic mutated-antigens. We examined the feasibility of using CXCR3 ligand mRNAs as sensors for detection of specific immune responses in human and murine systems. We further investigated whether this approach is applicable for the identification of immunogenic mutated-antigens by using murine sarcoma lines. Rapid synthesis of CXCR3 ligand mRNAs occurred shortly after specific immune responses in both human and murine immune systems. Particularly, in CMS5 tumor-bearing mice, we detected specific immune responses to mutated mitogen-activated protein kinase 2 (ERK2), which has previously been identified as an immunogenic mutated-antigen. Furthermore, by combining this approach with whole-exome and transcriptome sequencing analyses, we identified an immunogenic neo-epitope derived from mutated staphylococcal nuclease domain-containing protein 1 (Snd1) in CMS7 tumor-bearing mice. Most importantly, we successfully detected the specific immune response to this neo-epitope even without co-administration of anti-cytotoxic T-lymphocyte protein-4 (CTLA-4), anti-programmed cell death-1 (PD-1) and anti-glucocorticoid-induced TNFR-related protein (GITR) antibodies, which vigorously augmented the immune response and consequently enabled us to detect the specific immune response to this neo-epitope by conventional IFNγ intracellular staining method. Our data indicate the potential usefulness of this strategy for the identification of immunogenic mutated-antigens. We propose that this approach would be of great help for the development of personalized cancer vaccine therapies in future.

Reduced expression of APC-1B but not APC-1A by the deletion of promoter 1B is responsible for familial adenomatous polyposis.

Germline mutations in the tumor suppressor gene APC are associated with familial adenomatous polyposis (FAP). Here we applied whole-genome sequencing (WGS) to the DNA of a sporadic FAP patient in which we did not find any pathological APC mutations by direct sequencing. WGS identified a promoter deletion of approximately 10 kb encompassing promoter 1B and exon1B of APC. Additional allele-specific expression analysis by deep cDNA sequencing revealed that the deletion reduced the expression of the mutated APC allele to as low as 11.2% in the total APC transcripts, suggesting that the residual mutant transcripts were driven by other promoter(s). Furthermore, cap analysis of gene expression (CAGE) demonstrated that the deleted promoter 1B region is responsible for the great majority of APC transcription in many tissues except the brain. The deletion decreased the transcripts of APC-1B to 39-45% in the patient compared to the healthy controls, but it did not decrease those of APC-1A. Different deletions including promoter 1B have been reported in FAP patients. Taken together, our results strengthen the evidence that analysis of structural variations in promoter 1B should be considered for the FAP patients whose pathological mutations are not identified by conventional direct sequencing.

Detection of APC mosaicism by next-generation sequencing in an FAP patient.

Familial adenomatous polyposis (FAP) of the colon is characterized by multiple polyps in the intestine and extra-colonic manifestations. Most FAP cases are caused by a germline mutation in the tumor-suppressor gene APC, but some cases of adenomatous polyposis result from germline mutations in MUTYH, POLD1 or POLE. Although sequence analysis of APC by the Sanger method is routinely performed for genetic testing, there remain cases whose mutations are not detected by the analysis. Next-generation sequencing has enabled us to analyze the comprehensive human genome, improving the chance of identifying disease causative variants. In this study, we conducted whole-genome sequencing of a sporadic FAP patient in which we did not find any pathogenic APC mutations by the conventional Sanger sequencing. Whole-genome sequencing and subsequent deep sequencing identified a mosaic mutation of c.3175G>T, p.E1059X in ~12% of his peripheral leukocytes. Additional deep sequencing of his buccal mucosa, hair follicles, non-cancerous mucosa of the stomach and colon disclosed that these tissues harbored the APC mutation at different frequencies. Our data implied that genetic analysis by next-generation sequencing is an effective strategy to identify genetic mosaicism in hereditary diseases.

Attenuated familial adenomatous polyposis with desmoids caused by an APC mutation.

We present here a case of attenuated familial adenomatous polyposis (AFAP) with a family history of desmoids and thyroid tumors. This patient had no colonic polyps but did have multiple desmoids. Genetic analysis identified a 4-bp deletion in codon 2644 (c.7932_7935delTTAT: p.Tyr2645LysfsX14) of the adenomatous polyposis coli (APC) gene. In cases with limited numbers of colonic polyps and desmoids, AFAP may be caused by a mutation in the 3' region of APC.