ASLPrep: a platform for processing of arterial spin labeled MRI and quantification of regional brain perfusion.

Arterial spin labeled (ASL) magnetic resonance imaging (MRI) is the primary method for noninvasively measuring regional brain perfusion in humans. We introduce ASLPrep, a suite of software pipelines that ensure the reproducible and generalizable processing of ASL MRI data.

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.

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.

QSIPrep: an integrative platform for preprocessing and reconstructing diffusion MRI data.

Diffusion-weighted magnetic resonance imaging (dMRI) is the primary method for noninvasively studying the organization of white matter in the human brain. Here we introduce QSIPrep, an integrative software platform for the processing of diffusion images that is compatible with nearly all dMRI sampling schemes. Drawing on a diverse set of software suites to capitalize on their complementary strengths, QSIPrep facilitates the implementation of best practices for processing of diffusion images.

Interhemispheric Structural Connectivity Underlies Motor Recovery after Stroke.

OBJECTIVE: Although ample evidence highlights that the ipsilesional corticospinal tract (CST) plays a crucial role in motor recovery after stroke, studies on cortico-cortical motor connections remain scarce and provide inconclusive results. Given their unique potential to serve as structural reserve enabling motor network reorganization, the question arises whether cortico-cortical connections may facilitate motor control depending on CST damage. METHODS: Diffusion spectrum imaging (DSI) and a novel compartment-wise analysis approach were used to quantify structural connectivity between bilateral cortical core motor regions in chronic stroke patients. Basal and complex motor control were differentially assessed. RESULTS: Both basal and complex motor performance were correlated with structural connectivity between bilateral premotor areas and ipsilesional primary motor cortex (M1) as well as interhemispheric M1 to M1 connectivity. Whereas complex motor skills depended on CST integrity, a strong association between M1 to M1 connectivity and basal motor control was observed independent of CST integrity especially in patients who underwent substantial motor recovery. Harnessing the informational wealth of cortico-cortical connectivity facilitated the explanation of both basal and complex motor control. INTERPRETATION: We demonstrate for the first time that distinct aspects of cortical structural reserve enable basal and complex motor control after stroke. In particular, recovery of basal motor control may be supported via an alternative route through contralesional M1 and non-crossing fibers of the contralesional CST. Our findings help to explain previous conflicting interpretations regarding the functional role of the contralesional M1 and highlight the potential of cortico-cortical structural connectivity as a future biomarker for motor recovery post-stroke. ANN NEUROL 2023;94:785-797.

ModelArray: An R package for statistical analysis of fixel-wise data.

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.

Development of structure-function coupling in human brain networks during youth.

The protracted development of structural and functional brain connectivity within distributed association networks coincides with improvements in higher-order cognitive processes such as executive function. However, it remains unclear how white-matter architecture develops during youth to directly support coordinated neural activity. Here, we characterize the development of structure-function coupling using diffusion-weighted imaging and n-back functional MRI data in a sample of 727 individuals (ages 8 to 23 y). We found that spatial variability in structure-function coupling aligned with cortical hierarchies of functional specialization and evolutionary expansion. Furthermore, hierarchy-dependent age effects on structure-function coupling localized to transmodal cortex in both cross-sectional data and a subset of participants with longitudinal data (n = 294). Moreover, structure-function coupling in rostrolateral prefrontal cortex was associated with executive performance and partially mediated age-related improvements in executive function. Together, these findings delineate a critical dimension of adolescent brain development, whereby the coupling between structural and functional connectivity remodels to support functional specialization and cognition.

Retinal ganglion cell endowment is correlated with optic tract fiber cross section, not density.

There is substantial variation between healthy individuals in the number of retinal ganglion cells (RGC) in the eye, with commensurate variation in the number of axons in the optic tracts. Fixel-based analysis of diffusion MR produces estimates of fiber density (FD) and cross section (FC). Using these fixel measurements along with retinal imaging, we asked if individual differences in RGC tissue volume are correlated with individual differences in FD and FC measurements obtained from the optic tracts, and subsequent structures along the cortical visual pathway. We find that RGC endowment is correlated with optic tract FC, but not with FD. RGC volume had a decreasing relationship with measurements from subsequent regions of the visual system (LGN volume, optic radiation FC/FD, and V1 surface area). However, we also found that the variations in each visual area were correlated with the variations in its immediately adjacent visual structure. We only observed these serial correlations when FC is used as the measure of interest for the optic tract and radiations, but no significant relationship was found when FD represented these white matter structures. From these results, we conclude that the variations in RGC endowment, LGN volume, and V1 surface area are better predicted by the overall cross section of the optic tract and optic radiations as compared to the intra-axonal restricted signal component of these white matter pathways. Additionally, the presence of significant correlations between adjacent, but not distant, anatomical structures suggests that there are multiple, local sources of anatomical variation along the visual pathway.

Curation of BIDS (CuBIDS): A workflow and software package for streamlining reproducible curation of large BIDS datasets.

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.

Physiological acclimatization in high-latitude zooplankton.

How individual organisms adapt to nonoptimal conditions through physiological acclimatization is central to predicting the consequences of unusual abiotic and biotic conditions such as those produced by marine heat waves. The Northeast Pacific, including the Gulf of Alaska, experienced an extreme warming event (2014-2016, "The Blob") that affected all trophic levels and led to large-scale changes in the community. The marine copepod Neocalanus flemingeri is a key member of the subarctic Pacific pelagic ecosystem. During the spring phytoplankton bloom this copepod builds substantial lipid stores as it prepares for its nonfeeding adult phase. A 3-year comparison of gene expression profiles of copepods collected in Prince William Sound in the Gulf of Alaska between 2015 and 2017 included two high-temperature years (2015 and 2016) and one year with very low phytoplankton abundances (2016). The largest differences in gene expression were between high and low chlorophyll years, and not between warm and cool years. The observed gene expression patterns were indicative of physiological acclimatization. The predominant signal in 2016 was the down-regulation of genes involved in glycolysis and its incoming pathways, consistent with the modulation of metabolic rates in response to prolonged low food conditions. Despite the down-regulation of genes involved in metabolism, there was no evidence of suppression of protein synthesis based on gene expression or behavioural activity. Genes involved in muscle function were up-regulated, and the copepods were actively swimming and responsive to stimuli at collection. However, genes involved in fatty acid metabolism were down-regulated in 2016, suggesting reduced lipid accumulation.

Post-diapause transcriptomic restarts: insight from a high-latitude copepod.

BACKGROUND: Diapause is a seasonal dormancy that allows organisms to survive unfavorable conditions and optimizes the timing of reproduction and growth. Emergence from diapause reverses the state of arrested development and metabolic suppression returning the organism to an active state. The physiological mechanisms that regulate the transition from diapause to post-diapause are still unknown. In this study, this transition has been characterized for the sub-arctic calanoid copepod Neocalanus flemingeri, a key crustacean zooplankter that supports the highly productive North Pacific fisheries. Transcriptional profiling of females, determined over a two-week time series starting with diapausing females collected from > 400 m depth, characterized the molecular mechanisms that regulate the post-diapause trajectory. RESULTS: A complex set of transitions in relative gene expression defined the transcriptomic changes from diapause to post-diapause. Despite low temperatures (5-6 °C), the switch from a "diapause" to a "post-diapause" transcriptional profile occurred within 12 h of the termination stimulus. Transcriptional changes signaling the end of diapause were activated within one-hour post collection and included the up-regulation of genes involved in the 20E cascade pathway, the TCA cycle and RNA metabolism in combination with the down-regulation of genes associated with chromatin silencing. By 12 h, females exhibited a post-diapause phenotype characterized by the up-regulation of genes involved in cell division, cell differentiation and multiple developmental processes. By seven days post collection, the reproductive program was fully activated as indicated by up-regulation of genes involved in oogenesis and energy metabolism, processes that were enriched among the differentially expressed genes. CONCLUSIONS: The analysis revealed a finely structured, precisely orchestrated sequence of transcriptional changes that led to rapid changes in the activation of biological processes paving the way to the successful completion of the reproductive program. Our findings lead to new hypotheses related to potentially universal mechanisms that terminate diapause before an organism can resume its developmental program.

Spatial coherence of oriented white matter microstructure: Applications to white matter regions associated with genetic similarity.

We present a method to discover differences between populations with respect to the spatial coherence of their oriented white matter microstructure in arbitrarily shaped white matter regions. This method is applied to diffusion MRI scans of a subset of the Human Connectome Project dataset: 57 pairs of monozygotic and 52 pairs of dizygotic twins. After controlling for morphological similarity between twins, we identify 3.7% of all white matter as being associated with genetic similarity (35.1 k voxels, p<10-4, false discovery rate 1.5%), 75% of which spatially clusters into twenty-two contiguous white matter regions. Furthermore, we show that the orientation similarity within these regions generalizes to a subset of 47 pairs of non-twin siblings, and show that these siblings are on average as similar as dizygotic twins. The regions are located in deep white matter including the superior longitudinal fasciculus, the optic radiations, the middle cerebellar peduncle, the corticospinal tract, and within the anterior temporal lobe, as well as the cerebellum, brain stem, and amygdalae. These results extend previous work using undirected fractional anisotrophy for measuring putative heritable influences in white matter. Our multidirectional extension better accounts for crossing fiber connections within voxels. This bottom up approach has at its basis a novel measurement of coherence within neighboring voxel dyads between subjects, and avoids some of the fundamental ambiguities encountered with tractographic approaches to white matter analysis that estimate global connectivity.

Analytic tractography: A closed-form solution for estimating local white matter connectivity with diffusion MRI.

White matter structures composed of myelinated axons in the living human brain are primarily studied by diffusion-weighted MRI (dMRI). These long-range projections are typically characterized in a two-step process: dMRI signal is used to estimate the orientation of axon segments within each voxel, then these local orientations are linked together to estimate the spatial extent of putative white matter bundles. Tractography, the process of tracing bundles across voxels, either requires computationally expensive (probabilistic) simulations to model uncertainty in fiber orientation or ignores it completely (deterministic). Furthermore, simulation necessarily generates a finite number of trajectories, introducing "simulation error" to trajectory estimates. Here we introduce a method to analytically (via a closed-form solution) take an orientation distribution function (ODF) from each voxel and calculate the probabilities that a trajectory projects from a voxel into each directly adjacent voxels. We validate our method by demonstrating experimentally that probabilistic simulations converge to our analytically computed transition probabilities at the voxel level as the number of simulated seeds increases. We then show that our method accurately calculates the ground-truth transition probabilities from a publicly available phantom dataset. As a demonstration, we incorporate our analytic method for voxel transition probabilities into the Voxel Graph framework, creating a quantitative framework for assessing white matter structure, which we call "analytic tractography". The long-range connectivity problem is reduced to finding paths in a graph whose adjacency structure reflects voxel-to-voxel analytic transition probabilities. We demonstrate that this approach performs comparably to the current most widely-used probabilistic and deterministic approaches at a fraction of the computational cost. We also demonstrate that analytic tractography works on multiple diffusion sampling schemes, reconstruction method or parameters used to define paths. Open source software compatible with popular dMRI reconstruction software is provided.

Optimal trajectories of brain state transitions.

The complexity of neural dynamics stems in part from the complexity of the underlying anatomy. Yet how white matter structure constrains how the brain transitions from one cognitive state to another remains unknown. Here we address this question by drawing on recent advances in network control theory to model the underlying mechanisms of brain state transitions as elicited by the collective control of region sets. We find that previously identified attention and executive control systems are poised to affect a broad array of state transitions that cannot easily be classified by traditional engineering-based notions of control. This theoretical versatility comes with a vulnerability to injury. In patients with mild traumatic brain injury, we observe a loss of specificity in putative control processes, suggesting greater susceptibility to neurophysiological noise. These results offer fundamental insights into the mechanisms driving brain state transitions in healthy cognition and their alteration following injury.

Stimulation-Based Control of Dynamic Brain Networks.

The ability to modulate brain states using targeted stimulation is increasingly being employed to treat neurological disorders and to enhance human performance. Despite the growing interest in brain stimulation as a form of neuromodulation, much remains unknown about the network-level impact of these focal perturbations. To study the system wide impact of regional stimulation, we employ a data-driven computational model of nonlinear brain dynamics to systematically explore the effects of targeted stimulation. Validating predictions from network control theory, we uncover the relationship between regional controllability and the focal versus global impact of stimulation, and we relate these findings to differences in the underlying network architecture. Finally, by mapping brain regions to cognitive systems, we observe that the default mode system imparts large global change despite being highly constrained by structural connectivity. This work forms an important step towards the development of personalized stimulation protocols for medical treatment or performance enhancement.

Prediction of a neuropeptidome for the eyestalk ganglia of the lobster Homarus americanus using a tissue-specific de novo assembled transcriptome.

In silico transcriptome mining is a powerful tool for crustacean peptidome prediction. Using homology-based BLAST searches and a simple bioinformatics workflow, large peptidomes have recently been predicted for a variety of crustaceans, including the lobster, Homarus americanus. Interestingly, no in silico studies have been conducted on the eyestalk ganglia (lamina ganglionaris, medulla externa, medulla interna and medulla terminalis) of the lobster, although the eyestalk is the location of a major neuroendocrine complex, i.e., the X-organ-sinus gland system. Here, an H. americanus eyestalk ganglia-specific transcriptome was produced using the de novo assembler Trinity. This transcriptome was generated from 130,973,220 Illumina reads and consists of 147,542 unique contigs. Eighty-nine neuropeptide-encoding transcripts were identified from this dataset, allowing for the deduction of 62 distinct pre/preprohormones. Two hundred sixty-two neuropeptides were predicted from this set of precursors; the peptides include members of the adipokinetic hormone-corazonin-like peptide, allatostatin A, allatostatin B, allatostatin C, bursicon α, CCHamide, corazonin, crustacean cardioactive peptide, crustacean hyperglycemic hormone (CHH), CHH precursor-related peptide, diuretic hormone 31, diuretic hormone 44, eclosion hormone, elevenin, FMRFamide-like peptide, glycoprotein hormone α2, glycoprotein hormone ß5, GSEFLamide, intocin, leucokinin, molt-inhibiting hormone, myosuppressin, neuroparsin, neuropeptide F, orcokinin, orcomyotropin, pigment dispersing hormone, proctolin, pyrokinin, red pigment concentrating hormone, RYamide, short neuropeptide F, SIFamide, sulfakinin, tachykinin-related peptide and trissin families. The predicted peptides expand the H. americanus eyestalk ganglia neuropeptidome approximately 7-fold, and include 78 peptides new to the lobster. The transcriptome and predicted neuropeptidome described here provide new resources for investigating peptidergic signaling within/from the lobster eyestalk ganglia.

Feature interactions enable decoding of sensorimotor transformations for goal-directed movement.

Neurophysiology and neuroimaging evidence shows that the brain represents multiple environmental and body-related features to compute transformations from sensory input to motor output. However, it is unclear how these features interact during goal-directed movement. To investigate this issue, we examined the representations of sensory and motor features of human hand movements within the left-hemisphere motor network. In a rapid event-related fMRI design, we measured cortical activity as participants performed right-handed movements at the wrist, with either of two postures and two amplitudes, to move a cursor to targets at different locations. Using a multivoxel analysis technique with rigorous generalization tests, we reliably distinguished representations of task-related features (primarily target location, movement direction, and posture) in multiple regions. In particular, we identified an interaction between target location and movement direction in the superior parietal lobule, which may underlie a transformation from the location of the target in space to a movement vector. In addition, we found an influence of posture on primary motor, premotor, and parietal regions. Together, these results reveal the complex interactions between different sensory and motor features that drive the computation of sensorimotor transformations.

Regional white matter variation associated with domain-specific metacognitive accuracy.

The neural mechanisms that mediate metacognitive ability (the capacity to accurately reflect on one's own cognition and experience) remain poorly understood. An important question is whether metacognitive capacity is a domain-general skill supported by a core neuroanatomical substrate or whether regionally specific neural structures underlie accurate reflection in different cognitive domains. Providing preliminary support for the latter possibility, recent findings have shown that individual differences in metacognitive ability in the domains of memory and perception are related to variation in distinct gray matter volume and resting-state functional connectivity. The current investigation sought to build on these findings by evaluating how metacognitive ability in these domains is related to variation in white matter microstructure. We quantified metacognitive ability across memory and perception domains and used diffusion spectrum imaging to examine the relation between high-resolution measurements of white matter microstructure and individual differences in metacognitive accuracy in each domain. We found that metacognitive accuracy for perceptual decisions and memory were uncorrelated across individuals and that metacognitive accuracy in each domain was related to variation in white matter microstructure in distinct brain areas. Metacognitive accuracy for perceptual decisions was associated with increased diffusion anisotropy in white matter underlying the ACC, whereas metacognitive accuracy for memory retrieval was associated with increased diffusion anisotropy in the white matter extending into the inferior parietal lobule. Together, these results extend previous findings linking metacognitive ability in the domains of perception and memory to variation in distinct brain structures and connections.

Identification and developmental expression of the enzymes responsible for dopamine, histamine, octopamine and serotonin biosynthesis in the copepod crustacean Calanus finmarchicus.

Neurochemicals are likely to play key roles in physiological/behavioral control in the copepod crustacean Calanus finmarchicus, the biomass dominant zooplankton for much of the North Atlantic Ocean. Previously, a de novo assembled transcriptome consisting of 206,041 unique sequences was used to characterize the peptidergic signaling systems of Calanus. Here, this assembly was mined for transcripts encoding enzymes involved in amine biosynthesis. Using known Drosophila melanogaster proteins as templates, transcripts encoding putative Calanus homologs of tryptophan-phenylalanine hydroxylase (dopamine, octopamine and serotonin biosynthesis), tyrosine hydroxylase (dopamine biosynthesis), DOPA decarboxylase (dopamine and serotonin biosynthesis), histidine decarboxylase (histamine biosynthesis), tyrosine decarboxylase (octopamine biosynthesis), tyramine ß-hydroxylase (octopamine biosynthesis) and tryptophan hydroxylase (serotonin biosynthesis) were identified. Reverse BLAST and domain analyses show that the proteins deduced from these transcripts possess sequence homology to and the structural hallmarks of their respective enzyme families. Developmental profiling revealed a remarkably consistent pattern of expression for all transcripts, with the highest levels of expression typically seen in the early nauplius and early copepodite. These expression patterns suggest roles for amines during development, particularly in the metamorphic transitions from embryo to nauplius and from nauplius to copepodite. Taken collectively, the data presented here lay a strong foundation for future gene-based studies of aminergic signaling in this and other copepod species, in particular assessment of the roles they may play in developmental control.