Motivation: The human leukocyte antigen (HLA) gene cluster plays a crucial role in adaptive immunity and is thus relevant in many biomedical applications. While next-generation sequencing data are often available for a patient, deducing the HLA genotype is difficult because of substantial sequence similarity within the cluster and exceptionally high variability of the loci. Established approaches, therefore, rely on specific HLA enrichment and sequencing techniques, coming at an additional cost and extra turnaround time.Result: We present OptiType, a novel HLA genotyping algorithm based on integer linear programming, capable of producing accurate predictions from NGS data not specifically enriched for the HLA cluster. We also present a comprehensive benchmark dataset consisting of RNA, exome and whole-genome sequencing data. OptiType significantly outperformed previously published in silico approaches with an overall accuracy of 97% enabling its use in a broad range of applications.Contact: szolek@informatik.uni-tuebingen.deSupplementary information: Supplementary data are available at Bioinformatics online.

Immune checkpoint inhibitors (ICIs) belong to the therapeutic armamentarium in advanced hepatocellular carcinoma (HCC). However, only a minority of patients benefit from immunotherapy. Therefore, we aimed to identify indicators of therapy response. This multicenter analysis included 99 HCC patients. Progression-free (PFS) and overall survival (OS) were studied by Kaplan-Meier analyses for clinical parameters using weighted log-rank testing. Next-generation sequencing (NGS) was performed in a subset of 15 patients. The objective response (OR) rate was 19% median OS (mOS)16.7 months. Forty-one percent reached a PFS > 6 months; these patients had a significantly longer mOS (32.0 vs. 8.5 months). Child-Pugh (CP) A and B patients showed a mOS of 22.1 and 12.1 months, respectively. Ten of thirty CP-B patients reached PFS > 6 months, including 3 patients with an OR. Tumor mutational burden (TMB) could not predict responders. Of note, antibiotic treatment within 30 days around ICI initiation was associated with significantly shorter mOS (8.5 vs. 17.4 months). Taken together, this study shows favorable outcomes for OS with low AFP, OR, and PFS > 6 months. No specific genetic pattern, including TMB, could identify responders. Antibiotics around treatment initiation were associated with worse outcome, suggesting an influence of the host microbiome on therapy success.

Immune checkpoint inhibitors (ICI) have revolutionized the therapeutic landscape of metastatic melanoma. However, ICI are often associated with immune-related adverse events (IRAE) such as colitis, hepatitis, pancreatitis, hypophysitis, pneumonitis, thyroiditis, exanthema, nephritis, myositis, encephalitis, or myocarditis. Biomarkers associated with the occurrence of IRAE would be desirable. In the literature, there is only little data available and furthermore mostly speculative, especially in view of genetic alterations. Our major aim was to check for possible associations between NGS-based genetic alterations and IRAE. We therefore analyzed 95 melanoma patients with ICI and evaluated their NGS results. We checked the data in view of potential associations between copy number variations (CNVs), small variations (VARs), human leucocyte antigen (HLA), sex, blood count parameters, pre-existing autoimmune diseases and the occurrence of IRAE. We conducted a literature research on genetic alterations hypothesized to be associated with the occurrence of IRAE. In total, we identified 39 genes that have been discussed as hypothetical biomarkers. We compared the list of these 39 genes with the tumor panel that our patients had received and focused our study on those 16 genes that were also included in the tumor panel used for NGS. Therefore, we focused our analyses on the following genes: AIRE, TERT, SH2B3, LRRK2, IKZF1, SMAD3, JAK2, PRDM1, CTLA4, TSHR, FAN1, SLCO1B1, PDCD1, IL1RN, CD274, UNG. We obtained relevant results: female sex was significantly associated with the development of hepatitis, combined immunotherapy with colitis, increased total and relative monocytes at therapy initiation were significantly associated with the development of pancreatitis, the same, pre-existing autoimmune diseases. Further significant associations were as follows: HLA homozygosity (hepatitis), and VARs on SMAD3 (pancreatitis). Regarding CNVs, significant markers included PRDM1 deletions and IL1RN (IRAE), CD274 duplications and SLCO1B1 (hepatitis), PRDM1 and CD274 (encephalitis), and PRDM1, CD274, TSHR, and FAN1 (myositis). Myositis and encephalitis, both, were associated with alterations of PRDM1 and CD274, which might explain their joined appearance in clinical practice. The association between HLA homozygosity and IRAE was clarified by finding HLA-A homozygosity as determining factor. We identified several genetic alterations hypothesized in the literature to be associated with the development of IRAE and found significant results concerning pre-existing autoimmune diseases and specific blood count parameters. Our findings can help to better understand the development of IRAE in melanoma patients. NGS might be a useful screening tool, however, our findings have yet to be confirmed in larger studies.

No personalized therapy regimens could demonstrate a benefit in survival of intrahepatic cholangiocarcinoma (iCCA). Since genetic heterogeneity might influence single biopsy based targeted therapy or the outcome of clinical trials, aim of the present study was to investigate intratumoral heterogeneity of iCCA by whole exome sequencing. Therefore, samples from tumor center and tumor periphery of large iCCA lesions as well as a control from healthy liver tissue were obtained from four patients and whole exome sequencing was performed. Mutations that occurred only in the tumor center or periphery were defined as private, whereas mutations present in both samples were regarded as common. A mean of 3 non-synonymous private mutations (range 0–14) per sample compared to 33,3 common mutations per sample (range 24–41) was identified. Mean percentage of non-synonymous private mutations per sample was 12% (range 0–58). In all samples of patient 1-3 as well as the central sample of patient 4 ≤ 10% private mutations were found, whereas 58% of private mutations were identified in the peripheral sample of patient 4. In this sample a private mutation in the DNA mismatch repair protein MSH6 could be identified most likely causing the high amount of private mutations. No substantial intratumoral heterogeneity was found in copy number variation analysis. In conclusion, iCCA show a small but distinct intratumoral heterogeneity. Somatic mutations in mismatch repair proteins might contribute significantly to increased spatial tumor burden and thereby may influence clinical management.