Extreme deviations from the normative model reveal cortical heterogeneity and associations with negative symptom severity in first-episode psychosis from the OPTiMiSE and GAP studies

There is currently no quantifiable method to predict long-term clinical outcomes in patients presenting with a first episode of psychosis. A major barrier to developing useful markers for this is biological heterogeneity, where many different pathological mechanisms may underly the same set of symptoms in different individuals. Normative modelling has been used to quantify this heterogeneity in established psychotic disorders by identifying regions of the cortex which are thinner than expected based on a normative healthy population range. These brain atypicalities are measured at the individual level and therefore potentially useful in a clinical setting. However, it is still unclear whether alterations in individual brain structure can be detected at the time of the first psychotic episode, and whether they are associated with subsequent clinical outcomes. We applied normative modelling of cortical thickness to a sample of first-episode psychosis patients, with the aim of quantifying heterogeneity and to use any pattern of cortical atypicality to predict symptoms and response to antipsychotic medication at timepoints from baseline up to 95 weeks (median follow-ups = 4). T1-weighted brain magnetic resonance images from the GAP and OPTiMiSE samples were processed with Freesurfer V6.0.0 yielding 148 cortical thickness features. An existing normative model of cortical thickness (n = 37,126) was adapted to integrate data from each clinical site and account for effects of gender and site. Our test sample consisted of control participants (n = 149, mean age = 26, SD = 6.7) and patient data (n = 295, mean age = 26, SD = 6.7), this sample was used for estimating deviations from the normative model and subsequent statistical analysis. For each individual, the 148 cortical thickness features were mapped to centiles of the normative distribution and converted to z-scores reflecting the distance from the population mean. Individual cortical thickness metrics of +/- 2.6 standard deviations from the mean were considered extreme deviations from the norm. We found that no more than 6.4% of psychosis patients had extreme deviations in a single brain region (regional overlap) demonstrating a high degree of heterogeneity. Mann-Whitney U tests were run on z-scores for each region and significantly lower z-scores were observed in FEP patients in the frontal, temporal, parietal and occipital lobes. Finally, linear mixed-effects modelling showed that negative deviations in cortical thickness in parietal and temporal regions at baseline are related to more severe negative symptoms over the medium-term. This study shows that even at the early stage of symptom onset normative modelling provides a framework to identify individualised cortical markers which can be used for early personalised intervention and stratification.

The normative modeling framework for computational psychiatry

Normative modeling is an emerging and innovative framework for mapping individual differences at the level of a single subject or observation in relation to a reference model. It involves charting centiles of variation across a population in terms of mappings between biology and behavior, which can then be used to make statistical inferences at the level of the individual. The fields of computational psychiatry and clinical neuroscience have been slow to transition away from patient versus 'healthy' control analytic approaches, probably owing to a lack of tools designed to properly model biological heterogeneity of mental disorders. Normative modeling provides a solution to address this issue and moves analysis away from case-control comparisons that rely on potentially noisy clinical labels. Here we define a standardized protocol to guide users through, from start to finish, normative modeling analysis using the Predictive Clinical Neuroscience toolkit (PCNtoolkit). We describe the input data selection process, provide intuition behind the various modeling choices and conclude by demonstrating several examples of downstream analyses that the normative model may facilitate, such as stratification of high-risk individuals, subtyping and behavioral predictive modeling. The protocol takes ~1-3 h to complete.

Charting brain growth and aging at high spatial precision

Defining reference models for population variation, and the ability to study individual deviations is essential for understanding inter-individual variability and its relation to the onset and progression of medical conditions. In this work, we assembled a reference cohort of neuroimaging data from 82 sites (N=58,836; ages 2-100) and used normative modeling to characterize lifespan trajectories of cortical thickness and subcortical volume. Models are validated against a manually quality checked subset (N=24,354) and we provide an interface for transferring to new data sources. We showcase the clinical value by applying the models to a transdiagnostic psychiatric sample (N=1985), showing they can be used to quantify variability underlying multiple disorders whilst also refining case-control inferences. These models will be augmented with additional samples and imaging modalities as they become available. This provides a common reference platform to bind results from different studies and ultimately paves the way for personalized clinical decision-making.

Subcortical volumes across the lifespan: Data from 18,605 healthy individuals aged 3-90 years

Age has a major effect on brain volume. However, the normative studies available are constrained by small sample sizes, restricted age coverage and significant methodological variability. These limitations introduce inconsistencies and may obscure or distort the lifespan trajectories of brain morphometry. In response, we capitalized on the resources of the Enhancing Neuroimaging Genetics through Meta-Analysis (ENIGMA) Consortium to examine age-related trajectories inferred from cross-sectional measures of the ventricles, the basal ganglia (caudate, putamen, pallidum, and nucleus accumbens), the thalamus, hippocampus and amygdala using magnetic resonance imaging data obtained from 18,605 individuals aged 3-90 years. All subcortical structure volumes were at their maximum value early in life. The volume of the basal ganglia showed a monotonic negative association with age thereafter; there was no significant association between age and the volumes of the thalamus, amygdala and the hippocampus (with some degree of decline in thalamus) until the sixth decade of life after which they also showed a steep negative association with age. The lateral ventricles showed continuous enlargement throughout the lifespan. Age was positively associated with inter-individual variability in the hippocampus and amygdala and the lateral ventricles. These results were robust to potential confounders and could be used to examine the functional significance of deviations from typical age-related morphometric patterns.

Cortical thickness across the lifespan: Data from 17,075 healthy individuals aged 3-90 years

Delineating the association of age and cortical thickness in healthy individuals is critical given the association of cortical thickness with cognition and behavior. Previous research has shown that robust estimates of the association between age and brain morphometry require large-scale studies. In response, we used cross-sectional data from 17,075 individuals aged 3-90 years from the Enhancing Neuroimaging Genetics through Meta-Analysis (ENIGMA) Consortium to infer age-related changes in cortical thickness. We used fractional polynomial (FP) regression to quantify the association between age and cortical thickness, and we computed normalized growth centiles using the parametric Lambda, Mu, and Sigma method. Interindividual variability was estimated using meta-analysis and one-way analysis of variance. For most regions, their highest cortical thickness value was observed in childhood. Age and cortical thickness showed a negative association; the slope was steeper up to the third decade of life and more gradual thereafter; notable exceptions to this general pattern were entorhinal, temporopolar, and anterior cingulate cortices. Interindividual variability was largest in temporal and frontal regions across the lifespan. Age and its FP combinations explained up to 59% variance in cortical thickness. These results may form the basis of further investigation on normative deviation in cortical thickness and its significance for behavioral and cognitive outcomes.

The Normative Modeling Framework for Computational Psychiatry.. [PMID: 35650452 PMCID: PMC7613648 DOI: 10.1101/2021.08.08.455583]

Normative modeling is an emerging and innovative framework for mapping individual differences at the level of a single subject or observation in relation to a reference model. It involves charting centiles of variation across a population in terms of mappings between biology and behavior which can then be used to make statistical inferences at the level of the individual. The fields of computational psychiatry and clinical neuroscience have been slow to transition away from patient versus “healthy” control analytic approaches, likely due to a lack of tools designed to properly model biological heterogeneity of mental disorders. Normative modeling provides a solution to address this issue and moves analysis away from case-control comparisons that rely on potentially noisy clinical labels. In this article, we define a standardized protocol to guide users through, from start to finish, normative modeling analysis using the Predictive Clinical Neuroscience toolkit (PCNtoolkit). We describe the input data selection process, provide intuition behind the various modeling choices, and conclude by demonstrating several examples of down-stream analyses the normative model results may facilitate, such as stratification of high-risk individuals, subtyping, and behavioral predictive modeling. The protocol takes approximately 1-3 hours to complete.

Association of cannabis with glutamatergic levels in patients with early psychosis: Evidence for altered volume striatal glutamate relationships in patients with a history of cannabis use in early psychosis

The associative striatum, an established substrate in psychosis, receives widespread glutamatergic projections. We sought to see if glutamatergic indices are altered between early psychosis patients with and without a history of cannabis use and characterise the relationship to grey matter. 92 participants were scanned: Early Psychosis with a history of cannabis use (EPC = 29); Early Psychosis with minimal cannabis use (EPMC = 25); Controls with a history of cannabis use (HCC = 16) and Controls with minimal use (HCMC = 22). Whole brain T1 weighted MR images and localised proton MR spectra were acquired from head of caudate, anterior cingulate and hippocampus. We examined relationships in regions with known high cannabinoid 1 receptor (CB1R) expression (grey matter, cortex, hippocampus, amygdala) and low expression (white matter, ventricles, brainstem) to caudate Glutamine+Glutamate (Glx). Patients were well matched in symptoms, function and medication. There was no significant group difference in Glx in any region. In EPC grey matter volume explained 31.9% of the variance of caudate Glx (p = 0.003) and amygdala volume explained 36.9% (p = 0.001) of caudate Glx. There was no significant relationship in EPMC. The EPC vs EPMC interaction was significant (p = 0.042). There was no such relationship in control regions. These results are the first to demonstrate association of grey matter volume and striatal glutamate in the EPC group. This may suggest a history of cannabis use leads to a conformational change in distal CB1 rich grey matter regions to influence striatal glutamatergic levels or that such connectivity predisposes to heavy cannabis use.

Grey matter volume and thickness abnormalities in young people with a history of childhood abuse

BACKGROUND

Childhood abuse is associated with abnormalities in brain structure and function. Few studies have investigated abuse-related brain abnormalities in medication-naïve, drug-free youth that also controlled for psychiatric comorbidities by inclusion of a psychiatric control group, which is crucial to disentangle the effects of abuse from those associated with the psychiatric conditions.

METHODS

Cortical volume (CV), cortical thickness (CT) and surface area (SA) were measured in 22 age- and gender-matched medication-naïve youth (aged 13-20) exposed to childhood abuse, 19 psychiatric controls matched for psychiatric diagnoses and 27 healthy controls. Both region-of-interest (ROI) and whole-brain analyses were conducted.

RESULTS

For the ROI analysis, the childhood abuse group compared with healthy controls only, had significantly reduced CV in bilateral cerebellum and reduced CT in left insula and right lateral orbitofrontal cortex (OFC). At the whole-brain level, relative to healthy controls, the childhood abuse group showed significantly reduced CV in left lingual, pericalcarine, precuneus and superior parietal gyri, and reduced CT in left pre-/postcentral and paracentral regions, which furthermore correlated with greater abuse severity. They also had increased CV in left inferior and middle temporal gyri relative to healthy controls. Abnormalities in the precuneus, temporal and precentral regions were abuse-specific relative to psychiatric controls, albeit at a more lenient level. Groups did not differ in SA.

CONCLUSIONS

Childhood abuse is associated with widespread structural abnormalities in OFC-insular, cerebellar, occipital, parietal and temporal regions, which likely underlie the abnormal affective, motivational and cognitive functions typically observed in this population.

Test-retest reliability and longitudinal analysis of automated hippocampal subregion volumes in healthy ageing and Alzheimer's disease populations

The hippocampal formation is a complex brain structure that is important in cognitive processes such as memory, mood, reward processing and other executive functions. Histological and neuroimaging studies have implicated the hippocampal region in neuropsychiatric disorders as well as in neurodegenerative diseases. This highly plastic limbic region is made up of several subregions that are believed to have different functional roles. Therefore, there is a growing interest in imaging the subregions of the hippocampal formation rather than modelling the hippocampus as a homogenous structure, driving the development of new automated analysis tools. Consequently, there is a pressing need to understand the stability of the measures derived from these new techniques. In this study, an automated hippocampal subregion segmentation pipeline, released as a developmental version of Freesurfer (v6.0), was applied to T1-weighted magnetic resonance imaging (MRI) scans of 22 healthy older participants, scanned on 3 separate occasions and a separate longitudinal dataset of 40 Alzheimer's disease (AD) patients. Test-retest reliability of hippocampal subregion volumes was assessed using the intra-class correlation coefficient (ICC), percentage volume difference and percentage volume overlap (Dice). Sensitivity of the regional estimates to longitudinal change was estimated using linear mixed effects (LME) modelling. The results show that out of the 24 hippocampal subregions, 20 had ICC scores of 0.9 or higher in both samples; these regions include the molecular layer, granule cell layer of the dentate gyrus, CA1, CA3 and the subiculum (ICC > 0.9), whilst the hippocampal fissure and fimbria had lower ICC scores (0.73-0.88). Furthermore, LME analysis of the independent AD dataset demonstrated sensitivity to group and individual differences in the rate of volume change over time in several hippocampal subregions (CA1, molecular layer, CA3, hippocampal tail, fissure and presubiculum). These results indicate that this automated segmentation method provides a robust method with which to measure hippocampal subregions, and may be useful in tracking disease progression and measuring the effects of pharmacological intervention.

Brain structure in women at risk of postpartum psychosis: an MRI study

Postpartum psychosis (PP) is the most severe psychiatric disorder associated with childbirth. The risk of PP is very high in women with a history of bipolar affective disorder or schizoaffective disorder. However, the neurobiological basis of PP remains poorly understood and no study has evaluated brain structure in women at risk of, or with, PP. We performed a cross-sectional study of 256 women at risk of PP and 21 healthy controls (HC) in the same postpartum period. Among women at risk, 11 who developed a recent episode of PP (PPE) (n = 2 with lifetime bipolar disorder; n = 9 psychotic disorder not otherwise specified) and 15 at risk women who did not develop an episode of PP (NPPE) (n = 10 with lifetime bipolar disorder; n = 1 with schizoaffective disorder; n = 1 with a history of PP in first-degree family member; n = 3 with previous PP). We obtained T1-weighted MRI scans at 3T and examined regional gray matter volumes with voxel-based morphometry and cortical thickness and surface area with Freesurfer. Women with PPE showed smaller anterior cingulate gyrus, superior temporal gyrus and parahippocampal gyrus compared to NPPE women. These regions also showed decreased surface area. Moreover, the NPPE group showed a larger superior and inferior frontal gyrus volume than the HC. These results should be interpreted with caution, as there were between-group differences in terms of duration of illness and interval between delivery and MRI acquisition. Nevertheless, these are the first findings to suggest that MRI can provide information on brain morphology that characterize those women at risk of PP more likely to develop an episode after childbirth.

Cortical thickness, surface area and volume measures in Parkinson's disease, multiple system atrophy and progressive supranuclear palsy

Objective

Parkinson's disease (PD), Multiple System Atrophy (MSA) and Progressive Supranuclear Palsy (PSP) are neurodegenerative diseases that can be difficult to distinguish clinically. The objective of the current study was to use surface-based analysis techniques to assess cortical thickness, surface area and grey matter volume to identify unique morphological patterns of cortical atrophy in PD, MSA and PSP and to relate these patterns of change to disease duration and clinical features.

Methods

High resolution 3D T1-weighted MRI volumes were acquired from 14 PD patients, 18 MSA, 14 PSP and 19 healthy control participants. Cortical thickness, surface area and volume analyses were carried out using the automated surface-based analysis package FreeSurfer (version 5.1.0). Measures of disease severity and duration were assessed for correlation with cortical morphometric changes in each clinical group.

Results

Results show that in PSP, widespread cortical thinning and volume loss occurs within the frontal lobe, particularly the superior frontal gyrus. In addition, PSP patients also displayed increased surface area in the pericalcarine. In comparison, PD and MSA did not display significant changes in cortical morphology.

Conclusion

These results demonstrate that patients with clinically established PSP exhibit distinct patterns of cortical atrophy, particularly affecting the frontal lobe. These results could be used in the future to develop a useful clinical application of MRI to distinguish PSP patients from PD and MSA patients.

Diffusion tensor imaging of Parkinson's disease, multiple system atrophy and progressive supranuclear palsy: a tract-based spatial statistics study

Although often clinically indistinguishable in the early stages, Parkinson's disease (PD), Multiple System Atrophy (MSA) and Progressive Supranuclear Palsy (PSP) have distinct neuropathological changes. The aim of the current study was to identify white matter tract neurodegeneration characteristic of each of the three syndromes. Tract-based spatial statistics (TBSS) was used to perform a whole-brain automated analysis of diffusion tensor imaging (DTI) data to compare differences in fractional anisotropy (FA) and mean diffusivity (MD) between the three clinical groups and healthy control subjects. Further analyses were conducted to assess the relationship between these putative indices of white matter microstructure and clinical measures of disease severity and symptoms. In PSP, relative to controls, changes in DTI indices consistent with white matter tract degeneration were identified in the corpus callosum, corona radiata, corticospinal tract, superior longitudinal fasciculus, anterior thalamic radiation, superior cerebellar peduncle, medial lemniscus, retrolenticular and anterior limb of the internal capsule, cerebral peduncle and external capsule bilaterally, as well as the left posterior limb of the internal capsule and the right posterior thalamic radiation. MSA patients also displayed differences in the body of the corpus callosum corticospinal tract, cerebellar peduncle, medial lemniscus, anterior and superior corona radiata, posterior limb of the internal capsule external capsule and cerebral peduncle bilaterally, as well as the left anterior limb of the internal capsule and the left anterior thalamic radiation. No significant white matter abnormalities were observed in the PD group. Across groups, MD correlated positively with disease severity in all major white matter tracts. These results show widespread changes in white matter tracts in both PSP and MSA patients, even at a mid-point in the disease process, which are not found in patients with PD.