GENE EXPRESSION PROFILING OF B-CLL IN UKRAINIAN PATIENTS IN POST-CHERNOBYL PERIOD
Summary. Background: Genetic mechanisms that result in the development and progression of B-cell chronic lymphocytic leukemia (B-CLL) are mainly unknown. We have analyzed gene expression patterns in Ukrainian B-CLL patients with the aim of identifying B-CLL involved / associated genes in order to shed light on the biology of this pathological entity. Material and methods: The samples of the peripheral blood and bone marrow of 44 Ukrainian B-CLL patients with no characteristics indicative of unfavorable course of the disease such as CD38 were analyzed morphologically and immunocytochemically according to the new WHO classification. Total RNA was isolated, and gene expression levels were determined by microarray method comparing with the sample from 17 healthy donors. Results: We investigated interactions using the Ingenuity Pathway Analysis (IPA) software and found 1191 network eligible up-regulated genes and 3398 Functions/Pathways eligible up-regulated genes, 1225 network eligible down-regulated genes and 2657 Functions/Pathways eligible down-regulated genes. Conclusion: In B-CLL patients, gene networks around MYC, HNF1A and HNF4A, YWHAG, NF-κB1 and SP1 are identified as up-regulated; CEBPA, YWHAG, SATB1 and RB1 — as down-regulated. G protein coupled receptor signaling, arachidonic acid and linoleic acid metabolisms, calcium signaling, metabolism of xenobiotics by cytochrome P450 are found out as significant up-regulated pathways. EIF2 and Cdc42 signaling, regulation of eIF4 and p70S6k signaling, protein ubiquitination pathway and oxidative phosphorylation are the most significant down-regulated pathways obtained in our study. The involvement of NF-κB gene network and upregulated levels of G protein coupled receptor signaling pathway, which has an important role in transcription of NF-κB, are important and need further examination.
Received: February 2, 2012.
*Correspondence: E-mail: email@example.com
Abbreviations used: AML — acute myelogenous leukemia; B‑CLL — B-cell chronic lymphocytic leukemia; IPA — Ingenuity Pathway Analysis; IR — ionizing radiation.
Ionizing radiation (IR) is one of the most studied carcinogens in the development of multiple myeloma, primary myelofibrosis, polycythemia vera, non-Hodgkin’s lymphomas, myelodysplastic syndromes and some forms of acute and chronic leukemia, especially in acute myelogenous leukemia (AML) [1, 2]. Until recently, chronic lymphocytic leukemia (CLL) has not been considered as a radiation-associated leukemia. Nevertheless, current understanding of radiation-induced tumorigenesis and the etiology of lymphatic neoplasia show that IR exposure increases CLL risk .
After Chernobyl nuclear accident, people living in the contaminated areas of Ukraine are still exposed to low doses of IR. Analysis of the patients with various forms of the malignancies of hematopoietic and lymphoid tissues has not revealed the differences in B-CLL percentage among Chernobyl clean-up workers and Ukrainian population in whole. B-CLL was shown to be a predominant form of hematopoietic malignancies in clean-up workers as well as in general population . Genetic mechanisms that result in the development and progression of CLL are mainly unknown . Gene expression profiling by microarray is useful to understand B-CLL origin and development . The analysis of the molecular genetic features should be advantageous in elucidating the putative association of IR and B-CLL.
Earlier, we have studied gene expressions of several apoptosis related genes in different types of tumors of hematopoietic and lymphoid tissues in 189 patients including those with B-CLL living in areas of Ukraine contaminated with radionuclides in post-Chernobyl period . In the present study, we have analyzed gene expression patterns in samples from 44 B-CLL Ukrainian patients in post-Chernobyl period with the aim of identifying the genes related to or involved in this pathology in order to shed light on the biology of B-CLL.
MATERIAL AND METHODS
The samples of the peripheral blood of B-CLL patients were obtained from R.E. Kavetsky Institute of Experimental Pathology, Oncology and Radiobiology of the National Academy of Sciences of Ukraine. All the patients were referred to Reference Laboratory of Immunocytochemistry and Oncohematology Department of the Institute for verifying the diagnosis. Bone marrow and peripheral blood smears stained by May-Grunwald-Giemsa were studied morphologically. Immunocytochemical techniques (APAAP, LSAB-AP) and a broad panel of monoclonal antibodies against lineage specific, differentiation and activation antigens of leukocytes were employed for immunophenotyping pathological cells in blood and bone marrow. The main forms and cytological variants of hematological malignancies were diagnosed according to new WHO classification . All the samples were immunophenotyped, and only 44 samples from CD38-negative B-CLL patients out of 127 diagnosed patients with B-CLL/B-cell lymphoma from small lymphocytes  were included in the study. Control group comprised peripheral blood samples from 17 healthy donors. The ethic committees of both collaborating research institutions approved the design of the study.
Total RNA isolation. Total RNA was isolated from mononuclear cells for each patient using QIAamp RNA Blood Mini Kit (QIAGEN, Valencia, CA, USA) and treated with DNase I according to the manufacturer’s instructions. The quality of the RNA was assessed by loading 300 ng of total RNA onto an RNA Labchip (Agilent Technologies, Waldbronn, Germany), followed by analysis (A2100 Bioanalyzer; Agilent Technologies). An RNA integrity value (RIN) of 7.0 was considered acceptable.
RNAs from 44 B-CLL patients and 17 healthy donors were pooled seperately. Pooling process was performed in the way that 100 ng RNA sample was used from each B-CLL patient/healthy donor. Each RNA pool was prepared as three replicates.
Microarray analysis. Microarray analysis was performed using the Whole Human Genome Oligo Microarray (Agilent Technologies), encompassing more than 44,000 human DNA probes. The full list of cDNAs is available online (www.agilent.com). Protocols for sample preparation and hybridization of the mononuclear cells were adaptations of those in the Agilent Technical Manual. In short, first strand cDNA was transcribed from 300 ng of total RNA using T7-Oligo(dT) Promoter Primer. Samples were transcribed in vitro and Cy-3-labelled by using a Quick-AMP labeling kit (Agilent Technologies). Following a further clean-up round (Qiagen), cRNA was fragmented into pieces ranging from 35 to 200 bases in size. Fragmented cRNA samples (1.65 mg) were hybridized onto chips by means of 17 h of incubation at 65°C with constant rotation, followed by a two-step microarray wash of 1 min in two washing buffers (Agilent Technologies). Hybridized microarrays were scanned in a Agilent Technologies Scanner (model G2505B) and numerical results were extracted with Feature Extraction version 126.96.36.199 using 014850_D_F_20060807 grid, GE1-v5_95_Feb07 protocol and GE1_QCM_Feb07 QC metric set.
The microarray data were analyzed using GeneSpring software version 9.0 (Agilent Technologies, Santa Clara, CA). The fold changes were analyzed by filtering the dataset using P-value < 0.01 and a signal-to-noise ratio >2 for use in T-test statistical analysis. Additional filtering (minimum 2-fold change) was applied to extract the most these genes, which were analyzed using Ingenuity Pathway Analysis (IPA) software (Ingenuity Systems, Redwood City, CA). Those genes with known gene symbols (HUGO) and their corresponding expression values were uploaded into the software. Each gene symbol was mapped to its corresponding gene object in the Ingenuity Pathways Knowledge Base. Networks of these genes were algorithmically generated based on their connectivity and assigned a score. The score is a numerical value used to rank Networks according to how relevant they are to the genes in the input dataset but may not be an indication of the quality or significance of the network. The score takes into account the number of focus genes in the network and the size of the network to approximate how relevant this network is to the original list of focus genes. The network identified is then presented as a graph indicating the molecular relationships between genes/gene products. Genes are represented as nodes, and the biological relationship between two nodes is represented as an edge (line). The intensity of the node color indicated the degree of up- or down-regulation. The node shapes are disclosed in corresponding figure legends. Canonical pathway analysis identified the pathways from the IPA library of canonical pathways, which were most significant to the input data set. The significance of the association between the data set and the canonical pathway was determined based on two parameters: (1) A ratio of the number of genes from the data set that map to the pathway divided by the total number of genes that map to the canonical pathway and (2) a P value calculated using Fischer’s exact test determining the probability that the association between the genes in the data set and the canonical pathway is due to chance alone.
Quantitative real-time PCR (Q-RT-PCR). cDNA was synthesized using RevertAid First Strand cDNA Synthesis Kit (Fermentas Inc., Maryland, USA). Q-RT-PCR was performed as we described previously for determination of MYC, BAX, BCL-2 and FAS1 gene expressions [9, 10]. Standard curves were obtained using serial dilutions of the beta-globulin gene (DNA Control Kit, Roche). Gene-specific primers (Table 1) were obtained from Integrated DNA Technologies (Iowa, USA). Obtained gene expression values were normalized using a housekeeping gene of beta2 microglobulin. Gene expression ratios were compared in patient and control groups using REST (Relative Expression Software Tool).
Table 1. List of the primers used for the quantitative RT-PCR
Differentially expressed genes are shown in two separate tables. The 100 most up-regulated genes are shown in Table 2. The 100 most down-regulated genes are shown in Table 3. Both sets of results were obtained based on minimum 2-fold change using GeneSpring software version 9.0 (Agilent Technologies, Santa Clara, CA). In Table 4 the gene expression results of four genes (MYC, BAX, BCL-2 and FAS1) obtained by real-time PCR and microarray methods are compared. Real-time PCR results of MYC, BCL-2 and BAX are in a good agreement with microarray expression rates.
Table 2. The 100 most up-regulated genes in B-CLL
Table 3. The 100 most down-regulated genes in B-CLL
Table 4. Summarized real-time PCR confirmation results of the four genes
We investigated interactions using IPA software and found 1191 network eligible up-regulated genes and 3398 Functions/Pathways eligible up-regulated genes. Fig. 1 shows the most significant four gene networks of over-expressed genes in B-CLL samples. Top functions of these genes are related to hematopoiesis, lipid metabolism, small molecule biochemistry, cancer, infectious diseases, cell cycle, cardiovascular system development and function, gene expression, embryonic development, tissue morphology, inflammatory response. Up-regulated gene networks are identified around MYC, HNF1A and HNF4A, YWHAG, NF-κB1 and SP1.
We also found 1225 network eligible down-regulated genes, and 2657 Functions/Pathways eligible down-regulated genes. Fig. 2 shows four gene networks of down-regulated genes in B-CLL. The main functions of the genes are related to cellular functions and maintenance, protein synthesis, dermatological diseases and conditions, cell death, gene expression, inflammatory disease, cellular growth and proliferation, post-translational modification, cancer, infectious diseases, cell morphology, and development. Down-regulated gene networks are identified around CEBPA, YWHAG, SATB1 and RB1.
B-CLL is a heterogeneous disease and a predominant form of hematopoietic malignancies. Despite new molecular methods identifying important prognostic and diagnostic genetic markers, genetic mechanisms involved in B-CLL origin are mainly unknown. A number of novel prognostic markers such as Bcl-2, MAP-kinase, NF-κB, ZAP-70 were identified applying gene expression profiling before [11, 12].
We have analyzed gene expression patterns in samples from 44 B-CLL Ukrainian patients in post-Chernobyl period to identify genes associated with this form of lymphoproliferative malignancy. Our study has demonstrated new genetic networks and biological pathways in both up- and down-regulated gene expression levels.
Analysis using IPA software revealed 1191 network eligible up-regulated genes and 3398 Functions/Pathways eligible up-regulated genes. The individual genes are found in multiple categories of functions related to hematopoiesis, lipid metabolism, small molecule biochemistry, cancer, infectious diseases, cell cycle, development and function of cardiovascular system, gene expression, embryonic development, tissue morphology, inflammatory response.
One important gene network is identified around the up-regulated MYC and SP1 genes (Fig. 1, a). MYC, a strong proto-oncogene, plays very important roles in cellproliferation (by upregulating cyclins, downregulating p21), controlling cell growth (by upregulating ribosomal RNA and proteins), apoptosis (by downregulating BCL-2), differentiation and stem cell self-renewal. Mutations, overexpression, rearrangement and translocation of MYC have been associated with a variety of hematopoietic tumors, leukemias and lymphomas, including Burkitt lymphoma . High expression level of MYC has been reported in more aggressive and apoptosis resistant forms of B-CLL and might be used as molecular marker specific of resistant B-CLL subsets [14, 15]. SP1, a zinc finger transcription factor, is involved in cell differentiation, cell growth, apoptosis, immune response, response to DNA damage, and chromatin remodeling. SP1 and MYC are involved cooperatively in telomerase activation, which is a critical step in cellular immortalization and carcinogenesis. Kyo et al. have suggested that the level of SP1 expression might be a critical determinant of telomerase activity both in cancer and normal cells .
Another network is identified around NF-κB1 gene (Fig. 1, b). NF-κB regulates several genes that mediate tumorigenesis and metastasis and also plays an important role in pathogenesis of B-cell neoplasms. Carcinogens, tumor promoters, inflammatory cytokines, and chemotherapeutic agents activate NF-κB and this activation can suppress apoptosis, thus promoting chemoresistance and tumorigenesis. Bharti et al. suggested that NF-κB might be an ideal target for chemoprevention and chemosensitization [17, 18]. In addition, we have found NF-κB gene centered around two up- and down-regulated networks in our previous study on prostate cancer .
Canonical pathway analysis revealed that G-protein coupled receptor (GPCR) signaling is an important pathway modulated by the up-regulated genes in B-CLL. It is known that GPCRs regulate proliferation, differentiation, chemotaxis and also they play an important role in inflammatory diseases and cancer . GPCRs are involved in control of transcription factors such as STAT3, NF-κB and CREB by G protein subfamilies . Enhanced viability of CLL cells by the STAT3 phosphorylation and interaction between hepatocyte growth factor and its receptor (c-MET), which was found up-regulated in our study, was reported before . CREB (cAMP response element binding protein) had been found overexpressed in bone marrow samples from patients with acute lymphoid and myeloid leukemia and associated with a poor outcome in AML patients according to previous studies .
A network is also identified around HNF1A and HNF4A (Fig. 1, c). HNF1A is a transcription factor required for the expression of several liver-specific genes and the expression of this gene is controlled by HNF4A, which may play role in development of the liver, kidney and intestines.
Another significant signaling pathway is calcium signaling involved in many processes such as cell survival/apoptosis, cell cycle progression, differentiation, cross-talk between intracellular compartments (ER, mitochondria), general metabolism and telomerase activity. The calcineurin/NFAT signaling pathway is important in lymphoma/leukemogenesis . Deregulation of this signaling and/or abnormal expression of its components has been reported in solid tumors of epithelial origin, lymphoma and lymphoid leukemia. Mouse models of human T-ALL/lymphoma showed the pro-oncogenic effect of active calcineurin/NFAT signaling in vivo . NFAT transcription factors form four calcium signaling responsive members: NFATc1, NFATc2, NFATc3 and NFATc4. Among these members NFATc1 and NFATc2, which are found up-regulated in our study, were reported to be involved in the development, differentiation and functioning of multiple T-and B-cell subsets in previous studies. NFATc1 was found to be expressed in a majority of aggressive B-cell lymphomas. On the other hand, NFATc2 activation was shown to be responsible in B-CLL, in cooperation with STAT6, for the high expression of CD23 .
Metabolism of xenobiotics by cytochrome P450 pathway has been shown as highly significant in our study. The enhanced expression of several P450s like CYP1A, CYP2C and CYP3A, that are up-regulated in our study, was reported in tumor cells elsewhere .
Arachidonic acid and linoleic acid metabolisms are the other significant pathways modulated by the up-regulated genes in our study.
Analysis using IPA software revealed 1225 network eligible down-regulated genes, and 2657 Functions/Pathways eligible down-regulated genes. These individual genes are related to cell functions and maintenance, protein synthesis, dermatological diseases and conditions, cell death, gene expression, inflammatory disease, cell growth and proliferation, post-translational modification, cancer, infectious diseases, cell morphology, and development.
One important down-regulated network is identified around RB1 gene (Fig. 2, a). The role of RB1 in B-CLL has been reported based on cytogenetic data . RB1 deletions involved in 13q14 abnormalities have been reported in B-CLL before .
Another down-regulated network is identified around SATB1 gene (Fig. 2, b). SATB1 is a new type of gene regulator expressing in various human cancers and thought to be related to the malignant potential. Overexpression of this gene has been reported as a predictor of poor prognosis in lung and gastric cancers [29, 30].
An important network is identified around CEBPA gene (Fig. 2, c). CEBPA is a critical transcriptional factor and regulates the balance between cell proliferation and differentiation during early hematopoietic development and myeloid differentiation . CEBPA has a tumor-suppressor function in leukemogenesis and both loss of function and gain of function have leukemogenic potential. It was reported that overexpression of CEBPA could contribute to B-ALL and loss of function could contribute to AML . On the other hand, down-regulated CEBPA was found in acute promyelocytic leukemia stem cells in animal models .
Canonical pathway analysis revealed that oxidative phosphorylation is an important pathway modulated by the down-regulated genes in B-CLL. In fact, previous studies suggested that the oxidative phosphorylation (OXPHOS) system is severely compromised in various cancers .
EIF2 signaling is another significant pathway. Suppression of head and neck, colorectal carcinoma and multiple myeloma tumor growth and/or survival by phosphorylation of eIF2α was reported before .
IPA reveals regulation of eIF4 and p70s6K signaling pathway. eIF4E down-regulated in our study plays an important role in tumor initiation and progression when its overexpression cooperates with oncogenes to accelerate transformation in cell lines and animal models . p70s6K is a serine/threonine kinase and its target substrate is S6 ribosomal protein . Inhibition of p70s6K was related to cell cycle arrest at G0/G1 phase in human cancer cells before .
Protein ubiquitination is another pathway found significant in our study. Ubiquitination of key signaling molecules by E3 ubiquitin ligases forms an important regulatory mechanism for NF-κB signaling. Deubiquitinases (DUBs) counteract E3 ligases and play a substantial role in down-regulation of NF-κB signaling and homeostasis .
Cdc42 signaling is a highly significant pathway. Cdc42 promotes or inhibits tumor progression depending on the cellular context and contributes to cancer development through its different roles in intracellular trafficking, cell cycle regulation and survival, polarity, migration and transcriptional control . Cdc42 is also important in the development and progression of lymphoma. Genetic knockdown or pharmacological inhibition of Cdc42 resulting in a cell cycle arrest and apoptosis of anaplastic large cell lymphoma cells has been reported .
An important network is identified around both down-regulated and up-regulated YWHAG gene in our study (Fig. 1 and Fig. 2). This gene encoding for 14–3-3 protein gamma was found highly expressed in skeletal and heart muscles. It has been suggested that this protein has an important role in muscle tissue [42, 43]. 14–3-3 proteins play critical regulatory roles in signaling pathways in cell division and apoptosis . Further investigations are required to establish the function of YWHAG gene in B-CLL.
Fig. 1. Significant up-regulated gene networks identified around MYC and SP1 genes (a), NF-κB1 gene (b), HNF1A and HNF4A genes (c), YWHAG gene (d) in B-CLL samples. The node shapes denote enzymes (), phosphatases (), kinases (), peptidases
(), G-protein coupled receptor (), transmembrane receptor (), cytokines (), growth factor (), ion channel (), transporter (), translation factor (), nuclear receptor (), transcription factor () and other ().The intensity of the node color-(red) indicated the degree of up-regulation
Fig. 2. A significant down-regulated gene network identified around RB1 gene (a), SATB1 gene (b) in B-CLL samples, CEBPA gene (c), YWHAG gene (d) in B-CLL samples. The keys to the node shapes are the same as in Fig. 1. The intensity of the node color-(green) indicates the degree of down-regulation
Real-time PCR confirmation results of four genes (MYC, FAS1, BAX and BCL-2) show that only MYC, BAX and BCL-2 expressions are in agreement with microarray results. Up-regulation of MYC is compatible with our expectations.
Previous studies indicate that high ratio of Bcl-2 to Bax proteins confers a poor prognosis with decreased rates of complete remission and overall survival . In our study, BCL-2 upregulation level is superior to that of BAX in real-time PCR results but not in microarrays being analyzed.
FAS1 expression was found up-regulated in real-time PCR but down-regulated in microarrays in our study. It has been reported that Fas expression is not very high in B-CLL  that coincides with our findings of relatively small up-regulation by real-time PCR. It should be noted that Fas was mentioned as apoptosis regulator  in CLL cells exposed to IR.
NF-κB gene network was conspicuous in terms of being determined also in our previous studies of gene expression in prostate cancer. In addition, upregulated levels of G protein coupled receptor signaling pathway, which has an important role in transcription of NF-κB, need advanced examinations. In this sense, NF-κB gene which is important in both cell cycle regulation and cancer progression deserves further study.
Our study has presented the gene expression profiling in B-CLL patients of Ukrainian population as whole. We believe that the contribution of IR as the putative factor in the origin of B-CLL should be further evaluated using such molecular genetic approach.
The study was financed within the framework of the joint research project M/32–2008 “Cytomorphological, immunocytochemical and molecular biological features of leukemias in persons exposed to ionizing radiation” according to the Agreement between the Ministry of Education and Science of Ukraine and the Scientific and Technical Research Council of Turkey (TUBITAK).
1. Gluzman DF, Nadgornaya VA, Machilo V, et al. Malignant diseases of hematopoietic and lymphoid tissues in Chernobyl clean-up workers. Hematol J 2005; 5: 565–7.
No Comments » Add comments