Miracle Ozzoude1†
Joel Ramirez1*†
Pradeep Reddy Raamana2
Melissa F. Holmes1Kirstin Walker1Christopher J. M. Scott1Fuqiang Gao1Maged Goubran1,3Donna Kwan4
Maria C. Tartaglia5,6
Derek Beaton2Gustavo Saposnik7Ayman Hassan8Jane Lawrence-Dewar8
Dariush Dowlatshahi9
Stephen C. Strother2,3Sean Symons10
Robert Bartha11Richard H. Swartz12
Sandra E. Black1,12 and on behalf of the ONDRI Investigators- 1LC Campbell Cognitive Neurology Research, Hurvitz Brain Sciences Program, Sunnybrook Research Institute, University of Toronto, Toronto, ON, Canada
- 2Rotman Research Institute, Baycrest Health Sciences, Toronto, ON, Canada
- 3Department of Medical Biophysics, University of Toronto, Toronto, ON, Canada
- 4Centre for Neuroscience Studies, Queens University, Kingston, ON, Canada
- 5Tanz Centre for Research in Neurodegenerative Diseases, University of Toronto, Toronto, ON, Canada
- 6Division of Neurology, Toronto Western Hospital, University Health Network, Toronto, ON, Canada
- 7Stroke Outcomes and Decision Neuroscience Research Unit, Division of Neurology, St. Michael’s Hospital, University of Toronto, Toronto, ON, Canada
- 8Thunder Bay Regional Health Research Institute, Thunder Bay, ON, Canada
- 9Department of Medicine (Neurology), Ottawa Hospital Research Institute, University of Ottawa, Ottawa, ON, Canada
- 10Department of Medical Imaging, Sunnybrook Health Sciences Centre, University of Toronto, Toronto, ON, Canada
- 11Centre for Functional and Metabolic Mapping, Department of Medical Biophysics, Robarts Research Institute, University of Western Ontario, London, ON, Canada
- 12Department of Medicine (Neurology), Sunnybrook Health Sciences Centre, University of Toronto, Toronto, ON, Canada
Background: Regional changes to cortical thickness in individuals with neurodegenerative and cerebrovascular diseases (CVD) can be estimated using specialized neuroimaging software. However, the presence of cerebral small vessel disease, focal atrophy, and cortico-subcortical stroke lesions, pose significant challenges that increase the likelihood of misclassification errors and segmentation failures.
Purpose: The main goal of this study was to examine a correction procedure developed for enhancing FreeSurfer’s (FS’s) cortical thickness estimation tool, particularly when applied to the most challenging MRI obtained from participants with chronic stroke and CVD, with varying degrees of neurovascular lesions and brain atrophy.
Methods: In 155 CVD participants enrolled in the Ontario Neurodegenerative Disease Research Initiative (ONDRI), FS outputs were compared between a fully automated, unmodified procedure and a corrected procedure that accounted for potential sources of error due to atrophy and neurovascular lesions. Quality control (QC) measures were obtained from both procedures. Association between cortical thickness and global cognitive status as assessed by the Montreal Cognitive Assessment (MoCA) score was also investigated from both procedures.
Results: Corrected procedures increased “Acceptable” QC ratings from 18 to 76% for the cortical ribbon and from 38 to 92% for tissue segmentation. Corrected procedures reduced “Fail” ratings from 11 to 0% for the cortical ribbon and 62 to 8% for tissue segmentation. FS-based segmentation of T1-weighted white matter hypointensities were significantly greater in the corrected procedure (5.8 mL vs. 15.9 mL, p < 0.001). The unmodified procedure yielded no significant associations with global cognitive status, whereas the corrected procedure yielded positive associations between MoCA total score and clusters of cortical thickness in the left superior parietal (p = 0.018) and left insula (p = 0.04) regions. Further analyses with the corrected cortical thickness results and MoCA subscores showed a positive association between left superior parietal cortical thickness and Attention (p < 0.001).
Conclusion: These findings suggest that correction procedures which account for brain atrophy and neurovascular lesions can significantly improve FS’s segmentation results and reduce failure rates, thus maximizing power by preventing the loss of our important study participants. Future work will examine relationships between cortical thickness, cerebral small vessel disease, and cognitive dysfunction due to neurodegenerative disease in the ONDRI study.
Introduction
Cortical thickness quantification obtained from magnetic resonance imaging (MRI) can be used to examine regional variations of the cerebral cortex that have been associated with normal aging and dementia due to neurodegeneration (Fischl and Dale, 2000; Salat et al., 2004; Lerch and Evans, 2005; Raamana et al., 2015; Redolfi et al., 2015). Cortical thinning in specific topographical regions of the brain has been used to accurately determine patterns of neurodegeneration in mild cognitive impairment (MCI) (Raamana et al., 2014), Alzheimer’s disease (AD), frontotemporal dementia (FTD) (Du et al., 2007; Lerch et al., 2008; Bakkour et al., 2009; Richards et al., 2009; Cho et al., 2012; Hartikainen et al., 2012; Paternicó et al., 2016; Vuksanović et al., 2019), Parkinson’s disease (Uribe et al., 2016; Gao et al., 2018; Yau et al., 2018; Wilson et al., 2019), amyotrophic lateral sclerosis (ALS) (Verstraete et al., 2012; Mezzapesa et al., 2013; Schuster et al., 2014; Walhout et al., 2015), and vascular cognitive impairment (Seo et al., 2010; Kim et al., 2014; Jung et al., 2018).
FreeSurfer (FS) is a neuroimaging software package that includes a widely used surface-based analysis technique that is able to automatically estimate cortical thickness from T1-weighted MRI (Fischl, 2012). However, degraded image quality and subtle changes introduced by pathology makes it challenging for FS to achieve accurate and reliable brain extraction and white matter (WM) segmentation (Gronenschild et al., 2012; Moore et al., 2012; Iscan et al., 2015; Li et al., 2015a; Backhausen et al., 2016; Schwarz et al., 2016). Although FS provides manual intervention steps to troubleshoot its output (e.g., via control points, WM lesion edits, and pial edits), they are labor-intensive. Further, they may introduce user-bias (Klapwijk et al., 2019), especially in MRI from individuals with significant brain atrophy, cortical stroke lesions, and cerebral small vessel disease. Previous studies examining FS manual correction approaches found that while manual editing may result in differences in morphometrical estimation between the methods in some brain regions (Li et al., 2015b; McCarthy et al., 2015; Canna et al., 2018; Guenette et al., 2018; Waters et al., 2019; Beelen et al., 2020), sensitivity to cognitive-behavioral changes are inconsistent at individual or clinical group levels (McCarthy et al., 2015; Guenette et al., 2018; Waters et al., 2019).
Estimation of cortical thickness in patients with cerebrovascular disease (CVD) can be the most challenging due to cortico-subcortical chronic stroke lesions, significant volumes of white matter hyperintensities (WMH), lacunar infarcts, MRI-visible perivascular spaces (PVS), cortical microinfarcts, and the presence of focal brain atrophy. Given that the performance of FS’s tissue classification is highly dependent on a uniform intensity of voxels in a particular brain region and the integrity of the neighboring voxels, vascular lesions and focal brain atrophy often result in erroneous tissue segmentations, particularly in regions with high surface area and curvature (Fischl et al., 2002). These challenges reduce the accuracy of tissue segmentation, which in turn reduces the accuracy of cortical thickness estimation. Since many age-related neurodegenerative diseases have focal and diffuse brain atrophy that is further exacerbated by comorbid cerebrovascular pathology (Schneider and Bennett, 2010; Khlif et al., 2019a, b), additional procedures to account for these potentially challenging variations in image contrast are needed.
In this paper, we examined results from a FS correction procedure that was applied to MRI obtained from a heterogeneous CVD cohort with varying degrees cerebral small vessel disease, chronic cortico-subcortical stroke lesions, and brain atrophy.
Materials and Methods
Study Participants
Participants (N = 155) recruited to the CVD cohort of the Ontario Neurodegenerative Disease Research Initiative (ONDRI) (Farhan et al., 2017)1 were selected for methodological validation of the FS correction procedure for cortical thickness estimation. The ONDRI study is a multi-modal, multi-site observational research study investigating individuals with neurodegenerative diseases. Study participants were recruited at various health centers across Ontario, Canada: London Health Science Centre and Parkwood Institute in London; Hamilton General Hospital and McMaster Medical Centre in Hamilton; The Ottawa Civic Hospital in Ottawa; Thunder Bay Regional Health Sciences Centre in Thunder Bay; St. Michael’s Hospital, Sunnybrook Health Sciences Centre, Baycrest Health Sciences, Centre for Addiction and Mental Health, and Toronto Western Hospital (University Health Network) in Toronto.
Detailed inclusion and exclusion criteria for the ONDRI CVD participants are previously reported (Farhan et al., 2017; Sunderland et al., 2020). Briefly, participants who had experienced a mild to moderate ischemic stroke event, documented with MRI or CT, over 3 months prior to enrollment, a Modified Rankin Scale (MRS) score (van Swieten et al., 1988) ranging from 0 to 3, and a Montreal Cognitive Assessment (MoCA) score (Nasreddine et al., 2005) ranging 18–30 were included. Participants were excluded if they had severe cognitive impairment, aphasia, a non-vascular cause of symptoms, inability to write or had severe functional disability preventing them to perform assessments, a history of dementia prior to the stroke event, had severe claustrophobia or other contra-indications to MRI procedures. Ethics approval was obtained from all participating institutions and performed in accordance with the Declaration of Helsinki. All participants provided informed consent, and subsequently underwent clinical evaluation, MRI, and other assessments as part of the full ONDRI protocol (Farhan et al., 2017).
MRI Acquisition and Pre-processing
Magnetic resonance imaging protocols were harmonized with the Canadian Dementia Imaging Protocol (CDIP) (Duchesne et al., 2019), and were in compliance with the National Institute of Neurological Disorders and Stroke–Canadian Stroke Network Vascular Cognitive Impairment Harmonization Standards (Hachinski et al., 2006). Detailed MRI protocols are reported elsewhere (Haddad et al., 2020; Ramirez et al., 2020). In brief, the structural MRI used in the current study include: a high-resolution 3D T1-weighted (T1), an interleaved proton density (PD) and T2-weighted (T2), and a T2 fluid-attenuated inversion recovery (FLAIR) images.
Ontario Neurodegenerative Disease Research Initiative’s structural image processing pipeline (Ramirez et al., 2020) will be considered as the pre-processing step for the Corrected FS procedure. ONDRI’s neuroimaging platform used previously published and validated methods, where outputs were further subjected to comprehensive quality control (QC) measures from ONDRI’s neuroinformatic platform using a novel outlier detection algorithm for the identification of anomalous data (Beaton et al., 2019; Sunderland et al., 2019). Briefly, interleaved PD and T2 images and FLAIR images were co-registered to the T1 and a PD-T2 based intra-cranial vault mask was automatically generated and manually edited to create brain segmentation (skull stripping). Using this mask, the T1 was segmented using a multi-feature histogram method to generate a tissue segmentation containing gray matter, WM, and cerebrospinal fluid (CSF) volumes (Kovacevic et al., 2002), at which point ventricular CSF was manually relabeled and the cerebellum was manually removed from the image. WMH of presumed vascular origin and lacunes were automatically identified with FLAIR-based and PD-T2-based lesion segmentation pipelines, respectively and manual edits were applied to remove false positives and recover false negatives (Gibson et al., 2010; Ramirez et al., 2011). Furthermore, the PD-T2 lesion pipeline was used to capture enlarged PVS (Ramirez et al., 2011, 2015). Cortico-subcortical stroke lesions were identified and verified by an expert research neuroradiologist (FG) on T1 and FLAIR images, and manually traced by an experienced neuroimaging analyst (MH). The final output of the pipeline produced a skull-stripped brain mask with segmented voxels comprising of four different “normal tissue” classes and five different “lesion tissue” classes: normal appearing white matter (NAWM), normal appearing gray matter (NAGM), sulcal and ventricular cerebrospinal fluid (sCSF/vCSF), periventricular and deep WMH (pWMH/dWMH), lacunes, PVS, and cortico-subcortical stroke lesions. The skull stripped and lesion-labeled masks were introduced at different processing stages of the Corrected FS procedure described below.
FreeSurfer (FS) Processing Overview
All scans were processed using the stable version of FS (Linux FSv6.0). Two methods were applied to the same participant’s MRI: (a) Unmodified FS and (b) Corrected FS. After applying the two methods, visual inspection was performed by two experienced neuroimaging analysts (MH = rater1; KW = rater2). Raters were blinded to the individual and segmentation method used. Discordant ratings were resolved by a research neuroradiologist (FG). The images were either rated a “pass” or “fail” based on the overall cortical ribbon and tissue segmentation as described in the “Quality Control Assessment Procedures” in the following section.
Unmodified FS: The unmodified procedure involved the standard reconstruction steps in the FS pipeline with the default settings on all participants without any manual interventions. Briefly, the standard reconstruction steps included skull stripping, WM segmentation, intensity normalization, surface reconstruction, subcortical segmentation, cortical parcellation and thickness (Fischl, 2012).
Corrected FS: The corrected procedure involved dividing the reconstruction steps into the following three stages in order to incorporate the skull stripped brain and lesion masks from the ONDRI processing pipeline into FS’s pipeline:
Stage 1 (autorecon1) This involved replacing the “skull stripped mask” (brainmask.mgz) generated by FS’s standard skull stripping method with an improved skull stripped mask from the ONDRI skull stripping method.
Stage 2 (autorecon2) The second intervention (-autorecon2) involved the integration of lesion masks from ONDRI into the initial version of brain tissue segmentation file (aseg.presurf.mgz) generated by FS’s standard segmentation of the brain which includes subcortical structures, WM, GM, CSF, and WM hypointensities. The lesions were given an index value of “77” corresponding to the lesion value in the FS pipeline.
Stage 3 Lastly, the modified aseg.presurf.mgz and brain mask were used as inputs in the last stage of FS pipeline stage 3 (-autorecon2-noaseg -autorecon3) for automatic cortical parcellation and statistics.
Quality Control Assessment Procedures
The accuracy of the cortical ribbon and tissue segmentation from the Unmodified and Corrected FS procedures was evaluated using Freeview, a visualization tool that is packaged with FS. Using the T1 image as the reference, the cortical ribbon accuracy was assessed visually and given a “pass” or “fail” rating based on overestimation or underestimation. Overestimation was defined as a cortical ribbon segmentation that erroneously classified areas of non-brain matter, such as the dura mater or skull, as cortex. Underestimation was defined as a cortical ribbon segmentation that erroneously removed areas of the brain. A “pass” rating was given when the cortical ribbon segmentation showed minimal to no areas of over/underestimation. A “fail” rating was given when there were significant areas of overestimation and/or underestimation (Figures 1a,b), including the complete absence of a cortical ribbon for an entire brain lobule or hemisphere (Figure 1d). Additionally, using the T1 GM-WM intensity differences as a reference, the quality of the tissue segmentation (aseg.mgz) was given a pass/fail rating based on the accuracy of WM-GM boundary. For the tissue segmentation, a “pass” rating was given when the GM-WM boundary showed minimal to no errors in the GM-WM boundary. A “fail” rating was given when there were multiple brain regions where GM was misclassified as WM, and/or WM was misclassified as GM (Figure 1c). Two expert neuroimaging analysts performed the visual ratings and achieved strong inter-rater reliability results for the cortical ribbon segmentation (Cohen’s kappa (k) = 0.9, 95% C.I.: 0.7, 1.00, p < 0.001) and moderate inter-rater reliability results for the tissue segmentation (k = 0.7, 95% C.I.: 0.5, 0.9, p < 0.001) (McHugh, 2012).
Figure 1. (a) An example of a QC “fail” rating due to significant overestimation of the cortical ribbon segmentation. The red arrows point to areas where the cortical ribbon segmentation erroneously classified significant areas of non-brain matter, such as the dura mater or skull, as cortex. (b) An example of a QC “fail” rating due to significant underestimation of the cortical ribbon segmentation. The red arrows point to areas where the cortical ribbon segmentation erroneously removed significant areas of the brain. (c) An example of a QC “fail” rating due to significant tissue segmentation errors. The 1st blue arrow points to an area where GM was misclassified as WM, and the 2nd and 3rd blue arrows point to areas WM was misclassified as GM. (d) An example of a QC “fail” rating due the complete absence of a cortical ribbon for an entire brain lobule or hemisphere. The red arrows point to areas of the brain (temporal lobe and one entire hemisphere) that have been erroneously removed.
Statistical Analysis
Statistical analyses were conducted using Statistical Package for Social Sciences (SPSS v.24) and FS’s packaged analytic software when described. Paired sample t-tests were conducted to determine if the mean lesion volume was significantly different between the unmodified and corrected procedures. This was achieved using the “White matter Hypointensities” identified by FS, which was adjusted for head size using estimated total intracranial volume (eTIV) and log transformed.
A whole brain vertex-wise surface-based cortical thickness analysis was performed on both methods using the built-in general linear model (GLM). Thickness was calculated by the software as the distance between the GM and WM boundaries (also known as the pial surface boundaries) at every vertex in each hemisphere. Each participant’s cortex was anatomically parcellated with every sulcus and gyrus labeled, and resampled to FS’s default average surface map (fsaverage). A 10-mm full-width half-maximum (FWHM) Gaussian spatial smoothing kernel was applied to the surface to improve the signal-to-noise ratio. Age, stroke, and lacunar volumes were included as nuisance regressors. Stroke and lacunar volumes were head size corrected using eTIV.
Montreal Cognitive Assessment total score was included as a regressor of interest to determine the association between cortical thickness and global cognitive status in participants with CVD. Associations between cortical thickness and cognition were further explored using MoCA sub-scores (Visuospatial/Executive, Naming, Memory, Attention, Language, Abstraction, Delayed Recall, and Orientation). Monte Carlo simulations with 5000 iterations were used to correct for multiple comparisons. This method implemented a cluster-wise threshold of 2 (p = 0.01) and cluster-wise probability (p(cwp)) of p < 0.05 (two-sided). Bonferroni correction was applied across the two hemispheres.
Results
Study participant demographics and clinical characteristics are summarized in Table 1. QC results are summarized in Table 2.
Table 1. Study participant demographics and neuroimaging volumetrics (n = 155).
Table 2. Quality control results for FreeSurfer Unmodified and Corrected procedures.
For the cortical ribbon QC, compared to the Unmodified FS procedure, the Corrected “Acceptable” ratings increased from 18 to 76%. For tissue segmentation QC, compared to the Unmodified FS procedure, the Corrected procedure’s “Acceptable” ratings increased from 38 to 92%. For the cortical ribbon QC, the “Fail” ratings were reduced from 11% (Unmodified) to 0% (Corrected). While for the tissue segmentation QC, the “Fail” ratings were reduced from 62 to 8% for Unmodified and Corrected procedures, respectively (e.g., Figure 2). Thus, for the subsequent brain-behavior analysis, we excluded 24 individuals (∼15.5%) that failed due to processing and QC assessment. In order to be excluded due to QC failure, individuals had to fail both the cortical ribbon and tissue segmentation visual assessments. In order to be excluded due to FS processing, the FS pipeline simply did not proceed to completion on these individuals.
Figure 2. Overlaid on a skull-stripped T1, the left image shows the segmentation result from the FS unmodified procedure and the right image shows the segmentation result from the FS corrected procedure. The segmentation colors are as follows: red = GM; green/white = right/left WM; pink = lesion.
When comparing Unmodified and Corrected procedures, results from a paired sample t-test revealed a significant increase (∼63%) in eTIV-adjusted log WM hypointensity volumes, (5824.5 ± 6378.4 mm3 to 15877.1 ± 17964.2 mm3, p < 0.001).
Cortical thickness analyses based on Unmodified FS revealed no significant associations with MoCA total score after accounting for age, stroke, and lacunar volumes. However, the same analyses based on the Corrected FS revealed significant clusters in the left superior parietal and left insula regions were positively associated with MoCA total score (p(cwp) = 0.018; p(cwp) = 0.040, respectively) (Table 3). Further analyses with the Corrected data and MoCA sub-scores using the significant clusters showed a positive association between left superior parietal thickness and the Attention sub-score.
Table 3. Cortical thickness analyses showing significant clusters with Montreal Cognitive Assessment scores corrected for multiple comparisons.
Discussion
To our knowledge, this is the first article that has proposed an improved FS processing in individuals with chronic stroke lesion and cerebral small vessel disease using an automated-hybrid approach. Most published work on improving the FS pipeline have focused on comparing the fully automated pipeline with manual interventions within FS (i.e., editing erroneous skull stripping and control points) (Li et al., 2015a; McCarthy et al., 2015; Canna et al., 2018; Klapwijk et al., 2019; Waters et al., 2019), which are labor intensive, time consuming, and may introduce user-bias. Regional cortical thickness measures obtained from participants’ MRI using FS is a useful imaging biomarker of cortical atrophy, within and between the various disease cohorts represented in ONDRI. By minimizing data loss and increasing statistical power, our correction procedures have the opportunity to develop structural neuroimaging biomarkers for early diagnosis and treatment of cerebrovascular and neurodegenerative disease. This correction procedure enabled the investigation of participants with significant atrophy and cerebrovascular lesion burden, which can present significant challenges to cortical thickness estimation, cortical and subcortical volumetrics, and other downstream processes (e.g., connectivity analyses of functional and diffusion MRI). Moreover, the correction procedures may improve the sensitivity of estimated features that may have otherwise been undetectable.
In the unmodified procedure, a failure rate of more than 60% was reported for tissue segmentation. This is in line with the concept that most segmentation difficulties reported in individuals with CVD result from inaccurate identification of tissue boundaries, which is highly dependent on the homogenous intensity of voxels in a particular brain region, especially in those with high surface area and curvature (Wang et al., 2012). Accurate and reliable skull stripping is important for cortical thickness estimation, since false positive classification of non-brain tissue (e.g., skull, dura, and pial maters) could result in poor estimation of the GM-WM border, which in turn can result in erroneous patterns of cortical thickness. Skull stripping segmentation accuracy is particularly relevant in aging and neurodegenerative populations, where brain atrophy is accompanied by increased CSF volumes and a decreased separation between GM and WM intensities (Seo et al., 2010; Walhout et al., 2015; Yau et al., 2018; Vuksanović et al., 2019; Wilson et al., 2019).
While small acute strokes may have minimal effects on tissue segmentation, large chronic cortico-subcortical stroke lesions introduce alterations to brain morphometry resulting in failed segmentation in most brain segmentation algorithms (Wang et al., 2012; Yang et al., 2016; Siegel et al., 2017; Zavaliangos-Petropulu et al., 2020). Although this issue is particularly relevant in individuals with CVD, cerebral small vessel disease and brain atrophy that are commonly observed in patients with Alzheimer’s and other related dementias present similar challenges when estimating cortical thickness.
Incorporating more accurate brain extraction and lesion masks reduced the overall failure rate from 62% down to less than 8% when the corrected procedure was applied. This improvement could be attributed to the use of multi-modal imaging sequences in the ONDRI structural neuroimaging pipeline (Ramirez et al., 2020) that produces consistent and accurate brain extraction and lesion segmentation. Although imaging markers of small vessel disease, such as WMH, appear hyperintense (bright) on PD/T2 and FLAIR MRI, these lesions appear hypointense or isointense to GM on T1, thus overlapping in intensity with normal appearing GM (Wardlaw et al., 2013). If present in confluent patches, it can result in significant inflation of GM voxel misclassification when using only T1-based segmentation approaches (Levy-Cooperman et al., 2008). Considering the significant WMH burden and atrophy in our sample, it was helpful that the FS pipeline allowed for these types of interventions. In line with this, we found a significant increase in WM hypointensities burden (∼63%) after incorporating ONDRI’s lesion segmentation to the FS pipeline. This is consistent with Magon et al. (2014) that used a lesion filling approach in multiple sclerosis to reclassify the majority of WM lesion in the WM that were erroneous classified as GM or CSF. This method improved the accuracy of tissue segmentation and classification and cortical thickness estimation. Several studies have underscored the importance of optimal lesion segmentation in various clinical population (Habes et al., 2016; Dadar et al., 2017; Sudre et al., 2017; Heinen et al., 2019; Salvadó et al., 2019), particularly in populations at risk of developing small vessel disease (Longstreth et al., 1996; Jeerakathil et al., 2004; Jongen et al., 2007; Habes et al., 2016). A recent systematic review by Frey et al. (2019) provided a comprehensive overview of the importance of WMH segmentation in large-scale MRI studies. They proposed a clear need for developing a guideline to cover the description of WMH segmentation approach, as a way of optimizing the multitude of segmentation techniques available. This is crucial, especially in medium to large sample size studies with clinical populations that donate their time to research. Furthermore, the flexibility of the FS pipeline to allow for such modification supports the individualized imaging methods used in the ONDRI study. This increases the study’s statistical power whilst including participants with challenging pathologies that otherwise might have failed when processed using the default settings, and in turn, reduces sampling bias related to the imaging method requirements (Li et al., 2015a).
Only data that underwent the FS correction demonstrated a relationship with cognition, whereby greater corrected cortical thickness in the left superior parietal cortex and in the insula was associated with higher MoCA total scores. Further analysis with MoCA sub-scores revealed that corrected cortical thickness in the left superior parietal cortex was associated in particular with higher Attention sub-scores. Several studies have reported a significant association between cortical thickness and cognitive function in participants with SVD and other diseases associated with vascular risk factors (Righart et al., 2013; Fujishima et al., 2014; Jung et al., 2014; Tuladhar et al., 2015; Tchistiakova and MacIntosh, 2016). Across these studies the effect of cortical thickness varies, with some reporting relationships with executive function, processing speed, memory (Righart et al., 2013; Tchistiakova et al., 2014), whilst others reporting relationships with memory and attention (Fujishima et al., 2014; Tuladhar et al., 2015). A study by Hilal et al. (2015) demonstrated that WMH and microbleeds were associated with thinning in the temporal and insular regions and associated multi-domain cognitive dysfunction. The insula is an important structure with extensive connections to cortical and subcortical regions, and is involved in various processes, such as empathy, emotion, body sensation, and other aspects of social cognition (Craig, 2009; Couto et al., 2013). Thus, insular atrophy as a result of stroke could lead to significant cognitive dysfunction and socioemotional deficits in participants with cerebral small vessel disease and other comorbid neurodegenerative diseases (Ibañez et al., 2010; Limongi et al., 2014; Moon et al., 2014). Further, the observed association between superior parietal thickness and the Attention sub-score is consistent with recent work showing that smoking-related superior parietal thinning was associated with decreased global cognition, as well as decreased visuospatial and attentional functioning (Neth et al., 2020). This is in line with the concept that better vascular health is associated with increased superior parietal thickness in neurodegenerative diseases (Fortea et al., 2014; Johnson et al., 2014; Vemuri et al., 2018a), suggesting a compensatory response to early brain pathological changes (Vemuri et al., 2018b). Future analyses using our method will investigate the associations between vascular risk factors and cortical thickness in predicting cognitive decline in neurodegenerative diseases with comorbid cerebral small vessel disease.
The ability to decrease the failure rate was the key strength of this work. Although our correction procedures were derived from ONDRI’s imaging pipeline, similar correction procedures that can account for vascular lesions and brain atrophy could be applied in other studies using FS (or any number of cortical thickness estimation tools) (Magon et al., 2014) to study challenging clinical populations (Dalca et al., 2014; Dadar et al., 2017; Subbanna et al., 2019; Ntiri et al., 2020). Hence, the decision to validate and apply this method to individuals with CVD presenting with a range of various combined pathologies including focal and global atrophy, large and small cortico-subcortical chronic stroke lesions, diffuse and focal WMH, lacunar infarcts, cortical microinfarcts, and enlarged PVS (Gorelick et al., 2011; Wardlaw et al., 2013). This combination of brain pathologies brings a unique set of potential challenges for cortical thickness estimation.
The findings reported here should also be considered in light of several limitations. The cross-sectional analysis of this project limits our ability to examine the robustness of our method longitudinally. As ONDRI is a longitudinal study, future work will implement our method at several follow-up time points, within and between all disease cohorts, providing a unique opportunity to investigate relationships between cortical thickness and other neurodegenerative biomarkers for predicting disease progression. As outlined in numerous studies, cortico-subcortical stroke lesions and other neurovascular lesions associated with cerebral small vessel disease are optimally visualized with multi-sequence imaging (minimally with a 3DT1 and T2/FLAIR), in cases where only T1 is available, manual intervention may be the most practical choice for cortical thickness estimation. Another benefit to the FS correction is its potential to facilitate better understanding of brain-behavior relationships by maximizing power since it prevents the loss of our valuable study participants due to poor segmentation and processing failure. As demonstrated, only the corrected cortical estimations correlated with a measure of global cognitive status, which was likely due to the increase in sample size. Future work will examine cross-sectional and longitudinal relationships between cortical thickness, vascular risk factors, neurodegeneration, and associations with more comprehensive neuropsychological testing (McLaughlin et al., 2020).
Conclusion
Accurate topographic measurement of the cortical thickness may benefit early detection and treatment of dementia due to cerebrovascular and neurodegenerative disease. Our results suggest that the integration of the individualized accounting of brain atrophy and neurovascular lesions can significantly improve segmentation results, reduce failure rates to minimize biased samples, and maximize power to examine brain-behavior relationships. Although it is important to acknowledge that the availability of optimized neuroimaging acquisitions may be a limiting factor for implementing this correction procedure, this proof of concept study opens an opportunity for stroke and CVD researchers to potentially examine regional GM thinning in patients who previously could not be analyzed using FS’s cortical thickness tool. Most importantly, these correction efforts invested to reduce data loss and inaccuracies, recognize the significant time and effort our patients have donated to participate in the ONDRI research study. With additional funding, we seek to validate and implement this method in the Canadian Consortium on Neurodegeneration in Aging (CCNA) and other future datasets examining patients with cerebral small vessel disease, focal atrophy, and chronic stroke.
Members of the ONDRI Investigators
Michael Strong, Peter Kleinstiver, Natalie Rashkovan, Susan Bronskil, Michael Borrie, Elizabeth Finger, Corinne Fischer, Andrew Frank, Morris Freedman, Sanjeev Kumar, Stephen Pasternak, Bruce Pollock, Tarek Rajji, Dallas Seitz, David Tang-Wai, Brenda Varriano, Agessandro Abrahao, Marvin Chum, Christen Shoesmith, John Turnbull, Lorne Zinman, Julia Fraser, Bill McIlroy, Ben Cornish, Karen Van Ooteghem, Frederico Faria, Manuel Montero-Odasso, Yanina Sarquis-Adamson, Alanna Black, Barry Greenberg, Wendy Hatch, Chris Hudson, Elena Leontieva, Ed Margolin, Efrem Mandelcorn, Faryan Tayyari, Sherif Defrawy, Don Brien, Ying Chen, Brian Coe, Doug Munoz, Alisia Bonnick, Leanne Casaubon, Jennifer Mandzia, Demetrios Sahlas, David Breen, David Grimes, Mandar Jog, Anthony Lang, Connie Marras, Mario Masellis, Tom Steeves, Dennis Bulman, Allison Ann Dilliott, Mahdi Ghani, Rob Hegele, John Robinson, Ekaterina Rogaeva, Sali Farhan, Hassan Haddad, Nuwan Nanayakkara, Courtney Berezuk, Sabrina Adamo, Mojdeh Zamyadi, Stephen Arnott, Malcolm Binns, Wendy Lou, Kelly Sunderland, Athena Theyers, Abiramy Uthirakumaran, Guangyong (GY) Zou, Sujeevini Sujanthan, Mojdeh Zamyadi, David Munoz, Roger A. Dixon, John Woulfe, Brian Levine, Paula McLaughlin, JB Orange, Alicia Peltsch, Angela Roberts, and Angela Troyer.
Data Availability Statement
The datasets presented in this article are not readily available because the ONDRI data will be made publicly available through an application process on October, 2020. For more information on the ONDRI project, please visit: http://ondri.ca/. Requests to access the datasets should be directed to http://ondri.ca/.
Ethics Statement
The studies involving human participants were reviewed and approved by ONDRI. Study participants were recruited at various health centers across Ontario, Canada: London Health Science Centre and Parkwood Institute in London; Hamilton General Hospital and McMaster Medical Centre in Hamilton; The Ottawa Civic Hospital in Ottawa; Thunder Bay Regional Health Sciences Centre in Thunder Bay; St. Michael’s Hospital, Sunnybrook Health Sciences Centre, Baycrest Health Sciences, Centre for Addiction and Mental Health, and Toronto Western Hospital (University Health Network) in Toronto. Ethics approval was obtained from all participating institutions and performed in accordance with the Declaration of Helsinki. All participants and study partners provided informed consent. The patients/participants provided their written informed consent to participate in this study.
Author Contributions
MO: conceptualization, data curation, formal analysis, investigation, methodology, project administration, software, validation, visualization, and writing (draft, review, and editing). JR: conceptualization, data curation, formal analysis, investigation, methodology, software, validation, visualization, writing (draft, review, and editing), and supervision. PR: data curation, formal analysis, and writing (review and editing). MH and KW: data curation, validation, visualization, and writing (review and editing). CS: data curation, project administration, and writing (review and editing). FG: data curation, formal analysis, validation, resources, and supervision. MG: writing (review and editing). DK: data curation, project administration, writing (review and editing). MT: writing (review and editing), supervision, funding acquisition. DB: data curation, software, writing (review and editing). GS: resources, writing (review and editing). AH: resources, funding acquisition. JL-D: resources, data curation, writing (review and editing). DD and SCS: resources, data curation, funding acquisition. SS: data curation, supervision. RB: data curation, resources, supervision, funding acquisition. RS: data curation, writing (review and editing), resources, supervision, funding acquisition. SB: conceptualization, methodology, supervision, writing (review and editing), resources, funding acquisition. All authors contributed to the article and approved the submitted version.
Funding
This research was conducted with the support of the Ontario Brain Institute, an independent non-profit corporation, funded partially by the Ontario government. The opinions, results, and conclusions are those of the authors and no endorsement by the Ontario Brain Institute is intended or should be inferred. Matching funds were provided by participant hospital and research foundations, including the Baycrest Foundation, Bruyere Research Institute, Centre for Addiction and Mental Health Foundation, London Health Sciences Foundation, McMaster University Faculty of Health Sciences, Ottawa Brain and Mind Research Institute, Queen’s University Faculty of Health Sciences, the Thunder Bay Regional Health Sciences Centre, the University of Ottawa Faculty of Medicine, and the Windsor/Essex County ALS Association. The Temerty Family Foundation provided the major infrastructure matching funds.
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.
Acknowledgments
We would like to thank the ONDRI participants for the time, consent, and participation in our study. Thank you to the L.C. Campbell Foundation, and the analysts and software developers in the L.C. Campbell Cognitive Neurology research team who have contributed to the ONDRI imaging analysis, including Edward Ntiri, Sabrina Adamo, Parisa Mojiri, Rita Meena, and Pugaliya Puveendrakumaran.
Footnotes
References
Backhausen, L. L., Herting, M. M., Buse, J., Roessner, V., Smolka, M. N., and Vetter, N. C. (2016). Quality control of structural mri images applied using FreeSurfer–a hands-on workflow to rate motion artifacts. Front. Neurosci. 10:558. doi: 10.3389/fnins.2016.00558
Bakkour, A., Morris, J. C., and Dickerson, B. C. (2009). The cortical signature of prodromal AD: regional thinning predicts mild AD dementia. Neurology 72, 1048–1055. doi: 10.1212/01.wnl.0000340981.97664.2f
Beaton, D., Sunderland, K. M., Adni, L. B., Mandzia, J., Masellis, M., Swartz, R. H., et al. (2019). Generalization of the minimum covariance determinant algorithm for categorical and mixed data types. bioRxiv[Preprint] doi: 10.1101/333005
Beelen, C., Phan, T. V., Wouters, J., Ghesquière, P., and Vandermosten, M. (2020). Investigating the added value of FreeSurfer’s manual editing procedure for the study of the reading network in a pediatric population. Front. Hum. Neurosci. 14:143. doi: 10.3389/fnhum.2020.00143
Canna, A., Russo, A. G., Ponticorvo, S., Manara, R., Pepino, A., Sansone, M., et al. (2018). Automated search of control points in surface-based morphometry. Neuroimage 176, 56–70. doi: 10.1016/j.neuroimage.2018.04.035
Cho, Y., Seong, J. K., Jeong, Y., and Shin, S. Y. (2012). Individual subject classification for Alzheimer’s disease based on incremental learning using a spatial frequency representation of cortical thickness data. Neuroimage 59, 2217–2230. doi: 10.1016/j.neuroimage.2011.09.085
Couto, B., Sedeño, L., Sposato, L. A., Sigman, M., Riccio, P. M., Salles, A., et al. (2013). Insular networks for emotional processing and social cognition: Comparison of two case reports with either cortical or subcortical involvement. Cortex 49, 1420–1434. doi: 10.1016/j.cortex.2012.08.006
Craig, A. D. B. (2009). How do you feel – now? The anterior insula and human awareness. Nat. Rev. Neurosci. 10, 59–70. doi: 10.1038/nrn2555
Dadar, M., Pascoal, T. A., Manitsirikul, S., Misquitta, K., Fonov, V. S., Tartaglia, M. C., et al. (2017). Validation of a regression technique for segmentation of white matter hyperintensities in Alzheimer’s disease. IEEE Trans. Med. Imaging 36, 1758–1768. doi: 10.1109/TMI.2017.2693978
Dalca, A. V., Sridharan, R., Cloonan, L., Fitzpatrick, K. M., Kanakis, A., Furie, K. L., et al. (2014). Segmentation of cerebrovascular pathologies in stroke patients with spatial and shape priors. Med. Image Comput. Comput. Assist. Interv. 17(Pt 2), 773–780. doi: 10.1007/978-3-319-10470-6_96
Du, A. T., Schuff, N., Kramer, J. H., Rosen, H. J., Gorno-Tempini, M. L., Rankin, K., et al. (2007). Different regional patterns of cortical thinning in Alzheimer’s disease and frontotemporal dementia. Brain 130, 1159–1166. doi: 10.1093/brain/awm016
Duchesne, S., Chouinard, I., Potvin, O., Fonov, V. S., Khademi, A., Bartha, R., et al. (2019). The Canadian dementia imaging protocol: harmonizing national cohorts. J. Magn. Reson. Imaging 49, 456–465. doi: 10.1002/jmri.26197
Farhan, S. M. K., Bartha, R., Black, S. E., Corbett, D., Finger, E., Freedman, M., et al. (2017). The Ontario neurodegenerative disease research initiative (ONDRI). Can. J. Neurol. Sci. 44, 196–202. doi: 10.1017/cjn.2016.415
Fischl, B., and Dale, A. M. (2000). Measuring the thickness of the human cerebral cortex from magnetic resonance images. Proc. Natl. Acad. Sci. U.S.A. 97, 11050–11055. doi: 10.1073/pnas.200033797
Fischl, B., Salat, D. H., Busa, E., Albert, M., Dieterich, M., Haselgrove, C., et al. (2002). Whole brain segmentation: automated labeling of neuroanatomical structures in the human brain. Neuron 33, 341–355.
Fortea, J., Vilaplana, E., Alcolea, D., Carmona-Iragui, M., Sánchez-Saudinos, M. B., Sala, I., et al. (2014). Cerebrospinal fluid β-amyloid and phospho-tau biomarker interactions affecting brain structure in preclinical Alzheimer disease. Ann. Neurol. 76, 223–230. doi: 10.1002/ana.24186
Frey, B. M., Petersen, M., Mayer, C., Schulz, M., Cheng, B., and Thomalla, G. (2019). Characterization of white matter hyperintensities in large-scale MRI-studies. Front. Neurol. 10:238. doi: 10.3389/fneur.2019.00238
Fujishima, M., Maikusa, N., Nakamura, K., Nakatsuka, M., Matsuda, H., and Meguro, K. (2014). Mild cognitive impairment, poor episodic memory, and late-life depression are associated with cerebral cortical thinning and increased white matter hyperintensities. Front. Aging Neurosci. 6:306. doi: 10.3389/fnagi.2014.00306
Gao, Y., Nie, K., Mei, M., Guo, M., Huang, Z., Wang, L., et al. (2018). Changes in cortical thickness in patients with early Parkinson’s disease at different Hoehn and Yahr stages. Front. Hum. Neurosci. 12:469. doi: 10.3389/fnhum.2018.00469
Gibson, E., Gao, F., Black, S. E., and Lobaugh, N. J. (2010). Automatic segmentation of white matter hyperintensities in the elderly using FLAIR images at 3T. J. Magn. Reson. 31, 1311–1322. doi: 10.1002/jmri.22004
Gorelick, P. B., Scuteri, A., Black, S. E., DeCarli, C., Greenberg, S. M., Iadecola, C., et al. (2011). Vascular contributions to cognitive impairment and dementia: a statement for healthcare professionals from the American heart association/american stroke association. Stroke 42, 2672–2713. doi: 10.1161/str.0b013e3182299496
Gronenschild, E. H., Habets, P., Jacobs, H. I., Mengelers, R., Rozendaal, N., van, O. J., et al. (2012). The effects of FreeSurfer version, workstation type, and Macintosh operating system version on anatomical volume and cortical thickness measurements. PLoS One 7:e38234. doi: 10.1371/journal.pone.0038234
Guenette, J. P., Stern, R. A., Tripodis, Y., Chua, A. S., Schultz, V., Sydnor, V. J., et al. (2018). Automated versus manual segmentation of brain region volumes in former football players. NeuroImage Clin. 18, 888–896. doi: 10.1016/j.nicl.2018.03.026
Habes, M., Erus, G., Toledo, J. B., Zhang, T., Bryan, N., Launer, L. J., et al. (2016). White matter hyperintensities and imaging patterns of brain ageing in the general population. Brain 139(Pt 4), 1164–1179. doi: 10.1093/brain/aww008
Hachinski, V., Iadecola, C., Petersen, R. C., Breteler, M. M., Nyenhuis, D. L., Black, S. E., et al. (2006). National institute of neurological disorders and stroke-Canadian stroke network vascular cognitive impairment harmonization standards. Stroke 37, 2220–2241. doi: 10.1161/01.str.0000237236.88823.47
Haddad, S. M. H., Scott, C. J. M., Ozzoude, M., Holmes, M., Arnott, S. R., Nanayakkara, N. D., et al. (2020). Comparison of quality control methods for automated diffusion tensor imaging analysis pipelines. PLoS One 14:e0226715. doi: 10.1371/journal.pone.0226715
Hartikainen, P., Räsänen, J., Julkunen, V., Niskanen, E., Hallikainen, M., Kivipelto, M., et al. (2012). Cortical thickness in frontotemporal dementia, mild cognitive impairment, and Alzheimer’s disease. J. Alzheimers Dis. 30, 857–874. doi: 10.3233/JAD-2012-112060
Heinen, R., Steenwijk, M. D., Barkhof, F., Biesbroek, J. M., van der Flier, W. M., Kuijf, H. J., et al. (2019). Performance of five automated white matter hyperintensity segmentation methods in a multicenter dataset. Sci. Rep. 9:16742. doi: 10.1038/s41598-019-52966-0
Hilal, S., Xin, X., Ang, S. L., Tan, C. S., Venketasubramanian, N., Niessen, W. J., et al. (2015). Risk factors and consequences of cortical thickness in an Asian population. Medicine (Baltimore) 94:e852. doi: 10.1097/MD.960402960402000852
Ibañez, A., Gleichgerrcht, E., and Manes, F. (2010). Clinical effects of insular damage in humans. Brain Struct. Funct. 214, 397–410. doi: 10.1007/s00429-010-0256-y
Iscan, Z., Jin, T. B., Kendrick, A., Szeglin, B., Lu, H., Trivedi, M., et al. (2015). Test-retest reliability of freesurfer measurements within and between sites: Effects of visual approval process. Hum. Brain Mapp. 36, 3472–3485. doi: 10.1002/hbm.22856
Jeerakathil, T., Wolf, P. A., Beiser, A., Hald, J. K., Au, R., Kase, C. S., et al. (2004). Cerebral microbleeds: prevalence and associations with cardiovascular risk factors in the Framingham Study. Stroke 35, 1831–1835. doi: 10.1161/01.str.0000131809.35202.1b
Johnson, S. C., Christian, B. T., Okonkwo, O. C., Oh, J. M., Harding, S., Xu, G., et al. (2014). Amyloid burden and neural function in people at risk for Alzheimer’s disease. Neurobiol. Aging 35, 576–584. doi: 10.1016/j.neurobiolaging.2013.09.028
Jongen, C., van der Grond, J., Kappelle, L. J., Biessels, G. J., Viergever, M. A., and Pluim, J. P. W. (2007). Automated measurement of brain and white matter lesion volume in type 2 diabetes mellitus. Diabetologia 50, 1509–1516. doi: 10.1007/s00125-007-0688-y
Jung, N.-Y., Cho, H., Kim, Y. J., Kim, H. J., Lee, J. M., Park, S., et al. (2018). The impact of education on cortical thickness in amyloid-negative subcortical vascular dementia: cognitive reserve hypothesis. Alzheimers Res. Ther. 10:103. doi: 10.1186/s13195-018-0432-5
Jung, W.-B., Mun, C.-W., Kim, Y.-H., Park, J. M., Lee, B. D., Lee, Y. M., et al. (2014). Cortical atrophy, reduced integrity of white matter and cognitive impairment in subcortical vascular dementia of Binswanger type. Psychiatry Clin. Neurosci. 68, 821–832. doi: 10.1111/pcn.12196
Khlif, M. S., Egorova, N., Werden, E., Redolfi, A., Boccardi, M., DeCarli, C. S., et al. (2019a). A comparison of automated segmentation and manual tracing in estimating hippocampal volume in ischemic stroke and healthy control participants. NeuroImage Clin. 21:101581. doi: 10.1016/j.nicl.2018.10.019
Khlif, M. S., Werden, E., Egorova, N., Boccardi, M., Redolfi, A., Bird, L., et al. (2019b). Assessment of longitudinal hippocampal atrophy in the first year after ischemic stroke using automatic segmentation techniques. NeuroImage Clin. 24:102008. doi: 10.1016/j.nicl.2019.102008
Kim, H. J., Ye, B. S., Yoon, C. W., Noh, Y., Kim, G. H., Cho, H., et al. (2014). Cortical thickness and hippocampal shape in pure vascular mild cognitive impairment and dementia of subcortical type. Eur. J. Neurol. 21, 744–751. doi: 10.1111/ene.12376
Klapwijk, E. T., van de Kamp, F., van der Meulen, M., Peters, S., and Wierenga, L. M. (2019). Qoala-T: a supervised-learning tool for quality control of FreeSurfer segmented MRI data. Neuroimage 189, 116–129. doi: 10.1016/j.neuroimage.2019.01.014
Kovacevic, N., Lobaugh, N. J., Bronskill, M. J., Levine, B., Feinstein, A., and Black, S. E. (2002). A robust method for extraction and automatic segmentation of brain images. Neuroimage 17, 1087–1100. doi: 10.1006/nimg.2002.1221
Lerch, J. P., and Evans, A. C. (2005). Cortical thickness analysis examined through power analysis and a population simulation. Neuroimage 24, 163–173. doi: 10.1016/j.neuroimage.2004.07.045
Lerch, J. P., Pruessner, J., Zijdenbos, A. P., Collins, D. L., Teipel, S. J., Hampel, H., et al. (2008). Automated cortical thickness measurements from MRI can accurately separate Alzheimer’s patients from normal elderly controls. Neurobiol. Aging 29, 23–30. doi: 10.1016/j.neurobiolaging.2006.09.013
Levy-Cooperman, N., Ramirez, J., Lobaugh, N. J., and Black, S. E. (2008). Misclassified tissue volumes in Alzheimer disease patients with white matter hyperintensities: importance of lesion segmentation procedures for volumetric analysis. Stroke 39, 1134–1141. doi: 10.1161/strokeaha.107.498196
Li, Q., Pardoe, H., Lichter, R., Werden, E., Raffelt, A., Cumming, T., et al. (2015b). Cortical thickness estimation in longitudinal stroke studies: a comparison of 3 measurement methods. NeuroImage Clin. 8, 526–535. doi: 10.1016/j.nicl.2014.08.017
Li, Q., Pardoe, H., Lichter, R., Werden, E., Raffelt, A., Cumming, T., et al. (2015a). Cortical thickness estimation in longitudinal stroke studies: a comparison of 3 measurement methods. NeuroImage Clin. 8, 526–535. doi: 10.1016/j.nicl.2014.08.017
Limongi, R., Tomio, A., and Ibanez, A. (2014). Dynamical predictions of insular hubs for social cognition and their application to stroke. Front. Behav. Neurosci. 8:380. doi: 10.3389/fnbeh.2014.00380
Longstreth, W. T. Jr., Manolio, T. A., Arnold, A., Burke, G. L., Bryan, N., Jungreis, C. A., et al. (1996). Clinical correlates of white matter findings on cranial magnetic resonance imaging of 3301 elderly people. The cardiovascular health study. Stroke 27, 1274–1282. doi: 10.1161/01.str.27.8.1274
Magon, S., Gaetano, L., Chakravarty, M. M., Lerch, J. P., Naegelin, Y., Stippich, C., et al. (2014). White matter lesion filling improves the accuracy of cortical thickness measurements in multiple sclerosis patients: a longitudinal study. BMC Neurosci. 15:106. doi: 10.1186/1471-2202-15-106
McCarthy, C. S., Ramprashad, A., Thompson, C., Botti, J.-A., Coman, I. L., and Kates, W. R. (2015). A comparison of FreeSurfer-generated data with and without manual intervention. Front. Neurosci. 9:379. doi: 10.3389/fnins.2015.00379
McHugh, M. L. (2012). Interrater reliability: the kappa statistic. Biochem. Med. 22, 276–282. doi: 10.11613/bm.2012.031
McLaughlin, P. M., Sunderland, K. M., Beaton, D., Binns, M. A., Kwan, D., Levine, B., et al. (2020). The quality assurance and quality control protocol for neuropsychological data collection and curation in the Ontario neurodegenerative disease research initiative (ONDRI) study. Assessment doi: 10.1177/1073191120913933 [Epub ahead of print],
Mezzapesa, D. M., D’Errico, E., Tortelli, R., Distaso, E., Cortese, R., Tursi, M., et al. (2013). Cortical thinning and clinical heterogeneity in amyotrophic lateral sclerosis. PLoS One 8:e80748. doi: 10.1371/journal.pone.0080748
Moon, Y., Moon, W.-J., Kim, H., and Han, S.-H. (2014). Regional atrophy of the insular cortex is associated with neuropsychiatric symptoms in Alzheimer’s disease patients. Eur. Neurol. 71, 223–229. doi: 10.1159/000356343
Moore, D. W., Kovanlikaya, I., Heier, L. A., Raj, A., Huang, C., Chu, K.-W., et al. (2012). Pilot study of quantitative MRI measurements of ventricular volume and cortical atrophy for the differential diagnosis of normal pressure hydrocephalus. Neurol. Res. Int. 2012, 1–6. doi: 10.1155/2012/718150
Nasreddine, Z. S., Phillips, N. A., Bedirian, V., Charbonneau, S., Whitehead, V., Collin, I., et al. (2005). The montreal cognitive assessment, MoCA: a brief screening tool for mild cognitive impairment. J. Am. Geriatr. Soc. 53, 695–699. doi: 10.1111/j.1532-5415.2005.53221.x
Neth, B. J., Graff-Radford, J., Mielke, M. M., Przybelski, S. A., Lesnick, T. G., Schwarz, C. G., et al. (2020). Relationship between risk factors and brain reserve in late middle age: implications for cognitive aging. Front. Aging Neurosci. 11:355. doi: 10.3389/fnagi.2019.00355
Ntiri, E. E., Holmes, M. F., Forooshani, P. M., Ramirez, J., Gao, F., Ozzoude, M., et al. (2020). Neuroinformatics Improved segmentation of the intracranial and ventricular volumes in populations with cerebrovascular lesions and atrophy using 3D CNNs. bioRxiv[Preprint] doi: 10.1101/2020.03.23.000844
Paternicó, D., Manes, M., Premi, E., Cosseddu, M., Gazzina, S., Alberici, A., et al. (2016). Frontotemporal dementia and language networks: cortical thickness reduction is driven by dyslexia susceptibility genes. Sci. Rep. 6:30848. doi: 10.1038/srep30848
Raamana, P. R., Weiner, M. W., Wang, L., and Beg, M. F. (2015). Thickness network features for prognostic applications in dementia. Neurobiol. Aging 36, S91–S102. doi: 10.1016/j.neurobiolaging.2014.05.040
Raamana, P. R., Wen, W., Kochan, N. A., Brodaty, H., Sachdev, P. S., Wang, L., et al. (2014). The sub-classification of amnestic mild cognitive impairment using MRI-based cortical thickness measures. Front. Neurol. 5:76. doi: 10.3389/fneur.2014.00076
Ramirez, J., Berezuk, C., McNeely, A. A., Scott, C. J. M., Gao, F., and Black, S. E. (2015). Visible Virchow-Robin spaces on magnetic resonance imaging of Alzheimer’s disease patients and normal elderly from the Sunnybrook dementia study. J. Alzheimers Dis. 43, 415–424. doi: 10.3233/JAD-132528
Ramirez, J., Gibson, E., Quddus, A., Lobaugh, N. J., Feinstein, A., Levine, B., et al. (2011). Lesion explorer: a comprehensive segmentation and parcellation package to obtain regional volumetrics for subcortical hyperintensities and intracranial tissue. Neuroimage 54, 963–973. doi: 10.1016/j.neuroimage.2010.09.013
Ramirez, J., Holmes, M. F., Scott, C. J. M., Ozzoude, M., Adamo, S., Szilagyi, G. M., et al. (2020). Ontario neurodegenerative disease research initiative (ONDRI): structural MRI methods & outcome measures. Front. Neurol. 11:847. doi: 10.3389/fneur.2020.00847
Redolfi, A., Manset, D., Barkhof, F., Wahlund, L.-O., Glatard, T., Mangin, J.-F., et al. (2015). Head-to-head comparison of two popular cortical thickness extraction algorithms: a cross-sectional and longitudinal study. PLoS One 10:e0117692. doi: 10.1371/journal.pone.0117692
Richards, B. A., Chertkow, H., Singh, V., Robillard, A., Massoud, F., Evans, A. C., et al. (2009). Patterns of cortical thinning in Alzheimer’s disease and frontotemporal dementia. Neurobiol. Aging 30, 1626–1636. doi: 10.1016/j.neurobiolaging.2007.12.019
Righart, R., Duering, M., Gonik, M., Jouvent, E., Reyes, S., Hervé, D., et al. (2013). Impact of regional cortical and subcortical changes on processing speed in cerebral small vessel disease. NeuroImage Clin. 2, 854–861. doi: 10.1016/j.nicl.2013.06.006
Salat, D. H., Buckner, R. L., Snyder, A. Z., Greve, D. N., Desikan, R. S., Busa, E., et al. (2004). Thinning of the cerebral cortex in aging. Cereb. Cortex 14, 721–730. doi: 10.1093/cercor/bhh032
Salvadó, G., Brugulat-Serrat, A., Sudre, C. H., Grau-Rivera, O., Suárez-Calvet, M., Falcon, C., et al. (2019). Spatial patterns of white matter hyperintensities associated with Alzheimer’s disease risk factors in a cognitively healthy middle-aged cohort. Alzheimers Res. Ther. 11:12. doi: 10.1186/s13195-018-0460-1
Schneider, J. A., and Bennett, D. A. (2010). Where vascular meets neurodegenerative disease. Stroke 41, S144–S146.
Schuster, C., Kasper, E., Machts, J., Bittner, D., Kaufmann, J., Benecke, R., et al. (2014). Longitudinal course of cortical thickness decline in amyotrophic lateral sclerosis. J. Neurol. 261, 1871–1880. doi: 10.1007/s00415-014-7426-4
Schwarz, C. G., Gunter, J. L., Wiste, H. J., Przybelski, S. A., Weigand, S. D., Ward, C. P., et al. (2016). A large-scale comparison of cortical thickness and volume methods for measuring Alzheimer’s disease severity. NeuroImage Clin. 11, 802–812. doi: 10.1016/j.nicl.2016.05.017
Seo, S. W., Ahn, J., Yoon, U., Im, K., Lee, J.-M., Tae Kim, S., et al. (2010). Cortical thinning in vascular mild cognitive impairment and vascular dementia of subcortical type. J. Neuroimaging 20, 37–45. doi: 10.1111/j.1552-6569.2008.00293.x
Siegel, J. S., Shulman, G. L., and Corbetta, M. (2017). Measuring functional connectivity in stroke: approaches and considerations. J. Cereb. Blood Flow Metab. 37, 2665–2678. doi: 10.1177/0271678X17709198
Subbanna, N. K., Rajashekar, D., Cheng, B., Thomalla, G., Fiehler, J., Arbel, T., et al. (2019). Stroke lesion segmentation in FLAIR MRI datasets using customized markov random fields. Front. Neurol. 10:541. doi: 10.3389/fneur.2019.00541
Sudre, C. H., Cardoso, M. J., and Ourselin, S. (2017). Longitudinal segmentation of age-related white matter hyperintensities. Med. Image Anal. 38, 50–64. doi: 10.1016/j.media.2017.02.007
Sunderland, K. M., Beaton, D., Arnott, S. R., Kleinstiver, P., Kwan, D., Lawrence-Dewar, J. M., et al. (2020). The Ontario neurodegenerative disease research initiative. medRxiv[Preprint] doi: 10.1101/2020.07.30.20165456
Sunderland, K. M., Beaton, D., Fraser, J., Kwan, D., McLaughlin, P. M., Montero-Odasso, M., et al. (2019). The utility of multivariate outlier detection techniques for data quality evaluation in large studies: an application within the ONDRI project. BMC Med. Res. Methodol. 19:102. doi: 10.1186/s12874-019-0737-5
Tchistiakova, E., Anderson, N. D., Greenwood, C. E., and MacIntosh, B. J. (2014). Combined effects of type 2 diabetes and hypertension associated with cortical thinning and impaired cerebrovascular reactivity relative to hypertension alone in older adults. NeuroImage Clin. 5, 36–41. doi: 10.1016/j.nicl.2014.05.020
Tchistiakova, E., and MacIntosh, B. J. (2016). Summative effects of vascular risk factors on cortical thickness in mild cognitive impairment. Neurobiol. Aging 45, 98–106. doi: 10.1016/j.neurobiolaging.2016.05.011
Tuladhar, A. M., Reid, A. T., Shumskaya, E., de Laat, K. F., van Norden, A. G. W., van Dijk, E. J., et al. (2015). Relationship between white matter hyperintensities, cortical thickness, and cognition. Stroke 46, 425–432. doi: 10.1161/STROKEAHA.114.007146
Uribe, C., Segura, B., Baggio, H. C., Abos, A., Marti, M. J., Valldeoriola, F., et al. (2016). Patterns of cortical thinning in nondemented Parkinson’s disease patients. Mov. Disord. 31, 699–708. doi: 10.1002/mds.26590
van Swieten, J. C., Koudstaal, P. J., Visser, M. C., Schouten, H. J., and van Gijn, J. (1988). Interobserver agreement for the assessment of handicap in stroke patients. Stroke 19, 604–607. doi: 10.1161/01.str.19.5.604
Vemuri, P., Lesnick, T. G., Przybelski, S. A., Graff−Radford, J., Reid, R. I., Lowe, V. J., et al. (2018b). Development of a cerebrovascular magnetic resonance imaging biomarker for cognitive aging. Ann. Neurol. 84, 705–716. doi: 10.1002/ana.25346
Vemuri, P., Lesnick, T. G., Przybelski, S. A., Graff−Radford, J., Reid, R. I., Lowe, V. J., et al. (2018a). Development of a cerebrovascular magnetic resonance imaging biomarker for cognitive aging. Ann. Neurol. 84, 705–716. doi: 10.1002/ana.25346
Verstraete, E., Veldink, J. H., Hendrikse, J., Schelhaas, H. J., van den Heuvel, M. P., and van den Berg, L. H. (2012). Structural MRI reveals cortical thinning in amyotrophic lateral sclerosis. J. Neurol. Neurosurg. Psychiatry 83, 383–388. doi: 10.1136/jnnp-2011-300909
Vuksanović, V., Staff, R. T., Ahearn, T., Murray, A. D., and Wischik, C. M. (2019). Cortical thickness and surface area networks in healthy aging, Alzheimer’s disease and behavioral variant fronto-temporal dementia. Int .J. Neural Syst. 29:1850055. doi: 10.1142/S0129065718500557
Walhout, R., Westeneng, H.-J., Verstraete, E., Hendrikse, J., Veldink, J. H., van den Heuvel, M. P., et al. (2015). Cortical thickness in ALS: towards a marker for upper motor neuron involvement. J. Neurol. Neurosurg. Psychiatry 86, 288–294. doi: 10.1136/jnnp-2013-306839
Wang, Y., Catindig, J. A., Hilal, S., Soon, H. W., Ting, E., Wong, T. Y., et al. (2012). Multi-stage segmentation of white matter hyperintensity, cortical and lacunar infarcts. Neuroimage 60, 2379–2388. doi: 10.1016/j.neuroimage.2012.02.034
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.
Waters, A. B., Mace, R. A., Sawyer, K. S., and Gansler, D. A. (2019). Identifying errors in Freesurfer automated skull stripping and the incremental utility of manual intervention. Brain Imaging Behav. 13, 1281–1291. doi: 10.1007/s11682-018-9951-8
Wilson, H., Niccolini, F., Pellicano, C., and Politis, M. (2019). Cortical thinning across Parkinson’s disease stages and clinical correlates. J. Neurol. Sci. 398, 31–38. doi: 10.1016/j.jns.2019.01.020
Yang, X., Han, X., Park, E., Aylward, S., Kwitt, R., and Niethammer, M. (2016). Registration of pathological images. Simul. Synth. Med. Imaging 9968, 97–107. doi: 10.1007/978-3-319-46630-9_10
Yau, Y., Zeighami, Y., Baker, T. E., Larcher, K., Vainik, U., Dadar, M., et al. (2018). Network connectivity determines cortical thinning in early Parkinson’s disease progression. Nat. Commun. 9:12. doi: 10.1038/s41467-017-02416-0
Keywords: cortical thickness, cerebrovascular disease, stroke, cerebral small vessel disease, FreeSurfer, ONDRI, MRI
Citation: Ozzoude M, Ramirez J, Raamana PR, Holmes MF, Walker K, Scott CJM, Gao F, Goubran M, Kwan D, Tartaglia MC, Beaton D, Saposnik G, Hassan A, Lawrence-Dewar J, Dowlatshahi D, Strother SC, Symons S, Bartha R, Swartz RH and Black SE (2020) Cortical Thickness Estimation in Individuals With Cerebral Small Vessel Disease, Focal Atrophy, and Chronic Stroke Lesions. Front. Neurosci. 14:598868. doi: 10.3389/fnins.2020.598868
Received: 25 August 2020; Accepted: 24 November 2020;
Published: 14 December 2020.
Edited by:
Olivier Colliot, Centre National de la Recherche Scientifique (CNRS), FranceReviewed by:
Baptiste Couvy-Duchesne, University of Queensland, AustraliaAlberto Redolfi, Centro San Giovanni di Dio Fatebenefratelli (IRCCS), Italy
Junhao Wen, University of Pennsylvania, United States
Copyright © 2020 Ozzoude, Ramirez, Raamana, Holmes, Walker, Scott, Gao, Goubran, Kwan, Tartaglia, Beaton, Saposnik, Hassan, Lawrence-Dewar, Dowlatshahi, Strother, Symons, Bartha, Swartz and Black. 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: Joel Ramirez, joel.ramirez1@sunnybrook.ca; joelr@sri.utoronto.ca
†These authors have contributed equally to this work