ABSTRACT
Psychiatric disturbances are commonly associated with adult-onset isolated dystonia (AOID); however, the mechanisms underlying psychiatric abnormalities in AOID remain unknown. We aimed to investigate the structural and functional brain changes in AOID patients with anxiety, and identify imaging biomarkers for diagnosing anxiety. Structural and functional magnetic resonance was performed on 69 AOID patients and 35 healthy controls (HCs). The Hamilton Anxiety Scale (HAMA) was used to assess anxiety symptoms in AOID patients and assign patients to AOID with and without anxiety groups. Group differences in grey matter volume, amplitude of low-frequency fluctuations (ALFF), fractional ALFF, and regional homogeneity (ReHo) were evaluated. Area under the receiver operating characteristic curve (ROC AUC) was used as a metric to identify imaging biomarkers for diagnosing anxiety. AOID patients with anxiety exhibited an increased ALFF and ReHo in the left angular gyrus (ANG.L) compared with those without and HCs (voxel P<0.001 and cluster P<0.05, corrected using GRF). A significant positive correlation was observed between ALFF (r = 0.627, P<0.001) and ReHo (r = 0.515, P<0.001) in the ANG.L and HAMA scores in AOID patients. ALFF and ReHo in the ANG.L exhibited an ROC AUC of 0.904 and 0.851, respectively, in distinguishing AOID patients with anxiety from those without and an ROC AUC of 0.887 and 0.853, respectively, in distinguishing AOID patients with anxiety from HCs. These findings provide new insights into the pathophysiology of psychiatric disturbances and highlight potential candidate biomarkers for identifying anxiety in AOID patients.
ABSTRACT
BACKGROUND: Accumulating neuroimaging evidence indicates that patients with cervical dystonia (CD) have changes in the cortico-subcortical white matter (WM) bundle. However, whether these patients' WM structural networks undergo reorganization remains largely unclear. We aimed to investigate topological changes in large-scale WM structural networks in patients with CD compared to healthy controls (HCs), and explore the network changes associated with clinical manifestations. METHODS: Diffusion tensor imaging (DTI) was conducted in 30 patients with CD and 30 HCs, and WM network construction was based on the BNA-246 atlas and deterministic tractography. Based on the graph theoretical analysis, global and local topological properties were calculated and compared between patients with CD and HCs. Then, the AAL-90 atlas was used for the reproducibility analyses. In addition, the relationship between abnormal topological properties and clinical characteristics was analyzed. RESULTS: Compared with HCs, patients with CD showed changes in network segregation and resilience, characterized by increased local efficiency and assortativity, respectively. In addition, a significant decrease of network strength was also found in patients with CD relative to HCs. Validation analyses using the AAL-90 atlas similarly showed increased assortativity and network strength in patients with CD. No significant correlations were found between altered network properties and clinical characteristics in patients with CD. CONCLUSION: Our findings show that reorganization of the large-scale WM structural network exists in patients with CD. However, this reorganization is attributed to dystonia-specific abnormalities or hyperkinetic movements that need further identification.
Subject(s)
Diffusion Tensor Imaging , Torticollis , White Matter , Humans , Torticollis/diagnostic imaging , Torticollis/pathology , White Matter/diagnostic imaging , White Matter/pathology , Female , Male , Diffusion Tensor Imaging/methods , Middle Aged , Adult , Nerve Net/diagnostic imaging , Nerve Net/pathology , AgedABSTRACT
BACKGROUND: The thalamus has a central role in the pathophysiology of idiopathic cervical dystonia (iCD); however, the nature of alterations occurring within this structure remain largely elusive. Using a structural magnetic resonance imaging (MRI) approach, we examined whether abnormalities differ across thalamic subregions/nuclei in patients with iCD. METHODS: Structural MRI data were collected from 37 patients with iCD and 37 healthy controls (HCs). Automatic parcellation of 25 thalamic nuclei in each hemisphere was performed based on the FreeSurfer program. Differences in thalamic nuclei volumes between groups and their relationships with clinical information were analysed in patients with iCD. RESULTS: Compared to HCs, a significant reduction in thalamic nuclei volume primarily in central medial, centromedian, lateral geniculate, medial geniculate, medial ventral, paracentral, parafascicular, paratenial, and ventromedial nuclei was found in patients with iCD (P < 0.05, false discovery rate corrected). However, no statistically significant correlations were observed between altered thalamic nuclei volumes and clinical characteristics in iCD group. CONCLUSION: This study highlights the neurobiological mechanisms of iCD related to thalamic volume changes.
Subject(s)
Magnetic Resonance Imaging , Thalamus , Torticollis , Humans , Male , Female , Middle Aged , Torticollis/diagnostic imaging , Torticollis/pathology , Magnetic Resonance Imaging/methods , Thalamus/diagnostic imaging , Thalamus/pathology , Adult , Aged , Thalamic Nuclei/diagnostic imaging , Thalamic Nuclei/pathologyABSTRACT
The thalamus is considered a key region in the neuromechanisms of blepharospasm. However, previous studies considered it as a single, homogeneous structure, disregarding potentially useful information about distinct thalamic nuclei. Herein, we aimed to examine (i) whether grey matter volume differs across thalamic subregions/nuclei in patients with blepharospasm and blepharospasm-oromandibular dystonia; (ii) causal relationships among abnormal thalamic nuclei; and (iii) whether these abnormal features can be used as neuroimaging biomarkers to distinguish patients with blepharospasm from blepharospasm-oromandibular dystonia and those with dystonia from healthy controls. Structural MRI data were collected from 56 patients with blepharospasm, 20 with blepharospasm-oromandibular dystonia and 58 healthy controls. Differences in thalamic nuclei volumes between groups and their relationships to clinical information were analysed in patients with dystonia. Granger causality analysis was employed to explore the causal effects among abnormal thalamic nuclei. Support vector machines were used to test whether these abnormal features could distinguish patients with different forms of dystonia and those with dystonia from healthy controls. Compared with healthy controls, patients with blepharospasm exhibited reduced grey matter volume in the lateral geniculate and pulvinar inferior nuclei, whereas those with blepharospasm-oromandibular dystonia showed decreased grey matter volume in the ventral anterior and ventral lateral anterior nuclei. Atrophy in the pulvinar inferior nucleus in blepharospasm patients and in the ventral lateral anterior nucleus in blepharospasm-oromandibular dystonia patients was negatively correlated with clinical severity and disease duration, respectively. The proposed machine learning scheme yielded a high accuracy in distinguishing blepharospasm patients from healthy controls (accuracy: 0.89), blepharospasm-oromandibular dystonia patients from healthy controls (accuracy: 0.82) and blepharospasm from blepharospasm-oromandibular dystonia patients (accuracy: 0.94). Most importantly, Granger causality analysis revealed that a progressive driving pathway from pulvinar inferior nuclear atrophy extends to lateral geniculate nuclear atrophy and then to ventral lateral anterior nuclear atrophy with increasing clinical severity in patients with blepharospasm. These findings suggest that the pulvinar inferior nucleus in the thalamus is the focal origin of blepharospasm, extending to pulvinar inferior nuclear atrophy and subsequently extending to the ventral lateral anterior nucleus causing involuntary lower facial and masticatory movements known as blepharospasm-oromandibular dystonia. Moreover, our results also provide potential targets for neuromodulation especially deep brain stimulation in patients with blepharospasm and blepharospasm-oromandibular dystonia.
ABSTRACT
Blepharospasm is a common form of focal dystonia characterized by excessive and involuntary spasms of the orbicularis oculi. In addition to idiopathic blepharospasm, lesions in various brain regions can also cause acquired blepharospasm. Whether these two types of blepharospasm share a common brain network remains largely unknown. Herein, we performed lesion coactivation network mapping, based on meta-analytic connectivity modeling, to test whether lesions causing blepharospasm could be mapped to a common coactivation brain network. We then tested the abnormality of the network in patients with idiopathic blepharospasm (n = 42) compared with healthy controls (n = 44). We identified 21 cases of lesion-induced blepharospasms through a systematic literature search. Although these lesions were heterogeneous, they were part of a co-activated brain network that mainly included the bilateral supplementary motor areas. Coactivation of these regions defines a single brain network that encompasses or is adjacent to most heterogeneous lesions causing blepharospasm. Moreover, the bilateral supplementary motor area is primarily associated with action execution, visual motion, and imagination, and participates in finger tapping and saccades. They also reported decreased functional connectivity with the left posterior cingulate cortex in patients with idiopathic blepharospasm. These results demonstrate a common convergent abnormality of the supplementary motor area across idiopathic and acquired blepharospasms, providing additional evidence that the supplementary motor area is an important brain region that is pathologically impaired in patients with blepharospasm.
Subject(s)
Blepharospasm , Motor Cortex , Humans , Blepharospasm/physiopathology , Blepharospasm/etiology , Female , Motor Cortex/physiopathology , Motor Cortex/diagnostic imaging , Male , Aged , Middle Aged , Magnetic Resonance Imaging , Nerve Net/physiopathology , Nerve Net/diagnostic imagingABSTRACT
Selective impairment in recognizing facial expressions of disgust was reported in patients with focal dystonia several years ago, but the basic neural mechanisms remain largely unexplored. Therefore, we investigated whether dysfunction of the brain network involved in disgust recognition processing was related to this selective impairment in blepharospasm. Facial emotion recognition evaluations and resting-state functional magnetic resonance imaging were performed in 33 blepharospasm patients and 33 healthy controls (HCs). The disgust processing network was constructed, and modularity analyses were performed to identify sub-networks. Regional functional indexes and intra- and inter-functional connections were calculated and compared between the groups. Compared to HCs, blepharospasm patients demonstrated a worse performance in disgust recognition. In addition, functional connections within the sub-network involved in perception processing rather than recognition processing of disgust were significantly decreased in blepharospasm patients compared to HCs. Specifically, decreased functional connections were noted between the left fusiform gyrus (FG) and right middle occipital gyrus (MOG), the left FG and right FG, and the right FG and left MOG. We identified decreased functional activity in these regions, as indicated by a lower amplitude of low-frequency fluctuation in the left MOG, fractional amplitude of low-frequency fluctuation in the right FG, and regional homogeneity in the right FG and left MOG in blepharospasm patients versus HCs. Our results suggest that dysfunctions of the disgust processing network exist in blepharospasm. A deficit in disgust emotion recognition may be attributed to disturbances in the early perception of visual disgust stimuli in blepharospasm patients.
Subject(s)
Blepharospasm , Facial Recognition , Humans , Blepharospasm/diagnostic imaging , Emotions , Brain/diagnostic imaging , Magnetic Resonance Imaging , Brain Mapping , Facial ExpressionABSTRACT
Background: Structural changes occur in brain regions involved in cortico-basal ganglia networks in idiopathic blepharospasm (iBSP); whether these changes influence the function connectivity patterns of cortico-basal ganglia networks remains largely unknown. Therefore, we aimed to investigate the global integrative state and organization of functional connections of cortico-basal ganglia networks in patients with iBSP. Methods: Resting-state functional magnetic resonance imaging data and clinical measurements were acquired from 62 patients with iBSP, 62 patients with hemifacial spasm (HFS), and 62 healthy controls (HCs). Topological parameters and functional connections of cortico-basal ganglia networks were evaluated and compared among the three groups. Correlation analyses were performed to explore the relationship between topological parameters and clinical measurements in patients with iBSP. Results: We found significantly increased global efficiency and decreased shortest path length and clustering coefficient of cortico-basal ganglia networks in patients with iBSP compared with HCs, however, such differences were not observed between patients with HFS and HCs. Further correlation analyses revealed that these parameters were significantly correlated with the severity of iBSP. At the regional level, the functional connectivity between the left orbitofrontal area and left primary somatosensory cortex and between the right anterior part of pallidum and right anterior part of dorsal anterior cingulate cortex was significantly decreased in patients with iBSP and HFS compared with HCs. Conclusion: Dysfunction of the cortico-basal ganglia networks occurs in patients with iBSP. The altered network metrics of cortico-basal ganglia networks might be served as quantitative markers for evaluation of the severity of iBSP.
ABSTRACT
Blepharospasm is traditionally thought to be a movement disorder that results from basal ganglia dysfunction. Recently, accumulating morphometric studies have revealed structural alterations outside the basal ganglia, such as in the brainstem, cerebellum and sensorimotor cortex, suggesting that blepharospasm may result from network disorders. However, the temporal and causal relationships between structural alterations and whether there are disease duration-related hierarchical structural changes in these patients remain largely unknown. Structural MRI was performed in 62 patients with blepharospasm, 62 patients with hemifacial spasm and 62 healthy controls to assess the structural alterations using voxel-based morphology and structural covariance networks. The use of the causal structural covariance network, modularity analysis and functional decoding were subsequently performed to map the causal effect of grey matter change pattern, hierarchical topography and functional characterizations of the structural network throughout the disease duration of blepharospasm. Greater grey matter volume in the left and right supplementary motor areas was identified in patients with blepharospasm compared to that in patients with hemifacial spasm and healthy controls, whereas no significant difference was identified between patients with hemifacial spasm and healthy controls. In addition, increased grey matter volume covariance between the right supplementary motor area and right brainstem, left superior frontal gyrus, left supplementary motor area and left paracentral gyrus was found in patients with blepharospasm compared to healthy controls. Further causal structural covariance network, modularity analysis and functional decoding showed that the right supplementary motor area served as a driving core in patients with blepharospasm, extending greater grey matter volume to areas in the cortico-basal ganglia-brainstem motor pathway and cortical regions in the vision-motor integration pathway. Taken together, our results suggest that the right supplementary motor area is an early and important pathologically impaired region in patients with blepharospasm. With a longer duration of blepharospasm, increased grey matter volume extends from the right supplementary motor area to the cortico-basal ganglia motor and visual-motor integration pathways, showing a hierarchy of structural abnormalities in the disease progression of blepharospasm, which provides novel evidence to support the notion that blepharospasm may arise from network disorders and is associated with a wide range of grey matter abnormalities.
Subject(s)
Blepharospasm , Hemifacial Spasm , Motor Cortex , Humans , Motor Cortex/diagnostic imaging , Blepharospasm/diagnostic imaging , Brain , Gray Matter/diagnostic imaging , Magnetic Resonance ImagingABSTRACT
Background: Numerous factors are related to the prognosis of rectal cancer, including T stage, N stage, metastasis, extramural venous invasion (EMVI), circumferential resection margin (CRM), and tumor differentiation. However, it is still a challenge to precisely evaluate them before therapy; therefore, we investigate whether synthetic magnetic resonance imaging and apparent diffusion coefficient (ADC) values could help predict the prognostic factors of rectal cancer. Methods: Eighty-seven patients (55 men and 32 women; mean age, 59±11 years) with pathologically confirmed rectal cancer were enrolled. Preoperative quantitative metrics, including T1, T2, proton density (PD), and ADC values, were measured with diffusion-weighted imaging (DWI) acquired by a single-shot echo-planar sequence and synthetic magnetic resonance imaging acquired by a multi-dynamic multi-echo sequence at 3.0 T, in patients with rectal cancer by two radiologists. We evaluated the diagnostic performance of synthetic magnetic resonance imaging using the independent sample t-test or Mann-Whitney U test and receiver operating characteristic (ROC) curve and multivariate logistic regression analyses and compared the area under the ROC curve of quantitative values using the DeLong test. Results: The T2 and PD values showed a significant reduction among patients with poor differentiation and lymph node metastasis in rectal cancer. The area under the ROC curve values of T2 and PD values for predicting magnetic resonance imaging N stage and differentiation were 0.734, 0.682, and 0.673, 0.686, respectively. Moreover, combining T2 and PD values for magnetic resonance imaging N stage slightly improved the area under the ROC curve value of 0.774 (95% CI, 0.673-0.876). In the present study, the ADC and T1 values were not significant in the differentiation or clinical stage of rectal cancer (RC). Conclusions: Quantitative T2 and PD values obtained by synthetic magnetic resonance imaging might be used for evaluating prognostic factors of rectal cancer noninvasively. Furthermore, combining T2 and PD values further improved the diagnostic performance of magnetic resonance imaging N staging in rectal cancer. The ADC and T1 values were not significant in the differentiation or clinical stage of RC.
ABSTRACT
BACKGROUND AND PURPOSE: Accumulating evidence indicates that dynamic amplitude of low-frequency fluctuations (dALFF) or dynamic functional connectivity (dFC) can provide complementary information, distinct from static amplitude of low-frequency fluctuations (sALFF) or static functional connectivity (sFC), in detecting brain functional abnormalities in brain diseases. We aimed to examine whether dALFF and dFC can offer valuable information for the detection of functional brain abnormalities in patients with blepharospasm. METHODS: We collected resting-state functional magnetic resonance imaging data from 46 patients each of blepharospasm, hemifacial spasm (HFS), and healthy controls (HCs). We examined intergroup differences in sALFF and dALFF to investigate abnormal regional brain activity in patients with blepharospasm. Based on the dALFF results, we conducted seed-based sFC and dFC analyses to identify static and dynamic connectivity changes in brain networks centered on areas showing abnormal temporal variability of local brain activity in patients with blepharospasm. RESULTS: Compared with HCs, patients with blepharospasm displayed different brain functional change patterns characterized by increased sALFF in the left primary motor cortex (PMC) but increased dALFF variance in the right PMC. However, differences were not found between patients with HFS and HCs. Additionally, patients with blepharospasm exhibited decreased dFC strength, but no change in sFC, between right PMC and ipsilateral cerebellum compared with HCs; these findings were replicated when patients with blepharospasm were compared to those with HFS. CONCLUSIONS: Our findings highlight that dALFF and dFC are complementary to sALFF and sFC and can provide valuable information for detecting brain functional abnormalities in blepharospasm. Blepharospasm may be a network disorder involving the cortico-ponto-cerebello-thalamo-cortical circuit.
Subject(s)
Blepharospasm , Motor Cortex , Blepharospasm/diagnostic imaging , Brain/diagnostic imaging , Brain Mapping , Humans , Magnetic Resonance Imaging , Motor Cortex/diagnostic imagingABSTRACT
BACKGROUND: Accumulating evidence indicates regional structural changes in the white matter (WM) of brains in patients with blepharospasm (BSP); however, whether large-scale WM structural networks undergo widespread reorganization in these patients remains unclear. OBJECTIVE: We investigated topology changes and global and local features of large-scale WM structural networks in BSP patients compared with hemifacial spasm (HFS) patients or healthy controls (HCs). METHODS: This cross-sectional study applied graph theoretical analysis to assess deterministic diffusion tensor tractography findings in 41 BSP patients, 41 HFS patients, and 41 HCs. WM structural connectivity in 246 cortical and subcortical regions was assessed, and topological parameters of the resulting graphs were calculated. Networks were compared among BSP, HFS, and HCs groups. RESULTS: Compared to HCs, both BSP and HFS patients showed alterations in network integration and segregation characterized by increased global efficiency and modularity and reduced shortest path length. Moreover, increased nodal efficiency in multiple cortical and subcortical regions was found in BSP and HFS patients compared with HCs. However, these differences were not found between BSP and HFS patients. Whereas all participants showed highly similar hub distribution patterns, BSP patients had additional hub regions not present in either HFS patients or HCs, which were located in the primary head and face motor cortex and basal ganglia. CONCLUSIONS: Our findings suggest that the large-scale WM structural network undergoes an extensive reorganization in BSP, probably due to both dystonia-specific abnormalities and facial hyperkinetic movements. © 2021 International Parkinson and Movement Disorder Society.
Subject(s)
Blepharospasm , White Matter , Blepharospasm/diagnostic imaging , Brain/diagnostic imaging , Cross-Sectional Studies , Diffusion Tensor Imaging/methods , Humans , White Matter/diagnostic imagingABSTRACT
Accumulating diffusion tensor imaging (DTI) evidence suggests that white matter abnormalities evaluated by local diffusion homogeneity (LDH) or fractional anisotropy (FA) occur in patients with blepharospasm (BSP), both of which are significantly correlated with disease severity. However, whether the individual severity of BSP can be identified using these DTI metrics remains unknown. We aimed to investigate whether a combination of machine learning techniques and LDH or FA can accurately identify the individual severity of BSP. Forty-one patients with BSP were assessed using the Jankovic Rating Scale and DTI. The patients were assigned to non-functionally and functionally limited groups according to their Jankovic Rating Scale scores. A machine learning scheme consisting of beam search and support vector machines was designed to identify non-functionally versus functionally limited outcomes, with the input features being LDH or FA in 68 white matter regions. The proposed machine learning scheme with LDH or FA yielded an overall accuracy of 88.67 versus 85.19% in identifying non-functionally limited versus functionally limited outcomes. The scheme also identified a sensitivity of 91.40 versus 85.87% in correctly identifying functionally limited outcomes, a specificity of 83.33 versus 83.67% in accurately identifying non-functionally limited outcomes, and an area under the curve of 93.7 versus 91.3%. These findings suggest that a combination of LDH or FA measurements and a sophisticated machine learning scheme can accurately and reliably identify the individual disease severity in patients with BSP.
ABSTRACT
BACKGROUND: Previous studies have indicated that quantitative MRI (qMR) is beneficial for diagnosis of breast cancer. As a novel qMR technology, synthetic MRI (syMRI) may be advantageous by offering simultaneous generation of T1 and T2 mapping in one scan within a few minutes and without concern to the deposition of the gadolinium contrast agent in cell nucleus. In this study, the potential of quantitative mapping derived from Synthetic MRI (SyMRI) to diagnose breast cancer was investigated. METHODS: From April 2018 to May 2019, a total of 87 patients with suspicious breast lesions underwent both conventional and SyMRI before treatment. The quantitative metrics derived from SyMRI, including T1 and T2 values, were measured in breast lesions. The diagnostic performance of SyMRI was evaluated with unpaired Student's t-tests, receiver operating characteristic curve analysis and multivariate logistic regression analysis. The AUCs of quantitative values were compared using Delong test. RESULTS: Among 77 patients who met the inclusion criteria, 48 were diagnosed with histopathological confirmed breast cancers, and the rest had benign lesions. The breast cancers showed significantly higher T1 (1611.61 ± 215.88 ms) values and lower T2 (80.93 ± 7.51 ms) values than benign lesions. The area under the ROC curve (AUC) values were 0.931 (95% CI: 0.874-0.989) and 0.883 (95% CI: 0.810-0.956) for T1 and T2 maps, respectively, in diagnostic discrimination between breast cancers and benign lesions. A slightly increased AUC of 0.978 (95% CI: 0.915-0.993) was achieved by combining those two relaxation-based quantitative metrics. CONCLUSION: In conclusion, our preliminary study showed that the quantitative T1 and T2 values obtained by SyMRI could distinguish effectively between benign and malignant breast lesions, and T1 relaxation time showed the highest diagnostic efficiency. Furthermore, combining the two quantitative relaxation metrics further improved their diagnostic performance.
Subject(s)
Breast Neoplasms/diagnostic imaging , Breast/diagnostic imaging , Image Interpretation, Computer-Assisted/methods , Magnetic Resonance Imaging/methods , Adult , Aged , Area Under Curve , Contrast Media , Diagnosis, Differential , Female , Fibrocystic Breast Disease/diagnostic imaging , Humans , Middle Aged , Prospective Studies , ROC Curve , Regression Analysis , Sensitivity and SpecificityABSTRACT
White matter abnormalities in blepharospasm (BSP) have been evaluated using conventional intra-voxel metrics, and changes in patterns of cortical thickness in BSP remain controversial. We aimed to determine whether local diffusion homogeneity, an inter-voxel diffusivity metric, could be valuable in detecting white matter abnormalities for BSP; whether these changes are related to disease features; and whether cortical thickness changes occur in BSP patients. Diffusion tensor and structural magnetic resonance imaging were collected for 29 patients with BSP and 30 healthy controls. Intergroup diffusion differences were compared using tract-based spatial statistics analysis and measures of cortical thickness were obtained. The relationship among cortical thickness, diffusion metric in significantly different regions, and behavioral measures were further assessed. There were no significant differences in cortical thickness and fractional anisotropy between the groups. Local diffusion homogeneity was higher in BSP patients than controls, primarily in the left superior longitudinal fasciculus, corpus callosum, left posterior corona radiata, and left posterior thalamic radiata (P < 0.05, family-wise error corrected). The local diffusion homogeneity values in these regions were positively correlated with the Jankovic rating scale (r s = 0.416, P = 0.031) and BSP disability index (r s = 0.453, P = 0.018) in BSP patients. These results suggest that intra- and inter-voxel diffusive parameters are differentially sensitive to detecting BSP-related white matter abnormalities and that local diffusion homogeneity might be useful in assessing disability in BSP patients.
ABSTRACT
OBJECTIVE: To improve the assessment of primary tumor heterogeneity in magnetic resonance imaging (MRI) of non-small cell lung cancer (NSCLC), we proposed a method using basic measurements from T1- and T2-weighted MRI. METHODS: One hundred and four NSCLC patients with different T stages were studied. Fifty-two patients were analyzed as training group and another 52 as testing group. The ratios of standard deviation (SD)/mean signal value of primary tumor from T1-weighted (T1WI), T1-enhanced (T1C), T2-weighted (T2WI), and T2 fat suppression (T2fs) images were calculated. In the training group, correlation analyses were performed between the ratios and T stages. Then an ordinal regression model was built to generate the tumor heterogeneous index (THI) for evaluating the heterogeneity of tumor. The model was validated in the testing group. RESULTS: There were 11, 32, 40, and 21 patients with T1, T2, T3, and T4 disease, respectively. In the training group, the median SD/mean on T1WI, T1C, T2WI, and T2fs sequences was 0.11, 0.19, 0.16, and 0.15 respectively. The SD/mean on T1C (p=0.003), T2WI (p=0.000), and T2fs sequences (p=0.002) correlated significantly with T stages. Patients with more advanced T stage showed higher SD/mean on T2-weighted, T2fs, and T1C sequences. The median THI in the training group was 2.15. THI correlated with T stage significantly (p=0.000). In the testing group, THI was also significantly related to T stages (p=0.001). Higher THI had relevance to more advanced T stage. CONCLUSIONS: The proposed ratio measurements and THI based on MRI can serve as functional radiomic markers that correlated with T stages for evaluating heterogeneity of lung tumors.