Research Paper

Pdf file on Synergy

Download Chinese abstract

Acta Biochim Biophys Sin 2005,37:618-624

doi:10.1111/j.1745-7270.2005.00086.x

cDNA Expression Array Analysis of Gene Expression in Human Hepatocarcinoma­ Hep3B Cells Induced By BNIPL-1

 

Li XIE1,2, Wen-Xin QIN2,3*, Jin-Jun LI2,3, Xiang-Huo HE2, Hui-Qun SHU2, Gen-Fu YAO2, Da-Fang WAN2,3*, and Jian-Ren GU2,3

 

1 Shanghai Medical College, Fudan University, Shanghai 200032, China;

2 National Laboratory for Oncogenes and Related Genes, Shanghai Cancer Institute, Shanghai 200032, China;

3 Cancer Institute, Shanghai Jiaotong University, Shanghai 200032, China

 

Received: May 9, 2005

Accepted: June 14, 2005

This work was supported by the grants from the Major State Basic Research Key Program of China (No. 2002CB513100), the National High Technology Research and Development Program of China (No. 2001AA221141, No. 2002BA711A02 and No. 2004AA227090), and the NSFC-RGC Joint Fund (30218006)

*Corresponding authors:

Wen-Xin QIN: Tel/Fax, 86-21-64432142; E-mail, [email protected]

Da-Fang WAN: Tel/Fax, 86-21-64177401; E-mail, [email protected]

 

Abstract        Bcl-2/adenovirus E1B 19 kDa interacting protein 2 like-1 (BNIPL-1) is a novel human protein identified in our laboratory, which can interact with Bcl-2 and Cdc42GAP and induce apoptosis via the BNIP-2 and Cdc42GAP homology (BCH) domain. In the present study, we established the Hep3B-Tet-on stable cell line in which expression of BNIPL-1 can be induced by doxycycline. The cell proliferation activity assay showed that the overexpression of BNIPL-1 suppresses Hep3B cell growth in vitro. The differential expression profiles of 588 known genes from BNIPL-1-transfected Hep3B-Tet-on and vector control cells were determined using the Atlas human cDNA expression array. Fifteen genes were differentially­ expressed between these two cell lines, among which seven genes were up-regulated and eight genes were down-regulated by BINPL-1. Furthermore, the differential expression result was confirmed by semiquantitative RT-PCR. Among these differentially expressed genes, p16INK4, IL-12, TRAIL and the lymphotoxin b gene involved in growth suppression or cell apoptosis were up-regulated, and PTEN involved in cell proliferation was down-regulated by BNIPL-1. These results suggest that BNIPL-1 might inhibit cell growth though cell cycle arrest and/or apoptotic cell death pathway(s).

 

Key words        BNIPL-1; Tet-on; Atlas human cDNA array; apoptosis

 

Bcl-2/adenovirus E1B 19 kDa interacting protein 2 like (BNIPL) has been recently identified and characterized as an apoptosis-associated protein [13]. BNIPL shares 72% homology (46% identity) with BNIP-2 [4] and has two variants, BNIPL-1 (GenBank accession No. AF193056, also called BNIP-Sa) and BNIPL-2 (accession No. AY033000). The full-length amino acid sequence of BNIPL-1 is 82 amino acids shorter at the N-terminus compared with that of BNIPL-2. Previous­ studies have shown that BNIPL-1 interacts with Bcl-2 and Cdc42GAP, and induces­ apoptosis via the BNIP-2 and Cdc42GAP homology­ (BCH) domain [3]. However, the molecular mechanism by which BNIPL-1 promotes apoptosis is not well defined so far.

Artificially regulated transgene expression systems are indispensable tools for studying the molecular mechanism of the differentiation process. Tetracycline (Tc)-regulated transgene expression systems, known as Tet-off and Tet-on systems, have been widely applied to a variety of biological­ materials, including mammalian cells [5,6]. The Tet-on system uses a reverse Tc-regulated transactivator (rtTA), which binds a Tc-responsive element (TRE) and induces transcription of the transgene in the presence of doxycycline (Dox) [7]. Microarray technology enables the monitoring of thousands of genes, and measuring the relative­ abundance of mRNA transcripts. There are a number­ of different microarray platforms available, which include cDNAs [8] or oligonucleotides [9] that are spotted­ on nylon membranes [10] or glass slides [11]. The cDNA expression array is a novel and robust method for comparing the expression of a large number of genes simultaneously and has been successfully used to study differential­ gene expression profiles [12].

In the present study, we established a Hep3B cell line with the doxycycline-regulated BNIPL-1 transgene system­ and identified differentially expressed genes regulated by BNIPL-1 in the Hep3B cell line using cDNA microarray technology.

 

 

Materials and Methods

 

Plasmid construction

 

The human BNIPL-1 full-length coding sequence was amplified by PCR using the primers containing BamHI and ClaI restriction sites and the Myc-tag sequence (Table 1). The PCR product was subcloned into the Tet-on expression vector pTRE2hyg which contains the hygromycin-resistant gene. The recombinant plasmid construct, pTRE2hyg-BNIPL-1, was confirmed by BamHI/ClaI digestion­ (New England Biolabs Inc., Beverly, USA) and DNA sequencing.

 

Cell culture

 

The human hepatocarcinoma cell line Hep3B was cultured­ in DMEM (Gibco BRL, Carlsbad, USA) containing­ 2 mM L-glutamine solution and 10% heat-inactivated Tet system-approved fetal bovine serum (Clontech, Palo Alto, USA). Cells were maintained at 37 ºC in 5% CO2. They were plated in plates or dishes, and allowed to grow until the cultures were confluent.

 

Cell transfection

 

Hep3B-Tet-on cells constructed by Wang et al. [13] were grown in 6-well plates to 70% confluence and transfected­ with pTRE2hyg-BNIPL-1 plasmid DNA using lipofectamine reagents (Invitrogen, Carlsbad, USA) under conditions recommended by the manufacturer. After 48 h, the medium was replaced with fresh DMEM containing 100 mg/L hygromycin to select the stable cell lines. Two weeks later, 43 hygromycin-resistant clones were obtained and subjected to further selection in the medium containing­ 25 mg/L hygromycin for more than three weeks, and then analyzed by Western blot using an anti-Myc monoclonal antibody (Invitrogen). BNIPL-1 expression in the cell line was induced by incubating cells with Dox, which is a Tetracycline analog (Clontech).

 

Western blot

 

Whole-cell lysate was extracted by the T-PERTM tissue protein extraction reagent containing 1% (V/V) proteinase inhibitor cocktail. Then, 10 mg of total protein from each sample was fractionated on 10% SDS-PAGE gels and electroblotted onto PROTRAN nitrocellulose membranes. After being blocked in TBST (20 mM Tris-HCl, pH 7.6, 130 mM NaCl, 0.1% Tween-20) containing 5% non-fat milk, the membranes were incubated with the primary anti-Myc antibody at 1:2000 dilution or monocolonal antibody actin pan Ab-5 (Neomarker, Westinghouse Drive, Fremont, USA) and then with horseradish peroxidase-coupled anti-mouse IgG (Santa Cruz, USA), and developed­ using SuperSignal west femto maximum sensitivity substrate­ (Promega, Madison, USA).

 

Cell proliferation assay

 

Cell viability was determined by Cell-titer 96 aqueous­ one solution cell proliferation assay kit (Promega). Aliquots of 2103 cells/well were cultured in 96-well plates with or without 2000 mg/L Dox for 24, 48, 72, 96 and 120 h. The medium was removed at the end of the incubation. Then, 100 ml of DMEM and 20 ml of Cell-titer 96 aqueous­ one solution were added to each well and incubated­ for an additional 4 h. The absorbance at 490 nm was recorded with a 96-well plate reader. Each experiment­ was performed in triplicate and repeated at least three times.

 

RNA extraction

 

Hep3B-Tet-on-BNIPL-1/20 cells were induced with or without 2000 mg/L Dox for 24 h. Cells were washed once with PBS, and total RNA samples were extracted from 1´107 cells using Trizol reagent (Invitrogen) according to the manufacturer's instructions, followed by treatment with RNase-free DNase at 37 ºC for 20 min to avoid contamination of genomic DNA. The RNA quality and concentration were assessed using agarose gel electrophoresis and spectrophotometer readings.

 

Probe preparation

 

Five micrograms of total RNA was reverse transcribed using 1 ml of CDS primer mix, 3.5 ml of [a-32P]dATP (3000 Ci/mmol; Amersham Pharmacia Biotech, Piscata­way, USA), 0.5 ml of 0.1 M DTT, 1 ml of dNTP mixture (con­taining dCTP, dGTP and dTTP), 1 ml of MMLV reverse­ transcriptase­ and 2 ml of 5´reaction buffer provided­ in the Atlas human cDNA expression array kit (Clontech). The 32P-labeled­ cDNA probes (210106 cpm) were purified using­ a NucleoSpin extraction spin column to remove the unincorporated­ nucleotides. After purification, labeled cDNAs were denatured before use.

 

Hybridization and imaging

 

Broad-scale expression profiling was performed using Atlas human cDNA expression array kit. Each array­ is spotted with duplicate cDNA fragments from 588 known genes and nine housekeeping genes or control sequences. The 588 known genes spotted on the Atlas membrane consist­ of cDNA fragments for cell-cycle control­ proteins, stress response proteins, apoptosis-associated­ proteins, DNA transcription factors, cell receptors, extracellular cell signaling and communication proteins, etc.. A complete list of the genes spotted on the array with their array positions­ and GenBank accession numbers is available at Clontech's website (http://www.clontech.com). Hybridi­zation was carried out at 68 ºC overnight in a hybridization oven. The membranes were washed thrice with 250 ml of 2´SSC containing 1% SDS and once with 250 ml of 1´SSC containing 0.1% SDS solution at 68 ºC for 30 min. Then the membranes were wrapped and exposed to X-ray film for autoradiography at -70 ºC for 48 h.

 

Statistical analysis

 

The images for the hybridized Atlas arrays were scanned with the CycloneTM storage phosphor system (Packard Bioscience Company, Meriden, USA) and saved as *.tiff files. The images were imported into the Quantarray image­ system (Packard Bioscience Company) and analyzed. Housekeeping ubiquitin and b-actin genes were selected for normalization. The adjusted intensity was equal to the mean intensity of each duplicate gene minus the background­ value. The signal ratio of each given gene between two compared array membranes was then calculated. In the present study, the ratio threshold was set to be 2.0. Only those genes that showed a 2-fold or greater increase or decrease in number were considered to be differentially expressed.

 

RT-PCR analysis

 

The expression of some prominently changed genes was further examined by RT-PCR to confirm the array results. First-strand cDNA was synthesized from total RNA using the SuperScriptTM first-strand synthesis system­ (Invitrogen). The PCRs were performed with first-strand cDNA for 26 cycles of 30 s each at 94 ºC, 26 cycles of 30 s each at 60 ºC and 26 cycles of 30 s each at 72 ºC, followed by an extension for 3 min at 72 ºC. The primers for the PCR are listed in Table 1. The RT-PCR products were subjected to electrophoresis in 2% agarose gels. To ensure that equal amounts of reverse­-transcribed cDNA were applied to the PCR reaction, the b-actin gene was also included in the PCR as a control.

 

 

Results

 

 

Construction and identification of the BNIPL-1 expressing­ vector pTRE2hyg-BNIPL-1

 

The human BNIPL-1 gene (0.8 kb) was cloned into the multicloning site of the pTRE2hyg vector (5.3 kb), designated­ as pTRE2hyg-BNIPL-1 (6.1 kb), and verified by BamHI and ClaI restriction endonuclease digestion (Fig. 1). The complete sequence of BNIPL-1 was further verified by sequencing (data not shown).

 

Expression of BNIPL-1 by doxycycline induction

 

After screening in the medium containing hygromycin for more than six weeks, a total of 43 independent hygromycin-resistant cell lines were obtained from the Hep3B-Tet-on cells transfected with pTRE2hyg-BNIPL-1. To verify the induction of BNIPL-1 protein by doxycycline­ treatment, Western blot analysis was performed with equal amounts of total cell protein from different Hep3B-Tet-on-BNIPL-1 cell lines. The clone Hep3B-Tet-on-BNIPL-1/20 showed a low background and high Dox-induced expression of BNIPL-1, and was selected as the cell line for further investigations. The amount of BNIPL-1 induced in the clone Hep3B-Tet-on-BNIPL-1/20 cells treated with Dox (2000 mg/L) was 26-fold higher than that in the cells without Dox treatment [Fig. 2(A)]. Moreover, the BNIPL-1 expression increased in a time-dependent manner­ [Fig. 2(B)].

 

Cell proliferation assay

 

To assess cell proliferation, we used the Cell-titer 96 aqueous one solution cell proliferation assay that is based on the metabolic conversion of a tetrazolium compound MTS to a colored product by living cells. The absorbance intensity of the MTS product is directly proportional to the number of viable cells in culture. As shown in Fig. 3, overexpression of BNIPL-1 significantly inhibited Hep3B cell growth (P<0.05). At 120 h, BNIPL-1 inhibited the proliferation rate by 21.7%. Each point in Fig. 3 represents the mean+/- SD of triplicate experiments.

 

Identification of differentially expressed genes in clone Hep3B-Tet-on-BNIPL-1/20 cells

 

To identify differentially expressed genes regulated by BNIPL-1, we compared the gene expression profiles of 588 genes in clone Hep3B-Tet-on-BNIPL-1/20 cells treated with or without Dox (2000 mg/L) for 24 h. No signals were visible in the blank spots and negative control spots, indicating that the Atlas human cDNA array hybridization was highly specific. The density of housekeeping genes was very similar among the samples, indicating that the results were reliable. We used ubiquitin and b-actin genes to normalize­ the intensities among arrays. The comparison results analyzed by the Quantarray image system showed that the expression of seven genes was up-regulated (Table 2) and that of eight genes was down-regulated (Table 3) following overexpression of BNIPL-1 in Hep3B cells. The genes involved in growth inhibition (p16INK4 and IL-12) or cell apoptosis [TRAIL and lymphotoxin b (LTB) gene] were up-regulated, and the genes involved in pro­liferation (PTN) were down-regulated.

 

Semiquantitative RT-PCR

 

Although differences in gene expression regulated by BNIPL-1 were observed with the cDNA array, it was essential­ to verify these results using other methods. We monitored the mRNA level of four selected genes shown to be differentially expressed in the Atlas arrays by RT-PCR. The results were consistent with the hybridization data for each of the genes monitored (Fig. 4); that is, p16INK4 and Apo-2 ligand were up-regulated, and HBNF-1 and endothelin ET-2 were down-regulated in BNIPL-1 overexpressed cells.

 

 

Discussion

 

BNIPL-1 is homologous to BNIP-2, a previously known Bcl-2 and E1B-associated protein [4]. It contains a complete BNIP-2 and Cdc42GAP homology (BCH) domain. The expression of BNIPL-1 or the complete BCH domain can cause extensive apoptosis in cells. BNIPL-1 forms a homophilic complex via a unique sequence motif within its BCH domain, and deletion of this interacting motif removes­ its pro-apoptotic effect [3]. However, the downstream signaling molecules that are involved in BNIPL-1-induced apoptosis or cell growth suppression should be examined.

In the present study, we established a cell line, Hep3B-Tet-on-BNIPL-1/20, in which the BNIPL-1 gene can be conditionally induced by Dox treatment. Dox activates the Tet-response element and significantly enhances the expression­ of the BNIPL-1 gene in Hep3B cells in a dose-dependent manner. A stable cell line with a low background and high Dox-induced expression of BNIPL-1 was obtained, which enabled us to investigate BNIPL-1-induced gene expression efficiently and rapidly.

To define gene expression profiles and compare patterns­ of expression at different stages of inducement, we can use RT-PCR, RNase protection assays or Northern blot analysis, but these methods focus on only a few genes at a time. A more promising approach for analyzing multiple genes simultaneously is the hybridization of entire cDNA populations to nucleic acid arrays, a method that has been adopted for high-throughput analysis of gene expression [12]. It allows for the rapid detection of the gene expression­ profiles of hundreds or thousands of genes simultaneously. Here, we applied cDNA array technology to compare the expression profiles in the human hepatocarcinoma cell line Hep3B regulated by BNIPL-1. Following the overexpression of BNIPL-1 in Hep3B cells, the expression­ levels of 15 genes were altered, of which seven were up-regulated and eight were down-regulated. Genes involved­ in growth inhibition (p16INK4 and IL-12) or cell apoptosis (TRAIL and LTB) were up-regulated. P16 has been reported to induce apoptosis in K562 cells [14]. IL-12 inhibits proliferation and increases the apoptotic rate of IL-12Rb2-transfected B cell lines in vitro [15]. TRAIL is a member of the TNF gene superfamily that is able to induce apoptosis through the binding of two cell-surface DRs (DR4 and DR5). TRAIL directly engages the cell-extrinsic apoptosis pathway through recruitment of the death-inducing signal complex and activation of initiator caspases, which activate effector caspases that ultimately execute programmed cell death. TRAIL can also engage the intrinsic apoptosis pathway through caspase-dependent­ activation of proapoptotic Bcl-2 family members [16]. The protein of LTB gene which is localized to the major histocompatibility complex region on chromosome 6p21.3 has an important role in the formation of germinal center reactions and regulation­ of immune response and apoptosis [17]. Genes involved in proliferation (PTN) and cell apoptosis (ET-2) were found to be down-regulated. In serum-starved NIH3T3 fibroblasts, PTN prevents apoptosis and induces cell growth. Anti-apoptotic signaling of PTN in NIH3T3 fibroblasts takes place via the MAP kinase pathway [18]. Endothelins (ET-1, ET-2 and ET-3) are a group of small vasoactive peptides (21 amino acids) with diverse paracrine/autocrine actions. ETs are expressed in a broad range of tumors; adding the ET-2 peptide to a culture medium moderately reduces the tumor cell death associated­ with hypoxia [19].

In conclusion, our results show that the combination of a Tet-on inducible system and cDNA microarray is an effective­ and powerful strategy to identify downstream signaling molecules induced by BNIPL-1. Among these differentially expressed genes induced by BNIPL-1 in hepatoma Hep3B cells, p16INK4, IL-12, TRAIL and LTB were up-regulated, and PTN was down-regulated. These results suggest that BNIPL-1 might promote cell apoptosis and growth suppression through cell cycle arrest and/or apoptotic cell death pathway(s). However, the exact mechanisms need to be investigated further.

 

 

References

 

1     Qin W, Hu J, Guo M, Xu J, Li J, Yao G, Zhou X et al. BNIPL-2, a novel homologue of BNIP-2, interacts with Bcl-2 and Cdc42GAP in apoptosis. Biochem Biophys Res Commun 2003, 308: 379385

2     Xie L, Qin WX, He XH, Shu HQ, Yao GF, Wan DF, Gu JR. Differential gene expression in human hepatocellular carcinoma Hep3B cells induced by apoptosis-related gene BNIPL-2. World J Gastroenterol 2004, 10: 12861291

3     Zhou YT, Soh UJ, Shang X, Guy GR, Low BC. The BNIP-2 and Cdc42GAP homology/Sec14p-like domain of BNIP-Sa is a novel apoptosis-inducing sequence. J Biol Chem 2002, 277: 74837492

4     Boyd JM, Malstrom S, Subramanian T, Venkatesh LK, Schaeper U, Elangovan B, Da-Eipper C et al. Adenovirus E1B 19 kDa and Bcl-2 proteins interact with a common set of cellular proteins. Cell 1994, 79: 341351

5     Lewandoski M. Conditional control of gene expression in the mouse. Nat Rev Genet 2001, 2: 743755

6     Gossen M, Bujard H. Studying gene function in eukaryotes by conditional gene inactivation. Annu Rev Genet 2002, 36: 153173

7     Gossen M, Freundlieb S, Bender G, Muller G, Hillen W, Bujard H. Transcriptional activation by tetracyclines in mammalian cells. Science 1995, 268: 17661769

8     Mirnics K, Middleton FA, Marquez A, Lewis DA, Levitt P. Molecular characterization of schizophrenia viewed by microarray analysis of gene expression in prefrontal cortex. Neuron 2000, 28: 5367

9     Yeoh EJ, Ross ME, Shurtleff SA, Williams WK, Patel D, Mahfouz R, Behm FG et al. Classification, subtype discovery, and prediction of outcome in pediatric acute lymphoblastic leukemia by gene expression profiling. Cancer Cell 2002, 1: 133143

10 Bertucci F, Bernard K, Loriod B, Chang YC, Granjeaud S, Birnbaum D, Nguyen C et al. Sensitivity issues in DNA array-based expression measurements and performance of nylon microarrays for small samples. Hum Mol Genet 1999, 8: 17151722

11   Duggan DJ, Bittner M, Chen Y, Meltzer P, Trent JM. Expression profiling using cDNA microarrays. Nat Genet 1999, 21: 1014

12   Yu Z, Guo R, Ge Y, Ma J, Guan J, Li S, Sun X et al. Gene expression profiles in different stages of mouse spermatogenic cells during spermatogenesis. Biol Reprod 2003, 69: 3747

13   Wang JR, Qin WX, Shu HQ, Pan Y, Zhang YY, Wan DF. Establishment of Hep3B cell line controlled by the Tet-on regulatory system. Tumor 2002, 22: 191193

14   Minami R, Muta K, Umemura T, Motomura S, Abe Y, Nishimura J, Nawata H. p16INK4a induces differentiation and apoptosis in erythroid lineage cells. Exp Hematol 2003, 31: 355362

15   Airoldi I, di Carlo E, Banelli B, Moserle L, Cocco C, Pezzolo A, Sorrentino C et al. The IL-12Rb2 gene functions as a tumor suppressor in human B cell malignancies. J Clin Invest 2004, 113: 16511659

16   Kelley SK, Ashkenazi A. Targeting death receptors in cancer with Apo2L/TRAIL. Curr Opin Pharmacol 2004, 4: 333339

17   Nagy B, Ferrer A, Larramendy ML, Galimberti S, Aalto Y, Casas S, Vilpo J et al. Lymphotoxin b expression is high in chronic lymphocytic leukemia but low in small lymphocytic lymphoma: A quantitative real-time reverse transcriptase polymerase chain reaction analysis. Haematologica 2003, 88: 654658

18   Bowden ET, Stoica GE, Wellstein A. Anti-apoptotic signaling of pleiotrophin through its receptor, anaplastic lymphoma kinase. J Biol Chem 2002, 277: 3586235868

19   Grimshaw MJ, Naylor S, Balkwill FR. Endothelin-2 is a hypoxia-induced autocrine survival factor for breast tumor cells. Mol Cancer Ther 2002, 1: 12731281