TF-FVIIa PAR2-ß-Arrestin Signaling Sustains Organ Dysfunction in Coxsackievirus B3 Infection of Mice.

BACKGROUND: Accumulating evidence implicates the activation of G-protein-coupled PARs (protease-activated receptors) by coagulation proteases in the regulation of innate immune responses. METHODS: Using mouse models with genetic alterations of the PAR2 signaling platform, we have explored contributions of PAR2 signaling to infection with coxsackievirus B3, a single-stranded RNA virus provoking multiorgan tissue damage, including the heart. RESULTS: We show that PAR2 activation sustains correlates of severe morbidity-hemodynamic compromise, aggravated hypothermia, and hypoglycemia-despite intact control of the virus. Following acute viral liver injury, canonical PAR2 signaling impairs the restoration process associated with exaggerated type I IFN (interferon) signatures in response to viral RNA recognition. Metabolic profiling in combination with proteomics of liver tissue shows PAR2-dependent reprogramming of liver metabolism, increased lipid droplet storage, and gluconeogenesis. PAR2-sustained hypodynamic compromise, reprograming of liver metabolism, as well as imbalanced IFN responses are prevented in ß-arrestin coupling-deficient PAR2 C-terminal phosphorylation mutant mice. Thus, wiring between upstream proteases and immune-metabolic responses results from biased PAR2 signaling mediated by intracellular recruitment of ß-arrestin. Importantly, blockade of the TF (tissue factor)-FVIIa (coagulation factor VIIa) complex capable of PAR2 proteolysis with the NAPc2 (nematode anticoagulant protein c2) mitigated virus-triggered pathology, recapitulating effects seen in protease cleavage-resistant PAR2 mice. CONCLUSIONS: These data provide insights into a TF-FVIIa signaling axis through PAR2-ß-arrestin coupling that is a regulator of inflammation-triggered tissue repair and hemodynamic compromise in coxsackievirus B3 infection and can potentially be targeted with selective coagulation inhibitors.

Transcriptional reprogramming by mutated IRF4 in lymphoma.

Disease-causing mutations in genes encoding transcription factors (TFs) can affect TF interactions with their cognate DNA-binding motifs. Whether and how TF mutations impact upon the binding to TF composite elements (CE) and the interaction with other TFs is unclear. Here, we report a distinct mechanism of TF alteration in human lymphomas with perturbed B cell identity, in particular classic Hodgkin lymphoma. It is caused by a recurrent somatic missense mutation c.295 T > C (p.Cys99Arg; p.C99R) targeting the center of the DNA-binding domain of Interferon Regulatory Factor 4 (IRF4), a key TF in immune cells. IRF4-C99R fundamentally alters IRF4 DNA-binding, with loss-of-binding to canonical IRF motifs and neomorphic gain-of-binding to canonical and non-canonical IRF CEs. IRF4-C99R thoroughly modifies IRF4 function by blocking IRF4-dependent plasma cell induction, and up-regulates disease-specific genes in a non-canonical Activator Protein-1 (AP-1)-IRF-CE (AICE)-dependent manner. Our data explain how a single mutation causes a complex switch of TF specificity and gene regulation and open the perspective to specifically block the neomorphic DNA-binding activities of a mutant TF.

CoVigator-A Knowledge Base for Navigating SARS-CoV-2 Genomic Variants.

BACKGROUND: The outbreak of the severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2) resulted in the global COVID-19 pandemic. The urgency for an effective SARS-CoV-2 vaccine has led to the development of the first series of vaccines at unprecedented speed. The discovery of SARS-CoV-2 spike-glycoprotein mutants, however, and consequentially the potential to escape vaccine-induced protection and increased infectivity, demonstrates the persisting importance of monitoring SARS-CoV-2 mutations to enable early detection and tracking of genomic variants of concern. RESULTS: We developed the CoVigator tool with three components: (1) a knowledge base that collects new SARS-CoV-2 genomic data, processes it and stores its results; (2) a comprehensive variant calling pipeline; (3) an interactive dashboard highlighting the most relevant findings. The knowledge base routinely downloads and processes virus genome assemblies or raw sequencing data from the COVID-19 Data Portal (C19DP) and the European Nucleotide Archive (ENA), respectively. The results of variant calling are visualized through the dashboard in the form of tables and customizable graphs, making it a versatile tool for tracking SARS-CoV-2 variants. We put a special emphasis on the identification of intrahost mutations and make available to the community what is, to the best of our knowledge, the largest dataset on SARS-CoV-2 intrahost mutations. In the spirit of open data, all CoVigator results are available for download. The CoVigator dashboard is accessible via covigator.tron-mainz.de. CONCLUSIONS: With increasing demand worldwide in genome surveillance for tracking the spread of SARS-CoV-2, CoVigator will be a valuable resource of an up-to-date list of mutations, which can be incorporated into global efforts.

Characterization of BV6-Induced Sensitization to the NK Cell Killing of Pediatric Rhabdomyosarcoma Spheroids.

Although the overall survival in pediatric rhabdomyosarcoma (RMS) has increased over the last decades, the most aggressive subtype of alveolar RMS is in dire need of novel treatment strategies. RMS cells evade cell death induction and immune control by increasing the expression of inhibitors of apoptosis proteins (IAPs), which can be exploited and targeted with stimulation with Smac mimetics. Here, we used the Smac mimetic BV6 to re-sensitize RMS spheroids to cell death, which increased killing induced by natural killer (NK) cells. Single BV6 treatment of RMS spheroids did not reduce spheroidal growth. However, we observed significant spheroidal decomposition upon BV6 pre-treatment combined with NK cell co-cultivation. Molecularly, IAPs s are rapidly degraded by BV6, which activates NF-κB signal transduction pathways in RMS spheroids. RNA sequencing analysis validated NF-κB activation and identified a plethora of BV6-regulated genes. Additionally, BV6 released caspases from IAP-mediated inhibition. Here, caspase-8 might play a major role, as knockdown experiments resulted in decreased NK cell-mediated attack. Taken together, we improved the understanding of the BV6 mechanism of RMS spheroid sensitization to cytotoxic immune cells, which could be suitable for the development of novel combinatory cellular immunotherapy with Smac mimetics.

Large-scale analysis of SARS-CoV-2 spike-glycoprotein mutants demonstrates the need for continuous screening of virus isolates.

Due to the widespread of the COVID-19 pandemic, the SARS-CoV-2 genome is evolving in diverse human populations. Several studies already reported different strains and an increase in the mutation rate. Particularly, mutations in SARS-CoV-2 spike-glycoprotein are of great interest as it mediates infection in human and recently approved mRNA vaccines are designed to induce immune responses against it. We analyzed 1,036,030 SARS-CoV-2 genome assemblies and 30,806 NGS datasets from GISAID and European Nucleotide Archive (ENA) focusing on non-synonymous mutations in the spike protein. Only around 2.5% of the samples contained the wild-type spike protein with no variation from the reference. Among the spike protein mutants, we confirmed a low mutation rate exhibiting less than 10 non-synonymous mutations in 99.6% of the analyzed sequences, but the mean and median number of spike protein mutations per sample increased over time. 5,472 distinct variants were found in total. The majority of the observed variants were recurrent, but only 21 and 14 recurrent variants were found in at least 1% of the mutant genome assemblies and NGS samples, respectively. Further, we found high-confidence subclonal variants in about 2.6% of the NGS data sets with mutant spike protein, which might indicate co-infection with various SARS-CoV-2 strains and/or intra-host evolution. Lastly, some variants might have an effect on antibody binding or T-cell recognition. These findings demonstrate the continuous importance of monitoring SARS-CoV-2 sequences for an early detection of variants that require adaptations in preventive and therapeutic strategies.

A noninflammatory mRNA vaccine for treatment of experimental autoimmune encephalomyelitis.

The ability to control autoreactive T cells without inducing systemic immune suppression is the major goal for treatment of autoimmune diseases. The key challenge is the safe and efficient delivery of pharmaceutically well-defined antigens in a noninflammatory context. Here, we show that systemic delivery of nanoparticle-formulated 1 methylpseudouridine-modified messenger RNA (m1Ψ mRNA) coding for disease-related autoantigens results in antigen presentation on splenic CD11c+ antigen-presenting cells in the absence of costimulatory signals. In several mouse models of multiple sclerosis, the disease is suppressed by treatment with such m1Ψ mRNA. The treatment effect is associated with a reduction of effector T cells and the development of regulatory T cell (Treg cell) populations. Notably, these Treg cells execute strong bystander immunosuppression and thus improve disease induced by cognate and noncognate autoantigens.

Positive Role of the MHC Class-I Antigen Presentation Regulator m04/gp34 of Murine Cytomegalovirus in Antiviral Protection by CD8 T Cells.

Murine cytomegalovirus (mCMV) codes for MHC class-I trafficking modulators m04/gp34, m06/gp48, and m152/gp40. By interacting with the MHC class-Iα chain, these proteins disconnect peptide-loaded MHC class-I (pMHC-I) complexes from the constitutive vesicular flow to the cell surface. Based on the assumption that all three inhibit antigen presentation, and thus the recognition of infected cells by CD8 T cells, they were referred to as "immunoevasins." Improved antigen presentation mediated by m04 in the presence of m152 after infection with deletion mutant mCMV-Δm06W, compared to mCMV-Δm04m06 expressing only m152, led us to propose renaming these molecules "viral regulators of antigen presentation" (vRAP) to account for both negative and positive functions. In accordance with a positive function, m04-pMHC-I complexes were found to be displayed on the cell surface, where they are primarily known as ligands for Ly49 family natural killer (NK) cell receptors. Besides the established role of m04 in NK cell silencing or activation, an anti-immunoevasive function by activation of CD8 T cells is conceivable, because the binding site of m04 to MHC class-Iα appears not to mask the peptide binding site for T-cell receptor recognition. However, functional evidence was based on mCMV-Δm06W, a virus of recently doubted authenticity. Here we show that mCMV-Δm06W actually represents a mixture of an authentic m06 deletion mutant and a mutant with an accidental additional deletion of a genome region encompassing also gene m152. Reanalysis of previously published experiments for the authentic mutant in the mixture confirms the previously concluded positive vRAP function of m04.

A liposomal RNA vaccine inducing neoantigen-specific CD4+ T cells augments the antitumor activity of local radiotherapy in mice.

Antigen-encoding, lipoplex-formulated RNA (RNA-LPX) enables systemic delivery to lymphoid compartments and selective expression in resident antigen-presenting cells. We report here that the rejection of CT26 tumors, mediated by local radiotherapy (LRT), is further augmented in a CD8+ T cell-dependent manner by an RNA-LPX vaccine that encodes CD4+ T cell-recognized neoantigens (CD4 neoantigen vaccine). Whereas CD8+ T cells induced by LRT alone were primarily directed against the immunodominant gp70 antigen, mice treated with LRT plus the CD4 neoantigen vaccine rejected gp70-negative tumors and were protected from rechallenge with these tumors, indicating a potent poly-antigenic CD8+ T cell response and T cell memory. In the spleens of CD4 neoantigen-vaccinated mice, we found a high number of activated, poly-functional, Th1-like CD4+ T cells against ME1, the immunodominant CD4 neoantigen within the poly-neoantigen vaccine. LRT itself strongly increased CD8+ T cell numbers and clonal expansion. However, tumor infiltrates of mice treated with CD4 neoantigen vaccine/LRT, as compared to LRT alone, displayed a higher fraction of activated gp70-specific CD8+ T cells, lower PD-1/LAG-3 expression and contained ME1-specific IFNγ+ CD4+ T cells capable of providing cognate help. CD4 neoantigen vaccine/LRT treatment followed by anti-CTLA-4 antibody therapy further enhanced the efficacy with complete remission of gp70-negative CT26 tumors and survival of all mice. Our data highlight the power of combining synergistic modes of action and warrants further exploration of the presented treatment schema.

Multi-Omics Characterization of the 4T1 Murine Mammary Gland Tumor Model.

Background: Tumor models are critical for our understanding of cancer and the development of cancer therapeutics. The 4T1 murine mammary cancer cell line is one of the most widely used breast cancer models. Here, we present an integrated map of the genome, transcriptome, and immunome of 4T1. Results: We found Trp53 (Tp53) and Pik3g to be mutated. Other frequently mutated genes in breast cancer, including Brca1 and Brca2, are not mutated. For cancer related genes, Nav3, Cenpf, Muc5Ac, Mpp7, Gas1, MageD2, Dusp1, Ros, Polr2a, Rragd, Ros1, and Hoxa9 are mutated. Markers for cell proliferation like Top2a, Birc5, and Mki67 are highly expressed, so are markers for metastasis like Msln, Ect2, and Plk1, which are known to be overexpressed in triple-negative breast cancer (TNBC). TNBC markers are, compared to a mammary gland control sample, lower (Esr1), comparably low (Erbb2), or not expressed at all (Pgr). We also found testis cancer antigen Pbk as well as colon/gastrointestinal cancer antigens Gpa33 and Epcam to be highly expressed. Major histocompatibility complex (MHC) class I is expressed, while MHC class II is not. We identified 505 single nucleotide variations (SNVs) and 20 insertions and deletions (indels). Neoantigens derived from 22 SNVs and one deletion elicited CD8+ or CD4+ T cell responses in IFNγ-ELISpot assays. Twelve high-confidence fusion genes were observed. We did not observe significant downregulation of mismatch repair (MMR) genes or SNVs/indels impairing their function, providing evidence for 6-thioguanine resistance. Effects of the integration of the murine mammary tumor virus were observed at the genome and transcriptome level. Conclusions: 4T1 cells share substantial molecular features with human TNBC. As 4T1 is a common model for metastatic tumors, our data supports the rational design of mode-of-action studies for pre-clinical evaluation of targeted immunotherapies.

Bioinformatics for Cancer Immunotherapy.

Our immune system plays a key role in health and disease as it is capable of responding to foreign antigens as well as acquired antigens from cancer cells. Latter are caused by somatic mutations, the so-called neoepitopes, and might be recognized by T cells if they are presented by HLA molecules on the surface of cancer cells. Personalized mutanome vaccines are a class of customized immunotherapies, which is dependent on the detection of individual cancer-specific tumor mutations and neoepitope (i.e., prediction, followed by a rational vaccine design, before on-demand production. The development of next generation sequencing (NGS) technologies and bioinformatic tools allows a large-scale analysis of each parameter involved in this process. Here, we provide an overview of the bioinformatic aspects involved in the design of personalized, neoantigen-based vaccines, including the detection of mutations and the subsequent prediction of potential epitopes, as well as methods for associated biomarker research, such as high-throughput sequencing of T-cell receptors (TCRs), followed by data analysis and the bioinformatics quantification of immune cell infiltration in cancer samples.

An RNA vaccine drives expansion and efficacy of claudin-CAR-T cells against solid tumors.

Chimeric antigen receptor (CAR)-T cells have shown efficacy in patients with B cell malignancies. Yet, their application for solid tumors has challenges that include limited cancer-specific targets and nonpersistence of adoptively transferred CAR-T cells. Here, we introduce the developmentally regulated tight junction protein claudin 6 (CLDN6) as a CAR target in solid tumors and a strategy to overcome inefficient CAR-T cell stimulation in vivo. We demonstrate that a nanoparticulate RNA vaccine, designed for body-wide delivery of the CAR antigen into lymphoid compartments, stimulates adoptively transferred CAR-T cells. Presentation of the natively folded target on resident antigen-presenting cells promotes cognate and selective expansion of CAR-T cells. Improved engraftment of CAR-T cells and regression of large tumors in difficult-to-treat mouse models was achieved at subtherapeutic CAR-T cell doses.

HPV16 RNA-LPX vaccine mediates complete regression of aggressively growing HPV-positive mouse tumors and establishes protective T cell memory.

HPV16 infections are associated with a variety of cancers and there is compelling evidence that the transforming activity of HPV16 critically depends on the expression of the viral oncoproteins E6 and E7. Therapeutic cancer vaccines capable of generating durable and specific immunity against these HPV16 antigens hold great promise to achieve long-term disease control. Here we show in mice that HPV16 E7 RNA-LPX, an intravenously administered cancer vaccine based on immuno-pharmacologically optimized antigen-encoding mRNA, efficiently primes and expands antigen-specific effector and memory CD8+ T cells. HPV-positive TC-1 and C3 tumors of immunized mice are heavily infiltrated with activated immune cells and HPV16-specific T cells and are polarized towards a proinflammatory, cytotoxic and less immune-suppressive contexture. E7 RNA-LPX immunization mediates complete and durable remission of progressing tumors. Circulating memory T cells are highly cytotoxic and protect from tumor rechallenge. Moreover, E7 RNA-LPX immunization sensitizes anti-PD-L1 refractory tumors to checkpoint blockade. In conclusion, our data highlight the potential of HPV16 RNA-LPX for the treatment of HPV-driven cancers.

In Silico Typing of Classical and Non-classical HLA Alleles from Standard RNA-Seq Reads.

Next-Generation Sequencing (NGS) enables the rapid generation of billions of short nucleic acid sequence fragments (i.e., "sequencing reads"). Especially, the adoption of gene expression profiling using whole transcriptome sequencing (i.e., "RNA-Seq") has been rapid. Here, we describe an in silico method, seq2HLA, that takes standard RNA-Seq reads as input and determines a sample's (classical and non-classical) HLA class I and class II types as well as HLA expression. We demonstrate the application of seq2HLA using publicly available RNA-Seq data from the Burkitt's lymphoma cell line DAUDI and the choriocarcinoma cell line JEG-3.

HLA and proteasome expression body map.

BACKGROUND: The presentation of HLA peptide complexes to T cells is a highly regulated and tissue specific process involving multiple transcriptionally controlled cellular components. The extensive polymorphism of HLA genes and the complex composition of the proteasome make it difficult to map their expression profiles across tissues. METHODS: Here we applied a tailored gene quantification pipeline to 4323 publicly available RNA-Seq datasets representing 55 normal tissues and cell types to examine expression profiles of (classical and non-classical) HLA class I, class II and proteasomal genes. RESULTS: We generated the first comprehensive expression atlas of antigen presenting-related genes across 56 normal tissues and cell types, including immune cells, pancreatic islets, platelets and hematopoietic stem cells. We found a surprisingly heterogeneous HLA expression pattern with up to 100-fold difference in intra-tissue median HLA abundances. Cells of the immune system and lymphatic organs expressed the highest levels of classical HLA class I (HLA-A,-B,-C), class II (HLA-DQA1,-DQB1,-DPA1,-DPB1,-DRA,-DRB1) and non-classical HLA class I (HLA-E,-F) molecules, whereas retina, brain, muscle, megakaryocytes and erythroblasts showed the lowest abundance. In contrast, we identified a distinct and highly tissue-restricted expression pattern of the non-classical class I gene HLA-G in placenta, pancreatic islets, pituitary gland and testis. While the constitutive proteasome showed relatively constant expression across all tissues, we found the immunoproteasome to be enriched in lymphatic organs and almost absent in immune privileged tissues. CONCLUSIONS: Here, we not only provide a reference catalog of tissue and cell type specific HLA expression, but also highlight extremely variable expression of the basic components of antigen processing and presentation in different cell types. Our findings indicate that low expression of classical HLA class I molecules together with lack of immunoproteasome components as well as upregulation of HLA-G may be of key relevance to maintain tolerance in immune privileged tissues.

TCLP: an online cancer cell line catalogue integrating HLA type, predicted neo-epitopes, virus and gene expression.

Human cancer cell lines are an important resource for research and drug development. However, the available annotations of cell lines are sparse, incomplete, and distributed in multiple repositories. Re-analyzing publicly available raw RNA-Seq data, we determined the human leukocyte antigen (HLA) type and abundance, identified expressed viruses and calculated gene expression of 1,082 cancer cell lines. Using the determined HLA types, public databases of cell line mutations, and existing HLA binding prediction algorithms, we predicted antigenic mutations in each cell line. We integrated the results into a comprehensive knowledgebase. Using the Django web framework, we provide an interactive user interface with advanced search capabilities to find and explore cell lines and an application programming interface to extract cell line information. The portal is available at http://celllines.tron-mainz.de.

Immunomic, genomic and transcriptomic characterization of CT26 colorectal carcinoma.

BACKGROUND: Tumor models are critical for our understanding of cancer and the development of cancer therapeutics. Here, we present an integrated map of the genome, transcriptome and immunome of an epithelial mouse tumor, the CT26 colon carcinoma cell line. RESULTS: We found that Kras is homozygously mutated at p.G12D, Apc and Tp53 are not mutated, and Cdkn2a is homozygously deleted. Proliferation and stem-cell markers, including Top2a, Birc5 (Survivin), Cldn6 and Mki67, are highly expressed while differentiation and top-crypt markers Muc2, Ms4a8a (MS4A8B) and Epcam are not. Myc, Trp53 (tp53), Mdm2, Hif1a, and Nras are highly expressed while Egfr and Flt1 are not. MHC class I but not MHC class II is expressed. Several known cancer-testis antigens are expressed, including Atad2, Cep55, and Pbk. The highest expressed gene is a mutated form of the mouse tumor antigen gp70. Of the 1,688 non-synonymous point variations, 154 are both in expressed genes and in peptides predicted to bind MHC and thus potential targets for immunotherapy development. Based on its molecular signature, we predicted that CT26 is refractory to anti-EGFR mAbs and sensitive to MEK and MET inhibitors, as have been previously reported. CONCLUSIONS: CT26 cells share molecular features with aggressive, undifferentiated, refractory human colorectal carcinoma cells. As CT26 is one of the most extensively used syngeneic mouse tumor models, our data provide a map for the rationale design of mode-of-action studies for pre-clinical evaluation of targeted- and immunotherapies.

Cell contact-dependent priming and Fc interaction with CD32+ immune cells contribute to the TGN1412-triggered cytokine response.

Following inconspicuous preclinical testing, the superagonistic anti-CD28 mAb TGN1412 was applied to six study participants who all developed a devastating cytokine storm. We verified that TGN1412 treatment of fresh PBMCs induced only moderate responses, whereas restoration of tissue-like conditions by high-density preculture (HDC) allowed vigorous cytokine production. TGN1412 treatment of T cells isolated from HDC-PBMCs induced moderate cytokine responses, which upon additional anti-IgG crosslinking were significantly boosted. Moreover, coincubation of TGN1412-treated T cells with B cells expressing the intermediate affinity Fcγ receptor IIB (CD32B), or coincubation with CD32B(+) transfectants, resulted in robust T cell activation. This was surprising because TGN1412 was expressed as an Ig of the subclass 4 (IgG4), which was shown before to exhibit only minor affinity to FcγRs. Transcriptome analysis of TGN1412-treated T cells revealed that similar gene signatures were induced irrespective of whether T cells derived from fresh or HDC-PBMCs were studied. Collectively, these data indicate that HDC-PBMCs and HDC-PBMC-derived T cells mount rapid TGN1412 responses, which are massively boosted by FcγR crosslinking, in particular by CD32-expressing B cells. These results qualify HDC-PBMCs as a valuable in vitro test system for the analysis of complex mAb functions.

A catalog of HLA type, HLA expression, and neo-epitope candidates in human cancer cell lines.

Cancer cell lines are a tremendous resource for cancer biology and therapy development. These multipurpose tools are commonly used to examine the genetic origin of cancers, to identify potential novel tumor targets, such as tumor antigens for vaccine devel-opment, and utilized to screen potential therapies in preclinical studies. Mutations, gene expression, and drug sensitivity have been determined for many cell lines using next-generation sequencing (NGS). However, the human leukocyte antigen (HLA) type and HLA expression of tumor cell lines, characterizations necessary for the development of cancer vaccines, have remained largely incomplete and, such information, when available, has been distributed in many publications. Here, we determine the 4-digit HLA type and HLA expression of 167 cancer and 10 non-cancer cell lines from publically available RNA-Seq data. We use standard NGS RNA-Seq short reads from "whole transcriptome" sequencing, map reads to known HLA types, and statistically determine HLA type, heterozygosity, and expression. First, we present previously unreported HLA Class I and II genotypes. Second, we determine HLA expression levels in each cancer cell line, providing insights into HLA downregulation and loss in cancer. Third, using these results, we provide a fundamental cell line "barcode" to track samples and prevent sample annotation swaps and contamination. Fourth, we integrate the cancer cell-line specific HLA types and HLA expression with available cell-line specific mutation information and existing HLA binding prediction algorithms to make a catalog of predicted antigenic mutations in each cell line. The compilation of our results are a fundamental resource for all researchers selecting specific cancer cell lines based on the HLA type and HLA expression, as well as for the development of immunotherapeutic tools for novel cancer treatment modalities.

HLA typing from RNA-Seq sequence reads.

We present a method, seq2HLA, for obtaining an individual's human leukocyte antigen (HLA) class I and II type and expression using standard next generation sequencing RNA-Seq data. RNA-Seq reads are mapped against a reference database of HLA alleles, and HLA type, confidence score and locus-specific expression level are determined. We successfully applied seq2HLA to 50 individuals included in the HapMap project, yielding 100% specificity and 94% sensitivity at a P-value of 0.1 for two-digit HLA types. We determined HLA type and expression for previously un-typed Illumina Body Map tissues and a cohort of Korean patients with lung cancer. Because the algorithm uses standard RNA-Seq reads and requires no change to laboratory protocols, it can be used for both existing datasets and future studies, thus adding a new dimension for HLA typing and biomarker studies.