Trajectories of remitted psychotic depression: identification of predictors of worsening by machine learning.

BACKGROUND: Remitted psychotic depression (MDDPsy) has heterogeneity of outcome. The study's aims were to identify subgroups of persons with remitted MDDPsy with distinct trajectories of depression severity during continuation treatment and to detect predictors of membership to the worsening trajectory. METHOD: One hundred and twenty-six persons aged 18-85 years participated in a 36-week randomized placebo-controlled trial (RCT) that examined the clinical effects of continuing olanzapine once an episode of MDDPsy had remitted with sertraline plus olanzapine. Latent class mixed modeling was used to identify subgroups of participants with distinct trajectories of depression severity during the RCT. Machine learning was used to predict membership to the trajectories based on participant pre-trajectory characteristics. RESULTS: Seventy-one (56.3%) participants belonged to a subgroup with a stable trajectory of depression scores and 55 (43.7%) belonged to a subgroup with a worsening trajectory. A random forest model with high prediction accuracy (AUC of 0.812) found that the strongest predictors of membership to the worsening subgroup were residual depression symptoms at onset of remission, followed by anxiety score at RCT baseline and age of onset of the first lifetime depressive episode. In a logistic regression model that examined depression score at onset of remission as the only predictor variable, the AUC (0.778) was close to that of the machine learning model. CONCLUSIONS: Residual depression at onset of remission has high accuracy in predicting membership to worsening outcome of remitted MDDPsy. Research is needed to determine how best to optimize the outcome of psychotic MDDPsy with residual symptoms.

Expansion of a frontostriatal salience network in individuals with depression.

Hundreds of neuroimaging studies spanning two decades have revealed differences in brain structure and functional connectivity in depression, but with modest effect sizes, complicating efforts to derive mechanistic pathophysiologic insights or develop biomarkers. 1 Furthermore, although depression is a fundamentally episodic condition, few neuroimaging studies have taken a longitudinal approach, which is critical for understanding cause and effect and delineating mechanisms that drive mood state transitions over time. The emerging field of precision functional mapping using densely-sampled longitudinal neuroimaging data has revealed unexpected, functionally meaningful individual differences in brain network topology in healthy individuals, 2-5 but these approaches have never been applied to individuals with depression. Here, using precision functional mapping techniques and 11 datasets comprising n=187 repeatedly sampled individuals and >21,000 minutes of fMRI data, we show that the frontostriatal salience network is expanded two-fold in most individuals with depression. This effect was replicable in multiple samples, including large-scale, group-average data (N=1,231 subjects), and caused primarily by network border shifts affecting specific functional systems, with three distinct modes of encroachment occurring in different individuals. Salience network expansion was unexpectedly stable over time, unaffected by changes in mood state, and detectable in children before the subsequent onset of depressive symptoms in adolescence. Longitudinal analyses of individuals scanned up to 62 times over 1.5 years identified connectivity changes in specific frontostriatal circuits that tracked fluctuations in specific symptom domains and predicted future anhedonia symptoms before they emerged. Together, these findings identify a stable trait-like brain network topology that may confer risk for depression and mood-state dependent connectivity changes in frontostriatal circuits that predict the emergence and remission of depressive symptoms over time.

A somato-cognitive action network alternates with effector regions in motor cortex.

Motor cortex (M1) has been thought to form a continuous somatotopic homunculus extending down the precentral gyrus from foot to face representations1,2, despite evidence for concentric functional zones3 and maps of complex actions4. Here, using precision functional magnetic resonance imaging (fMRI) methods, we find that the classic homunculus is interrupted by regions with distinct connectivity, structure and function, alternating with effector-specific (foot, hand and mouth) areas. These inter-effector regions exhibit decreased cortical thickness and strong functional connectivity to each other, as well as to the cingulo-opercular network (CON), critical for action5 and physiological control6, arousal7, errors8 and pain9. This interdigitation of action control-linked and motor effector regions was verified in the three largest fMRI datasets. Macaque and pediatric (newborn, infant and child) precision fMRI suggested cross-species homologues and developmental precursors of the inter-effector system. A battery of motor and action fMRI tasks documented concentric effector somatotopies, separated by the CON-linked inter-effector regions. The inter-effectors lacked movement specificity and co-activated during action planning (coordination of hands and feet) and axial body movement (such as of the abdomen or eyebrows). These results, together with previous studies demonstrating stimulation-evoked complex actions4 and connectivity to internal organs10 such as the adrenal medulla, suggest that M1 is punctuated by a system for whole-body action planning, the somato-cognitive action network (SCAN). In M1, two parallel systems intertwine, forming an integrate-isolate pattern: effector-specific regions (foot, hand and mouth) for isolating fine motor control and the SCAN for integrating goals, physiology and body movement.

Functional Connectivity Mapping for rTMS Target Selection in Depression.

OBJECTIVE: Repetitive transcranial magnetic stimulation (rTMS) protocols increasingly use subgenual anterior cingulate cortex (sgACC) functional connectivity to individualize treatment targets. However, the efficacy of this approach is unclear, with conflicting findings and varying effect sizes across studies. Here, the authors investigated the effect of the stimulation site's functional connectivity with the sgACC (sgACC-StimFC) on treatment outcome to rTMS in 295 patients with major depression. METHODS: The reliability and accuracy of estimating sgACC functional connectivity were validated with data from individuals who underwent extensive functional MRI testing. Electric field modeling was used to analyze associations between sgACC-StimFC and clinical improvement using standardized assessments and to evaluate sources of heterogeneity. RESULTS: An imputation-based method provided reliable and accurate sgACC functional connectivity estimates. Treatment responses weakly but robustly correlated with sgACC-StimFC (r=-0.16), but only when the stimulated cortex was identified using electric field modeling. Surprisingly, this association was driven by patients with strong global signal fluctuations stemming from a specific periodic respiratory pattern (r=-0.49). CONCLUSIONS: Functional connectivity between the sgACC and the stimulated cortex was correlated with individual differences in treatment outcomes, but the association was weaker than those observed in previous studies and was accentuated in a subgroup of patients with distinct, respiration-related signal patterns in their scans. These findings indicate that in a large representative sample of patients with major depressive disorder, individual differences in sgACC-StimFC explained only ∼3% of the variance in outcomes, which may limit the utility of existing sgACC-based targeting protocols. However, these data also provide strong evidence for a true-albeit small-effect and highlight opportunities for incorporating additional functional connectivity measures to generate models of rTMS response with enhanced predictive power.

Precision mapping and transcranial magnetic stimulation of individual-specific functional brain networks in humans.

Spatial targeting in transcranial magnetic stimulation protocols does not typically account for the idiosyncratic functional organization of individual human brains. Here, we provide a protocol for implementing targeted functional network stimulation (TANS), which accounts for each individual's unique functional neuroanatomy and cortical folding patterns. Using an example dataset, we describe how to create a head model and estimate the best coil placement and stimulation intensity to minimize off-target effects. For complete details on the use and execution of this protocol, please refer to Lynch et al. (2022).1.

Automated optimization of TMS coil placement for personalized functional network engagement.

Transcranial magnetic stimulation (TMS) is used to treat multiple psychiatric and neurological conditions by manipulating activity in particular brain networks and circuits, but individual responses are highly variable. In clinical settings, TMS coil placement is typically based on either group average functional maps or scalp heuristics. Here, we found that this approach can inadvertently target different functional networks in depressed patients due to variability in their functional brain organization. More precise TMS targeting should be feasible by accounting for each patient's unique functional neuroanatomy. To this end, we developed a targeting approach, termed targeted functional network stimulation (TANS). The TANS approach improved stimulation specificity in silico in 8 highly sampled patients with depression and 6 healthy individuals and in vivo when targeting somatomotor functional networks representing the upper and lower limbs. Code for implementing TANS and an example dataset are provided as a resource.

An effort-based social feedback paradigm reveals aversion to popularity in socially anxious participants and increased motivation in adolescents.

We created a novel social feedback paradigm to study how motivation for potential social links is influenced in adolescents and adults. 88 participants (42F/46M) created online posts and then expended physical effort to show their posts to other users, who varied in number of followers and probability of positive feedback. We focused on two populations of particular interest from a social feedback perspective: adolescents relative to young adults (13-17 vs 18-24 years of age), and participants with social anxiety symptoms. Individuals with higher self-reported symptoms of social anxiety did not follow the typical pattern of increased effort to obtain social feedback from high status peers. Adolescents were more willing to exert physical effort on the task than young adults. Overall, participants were more likely to exert physical effort for high social status users and for users likely to yield positive feedback, and men were more likely to exert effort than women, findings that parallel prior results in effort-based tasks with financial rather than social rewards. Together the findings suggest social motivation is malleable, driven by factors of social status and the likelihood of a positive social outcome, and that age, sex, and social anxiety significantly impact patterns of socially motivated decision-making.

Rapid Precision Functional Mapping of Individuals Using Multi-Echo fMRI.

Resting-state functional magnetic resonance imaging (fMRI) is widely used in cognitive and clinical neuroscience, but long-duration scans are currently needed to reliably characterize individual differences in functional connectivity (FC) and brain network topology. In this report, we demonstrate that multi-echo fMRI can improve the reliability of FC-based measurements. In four densely sampled individual humans, just 10 min of multi-echo data yielded better test-retest reliability than 30 min of single-echo data in independent datasets. This effect is pronounced in clinically important brain regions, including the subgenual cingulate, basal ganglia, and cerebellum, and is linked to three biophysical signal mechanisms (thermal noise, regional variability in the rate of T2∗ decay, and S0-dependent artifacts) with spatially distinct influences. Together, these findings establish the potential utility of multi-echo fMRI for rapid precision mapping using experimentally and clinically tractable scan times and will facilitate longitudinal neuroimaging of clinical populations.

Prevalent and sex-biased breathing patterns modify functional connectivity MRI in young adults.

Resting state functional connectivity magnetic resonance imaging (fMRI) is a tool for investigating human brain organization. Here we identify, visually and algorithmically, two prevalent influences on fMRI signals during 440 h of resting state scans in 440 healthy young adults, both caused by deviations from normal breathing which we term deep breaths and bursts. The two respiratory patterns have distinct influences on fMRI signals and signal covariance, distinct timescales, distinct cardiovascular correlates, and distinct tendencies to manifest by sex. Deep breaths are not sex-biased. Bursts, which are serial taperings of respiratory depth typically spanning minutes at a time, are more common in males. Bursts share features of chemoreflex-driven clinical breathing patterns that also occur primarily in males, with notable neurological, psychiatric, medical, and lifespan associations. These results identify common breathing patterns in healthy young adults with distinct influences on functional connectivity and an ability to differentially influence resting state fMRI studies.

A Critical, Event-Related Appraisal of Denoising in Resting-State fMRI Studies.

This article advances two parallel lines of argument about resting-state functional magnetic resonance imaging (fMRI) signals, one empirical and one conceptual. The empirical line creates a four-part organization of the text: (1) head motion and respiration commonly cause distinct, major, unwanted influences (artifacts) in fMRI signals; (2) head motion and respiratory changes are, confoundingly, both related to psychological and clinical and biological variables of interest; (3) many fMRI denoising strategies fail to identify and remove one or the other kind of artifact; and (4) unremoved artifact, due to correlations of artifacts with variables of interest, renders studies susceptible to identifying variance of noninterest as variance of interest. Arising from these empirical observations is a conceptual argument: that an event-related approach to task-free scans, targeting common behaviors during scanning, enables fundamental distinctions among the kinds of signals present in the data, information which is vital to understanding the effects of denoising procedures. This event-related perspective permits statements like "Event X is associated with signals A, B, and C, each with particular spatial, temporal, and signal decay properties". Denoising approaches can then be tailored, via performance in known events, to permit or suppress certain kinds of signals based on their desirability.

Characteristics of respiratory measures in young adults scanned at rest, including systematic changes and "missed" deep breaths.

Breathing rate and depth influence the concentration of carbon dioxide in the blood, altering cerebral blood flow and thus functional magnetic resonance imaging (fMRI) signals. Such respiratory fluctuations can have substantial influence in studies of fMRI signal covariance in subjects at rest, the so-called "resting state functional connectivity" technique. If respiration is monitored during fMRI scanning, it is typically done using a belt about the subject's abdomen to record abdominal circumference. Several measures have been derived from these belt records, including the windowed envelope of the waveform (ENV), the windowed variance in the waveform (respiration variation, RV), and a measure of the amplitude of each breath divided by the cycle time of the breath (respiration volume per time, RVT). Any attempt to gauge respiratory contributions to fMRI signals requires a respiratory measure, but little is known about how these measures compare to each other, or how they perform beyond the small studies in which they were initially proposed. Here, we examine the properties of these measures in hundreds of healthy young adults scanned for an hour each at rest, a subset of the Human Connectome Project chosen for having high-quality physiological records. We find: 1) ENV, RV, and RVT are all correlated, and ENV and RV are more highly correlated to each other than to RVT; 2) respiratory events like deep breaths exhibit characteristic heart rate elevations, fMRI signal changes, head motions, and image quality abnormalities time-locked to large deflections in the belt traces; 3) all measures can "miss" deep breaths; 4) RVT "misses" deep breaths more than ENV or RV; 5) all respiratory measures change systematically over the course of a 14.4-min scan. We discuss the implications of these findings for the literature and ways to move forward in modeling respiratory influences on fMRI scans.

Distinctions among real and apparent respiratory motions in human fMRI data.

Head motion estimates in functional magnetic resonance imaging (fMRI) scans appear qualitatively different with sub-second image sampling rates compared to the multi-second sampling rates common in the past. Whereas formerly the head appeared still for much of a scan with brief excursions from baseline, the head now appears to be in constant motion, and motion estimates often seem to divulge little information about what is happening in a scan. This constant motion has been attributed to respiratory oscillations that do not alias at faster sampling rates, and investigators are divided on the extent to which such motion is "real" motion or only "apparent" pseudomotion. Some investigators have abandoned the use of motion estimates entirely due to these considerations. Here we investigate the properties of motion in several fMRI datasets sampled at rates between 720 and 1160 ms, and describe 5 distinct kinds of respiratory motion: 1) constant real respiratory motion in the form of head nodding most evident in vertical position and pitch, which can be very large; 2) constant pseudomotion at the same respiratory rate as real motion, occurring only in the phase encode direction; 3) punctate real motions occurring at times of very deep breaths; 4) a low-frequency pseudomotion in only the phase encode direction at and after very deep breaths; 5) slow modulation of vertical and anterior-posterior head position by the respiratory envelope. We reformulate motion estimates in light of these considerations and obtain good concordance between motion estimates, physiologic records, image quality measures, and events evident in the fMRI signals. We demonstrate how variables describing respiration or body habitus separately scale with distinct kinds of head motion. We also note heritable aspects of respiration and motion.

Customized head molds reduce motion during resting state fMRI scans.

Head motion causes artifacts in functional magnetic resonance imaging (fMRI) scans, a problem especially relevant for task-free resting state paradigms and for developmental, aging, and clinical populations. In a cohort spanning 7-28 years old (mean age 15) we produced customized head-anatomy-specific Styrofoam molds for each subject that inserted into an MRI head coil. We scanned these subjects under two conditions: using our standard procedure of packing the head coil with foam padding about the head to reduce head motion, and using the customized molds to reduce head motion. In 12 of 13 subjects, the molds reduced head motion throughout the scan and reduced the fraction of a scan with substantial motion (i.e., volumes with motion notably above baseline levels of motion). Motion was reduced in all 6 head position estimates, especially in rotational, left-right, and superior-inferior directions. Motion was reduced throughout the full age range studied, including children, adolescents, and young adults. In terms of the fMRI data itself, quality indices improved with the head mold on, scrubbing analyses detected less distance-dependent artifact in scans with the head mold on, and distant-dependent artifact was less evident in both the entire scan and also during only low-motion volumes. Subjects found the molds comfortable. Head molds are thus effective tools for reducing head motion, and motion artifacts, during fMRI scans.

Reward-related regions form a preferentially coupled system at rest.

Neuroimaging studies have implicated a set of striatal and orbitofrontal cortex (OFC) regions that are commonly activated during reward processing tasks. Resting-state functional connectivity (RSFC) studies have demonstrated that the human brain is organized into several functional systems that show strong temporal coherence in the absence of goal-directed tasks. Here we use seed-based and graph-theory RSFC approaches to characterize the systems-level organization of putative reward regions of at rest. Peaks of connectivity from seed-based RSFC patterns for the nucleus accumbens (NAcc) and orbitofrontal cortex (OFC) were used to identify candidate reward regions which were merged with a previously used set of regions (Power et al., 2011). Graph-theory was then used to determine system-level membership for all regions. Several regions previously implicated in reward-processing (NAcc, lateral and medial OFC, and ventromedial prefrontal cortex) comprised a distinct, preferentially coupled system. This RSFC system is stable across a range of connectivity thresholds and shares strong overlap with meta-analyses of task-based reward studies. This reward system shares between-system connectivity with systems implicated in cognitive control and self-regulation, including the fronto-parietal, cingulo-opercular, and default systems. Differences may exist in the pathways through which control systems interact with reward system components. Whereas NAcc is functionally connected to cingulo-opercular and default systems, OFC regions show stronger connectivity with the fronto-parietal system. We propose that future work may be able to interrogate group or individual differences in connectivity profiles using the regions delineated in this work to explore potential relationships to appetitive behaviors, self-regulation failure, and addiction.

Ridding fMRI data of motion-related influences: Removal of signals with distinct spatial and physical bases in multiecho data.

"Functional connectivity" techniques are commonplace tools for studying brain organization. A critical element of these analyses is to distinguish variance due to neurobiological signals from variance due to nonneurobiological signals. Multiecho fMRI techniques are a promising means for making such distinctions based on signal decay properties. Here, we report that multiecho fMRI techniques enable excellent removal of certain kinds of artifactual variance, namely, spatially focal artifacts due to motion. By removing these artifacts, multiecho techniques reveal frequent, large-amplitude blood oxygen level-dependent (BOLD) signal changes present across all gray matter that are also linked to motion. These whole-brain BOLD signals could reflect widespread neural processes or other processes, such as alterations in blood partial pressure of carbon dioxide (pCO2) due to ventilation changes. By acquiring multiecho data while monitoring breathing, we demonstrate that whole-brain BOLD signals in the resting state are often caused by changes in breathing that co-occur with head motion. These widespread respiratory fMRI signals cannot be isolated from neurobiological signals by multiecho techniques because they occur via the same BOLD mechanism. Respiratory signals must therefore be removed by some other technique to isolate neurobiological covariance in fMRI time series. Several methods for removing global artifacts are demonstrated and compared, and were found to yield fMRI time series essentially free of motion-related influences. These results identify two kinds of motion-associated fMRI variance, with different physical mechanisms and spatial profiles, each of which strongly and differentially influences functional connectivity patterns. Distance-dependent patterns in covariance are nearly entirely attributable to non-BOLD artifacts.

Neural correlates of preferred activities: development of an interest-specific go/nogo task.

The activities we choose to spend our leisure time with are intrinsically motivating and vary across individuals. Yet it is unknown how impulse control or neural activity changes when processing a preferred stimulus related to a hobby or interest. Developing a task that assesses the response to preferred interests is of importance as it would be relevant to a range of psychiatric disorders that have hyper- or hypo-arousal to such cues. During functional Magnetic Resonance Imaging (fMRI), 39 healthy adults completed a novel task to test approach behavior and cognitive control to cues that were personalized to the participants' interests compared to stimuli the participants identified as being of non-interest and colored shapes. fMRI results showed that cues of one's interest elicited activation in the anterior insula compared to colored shapes. Interests did not change inhibition compared to non-interests and colored shapes and all stimuli equally engaged a frontostriatal circuit. Together the results suggest that adults were sensitive to their interests but were effective at regulating their impulses towards these cues, a skill that is critical for navigating the temptations and distractions in our daily environment.