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.

Deconvolution of sarcoma methylomes reveals varying degrees of immune cell infiltrates with association to genomic aberrations.

BACKGROUND: Soft-tissue sarcomas (STS) are a heterogeneous group of mesenchymal tumors for which response to immunotherapies is not well established. Therefore, it is important to risk-stratify and identify STS patients who will most likely benefit from these treatments. RESULTS: To reveal shared and distinct methylation signatures present in STS, we performed unsupervised deconvolution of DNA methylation data from the TCGA sarcoma and an independent validation cohort. We showed that leiomyosarcoma can be subclassified into three distinct methylation groups. More importantly, we identified a component associated with tumor-infiltrating leukocytes, which suggests varying degrees of immune cell infiltration in STS subtypes and an association with prognosis. We further investigated the genomic alterations that may influence tumor infiltration by leukocytes including RB1 loss in undifferentiated pleomorphic sarcomas and ELK3 amplification in dedifferentiated liposarcomas. CONCLUSIONS: In summary, we have leveraged unsupervised methylation-based deconvolution to characterize the immune compartment and molecularly stratify subtypes in STS, which may benefit precision medicine in the future.

TelomereHunter - in silico estimation of telomere content and composition from cancer genomes.

BACKGROUND: Establishment of telomere maintenance mechanisms is a universal step in tumor development to achieve replicative immortality. These processes leave molecular footprints in cancer genomes in the form of altered telomere content and aberrations in telomere composition. To retrieve these telomere characteristics from high-throughput sequencing data the available computational approaches need to be extended and optimized to fully exploit the information provided by large scale cancer genome data sets. RESULTS: We here present TelomereHunter, a software for the detailed characterization of telomere maintenance mechanism footprints in the genome. The tool is implemented for the analysis of large cancer genome cohorts and provides a variety of diagnostic diagrams as well as machine-readable output for subsequent analysis. A novel key feature is the extraction of singleton telomere variant repeats, which improves the identification and subclassification of the alternative lengthening of telomeres phenotype. We find that whole genome sequencing-derived telomere content estimates strongly correlate with telomere qPCR measurements (r = 0.94). For the first time, we determine the correlation of in silico telomere content quantification from whole genome sequencing and whole genome bisulfite sequencing data derived from the same tumor sample (r = 0.78). An analogous comparison of whole exome sequencing data and whole genome sequencing data measured slightly lower correlation (r = 0.79). However, this is considerably improved by normalization with matched controls (r = 0.91). CONCLUSIONS: TelomereHunter provides new functionality for the analysis of the footprints of telomere maintenance mechanisms in cancer genomes. Besides whole genome sequencing, whole exome sequencing and whole genome bisulfite sequencing are suited for in silico telomere content quantification, especially if matched control samples are available. The software runs under a GPL license and is available at https://www.dkfz.de/en/applied-bioinformatics/telomerehunter/telomerehunter.html .

Tpbpa-Cre-mediated deletion of TFAP2C leads to deregulation of Cdkn1a, Akt1 and the ERK pathway, causing placental growth arrest.

Loss of TFAP2C in mouse leads to developmental defects in the extra-embryonic compartment with lethality at embryonic day (E)7.5. To investigate the requirement of TFAP2C in later placental development, deletion of TFAP2C was induced throughout extra-embryonic ectoderm at E6.5, leading to severe placental abnormalities caused by reduced trophoblast population and resulting in embryonic retardation by E8.5. Deletion of TFAP2C in TPBPA(+) progenitors at E8.5 results in growth arrest of the junctional zone. TFAP2C regulates its target genes Cdkn1a (previously p21) and Dusp6, which are involved in repression of MAPK signaling. Loss of TFAP2C reduces activation of ERK1/2 in the placenta. Downregulation of Akt1 and reduced activation of phosphorylated AKT in the mutant placenta are accompanied by impaired glycogen synthesis. Loss of TFAP2C led to upregulation of imprinted gene H19 and downregulation of Slc38a4 and Ascl2. The placental insufficiency post E16.5 causes fetal growth restriction, with 19% lighter mutant pups. Knockdown of TFAP2C in human trophoblast choriocarcinoma JAr cells inhibited MAPK and AKT signaling. Thus, we present a model where TFAP2C in trophoblasts controls proliferation by repressing Cdkn1a and activating the MAPK pathway, further supporting differentiation of glycogen cells by activating the AKT pathway.

Heterozygous FGF8 mutations in patients presenting cryptorchidism and multiple VATER/VACTERL features without limb anomalies.

BACKGROUND: The acronym VATER/VACTERL association describes the combination of at least three of the following cardinal features: vertebral defects, anorectal malformations, cardiac defects, tracheoesophageal fistula with or without esophageal atresia, renal malformations, and limb defects. Although fibroblast growth factor-8 (FGF8) mutations have mainly found in patients with Kallmann syndrome, mice with a hypomorphic Fgf8 allele or complete gene invalidation display, aside from gonadotropin-releasing hormone deficiency, parts or even the entire spectrum of human VATER/VACTERL association. METHODS: We performed FGF8 gene analysis in 49 patients with VATER/VACTERL association and 27 patients presenting with a VATER/VACTERL-like phenotype (two cardinal features). RESULTS: We identified two heterozygous FGF8 mutations in patients displaying either VATER/VACTERL association (p.Gly29_Arg34dup) or a VATER/VACTERL-like phenotype (p.Pro26Leu) without limb anomalies. Whereas the duplication mutation has not been reported before, p.Pro26Leu was once observed in a Kallmann syndrome patient. Both our patients had additional bilateral cryptorchidism, a key phenotypic feature in males with FGF8 associated Kallmann syndrome. Each mutation was paternally inherited. Besides delayed puberty in both and additional unilateral cryptorchidism in one of the fathers, they were otherwise healthy. Serum hormone levels downstream the gonadotropin-releasing hormone in both patients and their fathers were within normal range. CONCLUSION: Our results suggest FGF8 mutations to contribute to the formation of the VATER/VACTERL association. Further studies are needed to support this observation.

Isolated bladder exstrophy associated with a de novo 0.9 Mb microduplication on chromosome 19p13.12.

BACKGROUND: The exstrophy-epispadias complex (BEEC) is a urogenital birth defect of varying severity. The causes of the BEEC are likely to be heterogeneous, with individual environmental or genetic risk factors still being largely unknown. In this study, we aimed to identify de novo causative copy number variations (CNVs) that contribute to the BEEC. METHODS: Array-based molecular karyotyping was performed to screen 110 individuals with BEEC. Promising CNVs were tested for de novo occurrence by investigating parental DNAs. Genes located in regions of rearrangements were prioritized through expression analysis in mice to be sequenced in the complete cohort, to identify high-penetrance mutations involving small sequence changes. RESULTS: A de novo 0.9 Mb microduplication involving chromosomal region 19p13.12 was identified in a single patient. This region harbors 20 validated RefSeq genes, and in situ hybridization data showed specific expression of the Wiz gene in regions surrounding the cloaca and the rectum between GD 9.5 and 13.5. Sanger sequencing of the complete cohort did not reveal any pathogenic alterations affecting the coding region of WIZ. CONCLUSIONS: The present study suggests chromosomal region 19p13.12 as possibly involved in the development of CBE, but further studies are needed to prove a causal relation. The spatiotemporal expression patterns determined for the genes encompassed suggest a role for Wiz in the development of the phenotype. Our mutation screening, however, could not confirm that WIZ mutations are a frequent cause of CBE, although rare mutations might be detectable in larger patient samples. 19p13.12, microduplication, bladder exstrophy-epispadias complex, array-based molecular karyotyping, in situ hybridization analysis, copy number variations, WIZ gene.

Liddle syndrome in a Serbian family and literature review of underlying mutations.

UNLABELLED: Severe and reproducible low-renin hypertension responsive to salt restriction and amiloride-thiazide therapy in a 13-year-old otherwise asymptomatic boy suggested Liddle syndrome. This assumption was strengthened by a positive family history of hypertension poorly responsive to conventional treatment or sudden deaths under 40 years of age in four generations. DNA analysis of the beta and gamma subunits of the epithelial sodium channel revealed a heterozygous mutation c.C1852T (p.Pro618Ser) in the SCNN1B gene in the patient and in both his hypertensive mother and uncle. A PubMed search revealed 21 different disease-causing mutations reported to date, all but two clustering in the cytoplasmic C-terminal regions of either beta (16 mutations) or gamma (5) subunit, leading to a three- to eightfold increase in the amiloride-sensitive sodium current. Inter- and intrafamilial variability in both hypertension and hypokalemia were disclosed, which may not be obligatory among the subjects carrying a Liddle mutation. CONCLUSION: Liddle syndrome should be considered as a cause of hypertension in children or adolescents particularly with suppressed renin activity. Early diagnosis and appropriately tailored treatment avoid complications of long-term unrecognized or inappropriately managed hypertension.

Familial occurrence of the VATER/VACTERL association.

The acronym VATER/VACTERL association is used to describe the non-random co-occurrence of vertebral defects (V), anorectal malformations (A), cardiac defects (C), tracheo-esophageal fistula with or without esophageal atresia (TE), renal malformations (R), and limb defects (L). We report a familial case of VATER/VACTERL association in which both the index case and her maternal uncle displayed four major component features of the disorder. A systematic literature search identified 12 previously described familial cases. However, on comparison, both members fulfilled the diagnostic criteria for VATER/VACTERL association only in one instance, and ours is the second such report. Although, a SNP array-based analysis identified no causal genomic alteration, the findings in the present family suggest that genetic factors are implicated in the development of the disorder.

Gene expression-based prediction of pazopanib efficacy in sarcoma.

BACKGROUND: The multi-receptor tyrosine kinase inhibitor pazopanib is approved for the treatment of advanced soft-tissue sarcoma and has also shown activity in other sarcoma subtypes. However, its clinical efficacy is highly variable, and no reliable predictors exist to select patients who are likely to benefit from this drug. PATIENTS AND METHODS: We analysed the molecular profiles and clinical outcomes of patients with pazopanib-treated sarcoma enrolled in a prospective observational study by the German Cancer Consortium, DKTK MASTER, that employs whole-genome/exome sequencing and transcriptome sequencing to inform the care of young adults with advanced cancer across histology and patients with rare cancers. RESULTS: Among 109 patients with available whole-genome/exome sequencing data, there was no correlation between clinical parameters, specific genetic alterations or mutational signatures and clinical outcome. In contrast, the analysis of a subcohort of 62 patients who underwent molecular analysis before pazopanib treatment and had transcriptome sequencing data available showed that mRNA levels of NTRK3 (hazard ratio [HR] = 0.53, p = 0.021), IGF1R (HR = 1.82, p = 0.027) and KDR (HR = 0.50, p = 0.011) were independently associated with progression-free survival (PFS). Based on the expression of these receptor tyrosine kinase genes, i.e. the features NTRK3-high, IGF1R-low and KDR-high, we developed a pazopanib efficacy predictor that stratified patients into three groups with significantly different PFS (p < 0.0001). Application of the pazopanib efficacy predictor to an independent cohort of patients with pazopanib-treated sarcoma from DKTK MASTER (n = 43) confirmed its potential to separate patient groups with significantly different PFS (p = 0.02), whereas no such association was observed in patients with sarcoma from DKTK MASTER (n = 97) or The Cancer Genome Atlas sarcoma cohort (n = 256) who were not treated with pazopanib. CONCLUSION: A score based on the combined expression of NTRK3, IGF1R and KDR allows the identification of patients with sarcoma and with good, intermediate and poor outcome following pazopanib therapy and warrants prospective investigation as a predictive tool to optimise the use of this drug in the clinic.

Macrophages/Microglia Represent the Major Source of Indolamine 2,3-Dioxygenase Expression in Melanoma Metastases of the Brain.

The manifestation of brain metastases in patients with advanced melanoma is a common event that limits patient's survival and quality of life. The immunosuppressive properties of the brain parenchyma are very different compared to the rest of the body, making it plausible that the current success of cancer immunotherapies is specifically limited here. In melanoma brain metastases, the reciprocal interplay between immunosuppressive mediators such as indoleamine 2, 3-dioxygenase (IDO) or programmed cell death-ligand 1 (PD-L1) in the context of neoplastic transformation are far from being understood. Therefore, we analyzed the immunoreactive infiltrate (CD45, CD3, CD8, Forkhead box P3 [FoxP3], CD11c, CD23, CD123, CD68, Allograft Inflammatory factor 1[AIF-1]) and PD-L1 with respect to IDO expression and localization in melanoma brain metastases but also in matched metastases at extracranial sites to correlate intra- and interpatient data with therapy response and survival. Comparative tissue analysis identified macrophages/microglia as the major source of IDO expression in melanoma brain metastases. In contrast to the tumor infiltrating lymphocytes, melanoma cells per se exhibited low IDO expression levels paralleled by cell surface presentation of PD-L1 in intracranial metastases. Absolute numbers and pattern of IDO-expressing cells in metastases of the brain correlated with recruitment and localization of CD8+ T cells, implicating dynamic impact on the regulation of T cell function in the brain parenchyma. However, paired analysis of matched intra- and extracranial metastases identified significantly lower fractions of cytotoxic CD8+ T cells in intracranial metastases while all other immune cell populations remain unchanged. In line with the already established clinical benefit for PD-L1 expression in extracranial melanoma metastases, Kaplan-Meier analyses correlated PD-L1 expression in brain metastases with favorable outcome in advanced melanoma patients undergoing immune checkpoint therapy. In summary, our data provide new insights into the landscape of immunosuppressive factors in melanoma brain metastases that may be useful in the implication of novel therapeutic strategies for patients undergoing cancer immunotherapy.

Integrative genomic and transcriptomic analysis of leiomyosarcoma.

Leiomyosarcoma (LMS) is an aggressive mesenchymal malignancy with few therapeutic options. The mechanisms underlying LMS development, including clinically actionable genetic vulnerabilities, are largely unknown. Here we show, using whole-exome and transcriptome sequencing, that LMS tumors are characterized by substantial mutational heterogeneity, near-universal inactivation of TP53 and RB1, widespread DNA copy number alterations including chromothripsis, and frequent whole-genome duplication. Furthermore, we detect alternative telomere lengthening in 78% of cases and identify recurrent alterations in telomere maintenance genes such as ATRX, RBL2, and SP100, providing insight into the genetic basis of this mechanism. Finally, most tumors display hallmarks of "BRCAness", including alterations in homologous recombination DNA repair genes, multiple structural rearrangements, and enrichment of specific mutational signatures, and cultured LMS cells are sensitive towards olaparib and cisplatin. This comprehensive study of LMS genomics has uncovered key biological features that may inform future experimental research and enable the design of novel therapies.

PD-L1 (CD274) copy number gain, expression, and immune cell infiltration as candidate predictors for response to immune checkpoint inhibitors in soft-tissue sarcoma.

Soft-tissue sarcomas (STS) are rare malignancies that account for 1% of adult cancers and comprise more than 50 entities. Current therapeutic options for advanced-stage STS are limited. Immune checkpoint inhibitors targeting the PD-1/PD-L1 signaling axis are being explored as new treatment modality in STS; however, the determinants of response to these agents are largely unknown. Using the sarcoma data set of The Cancer Genome Altas (TCGA) and an independent cohort of untreated high-grade STS, we analyzed DNA copy number status and mRNA expression of PD-L1 in a total of 335 STS cases. Copy number gains (CNG) were detected in 54 TCGA cases (21.1%), of which 21 (8.2%) harbored focal PD-L1 CNG and that were most prevalent in myxofibrosarcoma (35%) and undifferentiated pleomorphic sarcoma (34%). In the untreated high-grade STS cohort, we detected CNG in six cases (7.6%). Analysis of co-amplified genes identified a 5.6-Mb core region comprising 27 genes, including JAK2. Patients with PD-L1 CNG had higher PD-L1 expression compared with STS without CNG (fold change, 1.8; p = 0.02), an effect that was most pronounced in the setting of focal PD-L1 CNG (fold change, 3.0; p = 0.0027). STS with PD-L1 CNG showed a significantly higher mutational load compared with tumors with a diploid PD-L1 locus (median number of mutated genes; 58 vs. 40; p = 3.6E-06), and PD-L1 CNG were associated with inferior survival (HR = 1.82; p = 0.025). In contrast, T-cell infiltrates quantified by mRNA expression of CD3Z were associated with improved survival (HR = 0.88; p = 0.024) and consequently influenced the prognostic power of PD-L1 CNG, with low CD3Z levels conferring poor survival in cases with PD-L1 CNG (HR = 1.8; p = 0.049). These data demonstrate that PD-L1 GNG and elevated expression of PD-L1 occur in a substantial proportion of STS, have prognostic impact that is modulated by T-cell infiltrates, and thus warrant investigation as response predictors for immune checkpoint inhibition.

Targeting Fibroblast Growth Factor Receptor 1 for Treatment of Soft-Tissue Sarcoma.

Purpose: Altered FGFR1 signaling has emerged as a therapeutic target in epithelial malignancies. In contrast, the role of FGFR1 in soft-tissue sarcoma (STS) has not been established. Prompted by the detection and subsequent therapeutic inhibition of amplified FGFR1 in a patient with metastatic leiomyosarcoma, we investigated the oncogenic properties of FGFR1 and its potential as a drug target in patients with STS.Experimental Design: The frequency of FGFR1 amplification and overexpression, as assessed by FISH, microarray-based comparative genomic hybridization and mRNA expression profiling, SNP array profiling, and RNA sequencing, was determined in three patient cohorts. The sensitivity of STS cell lines with or without FGFR1 alterations to genetic and pharmacologic FGFR1 inhibition and the signaling pathways engaged by FGFR1 were investigated using viability assays, colony formation assays, and biochemical analysis.Results: Increased FGFR1 copy number was detected in 74 of 190 (38.9%; cohort 1), 13 of 79 (16.5%; cohort 2), and 80 of 254 (31.5%; cohort 3) patients. FGFR1 overexpression occurred in 16 of 79 (20.2%, cohort 2) and 39 of 254 (15.4%; cohort 3) patients. Targeting of FGFR1 by RNA interference and small-molecule inhibitors (PD173074, AZD4547, BGJ398) revealed that the requirement for FGFR1 signaling in STS cells is dictated by FGFR1 expression levels, and identified the MAPK-ERK1/2 axis as critical FGFR1 effector pathway.Conclusions: These data identify FGFR1 as a driver gene in multiple STS subtypes and support FGFR1 inhibition, guided by patient selection according to the FGFR1 expression and monitoring of MAPK-ERK1/2 signaling, as a therapeutic option in this challenging group of diseases. Clin Cancer Res; 23(4); 962-73. ©2016 AACR.

Involvement of the WNT and FGF signaling pathways in non-isolated anorectal malformations: sequencing analysis of WNT3A, WNT5A, WNT11, DACT1, FGF10, FGFR2 and the T gene.

Anorectal malformations (ARMs) comprise a broad spectrum of anomalies, including anal atresia, congenital anal fistula and persistence of the cloaca. Research suggests that genetic factors play an important role in ARM development. However, few genetic variants have been identified. Embryogenesis is orchestrated by crosstalk of the wingless-type MMTV integration site family (WNT) and fibroblast growth factor (FGF) signaling pathways in a process that involves several intracellular cascades. Studies in mice have implicated several genes from these pathways in the etiology of ARMs. We performed sequencing analysis of seven of these previously reported genes in 78 patients with ARMs occurring within the context of at least one additional congenital anomaly. No associations were identified with variants in WNT3A, WNT5A, WNT11, DACT1, FGF10 or the T gene. In the FGFR2 gene, three novel heterozygous nucleotide substitutions were identified. Further investigations, including the study of family members, revealed that these variants were not causally related to the phenotype in the present ARM cohort. Mutations in the seven investigated genes may nonetheless be a cause of ARMs in rare cases. However, further studies should consider genes encoding other proteins in the WNT/FGF signaling pathways as possible candidates.