TarPan: an easily adaptable targeted sequencing panel viewer for research and clinical use.

BACKGROUND: The study of cancer genomics continually matures as the number of patient samples sequenced increases. As more data is generated, oncogenic drivers for specific cancer types are discovered along with their associated risks. This in turn leads to potential treatment strategies that pave the way to precision medicine. However, significant financial and analytical barriers make it infeasible to sequence the entire genome of every patient. In contrast, targeted sequencing panels give reliable information on relevant portions of the genome at a fiscally responsible cost. Therefore, we have created the Targeted Panel (TarPan) Viewer, a software tool, to investigate this type of data. RESULTS: TarPan Viewer helps investigators understand data from targeted sequencing data by displaying the information through a web browser interface. Through this interface, investigators can easily observe copy number changes, mutations, and structural events in cancer samples. The viewer runs in R Shiny with a robust SQLite backend and its input is generated from bioinformatic algorithms reliably described in the literature. Here we show the results from using TarPan Viewer on publicly available follicular lymphoma, breast cancer, and multiple myeloma data. In addition, we have tested and utilized the viewer internally, and this data has been used in high-impact peer-reviewed publications. CONCLUSIONS: We have designed a flexible, simple to setup viewer that is easily adaptable to any type of cancer targeted sequencing, and has already proven its use in a research laboratory environment. Further, we believe with deeper sequencing and/or more targeted application it could be of use in the clinic in conjunction with an appropriate targeted sequencing panel as a cost-effective diagnostic test, especially in cancers such as acute leukemia or diffuse large B-cell lymphoma that require rapid interventions.

Identification of novel mutational drivers reveals oncogene dependencies in multiple myeloma.

Understanding the profile of oncogene and tumor suppressor gene mutations with their interactions and impact on the prognosis of multiple myeloma (MM) can improve the definition of disease subsets and identify pathways important in disease pathobiology. Using integrated genomics of 1273 newly diagnosed patients with MM, we identified 63 driver genes, some of which are novel, including IDH1, IDH2, HUWE1, KLHL6, and PTPN11 Oncogene mutations are significantly more clonal than tumor suppressor mutations, indicating they may exert a bigger selective pressure. Patients with more driver gene abnormalities are associated with worse outcomes, as are identified mechanisms of genomic instability. Oncogenic dependencies were identified between mutations in driver genes, common regions of copy number change, and primary translocation and hyperdiploidy events. These dependencies included associations with t(4;14) and mutations in FGFR3, DIS3, and PRKD2; t(11;14) with mutations in CCND1 and IRF4; t(14;16) with mutations in MAF, BRAF, DIS3, and ATM; and hyperdiploidy with gain 11q, mutations in FAM46C, and MYC rearrangements. These associations indicate that the genomic landscape of myeloma is predetermined by the primary events upon which further dependencies are built, giving rise to a nonrandom accumulation of genetic hits. Understanding these dependencies may elucidate potential evolutionary patterns and lead to better treatment regimens.

Microhomology-mediated end joining drives complex rearrangements and overexpression of MYC and PVT1 in multiple myeloma.

MYC is a widely acting transcription factor and its deregulation is a crucial event in many human cancers. MYC is important biologically and clinically in multiple myeloma, but the mechanisms underlying its dysregulation are poorly understood. We show that MYC rearrangements are present in 36.0% of newly diagnosed myeloma patients, as detected in the largest set of next generation sequencing data to date (n=1,267). Rearrangements were complex and associated with increased expression of MYC and PVT1, but not other genes at 8q24. The highest effect on gene expression was detected in cases where the MYC locus is juxtaposed next to super-enhancers associated with genes such as IGH, IGK, IGL, TXNDC5/BMP6, FAM46C and FOXO3 We identified three hotspots of recombination at 8q24, one of which is enriched for IGH-MYC translocations. Breakpoint analysis indicates primary myeloma rearrangements involving the IGH locus occur through non-homologous end joining, whereas secondary MYC rearrangements occur through microhomology-mediated end joining. This mechanism is different to lymphomas, where non-homologous end joining generates MYC rearrangements. Rearrangements resulted in overexpression of key genes and chromatin immunoprecipitation-sequencing identified that HK2, a member of the glucose metabolism pathway, is directly over-expressed through binding of MYC at its promoter.

Daratumumab and dexamethasone is safe and effective for triple refractory myeloma patients: final results of the IFM 2014-04 (Etoile du Nord) trial.

Single agent daratumumab has shown clinical activity in relapsed, refractory multiple myeloma (RRMM). The Intergroupe Francophone du Myélome 2014-04 trial was designed to further investigate daratumumab in combination with dexamethasone in triple RRMM patients. Patients received daratumumab infusions in combination with weekly dexamethasone until disease progression or unacceptable toxicity. Fifty-seven patients were included in the trial and evaluable for response. The overall response rate and the clinical benefit rate were 33% (n = 19) and 48% (n = 27), respectively. Five (8·8%) patients achieved a very good partial response or better. The median time to response was 4 weeks. For responding patients, the median progression-free survival was 6·6 months, compared to 3·7 months (3·0-5·5) for those with a minimal or stable disease. The median overall survival (OS) for all patients was 16·7 months (11·2-24·0). For responding patients, the median OS was 23·23 months, whereas that of patients with progressive disease was 2·97 months. The incidence of infusion-related reactions was 37%; all cases were manageable and did not lead to dose reduction or permanent treatment discontinuation. These data demonstrate that treatment with daratumumab and dexamethasone results in a meaningful long-term benefit with an acceptable safety profile for patients with triple RRMM.

Genomic characterization of biliary tract cancers identifies driver genes and predisposing mutations.

BACKGROUND & AIMS: Biliary tract cancers (BTCs) are clinically and pathologically heterogeneous and respond poorly to treatment. Genomic profiling can offer a clearer understanding of their carcinogenesis, classification and treatment strategy. We performed large-scale genome sequencing analyses on BTCs to investigate their somatic and germline driver events and characterize their genomic landscape. METHODS: We analyzed 412 BTC samples from Japanese and Italian populations, 107 by whole-exome sequencing (WES), 39 by whole-genome sequencing (WGS), and a further 266 samples by targeted sequencing. The subtypes were 136 intrahepatic cholangiocarcinomas (ICCs), 101 distal cholangiocarcinomas (DCCs), 109 peri-hilar type cholangiocarcinomas (PHCs), and 66 gallbladder or cystic duct cancers (GBCs/CDCs). We identified somatic alterations and searched for driver genes in BTCs, finding pathogenic germline variants of cancer-predisposing genes. We predicted cell-of-origin for BTCs by combining somatic mutation patterns and epigenetic features. RESULTS: We identified 32 significantly and commonly mutated genes including TP53, KRAS, SMAD4, NF1, ARID1A, PBRM1, and ATR, some of which negatively affected patient prognosis. A novel deletion of MUC17 at 7q22.1 affected patient prognosis. Cell-of-origin predictions using WGS and epigenetic features suggest hepatocyte-origin of hepatitis-related ICCs. Deleterious germline mutations of cancer-predisposing genes such as BRCA1, BRCA2, RAD51D, MLH1, or MSH2 were detected in 11% (16/146) of BTC patients. CONCLUSIONS: BTCs have distinct genetic features including somatic events and germline predisposition. These findings could be useful to establish treatment and diagnostic strategies for BTCs based on genetic information. LAY SUMMARY: We here analyzed genomic features of 412 BTC samples from Japanese and Italian populations. A total of 32 significantly and commonly mutated genes were identified, some of which negatively affected patient prognosis, including a novel deletion of MUC17 at 7q22.1. Cell-of-origin predictions using WGS and epigenetic features suggest hepatocyte-origin of hepatitis-related ICCs. Deleterious germline mutations of cancer-predisposing genes were detected in 11% of patients with BTC. BTCs have distinct genetic features including somatic events and germline predisposition.

Whole genome sequencing discriminates hepatocellular carcinoma with intrahepatic metastasis from multi-centric tumors.

BACKGROUND & AIMS: Patients with hepatocellular carcinoma (HCC) have a high-risk of multi-centric (MC) tumor occurrence due to a strong carcinogenic background in the liver. In addition, they have a high risk of intrahepatic metastasis (IM). Liver tumors withIM or MC are profoundly different in their development and clinical outcome. However, clinically or pathologically discriminating between IM and MC can be challenging. This study investigated whether IM or MC could be diagnosed at the molecular level. METHODS: We performed whole genome and RNA sequencing analyses of 49 tumors including two extra-hepatic metastases, and one nodule-in-nodule tumor from 23 HCC patients. RESULTS: Sequencing-based molecular diagnosis using somatic single nucleotide variation information showed higher sensitivity compared to previous techniques due to the inclusion of a larger number of mutation events. This proved useful in cases, which showed inconsistent clinical diagnoses. In addition, whole genome sequencing offered advantages in profiling of other genetic alterations, such as structural variations, copy number alterations, and variant allele frequencies, and helped to confirm the IM/MCdiagnosis. Divergent alterations between IM tumors with sorafenib treatment, long time-intervals, or tumor-in-tumor nodules indicated high intra-tumor heterogeneity, evolution, and clonal switching of liver cancers. CONCLUSIONS: It is important to analyze the differences between IM tumors, in addition to IM/MC diagnosis, before selecting a therapeutic strategy for multiple tumors in the liver. LAY SUMMARY: Whole genome sequencing of multiple liver tumors enabled the accuratediagnosis ofmulti-centric occurrence and intrahepatic metastasis using somatic single nucleotide variation information. In addition, genetic discrepancies between tumors help us to understand the physical changes during recurrence and cancer spread.

Coexistent hyperdiploidy does not abrogate poor prognosis in myeloma with adverse cytogenetics and may precede IGH translocations.

The acquisition of the cytogenetic abnormalities hyperdiploidy or translocations into the immunoglobulin gene loci are considered as initiating events in the pathogenesis of myeloma and were often assumed to be mutually exclusive. These lesions have clinical significance; hyperdiploidy or the presence of the t(11;14) translocation is associated with a favorable outcome, whereas t(4;14), t(14;16), and t(14;20) are unfavorable. Poor outcomes are magnified when lesions occur in association with other high-risk features, del17p and +1q. Some patients have coexistence of both good and poor prognostic lesions, and there has been no consensus on their risk status. To address this, we have investigated their clinical impact using cases in the Myeloma IX study (ISRCTN68454111) and shown that the coexistence of hyperdiploidy or t(11;14) does not abrogate the poor prognosis associated with adverse molecular lesions, including translocations. We have also used single-cell analysis to study cases with coexistent translocations and hyperdiploidy to determine how these lesions cosegregate within the clonal substructure, and we have demonstrated that hyperdiploidy may precede IGH translocation in a proportion of patients. These findings have important clinical and biological implications, as we conclude patients with coexistence of adverse lesions and hyperdiploidy should be considered high risk and treated accordingly.

The spectrum of somatic mutations in monoclonal gammopathy of undetermined significance indicates a less complex genomic landscape than that in multiple myeloma.

Monoclonal gammopathy of undetermined significance is a pre-malignant precursor of multiple myeloma with a 1% risk of progression per year. Although targeted analyses have shown the presence of specific genetic abnormalities such as IGH translocations, RB1 deletion, 1q gain, hyperdiploidy or RAS gene mutations, little is known about the molecular mechanism of malignant transformation. We performed whole exome sequencing together with comparative genomic hybridization plus single nucleotide polymorphism array analysis in 33 flow-cytometry-separated abnormal plasma cell samples from patients with monoclonal gammopathy of undetermined significance to describe somatic gene mutations and chromosome changes at the genome-wide level. Non-synonymous mutations and copy-number alterations were present in 97.0% and in 60.6% of cases, respectively. Importantly, the number of somatic mutations was significantly lower in monoclonal gammopathy of undetermined significance than in myeloma (P<10-4) and we identified six genes that were significantly mutated in myeloma (KRAS, NRAS, DIS3, HIST1H1E, EGR1 and LTB) within the monoclonal gammopathy of undetermined significance dataset. We also found a positive correlation with increasing chromosome changes and somatic gene mutations. IGH translocations, comprising t(4;14), t(11;14), t(14;16) and t(14;20), were present in 27.3% of cases and in a similar frequency to myeloma, consistent with the primary lesion hypothesis. MYC translocations and TP53 deletions or mutations were not detected in samples from patients with monoclonal gammopathy of undetermined significance, indicating that they may be drivers of progression to myeloma. Data from this study show that monoclonal gammopathy of undetermined significance is genetically similar to myeloma, however overall genetic abnormalities are present at significantly lower levels in monoclonal gammopathy of undetermined significant than in myeloma.

Genomewide profiling of copy-number alteration in monoclonal gammopathy of undetermined significance.

Monoclonal gammopathy of undetermined significance (MGUS) is a benign condition with an approximate 1% annual risk of symptomatic plasma cell disorder development, mostly to multiple myeloma (MM). We performed genomewide screening of copy-number alterations (CNAs) in 90 MGUS and 33 MM patients using high-density DNA microarrays. We identified CNAs in a smaller proportion of MGUS (65.6%) than in MM (100.0%, P = 1.31 × 10-5 ) and showed median number of CNAs is lower in MGUS (3, range 0-22) than in MM (13, range 4-38, P = 1.82 × 10-10 ). In the MGUS cohort, the most frequent losses were located at 1p (5.6%), 6q (6.7%), 13q (30.0%), 14q (14.4%), 16q (8.9%), 21q (5.6%), and gains at 1q (23.3%), 2p (6.7%), 6p (13.3%), and Xq (7.8%). Hyperdiploidy was detected in 38.9% of MGUS cases, and the most frequent whole chromosome gains were 3 (25.6%), 5 (23.3%), 9 (37.8%), 15 (23.3%), and 19 (32.2%). We also identified CNAs such as 1p, 6q, 8p, 12p, 13q, 16q losses, 1q gain and hypodiploidy, which are potentially associated with an adverse prognosis in MGUS. In summary, we showed that MGUS is similar to MM in that it is a genetically heterogeneous disorder, but overall cytogenetic instability is lower than in MM, which confirms that genetic abnormalities play important role in monoclonal gammopathies.

A molecular diagnostic approach able to detect the recurrent genetic prognostic factors typical of presenting myeloma.

Risk stratification in myeloma requires an accurate assessment of the presence of a range of molecular abnormalities including the differing IGH translocations and the recurrent copy number abnormalities that can impact clinical behavior. Currently, interphase fluorescence in situ hybridization is used to detect these abnormalities. High failure rates, slow turnaround, cost, and labor intensiveness make it difficult and expensive to use in routine clinical practice. Multiplex ligation-dependent probe amplification (MLPA), a molecular approach based on a multiplex polymerase chain reaction method, offers an alternative for the assessment of copy number changes present in the myeloma genome. Here, we provide evidence showing that MLPA is a powerful tool for the efficient detection of copy number abnormalities and when combined with expression assays, MLPA can detect all of the prognostically relevant molecular events which characterize presenting myeloma. This approach opens the way for a molecular diagnostic strategy that is efficient, high throughput, and cost effective.

Global methylation analysis identifies prognostically important epigenetically inactivated tumor suppressor genes in multiple myeloma.

Outcome in multiple myeloma is highly variable and a better understanding of the factors that influence disease biology is essential to understand and predict behavior in individual patients. In the present study, we analyzed combined genomewide DNA methylation and gene expression data of patients treated in the Medical Research Council Myeloma IX trial. We used these data to identify epigenetically repressed tumor suppressor genes with prognostic relevance in myeloma. We identified 195 genes with changes in methylation status that were significantly associated with prognosis. Combining DNA methylation and gene expression data led to the identification of the epigenetically regulated tumor modulating genes GPX3, RBP1, SPARC, and TGFBI. Hypermethylation of these genes was associated with significantly shorter overall survival, independent of age, International Staging System score, and adverse cytogenetics. The 4 differentially methylated and expressed genes are known to mediate important tumor suppressive functions including response to chemotherapy (TGFBI), interaction with the microenvironment (SPARC), retinoic acid signaling (RBP1), and the response to oxidative stress (GPX3), which could explain the prognostic impact of their differential methylation. Assessment of the DNA methylation status of the identified genes could contribute to the molecular characterization of myeloma, which is prerequisite for an individualized treatment approach.

Characterization of IGH locus breakpoints in multiple myeloma indicates a subset of translocations appear to occur in pregerminal center B cells.

Translocations in myeloma are thought to occur solely in mature B cells in the germinal center through class switch recombination (CSR). We used a targeted captured technique followed by massively parallel sequencing to determine the exact breakpoints in both the immunoglobulin heavy chain (IGH) locus and the partner chromosome in 61 presentation multiple myeloma samples. The majority of samples (62%) have a breakpoint within the switch regions upstream of the IGH constant genes and are generated through CSR in a mature B cell. However, the proportion of CSR translocations is not consistent between cytogenetic subgroups. We find that 100% of t(4;14) are CSR-mediated; however, 21% of t(11;14) and 25% of t(14;20) are generated through DH-JH recombination activation gene-mediated mechanisms, indicating they occur earlier in B-cell development at the pro-B-cell stage in the bone marrow. These 2 groups also generate translocations through receptor revision, as determined by the breakpoints and mutation status of the segments used in 10% and 50% of t(11;14) and t(14;20) samples, respectively. The study indicates that in a significant number of cases the translocation-based etiological events underlying myeloma may arise at the pro-B-cell hematological progenitor cell level, much earlier in B-cell development than was previously thought.

Intraclonal heterogeneity and distinct molecular mechanisms characterize the development of t(4;14) and t(11;14) myeloma.

We have used whole exome sequencing to compare a group of presentation t(4;14) with t(11;14) cases of myeloma to define the mutational landscape. Each case was characterized by a median of 24.5 exonic nonsynonymous single-nucleotide variations, and there was a consistently higher number of mutations in the t(4;14) group, but this number did not reach statistical significance. We show that the transition and transversion rates in the 2 subgroups are similar, suggesting that there was no specific mechanism leading to mutation differentiating the 2 groups. Only 3% of mutations were seen in both groups, and recurrently mutated genes include NRAS, KRAS, BRAF, and DIS3 as well as DNAH5, a member of the axonemal dynein family. The pattern of mutation in each group was distinct, with the t(4;14) group being characterized by deregulation of chromatin organization, actin filament, and microfilament movement. Recurrent RAS pathway mutations identified subclonal heterogeneity at a mutational level in both groups, with mutations being present as either dominant or minor subclones. The presence of subclonal diversity was confirmed at a single-cell level using other tumor-acquired mutations. These results are consistent with a distinct molecular pathogenesis underlying each subgroup and have important impacts on targeted treatment strategies. The Medical Research Council Myeloma IX trial is registered under ISRCTN68454111.

Identification of a novel t(7;14) translocation in multiple myeloma resulting in overexpression of EGFR.

IGH translocations in myeloma are a primary event and determine the prognostic outcome of a patient. These events are characterized by FISH and classical cytogenetics, but in a small proportion of samples a translocation involving the IGH locus can be detected but the partner chromosome cannot be identified. These cases are usually genetically complex and are the result of cryptic events that cannot be discerned at the resolution of FISH. Here we analyzed a sample with an unidentified translocation partner using a targeted capture and massively parallel sequencing. We identified the partner chromosome as a t(7;14) with the breakpoint upstream of EGFR. This sample over-expresses the target oncogene, EGFR. This case represents a rare and novel translocation in myeloma, from which a targeted personalized treatment, in the form of EGFR inhibitors, which are commonly used in other cancer types, could be used.

Long-read sequencing for brain tumors.

Brain tumors and genomics have a long-standing history given that glioblastoma was the first cancer studied by the cancer genome atlas. The numerous and continuous advances through the decades in sequencing technologies have aided in the advanced molecular characterization of brain tumors for diagnosis, prognosis, and treatment. Since the implementation of molecular biomarkers by the WHO CNS in 2016, the genomics of brain tumors has been integrated into diagnostic criteria. Long-read sequencing, also known as third generation sequencing, is an emerging technique that allows for the sequencing of longer DNA segments leading to improved detection of structural variants and epigenetics. These capabilities are opening a way for better characterization of brain tumors. Here, we present a comprehensive summary of the state of the art of third-generation sequencing in the application for brain tumor diagnosis, prognosis, and treatment. We discuss the advantages and potential new implementations of long-read sequencing into clinical paradigms for neuro-oncology patients.

Aberrant global methylation patterns affect the molecular pathogenesis and prognosis of multiple myeloma.

We used genome-wide methylation microarrays to analyze differences in CpG methylation patterns in cells relevant to the pathogenesis of myeloma plasma cells (B cells, normal plasma cells, monoclonal gammopathy of undetermined significance [MGUS], presentation myeloma, and plasma cell leukemia). We show that methylation patterns in these cell types are capable of distinguishing nonmalignant from malignant cells and the main reason for this difference is hypomethylation of the genome at the transition from MGUS to presentation myeloma. In addition, gene-specific hypermethylation was evident at the myeloma stage. Differential methylation was also evident at the transition from myeloma to plasma cell leukemia with remethylation of the genome, particularly of genes involved in cell-cell signaling and cell adhesion, which may contribute to independence from the bone marrow microenvironment. There was a high degree of methylation variability within presentation myeloma samples, which was associated with cytogenetic differences between samples. More specifically, we found methylation subgroups were defined by translocations and hyperdiploidy, with t(4;14) myeloma having the greatest impact on DNA methylation. Two groups of hyperdiploid samples were identified, on the basis of unsupervised clustering, which had an impact on overall survival. Overall, DNA methylation changes significantly during disease progression and between cytogenetic subgroups.

Intracranial hematolymphoid malignancies: A case series with molecular characterization.

OBJECTIVE: Central nervous system (CNS) manifestations of hematologic malignancies are uncommon and often have a poor prognosis. As hematologic neoplasms are typically chemotherapy- and radiotherapy-sensitive, surgical resection is usually not indicated; thus, opportunities for in-depth characterization of CNS hematologic tumors are limited. Here, we report four cases of rare intracranial hematologic tumors requiring surgical intervention, allowing for histopathologic and genomic characterization. METHODS: The clinical course, genetic perturbations, and histopathological features are described for a case of 1) primary marginal zone B-cell lymphoma of the dura as well as cases of brain metastases of 2) cutaneous T-cell lymphoma, 3) acute myeloid leukemia/myeloid sarcoma, and 4) multiple myeloma. Targeted DNA sequencing, fluorescence in situ hybridization, cytogenetic analysis, flow cytometry and immunohistochemical staining were used to assess the lesions. RESULT: Molecular and histopathological characterizations of four unusual presentations of hematolymphoid diseases involving the CNS are presented. Genetic abnormalities were identified in each lesion, including chromosomal aberrations and single nucleotide variants resulting in missense or nonsense mutations in oncogenes. CONCLUSIONS: Our case series provides insight into unique pathological phenotypes of hematologic neoplasms with atypical CNS involvement. We offer targets for future studies by identifying potentially pathogenic genetic variants in these lesions, as the full implications of the novel molecular abnormalities described remain unclear.

The clinical impact and molecular biology of del(17p) in multiple myeloma treated with conventional or thalidomide-based therapy.

Hemizygous deletion of 17p (del(17p)) has been identified as a variable associated with poor prognosis in myeloma, although its impact in the context of thalidomide therapy is not well described. The clinical outcome of 85 myeloma patients with del(17p) treated in a clinical trial incorporating both conventional and thalidomide-based induction therapies was examined. The clinical impact of deletion, low expression, and mutation of TP53 was also determined. Patients with del(17p) did not have inferior response rates compared to patients without del(17p), but, despite this, del(17p) was associated with impaired overall survival (OS) (median OS 26.6 vs. 48.5 months, P < 0.001). Within the del(17p) group, thalidomide induction therapy was associated with improved response rates compared to conventional therapy, but there was no impact on OS. Thalidomide maintenance was associated with impaired OS, although our analysis suggests that this effect may have been due to confounding variables. A minimally deleted region on 17p13.1 involving 17 genes was identified, of which only TP53 and SAT2 were underexpressed. TP53 was mutated in <1% in patients without del(17p) and in 27% of patients with del(17p). The higher TP53 mutation rate in samples with del(17p) suggests a role for TP53 in these clinical outcomes. In conclusion, del(17p) defined a patient group associated with short survival in myeloma, and although thalidomide induction therapy was associated with improved response rates, it did not impact OS, suggesting that alternative therapeutic strategies are required for this group.