Aberrant levels of cortical myelin distinguish individuals with unipolar depression from healthy controls
=========================================================================================================

* David A.A. Baranger
* Yaroslav O. Halchenko
* Skye Satz
* Rachel Ragozzino
* Satish Iyengar
* Holly A. Swartz
* Anna Manelis

## ABSTRACT

The association of unipolar depression (UD), relative to healthy controls (HC), with cortical myelin is underexplored, despite growing evidence of associations with white matter tract integrity. We characterized cortical myelin in the 360 Glasser atlas regions using the T1w/T2w ratio in 39 UD and 47 HC participants (ages=19-44, 75% female). A logistic elastic net regularized regression with nested cross-validation and a subsequent linear discriminant analysis conducted on held-out samples were used to select brain regions and classify UD vs. HC. True-label model performance was compared against permuted-label model performance. Cortical myelin distinguished UD from HC with 68% accuracy (p<0.001; sensitivity=63.8%, specificity=71.5%). Brain regions contributing to this classification performance were located in the orbitofrontal cortex, anterior cingulate, extended visual, and auditory cortices, and showed statistically significant decreases and increases in myelin levels in UD vs. HC. The patterns of cortical myelin in these regions may be a biomarker of UD.

## 1 INTRODUCTION

Unipolar depression (UD) is a leading cause of disability worldwide (1), with an economic burden of $210 billion dollars in the United States alone (2). Despite its impact, treatments for the disorder remain ineffective for many patients (3). Thus, there is a pressing need to understand the neurobiological etiology of UD to facilitate the development of improved treatments and prevention strategies.

Unipolar depression is characterized by dysfunctional affective and cognitive processing (4), including reduced executive functioning (5), biased emotional processing (6), and impaired reward processing (7). Correspondingly, individuals with UD show aberrant activation during tasks which recruit these processes, including activation in the striatum, hippocampus, amygdala, orbitofrontal cortex, prefrontal cortex, insula, cingulate, and occipital cortex (8–11). In addition, a growing body of literature has reported structural abnormalities associated with depression in many of these same regions described above, both in grey matter (12–14) and in white matter (15,16). Meta-analyses of diffusion weighted imaging (DWI) studies have repeatedly found evidence for lower fractional anisotropy (FA) in depressed populations (17–19). More recently, studies using very large samples (i.e., the UK Biobank), as well as meta-analyses combining both published and unpublished data (i.e., the ENIGMA consortium), have observed widespread and replicable reductions in FA(15,16). Notably, white matter integrity in identified regions has been shown to correlate with the cognitive processes disrupted in depression, including processing speed (20,21), emotion regulation (22,23), and reward learning (24).

Emerging evidence suggests that cortical myelin may be impacted in individuals with UD and that it may partially mediate some of the cognitive processes that are impaired in affected individuals. For example, studies of *post mortem* brain tissue from donors with UD have observed reduced myelination, a reduction in the number of oligodendrocytes and other glia (cells whose functions include generating and maintaining myelin (25)), and reduced expression of oligodendrocyte lineage genes (26,27). A study of individuals with treatment resistant depression revealed a reduced magnetization transfer ratio (MTR), which is thought to reflect lower myelin levels, in the cingulate cortex and insula (28). A recent study using R1 (1/T1) as a measure of myelination observed reduced whole-brain myelin, but no significant difference in cortical myelin in a handful of *a priori* bilateral regions between individuals with depression and healthy controls (HC) (29),

Developments in magnetic resonance imaging (MRI) methodology permit the examination of cortical myelin via the T1w/T2w ratio (30,31). Studies in population-based samples using this metric have found that lower myelin in the cingulate, orbitofrontal cortex, and middle temporal cortex correlated with poor sleep quality (32), lower frontal-pole myelin and greater myelin in the occipital cortex correlated with neuroticism (33), and lower myelin in the motor and higher myelin in the insular, cingulate, prefrontal, and superior parietal cortices correlated with trait anxiety (34). While poor sleep, high neuroticism, and trait anxiety might represent concurrent symptoms of depression, prior studies have not systematically examined cortical myelin in participants with depression as compared to HC.

The goals of the present study were (1) to ascertain whether cortical myelin content (characterized by the T1w/T2w ratio) is predictive of unipolar depression (UD), and (2) to characterize the brain regions that are predictive of case/control status. Based on the prior studies mentioned above, we hypothesized that cortical myelin levels would distinguish individuals with UD from HC and that these differences will be especially pronounced in the prefrontal cortical (PFC), cingulate, parietal and occipital regions that support reward and emotional processing, which are dysregulated in UD (11,35).

## 2 METHODS

### 2.1 Participants

The study was approved by the University of Pittsburgh Institutional Review Board. Participants were recruited from the community, universities, and counseling and medical centers. They gave written informed consent, were right-handed, fluent in English, and were matched on age and sex. Individuals with unipolar depression (UD) met DSM-5 criteria for major depressive or persistent depressive disorders. Healthy controls (HC) had no personal or family history of psychiatric disorders. Exclusion criteria included a history of head injury, metal in the body, pregnancy, claustrophobia, neurodevelopmental disorders, systemic medical illness, premorbid IQ<85 per the National Adult Reading Test (36), current alcohol/drug abuse, Young Mania Rating Scale scores>10 (YMRS(37)) at scan, or meeting criteria for any psychotic-spectrum disorder. Data were drawn from an ongoing longitudinal study that includes neuroimaging sessions at baseline and 6-month follow-up and clinical evaluations at baseline, 6-months, and 12-months. The present report includes baseline data available for 55 HC and 50 UD. Participants were excluded from analyses due to (1) previously undetected brain abnormalities of potential clinical relevance: 2 UD, (2) diagnosis conversion during the course of the study: 1 HC was diagnosed with major depressive disorder, and 1 UD was diagnosed with bipolar disorder; (3) scanner or movement-related artifacts in MRI data (4 HC, 7 UD), and (4) poor-quality myelin maps (see Section 2.4.2.1 *Subject-level processing;* 3 HC, 1 UD). The final sample included 47 HC and 39 UD.

### 2.2. Clinical assessment

All diagnoses were made by a trained clinician and confirmed by a psychiatrist according to DSM-5 criteria using SCID-5 (38). Additional information collected included illness onset and duration, number of current episodes, comorbid psychiatric disorders, current depression symptoms using the Hamilton Rating Scale for Depression (HRSD-25) (39), current mania symptoms using the Young Mania Rating Scale (YMRS) (37), and lifetime depression and hypo/mania spectrum symptomatology using the Mood Spectrum Self Report (MOODS-SR) (40). A total psychotropic medication load was calculated for each participant with UD, with greater numbers and doses of medications corresponding to a greater medication load (41,42).

### 2.3 Neuroimaging data acquisition

The neuroimaging data were collected at the University of Pittsburgh/UPMC Magnetic Resonance Research Center using a 3T Siemens Prisma scanner with a 64-channel receiver head coil and named according to the ReproIn convention (43). The DICOM images were converted to BIDS dataset using *heudiconv* (44) and *dcm2niix* (45). High-resolution T1w images were collected using the MPRAGE sequence with TR=2400ms, resolution=0.8×0.8×0.8mm, 208 slices, FOV=256, TE=2.22ms, flip angle=8°. High-resolution T2w images were collected using TR=3200ms, resolution=0.8×0.8×0.8mm, 208 slices, FOV=256, TE=563ms. Field maps were collected in the AP and PA directions using the spin echo sequence (TR=8000, resolution=2×2×2mm, FOV=210, TE=66ms, flip angle=90°, 72 slices).

### 2.4 Data analyses

#### 2.4.1 Clinical data analysis

HC and UD groups were compared on demographic and clinical variables using t-tests and chi-square tests. All analyses were conducted in R ([https://www.r-project.org/](https://www.r-project.org/)).

#### 2.4.2 Neuroimaging data processing

##### 2.4.2.1 Subject-level preprocessing

Data quality was examined using *mriqc version 0*.*15*.*1* (46) and visually inspected (Supplemental Methods). Each participant’s cortical myelin was characterized with the T1w/T2w ratio (30,47,48) using the PreFreeSurfer, FreeSurfer, and PostFreeSurfer minimal preprocessing pipelines for the human connectome project (47). *Workbench v1*.*4*.*2* and *HCPpipelines-4*.*1*.*3* were installed system-wide on a workstation with GNU/Linux Debian 10 operating system. The spin echo field maps collected in AP and PA phase encoding directions were used for bias field correction in *PreFreeSurfer*. Registration to standard space was achieved via MSMSulc (49) in *PostFreeSurfer*. If FreeSurfer images and myelin maps had artifacts and gross errors (e.g., large regions of apparent low myelin in the occipital cortex due to the transverse sinus interfering with accurate identification of the pial surface), the data were removed from analyses. The resulting myelin maps were parcellated in *Workbench* using the 360 region Glasser Atlas (48). Further data quality assurance resulted in the removal of 11 outlier parcels (Supplemental Methods), leaving 349 parcels.

#### 2.4.3 Neuroimaging data analysis

##### 2.4.3.1 Elastic net and linear discriminant analysis

Neuroimaging data decoding studies can capitalize on complex relationships between variables, but their large numbers can present a challenge in deriving a generalizable model. The elastic net approach has emerged as a flexible tool for use with neuroimaging data (50) as it is able to reduce the influence of overly large coefficients and reducing the number of variables while generating multivariate models predictive of complex behaviors (51–53). Elastic net is a regularized regression which combines lasso and ridge regression (i.e., L1- and L2-norm regularization) (54). Ridge regression penalizes overly large coefficients, while lasso regression removes variables with small coefficients.

We used logistic elastic net regularized regression (55) to select variables (brain parcels) that were most predictive of case/control status. Elastic net has two parameters: alpha (*α)* controls the balance between the ridge and lasso regularizations, and lambda (*λ)* controls the strength of regularization. To provide equal contribution of each penalty to the loss function, we used *α*=0.5. To avoid model overfitting and bias, we implemented nested cross-validation to identify the optimal *λ* parameter. A linear discriminant analysis (LDA) (56) model was subsequently trained using selected variables to make out-of-sample prediction on held-out participants. This strategy is illustrated in Figure 1A and described in detail in Supplemental Methods. For each repetition of the nested cross-validation loop, two participants (1 UD and 1 HC) were held out. The rest of the sample was used to identify the optimal *λ* parameter which were then used to fit the elastic net model and select variables whose myelin levels were predictive of UD/HC status. These variables were then used to train an LDA model, which was tested on held-out participants. Results (model fit, variable selection, prediction accuracy) were evaluated against results from a permuted-label analysis (Figure 1B and Supplemental Methods).

![Figure 1.](http://medrxiv.org/http://medrxiv.stage.highwire.org/content/medrxiv/early/2021/03/12/2021.02.25.21252472/F1.medium.gif)

[Figure 1.](http://medrxiv.org/content/early/2021/03/12/2021.02.25.21252472/F1)

Figure 1. 
Diagram of analysis steps. Conceptual depiction of analysis steps including: (1) a unique pair of one UD and one HC participant is held-out; (2) an elastic net regression is used to select variables; (3) the retained variables are used an LDA model predicting case/control status; (4) the LDA model is tested on the held-out sample; (5) this process is repeated for each of the n=1833 pairs of subjects; (6) for each held-out pair, the training procedure is repeated with 100 unique permutations.

#### 2.4.3.2 Post-hoc analyses

##### 2.4.3.2.1 The relationship between true-label and permuted-label sample demographics

To ensure that the permuted results were not due to changes in the internal structure of the permuted samples in terms of demographic variables (i.e., age, sex and IQ), we compared the age, sex and IQ values in the permuted samples with that in the true-label samples.

##### 2.4.3.2.2 Association of LDA accuracy and cortical myelin with demographic and clinical characteristics

To further characterize the parcels selected by the logistic elastic net regression, we compared cortical myelin in UD vs. HC while controlling for age, sex, and IQ. To assess the potential influence of confounding variables on model fit, we tested whether demographic or clinical characteristics were predictive of classification accuracy. Regressions tested whether the joint effect of group, the variable in question, or their interaction, was associated with participant-wise accuracy (57). Additional models tested whether clinical characteristics of UD participants were associated with model accuracy. To assess the influence of confounding variables on variable selection, regressions similarly tested the association of demographic and clinical variables with cortical myelin, controlling for age, sex, and IQ. Results for each variable were separately corrected for multiple comparisons using false discovery rate (FDR).

#### 2.4.3.3 Exploratory analysis of a HC participant who was diagnosed with major depressive disorder 12 months from the baseline scan

One participant entered the study as a HC but was diagnosed with major depressive disorder sometime between 6 and 12 months after study onset. At the study visit at 12 months the participant had mild depressive symptoms. While this participant was excluded from all primary analyses described above, exploratory analyses investigated whether the myelin was predictive of the participant’s conversion from HC to UD. This analysis used the primary 86 participants and the variables selected in primary analyses (cortical myelin in 33 parcels and IQ) to train an LDA model. UD/HC status was then predicted at study onset and at the 6-month follow-up (both time points were prior to conversion).

## 3 Results

### 3.1 Sample demographics

Individuals with UD did not differ by age or sex but had higher IQ and current and lifetime depression severity compared to HC (Table 1).

View this table:
[Table 1.](http://medrxiv.org/content/early/2021/03/12/2021.02.25.21252472/T1)

Table 1. 
Demographic and clinical characteristics

### 3.2 UD vs. HC LDA nested cross-validation classification accuracy

Cortical myelin levels and IQ distinguished UD from HC in subjects held-out during cross-validation, with an average accuracy of 68% (Figure 2A; sensitivity (UD): 63.8%, specificity (HC): 71.5%). The mean participant-wise accuracy across 84 classification loops ranged from 0% to 100%, and 84% of participants were classified either quite accurately (>80% n=53) or inaccurately (<20% n=19). When demographic variables were excluded from the LDA classification, nested cross-validation achieved 69% accuracy (UD: 65.2%, HC: 73%). Notably, while the mean classification accuracy was higher in the HC group, this difference was not statistically significant (t=0.88, p=0.38).

![Figure 2.](http://medrxiv.org/http://medrxiv.stage.highwire.org/content/medrxiv/early/2021/03/12/2021.02.25.21252472/F2.medium.gif)

[Figure 2.](http://medrxiv.org/content/early/2021/03/12/2021.02.25.21252472/F2)

Figure 2. 
Classification accuracy and variable selection. A) Average subject-wise nested cross-validation classification accuracy in both healthy controls (HC) and participants with unipolar depression (UD). B) Variable selection frequency (percent of models where a given parcel was retained) in the true data (light blue) and permutations (green). The n=90 parcels which were retained at least once are shown. IQ, the only demographic variable selected, is labeled. Solid black lines represent the median of the variable selection frequency for all n=350 variables +/-3.5 x the interquartile range (IQR). See Supplemental Table 2 for the full list of variable selection frequencies.

In permutation analyses, no variable was selected in 66.9% of 183300 models. Within each participant, the proportion of models that did not select any variables ranged from 65.4% to 68.7%. Excluding instances when no variables were selected, the average participant-wise LDA accuracy in permutation analyses ranged from 41.2% to 59.2%, with an average accuracy of 50.5% (i.e., chance level). The test of whether age, sex and IQ in permuted samples differed from that in the sample with true labels showed that the demographics of the permuted groups differed (p<0.05, uncorrected) from the true-label groups only in 1.7% of cases.

View this table:
[Table 2.](http://medrxiv.org/content/early/2021/03/12/2021.02.25.21252472/T2)

Table 2. 
Association of unipolar depression with cortical myelin in selected parcels

### 3.3 Elastic net variable selection

*True-label* nested cross-validated elastic net models predicting diagnostic status (UD vs. HC) with 349 myelin parcels, age, sex and IQ, selected 90 myelin parcels (Supplemental Table 1) and IQ in at least one model (Figure 2B). Each model selected between 9 and 68 variables, with a median of 17 variables. *Permuted-label* nested cross-validated elastic net models selected all predictor variables at least once, but no variable was selected by more than 3% of the models. In addition, 66.9% of models with permuted labels did not select any variable at all. Within permuted-label models where at least one variable was selected, the number of selected variables ranged from 1 to 113, with a median of 17 variables selected.

Given that the variable selection frequency in the models with permuted labels likely represents noise, we applied the criterion of the median + 3.5*IQR across all permuted variables (3.77%) as the cutoff value to separate the potential noise variables from ‘signal’ variables in the nested cross-validated elastic net with true labels. This latter analysis identified 33 out of 90 myelin parcels, plus IQ, across the true-label nested cross-validated elastic net models that were above the cutoff line (Table 2, Figure 2B). These 33 parcels included multiple regions in the orbitofrontal cortex, insula, cingulate, and frontal operculum, as well as regions in the auditory and visual cortices (Supplemental Figure 3). No model retained the age and sex variables.

![Figure 3.](http://medrxiv.org/http://medrxiv.stage.highwire.org/content/medrxiv/early/2021/03/12/2021.02.25.21252472/F3.medium.gif)

[Figure 3.](http://medrxiv.org/content/early/2021/03/12/2021.02.25.21252472/F3)

Figure 3. 
Association of diagnostic status with cortical myelin in selected parcels. A) The t-statistic for the association of unipolar depression (UD) with cortical myelin in the n=33 selected parcels. Regions where UD participants had greater average myelin than HC participants are in red, and regions where HC participants had greater average myelin than UD participants are in blue. B) The percent of models in which each parcel was retained. Regions are individually labeled (note that FFC is not visible, as it lies on the ventral surface).

Post-hoc analyses tested the association of diagnostic status with cortical myelin levels in the 33 selected myelin parcels (Table 2, Figure 3). These parcels include both regions where UD participants have lower mean cortical myelin than HC and regions where UD participants have greater mean cortical myelin, as well as regions where the two groups do not differ in their average level of cortical myelin. After FDR-correction for multiple comparisons, fourteen of the parcels showed evidence of significant differences between UD and HC participants (Table 2). An additional 8 parcels showed nominally significant differences between UD and HC participants (Table 2; p<0.05 uncorrected). Parcels that showed a greater absolute mean difference between the groups were selected more frequently in the nested cross-validation analysis (*r*=0.7, *p*=4×10−6).

The 33 selected parcels included 13 parcels in regions important for executive functioning and cognitive control, located in the anterior cingulate (left a24pr, left 24dv, left p32, and left s32), orbitofrontal cortex (left 11l, left p10p, and left 47m), posterior cingulate (right RSC, and left d23ab), frontal operculum (right FOP1 and left FOP4), and parietal cortex (right POS2 and right 7Pm). UD was largely associated with reduced cortical myelin in these regions, except in one region of the anterior cingulate (left p32), and the two regions in the parietal cortex (right POS2 and right 7Pm), where it was associated with increased myelin. An additional 19 of the identified parcels play roles in visual, somatomotor, and auditory processing. Visual processing regions included extrastriate regions (right V3A, left V4t, and left MT), ventral stream regions (right FFC and right TE2p), and dorsal stream regions (left MST, right PH, left 7PC, right 7AL). Auditory processing regions included regions in the auditory cortex (left LBelt, right LBelt, left RI, and right RI) and regions implicated in language (left STSdp, left FOP2, and right 52). Somatomotor regions included regions implicated in somatosensation (right 5m and right OP1) and a region in the premotor cortex (left 6a). UD was largely associated with reduced cortical myelin in auditory and somatomotor regions, except for two auditory regions (right 52 and right RI) where it was associated with increased myelin. In contrast, participants with UD had greater myelin in 6 visual regions and had reduced myelin in 3 regions (left V4t, left FFC, and left 7AL).

### 3.5 Post-hoc analyses

#### 3.5.1 Association of cortical myelin with demographic and clinical variables

The association of cortical myelin in the 33 selected parcels with demographic and clinical variables (see Table 1) was not statistically significant (Supplemental Table 4). Across all participants neither sex, IQ, nor MOODs-SR score were correlated with cortical myelin in selected parcels, nor did they interact with group (HC/UD) to predict cortical myelin. There was evidence for associations between age and cortical myelin in two regions (p-fdr<0.05), driven by a positive correlation in right LBelt, and a negative correlation in left a24pr. In individuals with UD, no clinical variable was correlated with cortical myelin in selected parcels. However, taking antidepressants was nominally associated (p<0.05 uncorrected) with increased cortical myelin in right OP1, left 24dv, and left p32, and decreased cortical myelin in right LBelt (Supplemental Table 4, Supplemental Figure 3).

#### 3.5.2 Association of LDA accuracy with demographic and clinical variables

Analyses found that no demographic or clinical variable was predictive of classification accuracy (Table 3). Similarly, within the UD participant group, no clinical or medication variable was associated with classification accuracy (Table 3).

View this table:
[Table 3.](http://medrxiv.org/content/early/2021/03/12/2021.02.25.21252472/T3)

Table 3. 
Association of demographic and clinical variables with classification accuracy

#### 3.5.3 Exploratory analysis of classification in a participant who converted from HC to UD

The LDA trained on the whole sample of 86 participants, with IQ and the 33 parcels identified in the previous analyses as predictors classified this participant as ‘UD’ both times: 12 months and 6 months before illness onset.

## 4 DISCUSSION

Cortical myelin distinguished healthy control participants (HC) from individuals with unipolar depression (UD) with 68% accuracy. The elastic net regression selected regions implicated in executive function and reward processing, as well as those involved in visual, auditory, and somatomotor processing. The pattern of cortical myelin distinguishing UD from HC was associated with a myelin reduction in some regions and increased myelin in other regions. These results demonstrate that aberrant levels of cortical myelin may be a biomarker of UD.

The 33 selected parcels included 13 parcels in regions important for executive functioning, cognitive control, and reward processing, including the anterior and posterior cingulate, orbitofrontal cortex, frontal operculum, and parietal cortex. There is abundant evidence that these processes are disrupted in depression (5,11). The present results converge with mounting evidence of disrupted myelination (28), thickness (12), connectivity (58), and activation (11) of these regions in those with depression. In addition to regions traditionally reported in neuroimaging studies of depression, we found that parcels located along the dorsal and ventral visual pathways, somatosensory and auditory processing regions and regions implicated in language also contributed to UD vs. HC classification. While depression is not typically considered a disorder of dysfunctional sensory processing, some studies report visual (59) and auditory (60) processing deficits in depression. In addition, somatosensory symptoms (e.g., motor retardation) are among the diagnostic criteria for this disorder and have been well described (61). Our results add to a growing body of literature documenting associations of depression with altered structure (12) and functional connectivity (62,63) in sensory regions. Broadly, our results contribute to the consensus that the neurobiology of regions that play important roles in cognitive processes and information transfer are disrupted in depression, which may contribute to the etiology of the disorder.

It is notable that in contrast to studies of the major white matter tracts, where depression is associated with lower integrity (15), we observed both decreased and increased cortical myelin in individuals with UD relative to HC. Additionally, the anterior cingulate, the superior parietal cortex, the retroinsular cortex, and visual processing regions included a combination of proximal regions with myelin decreases and increases. While the relationship between cortical myelin and functional activation or connectivity remains under-explored, these observations suggest that cortical myelin imbalance, rather than a uniform reduction, could drive some of the observed functional differences in depression, such as disrupted network integration (63,64) and patterns of both hypo-and hyper-connectivity (65,66).

Intriguingly, while classification accuracy was bimodal, in permutations participants were classified at chance (50%). This suggests that participants who were consistently misclassified differed from those who were correctly classified in their demographic or clinical features. However, *post-hoc* analyses did not reveal any associations between classification accuracy and demographic or clinical variables, thus suggesting that misclassification could reflect the presence of environmental or genetic risk factors that are not captured by clinical measures. These factors could alter the myelin content in selected regions thus creating predisposition to depression at the brain level. As a preliminary exploration of this idea, we explored the myelin patterns in a participant who experienced the onset of UD during the course of the study and was not used in primary analyses. This participant completed two MRI scans: at baseline and 6 months follow-up with both time points occurring prior to UD onset. Remarkably, this participant was classified as UD on both scans, despite not yet meeting criteria for a UD diagnosis at the time of scan. This preliminary result suggests that the pattern of cortical myelin in frontal, sensorimotor and extended visual cortices may be a biological risk marker predictive of UD diagnosis in the future. Further longitudinal studies are needed to test this hypothesis.

While the present results demonstrate that the pattern of cortical myelin is disrupted in UD, the cause of myelin disruption remains unknown. *Post-hoc* analyses suggest that myelin disruption is not related to medication use, lifetime depression severity, or illness duration. However, disrupted cortical myelin may reflect other risk factors for UD. For instance, sleep disturbance, which is a well-established risk factor and symptom of UD (67), was recently shown to correlate with cortical myelin in several of the same regions found in the present study, including the cingulate and middle temporal cortex (32). Stress is also associated with reduced white matter integrity in depression (15,68,69). After stress, remyelination can occur (70) thus resulting into altered patterns of myelination across the cortex (71). This ‘remyelination’ hypothesis could potentially explain our observations of both decreased and increased cortical myelin in unipolar depression.

Limitations of this work includes the need to replicate our findings in an independent sample. To partially address this limitation, our analyses used robust machine learning methods involving model testing using held-out samples. The major strength of this approach is that it helps to reduce model bias, which occurs when the same participants are used to train and test a model (72). It has been suggested that model performance with nested cross-validation is close to the accuracy that would be achieved on fully independent data (73). The second limitation concerns measurement noise due to susceptibility artifacts. Several regions with documented relevance to depression, including the bilateral hippocampus, entorhinal cortex, and posterior orbitofrontal cortex complex, were not included in the present analyses, as these regions showed an excess of between-person variability (see Supplemental Methods*)*. Future research should explore the association of cortical myelin in these regions with unipolar depression.

In summary, cortical myelin can distinguish participants with UD from HC, even when clinical and demographic variables are not included in analyses. Regions that were most important for this classification include several that play key roles in reward and emotion processing as well as a host of regions important for sensory processing. This result highlights that the association of UD with sensory processing bears further investigation. Notably, UD was associated with both decreased and increased cortical myelin, suggesting that observations of reduced integrity of major white matter tracts in UD may not fully extend to the cortex. These results suggest that cortical myelin holds promise as a biomarker of unipolar depression and may be an early predictor of risk for this disorder.

## Supporting information

Primary Supplement [[supplements/252472_file02.pdf]](pending:yes)

Supplemental Tables [[supplements/252472_file03.xlsx]](pending:yes)

## Data Availability

Data will be made publicly available upon publication.

## AUTHORS CONTRIBUTION

D.B. – data quality assurance, analysis, visualization, and interpretation; drafted and critically evaluated the manuscript

Y.O.H. – curated data organization and analyses, drafted and critically evaluated the manuscript

S.S., R.R. – data acquisition and quality assurance; drafted and critically evaluated the manuscript

S.I. - curated data analyses, interpreted the data, critically evaluated the manuscript H.A.S., M.L.P. – curated participants’ recruitment, interpreted the data, critically evaluated the manuscript.

A.M. – obtained funding; designed the study; data acquisition, quality assurance, analysis and interpretation; drafted and critically evaluated the manuscript

All authors have approved the final version of the manuscript and agreed to be accountable for all aspects of this work.

## ACKNOWLEDGMENTS

Funding Acknowledgements: This work was supported by a grant from the National Institute of Health R01MH114870 to A.M., and Y.O.H was supported by P41EB019936 to the Center for Reproducible Neuroimaging Computation (PI: Kennedy). D.A.A.B, Y.O.H., S.S., R.R, S.I., and A.M.: declare no conflict of interest. H.A.S: receives royalties from Wolters Kluwer, royalties and an editorial stipend from APA Press, and honorarium from Novus Medical Education.

The authors thank participants for taking part in this research study. We also thank Dr. Mary L. Phillips for fruitful discussions of the study design.

*   Received February 25, 2021.
*   Revision received March 11, 2021.
*   Accepted March 12, 2021.


*   © 2021, Posted by Cold Spring Harbor Laboratory

This pre-print is available under a Creative Commons License (Attribution-NonCommercial-NoDerivs 4.0 International), CC BY-NC-ND 4.0, as described at [http://creativecommons.org/licenses/by-nc-nd/4.0/](http://creativecommons.org/licenses/by-nc-nd/4.0/)

## References

1.  1.World Health Organisation. Depression and other common mental disorders: global health estimates. World Health Organization. 2017;
    
    

2.  2.Greenberg PE, Fournier AA, Sisitsky T, Pike CT, Kessler RC. The economic burden of adults with major depressive disorder in the United States (2005 and 2010). Journal of Clinical Psychiatry. 2015;
    
    

3.  3.Cuijpers P, Stringaris A, Wolpert M. Treatment outcomes for depression: challenges and opportunities. The Lancet Psychiatry. 2020;7(11):925–7.
    
    

4.  4.Gotlib IH, Joormann J. Cognition and depression: Current status and future directions. Annual Review of Clinical Psychology. 2010.
    
    

5.  5.Snyder HR. Major depressive disorder is associated with broad impairments on neuropsychological measures of executive function: A meta-analysis and review. Psychological Bulletin. 2013;
    
    

6.  6.Suslow T, Hußlack A, Kersting A, Bodenschatz CM. Attentional biases to emotional information in clinical depression: A systematic and meta-analytic review of eye tracking findings. Journal of Affective Disorders. 2020.
    
    

7.  7.Halahakoon DC, Kieslich K, O’Driscoll C, Nair A, Lewis G, Roiser JP. Reward-Processing Behavior in Depressed Participants Relative to Healthy Volunteers: A Systematic Review and Meta-analysis. JAMA Psychiatry. 2020.
    
    

8.  8.Stuhrmann A, Suslow T, Dannlowski U. Facial emotion processing in major depression: A systematic review of neuroimaging findings. Biology of Mood and Anxiety Disorders. 2011;
    
    

9.  9.Hamilton JP, Etkin A, Furman DJ, Lemus MG, Johnson RF, Gotlib IH. Functional neuroimaging of major depressive disorder: A meta-analysis and new integration of baseline activation and neural response data. American Journal of Psychiatry. 2012.
    
    

10. 10.Wang XL, D. MY, Chen TL, Chen ZQ, Huang XQ, Luo Y, et al. Neural correlates during working memory processing in major depressive disorder. Progress in Neuro-Psychopharmacology and Biological Psychiatry. 2015;
    
    

11. 11.Ng TH, Alloy LB, Smith D V. Meta-analysis of reward processing in major depressive disorder reveals distinct abnormalities within the reward circuit. Translational Psychiatry. 2019;
    
    

12. 12.Schmaal L, Hibar DP, Sämann PG, Hall GB, Baune BT, Jahanshad N, et al. Cortical abnormalities in adults and adolescents with major depression based on brain scans from 20 cohorts worldwide in the ENIGMA Major Depressive Disorder Working Group. Molecular Psychiatry. 2017;
    
    

13. 13.Schmaal L, Veltman DJ,  Van Erp Tgm, Smann PG, Frodl T, Jahanshad N, et al. Subcortical brain alterations in major depressive disorder: Findings from the ENIGMA Major Depressive Disorder working group. Molecular Psychiatry. 2016;
    
    

14. 14.Ho TC, Gutman B, Pozzi E, Grabe HJ, Hosten N, Wittfeld K, et al. Subcortical shape alterations in major depressive disorder: Findings from the ENIGMA major depressive disorder working group. Human Brain Mapping. 2020;
    
    

15. 15.van Velzen LS, Kelly S, Isaev D, Aleman A, Aftanas LI, Bauer J, et al. White matter disturbances in major depressive disorder: a coordinated analysis across 20 international cohorts in the ENIGMA MDD working group. Molecular Psychiatry. 2020;
    
    

16. 16.Shen X, Reus LM, Cox SR, Adams MJ, Liewald DC, Bastin ME, et al. Subcortical volume and white matter integrity abnormalities in major depressive disorder: Findings from UK Biobank imaging data. Scientific Reports. 2017;
    
    

17. 17.Murphy ML, Frodl T. Meta-analysis of diffusion tensor imaging studies shows altered fractional anisotropy occurring in distinct brain areas in association with depression. Biology of Mood & Anxiety Disorders. 2011;
    
    

18. 18.Liao Y, Huang X, Wu Q, Yang C, Kuang W, Du M, et al. Is depression a disconnection syndrome? Meta-analysis of diffusion tensor imaging studies in patients with MDD. Journal of Psychiatry and Neuroscience. 2013;
    
    

19. 19.Wise T, Radua J, Nortje G, Cleare AJ, Young AH, Arnone D. Voxel-based meta-Analytical evidence of structural disconnectivity in major depression and bipolar disorder. Biological Psychiatry. 2016;
    
    

20. 20.Penke L, Maniega SM, Murray C, Gow AJ, Valdés Hernández MC, Clayden JD, et al. A general factor of brain white matter integrity predicts information processing speed in healthy older people. Journal of Neuroscience. 2010;
    
    

21. 21.Chopra S, Shaw M, Shaw T, Sachdev PS, Anstey KJ, Cherbuin N. More highly myelinated white matter tracts are associated with faster processing speed in healthy adults. NeuroImage. 2018;
    
    

22. 22.Eden AS, Schreiber J, Anwander A, Keuper K, Laeger I, Zwanzger P, et al. Emotion regulation and trait anxiety are predicted by the microstructure of fibers between amygdala and prefrontal cortex. Journal of Neuroscience. 2015;
    
    

23. 23.Welton T, Indja BE, Maller JJ, Fanning JP, Vallely MP, Grieve SM. Replicable brain signatures of emotional bias and memory based on diffusion kurtosis imaging of white matter tracts. Human Brain Mapping. 2020;
    
    

24. 24.de Boer L, Garzón B, Axelsson J, Riklund K, Nyberg L, Bäckman L, et al. Corticostriatal white matter integrity and dopamine D1 receptor availability predict age differences in prefrontal value signaling during reward learning. Cerebral Cortex. 2020;
    
    

25. 25.Simons M, Nave KA. Oligodendrocytes: Myelination and axonal support. Cold Spring Harbor Perspectives in Biology. 2016;8(1):1–16.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1101/cshperspect.a022053&link_type=DOI) 

26. 26.Boda E. Myelin and oligodendrocyte lineage cell dysfunctions: New players in the etiology and treatment of depression and stress-related disorders. European Journal of Neuroscience. 2019.
    
    

27. 27.Tham MW, Woon PS, Sum MY, Lee TS, Sim K. White matter abnormalities in major depression: Evidence from post-mortem, neuroimaging and genetic studies. Journal of Affective Disorders. 2011.
    
    

28. 28.Zhang TJ, Wu QZ, Huang XQ, Sun XL, Zou K, Lui S, et al. Magnetization transfer imaging reveals the brain deficit in patients with treatment-refractory depression. Journal of Affective Disorders. 2009;
    
    

29. 29.Sacchet MD, Gotlib IH. Myelination of the brain in major depressive disorder: An in vivo quantitative magnetic resonance imaging study. Scientific Reports. 2017;7(1).
    
    

30. 30.Glasser MF, van Essen DC. Mapping human cortical areas in vivo based on myelin content as revealed by T1-and T2-weighted MRI. Journal of Neuroscience. 2011;31(32):11597–616.
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6Njoiam5ldXJvIjtzOjU6InJlc2lkIjtzOjExOiIzMS8zMi8xMTU5NyI7czo0OiJhdG9tIjtzOjUwOiIvbWVkcnhpdi9lYXJseS8yMDIxLzAzLzEyLzIwMjEuMDIuMjUuMjEyNTI0NzIuYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) 

31. 31.Shafee R, Buckner RL, Fischl B. Gray matter myelination of 1555 human brains using partial volume corrected MRI images. NeuroImage. 2015;105:473–85.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.neuroimage.2014.10.054&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=25449739&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F03%2F12%2F2021.02.25.21252472.atom) 

32. 32.Toschi N, Passamonti L, Bellesi M. Sleep quality relates to emotional reactivity via intracortical myelination. Sleep. 2020;
    
    

33. 33.Toschi N, Passamonti L. Intra-cortical myelin mediates personality differences. Journal of Personality. 2019;
    
    

34. 34.Norbom LB, Doan NT, Alnæs D, Kaufmann T, Moberget T, Rokicki J, et al. Probing Brain Developmental Patterns of Myelination and Associations With Psychopathology in Youths Using Gray/White Matter Contrast. Biological Psychiatry. 2019;
    
    

35. 35.Dalili MN, Penton-Voak IS, Harmer CJ, Munafò MR. Meta-analysis of emotion recognition deficits in major depressive disorder. Psychological Medicine. 2015.
    
    

36. 36.Nelson HE. National Adult Reading Test (NART): For the Assessment of Premorbid Intelligence in Patients with Dementia: Test Manual. 1982. Nfer-Nelson; 1982.
    
    

37. 37.Young RC, Biggs JT, Ziegler VE, Meyer DA. A rating scale for mania: Reliability, validity and sensitivity. British Journal of Psychiatry. 1978;133:429–35.
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6MTA6ImJqcHJjcHN5Y2giO3M6NToicmVzaWQiO3M6OToiMTMzLzUvNDI5IjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjEvMDMvMTIvMjAyMS4wMi4yNS4yMTI1MjQ3Mi5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 

38. 38.First MB. Structured Clinical Interview for the DSM (SCID). In: The Encyclopedia of Clinical Psychology. 2015.
    
    

39. 39.Hamilton M. A rating scale for depression. Journal of neurology, neurosurgery, and psychiatry. 1960;23:56–62.
    
    [FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6MzoiUERGIjtzOjExOiJqb3VybmFsQ29kZSI7czo0OiJqbm5wIjtzOjU6InJlc2lkIjtzOjc6IjIzLzEvNTYiO3M6NDoiYXRvbSI7czo1MDoiL21lZHJ4aXYvZWFybHkvMjAyMS8wMy8xMi8yMDIxLjAyLjI1LjIxMjUyNDcyLmF0b20iO31zOjg6ImZyYWdtZW50IjtzOjA6IiI7fQ==) 

40. 40.Dell’Osso L, Armani A, Rucci P, Frank E, Fagiolini A, Corretti G, et al. Measuring mood spectrum: Comparison of interview (SCI-MOODS) and self-report (MOODS-SR) instruments. Comprehensive Psychiatry. 2002;43:69–73.
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=11788923&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F03%2F12%2F2021.02.25.21252472.atom) 

41. 41.Hassel S, Almeida JR, Kerr N, Nau S, Ladouceur CD, Fissell K, et al. Elevated striatal and decreased dorsolateral prefrontal cortical activity in response to emotional stimuli in euthymic bipolar disorder: no associations with psychotropic medication load. Bipolar disorders. 2008 Dec;10(8):916–27.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1111/j.1399-5618.2008.00641.x&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=19594507&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F03%2F12%2F2021.02.25.21252472.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000261054400008&link_type=ISI) 

42. 42.Manelis A, Almeida JRC, Stiffler R, Lockovich JC, Aslam HA, Phillips ML. Anticipation-related brain connectivity in bipolar and unipolar depression: A graph theory approach. Brain. 2016;139(9):2554–66.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/brain/aww157&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=27368345&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F03%2F12%2F2021.02.25.21252472.atom) 

43. 43.Visconti di Oleggio Castello M, Dobson JE, Sackett T, Kodiweera C, Haxby J V., Goncalves M, et al. ReproNim/reproin 0.6.0. 2020 Jan 22;
    
    

44. 44.Halchenko Y, Goncalves M,  Castello MV di O, Ghosh S, Hanke M, Dae, et al. nipy/heudiconv: v0.5.4 [0.5.4] - 2019-04-29. 2019 Apr 29;
    
    

45. 45.Li X, Morgan PS, Ashburner J, Smith J, Rorden C. The first step for neuroimaging data analysis: DICOM to NIfTI conversion. Journal of Neuroscience Methods. 2016 May 1;264:47–56.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.jneumeth.2016.03.001&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=26945974&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F03%2F12%2F2021.02.25.21252472.atom) 

46. 46.Esteban O, Birman D, Schaer M, Koyejo OO, Poldrack RA, Gorgolewski KJ. MRIQC: Advancing the automatic prediction of image quality in MRI from unseen sites. PLoS ONE. 2017;
    
    

47. 47.Glasser MF, Sotiropoulos SN, Wilson JA, Coalson TS, Fischl B, Andersson JL, et al. The minimal preprocessing pipelines for the Human Connectome Project. NeuroImage. 2013;
    
    

48. 48.Glasser MF, Coalson TS, Robinson EC, Hacker CD, Harwell J, Yacoub E, et al. A multi-modal parcellation of human cerebral cortex. Nature. 2016;
    
    

49. 49.Robinson EC, Garcia K, Glasser MF, Chen Z, Coalson TS, Makropoulos A, et al. Multimodal surface matching with higher-order smoothness constraints. NeuroImage. 2018;
    
    

50. 50.Jollans L, Boyle R, Artiges E, Banaschewski T, Desrivières S, Grigis A, et al. Quantifying performance of machine learning methods for neuroimaging data. NeuroImage. 2019;
    
    

51. 51.Manelis A, Iyengar S, Swartz HA, Phillips ML. Prefrontal cortical activation during working memory task anticipation contributes to discrimination between bipolar and unipolar depression. Neuropsychopharmacology. 2020;
    
    

52. 52.Acuff HE, Versace A, Bertocci MA, Ladouceur CD, Hanford LC, Manelis A, et al. Baseline and follow-up activity and functional connectivity in reward neural circuitries in offspring at risk for bipolar disorder. Neuropsychopharmacology. 2019;
    
    

53. 53.Bertocci MA, Bebko G, Versace A, Fournier JC, Iyengar S, Olino T, et al. Predicting clinical outcome from reward circuitry function and white matter structure in behaviorally and emotionally dysregulated youth. Molecular Psychiatry. 2016;
    
    

54. 54.Zou H, Hastie T. Regularization and variable selection via the elastic net. Journal of the Royal Statistical Society Series B: Statistical Methodology. 2005;
    
    

55. 55.Friedman J, Hastie T, Tibshirani R. Regularization paths for generalized linear models via coordinate descent. Journal of Statistical Software. 2010;
    
    

56. 56.Venables WN, Ripley BD. Modern Applied Statistics with S. Fourth edi. Springer; 2002.
    
    

57. 57.Dinga R, Schmaal L,  Penninx Bwjh, Veltman DJ, Marquand AF. Controlling for Effects of Confounding Variables on Machine Learning Predictions. bioRxiv. 2020;
    
    

58. 58.Weiskopf N, Hutton C, Josephs O, Deichmann R. Optimal EPI parameters for reduction of susceptibility-induced BOLD sensitivity losses: A whole-brain analysis at 3 T and 1.5 T. NeuroImage. 2006;33(2):493–504.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.neuroimage.2006.07.029&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=16959495&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F03%2F12%2F2021.02.25.21252472.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000241406800007&link_type=ISI) 

59. 59.Grange JA, Rydon-Grange M. Computational modelling of attentional selectivity in depression reveals perceptual deficits. Psychological Medicine. 2020;
    
    

60. 60.Zweerings J, Zvyagintsev M, Turetsky BI, Klasen M, König AA, Roecher E, et al. Fronto‐parietal and temporal brain dysfunction in depression: A fMRI investigation of auditory mismatch processing. Human Brain Mapping. 2019 May 12;40(12):hbm.24623.
    
    

61. 61.Simon GE, VonKorff M, Piccinelli M, Fullerton C, Ormel J. An International Study of the Relation between Somatic Symptoms and Depression. New England Journal of Medicine. 1999 Oct 28;341(18):1329–35.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1056/NEJM199910283411801&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=10536124&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F03%2F12%2F2021.02.25.21252472.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000083357800001&link_type=ISI) 

62. 62.Korgaonkar MS, Goldstein-Piekarski AN, Fornito A, Williams LM. Intrinsic connectomes are a predictive biomarker of remission in major depressive disorder. Molecular Psychiatry. 2020 Jul 1;25(7):1537–49.
    
    

63. 63.Chen H, Liu K, Zhang B, Zhang J, Xue X, Lin Y, et al. More optimal but less regulated dorsal and ventral visual networks in patients with major depressive disorder. Journal of Psychiatric Research. 2019 Mar 1;110:172–8.
    
    

64. 64.Zheng H, Li F, Bo Q, Li X, Yao L, Yao Z, et al. The dynamic characteristics of the anterior cingulate cortex in resting-state fMRI of patients with depression. Journal of Affective Disorders. 2018 Feb 1;227:391–7.
    
    

65. 65.Yu M, Linn KA, Shinohara RT, Oathes DJ, Cook PA, Duprat R, et al. Childhood trauma history is linked to abnormal brain connectivity in major depression. Proceedings of the National Academy of Sciences of the United States of America. 2019 Apr 23;116(17):8582–90.
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NDoicG5hcyI7czo1OiJyZXNpZCI7czoxMToiMTE2LzE3Lzg1ODIiO3M6NDoiYXRvbSI7czo1MDoiL21lZHJ4aXYvZWFybHkvMjAyMS8wMy8xMi8yMDIxLjAyLjI1LjIxMjUyNDcyLmF0b20iO31zOjg6ImZyYWdtZW50IjtzOjA6IiI7fQ==) 

66. 66.Kaiser RH, Andrews-Hanna JR, Wager TD, Pizzagalli DA. Large-scale network dysfunction in major depressive disorder: A meta-analysis of resting-state functional connectivity. JAMA Psychiatry. 2015 Jun 1;72(6):603–11.
    
    

67. 67.Franzen PL, Buysse DJ. Sleep disturbances and depression: Risk relationships for subsequent depression and therapeutic implications. Dialogues in Clinical Neuroscience. 2008;10(4):473–81.
    
    

68. 68.Ziegler G, Moutoussis M, Hauser TU, Fearon P, Bullmore ET, Goodyer IM, et al. Childhood socio-economic disadvantage predicts reduced myelin growth across adolescence and young adulthood. Human Brain Mapping. 2020;
    
    

69. 69.Ho TC, King LS, Leong JK, Colich NL, Humphreys KL, Ordaz SJ, et al. Effects of sensitivity to life stress on uncinate fasciculus segments in early adolescence. Social Cognitive and Affective Neuroscience. 2017;
    
    

70. 70.Bonnefil V, Dietz K, Amatruda M, Wentling M, Aubry A V., Dupree JL, et al. Region-specific myelin differences define behavioral consequences of chronic social defeat stress in mice. eLife. 2019 Aug 1;8.
    
    

71. 71.Orthmann-Murphy J, Call CL, Molina-Castro GC, Hsieh YC, Rasband MN, Calabresi PA, et al. Remyelination alters the pattern of myelin in the cerebral cortex. eLife. 2020 May 1;9:1–61.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.7554/eLife.52760&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=32338598&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F03%2F12%2F2021.02.25.21252472.atom) 

72. 72.Davatzikos C. Machine learning in neuroimaging: Progress and challenges. Vol. 197, NeuroImage. Academic Press Inc.; 2019. p. 652–6.
    
    

73. 73.Varma S, Simon R. Bias in error estimation when using cross-validation for model selection. BMC Bioinformatics. 2006 Feb 23;7(1):1–8.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1186/1471-2105-7-1&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=16393334&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F03%2F12%2F2021.02.25.21252472.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000236027800001&link_type=ISI)