| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
Molecular Biology, Pathobiology, and Genetics |
1 British Columbia Cancer Agency/Cancer Research Center and 2 Faculty of Dentistry, University of British Columbia, Vancouver, British Columbia, Canada
Requests for reprints: Martial Guillaud, Cancer Imaging Department, British Columbia Cancer Agency/Research Centre, 675 West 10th Avenue, Vancouver, BC, V5Z1L3, Canada. Phone: 604-675-8086; Fax: 604-675-8099; E-mail: mguillau{at}bccrc.ca.
| Abstract |
|---|
|
|
|---|
| Introduction |
|---|
|
|
|---|
Alternative approaches that facilitate the identification of high-risk OPLs, particularly those within hyperplasia or mild/moderate dysplasia, need to be developed. The foremost development in such approaches over the last 5 to 10 years has been aimed at finding molecular markers (7–9). Among these, microsatellite analysis of loss of heterozygosity (LOH) has shown promise as a powerful adjunct tool for risk prediction (7, 10–14). Retrospective studies from our research team have shown that LOH pattern can classify morphologically indistinguishable hyperplasia or mild/moderate dysplasia OPLs into different risk categories. OPLs without LOH at 3p and 9p have the lowest risk of progression. In comparison those with LOH at 3p and/or 9p but not in other arms have a 3.8-fold increase in relative cancer risk. Those with LOH at 3p and/or 9p plus additional losses (at 4q, 8p, 11q, or 17p, high risk) have a 33-fold increase in relative cancer risk (14).
Given that the largest proportion of OPLs are hyperplasia or mild/moderate dysplasia, there is a need for an inexpensive system to triage such samples for molecular analysis. There is some indication that imaging devices that characterize phenotype changes associated with cancer risk could play such a role by complementing and adding to our current knowledge of risk assessment (15–17). Moreover, because exploration of phenotype provides a window to the functional status of the cell, showing whether the genetic alterations within are associated with structural/functional changes, this approach could expand our understanding of the biology underlying malignant progression.
Histopathologic evaluation of presence and degree of dysplasia is an appraisal of the tissue (and cell) phenotype using multiple variables. Quantitative tissue pathology (QTP) breaks down the components of this phenotype into multiple quantifiable units that can be studied independently and in combination, allowing the investigator to examine associations of such change with progression risk (Fig. 1A and B ). For example, increased proportion of heterochromatin condensation is a high-risk nuclear characteristic for cancer (18). Traditional pathology judges this by one criterion, hyperchromatism, whereas QTP measures multiple chromatin features, such as whether the increased DNA is distributed around the edge of the nucleus or clustered in the center (Fig. 1B, left), whether the nucleus is dark with light areas or light with dark areas (Fig. 1B, middle), and whether the increased chromatin is evenly distributed (euchromatin) or clumped locally (heterochromatin; Fig. 1B, right).
|
| Materials and Methods |
|---|
|
|
|---|
The training set contained two groups of specimens representing the two extreme ends of the intraepithelium spectrum: 30 oral biopsies with areas of relatively normal oral mucosa (e.g., amalgam tattoo or melanotic macule) and a group of 29 SCC. Of the 29 SCC, the majority were from the high-risk oral site (floor of mouth and lateroventral tongue, 66%), were well to moderately differentiated (83%), and were stage 1 or stage 2 tumors (69%).
The test set contained 69 biopsies that can be separated into two subgroups: a group of 25 biopsies diagnosed as severe or CIS which was denoted as high-histology risk group because they are histologically at risk for progression (1, 2) and a group of 44 biopsies which were a subset of 83 well-characterized OPLs with sufficient material for QTP used in our previous study that explored the use of LOH to predict risk of progression for OPLs with hyperplasia or mild/moderate dysplasia (14). These 44 lesions (one biopsy per patient) diagnosed as hyperplasia, mild or moderate dysplasia are denoted the low-histology risk group. They were from patients without a prior history of head and neck cancer. The only prerequisite for inclusion in the present study was the availability of two adjacent unstained slides that were serial to slides used for the previous LOH analysis. Of these lesions, 15 have progressed to cancer and the remaining has not. Among the 15 progressing lesions, there were six hyperplasia, five mild dysplasia, and four moderate dysplasia. Among the nonprogressing lesions, there were 15 hyperplasia, 8 mild dysplasia, and 6 moderate dysplasia. There was no difference between the progressing and the nonprogressing lesions in terms of gender, age distribution, and smoking history (all P < 0.05). On average, the nonprogressing cases have been monitored for over twice the duration of progressing cases (71 mo versus 30 mo) to ensure that progression did not occur.
Sample preparation. All 128 samples were formalin fixed and paraffin embedded. The histologic diagnoses of the samples were confirmed by two oral pathologists (C.P. and L.Z.). Serial sections, 4-µm in thickness, were cut from each sample and placed on two glass slides; one was stained with H&E and the other with Feulgen-Thionin (19). The area on the H&E slide selected by the pathologist (C.P. or L.Z.) for conventional histopathologic diagnosis was the exact same area that was circled on the Feulgen-stained slide to be examined with the QTP imaging system.
QTP imaging system. This system is a modified version of the Cyto-Savant automated quantitative system (Cancer Imaging, BC Cancer Agency; ref. 19). The illumination wavelength was 600 ± 5 nm, corresponding to the absorption peak of the Thionin stain. The effective pixel sampling space within the plane of the sample was 0.34 µm2, and the effective pixel sampling area was 0.116 µm2. The software is specifically designed for interactive semiautomatic cellular and architectural analysis of tissue sections (19, 20). A strict quality control procedure was implemented to ensure the stability of the imaging system for each analysis (21). The imaging system characteristics followed the recommendation of the European Society of Analytical Cellular Pathology (22).
QTP image analysis. An experienced cytotechnician under the guidance of an oral pathologist (C.P. or L.Z.) delineated the region of interest (ROI) from the Feulgen-stained slide under a 20x objective. The area selected represented the most abnormal tissue according to the pathologist. The analysis of collected nuclei in each biopsy was done automatically by the system and consisted of four steps: (a) delineation and focusing of the ROI (Fig. 1A), (b) automatic thresholding and segmentation of the nuclei in the field (23), (c) interactive correction of segmentation errors and selection of only nonoverlapping cells, and (d) automatic collection of individually focused images of each selected cell (Fig. 1D). For normal epithelium and OPLs, the full width of the epithelium had to be present. Only in-focus, intact nuclei were collected by the cytotechnician with nuclei obtained from all layers of the epithelium. Basal and parabasal cells may be preferentially chosen since there are more of them, but they are also more likely to overlap or be in clusters (higher density) and so will be rejected from the analysis as being overlapping cells. Intermediate and superficial cells, although less numerous, are more often isolated, correctly segmented, and thus more often selected for analysis.
The regions used in each case of SCC were selected using three criteria: (a) the region of tumor selected needed to show tumor differentiation that was representative of the case; (b) the region had to have minimal amounts of inflammation and tissue; and (c) the tumor region needed to be reasonably large, i.e., not a single line or layer of tumor cells.
On average, >100 cell nuclei were collected from the delineated area for each biopsy. For each nucleus, 110 nuclear phenotype features were calculated by the system, measuring the size, the shape, the amount, and the distribution of stained DNA in the nucleus (24).
Calculation of NPS. To determine nuclear features that are predictors for OPLs with different cancer risk, we first trained the imaging system to separate normal specimens (30 samples) from cancer specimens (29 SCC samples) in the training set. This was done by establishing first a cell-by-cell NPS, which is a linear discriminant function of the nuclear features that are statistically significant, not collinear, that taken together, used as a predictor, can classify cells into separate categories. For this process, all epithelial cells (4,027 cells) from the 30 normal samples were grouped together and used to define the "normal-like" cell group, and all epithelial cells (4,298 cells) from the 29 SCC samples were grouped together to define the "cancer-like" cell group. A forward stepwise discriminant function analysis of the 110 nuclear phenotypic features was done to select which of these features could be used to generate a discriminant function that best separated the two groups of cells. All features were considered as continuous variables. We used the default values provided by STATISTICA package: tolerance value, 0.01; F to enter, 1.0; F to remove, 0.0. In stepwise discriminant function analysis, a model of discrimination is built step-by-step, automatically done by the software. Specifically, at each step, all variables are reviewed and evaluated to determine which one will contribute most to the discrimination between groups. That variable will then be included in the model, and the process repeats. The process stops when the resulting correct classification rate (percentage of cells correctly classified in each group) no longer increased with additional variable.
An interpretation of this discriminant function is that it represents a cell-by-cell phenotype classification procedure, which assigns a score to each nucleus indicating similarity to normal or cancer cells.
Similar to the method described in a companion paper (15), the range of this cellular score was divided into 10 continuous regions: the first region A represents values of phenotypically normal-like cells, and region J represents values of phenotypically cancer-like cells.
From the percentage of cells in each of these 10 regions, the NPS assigned to each biopsy was calculated as NPS = 1*A + 2*B + 3*C + 4*D + 5*E + 6*F + 7*G + 8*H + 9*I + 10*J, wherein A is the proportion of cells whose disciminant score falls in the region A, B is the proportion of cells whose disciminant score falls in the region B, and so forth.
This NPS represents the weighed sum of the 10 regions. The possible values for the NPS range from 1 to 10. An NPS with a value of 1 corresponds to a specimen containing 100% of normal-like cells, whereas an NPS with a value of 10 corresponds to a specimen containing 100% of cancer-like cells.
Statistical analysis. For comparison of difference between groups, a nonparametric Mann-Whitney U test was done. For comparison of difference between more than two groups, a nonparametric Kruskal-Wallis ANOVA was done. When the test showed a statistical difference across groups, a multiple comparison of mean ranks was done to confirm significant differences across group medians (multiple comparisons of mean ranks for all groups). Differences between group frequencies were assessed with the Yates corrected
2 test. Time-to-progression curves were estimated by the Kaplan-Meier method, and comparisons were done using the log-rank test. Relative risks were determined using Cox regression analysis. All analyses were done with Statistica6 (StatSoft, Inc.). All P values were two-sided. A two-sided P value of <0.05 was considered statistically significant.
| Results |
|---|
|
|
|---|
|
Correlation of NPS with pathology. Figure 3 shows the association between NPS and the pathology diagnoses. There was a monotonic increase in NPS with severity of pathology diagnosis [NPS median (range)]: 3.7 (2.01–5.03) for normal; 4.14 (3.03–7.90) for the low-histology risk group, i.e., the hyperplasia or mild/moderate dysplasia group; 5.24 (3.65–8.03) for high-histology risk group (severe dysplasia or CIS); and 7.25 (3.77–9.96) for SCC. There is a significant increase in NPS between the low-histology risk group (hyperplasia or mild/moderate dysplasia) and the high-histology risk group (P = 0.036). There are significant differences between all logical histologic groupings except normal versus hyperplasia or mild/moderate dysplasia and severe dysplasia/CIS. The increase in NPS between the severe dysplasia/CIS group and the SCC was close to significance (P = 0.08).
|
Correlation of NPS with progression to cancer. We sought to determine whether NPS would identify hyperplasia or mild/moderate dysplasia at high risk of developing into invasive SCC by comparing NPS values in progressing and nonprogressing lesions. Progressing cases showed significantly higher NPS value compared with nonprogressing lesions (NPS median, 5.7 for progressing lesions versus 3.8 for nonprogressing lesions; P < 0.0001). As shown in Fig. 4A , although the sample set was small, this trend remained statistically significant when the hyperplasia and dysplasia were examined separately [NPS median, 4.9 for progressing hyperplasia versus 3.8 for nonprogressing hyperplasia (P = 0.016) and 6.3 for progressing low-grade dysplasia versus 3.8 for nonprogressing low-grade dysplasia (P = 0.001)].
|
4.5 (high NPS) whereas 28 of the 29 SCCs (97%) had high NPS. When the hyperplasia or mild/moderate dysplasia lesions were categorized by this cutoff value, 17 of the 44 hyperplasia or mild/moderate dysplasia test study cases had high NPS. Thirteen of the 17 cases (76%) with high NPS progressed into invasive SCC in contrast to only 2 of the 27 cases (7%) with low NPS (P < 0.0001). In other words, with this cutoff, we could correctly classify 86% of nonprogressing cases (25 of 29) and 86% of progressing cases (13 of 15).
Time-to-progression curves were plotted as a function of the NPS classification into one of the two groups (high NPS versus low NPS; Fig. 4C). A significant difference was observed between the low-NPS group and the high-NPS group (P < 0.00015).
The proportion of lesions progressing to cancer at 5 years was 71% (12 of 17) for lesions with high-NPS and only 22% for lesions with low NPS. There was a 10-fold (relative risk, 10.3; confidence interval, 2.9–60.0) increase in the relative risk of progression to cancer for oral lesions with high NPS compared with those with a low NPS.
Correlation of NPS with LOH. Table 1 shows association between NPS and different LOH patterns, indicating the proportion of cases within each LOH pattern that are above and below the NPS cutoff. NPS values were consistently higher for samples with LOH at all seven chromosome regions examined. This increase was significant for 3p (P = 0.03), 4q (P = 0.003), 9p (P = 0.003), and 11q (P = 0.05). We also looked at associations with multiple losses and for LOH at 3p and/or 9p plus LOH at any of the arms 4q, 8p, 11q, 13q, and 17p, patterns previously identified as having markedly increased risk of progression (9, 11). NPS was strongly associated with the presence of these high-risk LOH patterns: multiple losses (P = 0.0005) and LOH at 3p and/or 9p plus LOH at any of the arms 4q, 8p, 11q, 13q and 17p (P < 0.0005).
|
Regression analysis with multiple covariates. We did a regression analysis with multiple covariates for the hyperplasia or mild/moderate dysplasia group, looking at histology, LOH, and NPS with cancer progression as primary end point. For each of these covariates, we assigned a value of either 0 or 1, denoting either low-risk or high-risk marker value, respectively. For histology, hyperplasia was coded 0 and mild or moderate dysplasia was coded 1. For LOH, cases with LOH at 3p and/or 9p plus LOH at any of the arms 4q, 8p, 11q, 13q, and 17p were coded 1 and 0 if otherwise. For NPS, cases with high NPS (NPS, >4.5) were coded 1 and cases with low NPS were coded 0. The results of the Cox proportional hazard regression are given in Table 2 . In this table, the first three models examine NPS, LOH, and histology separately and show that only LOH and NPS are significant in predicting cancer (with similar P value, P < 0.001). For hyperplasia or mild/moderate dysplasia, histology is not a predictor of cancer progression (P = 0.88). Categorization of hyperplasia with mild dysplasia versus moderate dysplasia was also not significant (data not shown). When LOH and NPS were entered together into the Cox regression analysis (model 4), NPS is no longer as significant as it was when used as a single predictor (model 1, P = 0.0007), although the borderline P value of 0.07 does suggest that a larger sample size may be required for a definitive conclusion. The different models clearly show that conventional histology has no predictive effect and that NPS is a better "phenotypic" predictive marker for cancer risk than histology (models 5 and 6).
|
| Discussion |
|---|
|
|
|---|
Our results showed that the in-house imaging system developed by our research team, although not yet fully automated and high throughput, is quantitative and can unravel subtle differences in OPLs. We used two well-defined groups of cells (normal and cancer) to generate an NPS for each specimen. This system was able to analyze 110 nuclear features and identify not only differences in nuclear size and chromatin density as noted by pathologists, but also differences in nuclear shape and chromatin texture, features not readily quantified by pathologists. The identification of different features of chromatin texture changes by the system is significant because such changes are now widely accepted as an indication of genetic or epigenetic changes (18). Furthermore, the system also quantified the changes measured, allowing an objective computation of an NPS for each sample.
In a previous study (15), we determined how sensitive the NPS process is to the interactive steps involved in the collection of the nuclei used to calculate NPS. In that study, 30 preneoplastic bronchial biopsies which spanned the ranges of pathology grades and NPS were selected and measured. All samples were measured twice by one experienced cytotechnologist and twice by another experienced cytotechnologist. The measurement of the NPS was highly reproducible between individual technologists. The coefficient of correlation between two histopathology technologists was 0.98. The intratechnologist variability was comparable.
The NPS was noted to continuously increase with severity of histologic grading (P < 0.0001). A significant difference was observed between the low-histology risk hyperplasia or mild/moderate dysplasia group and the high-histology risk group (severe dysplasia/CIS). In British Columbia, a histologic diagnosis of severe dysplasia/CIS results in the majority of cases receiving treatment by surgery. Thus, we were unable to examine a potential association of NPS in progression within this category. However, 72% of the severe dysplasia/CIS in this study had NPS of >4.5 indicating a strong association of NPS with this high-risk histologic category.
This is the first study to investigate the relationship between NPS and outcome for hyperplasia or mild/moderate dysplasia OPLs, which is ultimately the most important aspect of validation of any new test. High NPS was found to be strongly associated with cancer progression (P < 0.0001). This association was independent of histology because high NPS in oral lesions with low-grade or no dysplasia was still strongly associated with cancer progression, although the number of cases were small. These results were also assessed by multivariate Cox regression analysis and showed an absence of predictive value for histology classification among hyperplasia or mild/moderate dysplasia lesions. In a mixed model, LOH is still the most significant predictor, but a larger sample size is required to determine whether QTP adds significantly to the correlation in a mixed model. It does support the use of QTP for triaging samples that require LOH analysis.
This is the first study to investigate the relationship between NPS and molecular changes of premalignant lesions. We used the microsatellite analysis for LOH, a technique currently used by many research groups in the study of cancer risk of OPLs. Increased NPS was significantly associated with the presence of previously reported high-risk LOH patterns, including the presence of multiple losses and LOH at 3p and/or 9p plus additional losses (at 4q, 8p, 11q, or 17p; refs. 11–14).
We have determined an NPS cutoff value of 4.5 based on the training set (normal and SCC; Fig. 2). Seventeen of the 34 hyperplasia or mild/moderate dysplasia cases showed high NPS. Sixteen of these 17 cases had LOH data available, and all showed LOH at 3p and/or 9p, a loss believed to be essential although not sufficient for cancer progression (14). Eleven of these 16 cases progressed into invasive SCC. Individuals with high NPS have a 10-fold increase in relative risk of developing a cancer within 5 years compared with individuals with a low NPS (confidence interval, 2.05–46.6). Similar results were obtained when the NPS threshold used to define the high NPS and the low NPS groups was varied.
NPS could therefore offer clinicians additional information on the possible malignant potential of lesions. For example, 6 of the 21 hyperplasia had high NPS, five of these six had a molecular pattern available, and four of them showed the high-risk LOH pattern. All three of the hyperplasia that progressed into cancer had these characteristics. Of the 23 low-grade dysplasia, 11 (48%) had high NPS and 9 of these had LOH at 3p and/or 9p (7 had LOH at 3p and/or 9p plus another loss). Eight of the 11 progressed into cancer. These results suggest that the QTP image system could be a powerful adjunct tool and reproducible method to identify hyperplasia or mild/moderate dysplasia lesions for molecular assessment.
These findings are similar to results obtained by our team in bronchial preneoplastic lesions (15). Using the same approach as the one described in this paper, we showed that QTP correlates better with genetic damage and cancer progression than conventional pathology. Bronchial NPS is currently being used as a surrogate end-point biomarker in a chemoprevention trial in bronchial dysplasia (27). Similarly, using different morphometric features, quantitative analysis of endometrial hyperplasia correlates better with genetic clonality than conventional pathology (28). Similar results by other groups have been found in bladder, prostate, cervix, and breast cancer for the association of QTP with histologic grading (29–31). Regarding the nuclear features used in the calculation of our NPS, they are similar to the features used in other studies for the calculation of morphometric scores, such as quantitative nuclear grade by Veltri (29) in transitional cell carcinoma, discriminant score by Hanselaar et al. (32) in cervical cytologic specimens, and morphometric nuclear grade by Bacus et al. (16) in breast and cervical specimens. The main consistency between all these studies, different in terms of tissue types and specimens (cytology/histology), is that they all used a combination of different types of features, i.e., features describing the nuclear shape and size, as well as features measuring chromatin texture. Furthermore, nuclear feature calculations can easily be reproduced as described in Rodenacker (33).
Recently, Montinori et al. (17) has shown that these findings could be reproduced using H&E-stained material and fully automated analysis to predict recurrence in urothelial lesions. If this work is confirmed, it could be a definitive breakthrough for the implementation of quantitative histology in routine histopathology, as both conventional and quantitative analysis would be done on the same material, with the analysis fast and automated.
A large longitudinal study is needed to confirm the prognostic value of QTP, and such a longitudinal study is under way. The overall goal of that study is to evaluate three new innovative, potentially complementary devices for early detection and follow-up of OPLs and oral cancer and to determine their potential role in the triage of patients in a provincial network that is being established in British Columbia (34). One device is an in vivo clinical tissue evaluation tool; the others involve the development of quantitative pathology approaches.
In conclusion, the results presented in this paper show that NPS correlates strongly with genetic damage, cancer progression, and histology grading and could have a potential value for triaging OPL for molecular analysis.
| Acknowledgments |
|---|
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.
We thank Anita Carraro for her expert technical support and Dr. Michelle Follen (M. D. Anderson Medical Center) for her valuable criticism.
Received 6/ 6/07. Revised 12/28/07. Accepted 2/22/08.
| References |
|---|
|
|
|---|
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
| Cancer Research | Clinical Cancer Research |
| Cancer Epidemiology Biomarkers & Prevention | Molecular Cancer Therapeutics |
| Molecular Cancer Research | Cancer Prevention Research |
| Cancer Prevention Journals Portal | Cancer Reviews Online |
| Annual Meeting Education Book | Meeting Abstracts Online |