Branka Zukic

Systematic documentation and analysis of human genetic variation in hemoglobinopathies using the microattribution approach

We developed a series of interrelated locus-specific databases to store all published and unpublished genetic variation related to hemoglobinopathies and thalassemia and implemented microattribution to encourage submission of unpublished observations of genetic variation to these public repositories. A total of 1,941 unique genetic variants in 37 genes, encoding globins and other erythroid proteins, are currently documented in these databases, with reciprocal attribution of microcitations to data contributors. Our project provides the first example of implementing microattribution to incentivise submission of all known genetic variation in a defined system. It has demonstrably increased the reporting of human variants, leading to a comprehensive online resource for systematically describing human genetic variation in the globin genes and other genes contributing to hemoglobinopathies and thalassemias. The principles established here will serve as a model for other systems and for the analysis of other common and/or complex human genetic diseases.

Functional analysis of the role of the TPMT gene promoter VNTR polymorphism in TPMT gene transcription

AIMS Thiopurine S-methyltransferase (TPMT) activity is polymorphic, and a trimodal distribution has been demonstrated in Caucasians (low, intermediate and high methylator groups). The TPMT gene promoter contains a variable number of three GC-rich tandem repeats, namely A, B and C, ranging from three to nine in length in a A(n)B(m)C architecture. MATERIALS & METHODS Here, we investigated the influence of number and type of TPMT gene promoter tandem repeats on human TPMT gene transcription in K562 cells transiently transfected with reporter constructs bearing various variable number of tandem repeats (VNTR) and addressed the interaction of transcription factor binding to the VNTRs by electrophoretic mobility shift assays. RESULTS We found that the distribution patterns of VNTR alleles do not significantly differ among acute lymphoblastic leukemia patients, acute myeloid leukemia patients and normal individuals. We also demonstrated that the A repeat has a negative effect in TPMT gene transcription and that a positive regulatory element, identified immediately upstream to the VNTR region of the TPMT gene promoter, is indispensable for TPMT gene transcription. Our electrophoretic mobility shift assay analysis indicated that the Sp1 and Sp3 transcription factors bind to the VNTR repeats. CONCLUSION Overall, our data underline that both the number and type of VNTRs, as well as the upstream regulatory region of the TPMT gene promoter, determine the overall level of TPMT gene transcription. It remains to be seen whether these VNTRs can be employed as pharmacogenetic markers to individualize thiopurine therapy.

A European spectrum of pharmacogenomic biomarkers: Implications for clinical pharmacogenomics

Pharmacogenomics aims to correlate inter-individual differences of drug efficacy and/or toxicity with the underlying genetic composition, particularly in genes encoding for protein factors and enzymes involved in drug metabolism and transport. In several European populations, particularly in countries with lower income, information related to the prevalence of pharmacogenomic biomarkers is incomplete or lacking. Here, we have implemented the microattribution approach to assess the pharmacogenomic biomarkers allelic spectrum in 18 European populations, mostly from developing European countries, by analyzing 1,931 pharmacogenomics biomarkers in 231 genes. Our data show significant inter-population pharmacogenomic biomarker allele frequency differences, particularly in 7 clinically actionable pharmacogenomic biomarkers in 7 European populations, affecting drug efficacy and/or toxicity of 51 medication treatment modalities. These data also reflect on the differences observed in the prevalence of high-risk genotypes in these populations, as far as common markers in the CYP2C9, CYP2C19, CYP3A5, VKORC1, SLCO1B1 and TPMT pharmacogenes are concerned. Also, our data demonstrate notable differences in predicted genotype-based warfarin dosing among these populations. Our findings can be exploited not only to develop guidelines for medical prioritization, but most importantly to facilitate integration of pharmacogenomics and to support pre-emptive pharmacogenomic testing. This may subsequently contribute towards significant cost-savings in the overall healthcare expenditure in the participating countries, where pharmacogenomics implementation proves to be cost-effective.

Genomic variation in the MAP3K5 gene is associated with β-thalassemia disease severity and hydroxyurea treatment efficacy

AIM In this study we explored the association between genetic variations in MAP3K5 and PDE7B genes, residing on chromosome 6q23, and disease severity in β-hemoglobinopathy patients, as well as the association between these variants with response to hydroxyurea (HU) treatment. Furthermore, we examined MAP3K5 expression in the context of high fetal hemoglobin (HbF) and upon HU treatment in erythroid progenitor cells from healthy and KLF1 haploinsufficient individuals. MATERIALS & METHODS For this purpose, we genotyped β-thalassemia intermedia and major patients and healthy controls, as well as a cohort of compound heterozygous sickle cell disease/β-thalassemia patients receiving HU as HbF augmentation treatment. Furthermore, we examined MAP3K5 expression in the context of high HbF and upon HU treatment in erythroid progenitor cells from healthy and KLF1 haploinsufficient individuals. RESULTS A short tandem repeat in the MAP3K5 promoter and two intronic MAP3K5 gene variants, as well as a PDE7B variant, are associated with low HbF levels and a severe disease phenotype. Moreover, MAP3K5 mRNA expression levels are altered in the context of high HbF and are affected by the presence of HU. Lastly, the abovementioned MAP3K5 variants are associated with HU treatment efficacy. CONCLUSION Our data suggest that these MAP3K5 variants are indicative of β-thalassemia disease severity and response to HU treatment.

6-mercaptopurine influences TPMT gene transcription in a TPMT gene promoter variable number of tandem repeats-dependent manner

AIM TPMT activity is characterized by a trimodal distribution, namely low, intermediate and high methylator. TPMT gene promoter contains a variable number of GC-rich tandem repeats (VNTRs), namely A, B and C, ranging from three to nine repeats in length in an A(n)B(m)C architecture. We have previously shown that the VNTR architecture in the TPMT gene promoter affects TPMT gene transcription. MATERIALS, METHODS & RESULTS: Here we demonstrate, using reporter assays, that 6-mercaptopurine (6-MP) treatment results in a VNTR architecture-dependent decrease of TPMT gene transcription, mediated by the binding of newly recruited protein complexes to the TPMT gene promoter, upon 6-MP treatment. We also show that acute lymphoblastic leukemia patients undergoing 6-MP treatment display a VNTR architecture-dependent response to 6-MP. CONCLUSION These data suggest that the TPMT gene promoter VNTR architecture can be potentially used as a pharmacogenomic marker to predict toxicity due to 6-MP treatment in acute lymphoblastic leukemia patients.

Expression of Bcl2L12 in chronic lymphocytic leukemia patients: association with clinical and molecular prognostic markers

Dysregulation of apoptosis is a distinctive feature of chronic lymphocytic leukemia (CLL), although a unique mechanism underlying apoptosis resistance of CLL B lymphocytes has not been identified yet. Aberrant expression as well as genetic and epigenetic alterations of numerous genes involved in different pathways of apoptosis regulation has been described in CLL. Here, we report the expression analysis of Bcl2L12 (Bcl2-like 12), a novel apoptotic gene belonging to Bcl2 family, in 58 Serbian CLL patients. Quantitative reverse-transcriptase polymerase chain reaction (qRT-PCR) analysis revealed a significant overexpression of Bcl2L12 mRNA in CLL samples compared to non-leukemic samples, implying its role in the pathogenesis of the disease. Receiver operating characteristic (ROC) analysis showed that Bcl2L12 expression efficiently discriminates CLL cases from healthy controls. However, relatively homogenous Bcl2L12 mRNA expression among patients did not reflect their clinical characteristics (with the exception of lactate dehydrogenase status and time from diagnosis to treatment) and failed to show association with the most informative prognostic markers, namely the mutational status of rearranged immunoglobulin heavy chain variable region genes, CD38 and lipoprotein lipase gene (LPL) expression.

Transcriptional regulation and pharmacogenomics

Interindividual variable drug response is correlated with sequence alterations in genes encoding drug-metabolizing enzymes and transporters, affecting drug absorption, distribution, metabolism and excretion. This variable drug response may have an impact on disease therapeutic outcomes, tolerance to adverse drug reactions and even survival. Sequence alterations may occur not only within the coding region of a gene, but in its regulatory elements too, affecting gene transcription and gene-expression levels. Here, we provide a compilation of the current knowledge of pharmacogenomics related to transcription, with a focus on the effect of SNPs and short tandem repeats residing in cis-regulatory elements of 11 genes encoding for drug-metabolizing enzymes and drug transporters. In addition, we comment on two genes encoding enzymes that are drug targets themselves. Finally, we briefly discuss the currently available methodologies for clinically assessing pharmacogenomic profiles, which could potentially in the future facilitate drug treatment-individualization via the identification of molecular signatures in specific patient groups.

Novel transcriptional regulatory element in the phenylalanine hydroxylase gene intron 8

We present the first transcriptional regulatory element found in a PAH gene intron. The element is located in the PAH gene intron 8, acts as an enhancer specifically in the hepatoma cell line, and binds GATA-1 transcription factor. Herein the presented data could unlock a new area for the analysis of PAH gene expression and could contribute to refining genotype-phenotype correlation.

The influence of novel transcriptional regulatory element in intron 14 on the expression of Janus kinase 2 gene in myeloproliferative neoplasms

The expression of Janus kinase 2 (JAK2) gene is altered in myeloproliferative neoplasms (MPN) and the regulation of transcription could be a mechanism that modulates JAK2 gene expression. We analyzed the transcriptional potential of single-nucleotide polymorphism (SNP) rs12343867 T > C in JAK2 intron 14, tagging 46/1 haplotype, and its influence on JAK2 gene expression. Functional analysis of JAK2 intron 14 was performed using the pBLCAT5 reporter system in K562 cells. Identification of the proteins binding to the intron 14 regulatory element was accomplished by electrophoretic mobility shift assay (EMSA) and supershift assays. Quantification of the expression of JAK2 gene in a cohort of 51 MPN patients and 12 healthy controls was performed by real-time quantitative polymerase chain reaction (RQ-PCR). Functional analysis revealed that the intronic DNA element harboring SNP rs12343867 T > C acts as a transcriptional repressor in vitro. The repressor activity was significantly attenuated by the presence of nucleotide C. Supershift analysis showed the enrolment of transcriptional factor Meis1 in this process. RQ-PCR experiments showed increased JAK2 expression in patients with the JAK2V617F mutation, with a significant difference between essential thrombocythemia (ET), polycythemia vera (PV), and myelofibrosis (MF) patients. SNP rs12343867 showed no statistically significant influence on the expression of JAK2 gene in MNP patients.

Functional analysis of a novel KLF1 gene promoter variation associated with hereditary persistence of fetal hemoglobin

Hereditary persistence of fetal hemoglobin (HPFH) is a rare hereditary condition resulting in elevated levels of fetal hemoglobin (HbF) in adults. Typical HPFH is associated with promoter mutations or large deletions affecting the human fetal globin (HBG1 and HBG2) genes, while genetic defects in other genes involved in human erythropoiesis, e.g. KLF1, also result in atypical HPFH. Here, we report the first KLF1 gene promoter mutation (KLF1:g.-148G > A) that is associated with increased HbF level. This mutation was shown to result in drastically reduced CAT reporter gene expression in K562 cells, compared to the wild-type sequence (p = 0.009) and also in reduced KLF1 gene expression in vivo. Furthermore, consistent with in silico analysis, electrophoretic mobility shift analysis showed that the KLF1:g.-148G > A mutation resides in a Sp1 binding site and further that this mutation leads to the ablation of Sp1 binding in vitro. These data suggest that the KLF1:g-148G > A mutation could play a role in increasing HbF levels in adults and further underlines the role of KLF1 as one of the key transcription factors involved in human fetal globin gene switching.