Search

Summary: During adolescence, the brain undergoes extensive changes in white matter structure that support cognition. Data-driven approaches applied to cortical surface properties have led the field to understand brain development as a spatially and temporally coordinated mechanism that follows hierarchically organized gradients of change. Although white matter development also appears asynchronous, previous studies have relied largely on anatomical tract-based atlases, precluding a direct assessment of how white matter structure is spatially and temporally coordinated. Harnessing advances in diffusion modeling and machine learning, we identified 14 data-driven patterns of covarying white matter structure in a large sample of youth. Fiber covariance networks aligned with known major tracts, while also capturing distinct patterns of spatial covariance across distributed white matter locations. Most networks showed age-related increases in fiber network properties, which were also related to developmental changes in executive function. This study delineates data-driven patterns of white matter development that support cognition.

Summary: Human brain size changes dynamically through early development, peaks in adolescence, and varies up to 2-fold among adults. However, the molecular genetic underpinnings of interindividual variation in brain size remain unknown. Here, we leveraged postmortem brain RNA sequencing and measurements of brain weight (BW) in 2,531 individuals across three independent datasets to identify 928 genome-wide significant associations with BW. Genes associated with higher or lower BW showed distinct neurodevelopmental trajectories and spatial patterns that mapped onto functional and cellular axes of brain organization. Expression of BW genes was predictive of interspecies differences in brain size, and bioinformatic annotation revealed enrichment for neurogenesis and cell-cell communication. Genome-wide, transcriptome-wide, and phenome-wide association analyses linked BW gene sets to neuroimaging measurements of brain size and brain-related clinical traits. Cumulatively, these results represent a major step toward delineating the molecular pathways underlying human brain size variation in health and disease.

Background: Impairment in intrinsic motivation (IM), the drive to satisfy internal desires like mastery, may play a key role in disability in psychosis. However, we have limited knowledge regarding relative impairments in IM compared to extrinsic motivation (EM) or general motivation (GM), in part due to limitations in existing measures. Methods: Here we address this gap using a novel Trait Intrinsic and Extrinsic Motivation self-report scale in a sample of n = 243 participants including those with schizophrenia, psychosis-risk, and healthy controls. Each of the 7 IM and 6 EM items used a 7-point Likert scale assessing endorsement of dispositional statements. Bifactor analyses of these items yielded distinct IM, EM, and GM factor scores. Convergent and discriminant validity were examined in relation to General Causality Orientation Scale (GCOS-CP) and Quality of Life 3-item IM measure (QLS-IM). Utility was assessed in relation to psychosis-spectrum (PS) status and CAINS clinical amotivation. Results: IM and EM showed acceptable inter-item consistency (IM: α = 0.88; EM: α = 0.66); the bifactor model exhibited fit that varied from good to borderline to inadequate depending on the specific fit metric (SRMR = 0.038, CFI = 0.94, RMSEA = 0.106 ± 0.014). IM scores correlated with established IM measures: GCOS-CP Autonomy (rho = 0.38, p 0.05). Two-year stability in an available longitudinal subset (n = 35) was strong (IM: rho = 0.64, p

Delay discounting is a measure of impulsive choice relevant in adolescence as it predicts many real-life outcomes, including obesity and academic achievement. However, resting-state functional networks underlying individual differences in delay discounting during youth remain incompletely described. Here we investigate the association between multivariate patterns of functional connectivity and individual differences in impulsive choice in a large sample of children, adolescents, and adults. A total of 293 participants (9–23 years) completed a delay discounting task and underwent 3T resting-state fMRI. A connectome-wide analysis using multivariate distance-based matrix regression was used to examine whole-brain relationships between delay discounting and functional connectivity. These analyses revealed that individual differences in delay discounting were associated with patterns of connectivity emanating from the left dorsal prefrontal cortex, a default mode network hub. Greater delay discounting was associated with greater functional connectivity between the dorsal prefrontal cortex and other default mode network regions, but reduced connectivity with regions in the dorsal and ventral attention networks. These results suggest delay discounting in children, adolescents, and adults is associated with individual differences in relationships both within the default mode network and between the default mode and networks involved in attentional and cognitive control.

ABSTRACT: Diffusion MRI is the dominant non-invasive imaging method used to characterize white matter organization in health and disease. Increasingly, fiber-specific properties within a voxel are analyzed using fixels. While tools for conducting statistical analyses of fixel-wise data exist, currently available tools support only a limited number of statistical models. Here we introduce ModelArray, an R package for mass-univariate statistical analysis of fixel-wise data. At present, ModelArray supports linear models as well as generalized additive models (GAMs), which are particularly useful for studying nonlinear effects in lifespan data. In addition, ModelArray also aims for scalable analysis. With only several lines of code, even large fixel-wise datasets can be analyzed using a standard personal computer. Detailed memory profiling revealed that ModelArray required only limited memory even for large datasets. As an example, we applied ModelArray to fixel-wise data derived from diffusion images acquired as part of the Philadelphia Neurodevelopmental Cohort (n = 938). ModelArray revealed anticipated nonlinear developmental effects in white matter. Moving forward, ModelArray is supported by an open-source software development model that can incorporate additional statistical models and other imaging data types. Taken together, ModelArray provides a flexible and efficient platform for statistical analysis of fixel-wise data.

Background: Impulse control behaviors (ICB) are frequently observed in patients with Parkinson’s disease (PD) and are characterized by compulsive and repetitive behavior resulting from the inability to resist internal drives. Objectives: In this study, we aimed to provide a better understanding of structural and functional brain alterations and clinical parameters related to ICB in PD patients. Methods: We utilized a dataset from the Parkinson’s Progression Markers Initiative including 36 patients with ICB (PDICB+) compared to 76 without ICB (PDICB-) and 61 healthy controls (HC). Using multimodal MRI data we assessed gray matter brain volume, white matter integrity, and graph topological properties at rest. Results: Compared with HC, PDICB+ showed reduced gray matter volume in the bilateral superior and middle temporal gyrus and in the right middle occipital gyrus. Compared with PDICB-, PDICB+ showed volume reduction in the left anterior insula. Depression and anxiety were more prevalent in PDICB+ than in PDICB- and HC. In PDICB+, lower gray matter volume in the precentral gyrus and medial frontal cortex, and higher axial diffusivity in the superior corona radiata were related to higher depression score. Both PD groups showed disrupted functional topological network pattern within the cingulate cortex compared with HC. PDICB+ vs PDICB- displayed reduced topological network pattern in the anterior cingulate cortex, insula, and nucleus accumbens. Conclusions: Our results suggest that structural alterations in the insula and abnormal topological connectivity pattern in the salience network and the nucleus accumbens may lead to impaired decision making and hypersensitivity towards reward in PDICB+. Moreover, PDICB+ are more prone to suffer from depression and anxiety.

With the increasing availability of neuroimaging data from multiple modalities—each providing a different lens through which to study brain structure or function—new techniques for comparing, integrating, and interpreting information within and across modalities have emerged. Recent developments include hypothesis tests of associations between neuroimaging modalities, which can be used to determine the statistical significance of intermodal associations either throughout the entire brain or within anatomical subregions or functional networks. While these methods provide a crucial foundation for inference on intermodal relationships, they cannot be used to answer questions about where in the brain these associations are most pronounced. In this paper, we introduce a new method, called CLEAN-R, that can be used both to test intermodal correspondence throughout the brain and also to localize this correspondence. Our method involves first adjusting for the underlying spatial autocorrelation structure within each modality before aggregating information within small clusters to construct a map of enhanced test statistics. Using structural and functional magnetic resonance imaging data from a subsample of children and adolescents from the Philadelphia Neurodevelopmental Cohort, we conduct simulations and data analyses where we illustrate the high statistical power and nominal type I error levels of our method. By constructing an interpretable map of group-level correspondence using spatially-enhanced test statistics, our method offers insights beyond those provided by earlier methods.

The Brain Imaging Data Structure (BIDS) is a specification accompanied by a software ecosystem that was designed to create reproducible and automated workflows for processing neuroimaging data. BIDS Apps flexibly build workflows based on the metadata detected in a dataset. However, even BIDS valid metadata can include incorrect values or omissions that result in inconsistent processing across sessions. Additionally, in large-scale, heterogeneous neuroimaging datasets, hidden variability in metadata is difficult to detect and classify. To address these challenges, we created a Python-based software package titled “Curation of BIDS” (CuBIDS), which provides an intuitive workflow that helps users validate and manage the curation of their neuroimaging datasets. CuBIDS includes a robust implementation of BIDS validation that scales to large samples and incorporates DataLad––a version control software package for data––as an optional dependency to ensure reproducibility and provenance tracking throughout the entire curation process. CuBIDS provides tools to help users perform quality control on their images’ metadata and identify unique combinations of imaging parameters. Users can then execute BIDS Apps on a subset of participants that represent the full range of acquisition parameters that are present, accelerating pipeline testing on large datasets.

Summary: Polygenic scores (PGS) are commonly evaluated in terms of their predictive accuracy at the population level by the proportion of phenotypic variance they explain. To be useful for precision medicine applications, they also need to be evaluated at the individual level when phenotypes are not necessarily already known. We investigated the stability of PGS in European American (EUR) and African American (AFR)-ancestry individuals from the Philadelphia Neurodevelopmental Cohort and the Adolescent Brain Cognitive Development study using different discovery genome-wide association study (GWAS) results for post-traumatic stress disorder (PTSD), type 2 diabetes (T2D), and height. We found that pairs of EUR-ancestry GWAS for the same trait had genetic correlations >0.92. However, PGS calculated from pairs of same-ancestry and different-ancestry GWAS had correlations that ranged from

Environment (E) is pivotal in explaining variability in brain and behavior development, including suicidal ideation (SI) and behavior. It is therefore critical to systematically study relationships among environmental exposures (i.e., exposome) and suicidal phenotypes. Here, we evaluated the role of individual-level adversity and neighborhood environment and their interaction (E x E) in association with youth SI. Sample included youth (N = 7,054, ages 11–21) from the Philadelphia Neurodevelopmental Cohort, which investigated clinical phenotypes in a diverse US community population. We examined cross-sectional associations of environmental exposures with lifetime history of SI (n = 671), focusing on interactions between individual-level exposures to assaultive trauma (n = 917) and neighborhood-level socioeconomic status (SES) quantified using geocoded Census data. Models included potential confounds and overall psychopathology. Results showed that assaultive trauma was strongly associated with SI (OR = 3.3, 95%CI 2.7–4, p

Objective: The limbic system is involved in memory and in processing of emotional stimuli. We measured volume of the hippocampus, amygdala, and thalamus, and assessed their relative contribution to episodic memory and emotion identification in POMS. Method: Sixty-five POMS participants (Mage = 18.3 ± 3.9 years; 48 female (73.8%)), average disease duration = 3.8 ± 3.8 years) and 76 age- and sex-matched controls (Mage = 18.1 ± 4.6 years; 49 female (64.5%)) completed the Penn Computerized Neurocognitive Battery (PCNB); 59 of 65 POMS participants and 69 out of 76 controls underwent 3 T MRI scanning. We derived age-adjusted Z-scores on accuracy and response time (RT) measures of episodic memory and emotion identification of the PCNB. Magnetic resonance imaging (MRI) volumetrics were normalized using the scaling factor computed by SIENAx. On PCNB tests that differed between groups, we used multiple linear regression to assess relationships between regional brain volumes and either episodic memory or emotion identification outcomes controlling for age, sex, accuracy/RT, and parental education. Results: POMS participants were slower and less accurate than controls on the episodic memory domain but did not differ from controls on emotion outcomes. At the subtest level, POMS participants showed reduced accuracy on Word Memory (p = .002) and slower performance on Face Memory (p = .04) subtests. POMS participants had smaller total and regional brain volumes of the hippocampus, amygdala, and thalamus (p values ≤ 0.01). Collapsing across groups, both hippocampal and thalamic volume were significant predictors of Word Memory accuracy; hippocampal volume (B = 0.24, SE = 0.10, p = .02) was more strongly associated with Word Memory performance than thalamic volume (B = 0.16, SE = 0.05, p = .003), though the estimate with was less precise. Conclusions: POMS participants showed reduced episodic memory performance compared to controls. Aspects of episodic memory performance were associated with hippocampal and thalamic volume. Emotion identification was intact, despite volume loss in the amygdala.

Alterations in resting-state functional connectivity (rsFC) have been demonstrated in Posttraumatic Stress Disorder (PTSD). However, such reports have primarily focused on adult participants, whereas findings in adolescents with PTSD are mixed and not entirely consistent with the adult literature. Here, we examined rsFC in a non-treatment seeking adolescent sample with posttraumatic stress symptoms (PTSS; n = 59) relative to asymptomatic controls (n = 226). We also examined differences between trauma-exposed and non-exposed control subgroups (TEC n = 73 and Non-TEC n = 153) to examine alterations associated with more general trauma exposure. Finally, we compared the PTSS and TEC groups, to confirm that the reported alterations in PTSS were not driven by trauma exposure. Using a seed-based approach, we examined connectivity of default-mode (DMN) and salience (SN) networks, where alterations have been previously reported. Results suggest that PTSS are associated with less within-DMN connectivity and greater SN-DMN connectivity, as well as altered connectivity with attention regions. Trauma exposure is associated with greater within-SN connectivity. Additionally, we report findings from exploratory connectome-based analysis, which demonstrate a number of topological alterations within DMN in the PTSS group. Overall, our findings replicate prior reports of altered rsFC in PTSD and extend them to non-treatment seeking, trauma-exposed adolescents, who did or did not report PTSS. They specifically highlight SN-DMN desegregation, lower within-DMN and greater within-SN connectivity, as well as altered connectivity with attention regions, in trauma-exposed adolescents. Future research is required to confirm that adolescents with diagnosed PTSD have similar/exacerbated connectivity patterns.

To date, most transcriptome studies of schizophrenia focus on the analysis of protein-coding genes. Long noncoding RNAs (lncRNAs) are emerging as key tissue-specific regulators of cellular and disease processes. The amygdala brain region has been implicated in the pathophysiology of schizophrenia. We performed unbiased whole transcriptome profiling of amygdala tissues from 22 schizophrenia patients and 24 non-psychiatric controls using RNA-seq. We reconstructed amygdala transcriptome and employed systems biology approaches to annotating the functional roles of lncRNAs. As a result, we identified 839 novel lncRNAs in amygdala. We found in amygdala lncRNAs are more subtype-specific than protein-coding genes. We identified functional modules associated with “synaptic transmission”, “ribosome”, and “immune responses” which were related to schizophrenia pathophysiology that involved lncRNAs. Integrative functional analyses associating individual lncRNAs with specific pathways and functions further show that amygdala lncRNAs are connected with all of these pathways. Our study presents the first systematic landscape of lncRNAs in amygdala tissue from schizophrenia cases. Keywords: Schizophrenia, Next-generation sequencing, Transcriptome, Systems biology

With the advent of Big Data Imaging Analytics applied to neuroimaging, datasets from multiple sites need to be pooled into larger samples. However, heterogeneity across different scanners, protocols and populations, renders the task of finding underlying disease signatures challenging. The current work investigates the value of multi-task learning in finding disease signatures that generalize across studies and populations. Herein, we present a multi-task learning type of formulation, in which different tasks are from different studies and populations being pooled together. We test this approach in an MRI study of the neuroanatomy of schizophrenia (SCZ) by pooling data from 3 different sites and populations: Philadelphia, Sao Paulo and Tianjin (50 controls and 50 patients from each site), which posed integration challenges due to variability in disease chronicity, treatment exposure, and data collection. Some existing methods are also tested for comparison purposes. Experiments show that classification accuracy of multi-site data outperformed that of single-site data and pooled data using multi-task feature learning, and also outperformed other comparison methods. Several anatomical regions were identified to be common discriminant features across sites. These included prefrontal, superior temporal, insular, anterior cingulate cortex, temporo-limbic and striatal regions consistently implicated in the pathophysiology of schizophrenia, as well as the cerebellum, precuneus, and fusiform, middle temporal, inferior parietal, postcentral, angular, lingual and middle occipital gyri. These results indicate that the proposed multi-task learning method is robust in finding consistent and reliable structural brain abnormalities associated with SCZ across different sites, in the presence of multiple sources of heterogeneity. Keywords: Multi-task learning, Multi-site classification, Imaging heterogeneity, Sparsity, Schizophrenia, MRI

The 22q11.2 deletion syndrome (22q11DS) is an uncommon genetic disorder with an increased risk of psychosis. Although the neural substrates of psychosis and schizophrenia are not well understood, aberrations in cortical networks represent intriguing potential mechanisms. Investigations of anatomic networks within 22q11DS are sparse. We investigated group differences in anatomic network structure in 48 individuals with 22q11DS and 370 typically developing controls by analyzing covariance patterns in cortical thickness among 68 regions of interest using graph theoretical models. Subjects with 22q11DS had less robust geographic organization relative to the control group, particularly in the occipital and parietal lobes. Multiple global graph theoretical statistics were decreased in 22q11DS. These results are consistent with prior studies demonstrating decreased connectivity in 22q11DS using other neuroimaging methodologies.