NFATc1 Expression as a Prognosticator in Urothelial Carcinoma of the Upper Urinary Tract

We recently found that NFATc1, a member of the NFAT family and a key regulator of the immune response, could induce bladder carcinogenesis and cancer progression. In this study, we immunohistochemically stained for NFATc1 in upper urinary tract urothelial carcinoma (UUTUC) specimens and paired nonneoplastic urothelial tissues. NFATc1 was positive in 51 [52%; 40 (40%) weak (1+), 9 (9%) moderate (2+), and 2 (2%) strong (3+)] of 99 UUTUCs, which was significantly higher than in benign urothelium [30 (36%) of 83; 28 (34%) weak and 2 (2%) moderate] (0 vs 1+/2+/3+, P = .038; 0/1+ vs 2+/3+, P = .023). There were no significant associations between NFATc1 expression pattern and tumor grade or pT stage. However, the positive rates of NFATc1 expression tended to be higher in renal pelvic tumors (60%) than in ureteral tumors (42%; P = .080) as well as in pN+ tumors (75%) than in pN0 tumors (49%; P = .089). Kaplan-Meier and log-rank tests revealed that moderate (2+) to strong (3+) NFATc1 expression correlated with lower progression-free survival (P = .032) and cancer-specific survival (P = .005) rates in the 99 cases. Patients with high (2+/3+) NFATc1 muscle-invasive tumor (n = 9) also had a significantly higher risk of cancer-specific mortality (P = .021) compared to those with low (0/1+) NFATc1 muscle-invasive tumor (n = 53). Thus, compared with nonneoplastic urothelium, a significant increase in the expression of NFATc1 in UUTUC was seen, implying the involvement of NFATc1 signals in the development of UUTUC. The current results further suggest that NFATc1 overexpression serves as a predictor of poor prognosis in patients with UUTUC.     

A Long Noncoding RNA Signature That Predicts Pathological Complete Remission Rate Sensitively in Neoadjuvant Treatment of Breast Cancer

BACKGROUND: Mounting evidence suggests that long noncoding RNAs (lncRNAs) are closely related to pathological complete response (pCR) in neoadjuvant treatment of breast cancer. Here, we construct lncRNA associated models to predict pCR rate. METHODS: LncRNA expression profiles of breast cancer patients treated with neoadjuvant chemotherapy (NAC) were obtained from Gene Expression Omnibus by repurposing existing microarray data. The prediction model was firstly built by analyzing the correlation between pCR and lncRNA expression in the discovery dataset GSE 25066 (n = 488). Another three independent datasets, GSE20194 (n = 278), GSE20271 (n = 178), and GSE22093 (n = 97), were integrated as the validation cohort to assess the prediction efficiency. RESULTS: A novel lncRNA signature (LRS) consisting of 36 lncRNAs was identified. Based on this LRS, patients with NAC treatment were divided into two groups: LRS-high group and LRS-low group, with positive correlation of pCR rate in the discovery dataset. In the validation cohort, univariate and multivariate analyses both demonstrated that high LRS was associated with higher pCR rate. Subgroup analysis confirmed that this model performed well in luminal B [odds ratio (OR) = 5.4; 95% confidence interval (CI) = 2.7-10.8; P = 1.47e-06], HER2-enriched (OR = 2.5; 95% CI = 1.1-5.7; P = .029), and basal-like (OR = 5.5; 95% CI = 2.3-16.2; P = 5.32e-04) subtypes. Compared with other preexisting prediction models, LRS demonstrated better performance with higher area under the curve. Functional annotation analysis suggested that lncRNAs in this signature were mainly involved in cancer proliferation process. CONCLUSION: Our findings indicated that our lncRNA signature was sensitive to predict pCR rate in the neoadjuvant treatment of breast cancer, which deserves further evaluation.     

Prognostic Value of Molecular Breast Cancer Subtypes based on Her2, ESR1, PGR and Ki67 mRNA-Expression in Muscle Invasive Bladder Cancer

INTRODUCTION: Gene expression analyses have identified similarities between bladder and breast cancer, where clinical risk stratification is based on Her2, ESR1, PGR and Ki67 expression. The aim of the study was to assess the respective marker gene expression in patients treated with radical cystectomy for muscle-invasive bladder cancer (MIBC) and to evaluate the applicability of breast cancer subtypes for MIBC risk stratification. MATERIALS & METHODS: 102 patients treated with radical cystectomy for MIBC were assessed. Using routine FFPE tissue and an IVD validated kit, mRNA expression was measured by single step RT-qPCR. Partition test were employed to define cut-off values for high or low marker gene expression. Association of expression with outcome was assessed using Kaplan-Meier analysis and multivariate cox regression analysis. Finally, we performed validation of our results in the MD-Anderson cohort (n = 57). RESULTS: Cancer specific survival (CSS) was impaired in patients with high gene expression of Her2 (P = 0.0009) and ESR1 (P = 0.04). In the multivariate regression model Her2 expression remained significant for the prediction of CSS (HR = 2.11, CI 1.11-4.21, P = 0.024). Furthermore, molecular stratification by breast cancer subgroups was significant (P = 0.023) for CSS prediction. Especially the differentiation between Her2-positive and Luminal A (HR = 4.41, CI 1.53-18.71, P = 0.004) and Luminal B (HR = 1.96, CI 0.99-4.08, P = 0.053) respectively was an independent prognostic parameter for CSS. External validation resulted in comparable risk stratification with differences in fractional subgroups distribution. CONCLUSION: Gene expression of Her2, ESR1, PGR, Ki67 and corresponding breast cancer subtypes allow a risk-stratification in MIBC, whereby Her2 overexpressing tumors reveal a particularly poor prognosis.     

Proving of a Mathematical Model of Cell Calculation Based on Apparent Diffusion Coefficient

OBJECTIVES: Recently, Atuegwu et al. proposed a mathematical model based on ADC_mean and ADC_min to calculation of cellularity. Our purpose was to compare the calculated cellularity according to the formula with the estimated cell count by histopathology in different tumors. METHODS: For this study, we re-analyzed our previous data regarding associations between ADC parameters and histopathological findings. Overall, 134 patients with different tumors were acquired for the analysis. For all tumors, the number of tumor cells was calculated according to Atuegwu et al. 2013. We performed a correlation analysis between the calculated and estimated cellularity. Thereby, Pearson's correlation coefficient was used and P < .05 was taken to indicate statistical significance in all instances. RESULTS: The estimated and calculated cellularity correlated well together in HNSCC (r = 0.701, P = .016) and lymphomas (r = 0.661, P = .001), and moderately in rectal cancer (r = 0.510, P = .036). There were no statistically significant correlations between the estimated and calculated cellularity in uterine cervical cancer, meningiomas, and in thyroid cancer. CONCLUSION: The proposed formula for cellularity calculation does not apply for all tumors. It may be used for HNSCC, cerebral lymphomas and rectal cancer, but not for uterine cervical cancer, meningioma, and thyroid cancer. Furthermore, its usefulness should be proved for other tumors.     

Thymic tumors and results of radiotherapy

Aim

The aim of this study was to evaluate thymic epithelial tumors (TETs) for treatment outcomes and prognostic factors on survival.

Detection of KRAS Exon 2 Mutations in Circulating Tumor Cells Isolated by the ISET System from Patients with RAS Wild Type Metastatic Colorectal Cancer

INTRODUCTION: The presence of KRAS mutations in patients with metastatic colorectal cancer (mCRC) predicts poor response to agents targeting the EGFR. Even in patients with RAS wild type (WT) tumors, resistance eventually develops due to multiple mechanisms, including the expansion of previously undetected KRAS mutated clones. In this feasibility study, we aimed to detect KRAS exon 2 mutations in serial samples of circulating tumor cells (CTCs) of RAS WT patients with mCRC captured by the Isolation by Size of Epithelial Tumor cells (ISET) system. METHODS: CTC isolation using the ISET system was performed from prospectively collected blood samples obtained from patients with RAS and BRAF WT mCRC prior to first-line therapy initiation, at first imaging assessment and on disease progression. CTCs were enumerated using hematoxylin & eosin and CD45 double stain on a single membrane spot. DNA was extracted from 5 spots and KRAS exon 2 mutations were detected using a custom quantitative Polymerase Chain Reaction (qPCR) assay. RESULTS: Fifteen patients were enrolled and 28 blood samples were analyzed. In 9 (60%) patients, at least one sample was positive for the presence of a KRAS exon 2 mutation. In 11 out of 28 samples (39.2%) with detectable CTCs a KRAS mutation was detected; the corresponding percentages for baseline and on progression samples were 27% and 37.5%, respectively. The most commonly detected mutations were G13D and G12C (n = 3). The presence of KRAS mutated CTCs at baseline was not prognostic for either PFS (P = .950) or OS (P = .383). CTC kinetics did not follow tumor response patterns. CONCLUSION: The results demonstrate that using a qPCR-based assay, KRAS exon 2 mutations could be detected in CTCs captured by the ISET system from patients with RAS WT primary tumors. However, the clinical relevance of these CTCs remains to be determined in future studies.     

Anthropometric dimensions provide reliable estimates of abdominal adiposity: A validation study

Abdominal fat accumulation is a major risk factor for cardiometabolic morbidity and mortality. The purpose of the study is to assess the possibility of developing accurate estimation equations based on body measurements to determine total abdominal (TFA), subcutaneous (SFA) and visceral fat area (VFA). Hungarian volunteers (n = 198) aged between 20 and 81 years were enrolled in the study, which was conducted between July and November 2014. All persons underwent anthropometric measurements and computer tomographic (CT) scanning. Sex-specific multiple linear regression analyses were conducted in a subgroup of 98 participants to generate estimation models, then Bland–Altman's analyses were applied in the cross-validation group to compare their predictive efficiency. The variables best predicting VFA were hip circumference, calf circumference and waist-to-hip ratio (WHR) for males (R² = 0.713; SEE = 5602.1 mm²) and sagittal abdominal diameter (SAD), WHR, thigh circumference and triceps skinfold for females (R² = 0.845; SEE = 3835.6 mm²). The SFA prediction equation included SAD, thigh circumference and abdominal skinfold for males (R² = 0.848; SEE = 4124.1 mm²), body mass index and thigh circumference for females (R² = 0.861; SEE = 5049.7 mm²). Prediction accuracy was the highest in the case of TFA: hip circumference and WHR for males (R² = 0.910; SEE = 5637.2 mm²), SAD, thigh circumference and abdominal skinfold for females (R² = 0.915; SEE = 6197.5 mm²) were used in the equations. The results suggested that deviations in the predictions were independent of the amount of adipose tissue. Estimation of abdominal fat depots based on anthropometric traits could provide a cheap, reliable method in epidemiologic research and public health screening to evaluate the risk of cardiometabolic events.     

2D and 3D CT Radiomics Features Prognostic Performance Comparison in Non-Small Cell Lung Cancer

OBJECTIVE: To compare 2D and 3D radiomics features prognostic performance differences in CT images of non-small cell lung cancer (NSCLC). METHOD: We enrolled 588 NSCLC patients from three independent cohorts. Two sets of 463 patients from two different institutes were used as the training cohort. The remaining cohort with 125 patients was set as the validation cohort. A total of 1014 radiomics features (507 2D features and 507 3D features correspondingly) were assessed. Based on the dichotomized survival data, 2D and 3D radiomics indicators were calculated for each patient by trained classifiers. We used the area under the receiver operating characteristic curve (AUC) to assess the prediction performance of trained classifiers (the support vector machine and logistic regression). Kaplan–Meier and Cox hazard survival analyses were also employed. Harrell's concordance index (C-Index) and Akaike's information criteria (AIC) were applied to assess the trained models. RESULTS: Radiomics indicators were built and compared by AUCs. In the training cohort, 2D_AUC = 0.653, 3D_AUC = 0.671. In the validation cohort, 2D_AUC = 0.755, 3D_AUC = 0.663. Both 2D and 3D trained indicators achieved significant results (P < .05) in the Kaplan-Meier analysis and Cox regression. In the validation cohort, 2D Cox model had a C-Index = 0.683 and AIC = 789.047; 3D Cox model obtained a C-Index = 0.632 and AIC = 799.409. CONCLUSION: Both 2D and 3D CT radiomics features have a certain prognostic ability in NSCLC, but 2D features showed better performance in our tests. Considering the cost of the radiomics features calculation, 2D features are more recommended for use in the current study.     

The Pulmonary Fibrosis Associated MUC5B Promoter Polymorphism Is Prognostic of the Overall Survival in Patients with Non–Small Cell Lung Cancer (NSCLC) Receiving Definitive Radiotherapy

BACKGROUND: MUC5B is glycoprotein secreted by bronchial glands. A promoter variant in MUC5B, rs35705950, was previously found to be strongly associated with the incidence of idiopathic pulmonary fibrosis (IPF) and also the overall survival (OS) of such patients. Patients with IPF and patients with radiation pneumonitis (RP) have the similar pathologic process and clinical symptoms. However, the role of rs35705950 in patients receiving thoracic radiotherapy remains unclear. PATIENTS AND METHODS: In total, 664 patients with NSCLC receiving definitive radiotherapy (total dose ≥60 Gy) were included in our study. RP was scored via the Common Terminology Criteria for Adverse Events v3.0. OS was the second end point. MUC5B rs35705950 was genotyped, and Kaplan-Meier and Cox regression analyses were used to evaluate associations between MUC5B rs35705950 and the risk of RP or OS. RESULTS: The median patient age was 66 years (range 35-88); most (488 [73.2%]) had stage III of the disease. Until the last follow-up, 250 patients developed grade ≥ 2 RP, 82 patients developed grade ≥ 3 RP, and 440 patients died. The median mean lung dose was 17.9 Gy (range 0.15-32.74). No statistically significant associations were observed between genotypes of MUC5B rs35705950 and the incidence of RP ≥ grade 2 either in univariate analysis (hazard ratio [HR] 1.009, 95% confidence interval [CI] 0.728-1.399, P = .958) or in multivariate analysis (HR 0.921, 95% CI 0.645-1.315, P = .65). Similar results were also observed for RP ≥ grade 3, while TT/GT genotypes in MUC5B were significantly associated with poor OS in both univariate analysis (HR 1.287, 95% CI 1.009-1.640, P = .042) and multivariate analysis (HR 1.561, 95% CI 1.193-2.042, P = .001). CONCLUSION: MUC5B promoter polymorphism could be prognostic of the OS among NSCLC patients receiving definitive radiotherapy, although no significant associations were found with the risk of RP.     

Treatment outcomes and late toxicities of intensity-modulated radiation therapy for 1091 Japanese patients with localized prostate cancer

Aim

This study aimed to evaluate the treatment result of intensity-modulated radiation therapy (IMRT) in a large number of Japanese patients with prostate cancer.

Integrated Analysis Identifies Molecular Signatures and Specific Prognostic Factors for Different Gastric Cancer Subtypes

BACKGROUND: Gastric cancer (GC) is the fifth leading cause of cancer-related deaths worldwide. As an effective and easily performed method, microscopy-based Lauren classification has been widely accepted by gastrointestinal surgeons and pathologists for GC subtyping, but molecular characteristics of different Lauren subtypes were poorly revealed. METHODS: GSE62254 was used as a derivation cohort, and GSE15459 was used as a validation cohort. The difference between diffuse and intestinal GC on the gene expression level was measured. Gene ontology (GO) enrichment analysis was performed for both subgroups. Hierarchical clustering and heatmap exhibition were also performed. Kaplan-Meier plot and Cox proportional hazards model were used to evaluate survival grouped by the given genes or hierarchical clusters. RESULTS: A total of 4598 genes were found differentially expressed between diffuse and intestinal GC. Immunity- and cell adhesion–related GOs were enriched for diffuse GC, whereas DNA repair– and cell cycle–related GOs were enriched for intestinal GC. We proposed a 40-gene signature (χ² = 30.71, P < .001) that exhibits better discrimination for prognosis than Lauren classification (χ² = 12.11, P = .002). FRZB [RR (95% CI) = 1.824 (1.115-2.986), P = .017] and EFEMP1 [RR (95% CI) = 1.537 (0.969-2.437), P = .067] were identified as independent prognostic factors only in diffuse GC but not in intestinal GC patients. KRT23 [RR (95% CI) = 1.616 (0.938-2.785), P = .083] was identified as an independent prognostic factor only in intestinal GC patients but not in diffuse GC patients. Similar results were achieved in the validation cohort. CONCLUSION: We found that GCs with different Lauren classifications had different molecular characteristics and identified FRZB, EFEMP1, and KRT23 as subtype-specific prognostic factors for GC patients.     

Patterns of treatment failure in salivary gland cancers

Aim

The purpose of the study was to publish our experience of salivary gland cancer treatment with large number of patients treated at a single institution.

Identification of target proteins of clinical immunity to Plasmodium falciparum in a region of low malaria transmission

The target molecules of antibodies against falciparum malaria remain largely unknown. Recently we have identified multiple proteins as targets of immunity against Plasmodium falciparum using African serum samples. To investigate whether potential targets of clinical immunity differ with transmission intensity, we assessed immune responses in residents of low malaria transmission region in Thailand. Malaria asymptomatic volunteers (Asy: n = 19) and symptomatic patients (Sym: n = 21) were enrolled into the study. Serum immunoreactivity to 186 wheat germ cell-free system (WGCFS)-synthesized recombinant P. falciparum asexual-blood stage proteins were determined by AlphaScreen, and subsequently compared between the study groups. Forty proteins were determined as immunoreactive with antibody responses to 35 proteins being higher in Asy group than in Sym group. Among the 35 proteins, antibodies to MSP3, MSPDBL1, RH2b, and MSP7 were significantly higher in Asy than Sym (unadjusted p < 0.005) suggesting these antigens may have a protective role in clinical malaria. MSP3 reactivity remained significantly different between Asy and Sym groups even after multiple comparison adjustments (adjusted p = 0.033). Interestingly, while our two preceding studies using African sera were conducted differently (e.g., cross-sectional vs. longitudinal design, observed clinical manifestation vs. functional activity), those studies similarly identified MSP3 and MSPDBL1 as potential targets of protective immunity. This study further provides a strong rationale for the application of WGCFS-based immunoprofiling to malaria vaccine candidate and biomarker discovery even in low or reduced malaria transmission settings.     

Association Between Background Parenchymal Enhancement and Pathologic Complete Remission Throughout the Neoadjuvant Chemotherapy in Breast Cancer Patients

PURPOSE: To retrospectively investigate the quantitative background parenchymal enhancement (BPE) of the contralateral normal breast in patients with unilateral invasive breast cancer throughout multiple monitoring points of neoadjuvant chemotherapy (NAC) and to further determine whether BPE is associated with tumor response, especially at the early stage of NAC. MATERIALS AND METHODS: A total of 90 patients with unilateral breast cancer who then received six or eight cycles of NAC before surgery were analyzed retrospectively. BPE was measured in dynamic contrast-enhanced MRI at baseline and after 2nd, 4th, and 6th NAC, respectively. Correlation between BPE and tumor size was analyzed, and the association between pathologic complete remission (pCR) and BPE was also analyzed. RESULTS: The BPE of contralateral normal breast showed a constant reduction throughout NAC therapy regardless of the menopausal status (P < .001 in all). Both the BPEs and the changes of BPE in each of the three monitoring points were significantly correlated with those in tumor size (P < .05 in all), and the reduction of BPE after 2nd NAC had the largest diagnostic value for pCR (AUC = 0.726, P < .001), particularly in hormonal receptor (HR)-negative patients (OR = 0.243, 95%CI = 0.083 to 0.706, P = .009). CONCLUSION: The BPE of contralateral normal breast had a constant decreased tendency similar to the change of tumor size in NAC. Reduction of BPE at the early stage of NAC was positively associated with pCR, especially in HR-negative status.     

Analysis of Chemopredictive Assay for Targeting Cancer Stem Cells in Glioblastoma Patients

Introduction: The prognosis of glioblastoma (GBM) treated with standard-of-care maximal surgical resection and concurrent adjuvant temozolomide (TMZ)/radiotherapy remains very poor (less than 15 months). GBMs have been found to contain a small population of cancer stem cells (CSCs) that contribute to tumor propagation, maintenance, and treatment resistance. The highly invasive nature of high-grade gliomas and their inherent resistance to therapy lead to very high rates of recurrence. For these reasons, not all patients with similar diagnoses respond to the same chemotherapy, schedule, or dose. Administration of ineffective anticancer therapy is not only costly but more importantly burdens the patient with unnecessary toxicity and selects for the development of resistant cancer cell clones. We have developed a drug response assay (ChemoID) that identifies the most effective chemotherapy against CSCs and bulk of tumor cells from of a panel of potential treatments, offering great promise for individualized cancer management. Providing the treating physician with drug response information on a panel of approved drugs will aid in personalized therapy selections of the most effective chemotherapy for individual patients, thereby improving outcomes. A prospective study was conducted evaluating the use of the ChemoID drug response assay in GBM patients treated with standard of care. Methods: Forty-one GBM patients (mean age 54 years, 59% male), all eligible for a surgical biopsy, were enrolled in an Institutional Review Board–approved protocol, and fresh tissue samples were collected for drug sensitivity testing. Patients were all treated with standard-of-care TMZ plus radiation with or without maximal surgery, depending on the status of the disease. Patients were prospectively monitored for tumor response, time to recurrence, progression-free survival (PFS), and overall survival (OS). Odds ratio (OR) associations of 12-month recurrence, PFS, and OS outcomes were estimated for CSC, bulk tumor, and combined assay responses for the standard-of-care TMZ treatment; sensitivities/specificities, areas under the curve (AUCs), and risk reclassification components were examined. Results: Median follow-up was 8 months (range 3-49 months). For every 5% increase in in vitro CSC cell kill by TMZ, 12-month patient response (nonrecurrence of cancer) increased two-fold, OR = 2.2 (P = .016). Similar but somewhat less supported associations with the bulk tumor test were seen, OR = 2.75 (P = .07) for each 5% bulk tumor cell kill by TMZ. Combining CSC and bulk tumor assay results in a single model yielded a statistically supported CSC association, OR = 2.36 (P = .036), but a much attenuated remaining bulk tumor association, OR = 1.46 (P = .472). AUCs and [sensitivity/specificity] at optimal outpoints (>40% CSC cell kill and >55% bulk tumor cell kill) were AUC = 0.989 [sensitivity = 100/specificity = 97], 0.972 [100/89], and 0.989 [100/97] for the CSC only, bulk tumor only, and combined models, respectively. Risk categorization of patients was improved by 11% when using the CSC test in conjunction with the bulk test (risk reclassification nonevent net reclassification improvement [NRI] and overall NRI = 0.111, P = .030). Median recurrence time was 20 months for patients with a positive (>40% cell kill) CSC test versus only 3 months for those with a negative CSC test, whereas median recurrence time was 13 months versus 4 months for patients with a positive (>55% cell kill) bulk test versus negative. Similar favorable results for the CSC test were observed for PFS and OS outcomes. Panel results across 14 potential other treatments indicated that 34/41 (83%) potentially more optimal alternative therapies may have been chosen using CSC results, whereas 27/41 (66%) alternative therapies may have been chosen using bulk tumor results. Conclusions: The ChemoID CSC drug response assay has the potential to increase the accuracy of bulk tumor assays to help guide individualized chemotherapy choices. GBM cancer recurrence may occur quickly if the CSC test has a low in vitro cell kill rate even if the bulk tumor test cell kill rate is high.     

MRI Texture Analysis Reflects Histopathology Parameters in Thyroid Cancer – A First Preliminary Study

OBJECT: Thyroid cancer represents the most frequent malignancy of the endocrine system with an increasing incidence worldwide. Novel imaging techniques are able to further characterize tumors and even predict histopathology features. Texture analysis is an emergent imaging technique to extract extensive data from an radiology images. The present study was therefore conducted to identify possible associations between texture analysis and histopathology parameters in thyroid cancer. METHODS: The radiological database was retrospectively reviewed for thyroid carcinoma. Overall, 13 patients (3 females, 23.1%) with a mean age of 61.6 years were identified. The MaZda program was used for texture analysis. The T1-precontrast and T2-weighted images were analyzed and overall 279 texture feature for each sequence was investigated. For every patient cell count, Ki67-index and p53 count were investigated. RESULTS: Several significant correlations between texture features and histopathology were identified. Regarding T1-weighted images, S(0;1)Sum Averg correlated the most with cell count (r = 0.82). An inverse correlations with S(5;0)AngScMom, S(5;0)DifVarnc S(5;0), DiffEntrp and GrNonZeros (r = ?0.69, ?0.66, ?0.69 and ?0.63, respectively) was also identified. For T2-weighted images, Variance with r = 0.63 was the highest coefficient, WavEnLL_S3 correlated inversely with cell count (r = ?0.57). WavEnLL_S2 derived from T1-weighted images was the highest coefficient r = ?0.80, S(0;5)SumVarnc was positively with r = 0.74. Regarding T2-weighted images WavEnHL_s-1 was inverse correlated with Ki67 index (r = ?0.77). S(1;0)Correlat was with r = 0.75 the best correlation with Ki67 index. For T1-weighed images S(5;0)SumofSqs was the best with r = 0.65 with p53 count. For T2-weighted images S(1;?1)SumEntrp was the inverse correlation with r = ?0.72, whereas S(0;4)AngScMom correlated positively with r = 0.63. CONCLUSIONS: MRI texture analysis derived from conventional sequences reflects histopathology features in thyroid cancer. This technique might be a novel noninvasive modality to further characterize thyroid cancer in clinical oncology.