NDRG2 expression and mutation in human liver and pancreatic cancers

Abstract

AIM: To investigate the expression of NDRG2 and mutation of the entire coding region of NDRG2 in human liver and pancreatic cancers, and to further discuss the possible causes of NDRG2 distinct expression patterns.

METHODS: Reverse transcription-polymerase chain reaction (RT-PCR) was used to analyze the expression of NDRG2 mRNA in 37 fresh cancer specimens (including 8 cases of pancreatic cancer and 29 cases of liver cancer) and adjacent normal tissues collected from clinical operation. In addition, mutation analysis of the whole coding region of NDRG2 in these cancers was examined by polymerase chain reaction-single strand conformational polymorphism (PCR-SSCP).

RESULTS: Compared with adjacent normal tissues, the expression levels of NDRG2 mRNA in corresponding cancer tissues reduced significantly (pancreatic cancer: 0.680 ± 0.112 vs 2.089 ± 0.214, P < 0.01) (liver cancer: 0.894 ± 0.098 vs 1.345 ± 0.177, P < 0.05). Using PCR-SSCP, the mutation of the whole coding region of NDRG2 was not found in those cancer tissues where the expression of NDRG2 mRNA reduced markedly.

CONCLUSION: NDRG2 gene might express differently between normal tissues and cancer tissues, and might play an important role in the development of pancreatic cancer and liver cancer. Low expression of NDRG2 might be unrelated to the mutation of coding region of NDRG2.

Key Words: $[Keywords]

INTRODUCTION

NDRG2, a new member of the N-Myc downstream-regulated gene (NDRG) family including NDRG1-4, was first discovered and cloned from normal human whole brain cDNA library by subtractive hybridization in our laboratory[1]. NDRG2 is located on chromosome 14q11.2, including 16 exons and 15 introns, and its cDNA is 2024 bp in length[2]. Studies from our laboratory and others have shown that NDRG2 was highly expressed in many normal tissues, especially in the brain, heart, skeletal muscle, and kidney, whereas low expressed or not detected in various tumors and tumor cell lines[3,4]. This suggests that NDRG2 inactivation might play an important role in some tumor genesis or evolution, and NDRG2 might be a new cancer suppressor gene candidate. To date, no report on mutation analysis of NDRG2 in a variety of tumors has been presented. This study was conducted to determine whether the distinct expression of NDRG2 was related to gene mutation as most cancer suppressor genes. The expression of NDRG2 mRNA in 37 patients with liver cancer and pancreatic cancer was detected, and the mutation of the whole coding region of NDRG2 in normal and cancer specimens was analyzed by PCR-SSCP.

MATERIALS AND METHODS

Patients and tissue specimens

Thirty-seven fresh cancer specimens (including 8 cases of pancreatic cancer and 29 cases of liver cancer) and corresponding adjacent normal tissues were collected from patients who underwent surgery at First Affiliated Hospital of Xi’an Jiaotong University, Tangdu Hospital and Xijing Hospital between May 2003 and September 2003. All specimens were examined by microscope after haematoxylin and eosin (H&E) staining.

RNA isolation and RT-PCR

Total RNA was extracted from tissue specimens using TRIzol reagent (GIBCO BRL, USA) according to the manufacturer’s instructions. Reverse transcription reaction was set up according to Promega’s reverse transcription system protocol. cDNAs were then amplified by PCR using specific primers for human NDRG2: upstream primer 5’- GC GGA TCC ATG GCG GAG CTG CAG GAG GTC -3’; and downstream primer 5’- GC GAA TTC AAC AAG GGC CAT TCA ACA GGA GAC -3’, about 1200 bp. The primer sequences of glyceraldehyde-3-phosphate dehydrogenase (GAPDH) which was used as an internal quantitative control for the amplification were as follows: upstream primer 5’ -GCC TCA AGA TCA GCA AT- 3’; downstream primer 5’ -AGG TCC ACC ACT GAC ACG TT- 3’, 310 bp. Cycling condition for NDRG2 was at 94 °C for 5 min, followed by 32 cycles of 94 °C for 40 s, 59 °C for 50 s, and 72 °C for 60 s, and the cycling condition for GAPDH was at 94 °C for 5 min, followed by 30 cycles of 94 °C for 30 s, 57 °C for 30 s, and 72 °C for 30 s, and finally for 10 min at 72 °C. The PCR products were then analyzed and visualized on 12 g/L agarose gel containing 5 g/L ethidium bromides. Results were normalized by the ratio of band density of NDRG2 mRNA to GAPDH mRNA.

PCR-SSCP silver staining

According to cDNA sequence of NDRG2 gene, we designed by DNAstar analysis software and synthesized 6 pairs of primers covering the whole coding region of NDRG2 gene, with length of each PCR product around 150 to 370 bp, and PCR reaction conditions are described in Table 1. The PCR products were separated by electrophoresis on 12 g/L agarose gel to check the specificity of each PCR reaction. The products were mixed with the same amount of formamide-loading buffer (950 g/L formamide, 10 mmol/L EDTA, and 0.5 g/L bromophenol blue, 0.5 g/L xylene cyanol) and heated at 100 °C for 5 min, and loaded on 80 g/L non-denaturing polyacrylamide gel in 1 × TBE buffer for 8-12 h at 40 V. After electrophoresis, the gel was fixed, silver-stained, and finally developed, photographed and analyzed. The sample was considered normal if the band position was the same as that of the normal tissue. Otherwise, it was abnormal.

Table 1 Primer sequence, annealing temperature and PCR product size for PCR-SSCP.

	Primer	Sequence	Annealing temperature (°C)	PCR product size (bp)
A	Upstream	5’-ATGGCGGAGCTGCAGGAGGTC-3’	54	186
	Downstream	5’-ATTTATAGTTGAGTCCCACATC-3’
B	Upstream	5’-CACCATACGGCTCTGTCACTTTCACT-3’	58	188
	Downstream	5’-CTCTTCCATTCCAGGGGCATCCAC-3’
C	Upstream	5’-GCCAGCGATCCTTACCTACC-3’	58	365
	Downstream	5’-GGCTGCCCAATCCATCCAACC-3’
D	Upstream	5’-CCGGACACTGTTGAAGGT-3’	55	364
	Downstream	5’-GGGTCCAGTTTGAGTTACATT-3’
E	Upstream	5’-CTGTGATGCTGGTGGTAGGAGA-3’	58	286
	Downstream	5’-TGGGACAGGTGCGAGAG-3’
F	Upstream	5’-GTCCGGTCTCGTACAGCCTCTC-3’	56	153
	Downstream	5’-AACAAGGGCCATTCAACAGGAGAC-3’

Statistical analysis

Statistical analysis was performed by SPSS version 10.0. Data were expressed as mean±SD, and analyzed by Student’s t test. P < 0.05 was considered statistically significant.

RESULTS

RT-PCR data

We investigated the expression patterns of NDRG2 mRNA by RT-PCR in 37 cases of human cancers including 29 liver cancers and 8 pancreatic cancers. Figure 1A shows PCR products of representative agarose gels for NDRG2 and GAPDH in human liver cancer and pancreatic cancer. The size of the RT-PCR products was as predicted, 1200 bp for NDRG2 and 310 bp for GAPDH. The levels of NDRG2 mRNA expression were quantified by the ratio of band optical density of NDRG2/ GAPDH (Figure 1B). Compared with adjacent normal tissues, there were marked decreases in NDRG2 mRNA expression in liver cancer (0.894 ± 0.098 vs 1.345 ± 0.177, P < 0.05) and pancreatic cancer (0.680 ± 0.112 vs 2.089 ± 0.214, P < 0.01).

Figure 1 Expressions of NDRG2 mRNA in human liver cancer and pancreatic cancer. A: Representative agarose gels of RT-PCR-amplified NDRG2. M: DNA molecular weight marker; N: Normal tissue; C: Cancer tissue; B: Relative amount of NDRG2 mRNA was evaluated by the ratio of band optical density of NDRG2/GAPDH and the expression was significantly decreased in cancer (n = 29 for liver cancer, and 8 for pancreatic cancer respectively) vs normal tissues. ^aP<0.05, ^bP < 0.01 vs normal.

PCR-SSCP analysis

The PCR products, which were amplified by 6 pairs of primers covering the whole coding region of NDRG2 gene, were separated by electrophoresis on 12 g/L agarose gel to check the specificity of each PCR reaction (Figure 2). PCR products with only one specific band of right size in the gel can further be used for SSCP assay. Figure 3 shows representative results of SSCP silver staining of the PCR products amplified by 6 pairs of primers covering the whole coding region of NDRG2. By SSCP, all the PCR products from the cancer tissues of 27 patients of NDRG2 low-expression confirmed by RT-PCR, did not show any mutation.

Figure 2 Agarose gel analysis of PCR products amplified by 6 pairs of primers covering the whole coding region of NDRG2. M: DNA molecular weight marker, PCR product size was 186 bp for A, 188 bp for B, 365 bp for C, 364 bp for D, 286 bp for E and 153 bp for F.

Figure 3 Representative results of SSCP silver staining of the PCR products (186 bp for A, 188 bp for B, 365 bp for C, 364 bp for D, 286 bp for E, and 153 bp for F) amplified by 6 pairs of primers covering the whole coding region of NDRG2. C: cancer tissue; N: normal tissue.

DISCUSSION

NDRG2, together with NDRG1, NDRG3 and NDRG4, constitute NDRG gene family. The four members, which share 53%-65% amino acid identity, belong to the alpha/beta hydrolase super family[5]. Phylogenetic analysis of the family demonstrated that human NDRG1 and NDRG3 belong to a subfamily, and NDRG2 and NDRG4 to another. Expression of the fusion proteins showed that all of them were cytosolic proteins. Northern and dot blot analysis showed that NDRG gene family was evidently expressed in adult brain and almost not detected in some human cancer lines[6-8]. In addition, NDRG2 is highly expressed in adult skeletal muscle and brain, NDRG3 remarkably expresses in brain and testis, and NDRG4 evidently expresses in brain and heart, suggesting that they might display different specific functions in distinct tissues[9,10]. Although NDRG gene family is believed to be involved in cell differentiation, the exact molecular and cellular function of members of this family is still unknown.

In 2000, our laboratory first discovered and successfully cloned human NDRG2 cDNA from normal human whole brain cDNA library by subtractive hybridization. Our previous studies showed that the expression level of NDRG2 mRNA was very high in brain, salivary gland, skeletal muscle and mammary gland, and low in bone marrow, testis, peripheral blood and placenta, and not detectable in leukocyte, colon and some tumor cell lines. Other studies have also shown similar result[11,12], suggesting that NDRG2 inactivation might play an important role in some cancer genesis or evolution, and NDRG2 might be a new cancer suppressor gene candidate.

A recent study showed that over expression of NDRG1 in tumor cells decreased the proliferation rate, enhanced differentiation, and suppressed metastasis potency of cancer cells. More interestingly, a nonsense mutation in NDRG1 gene may be the cause of hereditary motor and sensory neuropathy-Lom (HMSNL)[13].

The activation of oncogene and inactivation of cancer suppressor gene can lead to uncontrolled cell proliferation and occurrence of cancer. The inactivation of most cancer suppressor genes has been confirmed to relate to gene mutation such as p53, p16 and Rb gene[14]. To the authors’ knowledge, we first investigated the expression patterns of NDRG2 mRNA by RT- PCR in 37 cases of human cancers including 29 liver cancers and 8 pancreatic cancers. The results showed that there was marked decrease in NDRG2 mRNA expression in human liver cancer (P < 0.05) and pancreatic cancer (P < 0.01), compared with adjacent normal tissue. We speculated that different expression pattern of NDRG2 might relate to gene mutation. The designed 6 pairs of primers covering the whole coding region of NDRG2 gene to amplify the cancer tissue of 27 patients of NDRG2 low-expression were confirmed by RT-PCR, but no mutation was found.

In conclusion, our results indicate that NDRG2 gene expresses differently between normal tissues and cancer tissues, and it might play an important role in the development of pancreatic cancer and liver cancer. The low expression of NDRG2 might be unrelated to mutation of the coding region of NDRG2. In addition, whether a mutation exists in the upstream sequence of transcription initial site of NDRG2, deserves further investigation.

ACKNOWLEDGEMENTS

We are grateful to Dai-Hua Yu, Bang-Yong Sun, Wei Lin, and Jiang-Wei Liu for the collection of tissue specimens and Jiang-Tian Yu, Li-Hu Wu, and Li-Feng Wang for digital photographs.