Preoperative Prediction of Ki-67 Labeling Index By Three-dimensional CT Image Parameters for Differential Diagnosis Of Ground-Glass Opacity (GGO)

Abstract

The aim of this study was to predict Ki-67 labeling index (LI) preoperatively by three-dimensional (3D) CT image parameters for pathologic assessment of GGO nodules. Diameter, total volume (TV), the maximum CT number (MAX), average CT number (AVG) and standard deviation of CT number within the whole GGO nodule (STD) were measured by 3D CT workstation. By detection of immunohistochemistry and Image Software Pro Plus 6.0, different Ki-67 LI were measured and statistically analyzed among preinvasive adenocarcinoma (PIA), minimally invasive adenocarcinoma (MIA) and invasive adenocarcinoma (IAC). Receiver operating characteristic (ROC) curve, Spearman correlation analysis and multiple linear regression analysis with cross-validation were performed to further research a quantitative correlation between Ki-67 labeling index and radiological parameters. Diameter, TV, MAX, AVG and STD increased along with PIA, MIA and IAC significantly and consecutively. In the multiple linear regression model by a stepwise way, we obtained an equation: prediction of Ki-67 LI=0.022*STD+0.001* TV+2.137 (R=0.595, R’s square=0.354, p<0.001), which can predict Ki-67 LI as a proliferative marker preoperatively. Diameter, TV, MAX, AVG and STD could discriminate pathologic categories of GGO nodules significantly. Ki-67 LI of early lung adenocarcinoma presenting GGO can be predicted by radiologic parameters based on 3D CT for differential diagnosis.

About the Authors:

Mingzheng Peng

Affiliation: Department of Thoracic Surgery, Shanghai First People’s Hospital Affiliated to The Shanghai Jiao Tong University School Of Medicine, Shanghai, China

Fei Peng

Affiliation: Department of Nephrology, People's Hospital of Hunan Province Affiliated to Hunan Normal University School Of Medicine, Changsha, Hunan Province, China

Chengzhong Zhang

Affiliation: Department of Radiology, Shanghai First People’s Hospital Affiliated to The Shanghai Jiao Tong University School Of Medicine, Shanghai, China

Qingguo Wang

Affiliation: Department of Radiology, Shanghai First People’s Hospital Affiliated to The Shanghai Jiao Tong University School Of Medicine, Shanghai, China

Zhao Li

Affiliation: Department of Thoracic Surgery, Shanghai First People’s Hospital Affiliated to The Shanghai Jiao Tong University School Of Medicine, Shanghai, China

Haiyang Hu

Affiliation: Department of Thoracic Surgery, Shanghai First People’s Hospital Affiliated to The Shanghai Jiao Tong University School Of Medicine, Shanghai, China

Sida Liu

Affiliation: Department of Thoracic Surgery, Shanghai First People’s Hospital Affiliated to The Shanghai Jiao Tong University School Of Medicine, Shanghai, China

Binbin Xu

Affiliation: Department of Thoracic Surgery, Shanghai First People’s Hospital Affiliated to The Shanghai Jiao Tong University School Of Medicine, Shanghai, China

Wenzhuo Zhu

Affiliation: Department of Thoracic Surgery, Shanghai First People’s Hospital Affiliated to The Shanghai Jiao Tong University School Of Medicine, Shanghai, China

Yudong Han

Affiliation: Department of Thoracic Surgery, Shanghai First People’s Hospital Affiliated to The Shanghai Jiao Tong University School Of Medicine, Shanghai, China

Qiang Lin

* E-mail: [email protected]

Affiliation: Department of Thoracic Surgery, Shanghai First People’s Hospital Affiliated to The Shanghai Jiao Tong University School Of Medicine, Shanghai, China

Introduction

Advances in high resolution CT (HRCT) scanning had increased the detection of ground-glass opacity (GGO) with data from many studies suggesting that localized GGO represents as a precursor of lung adenocarcinoma.[1–5] Since a new international multidisciplinary classification of lung adenocarcinoma had been proposed by International Association for the Study of Lung Cancer (IASLC), the American Thoracic Society (ATS), and the European Respiratory Society (ERS) in 2011, pathologic differentiation of GGO has been imperative and attractive for thoracic radiologists and surgeons.[6]

GGO is a finding on HRCT lung images, and has also been described as a hazy increase in lung attenuation without obscuring the underlying bronchial or vascular structures.[3, 7, 8] Furthermore, lesions without solid component in it were classified to pure GGO in comparison to mix GGO with solid component and ground glass attenuation in it as well. As a nonspecific finding that indicates a variety of disorders, it is generally difficult to differentiate with only two-dimensional CT image features, even at follow-up. [3, 9, 10] Advances in understanding of pathologic and radiologic features of GGO have led to changes in diagnostic and therapeutic strategies.[11] Especially, an awareness of the significance of CT attenuation number in assessing GGO has been reported recently.[12–14] Three-dimensional (3D) evaluation has been shown to be more sensitive and precise for quantifying small pulmonary nodules, particularly for asymmetric nodules, than one- or two-dimensional methods.[15, 16]

Besides, proliferation is a key feature for progressing of lung cancer, which is now widely estimated by the immunohistochemical assessment of the nuclear antigen Ki-67. Some authors[17–21] have demonstrated that proliferative activities determined by Ki-67 were correlated with the prognosis of lung cancer patients.

Thus, at this present study, we combined more objective and accuracy parameters obtained from 3D CT image and Ki-67 labeling index (LI) of GGO to analyze their correlation quantitatively and unprecedentedly to predict Ki-67 LI by 3D CT image parameters for preoperative assessment.

Materials and Methods

This study was reviewed and approved by Institutional Review Board of Shanghai First People’s Hospital with Certificate Number of 2014KY115 and written informed consent for patients to participate in this research was obtained before the retrospective study.

Patients

All of our selected cases were adenocarcinoma that recently has been the most frequent pathologic type of lung cancer.[22] All those patients detected out a GGO lesion by HRCT should have been treated by antibiotic for two weeks or so at first. Another CT scan would be performed after at least three months as a follow-up. Only the stable or size-increasing lesions after anti-inflammatory could be selected into this research with excluding those transient ones reckoned as inflammatory in most cases.[23] Only pure and mix focal GGO (part-solid) nodules with diameters less than 3 cm were included. Henschke et al reported that GGO nodules have a higher malignancy rate than full-solid nodules.[24] Feng Li et al also found that there were only 15 malignancy lesions in comparison to 122 benign ones in all 137 cases presenting small full-solid nodules on CT screen. Among full-solid nodules, a polygonal shape or a smooth or somewhat smooth margin was present less frequently in malignant than in benign lesions (polygonal shape: 7% vs. 38%, P = 0.02; smooth or somewhat smooth margin: 0% vs. 63%, P = 0.001), and 98% (46 of 47) of polygonal nodules and 100% (77 of 77) of nodules with a smooth or somewhat smooth margin were benign.[25] Hence, we did not include small full-solid nodules as GGO component attractive more our attention as well as other researchers’.

The exclusion criteria were as follows: (1) cases without pathologic diagnosis or preoperative CT scanning in our hospital, (2) mean diameters of three axes larger than 3 cm, (3) small cell lung cancer, squamous carcinoma and metastatic carcinoma, (4) adenocarcinoma exceed T1N0M0, (5) cases with chemo-radiotherapy or biopsy preoperatively, (6) cases with a history of previous primary lung cancer or extra-pulmonary malignancy. Eventually, 160 patients (54 men and 106 women, mean age 56.59±9.9 years) undergoing HRCT examinations preoperatively were selected. The interval between the latest preoperative HRCT and surgery was 11±4 days with range of 1 to 23 days. Pathologic diagnoses included atypical adenomatous hyperplasia (AAH, n = 26), adenocarcinoma in situ (AIS, n = 11), minimally invasive adenocarcinoma (MIA, n = 106) and invasive adenocarcinoma (IAC, n = 17). Table 1 shows the clinical, radiologic and pathologic characteristics of the selected cases.

[Figure omitted. See PDF.]

Table 1. Clinical, Radiologic and Pathologic Characteristics of All GGO Nodules with different pathologic categories (n = 160).

https://doi.org/10.1371/journal.pone.0129206.t001

Image Acquisition

CT images were obtained from thoracic inlet to lung base with patients at full inspiration. No intravenous contrast material was injected. CT images were acquired by using CT scanners (Somatom Plus 4, Siemens, Erlangen, Germany; LightSpeed Ultra, GE Medical Systems, Milwaukee, Wis; or Mx8000, Philips Medical System, Andover, Mass). Tube voltage ranged from 120 to 140 kV, and tube current ranged from 200 to 400 mA. Axial images of 0.625 mm thickness with 0.5 mm spacing were reconstructed with 512×512 matrix by using bone algorithm axial reconstruction and filtered back projection (FBP) algorithm and the smallest field of view that included both lungs. Two chest radiologists, with 8 and 12 years’ experience of reading CT images of the chest respectively, identified diagnoses of GGO on CT images by consensus.

Computerized Volumetry of GGO Nodules After 3D Measurement

In the present study, we analyzed 5 radiologic parameters of GGO: diameter, total volume (TV), the maximum CT number (MAX), average CT number (AVG) and standard deviation of CT number within the whole GGO nodule (STD) that were measured on a commercially available workstation (Advantage Workstation 4.3; GE Healthcare) with CT lung analysis software (Lung VCAR; GE Healthcare). This software can segment pulmonary nodules with ground-glass attenuation. Diameter, TV, MAX, AVE and STD are computed automatically after the operator placing a marker on the nodule. The CT lung analysis software system automatically identified the GGO nodules in all X-axis, Y-axis and Z-axis directions from the surrounding normal lung tissue. The elimination of normal structures within or around the nodule, such as vessels and bronchiole, was performed using several image-processing techniques[26]. Therefore, the nodule was identified as the lesion area without vessels and bronchiole. Some authors have described the concrete procedures and methods with the exact kind of software.[27–31] The judgment of successful segmentation was based on the observers’ visual assessment on axial CT images as well as sagittal and coronal multiplanar reconstructed images. (Fig 1)

[Figure omitted. See PDF.]

Fig 1. Example of nodule three-dimensional processed and measured on thin-section helical computed tomography (CT) images.

(A) Typical ground-glass opacity (GGO, black arrow) nodule on high-resolution CT. (B) The software automatically processed and measured the nodule that be placed a marker (green square) on it. (C) The magnified measurement list of relative image parameters from picture B including diameter, total volume (TV), the maximum CT value (MAX), average CT value (AVG) and standard deviation of CT value (STD).

https://doi.org/10.1371/journal.pone.0129206.g001

Pathologic Diagnoses

The surgically resected specimens were routinely fixed in 10% formalin, and processed into paraffin blocks for pathologic examination over the entire volumes of GGO nodules. Tissue sections were cut at 4μm thickness, including the largest cut surface of the tumor, and stained with both hematoxylin and eosin (H&E). Pathologic diagnoses were made by two experienced lung pathologists. According to the new IASLC/ATS/ERS classification,[6] all of GGO nodules were diagnosed as AAH, AIS, MIA and IAC. (Fig 2)

[Figure omitted. See PDF.]

Fig 2. Photomicrographs of ground-glass opacity (GGO) (hematoxylin-eosin stain, original magnification, ×400).

(A) A case of atypical adenomatous hyperplasia (AAH) with a well-defined boundary and pulmonary mesenchyme. (B) A case of adenocarcinoma in situ (AIS) with growth restricted to preexisting alveolar walls. No foci of invasion or scarring are seen. (C) A case of minimally invasive adenocarcinoma (MIA) consisting primarily of lepidic growth with a small area of acinar invasion.

https://doi.org/10.1371/journal.pone.0129206.g002

Detection of Ki-67 LI

Immunostaining was performed with the standard streptavidin-perosidase (SP) technique with the antibodies for Ki-67 (monoclonal mouse antibody MIB-1, 1: 100 dilution). Representative formalin fixed, paraffin-embedded tissues were selected for immunohistochemistry. 4 μm thick tumor tissue sections were cut and dried for 20mins at 68°C. Routine deparaffinage by dimethylbenzene and dehydration by gradient ethyl alcohol was performed. 3% hydrogen peroxide was used to block endogenous peroxidase for 10mins at 37°C and phosphate buffer solution (PBS) was used to douche lasting 5mins for 3 times. Then citrate buffer with pH of 6.0 was used to repair antigen for 3mins. After that, sections were incubated with primary monoclonal mouse antibodies overnight at 4°C. The next step is put them into horseradish peroxidase complex (HRP) and diaminobenzidine (DAB) for 30mins and 10mins respectively at 37°C. Finally, routine hematoxylin redyeing, dehydration, dimethylbenzene transparentizing and gum mounting had been performed step by step. Quantitative analyses of positive expression of Ki-67 indicated by tan staining particles located in nucleus were performed using Image Software Pro Plus 6.0. Not only verified its accuracy by some authors,[32] this software has been applied to numerous aspects in biomedicine recently.[33–36] (Fig 3)

[Figure omitted. See PDF.]

Fig 3. Representative immunohistochemical detection and measurement of Ki-67 labeling index (LI) (magnifyication ×400).

(A) A case with positive expression of Ki-67 presenting tan particles located in nucleus of proliferating cells. (B) Quantitative measurement of Ki-67 LI by Image Software Pro Plus 6.0 (red area means detection of the software).

https://doi.org/10.1371/journal.pone.0129206.g003

Statistic Methods

Numerical data were recorded as mean ± standard deviation (SD). Levene test was performed for evaluating homogeneity of variance of variables. Consequently, variances of all variables other than that of Ki-67 LI were heterogeneous. (Table 2) Kruskal-Wallis test was taken for analyzing significant difference of variables with heterogeneous variance and analysis of variance (ANOVA) for Ki-67 LI with homogeneous variance among groups. Furthermore, we made Least Significant Difference (LSD) test for pairwise comparison of Ki-67 LI and Tamhane’s T2 for pairwise comparison of diameter, TV, MAX, AVG, STD and prediction of Ki-67 LI, whose variances were heterogeneous. In receiver operating curve (ROC) analysis, optimal cut-off points were defined as those on the curves closest to the upper left-hand corner[37]. Spearman correlation and multiple linear regression analysis with 10-fold cross-validation were performed to research the quantitative relation between Ki-67 LI and radiological parameters. All analyses were conducted with a statistical software package (SPSS, version 20.0). P-values of less than 0.05 were considered to indicate statistical significance.

[Figure omitted. See PDF.]

Table 2. Levene test of variables for homogeneity of variance.

https://doi.org/10.1371/journal.pone.0129206.t002

Results

Radiological Features and Histological Diagnoses

Parameters of GGO were automatically calculated and displayed on the multi-display CT (MDCT) after 3D measurement by the CT lung analysis software. (Fig 4) The positive expression of Ki-67 presented tan particles located in nucleus of over proliferative cell. All of relevant data were recorded in Table 1.

[Figure omitted. See PDF.]

Fig 4. Graphs showing the mean diameter, total volume (TV), the maximum CT value (MAX), average CT value (AVG), standard deviation of CT value within the whole nodule (STD) and Ki-67 labeling index (LI) for preinvasive adenocarcinoma (PIA), minimally invasive adenocarcinoma (MIA) and invasive adenocarcinoma (IAC).

Ki-67 LI was analyzed by ANOVA and LSD test due to homogeneity of variance and other variables were analyzed by Kurskal-Wallis test and Tamhane’s T2 for heterogeneity of variance. While all variables were significantly different among PIA, MIA and IAC.

https://doi.org/10.1371/journal.pone.0129206.g004

ROC Analysis for Pathologically Discrimination

By ROC analysis, area under the curve (AUC) of diameter, TV, MAX, AVG, STD and Ki-67 LI is 0.801, 0.822, 0.890, 0.857, 0.901 and 0.907 respectively to differentiate PIA from MIA and 0.812, 0.793, 0.749, 0.731, 0.684 and 0.901 respectively to rule out MIA from IAC. When it comes to the maximum true-positive ratio (sensitivity) and the minimum false-positive ratio (1—specificity) simultaneously for parameters, we obtained the thresholds of all variables. The optimal cut-off point of diameter, TV, MAX, AVG, STD and Ki-67 LI is 10.55 mm, 217 mm3, -126.5 HU, -615.5 HU, 135.5 and 4.38% respectively for the comparison between PIA and MIA. Corresponding to discrimination between MIA and IAC, the thresholds were 21.8 mm, 1708.5 mm3, 189 HU, -464 HU, 169.4 and 9.88% for diameter, TV, MAX, AVG, STD and Ki-67 LI respectively. (Fig 5)

[Figure omitted. See PDF.]

Fig 5. Receiver operating characteristic (ROC) curves of diameter, total volume (TV), the maximum CT value (MAX), average CT value (AVG), standard deviation of CT value within the whole nodule (STD) and Ki-67 labeling index (LI) for differentiation among PIA, MIA and IAC.

According to area under the curve (AUC) of all variables, Ki-67 LI had a better differential capability than other parameters. It is easier to differential diagnose between PIA and MIA than between MIA and IAC.

https://doi.org/10.1371/journal.pone.0129206.g005

Correlation between 3D Image Parameters and Ki-67 LI

Based on Spearman correlative analysis to Ki-67 LI, the correlative coefficient was 0.575, 0.559, 0.605, 0.585 and 0.639 for diameter (p<0.001), TV (p<0.001), MAX (p<0.001), AVG (p<0.001) and STD (p<0.001) respectively. (Fig 6)

[Figure omitted. See PDF.]

Fig 6. Graphs showing the correlation between Ki-67 labeling index (LI) and diameter, total volume (TV), the maximum CT value (MAX), average CT value (AVG), standard deviation of CT value within the whole nodule (STD) respectively, and correlation between STD and percentage of solid component in ground-glass opacity (GGO) nodule.

(r means correlative coefficient) Diameter, TV, MAX, AVG and STD were positive correlation with Ki-67. Graph (S) shows that in the heterogeneous model of GGO, STD increase along with solid component growing in it instead of decrease after solid-part more than 50% of nodule.

https://doi.org/10.1371/journal.pone.0129206.g006

Furthermore, multiple linear regression analysis was performed to quantitatively explore the relation between 3D parameters and Ki-67 LI. We also made a 10-fold cross-validation. After dividing 160 cases into 10 equal groups randomly, we defined 1 group including 10 cases as testing set and the other 9 groups as training set every time by turns. Next, we established 10 multi-linear regression models by training sets by a stepwise way and then put respective testing sets into its model to get 10 predictive value of Ki-67 every time. In every model, we calculated mean absolute error (MAE), mean relative error (MRE), and root mean square error (RMSE) for evaluating accuracy of equations obtained above. (Table 3)

[Figure omitted. See PDF.]

Table 3. 10-fold cross-validation of multi-linear regression model for Ki-67 by CT image parameters of 160 cases.

https://doi.org/10.1371/journal.pone.0129206.t003

(1)(2)(3)

In above formulas, n means number of cases in testing set, means predictive value of Ki-67 LI, and means actual measured Ki-67 LI.

Among 10 regression equations, No. 5 owns higher predictive accuracy than others in spite that it was not the best fitting one. As a result of cross-validation, we established a more accurate regression model for predicting Ki-67 LI by CT image parameters based on three-dimensional process, stating as follows:(4)

The prediction of Ki-67 LI was 4.27±0.76 (range, 2.67~6.05%), 6.50±1.69 (range,. 3.00~11.74%) and 8.81±3.73 (range, 3.75~16.65%) for PIA, MIA and IAC respectively (ANOVA P<0.001, LSD PIA vs. MIA P<0.001, MIA vs. IAC P = 0.014, PIA vs. IAC P<0.001). (Fig 7) Additionally, based on ROC analysis, we compared the AUC of the predictions with actual Ki-67 LI. Although discriminating not as well as actual Ki-67 LI does, the prediction of Ki-67 has a relatively higher AUC in ROC analysis to differentiate pathology of GGO nodules both between PIA and MIA (AUC = 0.893, Threshold = 5.71%) and between MIA and IAC (AUC = 0.841, Threshold = 9.54%) than diameter, TV, MAX, AVG and STD separately. (Fig 8)

[Figure omitted. See PDF.]

Fig 7. Boxplots comparing the mean predictive Ki-67 labeling index (LI) with mean actual measured Ki-67 LI of preinvasive adenocarcinoma (PIA), minimally invasive adenocarcinoma (MIA) and invasive adenocarcinoma (IAC).

Predictive Ki-67 LI and actual measured Ki-67 LI were all significantly different among PIA, MIA and IAC.

https://doi.org/10.1371/journal.pone.0129206.g007

[Figure omitted. See PDF.]

Fig 8. Receiver operating characteristic (ROC) curves of predictive Ki-67 labeling index (LI) and actual Ki-67 LI for comparison.

Area under the curve (AUC) of predictive Ki-67 LI was lower than actual measured one, but approximated.

https://doi.org/10.1371/journal.pone.0129206.g008

Discussion

Recent advances in imaging technology promoted the more precise and integrated 3D measurement of volume and CT number of pulmonary nodules than traditional one-dimensional or two-dimensional measurement methods.[15, 38–40] As for quantitative and more accurate analysis, given GGO lesion was not homogeneous and no clear dividing line exists between pure GGO and mix GGO, we just regarded them as a continuous progression paralleling to pathological progression from AAH through AIS to MIA. From this point of view, pure GGO is in a very beginning stage and mix GGO may be a middle one depending on percent of solid. Pure GGO lesions most likely own a low MAX, AVG and STD. When the first percent solid appear in pure GGO, it becomes “juvenile” mix GGO and own a higher MAX and little higher STD but a nearly AVG. Thereafter, solid component increased and moderate mix GGO comes up which own a larger diameter, MAX, AVG and STD.

Furthermore, although extensive studies have been conducted on the CT features of GGO nodules, little was known about the relationship between 3D parameters of GGO nodules and quantitative pathologic index such as the proliferative marker Ki-67. It is recorded that Ki-67 protein was originally defined by the prototype monoclonal antibody Ki-67[41], which was generated by immunizing mice with nuclei of the Hodgkin lymphoma cell line L428. The fact that the Ki-67 protein is present during all active phases of the cell cycle (G1, S, G2, and mitosis), but is absent from resting cells (G0), makes it an excellent marker for determining the so-called growth fraction of a given cell population[42]. High proliferative rate is a hallmark of cancer. Thus proliferation is reported to predict poor survival in multiple myeloma,[43] prostate cancer,[42, 44] and breast cancer.[45, 46]

However, on account of parameters such as size and CT number of GGO manifesting the quantity of cancerous cells, to a certain extent, as the same with Ki-67 LI in lung cancer, we considered that the correlation and regression model between them are biological and pathophysiological more than just statistical and mathematical. Most importantly, Ki-67 LIs in lesions presenting GGO of different pathological categories are significant difference. If we could predict it before surgery, radiologic diagnosis of GGO will be more accurate than only through radiological features. Therefore, preoperative prediction of Ki-67 LI could contribute to differential diagnosis of GGO nodules from not only radiology but pathology as well.

The pathological diagnosis is significant for selection of surgical procedures and prognosis. As a precursor of lung adenocarcinoma, malignant lesions undergo a progression along with AAH, AIS, MIA and IAC. As no invasion of lung parenchyma and lymph nodes, limited resections such as wedge resection and segmentectomy are appropriate for PIA including AAH and AIS. As for MIA and IAC that some cases have been detected invasion of local lymph nodes, standard lobectomy would be necessary.[47–50] AIS and MIA should have 100% and near 100% five-year disease-free survival (DFS) respectively when completely resected.[51–54]

In this study, diameter, TV, MAX, AVG, STD and Ki-67 LI were significantly and consecutively increasing along with the progression from PIA, MIA to IAC. Diameter and TV characterize GGO nodules from the perspective of whole basic framework of proliferation in comparison to MAX and AVG describing its substantiality and internality. Ultimately, STD manifests the degree of heterogeneity within the whole nodule. However, all of them are preoperative variables. Pathological diagnosis couldn’t be obtained before surgery or biopsy. Only if there is a quantitative pathologic index that could be measured or predicted preoperatively, can we accurately judge pathological features of GGO nodules.

By ROC analysis, the threshold of diameter and AVG for differentiating PIA from MIA and MIA from IAC was 10.55 mm, 21.8 mm, -615 HU, -464 HU respectively. Lee HY[55] explored the threshold of diameter and AVG of pure GGO for discriminating invasive adenocarcinoma with 15 mm and -472 HU respectively. By ROC analysis as well, Mamoru Takahashi[56] got the threshold of diameter for differentiating noninvasive lesions from invasive with 10 mm.

By Spearman correlative analysis, all of our radiologic image parameters have a significant positive correlation with Ki-67 LI. The increasing of Ki-67 LI along with the progression means rising proliferation of cancer cells. It made the parenchyma more thick and substantial, which present a higher CT number on HRCT. As for the invasion of MIA and IAC presenting prompt proliferation of cancerous cells at scattered point in lung tissue, thicknesses within lesion will be out of consistent step resulting in higher STD. Therefore, it is reasonable to presume that Ki-67 LI from pathology corresponds to parameters with regard to size and CT number from radiology. Furthermore, aim to predict Ki-67 LI preoperatively, we established a multivariable linear regression equation to assess proliferative degree of GGO nodules approximately and compared the prediction Ki-67 LI with actual one. AUC and threshold of predictive Ki-67 LI were not only close to those of actual one, but also better than those of other parameters in this study separately. We obtained the predictable equation as follow: Prediction (Ki-67 LI) = 0.022*STD+0.001* TV+2.137.

Some papers had shown that DNA synthesis was inhibited by oligodeoxynucleotides complementary to Ki-67 mRNA. Moreover, the microinjection of antibodies directed against the murine Ki-67 or the human Ki-67 resulted in a decreased rate of cell division. The phosphorylation and dephosphorylation of the Ki-67 protein controlled by the key regulatory complex cyclin B/cdc 2 is parallel to the transit of cells through mitosis.[57–60] As for volume of lesion in the equation, benign lung parenchyma around fringe of lesion would proliferate when Ki-67 LI increased. The pixels there then would get higher density distinguishing parenchyma around nodule than more peripheral structure and making them become part of GGO lesion. It is obvious and understandable that volume of GGO would increase along with Ki-67 LI. STD means the degree of heterogeneity of nodule. (Fig 1) When a lesion appeared, it presented pure GGO judged by subjective observation. In AAH, cells all over the nodule just proliferate with a low expression of Ki-67 so that difference of proliferating capability among cells was small, resulting in low STD on the CT screen. In AIS, proliferation was constrained in alveolar epithelium where have a higher CT number. Interstitial tissue was clean and has not been invaded by carcinoma cells. So it presented lower CT number. Then STD increased due to increasing difference of CT numbers among scattered proliferating areas. As lung adenocarcinoma progressing, in MIA or IAC, more cells express Ki-67 and proliferate irregularly. Some fast growing areas in nodule produced solid and other slow ones remained ground glass attenuation. Invasive growing made some interstitium filled with proliferating carcinoma cells. Though, in the solid part, pixels would not manifest the same density in CT screen. In short, higher Ki-67 LI, more irregular proliferation and larger STD would accordingly appear along with progression of lung adenocarcinoma. Just as result of Spearman correlative analysis, TV and STD had a significant positive correlation with Ki-67LI.

In conclusion, image parameters, such as diameter, TV, MAX, AVG and STD based on 3D CT, could discriminate the different pathology in the progression of lung adenocarcinoma and contribute to prognosis relatively as well. Additionally, we established a regression formula to predict Ki-67 labeling index of GGO lesions based on those parameters so that we can estimate the GGO pathologic status more accuracy preoperatively.

Acknowledgments

It should be appreciated that Ms. Honghui Tang and Mr. Yifei Wang gave us statistical instruction.

Author Contributions

Conceived and designed the experiments: MP QL. Performed the experiments: MP CZ QW ZL BX. Analyzed the data: MP SL HH YH WZ. Contributed reagents/materials/analysis tools: QW QL. Wrote the paper: MP QL.

Citation: Peng M, Peng F, Zhang C, Wang Q, Li Z, Hu H, et al. (2015) Preoperative Prediction of Ki-67 Labeling Index By Three-dimensional CT Image Parameters for Differential Diagnosis Of Ground-Glass Opacity (GGO). PLoS ONE 10(6): e0129206. https://doi.org/10.1371/journal.pone.0129206

References

1. Jang HJ, Lee KS, Kwon OJ, Rhee CH, Shim YM, and Han J. Bronchioloalveolar carcinoma: focal area of ground-glass attenuation at thin-section CT as an early sign[J]. Radiology, 1996, 199(2): 485–488. pmid:8668800

2. Aoki T, Nakata H, Watanabe H, Nakamura K, Kasai T, Hashimoto H, et al. Evolution of peripheral lung adenocarcinomas: CT findings correlated with histology and tumor doubling time[J]. AJR Am J Roentgenol, 2000, 174(3): 763–768. pmid:10701622

3. Nakajima R, Yokose T, Kakinuma R, Nagai K, Nishiwaki Y, and Ochiai A. Localized pure ground-glass opacity on high-resolution CT: histologic characteristics[J]. J Comput Assist Tomogr, 2002, 26(3): 323–329. pmid:12016356

4. Kakinuma R, Ohmatsu H, Kaneko M, Kusumoto M, Yoshida J, Nagai K, et al. Progression of focal pure ground-glass opacity detected by low-dose helical computed tomography screening for lung cancer[J]. J Comput Assist Tomogr, 2004, 28(1): 17–23. pmid:14716227

5. Nakata M, Saeki H, Takata I, Segawa Y, Mogami H, Mandai K, et al. Focal ground-glass opacity detected by low-dose helical CT[J]. Chest, 2002, 121(5): 1464–1467. pmid:12006429

6. Travis WD, Brambilla E, Noguchi M, Nicholson AG, Geisinger KR, Yatabe Y, et al. International association for the study of lung cancer/american thoracic society/european respiratory society international multidisciplinary classification of lung adenocarcinoma[J]. J Thorac Oncol, 2011, 6(2): 244–285. pmid:21252716

7. Austin JH, Muller NL, Friedman PJ, Hansell DM, Naidich DP, Remy-Jardin M, et al. Glossary of terms for CT of the lungs: recommendations of the Nomenclature Committee of the Fleischner Society[J]. Radiology, 1996, 200(2): 327–331. pmid:8685321

8. Hansell DM, Bankier AA, MacMahon H, McLoud TC, Muller NL, and Remy J. Fleischner Society: glossary of terms for thoracic imaging[J]. Radiology, 2008, 246(3): 697–722. pmid:18195376

9. Collins J and Stern EJ. Ground-glass opacity at CT: the ABCs[J]. AJR Am J Roentgenol, 1997, 169(2): 355–367. pmid:9242736

10. Park CM, Goo JM, Lee HJ, Lee CH, Kim HC, Chung DH, et al. CT findings of atypical adenomatous hyperplasia in the lung[J]. Korean J Radiol, 2006, 7(2): 80–86. pmid:16799268

11. Nakata M, Sawada S, Saeki H, Takashima S, Mogami H, Teramoto N, et al. Prospective study of thoracoscopic limited resection for ground-glass opacity selected by computed tomography[J]. Ann Thorac Surg, 2003, 75(5): 1601–1605; discussion 1605–1606. pmid:12735586

12. Ikeda K, Awai K, Mori T, Kawanaka K, Yamashita Y, and Nomori H. Differential diagnosis of ground-glass opacity nodules: CT number analysis by three-dimensional computerized quantification[J]. Chest, 2007, 132(3): 984–990. pmid:17573486

13. Nomori H, Ohtsuka T, Naruke T, and Suemasu K. Differentiating between atypical adenomatous hyperplasia and bronchioloalveolar carcinoma using the computed tomography number histogram[J]. Ann Thorac Surg, 2003, 76(3): 867–871. pmid:12963218

14. Nomori H, Ohtsuka T, Naruke T, and Suemasu K. Histogram analysis of computed tomography numbers of clinical T1 N0 M0 lung adenocarcinoma, with special reference to lymph node metastasis and tumor invasiveness[J]. J Thorac Cardiovasc Surg, 2003, 126(5): 1584–1589. pmid:14666037

15. Kostis WJ, Yankelevitz DF, Reeves AP, Fluture SC, and Henschke CI. Small pulmonary nodules: reproducibility of three-dimensional volumetric measurement and estimation of time to follow-up CT[J]. Radiology, 2004, 231(2): 446–452. pmid:15128989

16. Revel MP, Lefort C, Bissery A, Bienvenu M, Aycard L, Chatellier G, et al. Pulmonary nodules: preliminary experience with three-dimensional evaluation[J]. Radiology, 2004, 231(2): 459–466. pmid:15128991

17. Ji Y, Zheng M, Ye S, Chen J, and Chen Y. PTEN and Ki67 expression is associated with clinicopathologic features of non-small cell lung cancer[J]. J Biomed Res, 2014, 28(6): 462–467. pmid:25469115

18. Berghoff AS, Ilhan-Mutlu A, Wohrer A, Hackl M, Widhalm G, Hainfellner JA, et al. Prognostic significance of Ki67 proliferation index, HIF1 alpha index and microvascular density in patients with non-small cell lung cancer brain metastases[J]. Strahlenther Onkol, 2014, 190(7): 676–685. pmid:24577133

19. Lei B, Liu S, Qi W, Zhao Y, Li Y, Lin N, et al. PBK/TOPK expression in non-small-cell lung cancer: its correlation and prognostic significance with Ki67 and p53 expression[J]. Histopathology, 2013, 63(5): 696–703. pmid:24025073

20. Ciancio N, Galasso MG, Campisi R, Bivona L, Migliore M, and Di Maria GU. Prognostic value of p53 and Ki67 expression in fiberoptic bronchial biopsies of patients with non small cell lung cancer[J]. Multidiscip Respir Med, 2012, 7(1): 29. pmid:22978804

21. Scagliotti GV, Micela M, Gubetta L, Leonardo E, Cappia S, Borasio P, et al. Prognostic significance of Ki67 labelling in resected non small cell lung cancer[J]. Eur J Cancer, 1993, 29A(3): 363–365. pmid:8398336

22. Siegel R, Naishadham D, and Jemal A. Cancer statistics, 2012[J]. CA Cancer J Clin, 2012, 62(1): 10–29. pmid:22237781

23. Lee HY and Lee KS. Ground-glass opacity nodules: histopathology, imaging evaluation, and clinical implications[J]. J Thorac Imaging, 2011, 26(2): 106–118. pmid:21508733

24. Henschke CI, Yankelevitz DF, Mirtcheva R, McGuinness G, McCauley D, Miettinen OS, et al. CT screening for lung cancer: frequency and significance of part-solid and nonsolid nodules[J]. AJR Am J Roentgenol, 2002, 178(5): 1053–1057. pmid:11959700

25. Li F, Sone S, Abe H, Macmahon H, and Doi K. Malignant versus benign nodules at CT screening for lung cancer: comparison of thin-section CT findings[J]. Radiology, 2004, 233(3): 793–798. pmid:15498895

26. Maksimovic R, Stankovic S, and Milovanovic D. Computed tomography image analyzer: 3D reconstruction and segmentation applying active contour models—'snakes'[J]. Int J Med Inform, 2000, 58–59: 29–37.

27. Oda S, Awai K, Murao K, Ozawa A, Yanaga Y, Kawanaka K, et al. Computer-aided volumetry of pulmonary nodules exhibiting ground-glass opacity at MDCT[J]. AJR Am J Roentgenol, 2010, 194(2): 398–406. pmid:20093602

28. Awai K, Murao K, Ozawa A, Nakayama Y, Nakaura T, Liu D, et al. Pulmonary nodules: estimation of malignancy at thin-section helical CT—effect of computer-aided diagnosis on performance of radiologists[J]. Radiology, 2006, 239(1): 276–284. pmid:16467210

29. Linning E and Daqing M. Volumetric measurement pulmonary ground-glass opacity nodules with multi-detector CT: effect of various tube current on measurement accuracy—a chest CT phantom study[J]. Acad Radiol, 2009, 16(8): 934–939. pmid:19409818

30. Matsumoto T, Kanzaki M, Amiki M, Shimizu T, Maeda H, Sakamoto K, et al. Comparison of three software programs for three-dimensional graphic imaging as contrasted with operative findings[J]. Eur J Cardiothorac Surg, 2012, 41(5): 1098–1103. pmid:22219443

31. Chen B, Barnhart H, Richard S, Robins M, Colsher J, and Samei E. Volumetric quantification of lung nodules in CT with iterative reconstruction (ASiR and MBIR)[J]. Med Phys, 2013, 40(11): 111902. pmid:24320435

32. Wang CJ, Zhou ZG, Holmqvist A, Zhang H, Li Y, Adell G, et al. Survivin expression quantified by Image Pro-Plus compared with visual assessment[J]. Appl Immunohistochem Mol Morphol, 2009, 17(6): 530–535. pmid:19407655

33. Blatt RJ, Clark AN, Courtney J, Tully C, and Tucker AL. Automated quantitative analysis of angiogenesis in the rat aorta model using Image-Pro Plus 4.1[J]. Comput Methods Programs Biomed, 2004, 75(1): 75–79. pmid:15158050

34. Francisco JS, Moraes HP, and Dias EP. Evaluation of the Image-Pro Plus 4.5 software for automatic counting of labeled nuclei by PCNA immunohistochemistry[J]. Braz Oral Res, 2004, 18(2): 100–104. pmid:15311310

35. Rawlinson A, Elcock C, Cheung A, Al-Buhairi A, Khanna S, Walsh TF, et al. An in-vitro and in-vivo methodology study of alveolar bone measurement using extra-oral radiographic alignment apparatus, Image Pro-Plus software and a subtraction programme[J]. J Dent, 2005, 33(9): 781–788. pmid:15922503

36. Neubauer AS, Rothschuh A, Ulbig MW, and Blum M. Digital fundus image grading with the non-mydriatic Visucam(PRO NM) versus the FF450(plus) camera in diabetic retinopathy[J]. Acta Ophthalmol, 2008, 86(2): 177–182. pmid:17944975

37. Akobeng AK. Understanding diagnostic tests 3: Receiver operating characteristic curves[J]. Acta Paediatr, 2007, 96(5): 644–647. pmid:17376185

38. Jennings SG, Winer-Muram HT, Tarver RD, and Farber MO. Lung tumor growth: assessment with CT—comparison of diameter and cross-sectional area with volume measurements[J]. Radiology, 2004, 231(3): 866–871. pmid:15163822

39. Yankelevitz DF, Reeves AP, Kostis WJ, Zhao B, and Henschke CI. Small pulmonary nodules: volumetrically determined growth rates based on CT evaluation[J]. Radiology, 2000, 217(1): 251–256. pmid:11012453

40. Marten K, Funke M, and Engelke C. Flat panel detector-based volumetric CT: prototype evaluation with volumetry of small artificial nodules in a pulmonary phantom[J]. J Thorac Imaging, 2004, 19(3): 156–163. pmid:15273611

41. Gerdes J, Schwab U, Lemke H, and Stein H. Production of a mouse monoclonal antibody reactive with a human nuclear antigen associated with cell proliferation[J]. Int J Cancer, 1983, 31(1): 13–20. pmid:6339421

42. Scholzen T and Gerdes J. The Ki-67 protein: from the known and the unknown[J]. J Cell Physiol, 2000, 182(3): 311–322. pmid:10653597

43. Alexandrakis MG, Passam FH, Kyriakou DS, Dambaki K, Niniraki M, and Stathopoulos E. Ki-67 proliferation index: correlation with prognostic parameters and outcome in multiple myeloma[J]. Am J Clin Oncol, 2004, 27(1): 8–13. pmid:14758126

44. Li R, Heydon K, Hammond ME, Grignon DJ, Roach M 3rd, Wolkov HB, et al. Ki-67 staining index predicts distant metastasis and survival in locally advanced prostate cancer treated with radiotherapy: an analysis of patients in radiation therapy oncology group protocol 86–10[J]. Clin Cancer Res, 2004, 10(12 Pt 1): 4118–4124. pmid:15217948

45. Urruticoechea A, Smith IE, and Dowsett M. Proliferation marker Ki-67 in early breast cancer[J]. J Clin Oncol, 2005, 23(28): 7212–7220. pmid:16192605

46. Stuart-Harris R, Caldas C, Pinder SE, and Pharoah P. Proliferation markers and survival in early breast cancer: a systematic review and meta-analysis of 85 studies in 32,825 patients[J]. Breast, 2008, 17(4): 323–334. pmid:18455396

47. Kodama K, Higashiyama M, Takami K, Oda K, Okami J, Maeda J, et al. Treatment strategy for patients with small peripheral lung lesion(s): intermediate-term results of prospective study[J]. Eur J Cardiothorac Surg, 2008, 34(5): 1068–1074. pmid:18760618

48. Yoshida J, Nagai K, Yokose T, Nishimura M, Kakinuma R, Ohmatsu H, et al. Limited resection trial for pulmonary ground-glass opacity nodules: fifty-case experience[J]. J Thorac Cardiovasc Surg, 2005, 129(5): 991–996. pmid:15867771

49. Nakao M, Yoshida J, Goto K, Ishii G, Kawase A, Aokage K, et al. Long-term outcomes of 50 cases of limited-resection trial for pulmonary ground-glass opacity nodules[J]. J Thorac Oncol, 2012, 7(10): 1563–1566. pmid:22878750

50. Tsutani Y, Miyata Y, Nakayama H, Okumura S, Adachi S, Yoshimura M, et al. Appropriate sublobar resection choice for ground glass opacity-dominant clinical stage IA lung adenocarcinoma: wedge resection or segmentectomy[J]. Chest, 2014, 145(1): 66–71. pmid:24551879

51. Borczuk AC, Qian F, Kazeros A, Eleazar J, Assaad A, Sonett JR, et al. Invasive size is an independent predictor of survival in pulmonary adenocarcinoma[J]. Am J Surg Pathol, 2009, 33(3): 462–469. pmid:19092635

52. Sakurai H, Maeshima A, Watanabe S, Suzuki K, Tsuchiya R, Maeshima AM, et al. Grade of stromal invasion in small adenocarcinoma of the lung: histopathological minimal invasion and prognosis[J]. Am J Surg Pathol, 2004, 28(2): 198–206. pmid:15043309

53. Suzuki K, Asamura H, Kusumoto M, Kondo H, and Tsuchiya R. "Early" peripheral lung cancer: prognostic significance of ground glass opacity on thin-section computed tomographic scan[J]. Ann Thorac Surg, 2002, 74(5): 1635–1639. pmid:12440622

54. Son JY, Lee HY, Lee KS, Kim JH, Han J, Jeong JY, et al. Quantitative CT analysis of pulmonary ground-glass opacity nodules for the distinction of invasive adenocarcinoma from pre-invasive or minimally invasive adenocarcinoma[J]. PLOS One, 2014, 9(8): e104066. pmid:25102064

55. Lee HY, Choi YL, Lee KS, Han J, Zo JI, Shim YM, et al. Pure ground-glass opacity neoplastic lung nodules: histopathology, imaging, and management[J]. AJR Am J Roentgenol, 2014, 202(3): W224–233. pmid:24555618

56. Takahashi M, Shigematsu Y, Ohta M, Tokumasu H, Matsukura T, and Hirai T. Tumor invasiveness as defined by the newly proposed IASLC/ATS/ERS classification has prognostic significance for pathologic stage IA lung adenocarcinoma and can be predicted by radiologic parameters[J]. J Thorac Cardiovasc Surg, 2014, 147(1): 54–59. pmid:24131786

57. Absmanner B, Schmeiser V, Kampf M, and Lehle L. Biochemical characterization, membrane association and identification of amino acids essential for the function of Alg11 from Saccharomyces cerevisiae, an alpha1,2-mannosyltransferase catalysing two sequential glycosylation steps in the formation of the lipid-linked core oligosaccharide[J]. Biochem J, 2010, 426(2): 205–217. pmid:19929855

58. Schluter C, Duchrow M, Wohlenberg C, Becker MH, Key G, Flad HD, et al. The cell proliferation-associated antigen of antibody Ki-67: a very large, ubiquitous nuclear protein with numerous repeated elements, representing a new kind of cell cycle-maintaining proteins[J]. J Cell Biol, 1993, 123(3): 513–522. pmid:8227122

59. Tian H, Qian GW, Li W, Chen FF, Di JH, Zhang BF, et al. A critical role of Sp1 transcription factor in regulating the human Ki-67 gene expression[J]. Tumour Biol, 2011, 32(2): 273–283. pmid:20963645

60. Kausch I, Jiang H, Ewerdwalbesloh N, Doehn C, Kruger S, Sczakiel G, et al. Inhibition of Ki-67 in a renal cell carcinoma severe combined immunodeficiency disease mouse model is associated with induction of apoptosis and tumour growth inhibition[J]. BJU Int, 2005, 95(3): 416–420. pmid:15679806