- Research article
- Open Access
Simultaneous down-regulation of tumor suppressor genes RBSP3/CTDSPL, NPRL2/G21 and RASSF1A in primary non-small cell lung cancer
BMC Cancer volume 10, Article number: 75 (2010)
The short arm of human chromosome 3 is involved in the development of many cancers including lung cancer. Three bona fide lung cancer tumor suppressor genes namely RBSP3 (AP20 region),NPRL2 and RASSF1A (LUCA region) were identified in the 3p21.3 region. We have shown previously that homozygous deletions in AP20 and LUCA sub-regions often occurred in the same tumor (P < 10-6).
We estimated the quantity of RBSP3, NPRL2, RASSF1A, GAPDH, RPN1 mRNA and RBSP3 DNA copy number in 59 primary non-small cell lung cancers, including 41 squamous cell and 18 adenocarcinomas by real-time reverse transcription-polymerase chain reaction based on TaqMan technology and relative quantification.
We evaluated the relationship between mRNA level and clinicopathologic characteristics in non-small cell lung cancer. A significant expression decrease (≥2) was found for all three genes early in tumor development: in 85% of cases for RBSP3; 73% for NPRL2 and 67% for RASSF1A (P < 0.001), more strongly pronounced in squamous cell than in adenocarcinomas. Strong suppression of both, NPRL2 and RBSP3 was seen in 100% of cases already at Stage I of squamous cell carcinomas. Deregulation of RASSF1A correlated with tumor progression of squamous cell (P = 0.196) and adenocarcinomas (P < 0.05). Most likely, genetic and epigenetic mechanisms might be responsible for transcriptional inactivation of RBSP3 in non-small cell lung cancers as promoter methylation of RBSP3 according to NotI microarrays data was detected in 80% of squamous cell and in 38% of adenocarcinomas. With NotI microarrays we tested how often LUCA (NPRL2, RASSF1A) and AP20 (RBSP3) regions were deleted or methylated in the same tumor sample and found that this occured in 39% of all studied samples (P < 0.05).
Our data support the hypothesis that these TSG are involved in tumorigenesis of NSCLC. Both genetic and epigenetic mechanisms contribute to down-regulation of these three genes representing two tumor suppressor clusters in 3p21.3. Most importantly expression of RBSP3, NPRL2 and RASSF1A was simultaneously decreased in the same sample of primary NSCLC: in 39% of cases all these three genes showed reduced expression (P < 0.05).
The short arm of human chromosome 3 (region 3p21.3) contains clusters of tumor suppressor genes (TSG) involved in multiple cancer types including lung, kidney, breast, cervical, nasopharyngeal and other carcinomas [1–5]. We performed a comprehensive deletion survey of 3p in more than 400 of major epithelial cancer samples and identified two most frequently affected regions - LUCA at the centromeric and AP20 at the telomeric border of 3p21.3 [6–10]. Aberrations in these loci were detected in more than 90% of studied tumors. Homozygous deletions (HD) were frequently detected in all tumors in both the LUCA and AP20 regions. The frequent chromosome losses in these regions suggested that they harbor several multiple TSG [9, 10]. More than 20 genes were localized in these two regions and among them at least three TSG were identified: RBSP3 (RB protein serine phosphatase from chromosome 3 gene or HYA22 or CTDSPL; CTD small phosphatase family), NPRL2 (nitrogen permease regulator-like 2 gene or G21 or TUSC4; NPR family) and RASSF1A (Ras association domain family member 1 gene). RBSP3 was mapped to AP20 and the others two to the LUCA region [1, 11–13].
The RBSP3 gene occupies more than 120 kb and contains at least 8 exons coding for a 4.8 kb mRNA that is ubiquitously expressed in normal tissues including lung. By sequence analysis RBSP3 belongs to a gene family of small C-terminal domain phosphatases that may control the RNA polymerase II transcription machinery . Two sequence splice variants of RBSP3 (A and B) were identified and an initial analysis of RBSP3 was performed in lung and other cancers . The expression of the gene was greatly decreased in several small cell lung cancer (SCLC) and NSCLC cell lines. RBSP3 showed growth suppression with regulated transgenes in cell culture and suppression of tumor formation in SCID mice. It was demonstrated that transient expression of variant A and B resulted in drastic reduction of the phosphorylated form of RB protein presumably leading to a block of the cell cycle at the G1/S boundary. In addition, frameshift, nonsense and missense mutations in RBSP3 have been reported . All these features are consistent with classical characteristics of a TSG.
The NPRL2/G21 gene covers 3.3 kb and contains 11 exons coding for the main 1.8 kb transcript with multiple splice isoforms that are expressed in all tested normal tissues including lung. By sequence analysis, the main product of NPRL2/G21 encodes a soluble protein that has a bipartite nuclear localization signal, a protein-binding domain, similarity to MutS core domain, and a newly identified nitrogen permease regulator 2 domain with unknown function. This information suggests that the nuclear protein NPRL2/G21 may be involved in DNA mismatch repair, cell cycle checkpoint signaling, and regulation of the apoptotic pathway. NPRL2 plays an important role in cisplatin-induced resistance in human non-small-cell lung cancer cells [16, 17]. Previously obtained results indicated that NPRL2/G21 is a multiple tumor suppressor gene [16, 18, 19].
The RASSF1 gene occupies 7.6 kb and contains 5 exons coding for 2 kb alternatively spliced mRNAs [6, 11, 20]. One of the major splicing forms is RASSF1A. Several studies have shown that loss of RASSF1A expression occurs in many different cancers because of tumor acquired promoter DNA methylation and the gene is able to suppress growth of lung cancer cells in culture and tumor formation in mice [13, 21–24]. For example, RASSF1A is silenced by promoter hypermethylation in 100% of SCLCs and in 63% of NSCLCs cell lines and in 50-100% SCLC and 21-58% NSCLC primary tumors [25–28]. As in the case of RBSP3, frameshift, nonsense and missense mutations in RASSF1A have been discovered [15, 29]. The amino acid sequence of RASSF1A (340 amino acids) contains a predicted diacylglycerol (DAG) binding domain and a Ras association domain. Association of human proteins RASSF1C and RASSF1A with Ras protein was demonstrated [30, 31]. RASSF1A can induce cell-cycle arrest by engaging the Rb-family cell cycle checkpoint . RASSF1A is involved in several growth regulating and apoptotic pathways and regulates cell proliferation, cellular integrity and cell death [24, 27]. These and other results strongly suggest that RASSF1A is an important human TSG involved in the development or progression of many epithelial tumors.
Previously only few studies were performed to compare expression of several 3p TSG in the same tumor sample [33, 34]. To investigate this further we chose RBSP3, NPRL2 and RASSF1A and analyzed their expression by qPCR in primary tumors: non-small cell lung cancer (NSCLC) - adenocarcinoma (AC) and squamous cell lung cancer (SCC).
For the first time we found that expression of all three genes was significantly decreased in 67-85% of tested NSCLC cases. Moreover, the simultaneous down-regulation of RBSP3, NPRL2 and RASSF1A in the same tumor sample was observed in 39% of all cases. Both genetic and epigenetic mechanisms contributed to deregulation of these three genes representing two TSG clusters in 3p21.3.
Paired specimens of non-small cell lung cancer (NSCLC) tissues including 41 squamous cell carcinomas (SCC), 18 adenocarcinomas AC) and adjacent morphologically normal tissues (conventional "normal" matched control samples) were obtained after surgical resection of primary lung cancer prior radiation or chemotherapy and stored in liquid nitrogen. "Normal" matched controls were obtained minimum at 2 cm distance from the tumor and confirmed histologically as normal lung epithelial cells. The diagnosis was verified by histopathology and only samples containing 70% or more tumor cells were used in the study. The samples were collected in accordance to the guidelines issued by the Ethics Committee of Blokhin Cancer Research Center, Russian Academy of Medical Sciences (Moscow). All patients gave written informed consent that is available upon request. The study was done in accordance with the principles outlined in the Declaration of Helsinki. All tumor specimens were characterized according to the International System of Clinico-Morphological Classification of Tumors (TNM), based on the tumor-node-metastasis and staging classification of 1989  and WHO criteria classification of 1999 . Relevant clinical and pathological characteristics of the patients with NSCLC included in this study are summarized in Table 1. Normal lung tissues (autopsy material) were obtained post mortally from ten healthy individuals (age 23-49 lacking cancer history with absence of chronic diseases).
DNA and total RNA extraction and reverse transcription reaction
DNA was extracted using the Dneasy Tissue kit (Qiagen, USA) and total RNA was isolated with Rneasy mini kit (Qiagen, USA) according to the manufacturer's recommendation. RNA quality was assessed with spectrophotometer NanoDrop ND-1000 (NanoDrop Technologies Inc. USA) and by gel electrophoresis. All RNA samples were treated with RNAse free DNase I (Fermentas, Lithuania) and cDNA was synthesized using MMLV reverse transcriptase and random hexamers according to standard manufacturer's protocol (Fermentas, Lithuania).
Analysis of mRNA and DNA copy number by qPCR
The sequences of primers and probes are shown in Table 2. All reactions were performed using ABI 7000 PRISM™ SDS (Applied Biosystems) with RQ software (PCR program: 10 min at 95°C, then 40 two-step cycles 15 s at 95°C and 60 s at 60°C) in total volume 25 μl in triplicate. All probes contained the dye FAM at 5'-end and RTQ1 at 3'-end. Final concentrations of primers and probes for target and reference genes were: RBSP3 cDNA primers- 350 nM, probe - 150 nM; RBSP3 DNA primers - 150 nM, probe - 100 nM; ACTB DNA primers- 200 nM, probe - 100 nM, NPRL2 cDNA primers - 500 nM, probe 300 nM; RASSF1A cDNA primers - 300 nM, probe 300 nM;GAPDH cDNA primers - 300 nM, probe - 150 nM; RPN1 cDNA primers - 350 nM, probe - 200 nM;RPN1 DNA primers - 200 nM, probe - 100 nM; GUSB cDNA primers - 350 nM, probe - 250 nM; GUSB DNA primers - 200 nM, probe - 200 nM. PCR products were analyzed in 1.8% agarose gels and nucleotide sequences of the amplicons were verified by sequencing with 3730 DNA Analyzer automated sequencer (Applied Biosystems).
qPCR data were analyzed using the relative quantification or ΔΔCt-method [37, 38] based on mRNA (or DNA) copy number ratio (R) of a target gene versus reference gene in a given tumor sample relative to matched normal control sample (see above Tissue specimens section) according to the formula:
where E - efficiency of reaction, C T - threshold cycle, ref - reference gene, tar - target gene.
All preliminary validation steps have been done: standardization of all assays, reproducibility of the qPCRs in parallel and in independent runs, selection of reference samples and testing of reference genes http://www.gene-quantification.info/.
and NotY: 5'-GGCCGCGAATTCGTCGGTATGCACTGTGTGTGACATTCAAA-3"
were used to create the NotI linker. Annealing was carried out in a final volume of 100 μl containing 20 μl of 100 μM NotX, 20 μl of 100 μM NotY, 10 μl of 10×M buffer (Roche Molecular Biochemicals) and 50 μl of H2O. Two micrograms of tumor and normal control DNA (50 μg/ml) were digested with 20 U of Sau3A (Roche Molecular Biochemicals) at 37°C for 5 h and then 0.4 μg of the digested DNAs were circularized overnight with the T4 DNA ligase (Roche Molecular Biochemicals) in the appropriate buffer in 1 ml reaction mixture. Then DNA was concentrated with ethanol, partially filled in and digested with 10 U of NotI at 37°C for 3 h. Following digestion, NotI was heat inactivated and DNAs were ligated overnight in the presence of a 50 M excess of NotI linker at room temperature. NotI- representation (NR) probes were labeled in a PCR reaction with NotX primer. The majority of products of the DNA amplification step were in the 0.2-1.0 kb range. Repeated PCR was conducted for labeling NR with fluorophores.
Hybridization of coupled normal/tumor NotI samples was carried out at 42°C for 15 h in a Lucidea Base device (Amersham Pharmacia Biotech) according to manufacturer's recommendations. Automatic washing of the microarrays was performed in the same device using manufacturer's protocol. The following solutions were sequentially used for the washing: 1) 0.2% SDS+1 SSC; 2) 0.2% SDS+0.1 SSC; 3) 0.1 SSC; 4) de-ionized water; 5) isopropyl alcohol. Then microarrays were scanned in the GenePix 4000 A and results were processed with GenePix Pro 6.0 software (Amersham Pharmacia Biotech).
Nonparametric Wilcoxon test was used to compare mRNA expression differences of target and reference genes for the same NSCLC sample. Then groups of samples were compared in respect to average level of mRNA decrease (LDav) and the frequency of decrease (FD). The LD was calculated as 1/R and reflects the n-fold factor by which the mRNA content decreased in the tumor compared to normal tissue. Nonparametric Kruskal-Wallis and Mann-Whitney rank-sum tests were used to test mRNA differences (both LDav and FD) for each target gene in NSCLC (AC, SCC) and with and without metastases. Nonparametric Spearmen's criterion was used to calculate the coefficient of correlation between the levels of mRNA decrease (LDav) for each set of pairs of target genes. P-values < 0.05 were considered statistically significant. All statistical procedures were performed using the BioStat software .
Expression of RBSP3, NPRL2 and RASSF1A genes in lung tissues of healthy donors
Genomic structure and location of primers and probes is shown in Figure 1.
It is known that normally looking cells surrounding tumors can already harbor genetic and epigenetic changes compared to the cells isolated from the same tissue from healthy individuals. Therefore we decided to check expression of these 3 genes in samples obtained from healthy donors and in biopsies of "normal" (matched) control samples obtained from patients with NSCLC. The mRNA levels were measured in lung tissues from healthy donors using different reference samples. RNA pool isolated from lung tissues of 10 healthy donors was used in one set of experiments and in another test the same samples were examined against RNA pool isolated from adjacent to tumor "normal" tissues of 30 NSCLC patients. The GAPDH was used as a reference gene and its expression did not vary more than 2-fold in all samples examined. Very similar results were obtained for target genes in both experiments. We concluded that expression of all three genes did not significantly vary (R = 1.5 ± 0.5) in normal lung tissues from healthy donors and patients with NSCLC ("normal" matched controls).
Expression of three TSGs in primary NSCLC - AC and SCC
The results of mRNA level quantification for three genes and statistical analysis are shown in Figures 2 and 3 and Table 3. The variability of GAPDH and RPN1 mRNA was less than 2-fold in studied samples, therefore in this study we considered significant greater than or equal to 2-fold decrease of target genes expression.
Significant (from 2 to 94-fold) decrease of the mRNA level (LD) was observed in 85% (50 of 59, P < 0.001) of all NSCLC cases. Frequency of mRNA decrease (FD) was high in both AC (78%) and SCC (88%). Interestingly at Stage I of the SCC cases expression of RBSP3 was down-regulated already in 100% cases and only in 71% of AC samples (Table 3). On the other hand the tendency of more severe deregulation of RBSP3 transcription during tumor progression was more evident for AC cases: in average a 3-fold declined expression in 70% of cases without metastases and a 6-fold in 88% of samples with metastases. The same was correct for cases at Stage I (LDav = 3, FD = 71%) and Stage III (LDav = 6, FD = 100%). Spearmen's coefficient of reverse correlation (rs) between tumor Stages (I, II, III) and expression R values was 0.36 (P = 0.2, i.e. non-significant).
Significant decrease of LD (from 2 to 100-fold) was seen in 73% (43 of 59, P < 0.001) of all NSCLC cases. Down-regulation of NPRL2 was more significant in SCC compared to AC: 9-fold and 85% vs. 4-fold and 44%, respectively, P = 0.002 (see Table 3). No expression was found in 12% (5 of 59) of SCC cases. The FD was comparable in SCC cases with metastases and without (82% and 88%, respectively). However FD in AC samples without metastases was less than with metastases (20% against 75%; P = 0.08). Again at Stage I of the SCC cases expression of NPRL2 was down-regulated already in 100% cases and in AC only 14% samples showed decreased expression. The tendency of more frequent decrease of NPRL2 mRNA level with tumor progression was seen only for AC cases, where Spearmen's coefficient of reverse correlation (rs) between tumor Stages (I, II, III) and expression R values was 0.47 (P = 0.06). Eighty three percent of AC cases had decreased expression at Stage III as compared to 14% at Stage I.
Significant decrease of LD (from 2-fold to 57-fold) was detected in 67% (24 of 36, P < 0.01) of all NSCLC cases, including 64% of AC and 68% of SCC samples (Figure 2, Table 3). More severe deregulation of RASSF1A expression was observed during tumor progression both in AC and SCC. Spearmen's coefficient of reverse correlation (rs) between tumor Stages (I, II, III) and R values was 0.54 (P = 0.05) for AC. Statistically significant decrease of LDav and increase of FD was observed in AC for patients with metastases and without (7-fold and 100% vs. 4-fold and 29%, P < 0.05). The same tendency was seen for LDav and FD in SCC cases with and without metastases: namely, 6-fold and 81% vs. 4-fold and 33%, respectively; P = 0.196. Moreover an average decrease of RASSF1A mRNA (LDav) in SCC at Stage I was 4-fold and at Stage III it was already 7-fold.
FD increased for both SCC (33% vs. 82%) and AC (29% vs. 100%) during tumor progression from Stage I to Stage III.
Surprisingly a 3-to-4 fold increase of mRNA level of two genes - NPRL2 and RASSF1A was found for several AC samples at Stage I (Figure 2, samples 1, 3, 5). These tumor samples were highly differentiated in contrast to samples 2, 4, 9, 10 and others that displayed low degree of differentiation and significant decrease of NPRL2 and RASSF1A. Some of them showed so high degree of anaplasia that the epithelial phenotype was almost undetectable. These undifferentiated cells were highly malignant and showed the strongest decrease of NPRL2 and RASSF1A expression.
Simultaneous down-regulated expression of three genes in the same tumor samples
As expression of all three target genes was frequently decreased in the same NSCLC samples we tested the probability that this decrease was not random. To do this we used 36 NSCLC samples and two control reference genes GAPDH and RPN1 (Figure 3, Table 4). Results demonstrated that this simultaneous down-regulation was very frequent and statistically valid: all three genes had reduced expression (≥2) in 39% (P < 0.001) of NSCLC specimens (36% in AC and 41% in SCC). Simultaneous decreased expression of RBSP3 and NPRL2 was observed in 61% of cases and for RBSP3 and RASSF1A this occurred in 50%. In the case of NPRL2 and RASSF1A simultaneous down-regulation was detected in 44% of tumors (Table 4). Spearmen's coefficient values rs for RBSP3, NPRL2/G21, RASSF1A gene pairs were high (0.63 - 1.0, P < 0.001) in different groups especially in AC and SCC with metastases.
Decreased expression of RBSP3, NPRL2 and RASSF1A in primary NSCLC can be caused by genetic and epigenetic factors
To understand the mechanism underlying the observed down-regulation of RBSP3 in NSCLC samples we used qPCR to test DNA copy number changes, i.e. genetic factors [9, 10] and NotI microarrays to examine methylation i.e. epigenetic factors . Methylation was detected in 38% of AC (3 of 8) and in 80% of SCC (8 of 10) cases with decreased expression of RBSP3. Deletions were detected in 25% of AC samples and in 30% of SCC tumors (Table 5). Both methylation and hemizygous deletions were observed in two SCC specimens. Methylation or/and deletions of RBSP3 were detected in 63% of AC and in 90% of SCC cases. The data suggested that both genetic and epigenetic mechanisms are important for transcriptional inactivation of RBSP3 in NSCLC.
Using NotI microarrays we cannot test methylation of NPRL2 and RASSF1A promoters because they don't contain NotI sites. However we can check genes SEMA3F and GNAI2 that are located in the same LUCA sub-region as NPRL2 and RASSF1A. We tested how often deletions and methylations occurred in the same NSCLC sample at both LUCA and AP20 sub-regions and found that this occurred in 58% of all studied cases (15 of 26). Thus, most likely both genetic and epigenetic mechanisms are responsible for simultaneous down-regulation of expression of these three genes.
Several candidate TSG from the 3p21.3 AP20 and LUCA sub-regions were examined in the gene inactivation test, GIT [2, 12, 13, 16, 34, 42, 43]. The test is based on the functional inactivation of analyzed genes that can be achieved in different ways: by mutation, deletion, methylation etc. According to these results, at least three genes can now be considered as bona fide lung TSG: NPRL2 and RASSF1A from LUCA and RBSP3 from AP20 sub-regions [12, 13, 16].
Earlier the decrease of RBSP3 expression was shown in SCLC, NSCLC, cervical, renal, breast, ovary, leukemia cell lines and primary tumors by Northern blot analysis, RT-PCR and qPCR [12, 44–46]. The decrease or absence of NPRL2/G21 expression was detected in some SCLC, NSCLC and renal cancer cell lines using Northern blot analysis [6, 16]. The decrease or absence of RASSF1A expression was found in SCLC, NSCLC and many other tumors and cancer cell lines [see [24, 27, 28]].
It was reported that promoter methylation was the main mechanism of RASSF1A loss of expression in lung cancer (see Introduction). Homozygous deletion of 3'-part of NPRL2 gene and rare mutations were found in NSCLC and SCLC cell lines [6, 16]. There are no methylation data explaining the loss of RBSP3 expression in lung cancer. However frequent deletions and mutations were reported [10, 12, 15]. In some leukemia cell lines (up to 98%) and acute leukemia lymphoma blood samples (24%) methylation of the promoter region of RBSP3 was reported . Methylation (up to 26%), deletions and decreased expression of RBSP3 were significantly associated with poor prognosis of cervical cancer . Thus, inactivation of RBSP3 might be one of the early events in cervical carcinogenesis.
Loss of heterozygosity and quantitative real-time PCR demonstrated that aberrations in both LUCA and AP20 sub-regions occurred simultaneously in the same tumor with high probability. Thus, it was suggested that aberrations in both LUCA and AP20 sub-regions could be linked [9, 10]. Indeed, homozygous deletions in both regions often occur in the same tumor (P < 3 × 10-7). The estimation of possible interdependency between all aberrations in the loci NLJ-003 (AP20) and NL3-001 (LUCA) as different events was carried out using a permutation test for four types of cancers: lung, renal, breast and ovarian. This test also revealed a significant correlation between different aberrations in these two loci (P < 10-6). The same results were obtained using Pearson correlation for numeric values of copy number changes of these loci. Indeed, proteins RBSP3 and RASSF1A could collaborate in cell cycle arrest: RASSF1A by inhibiting cyclin D1  and RBSP3 by dephosphorylating pRB . Thus functional collaboration of these two genes could result in activation of the RB1 gene.
In this study we tested the hypothesis that TSG in AP20 and LUCA regions were not only deleted but their expression could also be simultaneously down-regulated in NSCLC. This suggestion was indirectly supported by other studies that showed that genes over large chromosomal regions could be regulated in a coordinated fashion [33, 47–49].
First we found that expression of all three genes is rather uniform in lung samples isolated from healthy donors and from normally looking lung samples obtained from NSCLC patients ("normal" matched control samples). Thus adjacent morphologically normal tissues from the patients can be used as paired reference controls to tumor samples. In the study two parameters were analyzed - the level of mRNA decrease and frequency of mRNA decrease in two major NSCLC histological subtypes (AC and SCC) and their subgroups with different characteristics such as clinical Stage, grade, tumor localization, presence of metastases and others. Although both parameters reflect deregulation of gene expression, they are not randomly but rather functionally related.
Expression analysis of the three genes revealed the following main features.
1. Expression of the three studied TSG was significantly decreased in NSCLC: 85% for RBSP3, 67% for RASSF1A and 73% for NPRL2 (P < 0.001). It was statistically valid both for SCC and AC.
2. Down-regulation of the three genes was already evident at Stage I of NSCLC samples. Statistically significant down-regulation of both NPRL2 and RBSP3 was seen in 100% cases at Stage I of SCC.
3. The degree and frequency of the expression decrease for all three genes was more strongly pronounced in SCC than in AC samples (see Table 3). This difference was statistically valid in the case of NPRL2 (P = 0.002).
4. All studied genes were involved in progression of AC. The tendency of more severe expression down-regulation of the RBSP3 was evident during tumor progression of AC with respect to FD and LD (70% and 3-fold in cases without metastases in contrast to 88% and 6-fold decrease in cases with metastases, P = 0.13). For NPRL2, this tendency was also seen only for AC - 83% of cases had decreased expression at Stage III compared to 14% at Stage I and in 75% of AC cases with metastases vs. 20% of cases with metastases, P = 0.08, see Table 3). Expression of RASSF1A revealed the most strongly pronounced correlation between decrease of expression (FD and LD) and tumor progression both in SCC and AC. For example, difference in FD values was obvious between cases with and without metastases. For SCC cases this difference was 33% vs. 81% (P = 0.196) and for AC it was even more sharp, 29% vs. 100% (P < 0.05).
5. Expression of RBSP3 and RASSF1A was most seriously affected in respect to FD. For RBSPS3 it was detected in 85% of all NSCLC cases and for RASSF1A in 67%. However, regarding LD, the expression of NPRL2 was most strongly inhibited (LDav = 8), while RBSP3 showed weaker inhibition (LDav = 5).
6. Preliminary data suggested that no statistically significant difference was observed in cases with relation to age, smoking history and other cytological and pathological characteristics.
7. NotI microarrays and qPCR on genomic DNA we tested for possible mechanisms of the declined expression of RBSP3 in NSCLC. The data suggested that both genetic and epigenetic mechanisms were important for transcriptional inactivation of RBSP3 in NSCLC. Altogether deletions were detected in 25% of AC samples and in 30% of SCC patients. Methylation of RBSP3 was detected in 38% of AC and in 80% of SCC cases. With NotI microarrays we also tested how often LUCA (NPRL2 and RASSF1A) and AP20 (RBSP3) regions were deleted or methylated in the same tumor and found that this occurred in 58% of all studied cases (18 of 26). Thus, most likely both genetic and epigenetic mechanisms are responsible for simultaneous down-regulation of expression of these three TSG.
The detailed analysis of mRNA expression levels of three 3p TSG was performed in two histological subtypes of NSCLC - AC and SCC respectively. The most important finding was that expression of RBSP3, NPRL2 and RASSF1A decreased in the same samples of primary NSCLC: all 3 genes have reduced expression in 39% of cases (P < 0.05). Declined expression of RBSP3 and NPRL2 was observed in 61% of samples (P < 0.05), RBSP3 and RASSF1A in 50% (P < 0.05), and NPRL2 and RASSF1A in 44% (P < 0.05) of cases. RASSF1A and NPRL2 are located in the same locus and if close localization is the reason for simultaneous decrease of expression then such decrease should be significantly less frequent when genes from the LUCA region are compared with RBSP3. However this was not the case and thus close localization is not the main factor for simultaneous decrease of expression. At present expression of many thousands of genes in NSCLC is investigated using cDNA microarrays [50, 51]. These experiments already produced very valuable data. However, results of different studies varied significantly. They should be proved by independent methods and qPCR is one of the methods of choice that allows more detailed analysis of particular genes. This work is one of such studies. All presented data confirm and extend our previous results for lung, renal, breast, cervical cancers [9, 10] and support the hypothesis that two TSG clusters (in AP20 and LUCA) are very likely co-regulated by common mechanisms.
quantitative real-time reverse transcription-polymerase chain reaction
non-small cell lung cancer
squamous cell carcinoma
RB1 serine phosphatase from human chromosome 3
nitrogen permease regulator-like 2 gene
Ras association domain family 1A
Zabarovsky ER, Lerman MI, Minna JD: Tumor suppressor genes on chromosome 3p involved in the pathogenesis of lung and other cancers. Oncogene. 2002, 21: 6915-6935. 10.1038/sj.onc.1205835.
Imreh S, Klein G, Zabarovsky ER: Search for unknown tumor-antagonizing genes. Genes Chromosomes Cancer. 2003, 38: 307-321. 10.1002/gcc.10271.
Ji L, Minna JD, Roth JA: 3p21.3 tumor suppressor cluster: prospects for translational applications. Future Oncol. 2005, 1: 79-92. 10.1517/147966184.108.40.206.
Klein G, Imreh S, Zabarovsky ER: Why do we not all die of cancer at an early age?. Adv Cancer Res. 2007, 98: 1-16. full_text.
Angeloni D: Molecular analysis of deletions in human chromosome 3p21 and the role of resident cancer genes in disease. Brief Funct Genomic Proteomic. 2007, 6: 19-39. 10.1093/bfgp/elm007.
Lerman MI, Minna JD: The 630-kb lung cancer homozygous deletion region on human chromosome 3p21.3: identification and evaluation of the resident candidate tumor suppressor genes. The International Lung Cancer Chromosome 3p21.3 Tumor Suppressor Gene Consortium. Cancer Res. 2000, 60: 6116-6133.
Braga E, Senchenko V, Bazov I, Loginov W, Liu J, Ermilova V, Kazubskaya T, Garkavtseva R, Mazurenko N, Kisseljov F, Lerman MI, Klein G, Kisselev L, Zabarovsky ER: Critical tumor suppressor gene regions on chromosome 3 regions in major human epithelial malignancies: allelotyping and quantitative real time PCR. Int J Cancer. 2002, 100: 534-541. 10.1002/ijc.10511.
Protopopov A, Kashuba V, Zabarovska VI, Muravenko OV, Lerman MI, Klein G, Zabarovsky ER: An integrated physical and gene map of the 3.5-Mb chromosome 3p21.3 (AP20) region implicated in major human epithelial malignancies. Cancer Res. 2003, 63: 404-412.
Senchenko V, Liu J, Braga E, Mazurenko N, Loginov W, Seryogin Y, Bazov I, Protopopov A, Kisseljov FL, Kashuba V, Lerman MI, Klein G, Zabarovsky ER: Deletion mapping of cervical carcinomas using quantitative real time PCR identifies two frequently affected regions in 3p21.3. Oncogene. 2003, 22: 2984-2992. 10.1038/sj.onc.1206429.
Senchenko VN, Liu J, Loginov W, Bazov I, Angeloni D, Seryogin Y, Ermilova V, Kazubskaya T, Garkavtseva R, Zabarovska VI, Kashuba VI, Kisselev LL, Minna JD, Lerman MI, Klein G, Braga EA, Zabarovsky ER: Discovery of frequent homozygous deletions in chromosome 3p21.3 LUCA and AP20 regions in renal, lung and breast carcinomas. Oncogene. 2004, 23: 5719-5728. 10.1038/sj.onc.1207760.
Dammann R, Li C, Yoon JH, Chin PL, Bates S, Pfeifer GP: Epigenetic inactivation of a RAS association domain family protein from the lung tumour suppressor locus 3p21.3. Nat Genet. 2000, 25: 315-319. 10.1038/77083.
Kashuba VI, Li J, Wang F, Senchenko VN, Protopopov A, Malyukova A, Kutsenko AS, Kadyrova E, Zabarovska VI, Muravenko OV, Zelenin AV, Kisselev LL, Kuzmin I, Minna JD, Winberg G, Ernberg I, Braga E, Lerman MI, Klein G, Zabarovsky ER: RBSP3 (HYA22) is a tumor suppressor gene implicated in major epithelial malignancies. Proc Natl Acad Sci USA. 2004, 101: 4906-4911. 10.1073/pnas.0401238101.
Li J, Wang F, Protopopov A, Malyukova A, Kashuba V, Minna JD, Lerman MI, Klein G, Zabarovsky E: Inactivation of RASSF1C during in vivo tumor growth identifies it as a tumor suppressor gene. Oncogene. 2004, 23: 5941-5949. 10.1038/sj.onc.1207789.
Yeo M, Lin PS, Dahmus ME, Gill GN: A novel RNA polymerase II C-terminal domain phosphatase that preferentially dephosphorylates serine 5. J Biol Chem. 2003, 278: 26078-26085. 10.1074/jbc.M301791200.
Kashuba VI, Pavlova TV, Grigorieva EV, Kutsenko A, Yenamandra SP, Li J, Wang F, Protopopov AI, Zabarovska VI, Senchenko V, Haraldson K, Eshchenko T, Kobliakova J, Vorontsova O, Kuzmin I, Braga E, Blinov VM, Kisselev LL, Zeng YX, Ernberg I, Lerman MI, Klein G, Zabarovsky ER: High mutability of the tumor suppressor genes RASSF1 and RBSP3 (CTDSPL) in cancer. PLoS ONE. 2009, 4: e5231-10.1371/journal.pone.0005231.
Li J, Wang F, Haraldson K, Protopopov A, Duh FM, Geil L, Kuzmin I, Minna JD, Stanbridge E, Braga E, Kashuba VI, Klein G, Lerman MI, Zabarovsky ER: Functional characterization of the candidate tumor suppressor gene NPRL2/G21 located in 3p21.3C. Cancer Res. 2004, 64: 6438-6443. 10.1158/0008-5472.CAN-03-3869.
Ueda K, Kawashima H, Ohtani S, Deng WG, Ravoori M, Bankson J, Gao B, Girard L, Minna JD, Roth JA, Kundra V, Ji L: The 3p21.3 tumor suppressor NPRL2 plays an important role in cisplatin-induced resistance in human non-small-cell lung cancer cells. Cancer Res. 2006, 66: 9682-9690. 10.1158/0008-5472.CAN-06-1483.
Ji L, Nishizaki M, Gao B, Burbee D, Kondo M, Kamibayashi C, Xu K, Yen N, Atkinson EN, Fang B, Lerman MI, Roth JA, Minna JD: Expression of several genes in the human chromosome 3p21.3 homozygous deletion region by an adenovirus vector results in tumor suppressor activities in vitro and in vivo. Cancer Res. 2002, 62: 2715-2720.
Otani S, Takeda S, Yamada S, Sakakima Y, Sugimoto H, Nomoto S, Kasuya H, Kanazumi N, Nagasaka T, Nakao A: The tumor suppressor NPRL2 in hepatocellular carcinoma plays an important role in progression and can be served as an independent prognostic factor. J Surg Oncol. 2009, 100: 358-363. 10.1002/jso.21241.
Wei MN, Latif F, Bader S, Kashuba V, Chen JY, Duh FM, Sedico Y, Lee CC, Geil L, Kuzmin I, Zabarovsky E, Klein G, Zbar B, Minna JD, Lerman MI: Construction of a 600-kilobase cosmid clone contig and generation of a transcriptional map surrounding the lung cancer tumor suppressor gene (TSG) locus on human chromosome 3p21.3: progress toward the isolation of a lung cancer TSG. Cancer Res. 1996, 56: 1487-1492.
Dreijerink K, Braga E, Kuzmin I, Geil L, Duh FM, Angeloni D, Zbar B, Lerman MI, Stanbridge EJ, Minna JD, Protopopov A, Li J, Kashuba V, Klein G, Zabarovsky ER: The candidate tumor suppressor gene RASSF1A from human chromosome 3p21.3 is involved in kidney tumorigenesis. Proc Natl Acad Sci USA. 2001, 98: 7504-7509. 10.1073/pnas.131216298.
Kuzmin I, Gillespie JW, Protopopov A, Geil L, Dreijerink K, Yang Y, Vocke CD, Duh FM, Zabarovsky E, Minna JD, Rhim JS, Emmert-Buck MR, Linehan WM, Lerman MI: The RASSF1A tumor suppressor gene is inactivated in prostate tumors and suppresses growth of prostate carcinoma cells. Cancer Res. 2002, 62: 3498-3502.
Inda MM, Castresana JS: RASSF1A promoter is highly methylated in primitive neuroectodermal tumors of the central nervous system. Neuropathology. 2007, 27: 341-346. 10.1111/j.1440-1789.2007.00788.x.
Hesson LB, Cooper WN, Latif F: Evaluation of the 3p21.3 tumour-suppressor gene cluster. Oncogene. 2007, 26: 7283-301. 10.1038/sj.onc.1210547.
Burbee DG, Forgacs E, Zöchbauer-Müller S, Shivakumar L, Fong K, Gao B, Randle D, Kondo M, Virmani A, Bader S, Sekido Y, Latif F, Milchgrub S, Toyooka S, Gazdar AF, Lerman MI, Zabarovsky E, White M, Minna JD: Epigenetic inactivation of RASSFA1 in lung and breast cancers and malignant. J Natl Cancer Inst. 2001, 93: 691-699. 10.1093/jnci/93.9.691.
Honorio S, Agathanggelou A, Schuermann M, Pankow W, Viacava P, Maher ER, Latif F: Detection of RASSF1A aberrant promoter hypermethylation in sputum from chronic smokers and ductal carcinoma in situ from breast cancer patients. Oncogene. 2003, 22: 147-150. 10.1038/sj.onc.1206057.
Dammann R, Schagdarsurengin U, Seidel C, Strunnikova M, Rastetter M, Baier K, Pfeifer GP: The tumor suppressor RASSF1A in human carcinogenesis: an update. Histol Histopathol. 2005, 20: 645-663.
Agathanggelou A, Cooper WN, Latif F: Role of the Ras-association domain family 1 tumor suppressor gene in human cancers. Cancer Res. 2005, 65: 3497-3508. 10.1158/0008-5472.CAN-04-4088.
Pan ZG, Kashuba VI, Liu XQ, Shao JY, Zhang RH, Jiang JH, Guo C, Zabarovsky E, Ernberg I, Zeng YX: High frequency somatic mutations in RASSF1A in nasopharyngeal carcinoma. Cancer Biol Ther. 2005, 4: 1116-1122.
Vos MD, Ellis CA, Bell A, Birrer MJ, Clark GJ: Ras uses the novel tumor suppressor RASSF1 as an effector to mediate apoptosis. J Biol Chem. 2000, 275: 35669-35672. 10.1074/jbc.C000463200.
Ortiz-Vega S, Khokhlatchev A, Nedwidek M, Zhang XF, Dammann R, Pfeifer GP, Avruch J: The putative tumor suppressor RASSF1A homodimerizes and heterodimerizes with the Ras-GTP binding protein Nore1. Oncogene. 2002, 21: 1381-1390. 10.1038/sj.onc.1205192.
Shivakumar L, Minna JD, Sakamaki T, Pestell RG, White MA: The RASSF1A tumor suppressor blocks cell cycle progression and inhibits cyclin D1 accumulation. Mol Cell Biol. 2002, 22: 4309-4318. 10.1128/MCB.22.12.4309-4318.2002.
Sinha S, Singh RK, Alam N, Roy A, Roychoudhury S, Panda CK: Frequent alterations of hMLH1 and RBSP3/HYA22 at chromosomal 3p22.3 region in early and late-onset breast carcinoma: clinical and prognostic significance. Cancer Sci. 2008, 99: 1984-1991.
Wang F, Grigorieva EV, Li J, Senchenko VN, Pavlova TV, Anedchenko EA, Kudryavtseva AV, Tsimanis A, Angeloni D, Lerman MI, Kashuba VI, Klein G, Zabarovsky ER: HYAL1 and HYAL2 inhibit tumour growth in vivo but not in vitro. PLoS ONE. 2008, 3: e3031-10.1371/journal.pone.0003031.
Speissel B, Beahrs OH, Hermanek P, Hutter RVP, Scheibe O: TNM atlas: illustrated guide to the TNM/pTNM classification of malignant tumours. 1989, Springer-Verlag Berlin New York, 3
Travis WD, Coby TV, Corrin B, Shimosato Y, Brambilla E: World Health Organization International Histological Classification of Tumours; Histological typing of lung and pleural tumours. 1999, Springer Berlin, 3
Livak KJ, Schmittgen TD: Analysis of relative gene expression data using real-time quantitative PCR and the 2(-ΔΔCT) method. Methods. 2001, 25: 402-408. 10.1006/meth.2001.1262.
Relative Quantitation of Gene Expression. [http://docs.appliedbiosystems.com/pebiodocs/04303859.pdf]
Li J, Protopopov A, Wang F, Senchenko V, Petushkov V, Vorontsova O, Petrenko L, Zabarovska V, Muravenko O, Braga E, Kisselev L, Lerman MI, Kashuba V, Klein G, Ernberg I, Wahlestedt C, Zabarovsky ER: NotI subtraction and NotI-specific microarrays to detect copy number and methylation changes in whole genomes. Proc Natl Acad Sci USA. 2002, 99: 10724-10729. 10.1073/pnas.132271699.
Dunwell TL, Hesson LB, Pavlova T, Zabarovska V, Kashuba V, Catchpoole D, Chiaramonte R, Brini AT, Griffiths M, Maher ER, Zabarovsky E, Latif F: Epigenetic analysis of childhood acute lymphoblastic leukemia. Epigenetics. 2009, 4: 185-193.
Glantz SA: Primer of biostatistics. 1994, McGraw-Hill New York, 4
Li J, Protopopov AI, Gizatullin RZ, Kiss C, Kashuba VI, Winberg G, Klein G, Zabarovsky ER: Identification of new tumor suppressor genes based on in vivo functional inactivation of a candidate gene. FEBS Lett. 1999, 451: 289-294. 10.1016/S0014-5793(99)00598-0.
Protopopov AI, Li J, Winberg G, Gizatullin RZ, Kashuba VI, Klein G, Zabarovsky ER: Human cell lines engineered for tetracycline-regulated expression of tumor suppressor candidate genes from a frequently affected chromosomal region, 3p21. J Gene Med. 2002, 4: 397-406. 10.1002/jgm.283.
Shu J, Jelinek J, Chang H, Shen L, Qin T, Chung W, Oki Y, Issa JP: Silencing of bidirectional promoters by DNA methylation in tumorogenesis. Cancer Res. 2006, 66: 5077-5084. 10.1158/0008-5472.CAN-05-2629.
Mitra S, Mazumder Indra D, Bhattacharya N, Singh RK, Basu PS, Mondal RK, Roy A, Zabarovsky ER, Roychoudhury S, Panda CK: RBSP3 is frequently altered in premalignant cervical lesions: Clinical and prognostic significance. Genes Chromosomes Cancer. 2010, 49: 155-170.
Anedchenko EA, Kiseleva NP, Dmitriev AA, Kiselev FL, Zabarovsky ER, Senchenko VN: Tumor suppressor gene RBSP3 in cervical carcinoma: copy number and transcriptional level. Mol Biol (Moscow), Russian. 2007, 41: 86-95.
Frigola J, Song J, Stirzaker C, Hinshelwood RA, Peinado MA, Clark SJ: Epigenetic remodeling in colorectal cancer results in coordinate gene suppression across an entire chromosome band. Nat Genet. 2006, 38: 540-549. 10.1038/ng1781.
Gao C, Furge K, Koeman J, Dykema K, Su Y, Cutler ML, Werts A, Haak P, Woude Vande GF: Chromosome instability, chromosome transcriptome, and clonal evolution of tumor cell populations. Proc Natl Acad Sci USA. 2007, 104: 8995-9000. 10.1073/pnas.0700631104.
Smith IM, Glazer CA, Mithani SK, Ochs MF, Sun W, Bhan S, Vostrov A, Abdullaev Z, Lobanenkov V, Gray A, Liu C, Chang SS, Ostrow KL, Westra WH, Begum S, Dhara M, Califano J: Coordinated activation of candidate proto-oncogenes and cancer testes antigens via promoter demethylation in head and neck cancer and lung cancer. PLoS ONE. 2009, 4: e4961-10.1371/journal.pone.0004961.
Oyama T, Osaki T, Baba T, Nagata Y, Mizukami M, So T, Nakata S, Ichiki Y, Uramoto H, Sugaya M, Yoshimatsu T, Morita M, Hanagiri T, Sugio K, Kawamoto T, Yasumoto K: Molecular genetic tumor markers in non-small cell lung cancer. Anticancer Res. 2005, 25: 1193-1196.
Petty RD, Nicolson MC, Kerr KM, Collie-Duguid E, Murray GI: Gene expression profiling in non-small cell lung cancer: from molecular mechanisms to clinical application. Clin Cancer Res. 2004, 10: 3237-3248. 10.1158/1078-0432.CCR-03-0503.
The pre-publication history for this paper can be accessed here:http://www.biomedcentral.com/1471-2407/10/75/prepub
The authors are grateful to Ms. Leah V. Lerman, esq. for the help in preparation of the manuscript. This study was financially supported by the Russian Foundation for Basic Research (project no 05-04-49408), The Agency of Science and Innovations (State Contracts no 02.512.11.2241, Agreement 1-2008/07/07) and by research grants from the Swedish Cancer Society, the Swedish Research Council, the Swedish Foundation for International Cooperation in Research and Higher Education (STINT), the Swedish Institute, and Karolinska Institute.
The authors declare that they have no competing interests.
VNS designed the experiments, carried out the qPCR analysis, and wrote the manuscript; EAA collected and classified all lung tumors, performed the qPCR data; TTK - contributed with clinical information for these patients and drafting the manuscript; GSK - was responsible for primers and probes design and analyzed the data; AAD -was responsible for calculations and statistical analysis. VIZ, TVP, VIK performed the experiments and analyzed the NotI-data; MIL helped in the evaluation of the results and revised the manuscript critically for important intellectual content; ERZ was responsible for design study and coordination, wrote the manuscript. All authors read and approved the final manuscript.
About this article
Cite this article
Senchenko, V.N., Anedchenko, E.A., Kondratieva, T.T. et al. Simultaneous down-regulation of tumor suppressor genes RBSP3/CTDSPL, NPRL2/G21 and RASSF1A in primary non-small cell lung cancer. BMC Cancer 10, 75 (2010). https://doi.org/10.1186/1471-2407-10-75
- Small Cell Lung Cancer
- Squamous Cell Carcinoma
- NSCLC Sample
- Primary NSCLC
- NSCLC Case