- Open Access
Reduced striatal volumes in Parkinson’s disease: a magnetic resonance imaging study
Translational Neurodegenerationvolume 1, Article number: 17 (2012)
The presence and extent of structural changes in the brain as a consequence of Parkinson’s disease (PD) is still poorly understood.
High-resolution 3-tesla T1-weighted structural magnetic resonance images in sixty-five PD and 27 age-matched healthy control participants were examined. Putamen, caudate, and intracranial volumes were manually traced in the axial plane of 3D reconstructed images. Striatal nuclei volumes were normalized to intracranial volume for statistical comparison. Disease status was assessed using the Unified Parkinson’s Disease Rating Scale and Hoehn and Yahr scale. Cognitive status was assessed using global status tests and detailed neuropsychological testing.
Both caudate and putamen volumes were smaller in PD brains compared to controls after adjusting for age and gender. Caudate volumes were reduced by 11% (p = 0.001) and putamen volumes by 8.1% (p = 0.025). PD striatal volumes were not found to be significantly correlated with cognitive or motor decline.
Small, but significant reductions in the volume of both the caudate and putamen occur in PD brains. These reductions are independent of the effects of age and gender, however the relation of these reductions to the functional loss of dopamine, which is characteristic of PD, remains unclear.
The striatum (caudate and putamen) is the major input structure of the basal ganglia complex and is an essential part of neural networks involved in motor and non-motor function . Striatal function is severely impaired in Parkinson's disease (PD), which depletes the neuromodulatory influence of ventral midbrain dopamine-producing neurons on these circuits and disrupts the balance of multiple corticostriatal circuits.
It is, however, unclear whether PD also produces gross morphological (volumetric) changes in the striatum. PD manifests progressively worsening motor and non-motor (cognitive and behavioral) dysfunction, which may in part reflect anatomical changes at the level of the putamen (motor) and caudate (oculomotor, cognitive and behavioral). Previous magnetic resonance imaging (MRI) studies have variously reported decreased or non-significant volume differences for these striatal structures using manual or semi-automated tracing methods [2–8] (Table 1).
We employed 3-tesla MRI to estimate in vivo volumes of the striatum, using a large sample that covered a wide range of the motor and cognitive changes evident in PD.
A total of 116 participants, 82 with a diagnosis of probable PD according to the UK Brain Bank criteria  and 34 healthy age-related controls were recruited as part of a continuing research programme. PD participants were recruited through a movement disorder clinic (TJA) and healthy controls were volunteers from the general population. Exclusion criteria were a history of another central nervous system disorder, such as stroke or head injury or a major depressive episode within the last six months. All participants underwent structural magnetic resonance scanning, clinical and cognitive assessment. Six (4 PD and 2 controls) did not meet the requirements of a standard magnetic resonance safety screen and were excluded. A further eighteen were subsequently excluded, five (3 PD and 2 controls) following clinical review of scans, three controls who met the criteria for mild cognitive impairment (MCI) and ten (all PD participants) because of the presence of lacunar infarcts or large perivascular spaces in the striatal or internal capsule regions (n = 6) or movement artefacts (n = 4). Demographic and clinical characteristics of the final cohort (PD n = 65 and control n = 27) included in this study are provided in Table 2. The Upper South Regional Ethics Committees, New Zealand, approved all procedures and all participants gave informed written consent.
PD participants underwent assessment of PD symptom severity and disease stage using the Unified Parkinson’s Disease Rating Scale (UPDRS)  and the modified Hoehn and Yahr scale (H & Y) . Participants were classified as having dementia, mild cognitive impairment (MCI) or being of normal cognition using the following battery of tests. Attention and working memory: digits forward and backward, digit ordering, map search, Trails A, Stroop colour and word naming; Executive function: letter, action and category fluency, and category switching, Trails B and Stroop interference; Visuospatial and perception: judgement of line orientation, fragmented letters and Rey-Osterrieth Complex Figure copy; Learning and memory: California Verbal Learning Test–Short Form and the Rey-Osterrieth Complex Figure Test. A dementia diagnosis followed the Movement Disorders Task Force criteria , requiring significant impairment in activities of daily living  and severe impairment (two or more standard deviations below normative data) in two or more cognitive domains. Mild cognitive impairment was defined as a score of 1.5 or more standard deviations below the normative score in at least two tests in one neuropsychological domain  and meets the Movement Disorders Task Force guidelines for MCI . Global cognitive status was assessed with the Mini Mental State Examination (MMSE)  and Montreal Cognitive Assessment (MoCA)  (not available in 3 PD participants). The levodopa-equivalent dose (LED) for PD participants was calculated using previously published conversion factors [18–20].
Structural magnetic resonance images were acquired with a 3.0-tesla General Electric HDx scanner (GE Healthcare, Milwaukee, WI, USA) using a high-resolution 3D T1-weighted spoiled gradient recalled acquisition sequence: echo time = 2.8 ms; repetition time = 6.6 ms; inversion time = 400 ms; flip angle = 15°, slice thickness 1 mm, field of view = 250 x 250 mm2, acquisition matrix 256 x 256 and voxel size of 0.98 x 0.98 x 1 mm3. Reconstructed images with a matrix of 512 x 512 and voxel size 0.5 x 0.5 x 1 mm3 were used for the manual segmentation of regions of interest (ROI).
Region of interest (ROI) manual tracing
Manual tracing of ROIs was performed by a single rater (TLP) using reconstructed images in MRIcroN  (http://www.mricro.com) implemented on a Windows PC. The rater was blinded to group membership. ROIs were traced in the axial plane with sagittal and coronal planes used for reference using the guidelines below. Each ROI was saved as a single volume, example tracings are shown in Figure 1.
Tracing of the putamen began in the first slice in which the putamen was clearly separated from the head of the caudate by the anterior limb of the internal capsule. Tracing continued on sequential slices for as long as a clearly distinguishable patch of grey matter was visible below the corona radiata. The globus pallidus and anterior limb of the internal capsule served as the medial border and the external capsule as the lateral border. Care was taken to ensure the claustrum was excluded.
Caudate nucleus ROI
The caudate was divided into two parts, head and body, which were traced separately. Tracing of the head of the caudate began in the first slice in which the internal capsule clearly separated the caudate from the putamen and in which the cavity of the lateral ventricle was visible. All gray matter below this slice was considered part of the nucleus accumbens. Small protuberances were included in the measurements of the head of the caudate as described by Looi et al.. The boundary between the caudate head and body was defined, in the coronal view, as the first slice in which the interventricular foramen was present , i.e. when the cavities of the lateral and third ventricle were seen to be continuous. An ROI, covering the entire coronal slice at this level, was created and overlaid on the image, giving a horizontal line across the axial slice. Regions of caudate above to the line were designated the head, and regions below, designated the body. The medial border of the caudate was the cavity of the lateral ventricle. The lateral and superior borders of the caudate were the white matter of the internal capsule and corona radiata, respectively. The gray matter projections between the head of the caudate and the putamen, seen crossing the anterior limb of the internal capsule, were not included in the volume estimates of either nuclei. The caudate nucleus head and body were unable to be accurately delineated in six participants (PD n = 4, controls n = 2) due to the participants head position during scanning. Whole caudate nucleus volumes from these participants were included in the final analysis.
The intracranial volume (ICV) was manually traced in the axial view on every tenth section with intervening slices interpolated . The ICV was traced between the foramen magnum and brain vertex. A starting slice within these boundaries was randomly selected for each brain. The outline was traced following the inner layer of the dura where visible, and the contour of the brain where dura was not visible. The venous sinuses and pituitary fossa were excluded from the ICV calculation.
ROI volume calculation
A custom Matlab (The MathWorks Inc, Natick, USA) script was used to count the number of voxels in each ROI. This value was then multiplied by the reconstructed voxel size to give the ROI volume in mm3. Striatal volumes were normalized to the intracranial volume (bilateral structure volume/ICV x 100%) for statistical comparison.
Statistical analyses were conducted using SPSS 17 (SPSS Inc, Chicago) using normalized volumes (significance was set at p < 0.05). Analysis of covariance was used to compare ICV normalized volumes with group and gender as fixed factors and age as a covariate. Effect sizes were calculated from the estimated marginal means produced by the analysis of covariance. Associations between volumes, demographic and clinical variables were assessed using Pearson’s correlation statistic in a single matrix.
There was no difference in the age (F(1,91) = 1.39; p = 0.241) or years of education (F(1,91) = 0.947; p = 0.333) between PD and control groups (Table 2). There was a significant difference in MoCA scores (F(1,91) = 9.22; p = 0.003), but not MMSE scores (F(1,91) = 3.41; p = 0.068).
ROI measurement reliability
The reliability of volumetric measurements was assessed by randomly selecting 15 scans, duplicating them and including them within the larger sample. The rater was blinded to this process. The reliability of ROI tracing was calculated using raw voxel size measurements from the duplicate scans and the intraclass coefficient statistic. Reliability coefficients were: right putamen, 0.99; left putamen, 0.96; right caudate head, 0.97; right caudate body, 0.87; left caudate head, 0.99; left caudate body, 0.97; and ICV, 0.99.
Intracranial volume did not differ between PD and healthy controls. (Control; mean ± sd. 1524.5 ± 177 cm3, PD; 1546.7 ± 159 cm3; F(1,91) = 0.35, p = 0.56, effect size = −0.13).
The volumes (mm3) of each ROI and of the total (left + right) putamen and caudate are shown in Table 3. To control for global differences in head size, striatal volumes were normalized to the ICV (see Table 3).
Analysis of normalized striatal nuclei volumes indicated that both the caudate and putamen were smaller in PD compared to controls (Table 4), representing a 11% volume reduction in the caudate and an 8.1% volume reduction in the putamen. The model controlled for gender and age, with marginal means calculated assuming an age of 67 years, the average for the entire sample. Although both the caudate head and body showed reductions in volume only the head volume was significantly smaller in PD compared to controls (Table 4).
Both age and gender were found to have significant effects on striatal nuclei volume. Age effects were greatest: caudate; F(1,91) = 14.6, p < 0.001 and putamen; F(1,91) = 39.5, p < 0.001. Regression analysis estimated the volume losses to be 28 mm3/year in the caudate and 57 mm3/year in the putamen. The gender effect indicated that females had larger striatal volumes, with the caudate but not the putamen reaching significance (see Table 4). There were no group by gender interactions: caudate; F(1,91) = 0.325, p = 0.57, putamen; F(1,91) = 0.01, p = 0.92.
Antiparkinsonian medications were not found to have an effect on striatal size: caudate; F(1,64) = 0.30, p = 0.56, putamen; F(1,64) = 2.66, p = 0.11.
Effect of clinical and cognitive measures
Effect sizes slightly smaller than those obtained in the full sample comparison were obtained when comparing cognitively normal PD participants (n = 41) to healthy controls, although the putamen no longer reached significance (see Table 4). There was no significant difference in striatal volumes between the PD cognitive groups (Caudate: PD-normal, 0.419 ± 0.060; PD-MCI, 0.411 ± 0.058; PDD, 0.380 ± 0.057. F(1,64) = 0.64, p = 0.53. Putamen: PD-normal, 0.429 ± 0.068; PD-MCI, 0.393 ± 0.048; PDD, 0.407 ± 0.071. F(1,64) = 0.99, p = 0.38). There were no significant differences in striatal volumes between patients with the initial onset of motor symptoms on the left versus the right (analysis restricted to early stage PD; H & Y stage 1 to 2.5 and cognitively intact). There was no association between striatal volumes and total UPDRS part III motor score. There was also no association when the UPDRS scores were separated into left and right-sided items and correlated with the contralateral striatal volumes. There was an overall correlation between striatal volumes and H & Y stage (Figure 2: caudate: r = −0.36, p = 0.004 and putamen: r = −0.27, p = 0.03). The variability within each H & Y stage meant that there were no significant differences between individual stages. Striatal volumes were also correlated with MoCA scores (Figure 2: caudate: r = 0.32, p = 0.010 and putamen: r = 0.289, p = 0.02) but not MMSE or disease duration.
Given the significant influence of age and gender on striatal volumes the significant correlations were retested using partial correlational analyses controlling for these factors, following which, no significant associations were observed (caudate: H & Y stage; r = −0.20, p = 0.13. MoCA; r = 0.11, p = 039, and putamen: H & Y; r = 0.00, p = 1.0, MoCA; r = 0.03, p = 0.80). The association between striatal volumes and age is shown in Figure 2. Because age, cognitive status and disease severity are highly correlated with each other, it is difficult to determine which factors are primarily driving the association with striatal volume.
Decreased volumes of both the caudate and putamen were found in PD brains compared to healthy controls. This study provides the first 3-tesla evidence of caudate atrophy in PD brains, using manual tracing methods. The present findings contrast with prior studies utilising fewer patients and in which, reduced putamen volumes but not caudate volumes were reported [2–5, 7, 8]. A possible source of variation in the findings across previous studies is variation in imaging parameters, including magnet strength, scanning sequence, slice thickness and inter-slice gap. The present study has the advantage of enhanced image resolution offered by acquisition at 3-tesla compared to 1.5-tesla, and 1 mm-thick slices with no inter slice gap, enabling the full extent of the regions of interest to be viewed. Another potential source of variation may be the boundaries used to delineate each structure. Although the putamen and caudate have well-defined boundaries due to the surrounding white matter and lateral ventricles, these structures fuse anteriorly with the ventral striatum. Unfortunately, most previous studies do not give detailed descriptions of the boundaries used.
A volume reduction in the putamen has previously been accepted due to the reduced dopaminergic activity in the nucleus in early disease stages [3, 4, 7], but little attention has been given to potential mechanisms of such a reduction. Although MRI is not able to directly inform us on such mechanisms, post-mortem studies provide some insight. A reduction in the dendritic spine density and a shortening of dendritic length in the medium spiny neurons of the striatum [25, 26], and frank loss of striatal cells  have been described in post-mortem tissue. Such changes may contribute to a volume reduction and could significantly influence the functionality of the vast cortical projections to the striatum, which synapse on the medium spiny neuron dendritic arbour. Whether these events occur as a direct result of dopamine loss is unclear.
The volume reduction in the caudate was greatest in the head of the nucleus with non-significant volume changes in the body. The head of the caudate is more involved in cognitive function than the body, which is primarily involved in oculomotor function. However, our results do not suggest that the volume loss in the caudate is associated with cognitive status, as we also observed reduced caudate volume when comparing the cognitively unimpaired PD participants to healthy controls and no significant changes in volume between the PD cognitive groups. Whether caudate volume loss contributes to an increased susceptibility to later cognitive decline remains unknown.
Aside from Ghaemi et al. all earlier studies from which we were able to calculate effect sizes indicated the putamen and caudate volumes to be smaller in PD brains compared to control brains (Table 1). The effect size for the caudate volume decrease in the study by Geng et al. was of a reasonable size, but no significant reduction in PD volume was reported, suggesting that the study may have been underpowered to detect the difference in this measurement. Given the consistency in direction of effect sizes across previous studies, and the large sample used in our study, we believe our results reflect true changes in striatal volume in PD.
The absence of an effect of laterality of symptoms (whether side of initial onset or currently most affected side) upon striatal volume is consistent with previous studies . The lack of association between striatal volumes and disease stage indicates that continued loss of volume over time due to disease processes may be of a relatively small magnitude and unable to be teased apart from the losses associated with increasing age. This is in contrast to the previously reported negative association between putamen volume and H & Y stage . Since H & Y stage was significantly correlated with age in this current study, and the correlation between H & Y stage and volumes was not present after correcting for age, the association reported previously may have been, at least in part, due to age rather than more advanced stages of PD. Further investigation into the relationship between disease stage and striatal volume would require greater numbers of younger people with more advanced symptoms. Nigro-striatal neuronal loss is approximately 60% by the time of symptom onset  with rates of loss greater in earlier rather than late stages of the disease. Thus, it should not be surprising that we observed a significant reduction in striatal volumes in patients with normal cognition and no significant change with advancing disease (as measured by H & Y stage).
A limitation to the study is the unbalanced sample sizes in the cognitive sub-groups. Increased numbers of MCI and dementia PD participants would strengthen the analysis based on cognitive status. Also, improved delineation of the sub-regions of the caudate may be achieved in the future with the utilisation of diffusion tensor imaging tractography.
In summary, we have provided evidence at 3-tesla, that the caudate and putamen undergo volume loss in Parkinson’s disease, even early in the disease course. These new data, in conjunction with the existing literature, confirm the presence of atrophy in the striatum in PD brains over and above the influence of natural aging, although not with declining cognitive function and increasing motor impairment. Given the overlap in volumes across groups, such measures would not be useful as a diagnostic tool, however, knowledge of such volume reductions is important for understanding disease effects on the brain.
Mild cognitive impairment
Parkinson’s disease with dementia
Unified parkinson’s disease rating scale
- H & Y:
Hoehn and yahr scale
Region of interest
Mini mental state exam
Montreal cognitive assessment.
Middleton FA, Strick PL: Basal ganglia output and cognition: evidence from anatomical, behavioral, and clinical studies. Brain Cogn 2000, 42: 183-200. 10.1006/brcg.1999.1099
Almeida OP, Burton EJ, McKeith I, Gholkar A, Burn D, O'Brien JT: MRI study of caudate nucleus volume in parkinson's disease with and without dementia with lewy bodies and alzheimer's disease. Demen Geriatr Cogn Disord 2003, 16: 57-63. 10.1159/000070676
Geng DY, Li YX, Zee CS: Magnetic resonance imaging-based volumetric analysis of basal ganglia nuclei and substantia nigra in patients with parkinson's disease. Neurosurgery 2006, 58: 256-262. 10.1227/01.NEU.0000194845.19462.7B
Ghaemi M, Hilker R, Rudolf J, Sobesky J, Heiss WD: Differentiating multiple system atrophy from parkinson's disease: contribution of striatal and midbrain MRI volumetry and multi-tracer PET imaging. J Neurol Neurosurg Psychiatry 2002, 73: 517-523. 10.1136/jnnp.73.5.517
Krabbe K, Karlsborg M, Hansen A, Werdelin L, Mehlsen J, Larsson HB, Paulson OB: Increased intracranial volume in parkinson's disease. J Neurol Sci 2005, 239: 45-52. 10.1016/j.jns.2005.07.013
Lisanby SH, McDonald WM, Massey EW, Doraiswamy PM, Rozear M, Boyko OB, Krishnan KR, Nemeroff C: Diminished subcortical nuclei volumes in parkinson's disease by MR imaging. J Neural Transm Suppl 1993, 40: 13-21.
O'Neill J, Schuff N, Marks WJ, Feiwell R, Aminoff MJ, Weiner MW: Quantitative 1 H magnetic resonance spectroscopy and MRI of Parkinson's disease. Mov Disord 2002, 17: 917-927. 10.1002/mds.10214
Schulz JB, Skalej M, Wedekind D, Luft AR, Abele M, Voigt K, Dichgans J, Klockgether T: Magnetic resonance imaging-based volumetry differentiates idiopathic parkinson's syndrome from multiple system atrophy and progressive supranuclear palsy. Ann Neurol 1999, 45: 65-74. 10.1002/1531-8249(199901)45:1<65::AID-ART12>3.0.CO;2-1
Hughes AJ, Daniel SE, Kilford L, Lees AJ: Accuracy of clinical diagnosis of idiopathic parkinson's disease: a clinico-pathological study of 100 cases. J Neurol Neurosurg Psychiatry 1992, 55: 181-184. 10.1136/jnnp.55.3.181
Fahn S, Elton R, Members of the UPDRS Development Committee: Unified parkinson's disease rating scale. In Recent Developments in Parkinson's Disease. Vol 2. Edited by: Fahn S, Marsden C, Calne D, Goldstein M, Florham Park NJ. Macmillan Healthcare Information; 1987:153-163.
Hoehn MM, Yahr MD: Parkinsonism: onset, progression and mortality. Neurology 1967, 17: 427-442. 10.1212/WNL.17.5.427
Emre M, Aarsland D, Brown R, Burn DJ, Duyckaerts C, Mizuno Y, Broe GA, Cummings J, Dickson DW, Gauthier S, et al.: Clinical diagnostic criteria for dementia associated with parkinson's disease. Move Disord 2007, 22: 1689-1707. 10.1002/mds.21507
Reisberg B, Finkel S, Overall J, Schmidt-Gollas N, Kanowski S, Lehfeld H, Hulla F, Sclan SG, Wilms HU, Heininger K, et al.: The alzheimer's disease activities of daily living international scale (ADL-IS). Int Psychogeriatr 2001, 13: 163-181. 10.1017/S1041610201007566
Dalrymple-Alford JC, Livingston L, MacAskill MR, Graham C, Melzer TR, Porter RJ, Watts R, Anderson TJ: Characterizing mild cognitive impairment in parkinson's disease. Mov Disord 2011, 26: 629-636. 10.1002/mds.23592
Litvan I, Goldman JG, Troster AI, Schmand BA, Weintraub D, Petersen RC, Mollenhauer B, Adler CH, Marder K, Williams-Gray CH, et al.: Diagnostic criteria for mild cognitive impairment in parkinson's disease: movement disorder society task force guidelines. Mov Disord 2012, 27: 349-356. 10.1002/mds.24893
Folstein MF, Folstein SE, McHugh PR: "Mini-mental state". a practical method for grading the cognitive state of patients for the clinician. J Psychiatr Res 1975, 12: 189-198. 10.1016/0022-3956(75)90026-6
Nasreddine ZS, Phillips NA, Bedirian V, Charbonneau S, Whitehead V, Collin I, Cummings JL, Chertkow H: The montreal cognitive assessment, MoCA: a brief screening tool for mild cognitive impairment. J Am Geriatr Soc 2005, 53: 695-699. 10.1111/j.1532-5415.2005.53221.x
Evans AH, Katzenschlager R, Paviour D, O'Sullivan JD, Appel S, Lawrence AD, Lees AJ: Punding in parkinson's disease: its relation to the dopamine dysregulation syndrome. Mov Disord 2004, 19: 397-405. 10.1002/mds.20045
Parkin SG, Gregory RP, Scott R, Bain P, Silburn P, Hall B, Boyle R, Joint C, Aziz TZ: Unilateral and bilateral pallidotomy for idiopathic parkinson's disease: a case series of 115 patients. Mov Disord 2002, 17: 682-692. 10.1002/mds.10186
Tomlinson CL, Stowe R, Patel S, Rick C, Gray R, Clarke CE: Systematic review of levodopa dose equivalency reporting in parkinson's disease. Mov Disord 2010, 25: 2649-2653. 10.1002/mds.23429
Rorden C, Brett M: Stereotaxic display of brain lesions. Behav Neurol 2000, 12: 191-200.
Looi JC, Lindberg O, Liberg B, Tatham V, Kumar R, Maller J, Millard E, Sachdev P, Hogberg G, Pagani M, et al.: Volumetrics of the caudate nucleus: reliability and validity of a new manual tracing protocol. Psychiatry Res 2008, 163: 279-288. 10.1016/j.pscychresns.2007.07.005
Young PA, Young PH: Basic Clinical Neuroanatomy. Williams & Wilkins, Baltimore; 1997.
Eritaia J, Wood SJ, Stuart GW, Bridle N, Dudgeon P, Maruff P, Velakoulis D, Pantelis C: An optimized method for estimating intracranial volume from magnetic resonance images. Magn Reson Med 2000, 44: 973-977. 10.1002/1522-2594(200012)44:6<973::AID-MRM21>3.0.CO;2-H
Stephens B, Mueller AJ, Shering AF, Hood SH, Taggart P, Arbuthnott GW, Bell JE, Kilford L, Kingsbury AE, Daniel SE, Ingham CA: Evidence of a breakdown of corticostriatal connections in parkinson's disease. Neuroscience 2005, 132: 741-754. 10.1016/j.neuroscience.2005.01.007
Zaja-Milatovic S, Milatovic D, Schantz AM, Zhang J, Montine KS, Samii A, Deutch AY, Montine TJ: Dendritic degeneration in neostriatal medium spiny neurons in parkinson disease. Neurology 2005, 64: 545-547. 10.1212/01.WNL.0000150591.33787.A4
Bugiani O, Perdelli F, Salvarani S, Leonardi A, Mancardi GL: Loss of striatal neurons in parkinson's disease: a cytometric study. Eur Neurol 1980, 19: 339-344. 10.1159/000115172
Fearnley JM, Lees AJ: Ageing and parkinson's disease: substantia nigra regional selectivity. Brain 1991, 114(Pt 5):2283-2301.
We thank all the participants and their families for volunteering their time and Tony Harley for helpful discussions. Funding was received from Lottery Health Research, The Neurological Foundation, Canterbury Research Foundation and The Neurology Trust.
The authors declare no competing interests.
TLP - Project: conception and execution; Statistical analysis: design and execution; Manuscript: writing of first draft; TRM – Project: conception and execution; Manuscript: review and critique; MRM – Project: conception; Statistical analysis: review and critique; Manuscript: review and critique; CFG – Project: execution; Manuscript: review and critique; LL - Project: organization and execution; Manuscript: review and critique; RK - Project: execution; Manuscript: review and critique; RW- Project: organization and execution; Manuscript: review and critique; JCD-A - Project: conception, execution; Statistical analysis: review and critique; Manuscript: review and critique; TJA - Project: conception, execution; Manuscript: review and critique. All authors read and approved the final manuscript.