Patient cohort and collection of ultrasound images
The patient cohort in this study comprised consecutive patients with thyroid nodules who underwent regular ultrasound detection and ultrasound-guided FNA between June 2018 and October 2018 at the Ultrasound Department of the First Affiliated Hospital of Dalian Medical University. All ultrasound studies were prospectively generated by a trained sonographer using 5-18 MHz frequencies of the HITACHI HI VISION Ascendus (HITACHI, Japan). Ultrasound dynamic video records were obtained with patients positioned supine and with slightly hyperextended necks. Transverse and longitudinal videos and images of each nodule, using both grayscale and CDFI, were preserved by the ultrasonic machine for the subsequent analysis. Doppler amplification was set to a level at which normal thyroid tissue did not display any noise and just under the level at which random noise appeared. All studies were checked by expert radiologists at the time of the examination to ensure image quality and adequacy. After acquisition of the images, ultrasound-guided FNAs were prospectively performed on these patients by a trained sonographer with 10 years of experience in interventional ultrasound. From the initial patient cohort, all FNA results were evaluated by the expert pathologist for the clinical diagnosis, and patients with unclear pathological diagnosis (Bethesda I, III and IV) or non-papillary carcinoma were excluded. The final nodule population was 299, including 161 malignant nodules (PTCs) and 138 benign nodules. Figure 1 is the detailed diagram of this study group.
Evaluation of ultrasound images
We reviewed electronic medical records to obtain pertinent patient information including sex, age, FNA-proven pathologic diagnosis, and ultrasound dynamic video records. The videos and images of each nodule were evaluated on the ultrasonic machine with Hitachi color monitor (Hitachi, Ltd., Higashi-Ueno, Taito-ku, Tokyo, Japan) with 2-megapixel 1280 × 1024 resolution. The location, size, adjacency, composition, echogenicity, shape, margin and echogenic foci of each nodule were evaluated by two trained sonographers and documented for the following analyses.
To evaluate the nodular microvascular flow, each nodule was visualized with CDFI by the video mode in both transverse and longitudinal section. The video clips were about 5 seconds. In terms of the merging methods recommended by Shin et al(23), Chen et al(24) and our experience, nodular vascularity on CDFI was classified into four types (Fig. 2): type I, absence of nodule vascularity; type II, predominantly perinodular vascularity with continuous (II a) or discontinuous circumferential vascularity at the margin of a nodule (II b); type III, predominantly intranodular vascularity, homogenous linear (III a) or branching (III b) with or without perinodular vascularity; type IV, heterogeneous short-line, strip-like, or microbubble colorful signals, which penetrated the margin of the nodule (perinodular neovascularization, IV a) or dispersely distributed within the nodule (internal neovascularization, IV b). Specifically, the presence of tumor neovascularization-like pattern in CDFI was identified as vascularity of type IV. The vascularity of type I, II and III were excluded from the tumor neovascularization-like pattern. Nineteen freeze-frame shots of these nodules from 19 video clips are provided in Figure 2.
For the evaluation of the echogenic areas, each nodule was evaluated in a binary manner in both transverse and longitudinal sections. Specifically, the echogenic area features (present or absent) were defined as homogeneous mild hyperechoic areas (excluding microcalcifications) located within the relatively hypoechoic nodule. The morphology of the hyperechoic areas was either conglomerate (Fig. 3A) or flocculent (Fig. 3B).
All of the above evaluations were performed independently by two trained sonographers, and the discrepant cases were reassessment by an expert radiologist to reach a consensus. All evaluators were blinded to the pathologic results.
Tissue samples
Tissue samples from 10 patients with malignant thyroid nodules (presence of the echogenic areas on binary), 10 patients with malignant nodules (presence of the tumor neovascularization-like pattern on CDFI) and 10 patients with benign nodules (absence of the tumor neovascularization-like pattern on CDFI) treated at the First Affiliated Hospital of Dalian Medical University (Dalian, Liaoning) were randomly selected. All tissue slides were reviewed by an expert pathologist for verification of the clinical diagnosis. Each tissue sample had at least 3 independent tissue slides for immunohistochemistry (IHC) staining.
IHC staining
IHC was performed on the 5-μm formalin-fixed, paraffin-embedded tissue slides. The DAB chromogenic reagent kit (Cat. ZLI-9019, ZSGB-BIO, Beijing, China) was used for IHC. Slides were dewaxed, rehydrated, antigen retrieved and endogenous peroxidase blocked following the protocol of the kit. For immunolabeling of CD34, mouse monoclonal CD34 (Cat. Kit-0004, MXB Biotechnologies, Fujian, China) was applied. Thereafter, secondary antibody testing, DAB chromogenic reaction, counterstaining, dehydration and transparency were performed according to the protocol of the kit.
IHC scoring system
Immunoreactivity of the vascular endothelial cells was first evaluated for each tissue slide. For all cases, slides that showed distinct immunostaining in vascular endothelial cells were further evaluated. Thereafter, the intensity of CD34 expression for vascular endothelial cells in each field was classified into five grades: extremely high (score 4), high (score 3), moderate (score 2), low (score 1), and negative (score 0). Final scores for each case were calculated as the mean score of all the individual field scores of each slide.
Measuring microvessel density (MVD)
Microvessels of thyroid nodules were highlighted by anti-CD34 immunostaining in formalin-fixed, paraffin-embedded slides. The MVD quantification for each slide was performed according to the detailed method described previously(25, 26). Briefly, the stained slides were examined at low-power magnification (40 and 100 total magnification) to identify the areas of highest neovascularization of the tumor. In each slide, the three most vascular areas were chosen. The microvessel counts in a 200 field (Olympus microscope) in each of these 3 areas was counted. The average counts of the 3 fields in each slide were calculated and thereafter were referred to as the density counts. Any brown-staining endothelial cell clearly separated from adjacent microvessels, tumor cells, or other connective tissue elements was considered to be a single, countable microvessel. Large vessels with thick, muscular walls and large vessels with lumina greater than approximately eight red blood cells were excluded from the count. All measurements were performed independently by two observers.
Scoring system for the individual features
For multivariate analyses, covariates for prediction models included the newly identified features, as well as established risk factors. Specifically, the absence and presence of tumor neovascularization-like pattern were scored as 0 and 1 point, respectively, as well as the echogenic areas. The traditional grayscale features were scored according to the ACR TI-RADS guideline(4). In brief, the composition of cystic, spongiform, mixed and solid was scored as 0, 0, 1, and 2 points, respectively; the echogenicity of anechoic, hyperechoic or isoechoic, hypoechoic and very hypoechoic was scored as 0, 1, 2, and 3 points, respectively; the shapes of wider-than-tall and taller-than-wide were scored as 0 and 3 points, respectively; the margins of smooth, ill-defined, lobulated or irregular and extrathyroidal extension were scored as 0, 0, 2, and 3 points, respectively; and the echogenic foci of none or large comet-tail artifacts, macrocalcifications, peripheral calcifications and punctate echogenic foci were scored as 0, 1, 2, and 3 points, respectively.
Statistical analysis
Data included in this study were analyzed using SPSS version 23.0 (IBM Corporation, Armonk, NY). Chi squared test was used to analyze differences of clinical parameters between two groups of patients. FNA-proven diagnosis was used as the reference standard to determine sensitivity and specificity. Fisher's exact test was used to analyze associations between the sonographic features and malignancy. Student’s t test was used to analyze differences in IHC scores and MVD counts. Logistic regression analysis was performed to evaluate the sonographic features as the predictor for malignancy and the associated odds ratio (OR). Additionally, the relative importance of individual covariates in multivariate logistic regression models was estimated by examining the partial Wald Chi-squared statistic. P < 0.05 was considered statistically significant.