
Clinical Cancer Research 14, 3746, June 15, 2008. doi: 10.1158/1078-0432.CCR-07-4552
© 2008 American Association for Cancer Research
Association of p16 Homozygous Deletions with Clinicopathologic Characteristics and EGFR/KRAS/p53 Mutations in Lung Adenocarcinoma
Reika Iwakawa1,6,
Takashi Kohno1,
Yoichi Anami4,
Masayuki Noguchi4,
Kenji Suzuki2,
Yoshihiro Matsuno3,
Kazuhiko Mishima5,
Ryo Nishikawa5,
Fumio Tashiro6 and
Jun Yokota1
Authors' Affiliations: 1 Biology Division, National Cancer Center Research Institute, 2 Thoracic Surgery Division and 3 Diagnostic Pathology Division, National Cancer Center Hospital, Tokyo, Japan; 4 Department of Pathology, Institute of Basic Medical Sciences, Graduate School of Comprehensive Human Sciences, Tsukuba University, Ibaraki, Japan; 5 Department of Neuro-Oncology, Comprehensive Cancer Center, International Medical Center, Saitama Medical University, Saitama, Japan; and 6 Department of Biological Science and Technology, Faculty of Industrial Science and Technology, Tokyo University of Science, Chiba, Japan
Requests for reprints: Jun Yokota, Biology Division, National Cancer Center Research Institute, 1-1, Tsukiji 5-chome, Chuo-ku, Tokyo 104-0045, Japan. Phone: 81-33547-5272; Fax: 81-33542-0807; E-mail: jyokota{at}ncc.go.jp.
 |
Abstract
|
|---|
Purpose: The p16 gene is frequently inactivated in lung adenocarcinoma. In particular, homozygous deletions (HD) have been frequently detected in cell lines; however, their frequency and specificity is not well-established in primary tumors. The purpose of this study was to elucidate the prevalence and the timing for the occurrence of p16 HDs in lung adenocarcinoma progression in vivo.
Experimental Design: Multiple ligation-dependent probe amplification was used for the detection of p16 HDs in 28 primary small-sized lung adenocarcinomas and 22 metastatic lung adenocarcinomas to the brain. Cancer cells were isolated from primary adenocarcinoma specimens by laser capture microdissection. HDs were confirmed by quantitative real-time genomic PCR analysis.
Results: HDs were detected in 8 of 28 (29%) primary tumors, including 2 of 8 (25%) noninvasive bronchioloalveolar carcinomas, and 5 of 22 (26%) brain metastases, respectively. No significant associations were observed between p16 HDs and gender, age, smoking history, stage, and prognosis. HDs were detected with similar frequencies (17–29%) among adenocarcinomas with epidermal growth factor receptor (EGFR) mutations, with KRAS mutations, and without EGFR/KRAS mutations, and with similar frequencies (22–28%) between adenocarcinomas with and without p53 mutations.
Conclusions: p16 HDs occur early in the development of lung adenocarcinomas and with similar frequencies among EGFR type, KRAS type, and non-EGFR/KRAS type lung adenocarcinomas. Tobacco carcinogens would not be a major factor inducing p16 HDs in lung adenocarcinoma progression.
Adenocarcinoma is the most common histologic type of lung cancer. Recent molecular analyses have indicated that lung adenocarcinoma can be divided into at least three types; the epidermal growth factor receptor (EGFR) type, the KRAS type, and the non-EGFR/KRAS type, based on accumulated genetic alterations in adenocarcinoma cells (1, 2). The p16 tumor suppressor gene is frequently inactivated in lung adenocarcinomas, most prominently through promoter methylation and homozygous deletion (HD), and less frequently through intragenic mutation (3–6). In particular, p16 methylation is known to occur in close association with tobacco carcinogen exposure and preferentially in the KRAS type adenocarcinomas (1). Molecular analyses of small-sized adenocarcinomas revealed that p16 methylation occurs in the course of progression from noninvasive bronchioloalveolar carcinomas (BAC) to invasive ones, and its occurrence is associated with smoking history, staging, and prognosis (7). However, it is still unknown whether p16 HD also occurs preferentially in the KRAS type or not. In addition, the involvement of tobacco smoking in the occurrence of p16 HD is also unclear because both positive and negative associations between tobacco smoking and HD have been reported (5, 6). Indeed, p16 HDs have been reported to occur with a highly variable frequency in primary lung adenocarcinomas (0–40%; refs. 3–6, 8). HD can be easily masked if a large fraction of noncancerous cells are contaminated in tumor tissues, in particular, in the analysis of small tumors with a BAC component. Thus, such an inconsistency could come not only from the diversity and/or heterogeneity of lung adenocarcinomas but also from tumor tissues used for the analysis and also from the methods used for the detection of HDs. For instance, immunohistochemical analyses of p16 proteins have shown that a considerable fraction of adenocarcinomas are negative for p16 protein expression without clear evidence of an inactivating event by molecular analyses of the p16 gene (4, 5, 7). Thus, it has been assumed that p16 HD is a causative event for the absence of p16 protein in these adenocarcinoma cases. However, due to the lack of comprehensive analysis for p16 HDs, the prevalence, specificity, and the timing for the occurrence of p16 HDs in lung adenocarcinoma progression in vivo is still unclear.
In this study, we investigated the association of p16 HDs with clinicopathologic characteristics, including smoking history, of lung adenocarcinomas and with the status of EGFR, KRAS, and p53 mutations in lung adenocarcinomas. To obtain more critical information than previous studies for the timing of its occurrence in the progression of lung adenocarcinomas, we analyzed two typical stages of lung adenocarcinomas. One is a primary small-sized adenocarcinoma of
2 cm in maximum diameter as a representative of early stage lung adenocarcinomas, and the other is a brain metastasis as a representative of late stage lung adenocarcinomas. To exclude the possible overlooking of HDs in the analysis due to contamination of noncancerous cells in tumor tissues, a laser capture microdissection method was applied for the isolation of cancer cells from small-sized primary lung adenocarcinomas. Brain metastases are known to be relatively solid and generally contain a small amount of noncancerous cells in tumor tissues; therefore, macrodissected samples were used for the analysis. Recently, a simple and effective method, which is called multiplex ligation-dependent probe amplification (MLPA), to measure the copy number of up to 45 genomic loci in a single experiment, was developed and has been applied for the detection of large deletions in/of various genes in the human genome DNA sequence (9). For instance, this method was successfully applied for the detection of p16 hemizygous deletions in melanoma families (10) and p16 HDs in head and neck squamous cell carcinoma cell lines (11). Thus, we applied this method for the detection of p16 HDs in surgically resected lung adenocarcinomas. HDs were detected with similar frequencies of 23% to 40% in lung adenocarcinomas of any progression stage—from noninvasive carcinomas to advanced ones. The deletions were also detected with similar frequencies among EGFR type, KRAS type, and non-EGFR/KRAS type adenocarcinomas, and were not associated with tobacco smoking and poor prognosis.
 |
Materials and Methods
|
|---|
Patients and tissues. Twenty-eight primary small-sized (
2 cm in maximum diameter) adenocarcinomas, 22 brain metastases, and corresponding noncancerous tissues were obtained at surgery from lung adenocarcinoma patients treated at the National Cancer Center Hospital, Tokyo and at the Saitama Medical University. In 4 of the 22 brain metastasis cases, the corresponding primary tumors were also obtained at surgery. The surgically resected specimens were fixed routinely with 10% formalin and embedded in paraffin for histologic examination. All the sections were stained with H&E and examined by light microscopy. The tumors were pathologically diagnosed according to the tumor-node-metastasis classification of malignant tumors (12). The primary small-sized adenocarcinomas were further classified into three types according to the histologic classification of small-sized adenocarcinoma of the lung reported by Noguchi and colleagues (13), and there were 8 type B, 15 type C, and 5 type D adenocarcinomas. These 28 tumors and corresponding noncancerous tissues were fixed with methanol and embedded in paraffin. Cancer cells were then microdissected using the Pixcell Laser Capture Microdissection system (Arcturus Engineering), and their genomic DNAs were extracted as described previously (14). Representative figures of type B and C tumors before and after microdissection have been shown in our previous article (14). All of the 22 brain metastases, corresponding primary tumors, and noncancerous tissues were macrodissected and stored at –80°C until DNA extraction. Genomic DNAs were prepared as previously described (15).
Cell lines. Three lung adenocarcinoma cell lines, H2126, A549, and PC3, were used as positive controls of p16 HDs in MLPA and quantitative real-time genomic PCR (QRT-G-PCR) analysis. Genomic DNA were prepared as previously described (15). Incidence as well as regions of p16 HDs in 55 lung adenocarcinoma cell lines, including these three cell lines, were determined by multiplex PCR analysis in our previous studies (3).
MLPA analysis. MLPA was carried out using the P024B kit for the 9p21 CDKN2A/2B region (MRC-Holland) according to the manufacturer's protocol. The kit contains 12 probes covering the p15/p14ARF/p16 genes, 12 probes for 9 other genes on chromosome 9p, and 15 control probes for nonchromosome 9p loci. Genes on chromosome 9p included in this screen were TEK, ELAVL2, CDKN2B (p15), CDKN2A (p14ARF/p16), MTAP, IFNA1, KIAA1354, IFNW1, IFNB1, MLLT3, and DOCK8 (FLJ00026) in the order from centromere to telomere (Fig. 1A
). Experiments were done in a half volume until the ligation reaction step, and then by the supplied protocol. Briefly, 12.5 to 50 ng of genomic DNA in 2.5 to 5 µL of TE buffer were heat-denatured and hybridized to probes for 16 h at 60°C. The hybridized probes were then ligated and amplified by PCR of 35 cycles at 60°C for 30 s and 72°C for 60 s. PCR products were separated by capillary electrophoresis using the ABI3700 Automated Capillary DNA Sequencer with a 50-cm capillary array, ABI POP-6 polymer, and GeneScan-ROX 500 size standards (Applied Biosystems). Analysis was automated using ABI PRISM GeneScan Analysis software version 3.7, and Genotyper Analysis software version 3.7 (Applied Biosystems). PCR and electrophoresis, respectively, were done in duplicate, and the mean of four values was calculated and considered to be the DNA copy number ratio of each locus.

View larger version (35K):
[in this window]
[in a new window]
[Download PPT slide]
|
Fig. 1. HDs of the p16 gene detected by MLPA analysis. A, physical map of the chromosome 9p region. Top, positions of 24 probes from nos. 1 to 24 (red) in the order from centromere (Cen.) to telomere (Tel.). Bottom, positions of 12 genes in this region (light blue). Exons (numbers) and orientations (arrows) of the p15, p14ARF, p16, MTAP, and MLLT3 genes. B, definition of p16 HDs by MLPA analysis. Relative copy numbers of the 24 loci in two lung adenocarcinoma cell lines, H2126 and PC3 (blue dots/lines), and those in mixed samples of the same amounts of DNA from these cell lines and normal lung tissue (pink dots/lines; H2126 half and PC3 half, respectively). C, results of MLPA analysis for cases B-2P, C-10P, D-4P, and N2131M/N2133P showing HDs of the p16 gene. As for case C-10P, both the BAC component (BAC; blue dots/lines) and the invasive region (in; pink dots/lines) of the tumor were analyzed. N2131M (blue) and N2133P (pink) are a brain metastasis and the primary tumor from the same patient, respectively.
|
|
Data analysis. The relative DNA copy number of each locus was calculated as follows. First, the value for the sum of 15 autosomal control probe peak areas in a lung adenocarcinoma sample was adjusted to the mean value for those in two or three normal lung tissue samples in the same run. Each of the 39 probe peak areas was then divided by the sum of the peak areas of the 15 autosomal control probe peak areas. Finally, the relative DNA copy number ratio of each of the 39 chromosome loci in the lung adenocarcinoma samples against the normal lung tissue samples was calculated. Theoretically, ratios close to 1.0 indicated that two DNA copies were present (i.e., wild-type), a ratio of 0.5 indicated that one copy was absent (i.e., hemizygous deletion), and a ratio of 0.0 indicated that both copies were absent (i.e., homozygous deletion). The criteria for hemizygous deletions and HDs of the 39 loci by MLPA analysis were defined using MLPA data from three lung adenocarcinoma cell lines with p16 HDs as described in Results and Discussion. The p16 gene was then defined as homozygously deleted or not in each sample. Because the purpose of this study is to evaluate the prevalence of p16 HD, hemizygous deletions were not evaluated from MLPA data.
QRT-G-PCR analysis. TaqMan-MGB probes and primers were designed using Primer Express software (Applied Biosystems) and were optimized according to the manufacturer's guidelines. Target and reference locus probes were labeled with FAM and VIC, respectively. Probe sequences were as follows: no. 9, 5'-AACTCCTCCACTGATTAC-3'; and 2p14, 5'-CCAGCCTATTCCTGC-3'. Primer sequences were as follows: no. 9-F/R, 5'-GGGTCTCCTTCATTTGGTGAAA-3'/5'-GGATCCCAGGGAGGAGAGTCT-3'; and 2p14-F/R, 5'-AAGAAGACTGGAGTGGTGTTTGG-3'/5'-CACAATGCTGAATACTGTCAATGAAA-3'. PCR was carried out in duplicate using 1 ng of DNA as a template. Primer and probe concentrations were optimized for each target according to the manufacturer's instructions. The PCR program consisted of 50°C for 2 min and 95°C for 15 min followed by 45 cycles of 95°C for 15 s and 60°C for 1 min. Standard curves for the copy numbers of the target and reference genes were generated using serially diluted (0.04–25 ng) normal lung tissue DNA. Data analysis was carried out using ABI Prism 7900HT Sequence Detection Software. DNA copy number ratios were calculated as the average copy number of the target locus divided by the average copy number of the reference locus, and then normalized against the normal lung tissue DNA to give a normalized DNA copy number ratio.
Immunohistochemistry. Immunohistochemical analysis of p16 protein was performed as described previously (7) using 4-µm sections cut from methanol-fixed and paraffin-embedded specimens of 28 primary small-sized adenocarcinomas, which were subjected to MLPA analysis. Positivity for p16 staining was scored by the same criteria as previously described (7). In brief, only nuclear staining was scored, and was considered to be positive when it was more intense than the background cytoplasmic staining. If <10% of tumor cells displayed p16 protein staining, it was judged negative; and if >10% of them showed strong staining, it was judged positive.
Mutation analysis of the EGFR, KRAS, and p53 genes. Thirteen of the 28 primary tumors, and 16 of the 22 brain metastases were previously examined for mutations of exons 18 to 21 in the EGFR gene, of exons 1 and 2 in the KRAS gene, and of exons 4 to 8 in the p53 gene by genomic PCR and direct sequencing (14, 16). The remaining 15 primary tumors and 6 brain metastases were also examined for these mutations in this study using the same methods.
Statistical analysis. Fisher's exact test was used to assess the association of p16 HDs with clinicopathologic characteristics or mutations of the EGFR, KRAS, and p53 genes. P < 0.05 was considered to be statistically significant.
 |
Results and Discussion
|
|---|
Detection of HDs by MLPA analysis. We previously determined regions of p16 HDs in various lung cancer cell lines by multiplex genomic PCR analysis (3, 17). Thus, we first validated the sensitivity and specificity of the MLPA method to detect p16 HDs using several lung adenocarcinoma cell lines. No PCR amplification was observed for probes that hybridized to the sequences of HD regions in all the cell lines examined. Representative results of MLPA analysis for two cell lines, H2126 and PC3, are shown in Fig. 1B. We then applied MLPA analysis for the detection of p16 HDs in surgically resected adenocarcinoma samples. Cancer cells of small-sized primary adenocarcinomas were isolated by laser capture microdissection, and brain metastases generally contain a small fraction of noncancerous cells. However, a complete absence of PCR products was not observed for these samples at any chromosomal loci by MLPA analysis. Therefore, we next defined the criteria for p16 HDs by MLPA analysis for surgically resected samples. For this purpose, DNA samples with virtual hemizygous deletions were prepared by mixing the same amounts of DNA from normal lung tissue and from three lung cancer cell lines with p16 HDs, H2126, A549, and PC3. In total, 64 probe loci were homozygously deleted in these cell lines. The mean ± 3 SD of relative DNA copy number ratios for the 64 probe loci in the HD regions among these mixed samples was 0.54 ± 0.17; thus, the range of virtual hemizygous deletions was 0.37 to 0.71. Indeed, none of the 64 probe loci in the HD regions of the mixed samples showed a DNA copy number ratio of >0.71 or <0.37 (Fig. 1B). Therefore, if the DNA copy number ratio for a locus was <0.37, the locus was judged as homozygously deleted in surgically resected lung adenocarcinoma samples.
Under this criterion, one or more loci in the p15/p14ARF/p16 gene region were judged as homozygously deleted in 8 of 28 primary tumors and in 5 of 22 brain metastases (Table 1
). Representative adenocarcinoma cases judged as having HDs of this region by MLPA analysis are shown in Fig. 1C. Case B-2P showed a HD of only one locus (probe no. 13) in exon 2 of the p16 gene. On the other hand, most of the other 12 cases showed HDs of several genes, including the p14ARF and/or p16 genes (Table 2
). Both the invasive region (C-10P-in) and the BAC component (C-10P-BAC) of case C-10P showed HDs of the same loci from p15 exon 1 to IFNW1. In case D-4P, the p15, p14ARF, p16, and MTAP genes were homozygously deleted. A metastasis to the brain, case N2131M, showed a large HD of the region from the ELAVL2 gene to the IFNB1 gene including the p16 gene. Furthermore, the relative DNA copy number ratios of the corresponding loci were also decreased in the corresponding primary tumor, N2133P, although the ratios were underrepresented as HDs. This may be attributed to a contamination of noncancerous cells in the primary tumor sample because this sample was macrodissected but not microdissected. There was no probe loci deleted in all the cases, but six loci (probe nos. 4, 5, 6, 8, 9, and 13) were deleted in 12 of 13 cases, including exon 2 of the p16 gene (Table 2).
Confirmation of HDs by QRT-G-PCR analysis. QRT-G-PCR was performed to confirm p16 HDs detected by MLPA analysis. Among 15 control probes for nonchromosome 9p loci, a probe for a chromosome 2p14 locus showed the most consistent DNA copy number ratios of nearly 1.0 in surgically resected adenocarcinoma samples. Thus, two primer sets were designed for QRT-G-PCR analysis. One was for the amplification of an MLPA probe locus between the p14ARF gene and the p16 gene (probe no. 9 in Fig. 1A, and hereinafter referred to as the p14ARF/p16 locus), and the other was for the amplification of a control probe locus on chromosome 2p14.
QRT-G-PCR was carried out for 16 brain metastases and the corresponding 4 primary tumors. Among them, five brain metastases were judged as having HDs of the p14ARF/p16 locus by MLPA analysis. DNA from normal lung tissue was used as a negative control, and DNA from two lung cancer cell lines with p16 HDs, A549 and H2126, were used as positive controls. No PCR products were detected for the p14ARF/p16 locus in the A549 and H2126 cell lines (data not shown). The DNA copy number ratios of the p14ARF/p16 locus in all five cases that had been determined to have HDs of this locus by MLPA analysis were <0.45. Thus, these five cases were also judged as having less than one copy of the p14ARF/p16 gene by QRT-G-PCR analysis. We further analyzed the association between the results of MLPA analysis and those of QRT-G-PCR analysis among all the 20 cases. The correlation coefficient was 0.95, and a highly significant correlation was observed between them (P = 1.87 x 10–10; Fig. 2
). This result gave the agreement for the appropriateness of the criterion for p16 HDs in MLPA analysis.

View larger version (10K):
[in this window]
[in a new window]
[Download PPT slide]
|
Fig. 2. Correlation between the results of MLPA analysis and those of QRT-G-PCR analysis. Relative DNA copy numbers of the probe no. 9 locus in Fig. 1A on chromosome 9 against the control locus on chromosome 2 defined by the QRT-G-PCR (X-axis) and MLPA analyses (Y-axis). Sixteen brain metastases and four of the corresponding primary tumors were analyzed. Values <0.5 by QRT-G-PCR analysis and values <0.37 by MLPA analysis (indicated by lines) were considered as having HDs.
|
|
Occurrence of p16 HDs in early stage lung adenocarcinoma. Based on the results of MLPA analysis, together with the confirmation by QRT-G-PCR analysis, we concluded that p16 HDs were present in 8 of 28 (29%) small-sized primary tumors and in 5 of 22 (23%) brain metastases (Table 1). Among four pairs of brain metastases and the corresponding primary tumors, only one metastasis (N2131M) was judged as having a p16 HD as described above and shown in Fig. 1C, and none of the remaining three cases showed HD in either brain metastases or primary tumors. Small-sized adenocarcinomas were further classified histologically into types B, C, and D (Fig. 3
). Type B tumors are noninvasive BACs and type C tumors are invasive adenocarcinomas with noninvasive BAC components. p16 HDs were detected in 2 of 8 (25%) type B tumors and in 4 of 15 (27%) type C tumors, in particular, in BAC components of two type C tumors (cases C-10P and C-65P). Type D tumors are invasive, poorly differentiated adenocarcinomas and p16 HDs were detected in 2 of 5 (40%) of the type D tumors. Thus, p16 HDs were present with similar frequencies of 25% to 40% in noninvasive and invasive primary adenocarcinomas. We previously reported that p16 HDs were present in 20 of 55 (36%) lung adenocarcinoma cell lines (17). Thus, the frequency of p16 HDs in small-sized primary adenocarcinomas was not significantly lower than, but similar to, those in brain metastases and cultured cell lines. These results strongly indicate that most p16 HDs detected in the cell lines occurred in vivo during adenocarcinoma progression and were retained during cultivation in vitro of adenocarcinoma cell lines. Similar frequencies of p16 HDs between noninvasive and invasive adenocarcinomas and between primary and metastatic adenocarcinomas further indicate that p16 HD occurs early in the multistage carcinogenic process of lung adenocarcinomas.

View larger version (72K):
[in this window]
[in a new window]
[Download PPT slide]
|
Fig. 3. Histology (H&E; original magnification, x100 for B-2P/C-12P/C-64P and x40 for C-10P; top) of small-sized lung adenocarcinomas and immunohistochemical staining of p16 protein (original magnification, x400; bottom) in the corresponding tissues. B-2P, C-10P, and D-12P are representative cases of type B, C, and D tumors with p16-negative staining, and C-64P is a representative case of type C tumors with p16-positive staining. A BAC component (BAC) and an invasive region (INV) are shown separately for C-10P.
|
|
No p16 expression in tumors with p16 HDs. Because one of the purposes of this study was to confirm the immunohistochemical negativity of adenocarcinoma cells with p16 HDs, we further performed an immunohistochemical analysis on 28 cases of small-sized adenocarcinomas, all of which were subjected to MLPA analysis. Fifteen of the 28 cases (54%) showed negative immunoreactivity for p16 protein. As predicted, all eight cases with p16 HDs following MLPA analysis were negative for p16 protein expression (Table 3
). Representative results of immunohistochemical staining for type B, C, and D tumors are shown in Fig. 3. It was noted that most tumor cells were clearly negative for nuclear p16 staining in these eight cases, and that both BAC components and invasive regions were negative in all four type C tumors with p16 HDs. This result strongly supports the reliability of MLPA analysis for the detection of p16 HDs in primary lung adenocarcinomas. Thus, it was concluded that p16 protein is not expressed in a considerable fraction of small-sized lung adenocarcinomas due to HDs of the p16 gene. The consistency of the results of immunohistochemistry with that of MLPA analysis further supports the hypothesis that p16 HD occurs early in lung adenocarcinoma progression.
View this table:
[in this window]
[in a new window]
|
Table 3. Associations of p16 HDs with p16 protein expression, clinicopathologic characteristics, and EGFR/KRAS/p53 mutations in lung adenocarcinoma
|
|
Previously, Dr. Noguchi, who is one of the authors of this article, and his colleagues reported that 29 of 57 (51%) small-sized lung adenocarcinomas of types A to F were negative for p16 immunostaining (7). In their report, the frequency of p16 negativity was higher in smokers and in patients with non-BACs. Aberrant methylation of the p16 gene promoter was detected more frequently in advanced BACs (type C) and non-BAC (types D–F) than in BACs (types A and B). In this study, the overall frequency (54%) of p16 negativity was quite similar to their reports, and the frequency was also higher in smokers than in nonsmokers (9 of 13 versus 6 of 15), and in type D (5 of 5) than in type B (7 of 15) or type C (3 of 8). Thus, it was highly suggested that a majority of cases with negative p16 expression without p16 HDs could be due to methylation of the p16 gene promoter and the methylation was associated with tobacco smoking. However, methylation and HD are not likely to coexist with each other because the regions of HDs in most of the 13 cases with HDs included the promoter region of the p16 gene (from probe no. 8 to probe no. 11). Thus, either methylation or HD of the p16 gene was suggested to occur equally and independently in smokers, whereas HD dominates in nonsmokers.
Association of p16 HDs with clinicopathologic characteristics. We then investigated the association of p16 HDs with clinicopathologic characteristics, including smoking history, of patients with lung adenocarcinomas (Table 3). No significant associations were observed between p16 HDs and gender, age, and smoking history in all 50 cases analyzed (Table 3), as well as in 28 cases of primary adenocarcinomas or in 22 cases of brain metastases (data not shown). p16 HDs were not associated with pathologic stage nor with 5-year survival in 28 cases of primary adenocarcinomas (Table 3). Previously, p16 methylation was shown to be associated with smoking history (1, 7). Thus, no association of p16 HDs with smoking history was in contrast with the status of p16 methylation in lung adenocarcinomas. The results indicate that causative factors for p16 HD were different from those for p16 methylation, although both alterations result in the inactivation of the same gene. Previously, Kraunz and colleagues reported that p16 HD occurred at a higher frequency in never-smokers as compared with former and current smokers (6). Although such an association was not observed in this study, both studies indicated the absence of a positive association between p16 HD and smoking history, and the possible association of p16 HD with other causative factors for lung adenocarcinoma. Thus, further studies are needed for the elucidation of such factors because little is known about the environmental as well as genetic risk factors for lung adenocarcinoma. The absence of any association between p16 HDs and 5-year survival was also in contrast with the status of p16 methylation in lung adenocarcinomas. However, because the number of cases examined was small (28 cases), and the patients with poor prognosis had a higher frequency of p16 HDs than those with good prognosis (44% versus 21%), further studies will be required on this subject. The absence of any association between p16 HDs and pathologic stage further supports their occurrence in the early stages, rather than in the late stages, of lung adenocarcinoma progression.
Association of p16 HDs with EGFR, KRAS, and p53 mutations. We next evaluated the association of p16 HDs with the status of EGFR, KRAS, and p53 mutations in these samples (Tables 1 and 3). Then, based on the results of these molecular analyses and a previous p16 methylation analysis (7), a stepwise malignant progression model for small-sized lung adenocarcinoma was depicted as shown in Fig. 4
. As previously reported (14, 18), EGFR and KRAS mutations were detected in a mutually exclusive manner in lung adenocarcinomas. EGFR mutations were most frequently detected in type B tumors (8 of 8, 100%), suggesting the involvement of EGFR mutations in the formation of noninvasive BACs. Although KRAS mutations were detected only in type C and D tumors in this study, it was recently reported that the mutations were frequently detected in noninvasive adenocarcinomas and also in atypical adenomatous hyperplasias (AAH; ref. 19). Thus, it is likely that either EGFR or KRAS mutations occur prior to HDs and methylations of the p16 gene in the progression of BACs, although it is also possible that p16 alterations occur earlier than depicted in Fig. 4. Frequencies of p53 mutations in primary adenocarcinomas were the highest in type D (5 of 5, 100%), intermediate in type C (9 of 15, 60%), and the lowest in type B (2 of 8, 25%), suggesting the accumulation of the mutations during progression from noninvasive adenocarcinomas to invasive ones. Thus, it was likely that p53 mutations had accumulated in adenocarcinomas with p16 HDs and/or EGFR/KRAS mutations. It is indispensable to analyze a considerable number of AAHs as well as type A tumors to fully understand the timing of each genetic alteration in sequential progression of lung adenocarcinomas because AAH is a putative precursor of peripheral lung adenocarcinoma including BAC (20), and sequential progression from type A to type C tumors through type B tumors was strongly indicated in previous studies (7, 21). However, because AAH is not routinely resected by surgery and because type A tumors are usually very small, the number of tumors as well as the amount of DNA obtained was not enough for the present study.

View larger version (20K):
[in this window]
[in a new window]
[Download PPT slide]
|
Fig. 4. A stepwise malignant progression model of small-sized lung adenocarcinoma in association with accumulated genetic alterations in cells. BAC, bronchioloalveolar carcinoma; type C, localized bronchioloalveolar carcinoma with foci of active fibroblastic proliferation; type D, poorly differentiated adenocarcinoma.
|
|
There was no specificity for the occurrence of p16 HDs among adenocarcinomas with regard to other accompanied genetic alterations. In particular, p16 HDs were detected with similar frequencies (17–29%), irrespective of the presence or absence of EGFR, KRAS, and p53 mutations (Table 3). This result indicates that p16 HDs occur with similar frequencies in EGFR type as well as KRAS type and non-EGFR/KRAS type adenocarcinomas. In this model, type C tumors were considered to progress from BACs, as previously indicated by Aoyagi and colleagues (21). However, it was still unclear whether type D tumors arise de novo by a distinct pathway from tumors with BAC components or progress from these tumors. Low frequencies of EGFR mutations and high frequencies of p53 mutations in type D tumors indicate that type D tumors progress from either the KRAS types or non-EGFR/KRAS types. Haneda and colleagues also reported a low frequency of EGFR mutations in type D tumors (22), supporting the presence of a non-EGFR pathway for the development of type D poorly differentiated adenocarcinomas.
Conclusions. MLPA analysis of microdissected small-sized primary adenocarcinoma cells revealed that p16 HDs are present in 20% to 40% of adenocarcinomas irrespective of the presence of mutations in the EGFR, KRAS, and p53 genes, and occur early in the development of lung adenocarcinomas. HDs were not associated with smoking history of the patients. It has been indicated that smoking is a major factor inducing p16 methylation as well as KRAS mutation. In contrast, EGFR mutations frequently occur in female nonsmokers and are associated with bronchioloalveolar morphology of adenocarcinomas. Interestingly, p16 HDs did not coexist with specific genetic alterations in adenocarcinoma cells. Thus, causative factors for p16 HD would be different from those for p16 methylation, KRAS mutations, and EGFR mutations. To elucidate the causative role for the occurrence of p16 HDs in multistage lung carcinogenesis, further studies should focus on the identification of environmental and genetic factors for the induction of DNA double-strand breaks surrounding the p16 gene locus and their repair systems.
 |
Disclosure of Potential Conflicts of Interest
|
|---|
No potential conflicts of interest were disclosed.
 |
Acknowledgments
|
|---|
We thank Drs. Higashi and Yokoyama of FALCO Biosystems, Ltd., for their technical support and critical discussion in MLPA analysis. We also thank Drs. Nakanishi and Matsumoto for the laser capture microdissection of cancer cells from primary small-sized lung adenocarcinomas.
 |
Footnotes
|
|---|
Grant support: Grants-in-Aid from the Ministry of Health, Labor and Welfare for the Third-Term Comprehensive 10-Year Strategy for Cancer Control and for Cancer Research (16-1) and a Grant-in-Aid for the Program for Promotion of Fundamental Studies in Health Sciences of the National Institute of Biomedical Innovation (NiBio).
The costs of publication of this article were defrayed in part by the payment of page charges. This article must therefore be hereby marked advertisement in accordance with 18 U.S.C. Section 1734 solely to indicate this fact.
Received 10/ 5/07;
revised 2/24/08;
accepted 2/28/08.
 |
References
|
|---|
- Toyooka S, Tokumo M, Shigematsu H, et al. Mutational and epigenetic evidence for independent pathways for lung adenocarcinomas arising in smokers and never smokers. Cancer Res 2006;66:1371–5.[Abstract/Free Full Text]
- Sato M, Shames DS, Gazdar AF, Minna JD. A translational view of the molecular pathogenesis of lung cancer. J Thorac Oncol 2007;2:327–43.[Medline]
- Hamada K, Kohno T, Kawanishi M, Ohwada S, Yokota J. Association of CDKN2A(p16)/CDKN2B(p15) alterations and homozygous chromosome arm 9p deletions in human lung carcinoma. Genes Chromosomes Cancer 1998;22:232–40.[CrossRef][Medline]
- Sanchez-Cespedes M, Reed AL, Buta M, et al. Inactivation of the INK4A/ARF locus frequently coexists with TP53 mutations in non-small cell lung cancer. Oncogene 1999;18:5843–9.[CrossRef][Medline]
- Sanchez-Cespedes M, Decker PA, Doffek KM, et al. Increased loss of chromosome 9p21 but not p16 inactivation in primary non-small cell lung cancer from smokers. Cancer Res 2001;61:2092–6.[Abstract/Free Full Text]
- Kraunz KS, Nelson HH, Lemos M, Godleski JJ, Wiencke JK, Kelsey KT. Homozygous deletion of p16INK4a and tobacco carcinogen exposure in nonsmall cell lung cancer. Int J Cancer 2006;118:1364–9.[CrossRef][Medline]
- Tanaka R, Wang D, Morishita Y, et al. Loss of function of p16 gene and prognosis of pulmonary adenocarcinoma. Cancer 2005;103:608–15.[CrossRef][Medline]
- Chen JT, Chen YC, Wang YC, Tseng RC, Chen CY, Wang YC. Alterations of the p16(ink4a) gene in resected nonsmall cell lung tumors and exfoliated cells within sputum. Int J Cancer 2002;98:724–31.[CrossRef][Medline]
- Schouten JP, McElgunn CJ, Waaijer R, Zwijnenburg D, Diepvens F, Pals G. Relative quantification of 40 nucleic acid sequences by multiplex ligation-dependent probe amplification. Nucleic Acids Res 2002;30:e57.[Abstract/Free Full Text]
- Worsham MJ, Chen KM, Tiwari N, et al. Fine-mapping loss of gene architecture at the CDKN2B (p15INK4b), CDKN2A (p14ARF, p16INK4a), and MTAP genes in head and neck squamous cell carcinoma. Arch Otolaryngol Head Neck Surg 2006;132:409–15.[Abstract/Free Full Text]
- Mistry SH, Taylor C, Randerson-Moor JA, et al. Prevalence of 9p21 deletions in UK melanoma families. Genes Chromosomes Cancer 2005;44:292–300.[CrossRef][Medline]
- Sobin LH, Wittekind CH, editors. TNM classification of malignant tumours. 6th ed. New York: Wiley-Liss; 2002. p. 99–103.
- Noguchi M, Morikawa A, Kawasaki M, et al. Small adenocarcinoma of the lung. Histologic characteristics and prognosis. Cancer 1995;75:2844–52.[CrossRef][Medline]
- Matsumoto S, Iwakawa R, Kohno T, et al. Frequent EGFR mutations in noninvasive bronchioloalveolar carcinoma. Int J Cancer 2006;118:2498–504.[CrossRef][Medline]
- Sakamoto H, Mori M, Taira M, et al. Transforming gene from human stomach cancers and a noncancerous portion of stomach mucosa. Proc Natl Acad Sci U S A 1986;83:3997–4001.[Abstract/Free Full Text]
- Matsumoto S, Takahashi K, Iwakawa R, et al. Frequent EGFR mutations in brain metastases of lung adenocarcinoma. Int J Cancer 2006;119:1491–4.[CrossRef][Medline]
- Hamada K, Kohno T, Takahashi M, et al. Two regions of homozygous deletion clusters at chromosome band 9p21 in human lung cancer. Genes Chromosomes Cancer 2000;27:308–18.[CrossRef][Medline]
- Riely GJ, Politi KA, Miller VA, Pao W. Update on epidermal growth factor receptor mutations in non-small cell lung cancer. Clin Cancer Res 2006;12:7232–41.[Abstract/Free Full Text]
- Sakamoto H, Shimizu J, Horio Y, et al. Disproportionate representation of KRAS gene mutation in atypical adenomatous hyperplasia, but even distribution of EGFR gene mutation from preinvasive to invasive adenocarcinomas. J Pathol 2007;212:287–94.[CrossRef][Medline]
- Travis WD, Brambilla E, Muller-Hermelink HK, Harris CC, editors. Pathology and genetics: tumours of the lung, pleura, thymus and heart. Lyon: IARC Press; 2004. p. 73–5.
- Aoyagi Y, Yokose T, Minami Y, et al. Accumulation of losses of heterozygosity and multistep carcinogenesis in pulmonary adenocarcinoma. Cancer Res 2001;61:7950–4.[Abstract/Free Full Text]
- Haneda H, Sasaki H, Shimizu S, et al. Epidermal growth factor receptor gene mutation defines distinct subsets among small adenocarcinomas of the lung. Lung Cancer 2006;52:47–52.[CrossRef][Medline]
This article has been cited by other articles:

|
 |

|
 |
 
H. Nakanishi, S. Matsumoto, R. Iwakawa, T. Kohno, K. Suzuki, K. Tsuta, Y. Matsuno, M. Noguchi, E. Shimizu, and J. Yokota
Whole Genome Comparison of Allelic Imbalance between Noninvasive and Invasive Small-Sized Lung Adenocarcinomas
Cancer Res.,
February 15, 2009;
69(4):
1615 - 1623.
[Abstract]
[Full Text]
[PDF]
|
 |
|