Automated Detection of Lupus White Matter Lesions in MRI

1. Introduction

Several brain diseases present abnormalities in the white matter tissue, usually denoted as white matter lesions (WML). Segmenting these WML is important to diagnose and better understand these diseases as well as monitoring its progression. However, performing this task manually is tedious and very time consuming. Hence, several works have been proposed to tackle automatically this lesion segmentation problem. For instance, various approaches have been presented in multiple sclerosis lesion segmentation (Van Leemput et al., 2001; Lladó et al., 2012a,b; Schmidt et al., 2012; Cabezas et al., 2014a,b; Guizard et al., 2015; Jain et al., 2015; Roura et al., 2015; Brosch et al., 2016), stroke (Mitsias et al., 2002), vascular dementia (Mohamed et al., 2001; Yamashita et al., 2008), and other diseases (Kruggel et al., 2008; Schwarz et al., 2009). Instead, few attempts have been done on semiautomatic or automatic segmentation of Lupus lesions (Appenzeller et al., 2008a,b; Petri et al., 2008; Scully et al., 2010), which have the particularity of being very small and focal WML, and they are few and isolated.

Magnetic resonance imaging (MRI) is the gold standard technique for studying the brain in lupus (Sarbu et al., 2015c). The neuroimaging findings are classified as small or large vessel disease, and inflammatory-type lesions (Sarbu et al., 2015b). Small vessel disease is represented by white-matter hyperintensities/lesions, recent small subcortical infarcts, lacunes, microbleeds, and brain atrophy (Wardlaw et al., 2013). WML are the most common findings of small vessel disease seen in lupus, and represent small T2-hyperintensities following the distribution of the white matter (periventricular, deep, subcortical), and including also the white matter at the basal ganglia, and cerebellum (Sarbu et al., 2015b,c). During the last years, WML have been shown to function as an independent predictor for the neurolupus activity and injury, and quantitative methods are increasingly proposed for the quantification and follow-up of the WML in neurolupus. As stated in Sarbu et al. (2015b), the stratification by the number of lesions is also important for the diagnosis of neuropsychiatric lupus. In their work, the authors introduced an stratification in three different groups: (1) low lesion burden (< 5 lesions); (2) medium lesion burden (between 5 – 25 lesions); and (3) high lesion burden (> 25 lesions).

To deal with the lupus WML segmentation problem, previous approaches, such as the automated one of Scully et al. (2010), have used a supervised strategy. In their work, local morphometric features extracted from multiple sequences, including T1w, T2w, and fluid-attenuated inversion recovery (FLAIR) images, were used to train a supervised classifier that takes advantage of a different subset of the features to segment lesion voxels. With a different viewpoint, in our work we present an unsupervised approach to automatically segment WML in Lupus patients by using only T1w and FLAIR images. This work can be seen as an extension of the tool recently presented by Roura et al. (2015), in which the focus was the segmentation of multiple sclerosis lesions. The whole pipeline can be considered as a two step process: pre-processing and WML segmentation. The first step is focused on the image enhancement by performing different intensity corrections on the brain and co-aligning all the image modalities. The second one, performs the lesion segmentation by detecting outliers to the normal apparent gray matter brain tissue on the FLAIR image as was previously done by Souplet et al. (2008) and Roura et al. (2015). Given the specific properties of the Lupus WML, we introduce a set of post-processing steps to reduce possible false positive (FP) detections which are based on lesion size, lesion tissue neighborhood, and lesion location. The last one aims to eliminate the FP detections usually found in the posterior fossa due to frequent scanner artifacts, yet this is an uncommon location for WML in neurolupus (up to 7% of patients vs. 40–60% in frontal lobes) (Sarbu et al., 2015b). We introduced this constraint in the segmentation by using an in-house atlas created with the unbiased template creation algorithm proposed by Fonov et al. (2011), which was then segmented into 12 brain structures including the posterior fossa using the Computational Morphometry Toolkit software¹.

The evaluation of the Lupus WML segmentation has been done on a dataset of 20 patients comparing quantitatively the results obtained by our tool with the ones performed manually by an expert neuroradiologist. This ground truth (GT) has been used to compute quantitative measures in terms of detection, such as True Positive Rate (TPR) and Positive Prediction Value (PPV), and in terms of segmentation accuracy by using the Dice Similarity Coefficient (DSC). Both detection and segmentation results show the ability of the approach to automatically detect and segment focal WML in Lupus patients. The code of our approach is publicly available as a Statistical Parametric Mapping (SPM8/12) toolbox extension with a simple parameter configuration².

2. Materials and Methods

2.1. Data

This study included 20 Lupus patients. The brain MRIs were performed between 2014 and 2015 at Hospital Clínic, University of Barcelona, Spain, the main national referral institution for lupus. All scans were performed at three Tesla Siemens MAGNETOM TIM Trio scanner, using a 32-channel head coil, with the same protocol including 3D T1 and 3D FLAIR, with a voxel size = 1 × 1 × 1mm³. The lesions were semiautomatically annotated on FLAIR images by an expert neuroradiologist. They present a lesion volume mean and range (min-max) per patient of 0.217 [11−1459] mm³. This study was carried out in accordance with the ethical recommendations of the Hospital Clínic committee (IDIBAPS, Barcelona), with written informed consent from all subjects.

2.2. MRI Pre-processing

To deal with the Lupus WML segmentation, several pre-processing steps (see Figure 1) are required to optimize the overall performance, as seen in previous works (Schmidt et al., 2012; Cabezas et al., 2014a,b; Valverde et al., 2014; Guizard et al., 2015; Roura et al., 2015). Since, our aim is to provide a publicly available Lupus segmentation tool as an extension of the SPM8/12 all the required steps are performed within the Matlab environment.

FIGURE 1

Figure 1. Scheme of the full pipeline. T1w and FLAIR images are the original subject images, which require a co-registration. Once co-registered they undergone a separate pre-processing, T1w by SPM8/12 to obtain the tissue segmentation and the brain mask, while FLAIR is denoised and unbiased by separate methods. The atlas and its structures also belong to the SPM space, so they can be brought to the subject space as the tissue segmentation and the brain mask. The corrected FLAIR image, the tissue segmentation and the posterior fossa mask are the inputs of the WML segmentation tool.

The first step of the pre-processing consists in the intra-subject registration. For this, we follow a similar procedure as the one used in Roura et al. (2015). In this case, we register FLAIR to the T1w image, where the target space used [corresponding to the Montreal Neurological Institute (MNI) (Mazziotta et al., 1995)] as well as the co-registration software are provided by the SPM toolbox.

One of the most common pre-processing step is the skull stripping process (Shattuck et al., 2001; Smith, 2002; Park and Lee, 2009; Roura et al., 2014), which we incorporate into our pipeline using the SPM tissue segmentation algorithm (Ashburner and Friston, 2005), avoiding therefore the use of external libraries such as BET (Smith, 2002) or BSE (Shattuck et al., 2001). Given that this process provides the probability map of the three main brain tissues [white matter (WM), gray matter (GM), and cerebrospinal fluid (CSF)], we created a maximum probability map to determine the three tissue masks and then thresholded the result at 0.5 to determine directly the brain mask. A similar procedure was also applied in previous segmentation works (Boesen et al., 2004; Roura et al., 2014). Notice that a threshold of 0.5 could lead to some holes in the corresponding masks due to the signal changes induced by the presence of WML. While partial volume effects will have only two tissue classes competing (WM/GM or GM/CSF), in the presence of lesions, areas such as periventricular, are prone to misclassifications among the three main tissue classes. In this case, a morphological operation should be applied in order to fill the inside holes of the brain mask. However, this is very unlikely and we did not experience this when using the SPM12 platform. This process is performed on the T1w image, although the brain mask is then applied on the FLAIR image where the rest of the pre-processing steps are carried out, since intensity corrections of the T1w image are handled by the SPM tissue segmentation process itself.

It is well-known that MRI images obtained directly from the scanner present noise and undesired artifacts (movement, high signal value, blood, flow artifacts, etc) (Lemieux et al., 1998). These abnormalities such as inhomogeneities in the magnetic field (Sha and Sutton, 2001) must be attenuated by applying post-scanning processes. We first apply the anisotropic diffusion filter of Perona and Malik (1990) in order to enhance the image by smoothing its histogram with the 3D Matlab implementation³ of this algorithm. Given the reduced size of the lesions, we have carefully run this method over all the patients with a restrictive parameter configuration (1 iteration, K = 50, and high contrast edges), reducing the iterations and gradient modulus, and focusing on contrast instead of region size.

To correct the bias field we used the Matlab method proposed by Larsen et al. (2014), which is based on an expectation maximization model (EM) that relies on the same generative models and bias field estimation computations of the well-known non-parametric, non-uniform intensity normalization (N3) method (Sled et al., 1998). This approach requires to mask out the low intensity voxels, thus the brain mask obtained from T1w image is used when correcting the FLAIR image.

2.3. Lupus Lesion Segmentation

Lupus lesions, similarly to other WML such as multiple sclerosis lesions, are characterized by being hyperintense regions in the FLAIR images. Due to the fact that the GM is the highest intensity tissue in this image modality, we used its histogram distribution to identify the hyperintense outliers. In order to obtain the GM distribution, we used the same SPM tissue segmentation (Ashburner and Friston, 2005) applied in the skull stripping process. At this point the lesion detection can be performed as a thresholding process, commonly computed by μ+ασ, where the standard deviation (σ) is determined using the full width at half maximum (FWHM) of the main peak (μ). We can then adjust the number of detected candidate lesions via the α parameter, observing a good trade off when setting this parameter to 2.5, assuming more than 98% of the histogram belonging to GM. Afterwards, we apply a set of post-processing steps to remove FP lesions that remained after thresholding the FLAIR image: (1) Lesion size: we constraint the minimum size of the lupus lesion to be 3 mm³. Therefore, we eliminate hyperintense voxels or a group of voxels smaller than this size. (2) Lesion tissue neighborhood (λ): because the lupus lesions should appear in the WM, the surrounding voxels must strictly belong to WM. Therefore, we introduce a parameter to limit the proportion of the WM over GM and CSF in the lesion neighborhood. We will see in Section 3 that the best trade off was obtained when using λ = 0.7. Looking at Figure 2, one can see how the neighbors of the two higher hyperintense regions marked in green in Figure 2B, all belong to WM in the tissue segmentation, while other candidate regions seen in the centre (marked in red) are not considered lesions because the neighbors voxels belong to GM. This neighborhood operation is applied in 3D. Figures 2C,D shows the original image and the tissue segmentation result of two slices forward, where the candidates marked in red are attached to GM and therefore eliminated with the neighborhood constraint. (3) Lesion location: since Lupus lesions are rarely present in the posterior fossa (Sarbu et al., 2015b,c), and this particular area is highly prone to present hyperintense artifacts, we have decided to exclude this region when looking for possible lesion candidates. This is done automatically by registering an atlas with the corresponding structure segmentation to the T1w image. In particular, we use an in-house 3T template created over healthy subjects using the unbiased template creation approach proposed by Fonov et al. (2011). This procedure, as stated by the authors, converges after 20 iterations, meaning that 20 non-rigid registrations must be performed for each subject of the population. The nonlinear registration process relies on the Automatic Nonlinear Image Matching and Anatomical Labeling (ANIMAL) of Collins et al. (1995). In order to obtain the structure segmentation of the healthy template, we have re-arranged the 83 labels of the T1w atlas from Hammers et al. (2003) into 12 regions⁴. Subsequently, our template was segmented into these 12 regions using the Computational Morphometry Toolkit (CMTK)⁵. To register the in-house template T1w image to each patient we used the SPM registration module, similarly to the intra-subject registration process. Finally, using the deformation field obtained by the non-rigid registration, we are able to bring the structure corresponding to the posterior fossa to each of the patient's space and therefore remove FP in this area caused by artifacts. A summary of the full pipeline is illustrated in Figure 1.

FIGURE 2

Figure 2. Neighboring rule. FLAIR 2D axial slice (A) showing 2 lesions (in green) of 4 and 5 voxels, respectively (B), both completely surrounded by white matter. Original image and tissue segmentation result of two slices forward are shown in (C,D). The bottom row shows four zooms of the original image, ground truth (green), candidates regions (red), and final lesion segmentation (blue).

3. Results

We have quantitatively analyzed the obtained results evaluating three different measures, TPR and PPV for lesion detection, sensitivity and precision, respectively, and DSC in terms of segmentation accuracy. The evaluation of our approach has been carried out through a two-fold cross-validation method where 10 random samples have been chosen for the first set while the rest have been used in the second set. For each training step we have exhaustively assessed the parameter configuration by computing a trade-off between DSC, TPR, and PPV using the traditional F-score measure (Zhang et al., 2015), which gives the harmonic mean of the three measures abovementioned as follows:

$\begin{array}{ll}F\mathrm{\text{-}}score=3\times \frac{DSC\times TPR\times PPV}{DSC+TPR+PPV}& (1)\end{array}$

A logical range of parameters has been tested for both α and λ (x and y axes, respectively in Figure 3), on each training set of the 10 patients randomly chosen within the two-fold cross-validation. The best configuration was α = 2.5 and λ = 0.70 for both sets, therefore these parameters were used in all the 20 patients of the dataset obtaining the results presented here.

FIGURE 3

Figure 3. Joint evaluation of both α and λ parameters (x and y axes, respectively) for each training set within the two-fold cross-validation. Each position in the map represents the mean F-score for a specific tissue neighborhood ratio λ and the initial candidate lesions adjustment α. Reddish colors show higher mean F-score values.

Figure 4 shows the obtained results per patient. As done in Sarbu et al. (2015b), we have stratified the population according to three different groups depending on the number of GT lesions per patient: (1) low lesion burden (< 5 lesions); (2) medium lesion burden (between 5 – 25 lesions); (3) high lesion burden (>25 lesions). Notice that this stratification by number of lesions may help on the diagnosis of neuropsychiatric lupus. Regarding the obtained results, both group and total averages of all the measures are over 50%, specially highlighting the group with more than 25 lesions, where we obtain a TPR = 0.81±0.14, a PPV = 0.96, and a DSC = 0.95±0.1. Notice, however, that this particular group only contains two samples and therefore we cannot extract significant conclusions. When considering the whole dataset, these values are: TPR = 0.62±0.19, PPV = 0.80±0.25, and DSC = 0.72±0.22; this decrease is due to the lower performance obtained in the first group of the stratification, where a small error represents a big percentage in the total measure. These cases with <5 lesions have also a small lesion volume which as stated in previous studies (Roura et al., 2015) worsen the performance of the automated segmentation methods.

FIGURE 4

Figure 4. Bar plots of each patient representing the DSC, TPR, and PPV values. The population is stratified in three groups depending on the GT number of lesions, from left to right: < 5; [5−25]; >25 lesions per patient.

To better understand the results, we also show two correlation plots, one with the number of lesions and one with the lesion volume (see Figure 5). In order to evaluate these correlations we first tested the normality of the data using the Lilliefors test based on the Kolmogorov-Smirnov normality test (Lilliefors, 1967). Given that the distributions did not follow a normal distribution we applied the Spearman's linear correlation coefficient. We have fitted a linear polynomial curve and showing also the expected fit, which is basically the ideal correlation. Looking at the number of lesions correlation, the model lies under the expected fit meaning that the approach underestimates the number of lesions. However, we obtained a very high linear correlation coefficient (r = 0.87) with a p ≪ 0.01, meaning that the coefficient is significantly different from zero, and therefore the null hypothesis of no correlation is not satisfied. The whole dataset can be linearly explained because all the samples follow the same trend, except one outlier which also has a correct stratification. Besides, one can see how the stratification results fit for most of the patients with the expected groups, except for two cases which are close to the group limit. Regarding the lesion volume correlation, the model fitting shows also a very good correlation, with only one sample out of the confidence level. The model coincides almost perfectly with the expected fit and Spearman's coefficient is also high with r = 0.90 with a p ≪ 0.01. Notice that this high fitting illustrates that the FP and FN are not significant compared to the TP in terms of total affected tissue volume, where the FNR (0.19±0.18) in terms of volume was relatively low with respect to the TPR (0.81±0.18), i.e., the lesions missed by the approach presented a very small volume and therefore the overall correlation per patient was high. To better understand the behavior on the low lesion volume per patient, we also show the Bland-Altman plot (Figure 6), where one can better perceive the good correlation of the scatter points skewed to the left in the first volume correlation plot. For a better visualization, we have also removed the two scatter points over 500 mm³ and also the outlier, which in fact shows a difference volume of around 400 mm³ while the other two scatter points lies close to the mean volume difference (10 mm³). Even though the mean difference volume was −8.9 mm³, indicating an overall oversegmentation trend of our tool, this is mainly due to the outlier case in which we have a very big overestimation of the total volume. In fact, when removing this outlier from the analysis the mean difference volume value changes to be 12.1 mm³ (see doted blue line in Figure 6), which indicates a general trend of our tool to undersegment the lesion volume. This reflects also the trend with the volume correlation plot shown in Figure 5. Furthermore, we have also analyzed the behavior individually for the 3 different groups depending on the number of lesions per patient: (1) < 5 lesions: mean volume difference of −7 mm³ (oversegmentation) with a mean number of lesions per patient of 2.5; (2) between 5 – 25 lesions (excluding the outlier): mean volume difference of 28.5 mm³ (undersegmentation) with a mean number of lesions per patient of 12.5; (3) >25 lesions: mean volume difference of 14.5 mm³ (undersegmentation) with a mean number of lesions per patient = 33. Therefore, from this analysis, we can point out that our tool tends to oversegment the cases with a very small lesion volume and number of lesions, in fact the extreme cases in which there is more confusion with possible artifacts and partial volume effects, while is more conservative (undersegmentation) with the rest of cases in which there are larger number of lesions and volume.

FIGURE 5

Figure 5. Correlation with number of lesions (stratified by the three groups) on the left and lesion volume, in terms of voxels, on the right.

FIGURE 6

Figure 6. Bland-Altman plot of the lesion volume for those cases below 500 mm³ of average volume. There is also one outlier not represented with 400 mm³ of difference volume and 370 mm³ of average volume. Besides, two scatter points over 500 mm³ of average volume have also been removed, although the difference volume was 10 mm³ for both. The* represents the mean difference volume when removing the outlier from the analysis.

Some samples of qualitative results are shown in Figure 7, where we compare the results of our automated tool with the GT annotations. We have chosen different samples to illustrate the performance in patients with different lesion load. Notice that the total lesion volume is very small in all of them, but the automatic detection provides a good performance in terms of TP while having a reduced number of FP and FN. When illustrating the whole 3D volume in the figure, those FP and FN are inappreciable because they are smaller than (10 mm³). However, we show some FP and FN examples on the 2D slices for the second and third group, zooming also into these regions in the first group.

FIGURE 7

Figure 7. Qualitative results of the approach. First row of each patient shows the original FLAIR image and second row shows the automatic segmentation (green, TP; red, FP; and yellow, FN).

The three cases shown in Figure 7 are the most representative for the FNs found within the juxtacortical area. As explained in the methodology, contextual information, and lesion size are the restrictions to filter the potential candidates, and both have to be accomplished. Given the small size of these juxtacortical lesions, a high initial threshold that only considers the focus of the lesions may fail at any of both requirements if they are low intense with respect to the rest. Although it can be solved increasing the α parameter, this could lead to higher FP detections in other regions.

In order to compare our approach with the state-of-the-art (Scully et al., 2010), we have also analyzed the performance of our approach in terms of sensitivity and specificity for some specific points of the receiver operating characteristic (ROC) curves, given that the combination of two parameters within a ROC curve is not trivial. In the work of (Scully et al., 2010), the authors reported a sensitivity of 94.3% and a specificity of 93.9% in the hold-out set at the lowest threshold. On the other hand, they obtained 100% specificity and 2.6% sensitivity when the highest threshold shown in the ROC curve was used. At the 80% of sensitivity (the midpoint of the curve), the reported specificity was around 99.5%, and at 90% sensitivity was close to 98% specificity. In our experimental results, we have increased the specificity maintaining similar sensitivity values. In particular, we observed 100% specificity at 3% of sensitivity using our approach, while at 80% of sensitivity the specificity was 99.99%. Besides, at 90% sensitivity our specificity was still higher than 99%. Notice that a specificity of 99% might still have large FP detections, since the brain volume is much larger than lesions volume.

4. Discussion

WML are the most common radiological finding in neurolupus. They are non-specific findings, being frequently observed in older age groups, migraine, chronic diseases, heart diseases, diabetes, high-blood pressure, dyslipidemia, and other vascular risk factors, although they are also present in asymptomatic subjects without known diseases (Sarbu et al., 2015c). However, WML are found in 40–60% of neurolupus patients, even at the onset of the disease, and many previous reports showed a higher frequency of WML in neurolupus when compared with lupus without neurolupus and general population (Castellino et al., 2008; Wardlaw et al., 2013; Sarbu et al., 2015b,c).

The pathogenesis of WM hyperintensity is attributed to chronic small vessel disease, which is supported by a study with radiologic-pathologic correlation in patients with neurolupus (Sibbitt et al., 2010). The underlying mechanisms for small vessel disease in neurolupus are not well understood, although multiple factors are incriminated, including accelerated atherosclerosis, direct immune mediated alterations, microembolisms, intimal hyperplasia, erythrocytes extravasation, fibrin thrombi, and coagulopathy (Joseph and Scolding, 2009; Sarbu et al., 2015a,c).

In neurolupus, WML involve preferentially the frontal and parietal regions, different from primary autoimmune demyelinating diseases such as multiple sclerosis. WML were repeatedly correlated with lupus duration, cognitive dysfunction, cerebrovascular syndrome, seizures, antiphospholipid antibody and low complement (C3, C4, CH50) levels (Ainiala et al., 2005; Appenzeller et al., 2008b; Toledano et al., 2013). A quantitative WML analysis in lupus patients demonstrated that age, duration of neuropsychiatric manifestations and total corticosteroid dosage were independent predictors for WML (Appenzeller et al., 2008b). Importantly, there was demonstrated a positive association between the lesion burden and the score of lupus activity (Systemic Lupus Erythematosus Disease Activity Index-SLEDAI) and injury (Systemic Lupus International Collaborating Clinic-SLICC). This means that WML are an independent predictor for lupus activity and injury, and suggests that the quantification of WML (either by number or, maybe better, by volumetric methods) and their follow-up, could be used for monitoring the disease progression and response to therapy (Appenzeller et al., 2008a,b; Sarbu et al., 2015c).

We have proposed in this work an automated tool which presents a good correlation in both number of lesions and lesion volume, as seen in Figure 5. Even though the obtained results tend to underestimate the lesion detections, the number of lesions detected have shown a good correlation with the stratified population into three groups. Notice that the FN rate has a weak influence on the final lesion volume, since Lupus WML are small focal lesions, characteristic of this particular disease. Besides, in order to reduce this FP detection rate introduced by the artifacts of the scanner, we have set a high lesion neighborhood restriction to belong to WM.

The parameter configuration has been set up with an exhaustive analysis over both α and lesion tissue neighborhood parameters, testing values from 1 to 3 each 0.1 and from 0 to 1 each 0.05, respectively. The analysis showed that the optimal configuration was with the α of 2.5 and tissue neighborhood ratio of 0.7. We want to remark that this will be the default configuration of the tool, however, other configurations provide very similar results. After performing this evaluation, we observed how the approach was fairly steady for a wide range of parameters. Moreover, as shown in the work of Roura et al. (2015), the use of the same parameters in a similar intensity thresholding strategy, provided also promising results when tested with more than 100 multiple sclerosis patient scans acquired from different machines (3T and 1.5T).

Even though this study has been evaluated with a dataset of 20 patients, we observed promising results in both lesion detection and segmentation, highly comparable to the state of the art approach of Scully et al. (2010). In their work they used T1w, T2w, and FLAIR sequences in a total of 27 patients, 10 of which were used for training and 17 for testing. Regarding the pre-processing we have included the same steps (co-registration, brain extraction and bias correction) being able to avoid the intensity standardization required in supervised approaches. We have based our approach on three basic characteristics (lesion size, tissue neighborhood of candidate regions, and brain location) while they used 49 morphological features (tissue types and its distance map for each tissue, left-right flipped images to enhance the lesion location, spatial location (x, y, z) to inform about global location, and neighborhoods means and medians). However, only a subset of those features was used at each level of the segmentation by the support vector machines classifier. The comparison at specific points of the ROC curve for both methods showed that our approach obtained similar or even better results. Although we have used a different dataset to evaluate the performance our approach, we have shown an increase in the specificity maintaining similar sensitivity values.

One of the limitations of these automatic WML detection and segmentation approaches are the FN found within the juxtacortical lesions. Even though in our experiments the ratio found was very low, some illustrative examples of this problem are shown in Figure 7. One could address this issue by decreasing the adjustment of the initial threshold (α parameter), but then the FP detections will be higher. This is actually one of the common limitations within the WML segmentations approaches, in which all of them tend to provide a good trade-off between sensitivity and specificity, assuming some FP and FN detections. Notice also that the inter- intra-rater variability could not be assessed in this study since we only had one manual annotation from a single rater.

We believe that the benefits of an unsupervised approach, which allows to avoid the training stage and therefore having manually annotated cases by experts, will help to the community to quantify WML on Lupus patients, specially considering that we provide a public tool which is straightforward to use in SPM8/12.

5. Conclusions

In this work we have presented an approach to perform WML segmentation on Lupus patients. We have maintained the same pre-processing pipeline applied in Roura et al. (2015), but now implemented in MATLAB code, fact that facilitates the integration to SPM8/12, the installation and execution. The lesion segmentation process has been modified specially on the application of the refinement constraints due to the difference on the lesion features. Results shown in this manuscript demonstrate the good performance of the approach. The correlation results in both number of lesions and lesion volume, illustrates the validity of the approach as a tool for clinicians when diagnosing Lupus patients or evaluating the disease evolution in patients treated with different therapies.

Author Contributions

ER is the corresponder author who has contributed to develop this work and also to write the manuscript. NS has contributed to annotate the data and also to review the manuscript. AO has contributed to develop and also to write this work. SG has contributed to develop this work. SV has contributed to develop this work and also to review the manuscript. RC has contributed to the acquisition of the database. NB has contributed to the acquisition of the database and also to review the paper. XL has contributed to develop and also to write this work.

Conflict of Interest Statement

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.

Acknowledgments

ER holds a BR-UdG2013 Ph.D. grant. SV holds a FI-DGR2013 Ph.D. grant. This work has been supported by “La Fundació la Marató de TV3”, by Retos de Investigación TIN2014-55710-R, and by MPC UdG 2016/022 grant.

Footnotes

1. ^http://www.nitrc.org/projects/cmtk/

2. ^http://eia.udg.edu/salem/slsToolbox/lupus.html

3. ^http://www.mathworks.com/matlabcentral/fileexchange/14995-anisotropic-diffusion--perona---malik-

4. ^http://www.pmod.com/files/download/v35/doc/pneuro/5674.htm

5. ^https://www.nitrc.org/projects/cmtk/

References

Ainiala, H., Dastidar, P., Loukkola, J., Lehtimaki, T., Korpela, M., Peltola, J., and Hietaharju, A. (2005). Cerebral MRI abnormalities and their association with neuropsychiatric manifestations in SLE: a population – based study. Scand. J. Rheumatol. 34, 376–382. doi: 10.1080/03009740510026643

PubMed Abstract | CrossRef Full Text | Google Scholar

Appenzeller, S., Pike, G., and Clarke, A. (2008a). Magnetic resonance imaging in the evaluation of central nervous system manifestations in systemic lupus erythematosus. Clin. Rev. Allergy Immunol. 34, 361–366. doi: 10.1007/s12016-007-8060-z

PubMed Abstract | CrossRef Full Text | Google Scholar

Appenzeller, S., Vasconcelos Faria, A., Li, L., Costallat, L., and Cendes, F. (2008b). Quantitative magnetic resonance imaging analyses and clinical significance of hyperintense white matter lesions in systemic lupus erythematosus patients. Ann. Neurol. 64, 635–643. doi: 10.1002/ana.21483

PubMed Abstract | CrossRef Full Text | Google Scholar

Ashburner, J., and Friston, K. J. (2005). Unified segmentation. Neuroimage 26, 839–851. doi: 10.1016/j.neuroimage.2005.02.018

PubMed Abstract | CrossRef Full Text | Google Scholar

Boesen, K., Rehm, K., Schaper, K., Stoltzner, S., Woods, R., Lüders, E., et al. (2004). Quantitative comparison of four brain extraction algorithms. Neuroimage 22, 1255–1261. doi: 10.1016/j.neuroimage.2004.03.010

PubMed Abstract | CrossRef Full Text | Google Scholar

Brosch, T., Tang, L. Y. W., Yoo, Y., Li, D. K. B., Traboulsee, A., and Tam, R. (2016). Deep 3d convolutional encoder networks with shortcuts for multiscale feature integration applied to multiple sclerosis lesion segmentation. IEEE Trans. Med. Imaging 35, 1229–1239. doi: 10.1109/TMI.2016.2528821

CrossRef Full Text | Google Scholar

Cabezas, M., Oliver, A., Roura, E., Freixenet, J., Vilanova, J., Ramió-Torrentà, L., et al. (2014a). Automatic multiple sclerosis lesion detection in brain MRI by FLAIR thresholding. Comput. Methods Programs Biomed. 115, 147–161. doi: 10.1016/j.cmpb.2014.04.006

PubMed Abstract | CrossRef Full Text | Google Scholar

Cabezas, M., Oliver, A., Valverde, S., Freixenet, J., Beltran, B., Vilanova, J., et al. (2014b). Boost: a supervised approach for multiple sclerosis lesion segmentation. J. Neurosci. Methods 237, 108–117. doi: 10.1016/j.jneumeth.2014.08.024

PubMed Abstract | CrossRef Full Text | Google Scholar

Castellino, G., Govoni, M., Giacuzzo, S., and Trotta, F. (2008). Optimizing clinical monitoring of central nervous system involvement in sle. Autoimmun. Rev. 7, 297–304. doi: 10.1016/j.autrev.2007.11.022

PubMed Abstract | CrossRef Full Text | Google Scholar

Collins, L., Holmes, C. J., Peters, T. M., and Evans, A. C. (1995). Automatic 3-D model-based neuroanatomical segmentation. Hum. Brain Mapp. 3, 190–208. doi: 10.1002/hbm.460030304

CrossRef Full Text | Google Scholar

Fonov, V., Evans, A., Botteron, K., Almli, R., McKinstry, R., and Collins, L. (2011). Unbiased average age-appropriate atlases for pediatric studies. Neuroimage 54, 313–327. doi: 10.1016/j.neuroimage.2010.07.033

PubMed Abstract | CrossRef Full Text | Google Scholar

Guizard, N., Coupé, P., Fonov, V., Manjón, J., Arnold, D., and Collins, L. (2015). Rotation-invariant multi-contrast non-local means for MS lesion segmentation. Neuroimage Clin. 8, 376–389. doi: 10.1016/j.nicl.2015.05.001

PubMed Abstract | CrossRef Full Text | Google Scholar

Hammers, A., Allom, R., Koepp, M., Free, S., Myers, R., Lemieux, L., et al. (2003). Three-dimensional maximum probability atlas of the human brain, with particular reference to the temporal lobe. Hum. Brain Mapp. 19, 224–247. doi: 10.1002/hbm.10123

PubMed Abstract | CrossRef Full Text | Google Scholar

Jain, S., Sima, D. M., Ribbens, A., Cambron, M., Maertens, A., Hecke, W. V., et al. (2015). Automatic segmentation and volumetry of multiple sclerosis brain lesions from MR images. Neuroimage Clin. 8, 367–375. doi: 10.1016/j.nicl.2015.05.003

PubMed Abstract | CrossRef Full Text | Google Scholar

Joseph, F., and Scolding, N. (2009). Neurolupus. Pract. Neurol. 10, 4–15. doi: 10.1136/jnnp.2009.200071

PubMed Abstract | CrossRef Full Text | Google Scholar

Kruggel, F., Paul, J., and Gertz, H. (2008). Texture-based segmentation of diffuse lesions of the brain as white matter. Neuroimage 39, 987–996. doi: 10.1016/j.neuroimage.2007.09.058

PubMed Abstract | CrossRef Full Text | Google Scholar

Larsen, C., Iglesias, J., and Van Leemput, K. (2014). “N3 bias field correction explained as a Bayesian modeling method,” in Bayesian and Graphical Models for Biomedical Imaging, Vol. 8677 of Lecture Notes in Computer Science, eds M. J. Cardoso, I. Simpson, T. Arbel, D. Precup, and A. Ribbens (Cambridge, MA: Springer International Publishing), 1–12.

Lemieux, L., Wieshmann, N., Moran, N., Fish, D., and Shorvon, S. (1998). The detection and significance of subtle changes in mixed-signal brain lesions by serial MRI scan matching and spatial normalization. Med. Image Anal. 2, 227–242.

PubMed Abstract | Google Scholar

Lilliefors, H. W. (1967). On the kolmogorov-smirnov test for normality with mean and variance unknown. J. Am. Stat. Assoc. 62, 399–402.

Google Scholar

Lladó, X., Ganiler, O., Oliver, A., Martí, R., Freixenet, J., Valls, L., et al. (2012a). Automated detection of multiple sclerosis lesions in serial brain MRI. Neuroradiology 54, 787–807. doi: 10.1007/s00234-011-0992-6

PubMed Abstract | CrossRef Full Text | Google Scholar

Lladó, X., Oliver, A., Cabezas, M., Freixenet, J., Vilanova, J. C., and Quiles, A. (2012b). Segmentation of multiple sclerosis lesions in brain MRI: a review of automated approaches. Inf. Sci. 186, 164–185. doi: 10.1016/j.ins.2011.10.011

CrossRef Full Text | Google Scholar

Mazziotta, J. C., Toga, A. W., Evans, A., Fox, P., and Lancaster, J. (1995). A probablistic atlas of the human brain: Theory and rationale for its development. Neuroimage 2, 89–101.

PubMed Abstract | Google Scholar

Mitsias, P., Jacobs, M., Hammoud, R., Pasnoor, M., Santhakumar, S., Papamitsakis, N., et al. (2002). Multiparametric MRI ISODATA ischemic lesion analysis: correlation with the clinical neurological deficit and single-parameter MRI techniques. Stroke 33, 2839–2844.

PubMed Abstract | Google Scholar

Mohamed, F., Vinitski, S., González, C., Faro, S., Lublin, F., Knobler, R., et al. (2001). Increased differentiation of intracranial white matter lesions by multispectral 3D-tissue segmentation: preliminary results. Magn. Reson. Imaging 19, 207–218. doi: 10.1016/S0730-725X(01)00291-0

PubMed Abstract | CrossRef Full Text | Google Scholar

Park, J., and Lee, C. (2009). Skull stripping based on region growing for magnetic resonance brain images. Neuroimage 47, 1394–1407. doi: 10.1016/j.neuroimage.2009.04.047

PubMed Abstract | CrossRef Full Text | Google Scholar

Perona, P., and Malik, J. (1990). Scale-space and edge detection using anisotropic diffusion. IEEE Trans. Pattern Anal. Mach. Intell. 12, 629–639.

Google Scholar

Petri, M., Naqibuddin, M., Carson, K., Wallace, D., Weisman, M., Holliday, S., et al. (2008). Brain magnetic resonance imaging in newly diagnosed systemic lupus erythematosus. J. Rheumatol. 35, 2348–2354. doi: 10.3899/jrheum.071010

PubMed Abstract | CrossRef Full Text | Google Scholar

Roura, E., Oliver, A., Cabezas, M., Valverde, S., Pareto, D., Vilanova, J., et al. (2015). A toolbox for multiple sclerosis lesion segmentation. Neuroradiology 57, 1031–1043. doi: 10.1007/s00234-015-1552-2

PubMed Abstract | CrossRef Full Text | Google Scholar

Roura, E., Oliver, A., Cabezas, M., Vilanova, J., Rovira, À., Ramió-Torrentà, L., et al. (2014). MARGA: Multispectral adaptive region growing algorithm for brain extraction on axial MRI. Comput. Methods Programs Biomed. 113, 655–673. doi: 10.1016/j.cmpb.2013.11.015

PubMed Abstract | CrossRef Full Text | Google Scholar

Sarbu, M., Salman-Monte, T., Rubio Muoz, P., Lisbona, M., Bernabé, M., and Carbonell, J. (2015a). Differences between clinical and laboratory findings in patients with recent diagnosis of SLE according to the positivity of anti-dsDNA by the crithidia luciliae method. Lupus 24, 1198–1203.

PubMed Abstract | Google Scholar

Sarbu, N., Alobeidi, F., Toledano, P., Espinosa, G., Giles, I., Rahman, A., et al. (2015b). Brain abnormalities in newly diagnosed neuropsychiatric lupus: Systematic MRI approach and correlation with clinical and laboratory data in a large multicenter cohort. Autoimmun. Rev. 14, 153–159. doi: 10.1016/j.autrev.2014.11.001

PubMed Abstract | CrossRef Full Text | Google Scholar

Sarbu, N., Bargalló, N., and Cervera, R. (2015c). Advanced and conventional magnetic resonance imaging in neuropsychiatric lupus. F1000Research 4:162. doi: 10.12688/f1000research.6522.2

PubMed Abstract | CrossRef Full Text | Google Scholar

Schmidt, P., Gaser, C., Arsic, M., Buck, D., Förschler, A., Berthele, A., et al. (2012). An automated tool for detection of FLAIR-hyperintense white-matter lesions in Multiple Sclerosis. Neuroimage 59, 3774–3783. doi: 10.1016/j.neuroimage.2011.11.032

PubMed Abstract | CrossRef Full Text | Google Scholar

Schwarz, C., Fletcher, E., DeCarli, C., and Carmichael, O. (2009). “Fully-automated white matter hyperintensity detection with anatomical prior knowledge and without FLAIR,” in Information Processing in Medical Imaging, Vol. 5636 of Lecture Notes in Computer Science, eds J. Prince, D. Pham, and K. Myers (Williamsburg, VA: Springer Berlin Heidelberg), 239–251.

PubMed Abstract

Scully, M., Anderson, B., Lane, T., Gasparovic, C., Magnotta, V., Sibbitt, W., et al. (2010). An automated method for segmenting white matter lesions through multi-level morphometric feature classification with application to lupus. Front. Hum. Neurosci. 4:27. doi: 10.3389/fnhum.2010.00027

PubMed Abstract | CrossRef Full Text | Google Scholar

Sha, D., and Sutton, J. (2001). Towards automated enhancement, segmentation and classification of digital brain images using networks of networks. Inf. Sci. 138, 45–77. doi: 10.1016/S0020-0255(01)00130-X

CrossRef Full Text | Google Scholar

Shattuck, D., Sandor-Leahy, S., Schaper, K., Rottenberg, D., and Leahy, R. (2001). Magnetic resonance image tissue classification using a partial volume model. Neuroimage 13, 856–876. doi: 10.1006/nimg.2000.0730

PubMed Abstract | CrossRef Full Text | Google Scholar

Sibbitt, W. J., Brooks, W., Kornfeld, M., Hart, B., Bankhurst, A., and Roldan, C. (2010). Magnetic resonance imaging and brain histopathology in neuropsychiatric systemic lupus erythematosus. Semin Arthritis Rheum. 40, 35–52. doi: 10.1016/j.semarthrit.2009.08.005

PubMed Abstract | CrossRef Full Text | Google Scholar

Sled, J., Zijdenbos, A., and Evans, A. (1998). A nonparametric method for automatic correction of intensity nonuniformity in MRI data. IEEE Trans. Med. Imaging 17, 87–97.

PubMed Abstract | Google Scholar

Smith, S. (2002). Fast robust automated brain extraction. Hum. Brain Mapp. 17, 143–155. doi: 10.1002/hbm.10062

PubMed Abstract | CrossRef Full Text | Google Scholar

Souplet, J., Lebrun, C., Ayache, N., and Malandain, G. (2008). “An automatic segmentation of T2-FLAIR multiple sclerosis lesions,” in Multiple Sclerosis Lesion Segmentation Challenge Workshop (MICCAI-2008), (New York, NY), 1–8. Available online at: http://hdl.handle.net/10380/1451

PubMed Abstract

Toledano, P., Sarbu, N., Espinosa, G., Bargalló, N., and Cervera, R. (2013). Neuropsychiatric systemic lupus erythematosus: magnetic resonance imaging findings and correlation with clinical and immunological features. Autoimmun. Rev. 12, 1166–1170. doi: 10.1016/j.autrev.2013.07.004

PubMed Abstract | CrossRef Full Text | Google Scholar

Valverde, S., Oliver, A., and Lladó, X. (2014). A white matter lesion-filling approach to improve brain tissue volume measurements. Neuroimage Clin. 6, 86–92. doi: 10.1016/j.nicl.2014.08.016

PubMed Abstract | CrossRef Full Text | Google Scholar

Van Leemput, K., Maes, F., Vandermeulen, D., Colchester, A., and Suetens, P. (2001). Automated segmentation of multiple sclerosis lesions by model outlier detection. IEEE Trans. Med. Imaging 20, 677–688. doi: 10.1109/42.938237

PubMed Abstract | CrossRef Full Text | Google Scholar

Wardlaw, J., Smith, E., Biessels, G., 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

PubMed Abstract | CrossRef Full Text | Google Scholar

Yamashita, Y., Arimura, H., and Tsuchiya, K. (2008). Computer-aided detection of ischemic lesions related to subcortical vascular dementia on magnetic resonance images. Acad. Radiol. 15, 978–985. doi: 10.1016/j.acra.2008.03.004

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, M., Wu, T., and Bennett, K. M. (2015). Small blob identification in medical images using regional features from optimum scale. IEEE Trans. Biomed. Eng. 62, 1051–1062.

PubMed Abstract | Google Scholar