Changes in brain connectivity and neurovascular dynamics during dexmedetomidine-induced loss of consciousness.

Understanding the neurophysiological changes that occur during loss and recovery of consciousness is a fundamental aim in neuroscience and has marked clinical relevance. Here, we utilize multimodal magnetic resonance neuroimaging to investigate changes in regional network connectivity and neurovascular dynamics as the brain transitions from wakefulness to dexmedetomidine-induced unconsciousness, and finally into early-stage recovery of consciousness. We observed widespread decreases in functional connectivity strength across the whole brain, and targeted increases in structure-function coupling (SFC) across select networks- especially the cerebellum-as individuals transitioned from wakefulness to hypnosis. We also observed robust decreases in cerebral blood flow (CBF) across the whole brain-especially within the brainstem, thalamus, and cerebellum. Moreover, hypnosis was characterized by significant increases in the amplitude of low-frequency fluctuations (ALFF) of the resting-state blood oxygen level-dependent signal, localized within visual and somatomotor regions. Critically, when transitioning from hypnosis to the early stages of recovery, functional connectivity strength and SFC-but not CBF-started reverting towards their awake levels, even before behavioral arousal. By further testing for a relationship between connectivity and neurovascular alterations, we observed that during wakefulness, brain regions with higher ALFF displayed lower functional connectivity with the rest of the brain. During hypnosis, brain regions with higher ALFF displayed weaker coupling between structural and functional connectivity. Correspondingly, brain regions with stronger functional connectivity strength during wakefulness showed greater reductions in CBF with the onset of hypnosis. Earlier recovery of consciousness was associated with higher baseline (awake) levels of functional connectivity strength, CBF, and ALFF, as well as female sex. Across our findings, we also highlight the role of the cerebellum as a recurrent marker of connectivity and neurovascular changes between states of consciousness. Collectively, these results demonstrate that induction of, and emergence from dexmedetomidine-induced unconsciousness are characterized by widespread changes in connectivity and neurovascular dynamics. SIGNIFICANCE STATEMENT: Elucidating the neurophysiological changes underlying loss and recovery of consciousness is a fundamental question in neuroscience. Here, we analyze magnetic resonance imaging data collected across multiple time-points to characterize how the human brain's connectivity and neurovascular dynamics change as it transitions from wakefulness to dexmedetomidine-induced unconsciousness, and early-stage recovery of consciousness. During hypnosis, brain regions become less functionally synchronized to each other; they attain a smaller number of functional configurations compared to wakefulness, and display functional connectivity patterns that are more similar to the underlying structural connectivity. Furthermore, cerebral blood flow significantly decreases across the whole brain, and less metabolically demanding low frequency fluctuations in the hemodynamic signal become more prominent. Collectively, loss of consciousness is accompanied by widespread connectivity and neurovascular changes in the brain, characteristic of less metabolically demanding dynamics.

Association of Anxiety with Uncinate Fasciculus Lesion Burden in Multiple Sclerosis.

Importance: Multiple sclerosis (MS) is an immune-mediated neurological disorder that affects 2.4 million people world-wide, and up to 60% experience anxiety. Objective: We investigated how anxiety in MS is associated with white matter lesion burden in the uncinate fasciculus (UF). Design: Retrospective case-control study of participants who received research-quality 3-tesla (3T) neuroimaging as part of MS clinical care from 2010-2018. Analyses were performed from June 1 st to September 30 th , 2024. Setting: Single-center academic medical specialty MS clinic. Participants: Participants were identified from the electronic medical record. All participants were diagnosed by an MS specialist and completed research-quality MRI at 3T. After excluding participants with poor image quality, 372 were stratified into three groups which were balanced for age and sex: 1) MS without anxiety (MS+noA, n=99); 2) MS with mild anxiety (MS+mildA, n=249); and 3) MS with severe anxiety (MS+severeA, n=24). Exposure: Anxiety diagnosis and anxiolytic medication. Main Outcome and Measure: We first evaluated whether MS+severeA patients had greater lesion burden in the UF than MS+noA. Next, we examined whether increasing anxiety severity was associated with greater UF lesion burden. Generalized additive models were employed, with the burden of lesions (e.g. proportion of fascicle impacted) within the UF as the outcome measure and sex and spline of age as covariates. Results: UF burden was higher in MS+severeA as compared to MS+noA (T=2.02, P=0.045, Cohen's f 2 =0.19). A dose-response effect was also found, where higher mean UF burden was associated with higher anxiety severity (T=2.08, P=0.038, Cohen's f 2 =0.10). Conclusions and Relevance: We demonstrate that overall lesion burden in UF was associated with the presence and severity of anxiety in patients with MS. Future studies linking white matter lesion burden in UF with treatment prognosis are warranted. KEY POINTS: Question: Are white matter lesions that impact the uncinate fasciculus (UF) associated with anxiety in patients with multiple sclerosis (MS)?Findings: This retrospective, case-control study of 372 patients with MS included 3 anxiety severity groups: 1) MS without anxiety (MS+noA, n=99); 2) MS with mild anxiety (MS+mildA, n=249); and 3) MS with severe anxiety (MS+severeA, n=24). We identified associations between anxiety and UF lesion burden. Specifically, we showed that MS+severeA had higher UF lesion burden than MS+noA, and worsening anxiety severity increased with greater UF burden.Meaning: Lesion burden in the UF may contribute to anxiety comorbidity in MS.

Reproducible Sex Differences in Personalized Functional Network Topography in Youth.

Background: A key step towards understanding psychiatric disorders that disproportionately impact female mental health is delineating the emergence of sex-specific patterns of brain organization at the critical transition from childhood to adolescence. Prior work suggests that individual differences in the spatial organization of functional brain networks across the cortex are associated with psychopathology and differ systematically by sex. Aims: We aimed to evaluate the impact of sex on the spatial organization of person-specific functional brain networks. Method: We leveraged person-specific atlases of functional brain networks defined using nonnegative matrix factorization in a sample of n = 6437 youths from the Adolescent Brain Cognitive Development Study. Across independent discovery and replication samples, we used generalized additive models to uncover associations between sex and the spatial layout ("topography") of personalized functional networks (PFNs). Next, we trained support vector machines to classify participants' sex from multivariate patterns of PFN topography. Finally, we leveraged transcriptomic data from the Allen Human Brain Atlas to evaluate spatial correlations between sex differences in PFN topography and gene expression. Results: Sex differences in PFN topography were greatest in association networks including the fronto-parietal, ventral attention, and default mode networks. Machine learning models trained on participants' PFNs were able to classify participant sex with high accuracy. Brain regions with the greatest sex differences in PFN topography were enriched in expression of X-linked genes as well as genes expressed in astrocytes and excitatory neurons. Conclusions: Sex differences in PFN topography are robust, replicate across large-scale samples of youth, and are associated with expression patterns of X-linked genes. These results suggest a potential contributor to the female-biased risk in depressive and anxiety disorders that emerge at the transition from childhood to adolescence.

Impact of insurance status on MRI phenotypes in MS.

BACKGROUND: Health insurance in the United States varies in coverage of essential diagnostic tests, therapies, and specialists. Health disparities between privately and publicly insured patients with MS have not been comprehensively assessed. The objective of this study is to evaluate the impact of public versus private insurance on longitudinal brain outcomes in MS. METHODS: Lesional, thalamic, and gray and white matter volumes were extracted from longitudinal MRI of 710 MS patients. Baseline volumes and atrophy rates of lesional, thalamic, and gray and white matter volumes were compared across insurance groups. RESULTS: After image quality assessment, 376 (284 private / 92 public), 638 (499 / 139), and 331 (250 / 81), patients were in MS lesion, thalamic, gray and white matter analyses respectively. Baseline lesion volume was higher for publicly insured patients but increased at a slightly higher rate in those privately insured (p = 0.01). Baseline gray matter measurements were lower for patients with public insurance, but thalamic (p < 0.01) and gray matter (p < 0.01) atrophy rates were slightly higher in the private insurance group. CONCLUSION: Insurance type was associated with lesion, thalamic, and gray matter volumes. The results suggest that patients with public insurance may present with more advanced disease.

Quantifying brain development in the HEALthy Brain and Child Development (HBCD) Study: The magnetic resonance imaging and spectroscopy protocol.

The HEALthy Brain and Child Development (HBCD) Study, a multi-site prospective longitudinal cohort study, will examine human brain, cognitive, behavioral, social, and emotional development beginning prenatally and planned through early childhood. The acquisition of multimodal magnetic resonance-based brain development data is central to the study's core protocol. However, application of Magnetic Resonance Imaging (MRI) methods in this population is complicated by technical challenges and difficulties of imaging in early life. Overcoming these challenges requires an innovative and harmonized approach, combining age-appropriate acquisition protocols together with specialized pediatric neuroimaging strategies. The HBCD MRI Working Group aimed to establish a core acquisition protocol for all 27 HBCD Study recruitment sites to measure brain structure, function, microstructure, and metabolites. Acquisition parameters of individual modalities have been matched across MRI scanner platforms for harmonized acquisitions and state-of-the-art technologies are employed to enable faster and motion-robust imaging. Here, we provide an overview of the HBCD MRI protocol, including decisions of individual modalities and preliminary data. The result will be an unparalleled resource for examining early neurodevelopment which enables the larger scientific community to assess normative trajectories from birth through childhood and to examine the genetic, biological, and environmental factors that help shape the developing brain.

Human adolescent brain similarity development is different for paralimbic versus neocortical zones.

Adolescent development of human brain structural and functional networks is increasingly recognized as fundamental to emergence of typical and atypical adult cognitive and emotional proodal magnetic resonance imaging (MRI) data collected from N [Formula: see text] 300 healthy adolescents (51%; female; 14 to 26 y) each scanned repeatedly in an accelerated longitudinal design, to provide an analyzable dataset of 469 structural scans and 448 functional MRI scans. We estimated the morphometric similarity between each possible pair of 358 cortical areas on a feature vector comprising six macro- and microstructural MRI metrics, resulting in a morphometric similarity network (MSN) for each scan. Over the course of adolescence, we found that morphometric similarity increased in paralimbic cortical areas, e.g., insula and cingulate cortex, but generally decreased in neocortical areas, and these results were replicated in an independent developmental MRI cohort (N [Formula: see text] 304). Increasing hubness of paralimbic nodes in MSNs was associated with increased strength of coupling between their morphometric similarity and functional connectivity. Decreasing hubness of neocortical nodes in MSNs was associated with reduced strength of structure-function coupling and increasingly diverse functional connections in the corresponding fMRI networks. Neocortical areas became more structurally differentiated and more functionally integrative in a metabolically expensive process linked to cortical thinning and myelination, whereas paralimbic areas specialized for affective and interoceptive functions became less differentiated, as hypothetically predicted by a developmental transition from periallocortical to proisocortical organization of the cortex. Cytoarchitectonically distinct zones of the human cortex undergo distinct neurodevelopmental programs during typical adolescence.

A network control theory pipeline for studying the dynamics of the structural connectome.

Network control theory (NCT) is a simple and powerful tool for studying how network topology informs and constrains the dynamics of a system. Compared to other structure-function coupling approaches, the strength of NCT lies in its capacity to predict the patterns of external control signals that may alter the dynamics of a system in a desired way. An interesting development for NCT in the neuroscience field is its application to study behavior and mental health symptoms. To date, NCT has been validated to study different aspects of the human structural connectome. NCT outputs can be monitored throughout developmental stages to study the effects of connectome topology on neural dynamics and, separately, to test the coherence of empirical datasets with brain function and stimulation. Here, we provide a comprehensive pipeline for applying NCT to structural connectomes by following two procedures. The main procedure focuses on computing the control energy associated with the transitions between specific neural activity states. The second procedure focuses on computing average controllability, which indexes nodes' general capacity to control the dynamics of the system. We provide recommendations for comparing NCT outputs against null network models, and we further support this approach with a Python-based software package called 'network control theory for python'. The procedures in this protocol are appropriate for users with a background in network neuroscience and experience in dynamical systems theory.

Corticospinal premotor fibers facilitate complex motor control after stroke.

OBJECTIVE: The corticospinal tract (CST) is considered the most important motor output pathway comprising fibers from the primary motor cortex (M1) and various premotor areas. Damage to its descending fibers after stroke commonly leads to motor impairment. While premotor areas are thought to critically support motor recovery after stroke, the functional role of their corticospinal output for different aspects of post-stroke motor control remains poorly understood. METHODS: We assessed the differential role of CST fibers originating from premotor areas and M1 in the control of basal (single-joint muscle synergies and strength) and complex motor control (involving inter-joint coordination and visuomotor integration) using a novel diffusion imaging approach in chronic stroke patients. RESULTS: While M1 sub-tract anisotropy was positively correlated with basal and complex motor skills, anisotropy of PMd, PMv, and SMA sub-tracts was exclusively associated with complex motor tasks. Interestingly, patients featuring persistent motor deficits showed an additional positive association between premotor sub-tract integrity and basal motor control. INTERPRETATION: While descending M1 output seems to be a prerequisite for any form of upper limb movements, complex motor skills critically depend on output from premotor areas after stroke. The additional involvement of premotor tracts in basal motor control in patients with persistent deficits emphasizes their compensatory capacity in post-stroke motor control. In summary, our findings highlight the pivotal role of descending corticospinal output from premotor areas for motor control after stroke, which thus serve as prime candidates for future interventions to amplify motor recovery.

A sensorimotor-association axis of thalamocortical connection development.

Human cortical development follows a sensorimotor-to-association sequence during childhood and adolescence1-6. The brain's capacity to enact this sequence over decades indicates that it relies on intrinsic mechanisms to regulate inter-regional differences in the timing of cortical maturation, yet regulators of human developmental chronology are not well understood. Given evidence from animal models that thalamic axons modulate windows of cortical plasticity7-12, here we evaluate the overarching hypothesis that structural connections between the thalamus and cortex help to coordinate cortical maturational heterochronicity during youth. We first introduce, cortically annotate, and anatomically validate a new atlas of human thalamocortical connections using diffusion tractography. By applying this atlas to three independent youth datasets (ages 8-23 years; total N = 2,676), we reproducibly demonstrate that thalamocortical connections develop along a maturational gradient that aligns with the cortex's sensorimotor-association axis. Associative cortical regions with thalamic connections that take longest to mature exhibit protracted expression of neurochemical, structural, and functional markers indicative of higher circuit plasticity as well as heightened environmental sensitivity. This work highlights a central role for the thalamus in the orchestration of hierarchically organized and environmentally sensitive windows of cortical developmental malleability.

Functional connectivity development along the sensorimotor-association axis enhances the cortical hierarchy.

Human cortical maturation has been posited to be organized along the sensorimotor-association axis, a hierarchical axis of brain organization that spans from unimodal sensorimotor cortices to transmodal association cortices. Here, we investigate the hypothesis that the development of functional connectivity during childhood through adolescence conforms to the cortical hierarchy defined by the sensorimotor-association axis. We tested this pre-registered hypothesis in four large-scale, independent datasets (total n = 3355; ages 5-23 years): the Philadelphia Neurodevelopmental Cohort (n = 1207), Nathan Kline Institute-Rockland Sample (n = 397), Human Connectome Project: Development (n = 625), and Healthy Brain Network (n = 1126). Across datasets, the development of functional connectivity systematically varied along the sensorimotor-association axis. Connectivity in sensorimotor regions increased, whereas connectivity in association cortices declined, refining and reinforcing the cortical hierarchy. These consistent and generalizable results establish that the sensorimotor-association axis of cortical organization encodes the dominant pattern of functional connectivity development.

A practical evaluation of measures derived from compressed sensing diffusion spectrum imaging.

Diffusion Spectrum Imaging (DSI) using dense Cartesian sampling of q-space has been shown to provide important advantages for modeling complex white matter architecture. However, its adoption has been limited by the lengthy acquisition time required. Sparser sampling of q-space combined with compressed sensing (CS) reconstruction techniques has been proposed as a way to reduce the scan time of DSI acquisitions. However prior studies have mainly evaluated CS-DSI in post-mortem or non-human data. At present, the capacity for CS-DSI to provide accurate and reliable measures of white matter anatomy and microstructure in the living human brain remains unclear. We evaluated the accuracy and inter-scan reliability of 6 different CS-DSI schemes that provided up to 80% reductions in scan time compared to a full DSI scheme. We capitalized on a dataset of 26 participants who were scanned over eight independent sessions using a full DSI scheme. From this full DSI scheme, we subsampled images to create a range of CS-DSI images. This allowed us to compare the accuracy and inter-scan reliability of derived measures of white matter structure (bundle segmentation, voxel-wise scalar maps) produced by the CS-DSI and the full DSI schemes. We found that CS-DSI estimates of both bundle segmentations and voxel-wise scalars were nearly as accurate and reliable as those generated by the full DSI scheme. Moreover, we found that the accuracy and reliability of CS-DSI was higher in white matter bundles that were more reliably segmented by the full DSI scheme. As a final step, we replicated the accuracy of CS-DSI in a prospectively acquired dataset (n = 20, scanned once). Together, these results illustrate the utility of CS-DSI for reliably delineating in vivo white matter architecture in a fraction of the scan time, underscoring its promise for both clinical and research applications.

Diffusion MRI head motion correction methods are highly accurate but impacted by denoising and sampling scheme.

Head motion correction is particularly challenging in diffusion-weighted MRI (dMRI) scans due to the dramatic changes in image contrast at different gradient strengths and directions. Head motion correction is typically performed using a Gaussian Process model implemented in FSL's Eddy. Recently, the 3dSHORE-based SHORELine method was introduced that does not require shell-based acquisitions, but it has not been previously benchmarked. Here we perform a comprehensive evaluation of both methods on realistic simulations of a software fiber phantom that provides known ground-truth head motion. We demonstrate that both methods perform remarkably well, but that performance can be impacted by sampling scheme and the extent of head motion and the denoising strategy applied before head motion correction. Furthermore, we find Eddy benefits from denoising the data first with MP-PCA. In sum, we provide the most extensive known benchmarking of dMRI head motion correction, together with extensive simulation data and a reproducible workflow. PRACTITIONER POINTS: Both Eddy and SHORELine head motion correction methods performed quite well on a large variety of simulated data. Denoising with MP-PCA can improve head motion correction performance when Eddy is used. SHORELine effectively corrects motion in non-shelled diffusion spectrum imaging data.

XCP-D: A Robust Pipeline for the post-processing of fMRI data.

Functional neuroimaging is an essential tool for neuroscience research. Pre-processing pipelines produce standardized, minimally pre-processed data to support a range of potential analyses. However, post-processing is not similarly standardized. While several options for post-processing exist, they tend not to support output from disparate pre-processing pipelines, may have limited documentation, and may not follow BIDS best practices. Here we present XCP-D, which presents a solution to these issues. XCP-D is a collaborative effort between PennLINC at the University of Pennsylvania and the DCAN lab at the University at Minnesota. XCP-D uses an open development model on GitHub and incorporates continuous integration testing; it is distributed as a Docker container or Singularity image. XCP-D generates denoised BOLD images and functional derivatives from resting-state data in either NifTI or CIFTI files, following pre-processing with fMRIPrep, HCP, and ABCD-BIDS pipelines. Even prior to its official release, XCP-D has been downloaded >3,000 times from DockerHub. Together, XCP-D facilitates robust, scalable, and reproducible post-processing of fMRI data.

Development of white matter fiber covariance networks supports executive function in youth.

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.

Mapping the Relationship of White Matter Lesions to Depression in Multiple Sclerosis.

BACKGROUND: Multiple sclerosis (MS) is an immune-mediated neurological disorder, and up to 50% of patients experience depression. We investigated how white matter network disruption is related to depression in MS. METHODS: Using electronic health records, 380 participants with MS were identified. Depressed individuals (MS+Depression group; n = 232) included persons who had an ICD-10 depression diagnosis, had a prescription for antidepressant medication, or screened positive via Patient Health Questionnaire (PHQ)-2 or PHQ-9. Age- and sex-matched nondepressed individuals with MS (MS-Depression group; n = 148) included persons who had no prior depression diagnosis, had no psychiatric medication prescriptions, and were asymptomatic on PHQ-2 or PHQ-9. Research-quality 3T structural magnetic resonance imaging was obtained as part of routine care. We first evaluated whether lesions were preferentially located within the depression network compared with other brain regions. Next, we examined if MS+Depression patients had greater lesion burden and if this was driven by lesions in the depression network. Primary outcome measures were the burden of lesions (e.g., impacted fascicles) within a network and across the brain. RESULTS: MS lesions preferentially affected fascicles within versus outside the depression network (ß = 0.09, 95% CI = 0.08 to 0.10, p < .001). MS+Depression patients had more lesion burden (ß = 0.06, 95% CI = 0.01 to 0.10, p = .015); this was driven by lesions within the depression network (ß = 0.02, 95% CI = 0.003 to 0.040, p = .020). CONCLUSIONS: We demonstrated that lesion location and burden may contribute to depression comorbidity in MS. MS lesions disproportionately impacted fascicles in the depression network. MS+Depression patients had more disease than MS-Depression patients, which was driven by disease within the depression network. Future studies relating lesion location to personalized depression interventions are warranted.

Myelination and excitation-inhibition balance synergistically shape structure-function coupling across the human cortex.

Recent work has demonstrated that the relationship between structural and functional connectivity varies regionally across the human brain, with reduced coupling emerging along the sensory-association cortical hierarchy. The biological underpinnings driving this expression, however, remain largely unknown. Here, we postulate that intracortical myelination and excitation-inhibition (EI) balance mediate the heterogeneous expression of structure-function coupling (SFC) and its temporal variance across the cortical hierarchy. We employ atlas- and voxel-based connectivity approaches to analyze neuroimaging data acquired from two groups of healthy participants. Our findings are consistent across six complementary processing pipelines: 1) SFC and its temporal variance respectively decrease and increase across the unimodal-transmodal and granular-agranular gradients; 2) increased myelination and lower EI-ratio are associated with more rigid SFC and restricted moment-to-moment SFC fluctuations; 3) a gradual shift from EI-ratio to myelination as the principal predictor of SFC occurs when traversing from granular to agranular cortical regions. Collectively, our work delivers a framework to conceptualize structure-function relationships in the human brain, paving the way for an improved understanding of how demyelination and/or EI-imbalances induce reorganization in brain disorders.

Using network control theory to study the dynamics of the structural connectome.

Network control theory (NCT) is a simple and powerful tool for studying how network topology informs and constrains dynamics. Compared to other structure-function coupling approaches, the strength of NCT lies in its capacity to predict the patterns of external control signals that may alter dynamics in a desired way. We have extensively developed and validated the application of NCT to the human structural connectome. Through these efforts, we have studied (i) how different aspects of connectome topology affect neural dynamics, (ii) whether NCT outputs cohere with empirical data on brain function and stimulation, and (iii) how NCT outputs vary across development and correlate with behavior and mental health symptoms. In this protocol, we introduce a framework for applying NCT to structural connectomes following two main pathways. Our primary pathway focuses on computing the control energy associated with transitioning between specific neural activity states. Our second pathway focuses on computing average controllability, which indexes nodes' general capacity to control dynamics. We also provide recommendations for comparing NCT outputs against null network models. Finally, we support this protocol with a Python-based software package called network control theory for python (nctpy).

ReeBundle: A Method for Topological Modeling of White Matter Pathways Using Diffusion MRI.

Tractography can generate millions of complex curvilinear fibers (streamlines) in 3D that exhibit the geometry of white matter pathways in the brain. Common approaches to analyzing white matter connectivity are based on adjacency matrices that quantify connection strength but do not account for any topological information. A critical element in neurological and developmental disorders is the topological deterioration and irregularities in streamlines. In this paper, we propose a novel Reeb graph-based method "ReeBundle" that efficiently encodes the topology and geometry of white matter fibers. Given the trajectories of neuronal fiber pathways (neuroanatomical bundle), we re-bundle the streamlines by modeling their spatial evolution to capture geometrically significant events (akin to a fingerprint). ReeBundle parameters control the granularity of the model and handle the presence of improbable streamlines commonly produced by tractography. Further, we propose a new Reeb graph-based distance metric that quantifies topological differences for automated quality control and bundle comparison. We show the practical usage of our method using two datasets: (1) For International Society for Magnetic Resonance in Medicine (ISMRM) dataset, ReeBundle handles the morphology of the white matter tract configurations due to branching and local ambiguities in complicated bundle tracts like anterior and posterior commissures; (2) For the longitudinal repeated measures in the Cognitive Resilience and Sleep History (CRASH) dataset, repeated scans of a given subject acquired weeks apart lead to provably similar Reeb graphs that differ significantly from other subjects, thus highlighting ReeBundle's potential for clinical fingerprinting of brain regions.

A reproducible and generalizable software workflow for analysis of large-scale neuroimaging data collections using BIDS Apps.

Neuroimaging research faces a crisis of reproducibility. With massive sample sizes and greater data complexity, this problem becomes more acute. Software that operates on imaging data defined using the Brain Imaging Data Structure (BIDS) - BIDS Apps - have provided a substantial advance. However, even using BIDS Apps, a full audit trail of data processing is a necessary prerequisite for fully reproducible research. Obtaining a faithful record of the audit trail is challenging - especially for large datasets. Recently, the FAIRly big framework was introduced as a way to facilitate reproducible processing of large-scale data by leveraging DataLad - a version control system for data management. However, the current implementation of this framework was more of a proof of concept, and could not be immediately reused by other investigators for different use cases. Here we introduce the BIDS App Bootstrap (BABS), a user-friendly and generalizable Python package for reproducible image processing at scale. BABS facilitates the reproducible application of BIDS Apps to large-scale datasets. Leveraging DataLad and the FAIRly big framework, BABS tracks the full audit trail of data processing in a scalable way by automatically preparing all scripts necessary for data processing and version tracking on high performance computing (HPC) systems. Currently, BABS supports jobs submissions and audits on Sun Grid Engine (SGE) and Slurm HPCs with a parsimonious set of programs. To demonstrate its scalability, we applied BABS to data from the Healthy Brain Network (HBN; n=2,565). Taken together, BABS allows reproducible and scalable image processing and is broadly extensible via an open-source development model.

Transcriptomic analysis reveals responses to a polluted sediment in the Mediterranean copepod Acartia clausi.

Marine sediments are regarded as sinks for several classes of contaminants. Characterization and effects of sediments on marine biota now require a multidisciplinary approach, which includes chemical and ecotoxicological analyses and molecular biomarkers. Here, a gene expression study was performed to measure the response of adult females of the Mediterranean copepod Acartia clausi to elutriates of polluted sediments (containing high concentrations of polycyclic aromatic hydrocarbons, PAHs, and heavy metals) from an industrial area in the Southern Tyrrhenian Sea (Bagnoli-Coroglio). Functional annotation of the A. clausi transcriptome generated as reference here, showed a good quality of the assembly and great homology with other copepod and crustacean sequences in public databases. This is one of the few available transcriptomic resources for this widespread copepod species of great ecological relevance in temperate coastal areas. Differential expression analysis between females exposed to the elutriate and those in control seawater identified 1000 differentially expressed genes, of which 743 up- and 257 down-regulated. Within the up-regulated genes, the most represented functions were related to proteolysis (lysosomal protease, peptidase, cathepsin), response to stress and detoxification (heat-shock protein, superoxide dismutase, glutathione-S-transferase, cytochrome P450), and cytoskeleton structure (α- and ß-tubulin). Down-regulated genes were mostly involved with ribosome structure (ribosomal proteins) and DNA binding (histone proteins, transcription factors). Overall, these results suggest that processes such as transcription, translation, protein degradation, metabolism of biomolecules, reproduction, and xenobiotic detoxification were altered in the copepod in response to polluted elutriates. In conclusion, our results contribute to gaining information on the transcriptomic responses of copepods to polluted sediments. They will also prompt the selection of genes of interest to be used as biomarkers of exposure to PAHs and heavy metals in molecular toxicology studies on copepods, and in general, in comparative functional genomic studies on marine zooplankton.