<title>Introduction</title> <p>The proportion of breast cancers currently being diagnosed in the <italic>in situ</italic> stage (ductal carcinoma in situ, DCIS) has been increasing with the increased use of mammographic screening and now accounts for up to about 25% of breast cancers diagnosed in some centers. Nuclear grade is a major determinant of therapeutic decisions although there is no single accepted grading system (##REF##9744307##Badve et al. 1998##; ##REF##7831528##Holland et al. 1994##; ##REF##15199110##Leonard and Swain, 2004##; ##REF##11345835##Leong et al. 2001##; ##REF##8635094##Silverstein et al. 1996##; ##REF##10774450##Silverstein, 2000##; ##REF##10088542##Sneige et al. 1999##; ##REF##9504685##Tavassoli, 1998##; ##REF##1317200##van Dongen et al. 1992##; ##REF##12160327##Wärnberg et al. 2002##; ##REF##10505030##Wärnberg et al. 1999##).</p> <p>Nuclear grade has been associated with DCIS recurrence (##REF##9744307##Badve et al. 1998##; ##REF##15070793##Burstein et al. 2004##; ##REF##2536582##Lagios et al. 1989##) and progression to invasive carcinoma (##REF##2536582##Lagios et al. 1989##; ##REF##14625260##Kerilkowske et al. 2003##). Based on such results, nuclear grade has been recommended for the primary stratification of DCIS (Schwarz, 1997). However, the association between nuclear grade and recurrence of DCIS or development of invasive carcinoma is not as clear in some studies in which as many as 50% of patients are reported to have more than one nuclear grade (##REF##11906438##Miller et al. 2001##). In ##REF##11906438##Miller et al. 2001## attempts were made to reflect and evaluate heterogeneity by evaluating both most dominant (most prominent) grade and worst grade; however, nuclear grading was not associated with clinical outcome. In this context, a single grade may not be sufficient to characterize a patient and this may affect the medical therapeutic decisions being made. It is important to recognize and account for heterogeneity of nuclear grade within each patient since as many as 50% of patients with DCIS may have more than one nuclear grade (##REF##11906438##Miller et al. 2001##; ##REF##8816587##Goldstein and Murphy, 1996##).</p> <p>It would be advantageous to have a method to obtain reliable quantitative data that could be objectively analyzed to assist with the classification of patients who have breast DCIS in which more than one nuclear grade is present. Computer-aided image analysis is such a method. Image analysis has been used to extract quantitative nuclear information useful for diagnosis of biopsy specimens of many tissues, including breast DCIS, and invasive breast (##REF##14625260##Kerlikowske et al. 2003##;Carpenter et al. 1985; ##REF##11147303##Dey et al. 2000##; ##REF##11777277##Frank et al. 2001##; ##UREF##3##Hoque et al. 2001##; ##REF##8832747##Mariuzzi et al. 1996##; ##REF##11169513##Mommers et al. 2001##; ##REF##7734410##Susnik et al. 1995##; ##REF##8978872##Tuczek et al. 1996##; ##REF##7628853##Wolberg et al. 1995##). However, these previous reports did not focus on evaluating information for individual nuclei in the classification of patients with mixed nuclear grades. We previously examined the role of image analysis of individual nuclei for the ##REF##11906438##Miller et al. 2001## cohort with respect to the development of invasive breast cancer, and found that the reproducible, quantitated image analysis features were associated with development of invasive disease (##REF##17845726##Chapman et al. 2007##) whereas nuclear grade was not. In that study, we determined that indications for particular image features across fields were not consistent when there was heterogeneous grading.</p> <p>The purpose of this study was to determine, if after considering clinical and pathologic factors, whether nuclear features measured by digital image analysis of hematoxylin and eosin stained slides would be significantly associated with DCIS recurrence, and therefore could help pathologists classify patients whose DCIS was more likely to recur. The clinical cohort of patients whose DCIS exhibited substantive heterogeneity (##REF##11906438##Miller et al. 2001##) was used for these investigations. We restricted the examinations to pooled assessments across fields, and an overall patient basis, for better per patient representation of nuclear features.</p>
<title>Patients and Methods</title> <title>Patients and specimens</title> <p>The DCIS patients whose specimens were used in this study were a subset of the 124 previously described in detail (##REF##11906438##Miller et al. 2001##; ##REF##9869520##Fish et al. 1998##; ##REF##17845726##Chapman et al. 2007##). Of the 124, the cohort of 88 patients who underwent lumpectomy alone were used, as these patients had experienced most of the events: 17 of 19 recurrences of DCIS with median 5.0 years of follow-up. Of these 88 patients, three did not have slides available for assessment and for five patients the hematoxylin and eosin (H&amp;E) stained slides resulted in images of poor quality, unsuitable for further evaluation. Digital images of the slides of the 80 study patients were acquired under the supervision of a breast pathologist (NM). Worst nuclear grade of 1, 2, or 3, by the Van Nuys system (##REF##8635094##Silverstein et al. 1996##; ##REF##7723550##Silverstein et al. 1995##) was assigned by the pathologist by viewing the slides directly in a bright field microscope, yielding 1 patient with grade 1, 31 with grade 2, and 48 with grade 3 DCIS.</p> <p>The focus for this study was an investigation of the effects of heterogeneity on the assessment of image analysis features. When more than one grade was present, predominant grade, and all grades were recorded for each of two representative fields for each patient. There were 32 patients in Group A, whose DCIS was nuclear grade 1 only (1 patient), grades 1 and grade 2 (8 patients), or grade 2 only (23 patients); 31 patients in Group B with grades 1, 2, and 3 (2 patients), or grades 2 and 3 (29 patients); and 17 patients in Group C with grade 3 only.</p> <title>Segmentation and preparation of image analysis data</title> <p>As described above, DCIS heterogeneity was observed pathologically with H&amp;E slides. Digital images of areas of DCIS identified by the pathologist were captured with a CCD camera, bright field microscope, desktop computer, and NIH-Image software, as described in detail in the <xref ref-type="app" rid="app1-cin-6-0099">Appendix</xref>. The experimental design involved imaging five ducts in one field, followed by imaging five ducts in a second field. Those fields demonstrated DCIS of the nuclear grade(s) recorded for that patient. The fields were present as two separate fields on one slide or one field on each of two slides. Thus, ten or more images in different regions of the slides(s) were captured for each patient, 20 or more nuclei per image were segmented resulting typically in 200 nuclei analyzed for each of the 80 patients. The focus for these investigations are the two fields with about 100 nuclei assessed per field and the pooled data for both fields, to yield an overall assessment for a patient. A blank field was subtracted from each captured image. Images were printed and nuclei labeled with an index number to avoid duplication. Contrast enhanced images were viewed on a 17 inch monitor and enlarged 4:1. These changes did not affect the values extracted from the saved image. Each nuclear region of interest was segmented by an operator using a computer mouse. The operator was blinded to the grade of the specimen determined by the pathologist. All nuclear images were segmented by a single person (DA). The reproducibility of feature values extracted by operator guided segmentation is discussed in the Results section. Some image feature calculations and the merging of all nuclear image feature data to per field and per patient attribution were accomplished with StatView v 5.01 (Brain Power, Calabasas, CA, USA) software.</p> <title>Image analysis features</title> <p>For each nucleus, 39 features were determined in three categories. (i) Morphometry: area, perimeter, ellipse major axis, ellipse minor axis, ellipticity (major axis/minor axis), shape form factor (4 X pi X area/perimeter squared), and roundness b (4 X area/pi X ellipticity squared) (##UREF##6##Russ, 1990##). (ii) Densitometry: mean density, standard deviation of density, modal density, minimum density, maximum density, sum density (mean density X area, used instead of I.O.D. of NIH-Image), range density. (iii) Markovian texture features (##UREF##2##Haralick et al. 1973##; ##REF##56387##Pressman, 1976##) were calculated from the co-occurrence matrix of pixel densities with a step size of 2. They were angular second moment, contrast, correlation, variance, inverse difference moment, sum average, sum variance (corrected from ##REF##56387##Pressman, 1976##), difference average, difference variance, initial entropy, final entropy, entropy, sum entropy, difference entropy, coefficient of variation, peak transition probability, diagonal variance, diagonal moment, second diagonal moment, product moment, and triangular symmetry. Additional texture features, calculated from the binned histogram of pixel gray scale values, included histogram mean, histogram variance, histogram skewness, and histogram kurtosis. Similar features have previously been used to classify breast ductal carcinoma in situ specimens (##REF##11777277##Frank et al. 2001##; ##UREF##3##Hoque et al. 2001##).</p> <title>Prognostic factors</title> <p>The clinical factors recorded on these patients were age (in years) and type of presentation (mammographic, clinically palpable, bloody nipple discharge). The histologic factors previously evaluated (##REF##11906438##Miller et al. 2001##) were maximum DCIS size (cm), percentage of parenchyma involved with DCIS (&lt;10%, 10%–50%, &gt;50%), predominant architecture (0—cribriform/micropapillary/other, 1—solid), worst architecture (0—cribriform/ micropapillary/other, 1—solid), nuclear grade [by the Van Nuys Classification system worst (nuclear grade 1, 2, 3); also, predominant (nuclear grade 1, 2, 3)] and heterogeneous nuclear grading (Group A, Group B, Group C)], necrosis [none, confluent (comedo-like)], calcification (none, crystalline/ amorphous), measured margin (zero margin, &lt;1 mm, 1–5 mm, &gt;5 mm), presence of uninvolved intervening duct (not assessable, no, yes), Van Nuys Prognostic Index. In addition, 39 nuclear image features were determined for about 200 nuclei per patient.</p> <p>For each patient, the image data were pooled across i) all nuclei in a field (2 assessments), and ii) all nuclei for a patient (1 assessment) to yield a summary feature value [adjusted mean = mean/ (standard error of the mean)], for each of the 39 image features for nuclei of the 3 different assessments per patient: 2 fields, 1 overall. In addition, grading discriminant classification functions, that are weighted combinations of image features, described below in the Analysis section, were assessed as prognostic factors.</p> <p>Three different assessments, corresponding to the 3 different ways of pooling the image analysis feature data, were performed to examine the effects of DCIS heterogeneity on apparent associations with clinical outcome. In other contexts, investigations have been restricted to single ducts, fields, or pooled per person assessments without an examination of replicability.</p> <title>Events</title> <p>Recurrence of ipsilateral DCIS made more than 90 days after the initial surgery was designated as an event. There were no developments of contralateral DCIS. There were no deaths from breast cancer, or another cause, in this group of patients over the study period.</p> <title>Statistical analysis</title> <p>Statistical analyses were performed with BMDP PC Dynamic Version 7.0 (same as BMDP-XP, Statistical Solutions, Sagua, MA, USA). Analyses included for each image feature and each patient, 1) Levene’s tests for equality of variance between fields for each person and between people, 2) use of the mean/S.E.M. of image features on a field and patient basis due to highly significant evidence against assumption of equal variances, 3) forward step-wise Fisher linear discriminant analyses, using an entry p-value of p = 0.05, and 4) jackknifed (leave-one-out) classification of patients to find the number of patients who were correctly classified by the discriminant functions.</p> <p>The histologic, clinical, and image analysis factors were assessed with respect to whether they were associated with ipsilateral DCIS recurrence. Univariate assessments were with Kaplan-Meier plots and the Wilcoxon (Peto-Prentice) test statistic (##REF##393312##Prentice and Marek, 1979##); for each image feature, standard image analysis cut-points at the means of the data were utilized after confirmation that the data were approximately symmetric.</p> <p>Multivariate assessments utilized continuous data where possible, and were with Cox forward stepwise regressions, using the likelihood ratio criterion (~χ<sup>2</sup><sub>(1),</sub> p ≤ 0.05) as the test statistic to determine if a factor would be added to the model. Since we had no knowledge of which of the image analysis features assessed would best reflect a patient’s DCIS, or the extent to which differences in image features might relate to prognosis, we performed 3 sets of multivariate analyses, corresponding to the 3 generations of image feature factors per patient: per 2 fields, 1 pooled across 2 fields.</p>
<title>Results</title> <title>Suitability of H&amp;E stained slides for image analysis</title> <p>Archival H&amp;E slides were satisfactory for image analysis for the majority of this DCIS cohort (##REF##17845726##Chapman et al. 2007##). Only 5.7% (5/88) patients had slides that were of too poor quality.</p> <title>Reproducibility of segmentation and image measurements</title> <p>All nuclear regions of interest were segmented manually by the same operator (DA). Reproducibility of manual segmentation was determined by repeatedly segmenting the same nucleus (CV = 3.4%, n = 150). In order to determine the reproducibility of extracted feature values, independent measurements were made of the same nucleus in images captured at different times. Ten or more nuclei, identified from images of specimens of seven patients, were segmented at two different times without knowledge of the previously segmented region of interest, or of the extracted feature values. The differences between the pairs of measurements of the same nucleus were determined with a two-tailed t-test. (The null hypothesis was that the differences were equal to zero, df = n-1, at the 0.05 level of significance). The percent of feature values that were not statistically different in duplicate measurements ranged from 89.7% to 97.4% among nuclei of the seven patients. Further, there were no statistically significant differences for 23 of the 39 feature values in 7/7 pairs of images, and 37 of 39 feature values in at least 5/7 pairs of images. Among the features whose values had minimal differences in all pairs of images there were two that were selected by discriminant analysis, the morphological feature ellipse minor axis and the texture feature peak transition probability.</p> <title>Classification of patients by image analysis</title> <p>In order to take into account patients with mixed grades, patients were classified as defined above into Groups A, B, C, corresponding to low, intermediate and high grades. These grading groups were used in discriminant analyses for nuclei in field 1 (80 patients), field 2 (79 patients), and overall pooled across both fields (80 patients). The results of the three discriminant analyses are provided in ##TAB##0##Table 1##. In each instance (field 1, field 2, and overall across both fields), there were image features significantly associated with the grading classifications, p &lt; 0.001 for each. The discriminant function for the first field included one morphological feature reflecting the size of nuclei (minor ellipse axis) and one texture feature reflecting the arrangement of DNA in the nucleus (sum entropy). Discriminant analysis of the second field included one morphological feature (perimeter), one densitometric feature (range density) and one texture feature (angular second moment). The analyses for both fields indicated one morphometric feature (minor ellipse axis) and one texture feature (peak transition probability). Different image analysis features were obtained for the first field, second field, and both fields. Correct classification of the nuclear grading with the image features for each field was respectively, 65%, 67.1% and 65.0%, in ##TAB##0##Table 1##.</p> <p>Discrimination using both fields would be most representative for a patient. A larger minor ellipse axis, indicative of a rounder nucleus, (p &lt; 0.001) and lower peak transition probability, indicating more uniform nuclear staining, (p &lt; 0.001) were associated with higher grading. ##TAB##1##Table 2## indicates the accuracy of image analysis classification by nuclear grading group. Image analysis features correctly classified 78.1% of patients in Group A (grade 1 and grade 2 nuclei), 48.4% in Group B (grade 2 and grade 3 nuclei), and 70.6% in Group C (grade 3 nuclei). The discriminant function was optimized to separate patients into the three nuclear grading groups. The distribution of patients within each of the three groups is illustrated in ##FIG##0##Figure 1##. For each patient the canonical variable = 1.00475 X (minor ellipse axis) – 0.60149 X (peak transition probability) – 1.7285, where the canonical variable has zero mean and coefficients standardized by pooled within group variance, F-statistic = 12.304, df = 4,152, p &lt; 0.001. Although the discriminant function was optimized to provide the best classification of patients within the three grading groups, there was considerable overlap between patients in the three grading groups. ##FIG##1##Figure 2## shows the discriminant function values ranked according to the value of the patient canonical variable; there is a nearly continuous distribution.</p> <p>##TAB##2##Table 3## indicates the factors significantly (p ≤ 0.05) associated with DCIS recurrence in the multivariate analyses based on data for the 2 fields and overall pooled assessments. Each model contains both clinical and image features; however, the image features differ. The overall pooled data would represent the best available summary for the patients. The factors significantly associated with recurrence of DCIS, were those previously found by ##REF##11906438##Miller et al. 2001##, type of initial presentation (p = 0.03), and amount of parenchymal involvement (p = 0.05), along with the morphometry image feature ellipticity, p = 0.04. ##FIG##2##Figure 3## shows the Kaplan-Meier plot for the image feature; smaller ellipticity (less elongated and more rounded nuclei) was associated with higher DCIS recurrence Van Nuys nuclear grade, predominant grade, and the grading discriminant function were not significantly associated with recurrence of ipsilateral DCIS.</p>
<title>Discussion</title> <p>Nuclear grade, or all grades when there was more than one grade, was reported for each patient by the pathologist according to the Van Nuys system (##REF##7831528##Holland et al. 1994##; ##REF##7723550##Silverstein et al. 1995##). This system has been compared to other systems that have been proposed to predict development of infiltrating carcinoma (##REF##9744307##Badve et al. 1998##; ##REF##11345835##Leong et al. 2001##; ##UREF##1##Douglas-Jones et al. 1996##). Several of these grading systems have recently been reviewed (##REF##15199110##Leonard and Swain, 2004##). The greatest consistency among pathologists seems to be obtained with systems that are based, in large part, on nuclear grade (##UREF##4##Morrow et al. 2002##) and a consensus conference on classification of DCIS recommended that DCIS should be stratified primarily by nuclear grade (##REF##9351550##Schwartz, 1997##). Our results indicate that image analysis can provide a reproducible quantitative description of nuclei for this purpose since duplicate measurements were similar for 89.7% to 97.4% of features assessed.</p> <p>We investigated the ability of 39 image features describing nuclear morphology (size and shape), densitometry (amount of stain), and texture (arrangement of DNA) to quantitatively discriminate tumors with pathologically determined nuclear grade. The most representative grading for a patient utilized data obtained for nuclei in both fields. Features included in the discriminant function were one whose value was determined by the size and shape of the nuclei (minor ellipse axis) and one whose value was determined by the arrangement of DNA in the nuclei (peak transition probability). Similar rates of accurate classification of grade were obtained from the first field assessed (about 100 nuclei), the second field (about 100 nuclei), and both fields (about 200 nuclei), with correct classifications of respectively 65.0%, 67.1%, and 65.0% of the patients. However, the discriminant functions for the three situations differed by image features. Discrimination using both fields would be the most representative of the patient’s grading.</p> <p>Discriminant analysis has frequently been used to classify patients based on feature values determined by computer-aided image cytometry (##UREF##0##Baak et al. 1991##; ##REF##6990845##Bartels, 1980##; ##UREF##5##Patterson, 1995##). The features selected and the weights assigned to each feature are readily interpretable. Further, we examined the relevance of the continuous discriminant functions here in a censored survival analysis framework.</p> <p>In this study, image features were extracted from H&amp;E stained slides rather than Feulgen stained slides. The Feulgen reaction is useful because it is specific for DNA, the intensity of the reaction is proportional to the amount of DNA, and nuclear regions of interest can be automatically segmented (##UREF##7##Schulte et al. 1995##). However, conventional H&amp;E stained slides have also been used to extract morphometric and other nuclear feature values (##REF##11777277##Frank et al. 2001##; ##REF##8978872##Tuczek et al. 1996##; ##REF##7628853##Wolberg et al. 1995##; ##REF##7507673##Christen et al. 1993##; ##REF##1655233##Peinta and Coffey, 1991##; ##REF##11002269##Weyn et al. 2000##). The possibility of extracting nuclear features from H&amp;E stained slides offers several advantages. First, existing H&amp;E stained slides available from pathology archives can be used without recovering the original paraffin blocks, recutting tissue sections, and staining the new slides with the Feulgen procedure; second, digital images of each microscopic field can be acquired by the image analysis system simultaneously with review by a pathologist viewing familiar H&amp;E stained tissue. Although the density measurements extracted from H&amp;E stained slides are not directly proportional to amount of DNA as they are for Feulgen stain, we have found that nuclear area is proportional to amount of DNA as determined by Feulgen (r = 0.886). Nevertheless, measurement of nuclear area in H&amp;E stained specimens is not an adequate substitute for measurement of DNA ploidy as determined with the DNA-specific Feulgen stain. In the future, the labor intensive chore of manually segmenting nuclear regions of H&amp;E stained slides might be overcome by acquiring color images and segmenting them with appropriate algorithms (Ferr-Roca et al. 1998; ##REF##14714298##Latson et al. 2003##). In our study only 5.7% of the archival slides had too poor stain quality for image analysis.</p> <p>Patients with breast DCIS that had intermediate nuclear grades, and/or more than one grade, offer challenges to pathologists. Such patients may represent as much as 50% of DCIS patients (##REF##11906438##Miller et al. 2001##). In this study computer-aided image analysis was used to assess a cohort of patients who presented with DCIS, many with intermediate or mixed grades. H&amp;E stained slides were viewed by a pathologist who recorded the nuclear grades in each field and digital images were acquired simultaneously. Nuclear image feature values were extracted and patients classified by discriminant analysis using the pathologist’s grouping of patients to supervise the analysis. Patients were placed into three groups according to the grades assigned by the pathologist taking into account the large proportion of patients with intermediate and/or mixed grades. The discriminant function correctly classify 78.1% of patients with low nuclear grade, and 70.6% of patients with high nuclear grade. The same discriminant function had a poorer success rate of only 48.4% for the intermediate group. This suggests that the intermediate group could be subtyped by characteristics different than those separating the low and high grade groups. Such subtypes of the intermediate grade group might be relevant for therapeutic decisions.</p> <p>Since high nuclear grade is one of the factors that has been associated with recurrence it might be expected that high nuclear grade would also be associated with recurrence within this cohort of patients. However, neither the highest nuclear grade, nor the most prominent nuclear grade, was associated with recurrence (##REF##11906438##Miller et al. 2001##; ##REF##17845726##Chapman et al. 2007##). The large proportion of patients with mixed nuclear grades may have contributed to the lack of association. In order to take into account the patients with mixed grades, discriminant analysis used the pathologist’s nuclear grade(s) for each patient, the grouping of patients with mixed grades, and the quantitative nuclear feature values determined by image cytometry. A discriminant function was derived that optimized the classification of patients into grading groups. The value of the canonical variable was determined for each patient. The canonical variable was standardized across patient values to have a mean of zero and a standard deviation of one. The distribution of patient values by canonical variable appeared to be continuous rather than three discrete and separate groups. The discriminant function was not associated with ipsilateral DCIS recurrence, although it had been with some pooled assessments concerning the development of invasive breast cancer (##REF##17845726##Chapman et al. 2007##). One of the image features included in the discriminant function was minor ellipse axis (p &lt; 0.001); it is noteworthy that a related feature, ellipticity, was associated with DCIS recurrence (p = 0.04). Rounder nuclei were associated both with higher grade and DCIS recurrence. It should be noted that this association was found with image feature adjusted means, after accounting for the greater variability observed with high nuclear grade. The image features alone would not be expected to be sufficient for prognosis. However, the additional information, added to other pathologic and molecular features, may improve the prognostic classification of patients.</p>
<title>Conclusion</title> <p>A discriminant function was derived that optimized the classification of DCIS patients with mixed nuclear grades. The classification of patients into grading groups was significant (p &lt; 0.001), although the separation of patient grading groups was not complete. Grouping patients with mixed nuclear grades and measuring their nuclear image features may contribute to their classification and prognosis.</p>
<title>Background</title> <p>Nuclear grade has been associated with breast DCIS recurrence and progression to invasive carcinoma; however, our previous study of a cohort of patients with breast DCIS did not find such an association with outcome. Fifty percent of patients had heterogeneous DCIS with more than one nuclear grade. The aim of the current study was to investigate the effect of quantitative nuclear features assessed with digital image analysis on ipsilateral DCIS recurrence.</p> <title>Methods</title> <p>Hematoxylin and eosin stained slides for a cohort of 80 patients with primary breast DCIS were reviewed and two fields with representative grade (or grades) were identified by a Pathologist and simultaneously used for acquisition of digital images for each field. Van Nuys worst nuclear grade was assigned, as was predominant grade, and heterogeneous grading when present. Patients were grouped by heterogeneity of their nuclear grade: Group A: nuclear grade 1 only, nuclear grades 1 and 2, or nuclear grade 2 only (32 patients), Group B: nuclear grades 1, 2 and 3, or nuclear grades 2 and 3 (31 patients), Group 3: nuclear grade 3 only (17 patients). Nuclear fine structure was assessed by software which captured thirty-nine nuclear feature values describing nuclear morphometry, densitometry, and texture. Step-wise forward Cox regressions were performed with previous clinical and pathologic factors, and the new image analysis features.</p> <title>Results</title> <p>Duplicate measurements were similar for 89.7% to 97.4% of assessed image features. The rate of correct classification of nuclear grading with digital image analysis features was similar in the two fields, and pooled assessment across both fields. In the pooled assessment, a discriminant function with one nuclear morphometric and one texture feature was significantly (p = 0.001) associated with nuclear grading, and provided correct jackknifed classification of a patient’s nuclear grade for Group A (78.1%), Group B (48.4%), and Group C (70.6%). The factors significantly associated with DCIS recurrence were those previously found, type of initial presentation (p = 0.03) and amount of parenchymal involvement (p = 0.05), along with the morphometry image feature of ellipticity (p = 0.04).</p> <title>Conclusion</title> <p>Analysis of nuclear features measured by image cytometry may contribute to the classification and prognosis of breast DCIS patients with more than one nuclear grade.</p>
<title>Acknowledgments</title> <p>We thank Dr. Edward B. Fish and Marilyn Link for their extensive work on the clinical database, Drs. D. August, Y. Gusev, R. Sklarew, D. Foran, and W. Hanna for helpful discussions, Dr. W. Sofer for equipment, and A. Khokhar and A. Kagan for laboratory assistance. This work was funded by the New Jersey Commission for Cancer Research 1076-CCR-S0, the National Institutes of Health U56 CA113004, the Hyde and Watson Foundation, the Busch Memorial Fund, and the E.B. Fish Research Fund.</p> <title>Appendix</title> <title>Image analysis system</title> <title>Hardware</title> <p>A custom image cytometry system was assembled which consisted of a CCD camera attached to a bright field microscope and linked to a desktop computer with a frame grabber card. Images of nuclei were acquired and stored as follows: hematoxylin and eosin stained slides were viewed with a bright field microscope (Wild model M20), 40X N.A. 0.75 objective, 1.25X phototube, 530–590 nm band pass green filter, detected with an 8 bit monochrome CCD camera (Sonyo model VDC3874) connected to a video monitor (RCA TC1112) and a frame grabber card (60 HZ Data Translation Quick Capture model DT2255) in a desktop computer (Apple Macintosh model IIci, 12 MB RAM, 80 MB hard disk), and stored as uncompressed TIFF files on removable media (Zip 100 disks). Ten frames were averaged and acquired using NIH-Image software (v. 1.57, written by Wayne Rasband, obtained from the internet by anonymous FTP). Each TIFF formatted image was 640 × 480 pixels, with 256 gray levels. The resulting pixel images were isotropic, with an effective size of 0.25 microns × 0.25 microns. Segregated nuclear images were of modest resolution, typically containing 800 to 1600 square pixels. Sizes were calibrated with a B&amp;L stage micrometer. Optical density of pixel gray values was standardized and camera response calibrated with a set of neutral density filters (50, 25, and 12.5% transmission).</p> <title>Software</title> <p>NIH-Image v.1.62b34-Arnv software (modified from <ext-link ext-link-type="uri" xlink:href="http://rsb.info.nih.gov/">http://rsb.info.nih.gov/</ext-link>) and StatView v. 5.01 statistical package (BrainPower, Calabasas, CA, USA) were used to measure and calculate DNA densitometric and nuclear morphometric features using a Mac G4 computer. TextureCalc v. 1.1ax, software (written by W. C.-B.) was used to rebin 256 gray levels into 8 intervals and to calculate texture features from the Markovian gray level co-occurrence matrices. Programs written in SAS release 6.12 for the Macintosh (SAS Institute, Inc., Cary, NC, USA) were used to format data and BMDP PC Dynamic Version 7.0 (Statistical Solutions, Sagua, MA, USA) statistical software was used for discriminant analysis.</p>
<fig id="f1-cin-6-0099" position="float"><label>Figure 1</label><caption><p>Distribution of patients between groups. <bold>a</bold>, Group A; <bold>b</bold>, Group B; <bold>c</bold>, Group C. The value of the discriminant function for each patient is determined by a weighted combination of image features significantly associated (p &lt; 0.001) with the characteristics of the grading groups. Factors that dealt with a larger minor ellipse axis, indicative of a rounder nucleus, (p &lt; 0.001) and lower peak transition probability, indicating more uniform nuclear staining, (p &lt; 0.001) were associated with higher grading.</p></caption></fig> <fig id="f2-cin-6-0099" position="float"><label>Figure 2</label><caption><p>Distribution of all patients by rank according to the value of their canonical variable. There is a nearly continuous distribution of patients.</p></caption></fig> <fig id="f3-cin-6-0099" position="float"><label>Figure 3</label><caption><p>Kaplan-Meier plot of image feature significantly associated with recurrence of DCIS: Morphometry (ellipticity), p = 0.04. The recurrence rates at 5 years are 25% for patients with Ellipticity less than or equal to the mean, and 17% for those with Ellipticity greater than the mean; the numbers of patients remaining at risk at 5 years are respectively, 23 and 24.</p></caption></fig>
<table-wrap id="t1-cin-6-0099" position="float"><label>Table 1</label><caption><p>Image features discriminating grading groups.</p></caption><table frame="hsides" rules="groups"><thead><tr><th align="left" rowspan="1" colspan="1">Field<xref ref-type="table-fn" rid="tfn1-cin-6-0099">a</xref></th><th align="center" rowspan="1" colspan="1">Image features</th><th align="center" rowspan="1" colspan="1">p-value<xref ref-type="table-fn" rid="tfn2-cin-6-0099">b</xref></th><th align="center" rowspan="1" colspan="1">Correct classification<xref ref-type="table-fn" rid="tfn3-cin-6-0099">c</xref></th></tr></thead><tbody><tr><td align="left" rowspan="1" colspan="1">Field 1</td><td align="center" rowspan="1" colspan="1">Morphologic: minor ellipse axis</td><td align="center" rowspan="1" colspan="1">&lt;0.001</td><td align="center" rowspan="1" colspan="1">65.0%</td></tr><tr><td align="left" rowspan="1" colspan="1"/><td align="center" rowspan="1" colspan="1">Texture: sum entropy</td><td align="center" rowspan="1" colspan="1"/><td align="center" rowspan="1" colspan="1"/></tr><tr><td align="left" rowspan="1" colspan="1">Field 2</td><td align="center" rowspan="1" colspan="1">Morphologic: perimeter</td><td align="center" rowspan="1" colspan="1">&lt;0.001</td><td align="center" rowspan="1" colspan="1">67.1%</td></tr><tr><td align="left" rowspan="1" colspan="1"/><td align="center" rowspan="1" colspan="1">Texture: angular second moment</td><td align="center" rowspan="1" colspan="1"/><td align="center" rowspan="1" colspan="1"/></tr><tr><td align="left" rowspan="1" colspan="1"/><td align="center" rowspan="1" colspan="1">Densitometric: range density</td><td align="center" rowspan="1" colspan="1"/><td align="center" rowspan="1" colspan="1"/></tr><tr><td align="left" rowspan="1" colspan="1">Both fields</td><td align="center" rowspan="1" colspan="1">Morphometric: minor ellipse axis</td><td align="center" rowspan="1" colspan="1">&lt;0.001</td><td align="center" rowspan="1" colspan="1">65.0%</td></tr><tr><td align="left" rowspan="1" colspan="1"/><td align="center" rowspan="1" colspan="1">Texture: peak transition probability</td><td align="center" rowspan="1" colspan="1"/><td align="center" rowspan="1" colspan="1"/></tr></tbody></table></table-wrap> <table-wrap id="t2-cin-6-0099" position="float"><label>Table 2</label><caption><p>Classification of groups with image features.</p></caption><table frame="hsides" rules="groups"><thead><tr><th align="left" rowspan="1" colspan="1">Group</th><th align="center" rowspan="1" colspan="1">Patients in group</th><th colspan="3" align="center" rowspan="1">Number of patients classified into group
<hr/></th><th align="center" rowspan="1" colspan="1">Percent correct</th></tr><tr><th align="left" rowspan="1" colspan="1"/><th align="center" rowspan="1" colspan="1"/><th align="center" rowspan="1" colspan="1">A</th><th align="center" rowspan="1" colspan="1">B</th><th align="center" rowspan="1" colspan="1">C</th><th align="center" rowspan="1" colspan="1"/></tr></thead><tbody><tr><td align="left" rowspan="1" colspan="1">A</td><td align="center" rowspan="1" colspan="1">32</td><td align="center" rowspan="1" colspan="1">25</td><td align="center" rowspan="1" colspan="1">4</td><td align="center" rowspan="1" colspan="1">3</td><td align="center" rowspan="1" colspan="1">78.1</td></tr><tr><td align="left" rowspan="1" colspan="1">B</td><td align="center" rowspan="1" colspan="1">31</td><td align="center" rowspan="1" colspan="1">8</td><td align="center" rowspan="1" colspan="1">15</td><td align="center" rowspan="1" colspan="1">8</td><td align="center" rowspan="1" colspan="1">48.4</td></tr><tr><td align="left" rowspan="1" colspan="1">C</td><td align="center" rowspan="1" colspan="1">17</td><td align="center" rowspan="1" colspan="1">0</td><td align="center" rowspan="1" colspan="1">5</td><td align="center" rowspan="1" colspan="1">12</td><td align="center" rowspan="1" colspan="1">70.6</td></tr><tr><td align="left" rowspan="1" colspan="1">Total</td><td align="center" rowspan="1" colspan="1">80</td><td align="center" rowspan="1" colspan="1">33</td><td align="center" rowspan="1" colspan="1">24</td><td align="center" rowspan="1" colspan="1">23</td><td align="center" rowspan="1" colspan="1">65.0</td></tr></tbody></table></table-wrap> <table-wrap id="t3-cin-6-0099" position="float"><label>Table 3</label><caption><p>Clinical, histologic, and image analysis factors affecting DCIS recurrence by image analysis assessment.</p></caption><table frame="hsides" rules="groups"><thead><tr><th colspan="2" align="center" rowspan="1">Field 1
<hr/></th><th colspan="2" align="center" rowspan="1">Field 2
<hr/></th></tr><tr><th align="center" rowspan="1" colspan="1">Factors<xref ref-type="table-fn" rid="tfn5-cin-6-0099">a</xref></th><th align="center" rowspan="1" colspan="1">p-value</th><th align="center" rowspan="1" colspan="1">Factors<xref ref-type="table-fn" rid="tfn5-cin-6-0099">a</xref></th><th align="center" rowspan="1" colspan="1">p-value</th></tr></thead><tbody><tr><td align="left" rowspan="1" colspan="1">Texture (Histogram mean)</td><td align="center" rowspan="1" colspan="1">0.01</td><td align="left" rowspan="1" colspan="1">Densitometry (Range density)</td><td align="center" rowspan="1" colspan="1">&lt;0.001</td></tr><tr><td align="left" rowspan="1" colspan="1">Initial presentation</td><td align="center" rowspan="1" colspan="1">0.01</td><td align="left" rowspan="1" colspan="1">Measured margin</td><td align="center" rowspan="1" colspan="1">&lt;0.001</td></tr><tr><td align="left" rowspan="1" colspan="1">Parenchymal involvement</td><td align="center" rowspan="1" colspan="1">0.02</td><td align="left" rowspan="1" colspan="1">Densitometry (Sum density)</td><td align="center" rowspan="1" colspan="1">0.02</td></tr><tr><td align="left" rowspan="1" colspan="1">Architecture</td><td align="center" rowspan="1" colspan="1">0.05</td><td align="left" rowspan="1" colspan="1">Van Nuys Prognostic Index</td><td align="center" rowspan="1" colspan="1">0.04</td></tr><tr><td align="left" rowspan="1" colspan="1"><bold> Both fields – Overall</bold></td><td align="center" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="center" rowspan="1" colspan="1"/></tr><tr><td align="left" rowspan="1" colspan="1"> Morphometry (Ellipticity)</td><td align="center" rowspan="1" colspan="1">p = 0.04</td><td align="left" rowspan="1" colspan="1"/><td align="center" rowspan="1" colspan="1"/></tr><tr><td align="left" rowspan="1" colspan="1"> Initial presentation</td><td align="center" rowspan="1" colspan="1">p = 0.03</td><td align="left" rowspan="1" colspan="1"/><td align="center" rowspan="1" colspan="1"/></tr><tr><td align="left" rowspan="1" colspan="1"> Parenchymal involvement</td><td align="center" rowspan="1" colspan="1">p = 0.05</td><td align="left" rowspan="1" colspan="1"/><td align="center" rowspan="1" colspan="1"/></tr></tbody></table></table-wrap>
<fn-group><fn><p><bold>Ethics approval</bold></p><p>This project was reviewed and approved by the Institutional Review Board of Rutgers University, Piscataway, New Jersey, and the Research Ethics Board of Women’s College Hospital, University of Toronto, Toronto, Ontario.</p></fn><fn><p><bold>Competing interests</bold></p><p>WC-B is President and Chief Scientist at Equipoise Imaging LLC. The other authors have no competing interests.</p></fn><fn><p><bold>Authors Contributions</bold></p><p>DA acquired the data, was involved in design of the study, analysis and interpretation of results, drafting and preparing the final manuscript. JC was involved in design of the study, analysis and interpretation of results, drafting and preparing the final manuscript. NM was involved in design of the study, pathological review of specimens, interpretation of results, and preparing the final manuscript. WC-B wrote the software used to extract the image features and reviewed the manuscript. JQ, YY and YF were involved in analysis of data and reviewed the manuscript. HL was involved in acquisition of specimens through surgical management of patients and reviewed the manuscript.</p></fn></fn-group> <table-wrap-foot><fn id="tfn1-cin-6-0099"><label>a</label><p>Field 1 and Field 2 each had 5 different ducts.</p></fn><fn id="tfn2-cin-6-0099"><label>b</label><p>p-values are based on F-statistic of final discriminant models.</p></fn><fn id="tfn3-cin-6-0099"><label>c</label><p>Jackknifed (leave-one-out) classification.</p></fn></table-wrap-foot> <table-wrap-foot><fn id="tfn4-cin-6-0099"><p><bold>Notes:</bold> Classification is based on about 200 nuclei per patient in two fields. Jackknifed (leave-one-out) assessment of classification was used.</p></fn></table-wrap-foot> <table-wrap-foot><fn id="tfn5-cin-6-0099"><label>a</label><p>Factors significantly (p ≤ 0.05) associated with DCIS recurrence, in the order entered into the step-wise forward Cox regression models.</p></fn></table-wrap-foot>
<graphic xlink:href="cin-6-0099f1"/> <graphic xlink:href="cin-6-0099f2"/> <graphic xlink:href="cin-6-0099f3"/>
[{"surname": ["Baak", "Langley", "Hermans", "Baak"], "given-names": ["JPA", "FA", "J", "JPA"], "year": ["1991"], "article-title": ["Classification and prognosis for new cases: Some aspects of univariate and multivariate analysis"], "source": ["Manual of quantitative pathology in cancer diagnosis and prognosis"], "publisher-loc": ["Berlin"], "publisher-name": ["Springer-Verlag"], "fpage": ["189"], "lpage": ["209"]}, {"surname": ["Douglas-Jones", "Gupta", "Attanoos"], "given-names": ["AG", "SK", "RL"], "year": ["1996"], "article-title": ["A critical appraisal of six modern classifications of ductal carcinoma in situ of the breast (DCIS): Correlation with grade of associated invasive carcinoma"], "source": ["Histopathol"], "volume": ["29"], "fpage": ["397"], "lpage": ["409"]}, {"surname": ["Haralick", "Shaanmugam", "Dinstein"], "given-names": ["RM", "K", "I"], "year": ["1973"], "article-title": ["Texture features for image classification"], "conf-name": ["IEEE Trans. Systems Man and Cybernetics"], "volume": ["SMC-3"], "fpage": ["610"], "lpage": ["21"]}, {"surname": ["Hoque", "Lippman", "Boiko"], "given-names": ["A", "SM", "IV"], "year": ["2001"], "article-title": ["Quantitative nuclear morphometry by image analysis for prediction of recurrence of ductal carcinoma in situ of the breast"], "source": ["Cancer Epidemiol. Biomark. Prev"], "volume": ["10"], "fpage": ["249"], "lpage": ["59"]}, {"surname": ["Morrow", "Strom", "Bassett"], "given-names": ["M", "EA", "LW"], "year": ["2002"], "article-title": ["Standard for the management of ductal carcinoma in situ of the breast (DCIS)"], "source": ["CA: Cancer J. Clinicians"], "volume": ["52"], "fpage": ["256"], "lpage": ["76"]}, {"surname": ["Patterson", "Hamilton", "Allen"], "given-names": ["CC", "PW", "DC"], "year": ["1995"], "article-title": ["Statistical analysis of morphometric and clinical data"], "source": ["Quantitative Clinical Pathology"], "publisher-loc": ["Oxford"], "publisher-name": ["Blackwell Sci"], "fpage": ["316"], "lpage": ["32"]}, {"surname": ["Russ"], "given-names": ["JC"], "year": ["1990"], "source": ["Computer-assisted microscopy: The measurement and analysis of images "], "publisher-loc": ["New York"], "publisher-name": ["Plenum Press"]}, {"surname": ["Schulte", "Seigneurin", "Giroud", "Hamilton", "Allen"], "given-names": ["EKW", "D", "F", "PW", "DC"], "year": ["1995"], "article-title": ["DNA densitometry"], "source": ["Quantitative Clinical Pathology"], "publisher-loc": ["Oxford"], "publisher-name": ["Blackwell Sci"], "fpage": ["140"], "lpage": ["69"]}]
Cancer Inform. 2008 Mar 1; 6:99-109