INTRODUCTION

Studies on the suppression of malignancy of tumor cells are important for the understanding of the mechanism of tumor suppressors and for the establishment of novel therapies of tumors. Since the last decade of the 20th century, several 3′untranslated regions of eukaryotic mRNAs were found to have in vitro tumor suppressor activity. They caused phenotypic reversion (lowering of tumorigenicity) in cultured tumor cells or malignant cells1, 2, 3, 4, 5, 6. However, the molecular mechanisms of these tumor suppressor effects remain elusive.

Several years ago, a regulatory element with tumor suppressor activity, NF-IL6 3′UTR, was found, which caused tumor suppression in a part of mouse malignant DT cells1, 2. Recently, in order to test its tumor-suppressive effect in different human tumor cell lines, we transfected it into SMMC7721 cells, a human hepatocarcinoma cell line, and resulted in markedly reduced malignancy (reversion) in significant part of the transfected cell clones. Then, we determined the changes in gene expression profile of a typical revertant clone by using the cDNA array analysis, to pave the way for elucidating possible molecular mechanisms of this tumor suppression. In this report, we will present results of this analysis.

MATERIALS AND METHODS

Cell lines and plasmids

SMMC7721, purchased from Cell Bank of Chinese Academy of Sciences, Shanghai, is a typical hepatoma cell line7. p14-6 is a recombinant pcD2 plasmid harboring a middle part of the 3′untranslated region of the human NF-IL6 cDNA4, stored in our lab (GenBank accession number M98056). pcDmr Cl 2 is a plasmid containing the same genes as p14-6, but harbors no insert.

Transfection of cell lines and malignancy tests

Stable transfection was performed by using the Calcium phosphate transfection kit (Invitrogen) according to the instruction manual. The concentration of cytotoxic antibiotic G418 (Sigma) was 2 mg / ml (valence 1 mg/ml) for all cells, including controls. G418R clones appeared in 1-2 week after screening. Cell growth curves determination and soft agar tests were performed according to our previous methods1. Nude mice tumorigenicity tests were performed with the Balb/c nu/nu nude mice, 4 -6 week age (from the Shanghai Center for Experimental Animals, Chinese Academy of Sciences), which were kept in the sterile rooms of Animal house of our Institute. 5× 106 cells per point were injected subcutaneously into the back of animals. Positive tumorigenicity was defined as the appearance of palpable tumor nodules in the injection points one month postinjection.

cDNA array analysis

The cDNA array analysis were done in CASarray Co., Shanghai. The cells were directly lysed in culture bottles with TrizolR (Invitrogen), then subjected to further mRNA purification, cDNA synthesis and hybridizations. The cDNA array membanes used contained 11,557 known genes and ESTs, and every gene had two spots. The analysis was performed twice to guarantee reproducible results, and passed strict quality control.

Semi-quantitative RT-PCR

6 genes, including 3 up-regulated and 3 down-regulated as revealed by cDNA array analysis, were chosen to compare their cDNA amounts by using semi-quantitative RT-PCR. The primers were designed using OLIGO version 6.44 software. For control, cDNA of ribosomal protein gene S29, a well-known housekeeping gene, was used. Reverse transcription of total RNAs from Cl 1, SMMC7721 and pcDmr was done, respectively, using Superscript II reverse transcriptase (Life Technologies and equal amounts (1 μl) of RT mixture were used in a 50 μl PCR reaction containing various primer pairs and FD Taq DNA polymerase (Fuhua Co., Shanghai) under the same conditions. The intensities of bands were measured using Tanon gel electrophoresis analysis software, and the values of their strengths relative to S29 were calculated by dividing their measured band intensities by the average S29 band intensities of RT PCR products of corresponding cell lines. The results were shown in a vertical square graph.

RESULTS AND DISCUSSION

Establishment of the NF-IL6 transfectants of SMMC7721 cell line

The NF-IL6 3′UTR cDNA (plasmid p14-6) was stably transfected into SMMC7721 cells by using standard calcium phosphate procedure. It was found that this hepatoma cell line has a certain endogenous G418 resistance, and thus a higher concentration (2 mg/ml) of G418 was used for screening. This drug concentration led to the reduction in the number of G418R clones; however, all the transfectants contained plasmids, as proved by Southern blots.

Phenotypic characteristics of the NF-IL6 3′UTR transfectants of SMMC7721

By high concentration G418 screening, 14 independent clones were obtained, 6 of which were found to be contaminated by mycoplasma and excluded. The remaining 8 clones were used in experiments.

These p14-6 transfectants had morphology similar to their parent cells. The majority of clones grew slower than SMMC7721 (see Fig 1B-1) and had reduced ability to form colonies in soft agar medium (Tab 1 and Fig 1A; in the figure, only Cl 1 and two controls were shown). In the in vivo tumorigenicity tests, those transfectant clones that formed less colonies in soft agar also appeared to have reduced tumorigenicity (Tab 2). Therefore, these 7 transfectants were revertants. Southern blots with 32P-radiolabeled NF-IL6 3′UTR cDNA as probe showed that the NF-IL6 3′UTR was inserted into different positions of the cellular genome (Fig 2a). Northern blots with the same probe showed that the 3′UTR was highly expressed (Fig 2B). The remaining clone, clone 8, did not revert, possibly due to some cellular reasons. In contrast, the control pcDmr Cl 1-Cl 4 transfectants (4 clones) showed tumorigenicity non-distinctive from SMMC7721 cells (Tab 2), although two cell clones (clones 1 and 3) grew slower than the former (Fig1B-2). Interestingly, the expression of endogenous NF-IL6 gene had no significant changes at both transcriptional (Fig 2B) and translational levels as shown by Western blotting (data not shown).

Figure 1
figure 1

Growth characteristics of cell clones. A, morphology (upper) and soft agar colonies (lower) of SMMC7721 (1), revertant Cl 1 (2) and control pcDmr Cl 2 (3). Bars: 10 μm. B, growth curves of all cell clones. 1, growth curves of 8 NF-IL6 3′UTR transfectants. The signs in the curves are: , Cl 1; iõ iõ, Cl 2; •, Cl 3; ió ió, Cl 4; ið, Cl 5; , Cl 6; ♦, Cl 7; iö, Cl 8. 2, growth curves of control pcDmr clones. The signs in the curves are: ió, pcDmr Cl 1; iõ, pcDmr Cl 2; , pcDmr Cl 3; , pcDmr Cl 4. In both figures the thicker curves with signs are growth curves of SMMC7721.

Table 1 Soft agar colony test of cell clones
Table 2 Nude mice tumorigenicity* of cell clones
Figure 2
figure 2

Integration status and expression of NF-IL6 3′UTR. The probe for all hybridizations was 32P-labeled 3′UTR fragment. In both figures, 1-8 stand for transfectants Cl 1-8; S, SMMC7721; I - IV, pcDmr Cl 1-4. A. Southern hybridization of cellular genomic DNAs cut by EcoRI (which had no sites on p14-6 plasmid). B. Northern hybridization of cellular total RNAs. N, endogenous NFIL6 transcript ( 2 kb). U, the 3′UTR transcript ( 0.5 kb). GAPDH, the same membrane re-probed by GAPDH probe.

Transcriptome analysis of the NF-IL6 3′UTR revertant

To obtain results as clear as possible, one typical revertant, Cl 1, was chosen for the cDNA array analysis to compare with the SMMC7721 cells and with one of the control transfectant, pcDmr Cl 2. For analysis, the genes that were consistently up- and down-regulated in Cl 1 only, but were differently regulated in SMMC7721 and pcDmr Cl 2, were considered to be induced by NFIL6 3′UTR. By this criterion, the genes, whose different expression was induced by NF-IL6 3′UTR, are shown in Tab 3. For clarity, the average expression ratios between Cl 1 and SMMC7721 cells only are shown. The semi-quantitative RT-PCR results are shown in Fig 3, which are in agreement with cDNA array analysis qualitatively (Tab 3 and Fig 3).

Table 3 Genes whose regulation is induced by NF-IL6 3′UTR, found by cDNA array analysis
Figure 3
figure 3

RT-PCR results of 3 up-regulated genes: S59049, AJ133355 and NM_020233, and 3 down-regulated genes: AF012281, M26326 and D31885, in the Cl1, SMMC7721 and pcDmr cell lines, respectively. The ordinate shows ratios of the amplified band strengths of various gene transcripts to the average band strength of the amplified S29., S59049; B, AJ133355; C, NM_020233; D, AF012281; E, M26326; F, D31885.

Among up-regulated genes, first two genes (p53-in duced ribonucleotide reductase small subunit 2 homolog, and regulator of G protein signalling 1) were directly or indirectly related to tumor suppression, and their enhanced expression favored this suppression. Among those genes that were down-regulated, the first gene (PDZ domaincontaining protein 1) is reported to be overexpressed in some epithelial tumors12; therefore, its reduced expression may contribute to the tumor suppression. The functions of other down-regulated genes include resistance to apoptosis (keratin 18), activation of signal transduction (ADP-ribosylation factor 6-like interacting protein), and activation of cell motility (thymosin β4) etc; clearly, their down-regulated expression are unfavorable for tumor growth and metastasis. Therefore, the changes in the expression of above genes, induced by the transfected NF-IL6 3′UTR, all favored to the reversion of cellular phenotype. No oncogene, that was detectably up- or down-regulated, was found by the cDNA array analysis.

We have also performed analysis on expressed proteins of these cell lines by using 2D electrophoresis (will be published elsewhere). The majority of results of 2D electrophoresis were different with cDNA array results. For example, maspin18 and annexin III19 were found up-regulated in 2D analysis. Besides technical reasons, this difference may suggest that the expression of a gene from transcription to final protein product is more complicated than that anticipated by our present knowledge.

Possible mechanism(s) of NF-IL6 3′UTR tumor suppression

It might be thought that some oncogenes were inactivated by the insertion of transfected plasmids, leading to tumor suppression. However, this does not seem to be the case, because no oncogene was found differentially and significantly expressed.

The activation of dsRNA-dependent protein kinase (PKR) by 3′UTR RNA22 was used to explain the tumor suppressor effect of α-tropomyosin; however, this explanation is not satisfactory enough, because, first, some 3′UTRs of mRNA, like the 3′UTR of NF-IL6, do not have stable double stranded structure; second, as the activation of PKR leads to general suppression of protein synthesis, the up-regulation of protein expression (like the maspin and annexin III in our case) is hard to be explained.

To sum up, we are inclined to believe that the transfection of NF-IL6 3′UTR cDNA might possibly play some regulatory role for the expression of a group of genes, either in up- or down-regulated ways. It is well known that the 3′UTRs are sites for interactions with trans-acting factors, and these interactions, as a whole, may be responsible for the maintenance of normal cell life. Therefore, the exogenous 3′UTR might also act as a regulator of gene expression. If such supposition is correct, its molecular mechanism of tumor suppression should be found in the analysis of its regulatory functions.