Original Article Influence of Variants of the Drosha, Mir499a, and Mir938
Total Page:16
File Type:pdf, Size:1020Kb
Am J Transl Res 2020;12(12):8216-8224 www.ajtr.org /ISSN:1943-8141/AJTR0121284 Original Article Influence of variants of thedrosha , mir499a, and mir938 genes on susceptibility to acute lymphoblastic leukemia in an admixed population from the brazilian amazon Tatiane Piedade de Souza2, Darlen Cardoso de Carvalho1, Alayde Viera Wanderley1,4, Sweny Marinho Fernandes1,4, Juliana Carla Gomes Rodrigues1, Amanda Cohen-Paes1, Marianne Rodrigues Fernandes1, Fernando Augusto Rodrigues Mello Junior5, Lucas Favacho Pastana1, Lui Wallacy Morikawa Souza Vinagre1, Artur Luiz da Costa Silva3, Paulo Pimentel de Assumpção1, Sidney Santos1,2, André Salim Khayat1, Ney Pereira Carneiro dos Santos1,2 1Research Center of Oncology, Federal University of Pará, Belém, Pará, PA, Brazil; 2Laboratory of Human and Medical Genetics, Institute of Biological Sciences, Belém, Pará, PA, Brazil; 3PTC Laboratory of Biological Engineering-Guamá, Belém, Pará, PA, Brazil; 4Otávio Lobo Children’s Cancer Hospital-HOIOL-Belém, Pará, PA, Brazil; 5Molecular Biology Laboratory of The Ophir Loyola Hospital-HOL-Belém, Pará, PA, Brazil Received August 28, 2020; Accepted October 25, 2020; Epub December 15, 2020; Published December 30, 2020 Abstract: Acute Lymphoblastic Leukemia (ALL) is the most common type of cancer in children. Polymorphisms that alter the normal function of the microRNAs involved in the development of ALL have been widely investigated, although published data on these polymorphisms in admixed populations are scarce. We investigated the role of 10 polymorphisms in the microRNA and protein-coding genes of the microRNA synthesis complex in susceptibility to pediatric B-cell ALL. The study includes 100 pediatric ALL patients and 180 healthy individuals. The statistical analyses were run in SPSS v.25.0. In the case of the microRNA synthesizing genes, a significant pattern was found in only gene, that is, the rs3805500 polymorphism of DROSHA, in which the homozygous mutant (AA) genotype was associated with a threefold increase in the risk of developing ALL when compared to other genotypes (P=0.004, OR=2.913, CI=1.415-5.998). In the microRNA coding genes, the homozygous mutant rs3746444 genotype of the MIR499A gene was associated with a 17-fold increase in the risk of development of ALL (P<0.001, OR=17.797, CI=5.55-57.016). A protective effect against the development of ALL was also observed in the carriers of the wild homozygous rs2505901 genotype in the MIR938 gene. Our findings highlight the potential of these polymorphisms in the genes involving in the coding of microRNAs for the evaluation of the risk of contracting ALL in the population of the Brazilian Amazon region. These findings contribute to a more complete understanding of the complex etiology of ALL. Keywords: Acute lymphoblastic leukemia, microRNAs, single nucleotide polymorphism, susceptibility Introduction carcinogens, as well as chromosomal and molecular alterations [4, 5]. In particular, many Acute Lymphoblastic Leukemia (ALL) is the single nucleotide polymorphisms (SNPs) in key most common type of cancer in children and genes of the regulatory pathways have been adolescents of up to 19 years of age, repre- ascribed a fundamental role in the develop- senting 28% of all malignancies in this age ment of ALL [6]. group and 75% of all cases of leukemia [1-3]. The etiology of ALL is considered to be mul- The investigation of the SNPs found in the tifactorial, involving both environmental and microRNA genes has emerged as a promising genetic factors, which have been investigated new field of genomic research. The principal extensively. These factors include exposure to function of these molecules is to regulate gene Polymorphisms in microRNA genes and the ALL risk expression via post-transcriptional silencing, Selection of the study polymorphisms the cleaving of messenger RNAs (mRNA) or inhibiting the initiation of translation through The SNPs selected for the present study were base pairing between the microRNA and its chosen due to their known association with mRNA target [7]. The presence of SNPs in susceptibility to ALL and different types of microRNA genes or the protein-coding genes cancer, based on the published data. The involved in the synthesis of the microRNAs references of the journal articles used for the may affect their correct functioning, resulting selection of the polymorphisms are listed in in impacts on gene regulation processes [8]. the Table S1 [38-48]. A total of 10 polymor- Despite the important regulatory potential of phisms were selected, seven in microRNA these polymorphisms, few studies have fo- genes and three in the protein-coding genes cused on the influence of the SNPs located in essential to the synthesis of microRNAs: rs- non-coding genomic regions. In addition, evi- 636832 (AGO1), rs10035440 e rs3805500 dence on the occurrence of these polymor- (DROSHA), rs213210 and rs107822 (MIR- phisms in populations with high levels of mis- 219-1), rs2910164 (MIR146a), rs12894467 cegenation, such as that of the Brazilian (MIR300), rs3746444 (MIR499a), rs4919510 Amazon region, is particularly scarce. Given (MIR608), and rs2505901 (MIR938). this, the present study investigated the role of Extraction and quantification of the DNA 10 polymorphisms of the microRNA and microRNA synthesis genes in the susceptibility The DNA was extracted from peripheral blood of pediatric patients from the Brazilian Amazon using a commercial DNA extraction kit (Bio- region to B-cell ALL. pur Mini Spin Plus-250 extraction kit, Biopur, Brazil). The concentration of the genetic mate- Material and methods rial was quantified in a NanoDrop 1000 spec- trophotometer (Scientific Term NanoDrop Patients and controls 1000; NanoDrop Technologies Wilmington, The case group of the present study included DE). 100 patients diagnosed with B-cell ALL by Genotyping of polymorphisms immunophenotyping and/or molecular analy- sis. All these individuals were in treatment The 10 polymorphisms were genotyped by at the Otávio Lobo Hospital in Belém, a refer- allelic discrimination using the TaqMan Open- ence center for the treatment of pediatric can- Array Genotyping technology with a set of 32 cer in the northern region of Brazil. The pa- customized assays, which were run in a tients included in the study were between 1 QuantStudio™ 12K Flex Real-Time PCR sy- and 18 years old. Recurrent patients or those stem (Applied BiosystEM, Life Technologies, with comorbidities were not included in the Carlsbad, USA), according to the manufactur- group. The control group consisted of 180 er’s protocol. This method is based on real- individuals with no diagnosis for B-cell ALL or time Polymerase Chain Reaction (qPCR). The any other type of cancer. None of these indi- quality of the readings of the genotypes and viduals were related to any of the patients in other data were analyzed in the TaqMan the case group, and they were all over the peak Genotyper software. risk age for the development of ALL. Quality control Ethical aspects To ensure an adequate level of accuracy, poly- This study was approved by the Research morphisms were only included in the present Ethics Committee of the Research Center of study if they satisfied three criteria: (i) MAF≥ Oncology of the Federal University of Pará 1%; (ii) genotyping rate ≥80%, and (iii) were in (CAAE number 11433019.5.0000.5634). All Hardy-Weinberg (HWE). The HWE was per- participants (or their legal guardians) signed a formed using the Arlequin software (v.3.5.1.2). term of informed consent authorizing the col- The significance of the HWE test was adjust- lection of samples and data. ed for multiple comparisons by the Bonferroni 8217 Am J Transl Res 2020;12(12):8216-8224 Polymorphisms in microRNA genes and the ALL risk Table 1. Primers used for RT-PCR in the chromosomal USA) with the primers designed specifi- fusion analysis cally for the five fusions mentioned above, Cromossomic Fusion Genes Primer (5’-3’) also following the manufacturer’s in- structions. The primers used for RT-PCR t(1;19)(q23;p13) E2A CTACTCCCCGGATCACTCAA are listed in Table 1. PBX1 AGGCTTCATTCTGTGGCAGT t(4;11)(q21;q23) MLL CGCCCAAGTATCCCTGTAAA Analysis of genomic ancestry AF4 GAGCATGGATGACGTTCCTT t(9;22)(q34;q11) BCR TCGCAGAACTCGCAACAGT The genomic ancestry of the partici- ABL ACACCATTCCCCATTGTGAT pants was analyzed following Santos et t(12;21)(p13;q22) TEL TCTCTCATCGGGAAGACCTG al. (2010) and Ramos et al. (2016), us- ing a set of 61 Ancestry Informative AML1 TGCGGTAGCATTTCTCAGC Markers (AIMs). The individual and glo- del(1)(p32;p32) SIL TCCTACCCTGCAAACAGACC bal proportions of European, Amerindian, TAL1 AGGCGGAGGATCTCATTCTT and African genetic ancestry were esti- mated using STRUCTURE v.2.3.4 [9-11]. Table 2. Demographic parameters for the case (pa- Statistical analyses tients with B-cell ALL) and control groups analyzed in the present study The Chi-square test was applied for the Variable Case (100) Control (180) p-value pairwise comparisons of the categorical Sex (M/F) 60/40 53/127 <0.001a variables between the case and control Age* 5.53 ± 3.991 65.97 ± 16.021 <0.001b groups, while the quantitative variables Genetic ancestry* were compared using Student’s t. The European 0.429 ± 0.133 0.454 ± 0.170 0.189c multivariate analyses considered the sex of the participants and Amerindian African 0.203 ± 0.089 0.241 ± 0.138 0.213c ancestry as confounding variables. The Amerindian 0.361 ± 0.154 0.304 ± 0.149 0.004c Mann-Whitney test was used to com- aChi-square, bStudenta0 t, cMann-Whitney9trU. *Mean ± eStandard Deviation. pare the estimates of genetic ancestry between the groups. All statistical analy- ses were run in SPSS v.25.0 (SPSS, method (P≤0.001), and the HWE values are Chicago, IL, USA), considering a P≤0.05 signifi- shown in the Table S2. As all 10 of the polymor- cance level. phisms selected for the present study satisfied these criteria, they were all included in the Results association analyses (Table S2).