BCOR and BCORL1 mutations in pediatric acute myeloid leukemia

FUS-ERG induces late-onset azacitidine resistance in acute myeloid leukaemia cells

FUS-ERG is a chimeric gene with a poor prognosis, found in myelodysplastic syndromes (MDS) and acute myeloid leukaemia (AML). It remains unclear whether DNA hypomethylating agents, including azacitidine (Aza), are effective in FUS-ERG-harbouring AML and how FUS-ERG induces chemoresistance. Stable Ba/F3 transfectants with FUS-ERG were repeatedly exposed to Aza for 7 days of treatment and at 21-day intervals to investigate Aza sensitivity. Stable FUS-ERG transfectants acquired resistance acquired resistance after three courses of Aza exposure. RNA sequencing (RNA-seq) was performed when Aza susceptibility began to change; genes with altered expression or transcript variants were identified. Molecular signatures of these genes were analysed using gene ontology. RNA-seq analyses identified 74 upregulated and 320 downregulated genes involved in cell motility, cytokine production, and kinase activity. Additionally, 1321 genes with altered transcript variants were identified, revealing their involvement in chromatin organisation. In a clinical case of AML with FUS-ERG, we compared whole-genome alterations between the initial MDS diagnosis and AML recurrence after Aza treatment. Genes with non-synonymous or near mutations in transcription regulatory areas (TRAs), additionally detected in AML recurrence, were collated with the gene list from RNA-seq to identify genes involved in acquiring Aza resistance in the presence of FUS-ERG. Whole-genome sequencing of clinical specimens identified 29 genes with non-synonymous mutations, including BCOR, and 48 genes located within 20 kb of 54 TRA mutations in AML recurrence. These genes were involved in chromatin organisation and included NCOR2 as an overlapping gene with RNA-seq data. Transcription regulators involved in mutated TRAs were skewed and included RCOR1 in AML recurrence. We tested the efficacy of BH3 mimetics, including venetoclax and S63845, in primary Aza-resistant AML cells treated with FUS-ERG. Primary FUS-ERG-harbouring AML cells acquiring Aza resistance affected the myeloid cell leukaemia-1 (MCL1) inhibitor S63845 but not while using venetoclax, despite no mutations in BCL2. FUS-ERG promoted Aza resistance after several treatments. The disturbance of chromatin organisation might induce this by co-repressors, including BCOR, NCOR2, and RCOR1. MCL1 inhibition could partially overcome Aza resistance in FUS-ERG-harbouring AML cells.

Next-Generation Sequencing-Based Genomic Profiling of Children with Acute Myeloid Leukemia

Pediatric acute myeloid leukemia (AML) represents a major cause of childhood leukemic mortality, with only a limited number of studies investigating the molecular landscape of the disease. Here, we present an integrative analysis of cytogenetic and molecular profiles of 75 patients with pediatric AML from a multicentric, real-world patient cohort treated according to AML Berlin-Frankfurt-Münster protocols. Targeted next-generation sequencing of 54 genes revealed 17 genes that were recurrently mutated in >5% of patients. Considerable differences were observed in the mutational profiles compared with previous studies, as BCORL1, CUX1, KDM6A, PHF6, and STAG2 mutations were detected at a higher frequency than previously reported, whereas KIT, NRAS, and KRAS were less frequently mutated. Our study identified novel recurrent mutations at diagnosis in the BCORL1 gene in 9% of the patients. Tumor suppressor gene (PHF6, TP53, and WT1) mutations were found to be associated with induction failure and shorter event-free survival, suggesting important roles of these alterations in resistance to therapy and disease progression. Comparison of the mutational landscape at diagnosis and relapse revealed an enrichment of mutations in tumor suppressor genes (16.2% versus 44.4%) and transcription factors (35.1% versus 55.6%) at relapse. Our findings shed further light on the heterogeneity of pediatric AML and identify previously unappreciated alterations that may lead to improved molecular characterization and risk stratification of pediatric AML.

The MLL-Menin Interaction is a Therapeutic Vulnerability in NUP98-rearranged AML

Chromosomal translocations involving the NUP98 locus are among the most prevalent rearrangements in pediatric acute myeloid leukemia (AML). AML with NUP98 fusions is characterized by high expression of HOXA and MEIS1 genes and is associated with poor clinical outcome. NUP98 fusion proteins are recruited to their target genes by the mixed lineage leukemia (MLL) complex, which involves a direct interaction between MLL and Menin. Here, we show that therapeutic targeting of the Menin-MLL interaction inhibits the propagation of NUP98-rearrranged AML both ex vivo and in vivo. Treatment of primary AML cells with the Menin inhibitor revumenib (SNDX-5613) impairs proliferation and clonogenicity ex vivo in long-term coculture and drives myeloid differentiation. These phenotypic effects are associated with global gene expression changes in primary AML samples that involve the downregulation of many critical NUP98 fusion protein-target genes, such as MEIS1 and CDK6. In addition, Menin inhibition reduces the expression of both wild-type FLT3 and mutated FLT3-ITD, and in combination with FLT3 inhibitor, suppresses patient-derived NUP98-r AML cells in a synergistic manner. Revumenib treatment blocks leukemic engraftment and prevents leukemia-associated death of immunodeficient mice transplanted with NUP98::NSD1 FLT3-ITD-positive patient-derived AML cells. These results demonstrate that NUP98-rearranged AMLs are highly susceptible to inhibition of the MLL-Menin interaction and suggest the inclusion of AML patients harboring NUP98 fusions into the clinical evaluation of Menin inhibitors.

Molecular evaluation of gene mutation profiles and copy number variations in pediatric acute myeloid leukemia

BACKGROUND

The objectives of this study were to investigate the mutation profiles of targeted genes and copy number variations (CNVs) in normal cytogenetics (CN) pediatric acute myeloid leukemia (AML).

METHODS

This prospective study was conducted from October 2018 to December 2020. The next-generation sequencing (NGS) and chromosomal microarray analyses (CMA) were performed in pediatric CN-AML patients.

RESULTS

Out of 94 children (aged ≤18 years), 70 patients with AML (24 excluded) underwent conventional karyotyping/cytogenetic analyses. Forty-five (64.3%) of patients had abnormal/ recurrent cytogenetic abnormalities and 25 (35.7%) had normal cytogenetics. Twenty-three out of 25 CN-AML were further processed for gene mutation profile and CNVs using NGS and CMA, respectively. Twenty-two out of 23 (95.7%) patients were detected to have mutations in various genes. The common mutations were: NRAS, NPM1, CEBPA, KRAS, KIT, RUNX1, NOTCH1, WT1, GATA1, GATA2, FLT3, KMT2D, FLT3-TKD, and PHF6. Copy number variations (CNVs) were detected in nine patients (39%), and eight (34.8%) had a long contiguous stretch of homozygosity (LCSH) /loss of heterozygosity (LOH). An LCSH was detected on chromosomes 5, 7, 11, and 19. The gains were more common than losses (8 vs 2). The gains were observed on chromosomes 8, 9, 14, 19, 21, and 22, and the losses were detected on chromosomes 7 and 10. Monosomy was observed in three patients. Three patients (monosomy7, n = 2, and FLT-ITD, n = 1) were reclassified into the high-risk category. Post-induction, complete remission was achieved in all evaluable patients.

CONCLUSION

CN-AML patients have genetic abnormalities that can be detected by more advanced techniques like NGS and CMA. These genetic abnormalities play a role in risk stratification that may remain hidden in otherwise CN-AML.

BCORL1 S878G, GNB1 G116S, SH2B3 A536T, and KMT2D S3708R tetramutation co-contribute to a pediatric acute myeloid leukemia: Case report and literature review

Acute myeloid leukemia (AML) is a clinically, morphologically, and genetically heterogeneous group of malignancies characterized by a wide range of genomic alterations responsible for defective regulation of the differentiation and self-renewal programs of hematopoietic stem cells. Here, we report a 4-month-old boy who had acute onset with leukocytosis and abdominal mass. The morphological analysis of bone marrow (BM) smear revealed extremely marrow hyperplasia, large quantities of immature cells, and primary and immature monocytic hyperplasia accounting for 57.5% of nucleated cells. The chromosome karyotype of the case was complex, representing 48, XY, +13, +19[12]/48, idem, del (p12)[8]. After RNAs sequencing, a mutation (c.346G > A, p.G116S) of the GNB1 gene was detected and localized to the mutational hotspot in Exon 7. Meanwhile, the other three mutations were identified by next-generation sequencing (NGS) and whole-exome sequencing (WES) of DNA from the BM aspirate and oral swab, including BCORL1 mutation [c.2632A > G, p.S878G, mutation allele frequency (VAF): 99.95%], SH2B3 mutation (c.1606G > A, p.A536T, VAF: 51.17%), and KMT2D mutation (c.11124C > G, p.S3708R, VAF: 48.95%). BCORL1 mutations have been associated with the pathogenesis of AML, whereas other mutations have rarely been previously reported in pediatric AML. The patient did not undergo the combination chemotherapy and eventually died of respiratory failure. In conclusion, the concurrence of BCORL1, GNB1, SH2B3, and KMT2D mutations may be a mutationally detrimental combination and contribute to disease progression.

Murine Models of Acute Myeloid Leukemia

Acute myeloid leukemia (AML) is a phenotypically and genetically heterogeneous hematologic malignancy. Extensive sequencing efforts have mapped the genomic landscape of adult and pediatric AML revealing a number of biologically and prognostically relevant driver lesions. Beyond identifying recurrent genetic aberrations, it is of critical importance to fully delineate the complex mechanisms by which they contribute to the initiation and evolution of disease to ultimately facilitate the development of targeted therapies. Towards these aims, murine models of AML are indispensable research tools. The rapid evolution of genetic engineering techniques over the past 20 years has greatly advanced the use of murine models to mirror specific genetic subtypes of human AML, define cell-intrinsic and extrinsic disease mechanisms, study the interaction between co-occurring genetic lesions, and test novel therapeutic approaches. This review summarizes the mouse model systems that have been developed to recapitulate the most common genomic subtypes of AML. We will discuss the strengths and weaknesses of varying modeling strategies, highlight major discoveries emanating from these model systems, and outline future opportunities to leverage emerging technologies for mechanistic and preclinical investigations.

A Case Series of BCOR Sarcomas With a New Splice Variant of BCOR/CCNB3 Fusion Gene

BACKGROUND/AIM

Undifferentiated round cell sarcomas are a heterogeneous group of sarcomas. Identification of BCOR alterations, such as BCOR/CCNB3 and BCOR/MAML3 fusion genes and BCOR ITD has recently contributed in the precise diagnosis of these neoplasms, defining a new entity of the current classification of soft tissue and bone sarcomas. BCOR sarcomas share both morphological and genetic characteristics distinct from Ewing sarcomas. The scope of our study was to retrospectively identify BCOR sarcomas and find the correlations with the clinical outcome of these patients.

PATIENTS AND METHODS

Histopathology and immunohistochemistry of pediatric tumor samples were combined with molecular testing (PCR) and fluorescent in situ hybridization to find BCOR sarcomas.

RESULTS

We, herein, present our experience with BCOR sarcomas in a referral center of Greece. Moreover, we report in one case the detection of a variant BCOR/CCNB3 fusion not previously described.

CONCLUSION

We are the first to report a splice variant of BCOR/CCNB3 which reveals the central position of BCOR in the oncogenesis of these tumors, furthermore we highlight the importance of molecular diagnostics in Ewing-like sarcomas and discuss the current treatment options for this rare entity.

BCOR gene alterations in hematological diseases

The BCL6 co-repressor (BCOR) is a transcription factor involved in the control of embryogenesis, mesenchymal stem cells function, hematopoiesis and lymphoid development. Recurrent somatic clonal mutations of the BCOR gene and its homologue BCORL1 have been detected in several hematological malignancies and aplastic anemia. They are scattered across the whole gene length and mostly represent frameshifts (deletions, insertions), nonsense and missence mutations. These disruptive events lead to the loss of full-length BCOR protein and to the lack or low expression of a truncated form of the protein, both consistent with the tumor suppressor role of BCOR. BCOR and BCORL1 mutations are similar to those causing two rare X-linked diseases: the oculo-facio-cardio-dental (OFCD) and the Shukla-Vernon syndromes, respectively. Here, we focus on the structure and function of normal BCOR and BCORL1 in normal hematopoietic and lymphoid tissues and review the frequency and clinical significance of the mutations of these genes in malignant and non-malignant hematological diseases. Moreover, we discuss the importance of mouse models to better understand the role of Bcor loss, alone and combined with alterations of other genes (e.g. Dnmt3a and Tet2), in promoting hematological malignancies and in providing a useful platform for the development of new targeted therapies.

Loss-of-Function Mutations of BCOR Are an Independent Marker of Adverse Outcomes in Intensively Treated Patients with Acute Myeloid Leukemia

Acute myeloid leukemia (AML) is characterized by recurrent genetic events. The BCL6 corepressor (BCOR) and its homolog, the BCL6 corepressor-like 1 (BCORL1), have been reported to be rare but recurrent mutations in AML. Previously, smaller studies have reported conflicting results regarding impacts on outcomes. Here, we retrospectively analyzed a large cohort of 1529 patients with newly diagnosed and intensively treated AML. BCOR and BCORL1 mutations were found in 71 (4.6%) and 53 patients (3.5%), respectively. Frequently co-mutated genes were DNTM3A, TET2 and RUNX1. Mutated BCORL1 and loss-of-function mutations of BCOR were significantly more common in the ELN2017 intermediate-risk group. Patients harboring loss-of-function mutations of BCOR had a significantly reduced median event-free survival (HR = 1.464 (95%-Confidence Interval (CI): 1.005-2.134), p = 0.047), relapse-free survival (HR = 1.904 (95%-CI: 1.163-3.117), p = 0.01), and trend for reduced overall survival (HR = 1.495 (95%-CI: 0.990-2.258), p = 0.056) in multivariable analysis. Our study establishes a novel role for loss-of-function mutations of BCOR regarding risk stratification in AML, which may influence treatment allocation.

Bcor deficiency perturbs erythro-megakaryopoiesis and cooperates with Dnmt3a loss in acute erythroid leukemia onset in mice

Recurrent loss-of-function mutations of BCL6 co-repressor (BCOR) gene are found in about 4% of AML patients with normal karyotype and are associated with DNMT3a mutations and poor prognosis. Therefore, new anti-leukemia treatments and mouse models are needed for this combinatorial AML genotype. For this purpose, we first generated a Bcor^-/- knockout mouse model characterized by impaired erythroid development (macrocytosis and anemia) and enhanced thrombopoiesis, which are both features of myelodysplasia/myeloproliferative neoplasms. We then created and characterized double Bcor^-/-/Dnmt3a^-/- knockout mice. Interestingly, these animals developed a fully penetrant acute erythroid leukemia (AEL) characterized by leukocytosis secondary to the expansion of blasts expressing c-Kit+ and the erythroid marker Ter119, macrocytic anemia and progressive reduction of the thrombocytosis associated with loss of Bcor alone. Transcriptomic analysis of double knockout bone marrow progenitors revealed that aberrant erythroid skewing was induced by epigenetic changes affecting specific transcriptional factors (GATA1-2) and cell-cycle regulators (Mdm2, Tp53). These findings prompted us to investigate the efficacy of demethylating agents in AEL, with significant impact on progressive leukemic burden and mice overall survival. Information gained from our model expands the knowledge on the biology of AEL and may help designing new rational treatments for patients suffering from this high-risk leukemia.

BCOR Expression in Mullerian Adenosarcoma: A Potential Diagnostic Pitfall

Adenosarcoma can mimic high-grade endometrial stromal sarcoma with ZC3H7B-BCOR fusion that may show entrapped glands and often exhibits diffuse BCOR expression. We encountered diffuse BCOR expression in rare adenosarcomas and sought to define its frequency among a larger cohort of these tumors. BCOR immunohistochemistry was performed on archival formalin-fixed paraffin-embedded tumor tissue in 13 of 14 adenosarcomas with and without stromal overgrowth arising in the uterus or ovary. The staining intensity and percentage of positive tumor nuclei in the mesenchymal component were evaluated. Eleven cases with sufficient tumoral tissue were subjected to fluorescence in situ hybridization for the detection of BCOR, BCORL1, NUTM1, ZC3H7B, and JAZF1 rearrangement. Three cases were subjected to targeted RNA sequencing. BCOR was expressed in 9 of 13 (70%) tumors, including 6 with and 3 without stromal overgrowth. Moderate to strong staining in >70% of cells was seen throughout in 1 low-grade and 6 high-grade tumors, 5 of which had stromal overgrowth. No staining was seen in 3 low-grade and 1 high-grade tumors with stromal overgrowth. One tumor demonstrating extensive sex cord-like differentiation and diffuse BCOR expression harbored JAZF1 and BCORL1 rearrangements. No BCOR or BCORL1 rearrangement was identified in the remaining tumors. BCOR expression is seen in most adenosarcomas with and without stromal overgrowth. BCORL1 rearrangement is seen in rare tumors with diffuse BCOR expression. Assessment of BCOR or BCORL1 rearrangement status is required in adenosarcomas demonstrating BCOR expression.

Case of aggressive metastatic follicular variant papillary thyroid carcinoma with BRAF K601E and BCORL1 mutations

BCL6 corepressor like-1 (BCORL1) mutation has rarely been described in thyroid cancer or in association with BRAF mutations in any malignancy. However, we report a 49-year-old woman who had aggressive follicular variant papillary thyroid carcinoma (FV-PTC) with both the BRAF K601E and BCORL1 mutations. The patient underwent a total thyroidectomy for a 3.6 cm right thyroid nodule and a smaller lesion in the left lobe in 2007; both were FV-PTCs with no lymphovascular invasion or metastases. In 2015, a positron emission tomography-CT scan showed a small defect in the left posterior lateral fifth rib with mild increased hypermetabolic activity with standardised uptake value of 3.9 and another lesion in the right hip at the junction of the femoral neck and trochanter. Tumour biopsy and genetic analysis revealed an uncommon BRAF K601E and a rare BCORL1 mutation. While rare, we report a case of aggressive FV-PTC with both the BRAF K601E and BCORL1 mutations.

BCOR involvement in cancer

BCOR is a gene that encodes for an epigenetic regulator involved in the specification of cell differentiation and body structure development and takes part in the noncanonical polycomb repressive complex 1. This review provides a comprehensive summary of BCOR’s involvement in oncology, illustrating that various BCOR aberrations, such as the internal tandem duplications of the PCGF Ub-like fold discriminator domain and different gene fusions (mainly BCOR–CCNB3, BCOR–MAML3 and ZC3H7B–BCOR), represent driver elements of various sarcomas such as clear cell sarcoma of the kidney, primitive mesenchymal myxoid tumor of infancy, small round blue cell sarcoma, endometrial stromal sarcoma and histologically heterogeneous CNS neoplasms group with similar genomic methylation patterns known as CNS-HGNET-BCOR. Furthermore, other BCOR alterations (often loss of function mutations) recur in a large variety of mesenchymal, epithelial, neural and hematological tumors, suggesting a central role in cancer evolution.

Molecular characteristic of acute leukemias with t(16;21)/FUS-ERG

T(16;21)(p11;q22)/FUS-ERG is a rare but recurrent translocation in acute leukemias and in some types of solid tumors. Due to multiple types of FUS-ERG transcripts, PCR-based minimal residual disease detection is impeded. In this study, we evaluated a cohort of pediatric patients with t(16;21)(p11;q22)/FUS-ERG and revealed fusion gene breakpoints. We implemented next-generation sequencing (NGS) on long PCR amplicons for the detection of fusion genes with unknown partners or DNA breakpoints. That allowed us to describe different fusion variants of FUS/ERG in different patients and to detect MRD on both RNA and DNA levels. We also found several accompanying mutations in epigenetic regulators (DNMT3A, ASXL1, BCOR) by targeted NGS approach in AML cases. These mutations preceded full transformation by t(16;21)(p11;q22)/FUS-ERG and allowed us to trace clonal evolution on all steps of therapy. As a casual observation, the ASXL1 mutation was found in the unrelated donor hematopoietic cells.

Pediatric Anaplastic Embryonal Rhabdomyosarcoma: Targeted Therapy Guided by Genetic Analysis and a Patient-Derived Xenograft Study

Therapy for rhabdomyosarcoma (RMS) has generally been limited to combinations of conventional cytotoxic agents similar to regimens originally developed in the late 1960s. Recently, identification of molecular alterations through next-generation sequencing of individual tumor specimens has facilitated the use of more targeted therapeutic approaches for various malignancies. Such targeted therapies have revolutionized treatment for some cancer types. However, malignancies common in children, thus far, have been less amenable to such targeted therapies. This report describes the clinical course of an 8-year-old female with embryonal RMS having anaplastic features. This patient experienced multiple relapses after receiving various established and experimental therapies. Genomic testing of this RMS subtype revealed mutations in BCOR, ARID1A, and SETD2 genes, each of which contributes to epigenetic regulation and interacts with or modifies the activity of histone deacetylases (HDAC). Based on these findings, the patient was treated with the HDAC inhibitor vorinostat as a single agent. The tumor responded transiently followed by subsequent disease progression. We also examined the efficacy of vorinostat in a patient-derived xenograft (PDX) model developed using tumor tissue obtained from the patient’s most recent tumor resection. The antitumor activity of vorinostat observed with the PDX model reflected clinical observations in that obvious areas of tumor necrosis were evident following exposure to vorinostat. Histologic sections of tumors harvested from PDX tumor-bearing mice treated with vorinostat demonstrated induction of necrosis by this agent. We propose that the evaluation of clinical efficacy in this type of preclinical model merits further evaluation to determine if PDX models predict tumor sensitivity to specific agents and/or combination therapies.

Whole-exome sequencing reveals the spectrum of gene mutations and the clonal evolution patterns in paediatric acute myeloid leukaemia

Acute myeloid leukaemia (AML) is a molecularly and clinically heterogeneous disease. Targeted sequencing efforts have identified several mutations with diagnostic and prognostic values in KIT, NPM1, CEBPA and FLT3 in both adult and paediatric AML. In addition, massively parallel sequencing enabled the discovery of recurrent mutations (i.e. IDH1/2 and DNMT3A) in adult AML. In this study, whole-exome sequencing (WES) of 22 paediatric AML patients revealed mutations in components of the cohesin complex (RAD21 and SMC3), BCORL1 and ASXL2 in addition to previously known gene mutations. We also revealed intratumoural heterogeneities in many patients, implicating multiple clonal evolution events in the development of AML. Furthermore, targeted deep sequencing in 182 paediatric AML patients identified three major categories of recurrently mutated genes: cohesion complex genes [STAG2, RAD21 and SMC3 in 17 patients (8·3%)], epigenetic regulators [ASXL1/ASXL2 in 17 patients (8·3%), BCOR/BCORL1 in 7 patients (3·4%)] and signalling molecules. We also performed WES in four patients with relapsed AML. Relapsed AML evolved from one of the subclones at the initial phase and was accompanied by many additional mutations, including common driver mutations that were absent or existed only with lower allele frequency in the diagnostic samples, indicating a multistep process causing leukaemia recurrence.

BCORL1 is an independent prognostic marker and contributes to cell migration and invasion in human hepatocellular carcinoma

Background

The deregulation of E-cadherin has been considered as a leading cause of hepatocellular carcinoma (HCC) metastasis. BCL6 corepressor-like 1 (BCORL1) is a transcriptional corepressor and contributes to the repression of E-cadherin. However, the clinical significance of BCORL1 and its role in the metastasis of HCC remain unknown.

Methods

Differentially expressed BCORL1 between HCC and matched tumor-adjacent tissues, HCC cell lines and normal hepatic cell line were detected by Western blot. The expression of BCORL1 was altered by siRNAs or lentivirus-mediated vectors. Transwell assays were performed to determine HCC cell invasion and migration.

Results

Increased expression of BCORL1 protein was detected in HCC specimens and cell lines. Clinical association analysis showed that BCORL1 protein was expressed at significant higher levels in HCC patients with multiple tumor nodes, venous infiltration and advanced TNM tumor stage. Survival analysis indicated that high expression of BCORL1 protein conferred shorter overall survival (OS) and recurrence-free survival (RFS) of HCC patients. Multivariate Cox regression analysis disclosed that BCORL1 expression was an independent prognostic marker for predicting survival of HCC patients. Our in vitro studies demonstrated that BCORL1 prominently promoted HCC cell migration and invasion. Otherwise, an inverse correlation between BCORL1 and E-cadherin expression was observed in HCC tissues. BCORL1 inversely regulated E-cadherin abundance and subsequently facilitated epithelial-mesenchymal transition (EMT) in HCC cells. Notably, the effect of BCORL1 knockdown on HCC cells was abrogated by E-cadherin silencing.

Conclusions

BCORL1 may be a novel prognostic factor and promotes cell migration and invasion through E-cadherin repression-induced EMT in HCC.

Electronic supplementary material

The online version of this article (doi:10.1186/s12885-016-2154-z) contains supplementary material, which is available to authorized users.