Functional MRI (fMRI) has been widely used to examine changes in neuronal activity during cognitive tasks. Commonly used measures of gray matter macrostructure (e.g., cortical thickness, surface ...area, volume) do not consistently appear to serve as structural correlates of brain function. In contrast, gray matter microstructure, measured using neurite orientation dispersion and density imaging (NODDI), enables the estimation of indices of neurite density (neurite density index; NDI) and organization (orientation dispersion index; ODI) in gray matter. Our study explored the relationship among neurite architecture, BOLD (blood-oxygen-level-dependent) fMRI, and cognition, using a large sample (n = 750) of young adults of the human connectome project (HCP) and two tasks that index more cortical (working memory) and more subcortical (emotion processing) targeting of brain functions. Using NODDI, fMRI, structural MRI and task performance data, hierarchical regression analyses revealed that higher working memory- and emotion processing-evoked BOLD activity was related to lower ODI in the right DLPFC, and lower ODI and NDI values in the right and left amygdala, respectively. Common measures of brain macrostructure (i.e., DLPFC thickness/surface area and amygdala volume) did not explain any additional variance (beyond neurite architecture) in BOLD activity. A moderating effect of neurite architecture on the relationship between emotion processing task-evoked BOLD response and performance was observed. Our findings provide evidence that neuro-/social-affective cognition-related BOLD activity is partially driven by the local neurite organization and density with direct impact on emotion processing. In vivo gray matter microstructure represents a new target of investigation providing strong potential for clinical translation.
Autism spectrum disorder (ASD), obsessive-compulsive disorder (OCD) and attention-deficit/hyperactivity disorder (ADHD) are clinically and biologically heterogeneous neurodevelopmental disorders ...(NDDs). The objective of the present study was to integrate brain imaging and behavioral measures to identify new brain-behavior subgroups cutting across these disorders. A subset of the data from the Province of Ontario Neurodevelopmental Disorder (POND) Network was used including participants with different NDDs (aged 6-16 years) that underwent cross-sectional T1-weighted and diffusion-weighted magnetic resonance imaging (MRI) scanning on the same 3T scanner, and behavioral/cognitive assessments. Similarity Network Fusion was applied to integrate cortical thickness, subcortical volume, white matter fractional anisotropy (FA), and behavioral measures in 176 children with ASD, ADHD or OCD with complete data that passed quality control. Normalized mutual information was used to determine top contributing model features. Bootstrapping, out-of-model outcome measures and supervised machine learning were each used to examine stability and evaluate the new groups. Cortical thickness in socio-emotional and attention/executive networks and inattention symptoms comprised the top ten features driving participant similarity and differences between four transdiagnostic groups. Subcortical volumes (pallidum, nucleus accumbens, thalamus) were also different among groups, although white matter FA showed limited differences. Features driving participant similarity remained stable across resampling, and the new groups showed significantly different scores on everyday adaptive functioning. Our findings open the possibility of studying new data-driven groups that represent children with NDDs more similar to each other than others within their own diagnostic group. Future work is needed to build on this early attempt through replication of the current findings in independent samples and testing longitudinally for prognostic value.
•Neuroimaging research usually examines group averages, failing to capture critical variability.•Made use of task fMRI data (n = 822) from the Human Connectome Project to explore the range of ...variability in brain.•Extensive variability which was not well captured in the group averages, along a ‘positive’ to ‘negative’ pattern.•Participants did not segregate into distinct sub-groups, appearing to fall along multidimensional spectrums.•Variability related to out-of-scanner cognitive abilities.
Human neuroimaging during cognitive tasks has provided unique and important insights into the neurobiology of cognition. However, the vast majority of research relies on group aggregate or average statistical maps of activity, which do not fully capture the rich intersubject variability in brain function. In order to fully understand the neurobiology of cognitive processes, it is necessary to explore the range of variability in activation patterns across individuals. To better characterize individual variability, hierarchical clustering was performed separately on six fMRI tasks in 822 participants from the Human Connectome Project. Across all tasks, clusters ranged from a predominantly ‘deactivating’ pattern towards a more ‘activating’ pattern of brain activity, with significant differences in out-of-scanner cognitive test scores between clusters. Cluster stability was assessed via a resampling approach; a cluster probability matrix was generated, as the probability of any pair of participants clustering together when both were present in a random subsample. Rather than forming distinct clusters, participants fell along a spectrum or into pseudo-clusters without clear boundaries. A principal components analysis of the cluster probability matrix revealed three components explaining over 90% of the variance in clustering. Plotting participants in this lower-dimensional ‘similarity space’ revealed manifolds of variations along an S ‘snake’ shaped spectrum or a folded circle or ‘tortilla’ shape. The ‘snake’ shape was present in tasks where individual variability related to activity along covarying networks, while the ‘tortilla’ shape represented multiple networks which varied independently.
Myelinated axons form long-range connections that enable rapid communication between distant brain regions, but how genetics governs the strength and organization of these connections remains ...unclear. We perform genome-wide association studies of 206 structural connectivity measures derived from diffusion magnetic resonance imaging tractography of 26,333 UK Biobank participants, each representing the density of myelinated connections within or between a pair of cortical networks, subcortical structures or cortical hemispheres. We identify 30 independent genome-wide significant variants after Bonferroni correction for the number of measures studied (126 variants at nominal genome-wide significance) implicating genes involved in myelination (SEMA3A), neurite elongation and guidance (NUAK1, STRN, DPYSL2, EPHA3, SEMA3A, HGF, SHTN1), neural cell proliferation and differentiation (GMNC, CELF4, HGF), neuronal migration (CCDC88C), cytoskeletal organization (CTTNBP2, MAPT, DAAM1, MYO16, PLEC), and brain metal transport (SLC39A8). These variants have four broad patterns of spatial association with structural connectivity: some have disproportionately strong associations with corticothalamic connectivity, interhemispheric connectivity, or both, while others are more spatially diffuse. Structural connectivity measures are highly polygenic, with a median of 9.1 percent of common variants estimated to have non-zero effects on each measure, and exhibited signatures of negative selection. Structural connectivity measures have significant genetic correlations with a variety of neuropsychiatric and cognitive traits, indicating that connectivity-altering variants tend to influence brain health and cognitive function. Heritability is enriched in regions with increased chromatin accessibility in adult oligodendrocytes (as well as microglia, inhibitory neurons and astrocytes) and multiple fetal cell types, suggesting that genetic control of structural connectivity is partially mediated by effects on myelination and early brain development. Our results indicate pervasive, pleiotropic, and spatially structured genetic control of white-matter structural connectivity via diverse neurodevelopmental pathways, and support the relevance of this genetic control to healthy brain function.
Repetitive transcranial magnetic stimulation (rTMS) to the dorsolateral prefrontal cortex (DLPFC) is effective in alleviating treatment-resistant depression (TRD). It has been proposed that regions ...within the left DLPFC that are anti-correlated with the right subgenual anterior cingulate cortex (sgACC) may represent optimal individualized target sites for high-frequency left rTMS (HFL).
This study aimed to explore the effects of low-frequency right rTMS (LFR) on left sgACC connectivity during concurrent TMS-fMRI.
34 TRD patients underwent an imaging session that included both a resting-state fMRI run (rs-fMRI
) and a run during which LFR was applied to the right DLPFC (TMS-fMRI). Participants subsequently completed four weeks of LFR treatment. The left sgACC functional connectivity was compared between the rs-fMRI
run and TMS-fMRI run. Personalized e-fields and a region-of-interest approach were used to calculate overlap of left sgACC functional connectivity at the TMS target and to assess for a relationship with treatment effects.
TMS-fMRI increased left sgACC functional connectivity to parietal regions within the ventral attention network; differences were not significantly associated with clinical improvements. Personalized e-fields were not significant in predicting treatment outcomes (p = 0.18).
This was the first study to examine left sgACC anti-correlation with the right DLPFC during an LFR rTMS protocol. In contrast to studies that targeted the left DLPFC, we did not find that higher anti-correlation was associated with clinical outcomes. Our results suggest that the antidepressant mechanism of action of LFR to the right DLPFC may be different than for HFL.
•Active single-pulse TMS generates similar waveform as auditory response but with larger amplitude.•Only active paired-pulse protocol induced significant cortical inhibition in dorsolateral ...prefrontal cortex.•Cortical inhibition is observed over all frequency bands in active long interval cortical inhibition.
We measured the neurophysiological responses of both active and sham transcranial magnetic stimulation (TMS) for both single pulse (SP) and paired pulse (PP; long interval cortical inhibition (LICI)) paradigms using TMS-EEG (electroencephalography).
Nineteen healthy subjects received active and sham (coil 90° tilted and touching the scalp) SP and PP TMS over the left dorsolateral prefrontal cortex (DLPFC). We measured excitability through SP TMS and inhibition (i.e., cortical inhibition (CI)) through PP TMS.
Cortical excitability indexed by area under the curve (AUC(25-275ms)) was significantly higher in the active compared to sham stimulation (F(1,18) = 43.737, p < 0.001, η2 = 0.708). Moreover, the amplitude of N100-P200 complex was significantly larger (F(1,18) = 9.118, p < 0.01, η2 = 0.336) with active stimulation (10.38 ± 9.576 µV) compared to sham (4.295 ± 2.323 µV). Significant interaction effects were also observed between active and sham stimulation for both the SP and PP (i.e., LICI) cortical responses. Finally, only active stimulation (CI = 0.64 ± 0.23, p < 0.001) resulted in significant cortical inhibition.
The significant differences between active and sham stimulation in both excitatory and inhibitory neurophysiological responses showed that active stimulation elicits responses from the cortex that are different from the non-specific effects of sham stimulation.
Our study reaffirms that TMS-EEG represents an effective tool to evaluate cortical neurophysiology with high fidelity.
Combined transcranial magnetic stimulation and electroencephalography (TMS-EEG) is an effective way to evaluate neurophysiological processes at the level of the cortex. To further characterize the ...TMS-evoked potential (TEP) generated with TMS-EEG, beyond the motor cortex, we aimed to distinguish between cortical reactivity to TMS versus non-specific somatosensory and auditory co-activations using both single-pulse and paired-pulse protocols at suprathreshold stimulation intensities over the left dorsolateral prefrontal cortex (DLPFC). Fifteen right-handed healthy participants received six blocks of stimulation including single and paired TMS delivered as active-masked (i.e., TMS-EEG with auditory masking and foam spacing), active-unmasked (TMS-EEG without auditory masking and foam spacing) and sham (sham TMS coil). We evaluated cortical excitability following single-pulse TMS, and cortical inhibition following a paired-pulse paradigm (long-interval cortical inhibition (LICI)). Repeated measure ANOVAs revealed significant differences in mean cortical evoked activity (CEA) of active-masked, active-unmasked, and sham conditions for both the single-pulse (F(1.76, 24.63) = 21.88, p < 0.001, η
= 0.61) and LICI (F(1.68, 23.49) = 10.09, p < 0.001, η
= 0.42) protocols. Furthermore, global mean field amplitude (GMFA) differed significantly across the three conditions for both single-pulse (F(1.85, 25.89) = 24.68, p < 0.001, η
= 0.64) and LICI (F(1.8, 25.16) = 14.29, p < 0.001, η
= 0.5). Finally, only active LICI protocols but not sham stimulation (active-masked (0.78 ± 0.16, P < 0.0001), active-unmasked (0.83 ± 0.25, P < 0.01)) resulted in significant signal inhibition. While previous findings of a significant somatosensory and auditory contribution to the evoked EEG signal are replicated by our study, an artifact attenuated cortical reactivity can reliably be measured in the TMS-EEG signal with suprathreshold stimulation of DLPFC. Artifact attenuation can be accomplished using standard procedures, and even when masked, the level of cortical reactivity is still far above what is produced by sham stimulation. Our study illustrates that TMS-EEG of DLPFC remains a valid investigational tool.
This study examined BOLD changes prior to interictal discharges in the EEG of patients with epilepsy. From a database of 143 EEG–fMRI studies, we selected the 16 data sets that showed both strong ...fMRI activation in the original analysis and only a single spike type in the EEG. Scans were then analyzed using seven model HRFs, peaking 3 or 1 s before the event, or 1, 3, 5, 7, or 9 s after it. An HRF was calculated using a deconvolution method for all activations seen in each analysis. The results showed that seven data sets had HRFs that peaked 1 s after the event or earlier, indicating a BOLD change starting prior to the spike seen on the scalp EEG. This is surprising since the BOLD change is expected to result from the spike. For most of the data sets with early peaking HRFs, the maximum activation in all of the statistical maps was when the model HRF peaked 1 s after the event, suggesting that the early activation was at least as important as any later activation. We suggest that this early activity is the result of neuronal changes occurring several seconds prior to a surface EEG event, but that these changes are not visible on the scalp. This is the first report of a BOLD response occurring several seconds prior to an interictal event seen on the scalp and could have important implications for our understanding of the generation of epileptic discharges.
Convergent data from imaging and postmortem brain transcriptome studies implicate corticolimbic circuit (CLC) dysregulation in the pathophysiology of depression. To more directly bridge these lines ...of work, we generated a novel transcriptome-based polygenic risk score (T-PRS), capturing subtle shifts toward depression-like gene expression patterns in key CLC regions, and mapped this T-PRS onto brain function and related depressive symptoms in a nonclinical sample of 478 young adults (225 men; age 19.79 +/- 1.24) from the Duke Neurogenetics Study. First, T-PRS was generated based on common functional SNPs shifting CLC gene expression toward a depression-like state. Next, we used multivariate partial least squares regression to map T-PRS onto whole-brain activity patterns during perceptual processing of social stimuli (i.e., human faces). For validation, we conducted a comparative analysis with a PRS summarizing depression risk variants identified by the Psychiatric Genomics Consortium (PGC-PRS). Sex was modeled as moderating factor. We showed that T-PRS was associated with widespread reductions in neural response to neutral faces in women and to emotional faces and shapes in men (multivariate p < 0.01). This female-specific reductions in neural response to neutral faces was also associated with PGC-PRS (multivariate p < 0.03). Reduced reactivity to neutral faces was further associated with increased self-reported anhedonia. We conclude that women with functional alleles mimicking the postmortem transcriptomic CLC signature of depression have blunted neural activity to social stimuli, which may be expressed as higher anhedonia.
•Overall, there are differences in the participants excluded from four different quality control approaches across three pediatric datasets.•In clinically enriched samples, the greatest ...correspondence of excluded participants was between automated and visual quality control procedures.•Implementing quality control led to the exclusion of younger participants and those with greater clinical impairments.•Specific QC approach implemented did not lead to measurable differences in clinical or brain metric characteristics.
Poor quality T1-weighted brain scans systematically affect the calculation of brain measures. Removing the influence of such scans requires identifying and excluding scans with noise and artefacts through a quality control (QC) procedure. While QC is critical for brain imaging analyses, it is not yet clear whether different QC approaches lead to the exclusion of the same participants. Further, the removal of poor-quality scans may unintentionally introduce a sampling bias by excluding the subset of participants who are younger and/or feature greater clinical impairment. This study had two aims: (1) examine whether different QC approaches applied to T1-weighted scans would exclude the same participants, and (2) examine how exclusion of poor-quality scans impacts specific demographic, clinical and brain measure characteristics between excluded and included participants in three large pediatric neuroimaging samples.
We used T1-weighted, resting-state fMRI, demographic and clinical data from the Province of Ontario Neurodevelopmental Disorders network (Aim 1: n = 553, Aim 2: n = 465), the Healthy Brain Network (Aim 1: n = 1051, Aim 2: n = 558), and the Philadelphia Neurodevelopmental Cohort (Aim 1: n = 1087; Aim 2: n = 619). Four different QC approaches were applied to T1-weighted MRI (visual QC, metric QC, automated QC, fMRI-derived QC). We used tetrachoric correlation and inter-rater reliability analyses to examine whether different QC approaches excluded the same participants. We examined differences in age, mental health symptoms, everyday/adaptive functioning, IQ and structural MRI-derived brain indices between participants that were included versus excluded following each QC approach.
Dataset-specific findings revealed mixed results with respect to overlap of QC exclusion. However, in POND and HBN, we found a moderate level of overlap between visual and automated QC approaches (rtet=0.52–0.59). Implementation of QC excluded younger participants, and tended to exclude those with lower IQ, and lower everyday/adaptive functioning scores across several approaches in a dataset-specific manner. Across nearly all datasets and QC approaches examined, excluded participants had lower estimates of cortical thickness and subcortical volume, but this effect did not differ by QC approach.
The results of this study provide insight into the influence of QC decisions on structural pediatric imaging analyses. While different QC approaches exclude different subsets of participants, the variation of influence of different QC approaches on clinical and brain metrics is minimal in large datasets. Overall, implementation of QC tends to exclude participants who are younger, and those who have more cognitive and functional impairment. Given that automated QC is standardized and can reduce between-study differences, the results of this study support the potential to use automated QC for large pediatric neuroimaging datasets.