Skip to main content

Identifying novel genetic alterations in pediatric acute lymphoblastic leukemia based on copy number analysis

Abstract

Copy number variations (CNVs) analysis may reveal molecular biomarkers and provide information on the pathogenesis of acute lymphoblastic leukemia (ALL). We investigated the gene copy number in childhood ALL by microarray and select three new recurrent CNVs to evaluate by real-time PCR assay: DMBT1, KIAA0125 and PRDM16 were selected due to high frequency of CNVs in ALL samples and based on their potential biological functions in carcinogenesis described in the literature. DBMT1 deletion was associated with patients with chromosomal translocations and is a potential tumor suppressor; KIAA0125 and PRDM16 may act as an oncogene despite having a paradoxical behavior in carcinogenesis. This study reinforces that microarrays/aCGH is it is a powerful tool for detection of genomic aberrations, which may be used in the risk stratification.

Introduction

Acute lymphoblastic leukemia (ALL) is the most common cancer in children [1]. Leukemia represents the ninth most common cancer in Brazil and the fifth most frequent in the north region [2]. Advances in cytogenetics and molecular cytogenetics has allowed the identification of genetic aberrations in more than 80% of ALL cases [3]. Establishing genetic background in ALL patients is important for the diagnosis, risk classification and therapeutic interventions [3]. However, some patients do not have an established chromosomal aberration, which complicates the risk classification.

Recent analysis has shown that copy number variations (CNVs) are common in ALL and leukemia in general, especially in genes involved in transcription, cell cycle regulation and B-cell differentiation, (e.g., CDKN2A/B, IKZF1, ETV6, EBF1, PAX5, BTG1 and PAR1) [4]. Additional CNVs could be helpful to refine ALL prognostic. The prognostic effect of CNVs depends on the other factors, such as the presence of additional molecular or cytogenetic aberrations; this situation reinforces the need to analyze these combined alterations [5].

The aim of this report is to assess and evaluate CNVs identified by aCGH from a cohort of Brazilian children with ALL. Three new recurrent CNVs were further evaluated by qPCR. We highlight that DMBT1, KIAA0125 and PRDM16 were chosen due to high frequency of aberrations in ALL samples and based on their biological functions as well the data present in the literature.

Methods

Patients

A total of 16 ALL pediatric patients (5 ± 3 years) treated at Octávio Lobo Children’s Cancer Hospital were selected for aCGH analysis. Additional 84 ALL pediatric samples were used as validation group in copy number qPCR assay. These patients were classified by immunophenotyping and morphology. Gene fusions were investigated by reverse transcription polymerase chain reaction (RTq-PCR) (Tables 1 and 2). The samples were collected before cancer treatment between 2017 and 2019.

Table 1 Characteristics of pediatric patients with ALL included in the study
Table 2 Nucleotide sequence of RTq-PCR primers

The age at diagnosis and white blood cell (WBC) count were the criteria for assigning prognostic risk of ALL, according to the National Cancer Institute (NCI): 1) high risk, WBC count greater than 50 × 109 cells/μL, age 1 year or less, or age 10 years or more; and 2) standard risk, WBC count 50 × 109 cells/μL or less, or between 1 and 10 years of age. The patients with BCR-ABL1 or MLL-AF4 also were assigned to the NCI high risk group. Written consent forms were obtained from all parents of patients. This study was approved by hospital ethics committee (CAAE: 00905812.1.0000.00.18).

Array comparative genomic hybridization

Genomic DNA was extracted from peripheral blood by Pure Link Genomic DNA Mini Kit (Invitrogen, California, USA). aCGH was performed using Agilent 4x180k CGH + SNP microarray (Santa Clara, USA). After DNA extraction, a restriction enzyme digestion step and labeling with fluorochrome cyanine 5 were performed using random primers and exo-Klenow fragment DNA polymerase. DNA control was labeled with fluorochrome cyanine 3. DNA samples from the patient and control were combined and hybridized on the microarray. Data were analyzed using the software Agilent’s CytoGenomics v5.0.

Real-time quantitative PCR

TaqMan Copy Number Assay (Applied Biosystems, California, USA) was used to assess copy number for DMBT1, KIAA0125 and PRDM16. Briefly, 1 μL of 10 ng DNA was added to 5 μL of TaqMan Universal Master Mix no UNG, with 0.5 μL of each probe and 3 μL of ultra-pure water. RNase P was used as a control. The amplification protocol consisted of: denaturation at 95 °C for 10 min, followed by 40 cycles of 95 °C for 15 s and 60 °C for 1 min. Relative quantification was determined using the 7500 Real-time PCR system and all samples were analyzed in quadruplicate. After amplification, we imported the experiment results containing threshold-cycle values for the copy number and reference assay into the Copy Caller Software v2.0 for post-PCR data analysis as previously described [6].

Statistical analysis

Fisher’s exact test was used to compare the distribution of aberrations between subgroups (high or standard risk; positive or negative for chromosomal translocation) and pathological features of the patients; Odds ratio (OR) with a 95% confidence interval (CI) were also calculated through the statistical program BioEstat® v5.0 [7]. p-values less than 0.05 were considered significant.

Results

aCGH profiling identifies recurrent alterations

aCGH date were available for all 16 cases, the average of copy number variations (CNVs) was 8.3 per sample. Gains were the most frequent event, the most frequently gained regions were on chromosomes 14 (q32.33) and 10 (q26.13), these regions include KIAA0125 (a lncRNA) and DMBT1 genes, respectively. Frequent losses were identified on chromosome 7 (7p12.2) and 9 (p21.3), which includes IKZF1 and CDKN2A/B genes, respectively.

Recurrent CNVs included gains of cytobands 14q (93.75%), 2p (68.75%), 17q (62.5%), 9q (56.25%), 10q (56.25%), 19q (56.25%), 22q (56.25%), 1p (50%), 7q (50%), 8p (50%) and 21q (50%); losses involving 7p (62.5%), 9p (56.25%), 15q (47.75%), 4q (37.5%) and 12q (31.25%). The list of recurrent CNVs found in at least two samples is provided in Table 3.

Table 3 The most frequent copy number variations found in pediatric ALL by aCGH

All patients have alteration in at least one of the main genes associated with ALL; ETV6, RUNX1, IKZF1, KMT2A (MLL) and BTG1 (Table 3.). The median of alterations in standard (SR) and high risk (HR) group were 56.6 (±15.4) and e 52.2 (±14.2), respectively. We confirmed the association of CDKN2A/B losses with positive cases for TCF3-PBX1 or BCR-ABL1 (p < 0.05). There was no statistically significant difference in the number of CNVs between patients with (CT+) or without (CT-) chromosomal translocation.

CNV evaluation by real-time qPCR

To validate aCGH results DMBT1, KIAA0125 and PRDM16 genes were analyzed by qPCR. Genes were chosen due to the high frequency of aberrations found in samples and based on their biological function (mainly transcriptional regulation) described in literature. It is noted that the CNV found in these genes are described here for the first time in leukemia, especially in ALL. The aberrations of the three selected genes identified from aCGH and qPCR were illustrated in Fig. 1.

Fig. 1
figure 1

Frequency of copy number variation of PRDM16, KIAA0125 and DMBT1 identified by aCGH and qPCR. aCGH for 16 samples; qPCR for 48 samples. AMP: amplification; DEL: deletion

The results of qPCR were compared between positive (CT+) or negative (CT-) for gene fusions subgroups. DMBT1 deletion was observed in 74% of patients (n = 62; 97.4% of CT+ and 53% of CT-); KIAA0125 amplification was detected in 59% of cases (n = 50; 95% of CT+ and 29% of CT-), these amplifications were more frequent in cases CT+; while PRDM16 was deleted in 50% of patients (n = 42; 87% CT- and 8% CT+), amplifications were observed in 42% of samples, which only correspond to CT+ cases (Fig. 2).

Fig. 2
figure 2

Copy number of PRDM16, KIAA0125 and DMBT1 in ALL samples identified by qPCR

Curiously, 50% (25/50) and 69% (24/35) of KIAA0125 and PRDM16 amplifications, respectively, were high-level amplifications (> 10 copies), however, classification of cases according to the level of amplification did not result in any significant association.

The frequency of aberrations in DMBT1, KIAA0125 and PRDM16 according to NCI risk group, gender, age and cytogenetic findings are show in Tables 4 and 5 Statistical analysis showed that DMBT1 deletion was more common in patients with > 1 to ≤10 years (OR = 3.38; 95% IC = 1.15–9.89) and more common in NCI-SR cases (OR = 0.198; 95% IC = 0.07–056). DMBT1 deletion also was associated with CT+ samples (OR = 33.2; 95% IC = 4.19–263.55) (Tables 4 and 5).

Table 4 Frequency of alterations according to characteristics of patients
Table 5 Frequency of alterations according to cytogenetic subgroups

KIAA0125 amplification was associated with CT+ cases (OR = 45.5; 95% IC = 9.54–217-16). PRDM16 gene deletion was associated with NCI-HR patients (OR = 91.4; 95% IC = 11.32–738.6) and CT- cases (OR = 0.01; 95% IC = 0.00–0.05) (Tables 4 and 5), while amplification was related to CT+ samples (p-value = < 0.001), data not shown in Tables 4 and 5).

Discussion

All patients analyzed by aCGH showed a heterogeneous copy-number pattern. We identified 133 CNVs, 18 them involved the most frequent changes already known or not yet related to ALL (Table 3). Unlike previous studies, here, amplifications were more frequent than deletions, possibly due the small sample number and the presence of hyperdiploid cases. On the other hand, similar to antecedent studies [4, 8, 9], the more frequently altered genes were related to cell cycle regulation (ETV6), tumor suppression (CDKN2A/B), apoptosis regulation (BTG1) and others (Table 3).

In agreement with the literature, in our study deletions of CDKN2A/B were associated with positive cases for TCF3-PBX1 or BCR-ABL1. CDKN2A/B are tumor suppressor genes acting in cell growth regulation and apoptosis [10]. The deletion of these genes are associated with poor prognosis, high white blood cell count and older age at diagnosis and BCR-ABL1 or TCF3-PBX1 translocations [11,12,13]; all characteristics found in our study group.

The aCGH study also identified for the first-time recurrent alterations of DMBT1, KIAA0125 and PRDM16 in ALL (Table 3). These genes were verified by qPCR in a larger sample number.

High amplification frequencies observed in aCGH was confirmed by qPCR just for KIAA0125. For the DMBT1 and PRDM16 deletions were prevalent in qPCR assays. This divergence is probably due to differences in sample size and by the presence of trisomy of chromosomes 1 and 10 in cases with copy number variation in PRDM16 and DMBT1, respectively. But new significant associations were observed for the three genes.

The high frequency of DMBT1 deletions observed here support aCGH results. DMBT1 encoding protein belongs to the scavenger receptor cysteine rich (SRCR) super family involved in mucosal immune defense, epithelial differentiation and tumor suppression [14, 15]. Many studies have showed that DMBT1 deletion or inactivation lead to tumorigenesis by regulating infiltration and metastasis of tumor cells [16]. Altered expression in certain stages of carcinogenesis was identified in different tumor types [17,18,19]. We found DMBT1 deletion associated with standard risk and CT+ cases. It is possible that DMBT1 deletion have a more specific function in development of ALL cases without a high risk chromosomal abnormality (which are mostly classified as standard risk), since only 14% of CT+ cases have high risk biomarker (BCR-ABL1 or MLL-AF4). Thus, DMBT1 loss collaborate as a secondary event in the progression of disease in CT+ patients, since it is know that chromosomal translocations are primary aberrations [13]. Although DMBT1 absence is considered a malignancy marker in many epithelial cancers, we reported for the first time DMBT1 deletion in ALL and we suggest that DMBT1 may be also involved in hematologic malignancies development.

LncRNAs are involved in gene expression at epigenetic, transcriptional and post-transcriptional level and are considered a strong promise as a biomarker and therapeutic target [20]. In this study, we found that KIAA0125 amplifications were more common in CT+ patients while in CT- cases, deletions were more prevalent. Recurrent KIAA0125 amplifications were statistically associated with CT+ cases. CNV or abnormal expression of KIAA0125 were observed in many tumor types [21,22,23,24,25,26]. Several recent studies in lncRNAs have shown that they have a critical role in different cancers acting as an oncogene or suppressor, in this sense, the role of KIAA0125 in carcinogenesis may be cell-type dependent [27]. In colon cancer development, KIAA0125 may contribute via the regulation of BCL2 expression by sponging hsa-miR-29b-3p or regulating PI3K-Akt signaling [28]. In addition, Forero-Castro et al. [4] identified losses on 14q32.33 (where KIAA0125 is located) related to overall survival of children ALL with leukocytosis. In 14q32 there are miRNA clusters that may influence the genes expression levels involved in lymphoid B-cell transformation and differentiation, suggesting that 14q32 losses could be used as a diagnostic marker [4, 29]. Hornung R. et al. (2018) have recently shown that KIAA0125 could play a mediating role in the influence RUNX1 gene fusions have on survival of LMA [30].

It is also presumable that KIAA0125 may act as miRNA sponges regulates mRNAs expression levels also in ALL, however the exact mechanism of action and possible target genes need to be further investigated. These findings along with our data leading to the assumption that KIAA0125 plays important role in development of leukemia and reinforce previous studies that suggested that lncRNAs may be utilized as diagnostic and prognostic markers in leukemia [20].

PRDM16 is characterized by the combination of a conserved N-terminal PR domain and a variable number of zinc fingers [31], it encodes a SMAD binding protein that may repress SMAD-mediated transcription, also functions as a modulator of TGF-beta signaling and exhibit methyltransferase activity [32, 33]. PRDM16 is involved in various biological processes including maintenance of brown adipocytes and hematopoiesis [34, 35]. Two main PRDM16 isoforms are the full-length and the PR-lacking generated by alternative splicing or alternative use of different promoters [36, 37]. Notably, PRDM proteins sometimes exert opposing effects on tumor development [38, 39].

In the present study most cases have PRMD16 deletions (50%), however in 90% of CT+ patients this gene is highly amplified (21 samples with > 10 copies) and significantly related to presence of gene fusions. Overexpression of PRDM16 in AML is associated with worse overall survival [39, 40] and is considered a risk factor for primary induction failure [41]. In addition it is associated with other gene fusions not investigated here [42]. Hu et al. [43] reported that PRDM16 transforming megakaryocyte-erythroid progenitors into myeloid leukemia stem cells. In another study, PRDM16 knockdown induced cell proliferation in rhabdoid tumor cells [44], suggesting that PRDM16 may be an oncogene in leukemia development, although in other tumor types PRDM16 has a controversial role [45, 46]. Thus, the role of PRDM16 in cancer biology has been poorly studied and remains to be fully elucidated.

A limitation of this study was the small sample size. However, this is one of the few studies from the northern region of Brazil with genomic analysis in leukemia. This region has a large territorial extension, which makes the diagnosis of cancer a challenge due to its financial viability and the difficult access to geographically isolated regions of cancer treatment centers [47].

In conclusion, this study reinforces that aCGH it is a powerful tool for to identify regions of copy number variations in childhood ALL patients and to identify new genes associated to leukemia. Through this technique, we identified recurrent alterations in genes DMBT1, KIAA0125 and PRDM16; these alterations were verified by qPCR and confirmed the possible involvement of these genes in the development of leukemia, especially in ALL. DMBT1 probably is also a tumor suppressor in leukemia and is associated with standard risk and cases with gene fusions. Although both have a paradoxical behavior in tumorigenesis our data indicates that KIAA0125 and PRDM16 may act as oncogene, once amplifications in these two genes were related to gene fusions and leukocytosis, respectively. The combination of two molecular cytogenetics techniques has identified three genes that may be targets for further biological analysis of acute lymphoblastic leukemia.

Availability of data and materials

The data will not be shared because some analyses of another genes study by aCGH are still being analyzed together with other data and have not yet been published.

Abbreviations

aCGH:

Array Comparative Genomic Hybridization

ALL:

Acute lymphoblastic leukemia

CNV:

Copy Number Variation

CT+:

Chromosome translocation positive

CT-:

Chromosome translocation negative

HR:

High Risk

NCI:

National Cancer Institute

qPCR:

quantitative Polymerase Chain Reaction

RT-PCR:

Reverse Transcription Polymerase Chain Reaction

SR:

Standard Risk

WBC:

White Blood Count

References

  1. Scheurer ME, Bondy ML, Gourney JG. Epidemiology of childhood cancer. In: Pizzo PA, Poplack DG, editors. Principles and practice of pediatric oncology. Pennsylvania: Williams and Wilkins; 2011. p. 2–16.

    Google Scholar 

  2. INCA - Instituto Nacional do Câncer. Estimativa 2020 - Incidência de câncer no Brasil. http://www.inca.gov.br. Acessed 17 fev 2020.

  3. Bacher U, Kohlmann A, Haferlach T. Gene expression profiling for diagnosis and therapy in acute leukaemia and other haematologic malignancies. Cancer Treat Rev. 2010. https://doi.org/10.1016/j.ctrv.2010.05.002.

  4. Forero-Castro M, Robledo C, Benito R, Abáigar M, África Martín A, Arefi M, et al. Genome-wide DNA copy number analysis of acute lymphoblastic leukemia identifies new genetic markers associated with clinical outcome. PLoS One. 2016. https://doi.org/10.1371/journal.pone.0148972.

  5. Hamadeh L, Enshaei A, Schwab C, Alonso CN, Attarbaschi A, Barbany G, et al. Validation of the United Kingdom copy-number alteration classifier in 3239 children with B-cell precursor ALL. Blood Adv. 2019. https://doi.org/10.1182/bloodadvances.2018025718.

  6. Graziano F, Galluccio N, Lorenzini P, Ruzzo A, Canestrari E, D’Emidio S, et al. Genetic activation ofthe MET pathway and prognosis of patients with high-risk, radicallyresected gastric cancer. J Clin Oncol. 2011. https://doi.org/10.1200/JCO.2011.36.7706.

  7. Ayres M, Ayres JR, Ayres DL, Santos AAS. BioEstat 5.0: statistical applications in the biological and medical sciences. Belém: CNPq; 2007.

    Google Scholar 

  8. Bungaro S, Dell'Orto MC, Zangrando A, Basso D, Gorletta T, Lo Nigro L, et al. Integration of genomic and gene expression data of childhood ALL without known aberrations identifies subgroups with specific genetic hallmarks. Genes Chromosomes Cancer. 2009. https://doi.org/10.1002/gcc.20616.

  9. Dirse V, Bertasiute A, Gineikiene E, Zvirblis T, Dambrauskiene R, Gerbutavicius R, et al. A population-based single nucleotide polymorphism array analysis of genomic aberrations in younger adult acute lymphoblastic leukemia patients. Genes Chromosomes Cancer. 2015. https://doi.org/10.1002/gcc.22246.

  10. Kim WY, Sharpless NE. The regulation of INK4/ARF in cancer and aging. Cell. 2006;127:265–75.

    Article  CAS  Google Scholar 

  11. Usvasalo A, Savola S, Räty R, Vettenrantaa K, Harila-Saarid A, Koistinene P, et al. CDKN2A deletions in acute lymphoblastic leukemia of adolescents and young adults: an array CGH study. Leuk Res. 2008. https://doi.org/10.1016/j.leukres.2008.01.014.

  12. Schwab CJ, Chilton L, Morrison H, Jones L, Al-Shehhi H, Erhorn A, et al. Genes commonly deleted in childhood B-cell precursor acute lymphoblastic leukemia: association with cytogenetics and clinical features. Haematologica. 2013. https://doi.org/10.3324/haematol.2013.085175.

  13. Moorman AV. Haematologica. New and emerging prognostic and predictive genetic biomarkers in B-cell precursor acute lymphoblastic leukemia. Haematologica. 2016. https://doi.org/10.3324/haematol.2015.141101.

  14. Mollenhauer J, Herbertz S, Holmskov U, Tolnay M, Krebs I, Merlo A, et al. DMBT1 encodes a protein involved in the immune defense and in epithelial differentiation and is highly unstable in cancer. Cancer Res. 2000;60:1704–10.

    CAS  PubMed  Google Scholar 

  15. Singh N, Gupta DK, Sharma S, Sahu DK, Mishra A, Yadav DK, et al. Single-nucleotide and copy-number variance related to severity of hypospadias. Pediatr Surg Int. 2018. https://doi.org/10.1007/s00383-018-4330-5.

  16. Takeshita H, Sato M, Shiwaku HO, Semba S, Sakurada A, Hoshi M, et al. Expression of the DMBT1 gene is frequently suppressed in human lung cancer. Jpn J Cancer Res. 1999;90(9):903–8.

    Article  CAS  Google Scholar 

  17. Cheng J, Greshoc J, Shi L, Zheng S, Menius A, Lee K. Good practice guidelines for biomarker discovery from array data: A case study for breast cancer prognosis. BMC Syst Bio. 2013. https://doi.org/10.1186/1752-0509-7-S4-S2.

  18. Muller H, Renner M, Helmke BM, Mollenhauer J, Felderhoff-Muser U. Elevated DMBT1 levels in neonatal gastrointestinal diseases. Histochem Cell Biol. 2016. https://doi.org/10.1007/s00418-015-1381-8.

  19. Muller H, Weiss C, Renner M, Felderhoff-Muser U, Mollenhauer J. DMBT1 promotes basal and meconium-induced nitric oxide production in human lung epithelial cells in vitro. HistochemCellBio. 2017. https://doi.org/10.1007/s00418-016-14939.

  20. Cruz-Miranda GM, Hidalgo-Miranda A, Bárcenas-López DA, Núñez-Enríquez JC, Ramírez-Bello J, Mejía-Aranguré JM, et al. Long non-coding RNA and acute leukemia. Int J Mol Sci. 2019. https://doi.org/10.3390/ijms20030735.

  21. Seabra AD, Araújo TM, Mello Junior FA, Di Felipe Ávila Alcântara D, De Barros AP, De Assumpção PP, et al. High-density array comparative genomic hybridization detects novel copy number alterations in gastric adenocarcinoma. Anticancer Res. 2014;34(11):6405–15.

    CAS  PubMed  Google Scholar 

  22. Chiu CG, Nakamura Y, Chong KK, Huang SK, Kawas NP, Triche T, et al. Genome-wide characterization of circulating tumor cells identifies novel prognostic genomic alterations in systemic melanoma metastasis. Clin Chem. 2014. https://doi.org/10.1373/clinchem.2013.213611.

  23. Singchat W, Hitakomate E, Rerkarmnuaychoke B, Suntronpong A, Fu B, Bodhisuwan W, et al. Genomic alteration in head and neck squamous cell carcinoma (HNSCC) cell lines inferred from karyotyping, molecular Cytogenetics, and Array comparative genomic hybridization. PLoS One. 2016. https://doi.org/10.1371/journal.pone.0160901.

  24. Li L, Qiang L, Yufeng L, Zhang L, Wu Y. Long non-coding RNA expression profiles predict clinical phenotypes of seminoma and yolk sac tumor. RSC Adv. 2017. https://doi.org/10.1039/c7ra12131h.

  25. Diniz MG, França JA, Vilas-Boas FAS, de Souza FTA, Calin GA, Gomez RS, et al. The long noncoding RNA KIAA0125 is upregulated in ameloblastomas. Pathol Res Pract. 2019. https://doi.org/10.1016/j.prp.2018.12.030.

  26. Yao K, Wang Q, Jia J, Zhao H. A competing endogenous RNA network identifies novel mRNA, miRNA and lncRNA markers for the prognosis of diabetic pancreatic cancer. Tumor Biol. 2017. https://doi.org/10.1177/1010428317707882.

  27. Yang Y, Zhao Y, Hu N, Zhao J, Bai Y. lncRNA KIAA0125 functions as a tumor suppressor modulating growth and metastasis of colorectal cancer via Wnt/β-catenin pathway. Cell Biol Int. 2019. https://doi.org/10.1002/cbin.11196.

  28. Yang Y, Zhao Y, Zhang W, Bai Y. Whole transcriptome sequencing identifies crucial genes associated with colon cancer and elucidation of their possible mechanisms of action. Onco Targets Ther. 2019. https://doi.org/10.2147/OTT.S195235.

  29. Agueli C, Cammarata G, Salemi D, Dagnino L, Nicoletti R, La Rosa M, et al. 14q32/miRNA clusters loss of heterozygosity in acute lymphoblastic leukemia is associated with up-regulation of BCL11a. Am J Hematol. 2010. https://doi.org/10.1002/ajh.21758.

  30. Hornung R, Jurinovic V, Batcha AMN, Bamopoulos SA, Rothenberg-Thurley M, Amler S, et al. Mediation analysis reveals common mechanisms of RUNX1 point mutations and RUNX1/RUNX1T1 fusions influencing survival of patients with acute myeloid leukemia. Sci Rep. 2018. https://doi.org/10.1038/s41598-018-29593-2.

  31. Di Zazzo E, De Rosa C, Abbondanza C, Moncharmont B. PRDM proteins: molecular mechanisms in signal transduction and transcriptional regulation. Biology (Basel). 2013. https://doi.org/10.3390/biology2010107.

  32. Warner DR, Greene RM, Pisano MM. PRDM16 in development and disease. Human Genet Embryol. 2014. https://doi.org/10.4172/2161-0436.1000121.

  33. Mehrian-Shai R, Yalon M, Moshe I, Barshack I, Nass D, Jacob J, et al. Identification of genomic aberrations in hemangioblastoma by droplet digital PCR and SNP microarray highlights novel candidate genes and pathways for pathogenesis. BMC Genomics. 2016. https://doi.org/10.1186/s12864-016-2370-6.

  34. Yin H, Pasut A, Soleimani VD, Bentzinger CF, Antoun G, Thorn S, et al. MicroRNA-133 controls brown adipose determination in skeletal muscle satellite cells by targeting Prdm16. Cell Metab. 2013. https://doi.org/10.1016/j.cmet.2013.01.004.

  35. Aguilo F, Avagyan S, Labar A. Prdm16 is a physiologic regulator of hematopoietic stem cells. Blood. 2011. https://doi.org/10.1182/blood-2010-08-300145.

  36. Sorrentino A, Federico A, Rienzo M, Gazzerro P, Bifulco M, Ciccodicola A. PR/SET domain family and Cancer: novel insights from the Cancer genome atlas. Int J Mol Sci. 2018. https://doi.org/10.3390/ijms19103250.

  37. Mzoughi S, Tan YX, Low D, Guccione E. The role of PRDMs in cancer: one family, two sides. Curr Opin Genet Dev. 2016. https://doi.org/10.1016/j.gde.2016.03.009.

  38. Fog CK, Galli GG, Lund AH. PRDM proteins: important players in differentiation and disease. BioEssays. 2012. https://doi.org/10.1002/bies.201100107.

  39. Yamato G, Yamaguchi H, Handa H, Shiba N, Kawamura M, Wakita S, et al. Clinical features and prognostic impact of PRDM16 expression in adult acute myeloid leukemia. Genes Chromosomes Cancer. 2017. https://doi.org/10.1002/gcc.22483.

  40. Duhoux FP, Ameye G, Montano-Almendras CP, Bahloula K, Mozziconacci MJ, Laibe S, et al. PRDM16 [1p36] translocations define a distinct entity of myeloid malignancies with poor prognosis but may also occur in lymphoid malignancies. Br J Haematol. 2012. https://doi.org/10.1111/j.1365-2141.2011.08918.x.

  41. Miyamura T, Moritake H, Nakayama H, Tanaka S, Tomizawa D, Shiba N, et al. Clinical and biological features of paediatric acute myeloid leukaemia (AML) with primary induction failure in the Japanese PaediatricLeukaemia/lymphoma study group AML-05 study. Br J Haematol. 2019. https://doi.org/10.1111/bjh.15799.

  42. Shiba N, Ohki K, Kobayashi T, Hara Y, Yamato G, Tanoshima R, et al. High PRDM16 expression identifies a prognostic subgroup of pediatric acute myeloid leukaemia correlated to FLT3-ITD, KMT2A-PTD, and NUP98-NSD1: the results of the Japanese Paediatric Leukaemia/lymphoma study group AML-05 trial. Br J Haematol. 2016. https://doi.org/10.1111/bjh.13869.

  43. Hu T, Morita K, Hill MC, Jiang Y, Kitano A, Saito Y, et al. PRDM16s transforms megakaryocyte-erythroid progenitors into myeloid leukemia-initiating cells. Blood. 2019. https://doi.org/10.1182/blood.2018888255.

  44. Tegeder I, Thiel K, Erkek S, Johann PD, Berlandi J, Thatikonda V, et al. Functional relevance of genes predicted to be affected by epigenetic alterations in atypical teratoid/rhabdoid tumors. J Neuro-Oncol. 2019. https://doi.org/10.1007/s11060-018-03018-6.

  45. Fei LR, Huang WJ, Wang Y, Lei L, Li ZH, Zheng YW, et al. PRDM16 functions as a suppressor of lung adenocarcinoma metastasis. J Exp Clin Cancer Res. 2019. https://doi.org/10.1186/s13046-019-1042-1.

  46. Deng J, Kong W, Mou X, Wang S, Zeng W. Identifying novel candidate biomarkers of RCC based on WGCNA analysis. Per Med. 2018. https://doi.org/10.2217/pme-2017-0091.

  47. Silva-Junior AL, Alves FS, Kerr MWA, Xabregas LA, Gama FM, Rodrigues MGA, et al. Acute lymphoid and myeloid leukemia in a Brazilian Amazon population: epidemiology and predictors of comorbidity and deaths. PLoS One. 2019. https://doi.org/10.1371/journal.pone.0221518.

Download references

Acknowledgements

We acknowledge all patients who participated in this study, Federal University of Pará, Evandro Chagas Institute for technical support and Coordination for Enhancement of Higher Education Personnel (CAPES) for fellowship support (Code 001).

Funding

This work was partially supported by the National Counsel of Technological and Scientific Development (CNPq), grant n° 460185/2014–4 and Amazon Foundation for Research Support (FAPESPA), grant n° PPSUS/2013.

Author information

Authors and Affiliations

Authors

Contributions

JG, FM and MS performed the experiments. AK, EO, NS and BK participated in coordination and helped to draft the manuscript. AW and LP were responsible for sample collection. All authors read and approved the final manuscript.

Corresponding author

Correspondence to André Salim Khayat.

Ethics declarations

Ethics approval and consent to participate

Octávio Lobo Children’s Cancer Hospital ethics committee approved this study (CAAE: 00905812.1.0000.00.18). Written consent forms were obtained from all parents of patients and enabled the collection of biological samples and publication of results.

Consent for publication

Consent forms from the parents of patients were obtained to publish and to report individual data.

Competing interests

The authors declared that they have no conflict of interest.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Batista-Gomes, J.A., Mello, F.A.R., de Oliveira, E.H.C. et al. Identifying novel genetic alterations in pediatric acute lymphoblastic leukemia based on copy number analysis. Mol Cytogenet 13, 25 (2020). https://doi.org/10.1186/s13039-020-00491-5

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s13039-020-00491-5

Keywords