- 1Department of Radiology, The First Affiliated Hospital of Shenzhen University, Shenzhen University, Shenzhen Second People’s Hospital, Shenzhen, Guangdong, China
- 2School of Information Science and Engineering, Shandong Normal University, Shandong, China
- 3Five-Year Clinical Medicine, Xiangya School of Medicine, Central South University, Changsha, Hunan, China
- 4Department of Biomedical Engineering, Sun Yat-sen University, Shenzhen, Guangdong, China
Purpose: Studying perivascular spaces (PVSs) is important for understanding the pathogenesis and pathological changes of neurological disorders. Although some methods for automated segmentation of PVSs have been proposed, most of them were based on 7T MR images that were majorly acquired in healthy young people. Notably, 7T MR imaging is rarely used in clinical practice. Herein, we propose a deep-learning-based method that enables automatic segmentation of PVSs on T2-weighted 3T MR images.
Method: Twenty patients with Parkinson’s disease (age range, 42–79 years) participated in this study. Specifically, we introduced a multi-scale supervised dense nested attention network designed to segment the PVSs. This model fosters progressive interactions between high-level and low-level features. Simultaneously, it utilizes multi-scale foreground content for deep supervision, aiding in refining segmentation results at various levels.
Result: Our method achieved the best segmentation results compared with the four other deep-learning-based methods, achieving a dice similarity coefficient (DSC) of 0.702. The results of the visual count of the PVSs in our model correlated extremely well with the expert scoring results on the T2-weighted images (basal ganglia: rs = 0.845, P < 0.001; rs = 0.868, P < 0.001; centrum semiovale: rs = 0.845, P < 0.001; rs = 0.823, P < 0.001 for raters 1 and 2, respectively). Experimental results show that the proposed method performs well in the segmentation of PVSs.
Conclusion: The proposed method can accurately segment PVSs; it will facilitate practical clinical applications and is expected to replace the method of visual counting directly on T1-weighted images or T2-weighted images.
1 Introduction
Perivascular spaces (PVSs), also known as Virchow–Robin spaces, are fluid-filled spaces that surround small blood vessels in the brain and traverse through the brain substance (Ballerini et al., 2018). PVSs are primarily distributed in areas such as the white matter, basal ganglia (BG), and brainstem (Salzman et al., 2005). They play a role in cerebrospinal fluid (CSF) circulation and clearance of metabolic waste from the brain (Iliff et al., 2013; Rangroo Thrane et al., 2013; Yang et al., 2013). In magnetic resonance imaging (MRI), PVSs tend to appear as longer linear structures when the direction of penetration is parallel to the scanning plane and as small dot-like structures when they are perpendicular to the scanning plane; PVSs usually have a signal intensity similar to that of CSF (Figure 1). It is important to note that the size and number of PVSs are related to aging (Francis et al., 2019) and various diseases, including Alzheimer’s disease (AD) (Hansen et al., 2015; Boespflug et al., 2018), Parkinson’s disease (PD) (Shibata et al., 2019; Shen et al., 2021), multiple sclerosis (MS) (Wuerfel et al., 2008; Conforti et al., 2014; George et al., 2021), and small vessel disease (SVD) (Doubal et al., 2010; Zhu et al., 2010). PVS changes may reflect the pathogenesis of certain neurodegenerative diseases (Charidimou et al., 2013). Therefore, accurate identification and quantification of PVSs are crucial when studying these diseases.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1457405/fnagi-16-1457405-HTML/image_m/fnagi-16-1457405-g001.jpg)
Figure 1. Examples of PVSs on magnetic resonance images. T1WI represents T1-weighted images, and T2WI represents T2-weighted images; (A,B) Indicate that the PVSs appear as lines and dots, respectively. PVSs, perivascular spaces.
Visual rating scales remain the gold standard for assessing PVSs burden, and numerous scales have been developed for this purpose (Patankar et al., 2005; Chen et al., 2011; Rowley, 2013; Wardlaw et al., 2013; Laveskog et al., 2018). Visual scoring is a simple and intuitive method; however, it only provides information about the presence and degree of PVSs. It does not offer more detailed quantitative information regarding parameters such as volume, morphological features, and distribution of features in the brain (Bouvy et al., 2016). Additionally, when using visual rating scales, scans are compared based on the category and not on the number of PVS. If the number of PVS increases over time but remains within the same category, visual scoring will not detect changes over time (Moses et al., 2023). Consequently, fully automated methods should be developed to compensate for and improve the shortcomings of the visual scoring method for rating PVSs.
In recent years, some researchers have proposed methods for automatic PVSs segmentation based on deep learning. However, most studies focused on 7T MRI (Park et al., 2016; Hou et al., 2017; Zhang et al., 2017; Lian et al., 2018; Spijkerman et al., 2022), with a relatively limited number of studies conducted on 3T MRI (Boutinaud et al., 2021; Rashid et al., 2023). Additionally, most of the study participants were healthy young people (Park et al., 2016; Hou et al., 2017; Zhang et al., 2017; Lian et al., 2018; Boutinaud et al., 2021; Spijkerman et al., 2022). Moreover, 7T MRI is rarely used in routine clinical practice. Therefore, we aimed to develop an automatic segmentation model for PVSs in T2-weighted 3T MRI based on patients with PD.
2 Materials and methods
2.1 Data
We prospectively collected cranial MR images of 20 patients with PD. The average age of the patients was 63.0 ± 10.3 years (mean ± SD, range: [42–79], median = 64 years), with a sex ratio of 1:1. As recommended by the 2023 STRIVE-2 guidelines (Duering et al., 2023), for all participants, 3D T1-weighted and high-resolution T2-weighted images were collected using a 3T MR scanner equipped with a 64-channel head coil (Magnetom Prisma; Siemens Healthineers, Germany). Detailed information on the MRI acquisition parameters is provided in Table 1.
2.2 Production of ground truth for PVSs
To develop and assess the proposed machine learning-based method for extracting PVSs, we initially created ground truth PVSs masks through a comprehensive process involving visual inspection. An expert manually segmented the PVSs of the 20 subjects using ITK-SNAP software (version 3.8.0). These segments were then refined by two experienced neuroradiologists, each with over 5 years of clinical practice experience. The experts meticulously reviewed the preliminary ground truth that was formulated by combining the T2-weighted images. Disagreements were resolved through discussion to ensure the precision of the ground truth for effective model training.
2.3 Visual counting of PVSs burden
Three independent raters counted the PVSs on T2-weighted images and the model segmentation results for all participants. Among the raters, one was a senior neurologist and the other two were senior neuroradiologists, each with over 20 years of clinical practice experience. Both the senior neurologist and neuroradiologist were provided the T1-weighted and T2-weighted images for all participants. They independently counted the PVSs in selected slices on T2-weighted images. Another senior neuroradiologist counted the PVSs in the same slices as those in the model segmentation results. Slice selection: for the BG, the slice shows the anterior commissure; for the centrum semiovale CSO, the slice was 1 cm above the uppermost part of the lateral ventricles. The number of PVSs in these slices correlates well with the number of PVSs in the entire volume of the region (Adams et al., 2013).
2.4 Network architecture
We assembled a series of U-shaped sub-networks to create a densely nested structure. This design caters to the varying optimal receptive fields required for targets of different sizes, with the depth of each U-shaped subnetwork being ideal for specific target sizes. This concept was further enhanced by placing multiple nodes along the encoder-to-decoder pathways. These nodes were interconnected and formed nested networks. As depicted in Figure 2, each node processes the features of its own layer as well as those of adjacent layers, enabling a thorough multilayered fusion of features. Consequently, this architecture effectively preserves the representation of smaller targets in the deeper layers, leading to improved outcomes, and we propose the incorporation of multi-scale Highlighting Foregrounds (HFs) to enhance deep supervision within the densely nested U-Net framework. Based on the above characteristics, we named the model MfNS_De. Our modified network, depicted in Figure 2, is based on an encoder-decoder structure. In both the training and testing phases, the network processed individual slices of T2-weighted (T2w) modalities as inputs. We chose axial slices for the 2D input because of the discontinuities in the third dimension, a characteristic inherent to our images.
In our network decoder, we integrated modified label images processed using the multi-scale HF approach at intermediate layers (as shown in Figure 2, right). These intermediate output convolutional layers transform feature maps from different decoder layers into multi-scale segmentation probability maps. Concurrently, multi-scale HFs were applied to down sample the label or ground truth images to various scales (Figure 2, right). The foreground/background label images created using multi-scale HFs were then utilized to generate losses by comparison with their corresponding outputs. For the loss function, we used a soft dice score.
2.5 Training details
We have implemented the proposed network on an NVIDIA Tesla V100-SXM2 computer using PyTorch. To compare the model performance between our network and others, we employed five-fold cross-validation. All subjects were initially categorized into five groups, with the age distribution being fairly uniform across these groups. During each cycle of model training, one group was designated as the test set, while another one was randomly chosen as the validation set, and the remaining three groups were used as the training sets. The composition of the training, validation, and test sets remained the same in all comparison networks. The hyper-parameters of the networks are set as follows: mini-batch = 4, optimiser = Adagrad, learn rate = 0.05, epoch = 1000. Furthermore, the version of nnU-Net at the time is nnU-Net v2, which is configured as 3D full-resolution, to ensure a fair comparison. The model was trained for 1,000 epochs, and good convergence was achieved after 50 epochs. Therefore, we selected the model at epoch 50 as the well-trained model for the test data. When calculating the performance metrics (e.g., DSC, SEN, and PPV), we used the same computational method for each model to ensure the consistency and comparability of the results.
2.6 Evaluation metrics
Segmentation performance was evaluated using the dice similarity coefficient (DSC), sensitivity (SEN), and positive predictive value (PPV), as defined below:
Where TP, FP, and FN denote the true positive, false positive, and false negative, respectively. DSC reflects the overall segmentation performance, SEN indicates the capability of detecting the PVSs voxels, and PPV represents the capability of discarding the confounding background voxels. The correlation of PVSs count between the model segmentation results and the original T2-weighted images was obtained by calculating the Spearman’s correlation coefficient.
3 Results
3.1 Segmentation performance
Table 2 and Figure 3 present the segmentation results obtained using our MfNS_De method and the four other comparison methods (IAANet, TriSegNet, U-Net, and nnU-Net). We could obtain the following observations: first, our method achieved approximately 23%, 11%, 4%, and 2% average DSC enhancement, respectively. Second, the proposed MfNS_De outperformed the original U-Net mainly because of the use of two key modules in the proposed method: the densely nested U-Net structure and the multi-scale feature learning strategy supervised by multi-scale down sampling of the label images.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1457405/fnagi-16-1457405-HTML/image_m/fnagi-16-1457405-t002.jpg)
Table 2. Segmentation results obtained by five different models: the best scores are highlighted as boldface.
Figure 4 shows the segmentation results for several typical PVSs. Clearly, IAANet and TriSegNet could not detect many PVSs of normal signals (Figure 4A). Additionally, IAANet, TriSegNet, and U-Net faced difficulty in distinguishing some of the fine sulci from the PVSs (Figures 4B, C), resulting in many sulci being incorrectly categorized as PVSs. Moreover, the IAANet was unable to distinguish some lacunar infarctions (Figure 4C). Similarly, TriSegNet and nnU-Net were prone to detecting some slightly higher-signal non-PVS voxels as PVSs (Figures 4C, D), leading to a large number of misreported voxels. Notably, MfNS_De effectively resolves the above problems, and its segmentation results are more consistent with the ground truth.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1457405/fnagi-16-1457405-HTML/image_m/fnagi-16-1457405-g004.jpg)
Figure 4. Illustration of typical PVSs segmentation by five different models. (A–D) Indicate different subjects, respectively. Red arrows indicate PVSs that are not successfully segmented; blue arrows indicate the identification of the sulcus as PVSs; red circle indicates lacunar infarction; and blue circles indicate the identification of some of these slightly higher signaling voxels as PVSs. PVSs, perivascular spaces.
Figure 5 shows a 3D view of different PVSs burdens, illustrating the distribution of whole-brain PVSs. In Figures 5B, C, IAANet and TriSegNet demonstrate a significantly higher whole-brain PVSs burden than the ground truth, and the PVSs morphology appears more irregular, which is consistent with the results shown in Figure 4. In contrast, the 3D view demonstrated by MfNS_De and nnU-Net closely resembles the ground truth. However, in Figure 5B, we can see that MfNS_De performs better than nnU-Net for detecting PVSs at the brain’s edges.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1457405/fnagi-16-1457405-HTML/image_m/fnagi-16-1457405-g005.jpg)
Figure 5. 3D results of the segmentation of large amounts of PVSs (A), medium amounts of PVSs (B), and small amounts of PVSs (C) using five different models, shown in axial view. PVSs, perivascular spaces.
3.2 Visual count comparison
We evaluated the ability of our MfNS_De model to detect PVSs by comparing them with expert visual counts. Figure 6 shows the scatter plots and Spearman’s correlation coefficients for PVSs counts in the BG and CSO. There was a very high correlation between counts by visual raters and our model’s detection of PVSs in the same section (BG:rs = 0.845, P < 0.001; rs = 0.868, P < 0.001; CSO:rs = 0.845, P < 0.001; rs = 0.823, P < 0.001 for raters 1 and 2, respectively), and this result almost reached an inter-rater correlation (rs = 0.920, P < 0.001; rs = 0.915, P < 0.001 for BG and CSO, respectively). Overall, the PVSs count from our model’s automatic segmentation was lower than that from T2-weighted MRI, which could be because the model ignored very small and low-signal PVSs. The data collected had a heavy PVSs burden, with counts reaching 132 (CSO) and 38 (BG) based on our model’s segmentation results. As shown in Figure 6, one can notice that even for scans with a heavy PVSs burden, the PVSs counts from the automatic segmentation of the model were satisfactory. If humans were to visually rate scans with many PVSs on the original T2-weighted images alone, it would be a very time-consuming and labor-intensive task, and the results might also be more variable owing to subjective judgment.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1457405/fnagi-16-1457405-HTML/image_m/fnagi-16-1457405-g006.jpg)
Figure 6. Scatter plots of visual counts of the PVSs between T2-weighted images and our MfNS_De segmentation results. Purple and blue denote centers of semiovales (CSO) and basal ganglia (BG), respectively; rs denotes Spearman’s correlation coefficient.
4 Discussion
An increasing number of studies have indicated that PVSs are associated with cerebrospinal fluid circulation, clearance of metabolic waste from the brain, and certain neurodegenerative diseases (Iliff et al., 2013; Rangroo Thrane et al., 2013; Yang et al., 2013; Francis et al., 2019). Accurate and convenient evaluation of PVSs has been widely discussed in recent years. In this study, we proposed an automatic PVSs segmentation method of the whole brain in T2-weighted 3T MRI based on deep learning and achieved excellent segmentation results. We compared the segmentation accuracy between the MfNS_De method and the other four methods (IAANet, TriSegNet, U-Net, and nnU-Net). Our method achieved about 23%, 11%, 4%, and 2% average DSC improvement, respectively.
nnU-Net utilizes the original U-Net structure to optimize segmentation results through pre-processing and post-processing. It achieved the best automatic segmentation performance on 33 of the 53 anatomical structures evaluated, demonstrating strong generalization characteristics. nnU-Net does not require expert knowledge or computational resources beyond standard network training, and no manual task-specific adaptation is necessary, which is why we chose nnU-Net for comparison. Our results showed that nnU-Net exhibited higher sensitivity than all the other models, including MfNS_De, but at the cost of a lower PPV. In other words, nnU-Net accurately detects PVS but also includes more false positive regions, which can be clinically harmful. Therefore, we still consider DSC to be a more reasonable metric. It is noteworthy that nnU-Net’s DSC of 0.68 does not reach MfNS_De’s value of 0.702, indicating that we outperform nnU-Net in the PVS segmentation task.
The nested design of our network architecture allows multi-scale feature learning, which is crucial in image segmentation tasks, where objects of interest vary in size and shape. This network can effectively capture features at different scales, thereby improving its accuracy and robustness in segmenting small objects. The densely nested U-Net architecture enhances feature fusion from different layers. By combining features across various depths of the network, the model can leverage both low-level texture and edge details, as well as high-level contextual information. This comprehensive feature integration is particularly effective in achieving precise and detailed segmentation results, making it a powerful tool for various image segmentation applications. Moreover, multi-scale HFs facilitate the detection and segmentation of small objects by focusing on features at various scales. This is particularly useful for capturing the nuances of smaller objects that may be lost on a single scale or using a conventional U-Net model.
As shown in Figures 4, 5, IAANet and TriSegNet either frequently detect many anomalous voxels as PVSs or miss many true PVSs. The U-Net cannot distinguish between a small number of blended boundaries. Similarly, nnU-Net tends to detect many slightly higher-signal non-PVS voxels as PVSs, which results in nnU-Net’s SEN scores being slightly higher than those in our method; however, its DSC and PPV scores were lower than those in our method owing to a higher false positive rate. Our method clearly outperforms these four comparison models. Specifically, we proposed the MfNS_De model, which introduces a novel approach to image segmentation that is particularly effective for small objects, such as PVS. Its unique architectural feature, namely densely nested layers, facilitates a comprehensive understanding of both high-level and low-level image details, which is crucial for accurately segmenting small objects. The integration of features from different network depths yields precise and detailed segmentation results. Additionally, deep supervision with multi-scale HFs introduces a novel approach for segmenting small objects, combining the strengths of multiple layers of models with the precision of multi-scale analysis. The integration of multi-scale HFs allows the effective capture of detailed features at various scales, which is crucial for the precise identification of smaller objects that are often missed when using traditional methods. This combination significantly reduces false positives and improves segmentation accuracy. Overall, the proposed MfNS_De model effectively improves the segmentation accuracy by integrating different network depth features as well as multi-scale salient foreground depth supervision methods. The segmentation results show that the proposed strategy improves the segmentation performance, and the MfNS_De method detects PVSs in the whole brain so that a 3D view of the PVSs can be obtained (Figure 5), which helps doctors visualize the morphology, number, and distribution of PVSs in the brain.
Numerous studies have investigated the automatic segmentation of PVSs (Hou et al., 2017; Zhang et al., 2017; Lian et al., 2018; Boutinaud et al., 2021; Spijkerman et al., 2022; Rashid et al., 2023). Lian et al. (2018) used a fully convolutional neural network (FCNN) machine learning approach on a dataset of T2-weighted magnetic MRI acquired using a 7 T scanner. The DSC of their method was 0.77, which is the highest result among the articles we retrieved thus far. Given that 7TMR imaging is rarely used in clinical practice, the segmentation of PVSs in 3T scanners is beneficial for studying various diseases. However, the number of current research on automatic PVSs segmentation in 3T scanners is limited (Boutinaud et al., 2021; Rashid et al., 2023); in the literature we have retrieved, the reported DSC has not exceeded 0.70, except for some specific PVS (PVS clusters and large PVSs) (Boutinaud et al., 2021). Most articles have not disclosed the dice score, a measure of overall segmentation performance and the most rigorous metric for evaluating the performance of PVSs segmentation algorithms (Pham et al., 2022). This also reinforces the clinical value of our 3T MRI model.
We found some literature that showcases several excellent methods for brain tumor segmentation. Liu et al. (2022) presented a glioma segmentation method based on adversarial learning. It fuses contrast-enhanced T1-weighted and Flair MRI images, using a semantic segmentation network as a discriminator to extract tumor-related information. Zhu et al. (2023) proposed a brain tumor segmentation method based on the fusion of deep semantics and edge information. It primarily utilizes an improved Swin Transformer for semantic segmentation, combines convolutional neural network for edge detection, and employs graph convolution for feature fusion. The method aims to enhance segmentation accuracy using multimodal MRI data by fully leveraging deep semantic features and edge features. Afterward, Zhu et al. (2024) proposed a new brain tumor segmentation method that consists of three modules, the Modality Information Extraction Module (MIE) for weighting different modality information, the Spatial Information Enhancement Module (SIE) for enhancing spatial information extraction through dilated convolutions, and the Boundary Shape Correction Module (BSC) for improving segmentation accuracy by selecting and constraining critical boundary points. I believe our upcoming PVS segmentation work can draw new ideas from these articles to build a more efficient PVS segmentation model.
Our study has some limitations. First, as shown in Figure 7, failed segmentation may occur in a few cases. The MfNS_De method may not detect PVSs with very low signals on T2-weighted images, leading to a lack of continuity in the segmentation results for long PVSs with uneven signal strengths along the extension direction. Moreover, a possible impact to model generalizability since the images were from PD patients. For future work, we need to consider a large sample to improve the model’s generalizability, including children, normal persons, and patients with other diseases. Finally, high-resolution T2 (slice thickness: 2 mm/slice gap: 0 mm) can reduce motion artifacts and display the anatomical structures of the BG more clearly (Rasouli et al., 2018; Vos et al., 2018). Therefore, we set the slice thickness of the T2-weighted images to 2mm, but this may also cause us to miss some unscanned PVSs.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1457405/fnagi-16-1457405-HTML/image_m/fnagi-16-1457405-g007.jpg)
Figure 7. Typical cases of MfNS_De method segmentation failure. Blue arrows indicate low-signal PVSs. PVSs, perivascular spaces.
5 Conclusion
In this study, we introduced a multi-scale supervised dense nested attention network designed to segment the PVS based on T2-weighted 3T MRI. Our method achieved the best segmentation results compared with the four other deep-learning-based methods, reaching a DSC of 0.702. The results of the visual count of PVSs in our model showed an extremely high correlation with those of experts on T2-weighted images (BG: rs = 0.845, P < 0.001; rs = 0.868, P < 0.001; CSO: rs = 0.845, P < 0.001; rs = 0.823, P < 0.001 for raters 1 and 2, respectively). We believe that this method will facilitate practical clinical applications and is expected to replace the method of direct visual counting directly on T1-weighted images or T2-weighted images.
Data availability statement
The original contributions presented in the study are included in the article/supplementary material, further inquiries can be directed to the corresponding authors.
Ethics statement
The studies involving humans were approved by the Shenzhen Second People’s Hospital’s medical ethics committee. The studies were conducted in accordance with the local legislation and institutional requirements. The participants provided their written informed consent to participate in this study.
Author contributions
DC: Data curation, Investigation, Methodology, Writing−original draft. MP: Methodology, Writing−original draft, Software. CL: Methodology, Writing−original draft, Formal analysis. WH: Formal analysis, Data curation, Writing−original draft. XG: Software, Visualization, Writing−original draft. JL: Formal analysis, Writing−original draft. RL: Formal analysis, Writing−original draft. ML: Conceptualization, Funding acquisition, Project administration, Writing−review and editing. JX: Conceptualization, Funding acquisition, Project administration, Writing−review and editing.
Funding
The author(s) declare financial support was received for the research, authorship, and/or publication of the article. This work was supported by the National Natural Science Foundation of China (grant number 82171913); the Shenzhen Science and Technology Program (grant number JCYJ20220818101816036); the Key Project of Clinical Research of Shenzhen Second People’s Hospital (grant number 20243357010); the Shenzhen Science and Technology Innovation Program (grant number JCYJ20220818102414031) and the “Young Innovation Team Program” of Shandong Provincial University (grant number 2022KJ250).
Acknowledgments
We would like to extend our gratitude to all study participants.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
References
Adams, H. H. H., Cavalieri, M., Verhaaren, B. F. J., Bos, D., Van Der Lugt, A., Enzinger, C., et al. (2013). Rating method for dilated Virchow-Robin spaces on magnetic resonance imaging. Stroke 44, 1732–1735. doi: 10.1161/STROKEAHA.111.000620
Ballerini, L., Lovreglio, R., Valdés Hernández, M. D. C., Ramirez, J., MacIntosh, B. J., Black, S. E., et al. (2018). Perivascular spaces segmentation in brain MRI using optimal 3D filtering. Sci. Rep. 8:2132. doi: 10.1038/s41598-018-19781-5
Boespflug, E. L., Simon, M. J., Leonard, E., Grafe, M., Woltjer, R., Silbert, L. C., et al. (2018). Targeted assessment of enlargement of the perivascular space in Alzheimer’s disease and vascular dementia subtypes implicates astroglial involvement specific to Alzheimer’s disease. J. Alzheimers Dis. 66, 1587–1597. doi: 10.3233/JAD-180367
Boutinaud, P., Tsuchida, A., Laurent, A., Adonias, F., Hanifehlou, Z., Nozais, V., et al. (2021). 3D segmentation of perivascular spaces on t1-weighted 3 tesla MR images with a convolutional Autoencoder and a U-Shaped neural network. Front. Neuroinform. 15:641600. doi: 10.3389/fninf.2021.641600
Bouvy, W. H., Zwanenburg, J. J. M., Reinink, R., Wisse, L. E. M., Luijten, P. R., Kappelle, L. J., et al. (2016). Perivascular spaces on 7 Tesla brain MRI are related to markers of small vessel disease but not to age or cardiovascular risk factors. J. Cereb. Blood Flow Metab. 36, 1708–1717. doi: 10.1177/0271678X16648970
Charidimou, A., Meegahage, R., Fox, Z., Peeters, A., Vandermeeren, Y., Laloux, P., et al. (2013). Enlarged perivascular spaces as a marker of underlying arteriopathy in intracerebral haemorrhage: A multicentre MRI cohort study. J. Neurol. Neurosurg. Psychiatry 84, 624–629. doi: 10.1136/jnnp-2012-304434
Chen, W., Song, X., and Zhang, Y., Alzheimer’s Disease, and Neuroimaging Initiative. (2011). Assessment of the Virchow-Robin spaces in Alzheimer disease, mild cognitive impairment, and normal aging, using high-field MR imaging. AJNR Am. J. Neuroradiol. 32, 1490–1495. doi: 10.3174/ajnr.A2541
Conforti, R., Cirillo, M., Saturnino, P. P., Gallo, A., Sacco, R., Negro, A., et al. (2014). Dilated virchow-robin spaces and multiple sclerosis: 3 T magnetic resonance study. Radiol. Med. 119, 408–414. doi: 10.1007/s11547-013-0357-9
Doubal, F. N., MacLullich, A. M. J., Ferguson, K. J., Dennis, M. S., and Wardlaw, J. M. (2010). Enlarged perivascular spaces on MRI are a feature of cerebral small vessel disease. Stroke 41, 450–454. doi: 10.1161/STROKEAHA.109.564914
Duering, M., Biessels, G. J., Brodtmann, A., Chen, C., Cordonnier, C., de Leeuw, F. E., et al. (2023). Neuroimaging standards for research into small vessel disease-advances since 2013. Lancet Neurol. 22, 602–618. doi: 10.1016/S1474-4422(23)00131-X
Francis, F., Ballerini, L., and Wardlaw, J. M. (2019). Perivascular spaces and their associations with risk factors, clinical disorders and neuroimaging features: A systematic review and meta-analysis. Int. J. Stroke 14, 359–371. doi: 10.1177/1747493019830321
George, I. C., Arrighi-Allisan, A., Delman, B. N., Balchandani, P., Horng, S., and Feldman, R. (2021). A novel method to measure venular perivascular spaces in patients with MS on 7T MRI. AJNR Am. J. Neuroradiol. 42, 1069–1072. doi: 10.3174/ajnr.A7144
Hansen, T. P., Cain, J., Thomas, O., and Jackson, A. (2015). Dilated perivascular spaces in the Basal Ganglia are a biomarker of small-vessel disease in a very elderly population with dementia. AJNR Am. J. Neuroradiol. 36, 893–898. doi: 10.3174/ajnr.A4237
Hou, Y., Park, S. H., Wang, Q., Zhang, J., Zong, X., Lin, W., et al. (2017). Enhancement of perivascular spaces in 7 T MR image using haar transform of non-local cubes and block-matching filtering. Sci. Rep. 7:8569. doi: 10.1038/s41598-017-09336-5
Iliff, J. J., Wang, M., Zeppenfeld, D. M., Venkataraman, A., Plog, B. A., Liao, Y., et al. (2013). Cerebral arterial pulsation drives paravascular CSF-interstitial fluid exchange in the murine brain. J. Neurosci. 33, 18190–18199. doi: 10.1523/JNEUROSCI.1592-13.2013
Laveskog, A., Wang, R., Bronge, L., Wahlund, L. O., and Qiu, C. (2018). Perivascular spaces in old age: Assessment, distribution, and correlation with white matter hyperintensities. AJNR Am. J. Neuroradiol. 39, 70–76. doi: 10.3174/ajnr.A5455
Lian, C., Zhang, J., Liu, M., Zong, X., Hung, S. C., Lin, W., et al. (2018). Multi-channel multi-scale fully convolutional network for 3D perivascular spaces segmentation in 7-T MR images. Med. Image Anal. 46, 106–117. doi: 10.1016/j.media.2018.02.009
Liu, Y., Shi, Y., Mu, F., Cheng, J., and Chen, X. (2022). Glioma segmentation-oriented multi-modal MR image fusion with adversarial learning. IEEE CAA J. Automat. Sin. 9, 1528–1531. doi: 10.1109/JAS.2022.105770
Moses, J., Sinclair, B., Law, M., O’Brien, T. J., and Vivash, L. (2023). Automated methods for detecting and quantitation of enlarged perivascular spaces on MRI. J. Magn. Reson. Imaging 57, 11–24. doi: 10.1002/jmri.28369
Park, S. H., Zong, X., Gao, Y., Lin, W., and Shen, D. (2016). Segmentation of perivascular spaces in 7-T MR image using auto-context model with orientation-normalised features. Neuroimage 134, 223–235. doi: 10.1016/j.neuroimage.2016.03.076
Patankar, T. F., Mitra, D., Varma, A., Snowden, J., Neary, D., and Jackson, A. (2005). Dilatation of the Virchow-Robin space is a sensitive indicator of cerebral microvascular disease: Study in elderly patients with dementia. AJNR Am. J. Neuroradiol. 26, 1512–1520.
Pham, W., Lynch, M., Spitz, G., O’Brien, T., Vivash, L., Sinclair, B., et al. (2022). A critical guide to the automated quantification of perivascular spaces in magnetic resonance imaging. Front. Neurosci. 16:1021311. doi: 10.3389/fnins.2022.1021311
Rangroo Thrane, V., Thrane, A. S., Plog, B. A., Thiyagarajan, M., Iliff, J. J., Deane, R., et al. (2013). Paravascular microcirculation facilitates rapid lipid transport and astrocyte signaling in the brain. Sci. Rep. 3:2582. doi: 10.1038/srep02582
Rashid, T., Liu, H., Ware, J. B., Li, K., Romero, J. R., Fadaee, E., et al. (2023). Deep learning based detection of enlarged perivascular spaces on brain MRI. Neuroimage Rep. 3:100162. doi: 10.1016/j.ynirp.2023.100162
Rasouli, J., Ramdhani, R., Panov, F. E., Dimov, A., Zhang, Y., Cho, C., et al. (2018). Utilization of quantitative susceptibility mapping for direct targeting of the subthalamic nucleus during deep brain stimulation surgery. Oper. Neurosurg. 14, 412–419. doi: 10.1093/on
Rowley, H. A. (2013). The alphabet of imaging in acute stroke: Does it spell improved selection and outcome? Stroke 44, S53–S54. doi: 10.1161/STROKEAHA.113.001939
Salzman, K. L., Osborn, A. G., House, P., Jinkins, J. R., Ditchfield, A., Cooper, J. A., et al. (2005). Giant tumefactive perivascular spaces. AJNR Am. J. Neuroradiol. 26, 298–305.
Shen, T., Yue, Y., Zhao, S., Xie, J., Chen, Y., Tian, J., et al. (2021). The role of brain perivascular space burden in early-stage Parkinson’s disease. NPJ Parkinsons Dis. 7:12. doi: 10.1038/s41531-021-00155-0
Shibata, K., Sugiura, M., Nishimura, Y., and Sakura, H. (2019). The effect of small vessel disease on motor and cognitive function in Parkinson’s disease. Clin. Neurol. Neurosurg. 182, 58–62. doi: 10.1016/j.clineuro.2019.04.029
Spijkerman, J. M., Zwanenburg, J. J. M., Bouvy, W. H., Geerlings, M. I., Biessels, G. J., Hendrikse, J., et al. (2022). Automatic quantification of perivascular spaces in T2-weighted images at 7 T MRI. Cereb. Circ. Cogn. Behav. 3:100142. doi: 10.1016/j.cccb.2022.100142
Vos, S. B., Micallef, C., Barkhof, F., Hill, A., Winston, G. P., Ourselin, S., et al. (2018). Evaluation of prospective motion correction of high-resolution 3D-T2-FLAIR acquisitions in epilepsy patients. J. Neuroradiol. 45, 368–373. doi: 10.1016/j.neurad.2018.02.007
Wardlaw, J. M., Smith, E. E., Biessels, G. J., Cordonnier, C., Fazekas, F., Frayne, R., et al. (2013). Neuroimaging standards for research into small vessel disease and its contribution to ageing and neurodegeneration. Lancet Neurol. 12, 822–838. doi: 10.1016/S1474-4422(13)70124-8
Wuerfel, J., Haertle, M., Waiczies, H., Tysiak, E., Bechmann, I., Wernecke, K. D., et al. (2008). Perivascular spaces–MRI marker of inflammatory activity in the brain? Brain 131, 2332–2340. doi: 10.1093/brain/awn171
Yang, L., Kress, B. T., Weber, H. J., Thiyagarajan, M., Wang, B., Deane, R., et al. (2013). Evaluating glymphatic pathway function utilising clinically relevant intrathecal infusion of CSF tracer. J. Transl. Med. 11:107. doi: 10.1186/1479-5876-11-107
Zhang, J., Gao, Y., Park, S. H., Zong, X., Lin, W., and Shen, D. (2017). Structured learning for 3-D perivascular spaces segmentation using vascular features. IEEE Trans. Bio Med. Eng. 64, 2803–2812. doi: 10.1109/TBME.2016.2638918
Zhu, Y. C., Tzourio, C., Soumaré, A., Mazoyer, B., Dufouil, C., and Chabriat, H. (2010). Severity of dilated Virchow-Robin spaces is associated with age, blood pressure, and MRI markers of small vessel disease: A population-based study. Stroke 41, 2483–2490. doi: 10.1161/STROKEAHA.110.591586
Zhu, Z. Q., He, X. Y., Qi, G. Q., Li, Y. Y., Cong, B. S., and Liu, Y. (2023). Brain tumor segmentation based on the fusion of deep semantics and edge information in multimodal MRI. Inf. Fusion 91, 376–387. doi: 10.1016/j.inffus.2022.10.022
Keywords: perivascular spaces, Virchow-Robin spaces, deep learning, multiscale supervised, dense nesting, 3T MR image
Citation: Cai D, Pan M, Liu C, He W, Ge X, Lin J, Li R, Liu M and Xia J (2024) Deep-learning-based segmentation of perivascular spaces on T2-Weighted 3T magnetic resonance images. Front. Aging Neurosci. 16:1457405. doi: 10.3389/fnagi.2024.1457405
Received: 30 June 2024; Accepted: 12 August 2024;
Published: 29 August 2024.
Edited by:
Guanghui Wang, Soochow University, ChinaReviewed by:
Weiwei Chen, Huazhong University of Science and Technology, ChinaRongchang Zhao, School of Computer Science, Central South University, China
Zhiqin Zhu, Chongqing University of Posts and Telecommunications, China
Copyright © 2024 Cai, Pan, Liu, He, Ge, Lin, Li, Liu and Xia. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Jun Xia, xiajun@email.szu.edu.cn; Mengting Liu, liumt55@mail.sysu.edu.cn
†These authors have contributed equally to this work