ZygosityPredictor.

Summary: ZygosityPredictor provides functionality to evaluate how many copies of a gene are affected by mutations in next generation sequencing data. In cancer samples, the tool processes both somatic and germline mutations. In particular, ZygosityPredictor computes the number of affected copies for single nucleotide variants and small insertions and deletions (Indels). In addition, the tool integrates information at gene level via phasing of several variants and subsequent logic to derive how strongly a gene is affected by mutations and provides a measure of confidence. This information is of particular interest in precision oncology, e.g. when assessing whether unmutated copies of tumor-suppressor genes remain. Availability and implementation: ZygosityPredictor was implemented as an R-package and is available via Bioconductor at https://bioconductor.org/packages/ZygosityPredictor. Detailed documentation is provided in the vignette including application to an example genome.

Ex vivo drug response profiling for response and outcome prediction in hematologic malignancies: the prospective non-interventional SMARTrial.

Ex vivo drug response profiling is a powerful tool to study genotype-drug response associations and is being explored as a tool set for precision medicine in cancer. Here we conducted a prospective non-interventional trial to investigate feasibility of ex vivo drug response profiling for treatment guidance in hematologic malignancies (SMARTrial, NCT03488641 ). The primary endpoint to provide drug response profiling reports within 7 d was met in 91% of all study participants (N = 80). Secondary endpoint analysis revealed that ex vivo resistance to chemotherapeutic drugs predicted chemotherapy treatment failure in vivo. We confirmed the predictive value of ex vivo response to chemotherapy in a validation cohort of 95 individuals with acute myeloid leukemia treated with daunorubicin and cytarabine. Ex vivo drug response profiles improved ELN-22 risk stratification in individuals with adverse risk. We conclude that ex vivo drug response profiling is clinically feasible and has the potential to predict chemotherapy response in individuals with hematologic malignancies beyond clinically established genetic markers.

A T-cell antigen atlas for meningioma: novel options for immunotherapy.

Meningiomas are the most common primary intracranial tumors. Although most symptomatic cases can be managed by surgery and/or radiotherapy, a relevant number of patients experience an unfavorable clinical course and additional treatment options are needed. As meningiomas are often perfused by dural branches of the external carotid artery, which is located outside the blood-brain barrier, they might be an accessible target for immunotherapy. However, the landscape of naturally presented tumor antigens in meningioma is unknown. We here provide a T-cell antigen atlas for meningioma by in-depth profiling of the naturally presented immunopeptidome using LC-MS/MS. Candidate target antigens were selected based on a comparative approach using an extensive immunopeptidome data set of normal tissues. Meningioma-exclusive antigens for HLA class I and II are described here for the first time. Top-ranking targets were further functionally characterized by showing their immunogenicity through in vitro T-cell priming assays. Thus, we provide an atlas of meningioma T-cell antigens which will be publicly available for further research. In addition, we have identified novel actionable targets that warrant further investigation as an immunotherapy option for meningioma.

PyOncoPrint: a python package for plotting OncoPrints.

OncoPrint, the plot to visualize an overview of genetic variants in sequencing data, has been widely used in the field of cancer genomics. However, still, there have been no Python libraries capable to generate OncoPrint yet, a big hassle to plot OncoPrints within Python-based genetic variants analysis pipelines. This paper introduces a new Python package PyOncoPrint, which can be easily used to plot OncoPrints in Python. The package is based on the existing widely used scientific plotting library Matplotlib, the resulting plots are easy to be adjusted for various needs.

Targeted therapy of advanced parathyroid carcinoma guided by genomic and transcriptomic profiling.

Parathyroid carcinoma (PC) is an ultra-rare malignancy with a high risk of recurrence after surgery. Tumour-directed systemic treatments for PC are not established. We used whole-genome and RNA sequencing in four patients with advanced PC to identify molecular alterations that could guide clinical management. In two cases, the genomic and transcriptomic profiles provided targets for experimental therapies that resulted in biochemical response and prolonged disease stabilization: (a) immune checkpoint inhibition with pembrolizumab based on high tumour mutational burden and a single-base substitution signature associated with APOBEC (apolipoprotein B mRNA editing enzyme, catalytic polypeptide-like) overactivation; (b) multi-receptor tyrosine kinase inhibition with lenvatinib due to overexpression of FGFR1 (Fibroblast Growth Factor Receptor 1) and RET (Ret Proto-Oncogene) and, (c) later in the course of the disease, PARP (Poly(ADP-Ribose) Polymerase) inhibition with olaparib prompted by signs of defective homologous recombination DNA repair. In addition, our data provided new insights into the molecular landscape of PC with respect to the genome-wide footprints of specific mutational processes and pathogenic germline alterations. These data underscore the potential of comprehensive molecular analyses to improve care for patients with ultra-rare cancers based on insight into disease biology.

Comprehensive genomic and epigenomic analysis in cancer of unknown primary guides molecularly-informed therapies despite heterogeneity.

The benefit of molecularly-informed therapies in cancer of unknown primary (CUP) is unclear. Here, we use comprehensive molecular characterization by whole genome/exome, transcriptome and methylome analysis in 70 CUP patients to reveal substantial mutational heterogeneity with TP53, MUC16, KRAS, LRP1B and CSMD3 being the most frequently mutated known cancer-related genes. The most common fusion partner is FGFR2, the most common focal homozygous deletion affects CDKN2A. 56/70 (80%) patients receive genomics-based treatment recommendations which are applied in 20/56 (36%) cases. Transcriptome and methylome data provide evidence for the underlying entity in 62/70 (89%) cases. Germline analysis reveals five (likely) pathogenic mutations in five patients. Recommended off-label therapies translate into a mean PFS ratio of 3.6 with a median PFS1 of 2.9 months (17 patients) and a median PFS2 of 7.8 months (20 patients). Our data emphasize the clinical value of molecular analysis and underline the need for innovative, mechanism-based clinical trials.

Recurrent Germline Variant in RAD21 Predisposes Children to Lymphoblastic Leukemia or Lymphoma.

Somatic loss of function mutations in cohesin genes are frequently associated with various cancer types, while cohesin disruption in the germline causes cohesinopathies such as Cornelia-de-Lange syndrome (CdLS). Here, we present the discovery of a recurrent heterozygous RAD21 germline aberration at amino acid position 298 (p.P298S/A) identified in three children with lymphoblastic leukemia or lymphoma in a total dataset of 482 pediatric cancer patients. While RAD21 p.P298S/A did not disrupt the formation of the cohesin complex, it altered RAD21 gene expression, DNA damage response and primary patient fibroblasts showed increased G2/M arrest after irradiation and Mitomycin-C treatment. Subsequent single-cell RNA-sequencing analysis of healthy human bone marrow confirmed the upregulation of distinct cohesin gene patterns during hematopoiesis, highlighting the importance of RAD21 expression within proliferating B- and T-cells. Our clinical and functional data therefore suggest that RAD21 germline variants can predispose to childhood lymphoblastic leukemia or lymphoma without displaying a CdLS phenotype.

Whole exome sequencing identifies novel germline variants of SLC15A4 gene as potentially cancer predisposing in familial colorectal cancer.

About 15% of colorectal cancer (CRC) patients have first-degree relatives affected by the same malignancy. However, for most families the cause of familial aggregation of CRC is unknown. To identify novel high-to-moderate-penetrance germline variants underlying CRC susceptibility, we performed whole exome sequencing (WES) on four CRC cases and two unaffected members of a Polish family without any mutation in known CRC predisposition genes. After WES, we used our in-house developed Familial Cancer Variant Prioritization Pipeline and identified two novel variants in the solute carrier family 15 member 4 (SLC15A4) gene. The heterozygous missense variant, p. Y444C, was predicted to affect the phylogenetically conserved PTR2/POT domain and to have a deleterious effect on the function of the encoded peptide/histidine transporter. The other variant was located in the upstream region of the same gene (GRCh37.p13, 12_129308531_C_T; 43 bp upstream of transcription start site, ENST00000266771.5) and it was annotated to affect the promoter region of SLC15A4 as well as binding sites of 17 different transcription factors. Our findings of two distinct variants in the same gene may indicate a synergistic up-regulation of SLC15A4 as the underlying genetic cause and implicate this gene for the first time in genetic inheritance of familial CRC.

The genomic and transcriptional landscape of primary central nervous system lymphoma.

Primary lymphomas of the central nervous system (PCNSL) are mainly diffuse large B-cell lymphomas (DLBCLs) confined to the central nervous system (CNS). Molecular drivers of PCNSL have not been fully elucidated. Here, we profile and compare the whole-genome and transcriptome landscape of 51 CNS lymphomas (CNSL) to 39 follicular lymphoma and 36 DLBCL cases outside the CNS. We find recurrent mutations in JAK-STAT, NFkB, and B-cell receptor signaling pathways, including hallmark mutations in MYD88 L265P (67%) and CD79B (63%), and CDKN2A deletions (83%). PCNSLs exhibit significantly more focal deletions of HLA-D (6p21) locus as a potential mechanism of immune evasion. Mutational signatures correlating with DNA replication and mitosis are significantly enriched in PCNSL. TERT gene expression is significantly higher in PCNSL compared to activated B-cell (ABC)-DLBCL. Transcriptome analysis clearly distinguishes PCNSL and systemic DLBCL into distinct molecular subtypes. Epstein-Barr virus (EBV)+ CNSL cases lack recurrent mutational hotspots apart from IG and HLA-DRB loci. We show that PCNSL can be clearly distinguished from DLBCL, having distinct expression profiles, IG expression and translocation patterns, as well as specific combinations of genetic alterations.

Germline Variants of CYBA and TRPM4 Predispose to Familial Colorectal Cancer.

Familial colorectal cancer (CRC) is only partially explained by known germline predisposing genes. We performed whole-genome sequencing in 15 Polish families of many affected individuals, without mutations in known CRC predisposing genes. We focused on loss-of-function variants and functionally characterized them. We identified a frameshift variant in the CYBA gene (c.246delC) in one family and a splice site variant in the TRPM4 gene (c.25-1 G > T) in another family. While both variants were absent or extremely rare in gene variant databases, we identified four additional Polish familial CRC cases and two healthy elderly individuals with the CYBA variant (odds ratio 2.46, 95% confidence interval 0.48-12.69). Both variants led to a premature stop codon and to a truncated protein. Functional characterization of the variants showed that knockdown of CYBA or TRPM4 depressed generation of reactive oxygen species (ROS) in LS174T and HT-29 cell lines. Knockdown of TRPM4 resulted in decreased MUC2 protein production. CYBA encodes a component in the NADPH oxidase system which generates ROS and controls, e.g., bacterial colonization in the gut. Germline CYBA variants are associated with early onset inflammatory bowel disease, supported with experimental evidence on loss of intestinal mucus barrier function due to ROS deficiency. TRPM4 encodes a calcium-activated ion channel, which, in a human colonic cancer cell line, controls calcium-mediated secretion of MUC2, a major component of intestinal mucus barrier. We suggest that the gene defects in CYBA and TRPM4 mechanistically involve intestinal barrier integrity through ROS and mucus biology, which converges in chronic bowel inflammation.

Whole-Exome Sequencing Identifies a Novel Germline Variant in PTK7 Gene in Familial Colorectal Cancer.

Colorectal cancer (CRC) is the third most frequently diagnosed malignancy worldwide. Only 5% of all CRC cases are due to germline mutations in known predisposition genes, and the remaining genetic burden still has to be discovered. In this study, we performed whole-exome sequencing on six members of a Polish family diagnosed with CRC and identified a novel germline variant in the protein tyrosine kinase 7 (inactive) gene (PTK7, ENST00000230419, V354M). Targeted screening of the variant in 1705 familial CRC cases and 1674 healthy elderly individuals identified the variant in an additional familial CRC case. Introduction of this variant in HT-29 cells resulted in increased cell proliferation, migration, and invasion; it also caused down-regulation of CREB, p21 and p53 mRNA and protein levels, and increased AKT phosphorylation. These changes indicated inhibition of apoptosis pathways and activation of AKT signaling. Our study confirmed the oncogenic function of PTK7 and supported its role in genetic predisposition of familial CRC.

Comprehensive genomic analysis of refractory multiple myeloma reveals a complex mutational landscape associated with drug resistance and novel therapeutic vulnerabilities.

The outcomes of patients with multiple myeloma (MM) refractory to immunomodulatory agents (IMiDs) and proteasome inhibitors (PIs) remain poor. In this study, we performed whole genome and transcriptome sequencing of 39 heavily pretreated relapsed/refractory MM (RRMM) patients to identify mechanisms of resistance and potential therapeutic targets. We observed a high mutational load and indications of increased genomic instability. Recurrently mutated genes in RRMM, which had not been previously reported or only observed at a lower frequency in newly diagnosed MM, included NRAS, BRAF, TP53, SLC4A7, MLLT4, EWSR1, HCFC2, and COPS3. We found multiple genomic regions with bi-allelic events affecting tumor suppressor genes and demonstrated a significant adverse impact of bi-allelic TP53 alterations on survival. With regard to potentially resistance conferring mutations, recurrently mutated gene networks included genes with relevance for PI and IMiD activity; the latter particularly affecting members of the Cereblon and the COP9 signalosome complex. We observed a major impact of signatures associated with exposure to melphalan or impaired DNA double-strand break homologous recombination repair in RRMM. The latter coincided with mutations in genes associated with PARP inhibitor sensitivity in 49% of RRMM patients; a finding with potential therapeutic implications. In conclusion, this comprehensive genomic characterization revealed a complex mutational and structural landscape in RRMM and highlights potential implications for therapeutic strategies.

Investigation of Rare Non-Coding Variants in Familial Multiple Myeloma.

Multiple myeloma (MM) is a plasma cell malignancy whereby a single clone of plasma cells over-propagates in the bone marrow, resulting in the increased production of monoclonal immunoglobulin. While the complex genetic architecture of MM is well characterized, much less is known about germline variants predisposing to MM. Genome-wide sequencing approaches in MM families have started to identify rare high-penetrance coding risk alleles. In addition, genome-wide association studies have discovered several common low-penetrance risk alleles, which are mainly located in the non-coding genome. Here, we further explored the genetic basis in familial MM within the non-coding genome in whole-genome sequencing data. We prioritized and characterized 150 upstream, 5' untranslated region (UTR) and 3' UTR variants from 14 MM families, including 20 top-scoring variants. These variants confirmed previously implicated biological pathways in MM development. Most importantly, protein network and pathway enrichment analyses also identified 10 genes involved in mitogen-activated protein kinase (MAPK) signaling pathways, which have previously been established as important MM pathways.

Identification of ZBTB26 as a Novel Risk Factor for Congenital Hypothyroidism.

Congenital primary hypothyroidism (CH; OMIM 218700) is characterized by an impaired thyroid development, or dyshormonogenesis, and can lead to intellectual disability and growth retardation if untreated. Most of the children with congenital hypothyroidism present thyroid dysgenesis, a developmental anomaly of the thyroid. Various genes have been associated with thyroid dysgenesis, but all known genes together can only explain a small number of cases. To identify novel genetic causes for congenital hypothyroidism, we performed trio whole-exome sequencing in an affected newborn and his unaffected parents. A predicted damaging de novo missense mutation was identified in the ZBTB26 gene (Zinc Finger A and BTB Domain containing 26). An additional cohort screening of 156 individuals with congenital thyroid dysgenesis identified two additional ZBTB26 gene variants of unknown significance. To study the underlying disease mechanism, morpholino knock-down of zbtb26 in Xenopus laevis was carried out, which demonstrated significantly smaller thyroid anlagen in knock-down animals at tadpole stage. Marker genes expressed in thyroid tissue precursors also indicated a specific reduction in the Xenopus ortholog of human Paired-Box-Protein PAX8, a transcription factor required for thyroid development, which could be rescued by adding zbtb26. Pathway and network analysis indicated network links of ZBTB26 to PAX8 and other genes involved in thyroid genesis and function. GWAS associations of ZBTB26 were found with height. Together, our study added a novel genetic risk factor to the list of genes underlying congenital primary hypothyroidism and provides additional support that de novo mutations, together with inherited variants, might contribute to the genetic susceptibility to CH.

Molecular profiling of pediatric meningiomas shows tumor characteristics distinct from adult meningiomas.

In contrast to adults, meningiomas are uncommon tumors in childhood and adolescence. Whether adult and pediatric meningiomas differ on a molecular level is unclear. Here we report detailed genomic analyses of 37 pediatric meningiomas by sequencing and DNA methylation profiling. Histologically, the series was dominated by meningioma subtypes with aggressive behavior, with 70% of patients suffering from WHO grade II or III meningiomas. The most frequent cytogenetic aberrations were loss of chromosomes 22 (23/37 [62%]), 1 (9/37 [24%]), 18 (7/37 [19%]), and 14 (5/37 [14%]). Tumors with NF2 alterations exhibited overall increased chromosomal instability. Unsupervised clustering of DNA methylation profiles revealed separation into three groups: designated group 1 composed of clear cell and papillary meningiomas, whereas group 2A comprised predominantly atypical meningiomas and group 2B enriched for rare high-grade subtypes (rhabdoid, chordoid). Meningiomas from NF2 patients clustered exclusively within groups 1 and 2A. When compared with a dataset of 105 adult meningiomas, the pediatric meningiomas largely grouped separately. Targeted panel DNA sequencing of 34 tumors revealed frequent NF2 alterations, while other typical alterations found in adult non-NF2 tumors were absent. These data demonstrate that pediatric meningiomas are characterized by molecular features distinct from adult tumors.

Combinations of Low-Frequency Genetic Variants Might Predispose to Familial Pancreatic Cancer.

Familial pancreatic cancer (FPC) is an established but rare inherited tumor syndrome that accounts for approximately 5% of pancreatic ductal adenocarcinoma (PDAC) cases. No major causative gene defect has yet been identified, but germline mutations in predisposition genes BRCA1/2, CDKN2A and PALB2 could be detected in 10-15% of analyzed families. Thus, the genetic basis of disease susceptibility in the majority of FPC families remains unknown. In an attempt to identify new candidate genes, we performed whole-genome sequencing on affected patients from 15 FPC families, without detecting BRCA1/2, CDKN2A or PALB2 mutations, using an Illumina based platform. Annotations from CADD, PolyPhen-2, SIFT, Mutation Taster and PROVEAN were used to assess the potential impact of a variant on the function of a gene. Variants that did not segregate with pancreatic disease in respective families were excluded. Potential predisposing candidate genes ATM, SUFU, DAB1, POLQ, FGFBP3, MAP3K3 and ACAD9 were identified in 7 of 15 families. All identified gene mutations segregated with pancreatic disease, but sometimes with incomplete penetrance. An analysis of up to 46 additional FPC families revealed that the identified gene mutations appeared to be unique in most cases, despite a potentially deleterious ACAD9 Ala326Thr germline variant, which occurred in 4 (8.7%) of 46 FPC families. Notably, affected PDAC patients within a family carried identical germline mutations in up to three different genes, e.g., DAB1, POLQ and FGFBP3. These results support the hypothesis that FPC is a highly heterogeneous polygenetic disease caused by low-frequency or rare variants.

Evidence That Non-Syndromic Familial Tall Stature Has an Oligogenic Origin Including Ciliary Genes.

Human growth is a complex trait. A considerable number of gene defects have been shown to cause short stature, but there are only few examples of genetic causes of non-syndromic tall stature. Besides rare variants with large effects and common risk alleles with small effect size, oligogenic effects may contribute to this phenotype. Exome sequencing was carried out in a tall male (height 3.5 SDS) and his parents. Filtered damaging variants with high CADD scores were validated by Sanger sequencing in the trio and three other affected and one unaffected family members. Network analysis was carried out to assess links between the candidate genes, and the transcriptome of murine growth plate was analyzed by microarray as well as RNA Seq. Heterozygous gene variants in CEP104, CROCC, NEK1, TOM1L2, and TSTD2 predicted as damaging were found to be shared between the four tall family members. Three of the five genes (CEP104, CROCC, and NEK1) belong to the ciliary gene family. All genes are expressed in mouse growth plate. Pathway and network analyses indicated close functional connections. Together, these data expand the spectrum of genes with a role in linear growth and tall stature phenotypes.

Mutational mechanisms shaping the coding and noncoding genome of germinal center derived B-cell lymphomas.

B cells have the unique property to somatically alter their immunoglobulin (IG) genes by V(D)J recombination, somatic hypermutation (SHM) and class-switch recombination (CSR). Aberrant targeting of these mechanisms is implicated in lymphomagenesis, but the mutational processes are poorly understood. By performing whole genome and transcriptome sequencing of 181 germinal center derived B-cell lymphomas (gcBCL) we identified distinct mutational signatures linked to SHM and CSR. We show that not only SHM, but presumably also CSR causes off-target mutations in non-IG genes. Kataegis clusters with high mutational density mainly affected early replicating regions and were enriched for SHM- and CSR-mediated off-target mutations. Moreover, they often co-occurred in loci physically interacting in the nucleus, suggesting that mutation hotspots promote increased mutation targeting of spatially co-localized loci (termed hypermutation by proxy). Only around 1% of somatic small variants were in protein coding sequences, but in about half of the driver genes, a contribution of B-cell specific mutational processes to their mutations was found. The B-cell-specific mutational processes contribute to both lymphoma initiation and intratumoral heterogeneity. Overall, we demonstrate that mutational processes involved in the development of gcBCL are more complex than previously appreciated, and that B cell-specific mutational processes contribute via diverse mechanisms to lymphomagenesis.

A Novel Low-Risk Germline Variant in the SH2 Domain of the SRC Gene Affects Multiple Pathways in Familial Colorectal Cancer.

Colorectal cancer (CRC) shows one of the largest proportions of familial cases among different malignancies, but only 5-10% of all CRC cases are linked to mutations in established predisposition genes. Thus, familial CRC constitutes a promising target for the identification of novel, high- to moderate-penetrance germline variants underlying cancer susceptibility by next generation sequencing. In this study, we performed whole genome sequencing on three members of a family with CRC aggregation. Subsequent integrative in silico analysis using our in-house developed variant prioritization pipeline resulted in the identification of a novel germline missense variant in the SRC gene (V177M), a proto-oncogene highly upregulated in CRC. Functional validation experiments in HT-29 cells showed that introduction of SRCV177M resulted in increased cell proliferation and enhanced protein expression of phospho-SRC (Y419), a potential marker for SRC activity. Upregulation of paxillin, ß-Catenin, and STAT3 mRNA levels, increased levels of phospho-ERK, CREB, and CCND1 proteins and downregulation of the tumor suppressor p53 further proposed the activation of several pathways due to the SRCV177M variant. The findings of our pedigree-based study contribute to the exploration of the genetic background of familial CRC and bring insights into the molecular basis of upregulated SRC activity and downstream pathways in colorectal carcinogenesis.

Whole Exome Sequencing Identifies APCDD1 and HDAC5 Genes as Potentially Cancer Predisposing in Familial Colorectal Cancer.

Germline mutations in predisposition genes account for only 20% of all familial colorectal cancers (CRC) and the remaining genetic burden may be due to rare high- to moderate-penetrance germline variants that are not explored. With the aim of identifying such potential cancer-predisposing variants, we performed whole exome sequencing on three CRC cases and three unaffected members of a Polish family and identified two novel heterozygous variants: a coding variant in APC downregulated 1 gene (APCDD1, p.R299H) and a non-coding variant in the 5' untranslated region (UTR) of histone deacetylase 5 gene (HDAC5). Sanger sequencing confirmed the variants segregating with the disease and Taqman assays revealed 8 additional APCDD1 variants in a cohort of 1705 familial CRC patients and no further HDAC5 variants. Proliferation assays indicated an insignificant proliferative impact for the APCDD1 variant. Luciferase reporter assays using the HDAC5 variant resulted in an enhanced promoter activity. Targeting of transcription factor binding sites of SNAI-2 and TCF4 interrupted by the HDAC5 variant showed a significant impact of TCF4 on promoter activity of mutated HDAC5. Our findings contribute not only to the identification of unrecognized genetic causes of familial CRC but also underline the importance of 5'UTR variants affecting transcriptional regulation and the pathogenesis of complex disorders.