Abstract
Background Dysmyelination could be part of the pathophysiology of schizophrenia spectrum (SCZ) and bipolar disorders (BPD), yet few studies have examined myelination of the cerebral cortex. The ratio of T1- and T2-weighted magnetic resonance images (MRI) correlates with intracortical myelin. We investigated the T1w/T2w-ratio and its age trajectories in patients and healthy controls (CTR) and explored associations with antipsychotic medication use and psychotic symptoms.
Methods Patients with SCZ (n=64; mean age = 30.4 years, SD=9.8), BPD (n=91; mean age 31.0 years, SD=10.2), and CTR (n=155; mean age = 31.9 years, SD=9.1) who participated in the TOP study (NORMENT, University of Oslo, Norway) were clinically assessed and scanned using a General Electric 3T MRI system. T1w/T2w-ratio images were computed using an optimized pipeline with intensity normalization and field inhomogeneity correction. Vertex-wise regression models were used to compare groups and examine group × age interactions. In regions showing significant differences, we explored associations with antipsychotic medication use and psychotic symptoms.
Results No main effect of diagnosis was found. However, age slopes of the T1w/T2w-ratio differed significantly between SCZ and CTR, predominantly in frontal and temporal lobe regions: Lower T1w/T2w-ratio values with higher age were found in CTR, but not in SCZ. Follow-up analyses revealed a more positive age slope in patients who were using antipsychotics and patients using higher chlorpromazine-equivalent doses.
Conclusions While we found no evidence of reduced intracortical myelin in SCZ or BPD relative to CTR, different regional age trajectories in SCZ may suggest a promyelinating effect of antipsychotic medication.
1. Introduction
Schizophrenia spectrum and bipolar disorders are severe mental disorders that affect more than 1% of the population (Perälä et al., 2007). These disorders are proposed to exist along a psychosis continuum (Pearlson, 2015; Tamminga et al., 2014). Psychotic symptoms are also common in bipolar disorders, with an estimated lifetime prevalence above 60% in bipolar I disorder and above 20% in bipolar II disorder (Aminoff et al., 2022; van Bergen et al., 2019). Current evidence favors the view that both disorders reflect dysconnectivity within and across several brain circuitries rather than having focal origin (Friston et al., 2016; Friston and Frith, 1995; Kelly et al., 2018; Xia et al., 2019). Dysmyelination has been proposed as one possible mechanism (Bartzokis, 2002; Whitford et al., 2012). This hypothesis is supported by genetic (Goudriaan et al., 2014; Hakak et al., 2001) and post-mortem studies (Kolomeets and Uranova, 2018; Uranova et al., 2011; Vikhreva et al., 2016) indicating lower myelin content as well as lower density and altered morphology of myelinating oligodendrocytes in the prefrontal cortex (Kolomeets and Uranova, 2018; Uranova et al., 2011; Vikhreva et al., 2016). Notably, the maturation of myelin in association fiber tracts and frontal regions of the cerebral cortex, regions known to be involved in psychotic disorders, extends into early adulthood, which is a period of heightened incidence of psychosis (Paus et al., 2008; Whitford et al., 2012).
Magnetic resonance imaging (MRI) studies on dysmyelination in psychotic disorders have mainly examined white matter, whereas few have examined intracortical myelination (Ganzetti et al., 2015; Iwatani et al., 2015; Wei et al., 2022; Wei et al., 2020). In a previous study, we examined the cortical gray-white matter contrast (GWC), i.e., the contrast between T1-weighted (T1w) intensities in gray matter and adjacent superficial white matter, which is inversely correlated with intracortical myelin. We found higher GWC values in sensory and motor regions in patients with schizophrenia spectrum disorders and, to a lesser extent, in bipolar disorders compared to healthy controls (Jorgensen et al., 2016). More recently, Makowski et al (2019) used structural covariance and principal component analysis on the GWC in patients with first-episode psychosis. They reported a similar trend-level difference in a GWC component representing sensory and motor regions in patients relative to healthy controls. However, the GWC is an indirect measure, and it is necessary to validate these findings using other measures of intracortical myelin.
In bipolar I disorder, Sehmbi et al (2018) found positive associations between a T1w intensity-based measure of intracortical myelin and verbal memory. They also reported an inverted U-shaped relationship with age in healthy controls but not in patients with bipolar I disorder (Sehmbi et al., 2019). Higher T1 relaxation times in bipolar disorders have been found using quantitative MRI (Rangel-Guerra et al., 1983) and in one study this was found in regions including the somatosensory and temporal cortices (Necus et al., 2019).
Myelin may also be involved in the therapeutic action of antipsychotic medications (Bartzokis, 2012; Kroken et al., 2014). In animal studies, a lipogenic effect of antipsychotic agents has been demonstrated (Ersland et al., 2017; Ferno et al., 2011), and in clinical studies an association between serum lipid levels and treatment response was reported (Gjerde et al., 2018a; Kim et al., 2019; Procyshyn et al., 2007). Notably, Tishler et al (2018) investigated a measure of frontal lobe intracortical myelin volume and found a significant association within the first year of antipsychotic medication exposure that declined with prolonged exposure. In a previous study, users of second-generation antipsychotics had a higher intracortical myelin volume compared with users of first-generation antipsychotics (Bartzokis et al., 2007). While these findings are based on an indirect measure, they provide in vivo evidence suggestive of a promyelinating effect of antipsychotic medication in the frontal lobe of the cerebral cortex.
The T1w/T2w-ratio has been proposed as a measure of intracortical myelin (Glasser et al., 2016; Glasser and Van Essen, 2011). This interpretation is based on observations that the T1- and T2-weighted (T2w) MRI signals have positive and negative correlations with myelin content, respectively, such that the contrast due to myelin is enhanced in the ratio (Koenig, 1991; Koenig et al., 1990). Furthermore, shared field inhomogeneities in the T1w and T2w images are attenuated (Glasser and Van Essen, 2011). The T1w/T2w-ratio has been studied in neurological disorders such as multiple sclerosis (Beer et al., 2016), Huntington’s disease (Rowley et al., 2018), and Alzheimer’s disease (Pelkmans et al., 2019). In patients with schizophrenia spectrum disorders, Iwatani et al (2015) reported lower global T1w/T2w-ratio means both in gray and white matter compared to healthy controls, but no voxel-wise group differences in the cerebral cortex. Ganzetti et al (2015) used non-brain intensities to calibrate the T1w/T2w-ratio and found lower regional gray matter values in patients with schizophrenia spectrum disorders, particularly in the temporal lobe, frontal lobe, and the insula. However, two recent studies have indicated a more complex layer-dependent pattern of changes in first-episode psychosis (Wei et al., 2022; Wei et al., 2020). In bipolar disorders, Ishida et al (2017) reported lower T1w/T2w-ratio values in white matter regions relative to healthy controls, but no significant differences in gray matter.
Between-subject comparisons of the T1w/T2w-ratio are challenging due to the non-quantitative nature of MRI signal intensities. In a previous study, we evaluated the measurement properties of 33 T1w/T2w-ratio processing pipelines (Nerland et al., 2021). Correction for field inhomogeneities improved the agreement with the expected myeloarchitecture (i.e., the expected distribution of myelin across cortical areas). Furthermore, intensity normalization ensured acceptable test-retest reliability, which is of particular importance for between-subject comparisons.
In the present study, we investigated cortical T1w/T2w-ratio values in patients with schizophrenia spectrum and bipolar disorders relative to healthy controls. An optimized intensity normalized pipeline was used for computing the T1w/T2w-ratio maps with corrections for partial volume effects, surface outliers, and field inhomogeneities (Nerland et al., 2021). We examined if T1w/T2w-ratio values or age trajectories differed between each patient group and healthy controls with the following two hypotheses: First, that T1w/T2w-ratio values would be lower in primary sensory and motor regions in both patient groups. Second, that use of antipsychotic medication would, particularly in frontal regions of the cerebral cortex, be positively associated with the T1w/T2w-ratio. In exploratory analyses, we examined associations with psychotic symptoms.
2. Methods
2.1 Study design
Participants were recruited from hospitals in the greater Oslo region to the Thematically Organized Psychosis (TOP) study conducted by the Norwegian Centre for Mental Disorders Research (NORMENT). Patients who met the criteria for a DSM-IV schizophrenia spectrum disorder (SCZ), including schizophrenia, schizophreniform disorder and schizo-affective disorder, or bipolar disorders (BPD), including bipolar I disorder, bipolar II disorder and bipolar disorder not otherwise specified, were included in the current study. Healthy controls were randomly drawn from the national population registry in the same geographical region and asked to participate. The study complied with the Helsinki Declaration and was approved by the Regional Committee for Medical Research Ethics (REC South-East Norway) and the Norwegian Data Inspectorate. All participants gave informed consent.
The inclusion criteria for the TOP study include age between 18-65 years, no mental disability (defined as IQ<70), no history of head trauma with loss of consciousness, and no neurological disorder or other organic disorder thought to affect brain function.
Healthy controls were screened using the PRIME-MD (Spitzer et al., 1994). The absence of a mental disorder, substance use disorder, and history of severe mental disorders among first-degree relatives were criteria for inclusion. Healthy controls were selected from a larger pool (n=278) based on age- and sex-matching to the patient sample (SCZ and BPD). Matching was performed with the MatchIt package in R (version 4.2.3; R Core Team, 2018) (Ho et al., 2011). One-to-one matching was performed with the nearest neighbor method and quantile-quantile (Q-Q) plots were inspected to ensure adequate matching.
2.2 Diagnostic and clinical assessment
Clinical assessments were conducted by trained physicians, psychiatrists, or clinical psychologists. Diagnoses were verified using the Structured Clinical Interview for the DSM-IV Axis I disorders (SCID-IV) (Spitzer et al., 1992). Current symptoms were rated using the Positive and Negative Syndrome Scales (PANSS) (Kay et al., 1987). Scores for positive, negative, disorganized, excited, and depressed symptoms were calculated according to the five-factor model by Wallwork et al (2012). Current medication use was obtained by interview or chart review and included information on antipsychotic, antiepileptic, antidepressant, and anxiolytics/hypnotic medication. For each medication, type and dose were recorded. Antipsychotic medication dosages were converted to chlorpromazine equivalent doses (CPZ) (Andreasen et al., 2010). Intelligence quotient (IQ) was assessed with the Wechsler Abbreviated Scale of Intelligence (WASI-II) (Wechsler, 2007), and general psychosocial functioning was rated using the split version of the Global Assessment of Functioning Scale (GAF) (Pedersen et al., 2007). The median time from clinical assessment (defined as the day of the PANSS interview) to MRI acquisition among patients was 12 days, with an interquartile range of 7 – 25 days.
2.3 MRI acquisition
Patients and healthy controls were scanned using a 3T General Electric Discovery MR 750 system, equipped with a 32-channel head coil, between 2015 and 2019. T1w and T2w sequences were both acquired with 1 mm isotropic resolution. The T1w sequence was a 3D inversion recovery-prepared fast spoiled gradient echo recall (BRAVO) sequence with the following parameters: Repetition time (TR) = 8.16 ms; Echo time (TE) = 3.18 ms; Inversion time (TI) = 450 ms; Flip angle = 12°; Bandwidth = 244 Hz/px; ARC = 2; Acquisition time (TA) = 04:43. The T2w sequence was a 3D fast spin echo (CUBE) sequence with the following parameters: TR = 2500 ms; TE = 71.68 ms; FA = 90°; Bandwidth = 488 Hz/px; Echo train length (ETL) = 100; ARC = 2x2; TA = 04:23. Phased array uniformity enhancement (PURE) was enabled for both sequences. MRI images were inspected by a neuroradiologist and excluded if pathological findings were present.
2.4 MRI post-processing
FreeSurfer (v6.0.0; https://surfer.nmr.mgh.harvard.edu/) was used to reconstruct cortical surfaces, representing the boundary between gray and white matter (i.e., the inner gray-white surface of the cortex) and between gray matter and cerebrospinal fluid (i.e., the outer surface of the cortex or ‘pial surface’), based on T1w images.
FreeSurfer is open source and has been described in detail previously (Fischl, 2012). Reconstructed surfaces were visually inspected and edited according to standard guidelines. Images were excluded in the event of substantial motion artifacts or otherwise poor image quality.
2.5 Calculation of the T1w/T2w-ratio
To compute the T1w/T2w-ratio, we rigidly registered T2w images to T1w images using bbregister in FreeSurfer with FSL initialization (Greve and Fischl, 2009). We then applied N4ITK field bias correction and normalized intensities with the WhiteStripe algorithm. The T1w image was then divided by the T2w image to form the T1w/T2w-ratio, which was corrected for partial volume effects (Shafee et al., 2015). Next, T1w/T2w-ratio voxel values were projected onto the gray-white surface by sampling along layers representing equivolumetric distances of 10% to 80% of the vertex-wise cortical thickness (Waehnert et al., 2016). Finally, we performed surface-based outlier correction based on a previously published approach (Glasser and Van Essen, 2011). This pipeline was shown in a previous study to be robust to the presence of field inhomogeneities and to improve test-retest reliability whilst preserving inter-individual variation (Nerland et al., 2021).
We visually inspected each T1w/T2w-ratio map. If the maps deviated from known myeloarchitecture, the T1w and T2w volumes were inspected. If artefacts or low image quality were found in either of the scans, the participant was excluded.
2.6 Statistical analyses
Group differences in demographical and clinical variables were assessed using analysis of variance (ANOVA) or Χ2 tests with post hoc Bonferroni tests where appropriate.
In the primary analyses, we fitted age- and sex-adjusted vertex-wise (i.e., one model per vertex on the fsaverage gray-white surface) general linear models (GLMs) to examine 1) main effects of diagnosis (SCZ versus CTR, BPD versus CTR) and 2) diagnosis × age interaction effects. T1w/T2w-ratio maps were concatenated and smoothed (10 mm FWHM) before running the analyses. The models were then specified using the FreeSurfer command mris_glmfit with two categorical factors, diagnosis (three levels) and sex (two levels), and age as a continuous variable. Cluster-wise correction for multiple testing was employed with a cluster-forming threshold of 0.001 and cluster-wise probability of 5%. Correction for analysis across two hemispheres was applied. If significant diagnosis × age interaction effects were found, we further examined if the age slope differed from zero within each diagnostic group separately. The latter analyses were regarded as follow-up examinations, and we chose a liberal p-value threshold of p<0.01.
Further, we aimed to examine whether T1w/T2w-ratio values were associated with antipsychotic medication use or psychotic symptoms in regions where significant differences were found in the primary analyses. We first extracted the mean T1w/T2w-ratio values for each significant cluster. These were defined as dependent variables. To examine associations with antipsychotic medication use, antipsychotic medication status (current use/no use) as well as medication status × age interaction terms were entered as predictors of interest in the first set of models conducted in the patient sample (n=155). In the second set of models, we examined only patients who were using antipsychotic medication (n=86) and entered CPZ and CPZ × age interaction terms as predictors. We then performed separate analyses among all patients (n=155) using the PANSS total score and each of the five Wallwork factor scores as predictors of interest. All models were adjusted for age, sex, and diagnosis (SCZ/BPD). The analyses were conducted in SPSS version 28.
3. Results
3.1 Description of the study sample
The study sample consisted of 64 patients with SCZ, 91 patients with BPD and 155 CTR. The mean age did not differ between the groups. The sex distribution differed between groups, and post hoc tests indicated a non-significant trend towards more males in the SCZ group and more females in the BPD group.
The distributions of other demographic and clinical variables are shown in Table 1. Briefly, years of education and estimated IQ differed between the groups. Patients with SCZ had higher current PANSS positive, negative, and disorganized symptom scores, but not excited or depressive symptom scores compared to patients with BPD. Patients with SCZ also had lower GAF scores compared to patients with BPD.
Sample characteristics.
3.2 Medication use in the patient sample
The use of antipsychotics was more prevalent in SCZ than in BPD, whereas BPD had a more frequent use of antiepileptic drugs. Other medication categories did not differ between groups (Table 1).
Patients with SCZ used higher doses and were more often treated with multiple antipsychotic agents or long-acting injectables than patients with BPD (Table 1). Further information is found in Supplementary Table 1.
The association between age and current antipsychotic dose (CPZ) was not significant (ρ = 0.18, p = 0.09).
3.3 No group differences in regional T1w/T2w-ratio values
When examining the main effects of diagnosis (SCZ versus CTR, BPD versus CTR), we found no significant differences in regional T1w/T2w-ratio values between either of the patient groups and CTR.
3.4 Different age trajectories of regional T1w/T2w-ratio values
In the group-wise comparison of T1w/T2w-ratio age slopes (i.e., diagnosis x age interaction terms), patients with SCZ had more positive age slopes compared to CTR in 22 clusters. These included clusters in frontal and temporal regions, e.g., bilateral regions of the superior frontal and insular cortices, as well as parietal and occipital regions. There were no significant differences in age slopes between patients with BPD and CTR. An overview of significant clusters is shown in Table 2 and Figure 1. See Supplementary Figure 1 for further details.
Clusters shown in red to yellow colors on the inflated cortical surfaces signify cortical regions where associations with age differed between SCZ and CTR. The cluster-forming threshold was 0.001 and the cluster-wise probability set at p<0.05 with correction for two hemispheres. In the left and right boxes, scatterplots illustrate the association with age among SCZ (red line), BPD (blue line) and CTR (green line) with age on the x-axes and mean T1w/T2-ratio values within the cluster on the y-axes.
Significant clusters based on vertex-wise analysis of schizophrenia × age interaction effect
Follow-up analyses of the linear age slopes in each group showed that CTR had predominantly negative age slopes in medial frontal and temporal regions, with positive age slopes only in the central sulcus. In contrast, patients with SCZ had several regions with positive age slopes, including frontal lobe regions (Figure 2).
The figure displays results from a vertex-wise general linear model adjusted for age and sex. Contrasts were specified to examine age slopes within each group separately. Colored regions indicate where the association with age deviated from the null hypothesis (i.e., no association) at the threshold p<0.01, uncorrected. Blue to light blue colors denote negative age slopes (i.e., lower T1w/T2w-ratio in older individuals). Red to yellow colors denote positive age slopes (i.e., higher T1w/T2w-ratio values in older individuals).
3.5 Age trajectories of the T1w/T2w-ratio and antipsychotic medication use
We found significant interaction effects indicating a more positive age slope in patients using antipsychotic medication compared to patients who did not. This was found in temporal lobe regions, insular regions, the precuneus bilaterally, the left precentral gyrus, the right superior and middle frontal lobe, and the postcentral gyrus (Table 3).
Analysis of cluster-wise T1w/T2w values, antipsychotic medication use and interaction with age.
When analyses were further restricted to patients using antipsychotic medication only (n=86), we found more positive age slopes of T1w/T2w-ratio values in patients using a higher current dose (CPZ). Significant CPZ × age interaction effects were found in all except two of the 22 significant clusters (Table 3).
3.6 No association with clinical symptoms
We found no associations between T1w/T2w-ratio values and PANSS total scores, nor with any of the five symptom factors after correction for multiple testing. See Supplementary Table 2 for details.
4. Discussion
We did not find lower T1w/T2w-ratio values in either of the patient groups compared to healthy controls. Thus, insofar as the T1w/T2w-ratio is a measure of intracortical myelin, our results provide little support for intracortical myelin deficits in these disorders. However, we observed divergent age trajectories in patients with schizophrenia spectrum disorders. Antipsychotic medication status and dose were both associated with divergent age slopes within the patient sample, which is consistent with a possible promyelinating effect of antipsychotic medication.
The absence of lower regional T1w/T2w-ratio values in patients contrasts with our previous study on the GWC (Jorgensen et al., 2016). While these are different measures, they show moderate to high correlations, with a reported overall correlation of 0.73 (Parent et al., 2023). Our results also differed from previous findings of lower global and regional T1w/T2w-ratio values in patients with schizophrenia spectrum disorders (Ganzetti et al., 2015; Iwatani et al., 2015). However, it is worth noting that only the study by Ganzetti et al (2015) reported lower regional T1w/T2w-ratio values. In this study, data was pooled from three different sites, which may have influenced the results given the known effects of scanner on the T1w/T2w-ratio (Nerland et al., 2021). Furthermore, the calibration method based on small masks covering the eyes and the temporal muscles may be unreliable, especially for low-resolution data. In the study by Iwatani et al (2015), a large smoothing kernel was used for computing the T1w/T2w-ratio, and a large portion of the sensorimotor cortices was excluded from the analyses. Notably, both these previous studies used low-resolution T2w images, which may introduce partial volume effects (Shafee et al., 2015).
In two recent studies on first-episode treatment-naïve patients with schizophrenia spectrum disorders (FEP) by Wei and colleagues (2022; 2020), a layer-dependent regional pattern was reported, with lower T1w/T2w-ratio values in the left cingulate and insula and higher values in the left superior temporal gyrus. Notably, these regions belong to the salience network and the language and auditory processing circuitry, respectively, both thought to be affected in psychotic disorders. Interestingly, the patients with FEP differed from healthy controls in the superficial and middle layers of the cortex, but not in the deep layer. In the present study, T1w/T2w-ratio values were sampled at distances of 10-80% of cortical thickness from the gray-white surface, and depth-dependent analyses were not performed. Furthermore, Wei et al (2020) normalized each individual T1w/T2w-ratio map by subtracting the subject-wise mean T1w/T2w-ratio and dividing by the variance, which makes the comparison to the present study difficult.
We observed more positive age trajectories of the T1w/T2w-ratio in patients with schizophrenia spectrum disorders. Furthermore, patients currently treated with antipsychotics showed more positive age trajectories compared with patients not currently on antipsychotics, and there was evidence of a dose-response relationship (higher CPZ associated with a more positive age slope). It has been proposed that some antipsychotic medications have myelin-promoting properties (Bartzokis, 2012) and that their lipogenic effects are related to clinical efficacy (Kim et al., 2019; Leucht et al., 2013; Procyshyn et al., 2007). Indeed, many psychotropic drugs are known to induce cholesterol synthesis (Ferno et al., 2011), and there is evidence to suggest that similar effects are present in human settings (Barth et al., 2020; Bartzokis et al., 2007; Gjerde et al., 2018b; Tishler et al., 2018). Thus, the interpretation that divergent age trajectories in patients reflect antipsychotic treatment is plausible. To directly test this hypothesis, we encourage future studies to assess intracortical myelination longitudinally in patients with FEP who are drug-naïve at baseline.
4.1 Strengths and limitations
Strengths of this study include the use of a clinically well-characterized sample where all participants were scanned on the same MRI system. We used a well-tested pipeline for computing the T1w/T2w-ratio, which showed good test-retest reliability and agreement with known myeloarchitecture using scan acquisitions from the same MRI scanner system and pulse sequence parameters.
To compute the T1w/T2w-ratio, we used an intensity normalization procedure, WhiteStripe, based on intensities in normal-appearing white matter (NAWM). While this was previously shown to improve test-retest reliability whilst preserving individual variation in T1w/T2w-ratio distributions, it may introduce dependencies between cortical T1w/T2w-ratio values and T1w and T2w intensity values in NAWM. We cannot rule out the possibility that the observed age-by-diagnosis interactions reflect age trajectories of NAWM rather than gray matter. Quantitative MRI pulse sequences, such as inversion recovery imaging, may be used to rule out this possibility. Such methods estimate biophysically meaningful properties of the MRI measurements and can be used for between-subject comparisons without the need for intensity normalization or calibration.
The T1w/T2w-ratio shows spatial correlation with cortical myeloarchitecture (Glasser et al., 2014) but is based on T1w and T2w image intensities which are inherently non-dimensional measures. Recent studies have indicated that the correlations between the T1w/T2w-ratio and other indices of myelination vary between brain regions. For instance, low correlations have been found with the myelin-water fraction (MWF) in densely myelinated regions in white matter (Sandrone et al., 2023; Uddin et al., 2019). High correlations have, however, been found between the T1w/T2w-ratio and T1 relaxation time mapping of the cerebral cortex (Parent et al., 2023; Shams et al., 2019). Still, strong conclusions regarding microstructural tissue properties should be avoided since the T1w/T2w-ratio remains a complex measure and other tissue properties than myelin content, such as iron content or dendritic density, may also influence it (Righart et al., 2017).
Since we employed a cross-sectional and naturalistic study design, it was not possible to identify causal mechanisms or investigate longitudinal T1w/T2w-ratio change. Although our study sample was large, given the clinical context, it may have been underpowered to detect subtle alterations in the T1w/T2w-ratio. The T1w and T2w sequences we used both had 1mm isotropic resolution; however, for cortical myelin mapping, images with submillimeter resolution would be optimal. Furthermore, while our optimized pipeline improved reliability (Nerland et al., 2021), correction for field inhomogeneities through the acquisition of B1+ field maps is an alternative approach (Glasser et al., 2022). Lastly, although our findings did not confirm an intracortical myelin deficit in schizophrenia spectrum disorders, such deficits could be present in early illness phases or in treatment-naïve individuals.
4.2 Conclusions
While our findings did not support the hypothesis of intracortical myelin deficits in schizophrenia spectrum or bipolar disorders, they were consistent with the hypothesized promyelinating effect of antipsychotic medication. The possibility that this effect could also mask an intracortical myelin deficit in patients cannot be ruled out. The findings should be followed up by applying quantitative MRI measures and assessing longitudinal trajectories of intracortical myelination in patients who are drug-naïve at baseline.
5. Conflict of interest
OAA has received a speaker’s honorarium from Lundbeck and Sunovion and is a consultant for HealthLytix. IA has received a speaker’s honorarium from Lundbeck. The other authors report no conflict of interest.
Data Availability
The data was collected by the Norwegian Research Centre for Mental Disorders (NORMENT) at the Oslo University Hospital (OUS). The data is subject to restrictions and is not publicly available but may be made available given reasonable request. Data can only be made available following permission from OUS, and insofar requests are in line with the relevant ethical agreements and the consent of the participants.
6. Author contributions
Authors IA and OAA took part in designing the project and data acquisition. Authors KNJ and SN were responsible for the MRI post-processing and quality control. Authors KNJ and SN undertook the statistical analyses. All authors contributed to the interpretation of results. Authors IA, EGJ, OMG and IIM supervised the study. Authors KNJ and SN wrote the initial draft of the manuscript. All authors critically reviewed and provided input towards the final version. All authors have contributed to and approved the manuscript.
7. Role of the funding source
This work was supported by The Research Council of Norway (grant numbers 223273, 274359) and the South-Eastern Norway Regional Health Authority (grant number 2019-104). The funding sources had no further role in study design, data collection, analysis and interpretation of data, writing the manuscript, or submitting the paper for publication.
Acknowledgements
We thank the study participants and the clinicians responsible for recruitment and assessment at the Norwegian Research Centre for Mental Disorders (NORMENT). We also thank the scientific assistants who performed quality assurance and editing of reconstructed surfaces. The work was partly conducted on a platform provided by the Services for sensitive data (TSD), operated and developed at the University of Oslo IT Department (USIT).
References
