Evaluating phase synchronization methods in fMRI: A comparison study and new approaches.

In recent years there has been growing interest in measuring time-varying functional connectivity between different brain regions using resting-state functional magnetic resonance imaging (rs-fMRI) data. One way to assess the relationship between signals from different brain regions is to measure their phase synchronization (PS) across time. There are several ways to perform such analyses, and we compare methods that utilize a PS metric together with a sliding window, referred to here as windowed phase synchronization (WPS), with those that directly measure the instantaneous phase synchronization (IPS). In particular, IPS has recently gained popularity as it offers single time-point resolution of time-resolved fMRI connectivity. In this paper, we discuss the underlying assumptions required for performing PS analyses and emphasize the importance of band-pass filtering the data to obtain valid results. Further, we contrast this approach with the use of Empirical Mode Decomposition (EMD) to achieve similar goals. We review various methods for evaluating PS and introduce a new approach within the IPS framework denoted the cosine of the relative phase (CRP). We contrast methods through a series of simulations and application to rs-fMRI data. Our results indicate that CRP outperforms other tested methods and overcomes issues related to undetected temporal transitions from positive to negative associations common in IPS analysis. Further, in contrast to phase coherence, CRP unfolds the distribution of PS measures, which benefits subsequent clustering of PS matrices into recurring brain states.

Investigating the impact of autocorrelation on time-varying connectivity.

In recent years, a number of studies have reported on the existence of time-varying functional connectivity (TVC) in resting-state functional magnetic resonance imaging (rs-fMRI) data. The sliding-window technique is currently one of the most commonly used methods to estimate TVC. Although previous studies have shown that autocorrelation can negatively impact estimates of static functional connectivity, its impact on TVC estimates is not well known at this time. In this paper, we show both theoretically and empirically that the existence of autocorrelation within a time series can inflate the sampling variability of TVC estimated using the sliding-window technique. This can in turn increase the risk of misinterpreting noise as true TVC and negatively impact subsequent estimation of whole-brain time-varying FC profiles, or "brain states". The latter holds as more variable input measures lead to more variable output measures in the state estimation procedure. Finally, we demonstrate that prewhitening the data prior to analysis can lower the variance of the estimated TVC and improve brain state estimation. These results suggest that careful consideration is required when making inferences on TVC.

Tests of cortical parcellation based on white matter connectivity using diffusion tensor imaging.

The cerebral cortex is conventionally divided into a number of domains based on cytoarchitectural features. Diffusion tensor imaging (DTI) enables noninvasive parcellation of the cortex based on white matter connectivity patterns. However, the correspondence between DTI-connectivity-based and cytoarchitectural parcellation has not been systematically established. In this study, we compared histological parcellation of New World monkey neocortex to DTI- connectivity-based classification and clustering in the same brains. First, we used supervised classification to parcellate parieto-frontal cortex based on DTI tractograms and the cytoarchitectural prior (obtained using Nissl staining). We performed both within and across sample classification, showing reasonable classification performance in both conditions. Second, we used unsupervised clustering to parcellate the cortex and compared the clusters to the cytoarchitectonic standard. We then explored the similarities and differences with several post-hoc analyses, highlighting underlying principles that drive the DTI-connectivity-based parcellation. The differences in parcellation between DTI-connectivity and Nissl histology probably represent both DTI's bias toward easily-tracked bundles and true differences between cytoarchitectural and connectivity defined domains. DTI tractograms appear to cluster more according to functional networks, rather than mapping directly onto cytoarchitectonic domains. Our results show that caution should be used when DTI-tractography classification, based on data from another brain, is used as a surrogate for cytoarchitectural parcellation.

Improved estimation of subject-level functional connectivity using full and partial correlation with empirical Bayes shrinkage.

Reliability of subject-level resting-state functional connectivity (FC) is determined in part by the statistical techniques employed in its estimation. Methods that pool information across subjects to inform estimation of subject-level effects (e.g., Bayesian approaches) have been shown to enhance reliability of subject-level FC. However, fully Bayesian approaches are computationally demanding, while empirical Bayesian approaches typically rely on using repeated measures to estimate the variance components in the model. Here, we avoid the need for repeated measures by proposing a novel measurement error model for FC describing the different sources of variance and error, which we use to perform empirical Bayes shrinkage of subject-level FC towards the group average. In addition, since the traditional intra-class correlation coefficient (ICC) is inappropriate for biased estimates, we propose a new reliability measure denoted the mean squared error intra-class correlation coefficient (ICCMSE) to properly assess the reliability of the resulting (biased) estimates. We apply the proposed techniques to test-retest resting-state fMRI data on 461 subjects from the Human Connectome Project to estimate connectivity between 100 regions identified through independent components analysis (ICA). We consider both correlation and partial correlation as the measure of FC and assess the benefit of shrinkage for each measure, as well as the effects of scan duration. We find that shrinkage estimates of subject-level FC exhibit substantially greater reliability than traditional estimates across various scan durations, even for the most reliable connections and regardless of connectivity measure. Additionally, we find partial correlation reliability to be highly sensitive to the choice of penalty term, and to be generally worse than that of full correlations except for certain connections and a narrow range of penalty values. This suggests that the penalty needs to be chosen carefully when using partial correlations.

Comparing test-retest reliability of dynamic functional connectivity methods.

Due to the dynamic, condition-dependent nature of brain activity, interest in estimating rapid functional connectivity (FC) changes that occur during resting-state functional magnetic resonance imaging (rs-fMRI) has recently soared. However, studying dynamic FC is methodologically challenging, due to the low signal-to-noise ratio of the blood oxygen level dependent (BOLD) signal in fMRI and the massive number of data points generated during the analysis. Thus, it is important to establish methods and summary measures that maximize reliability and the utility of dynamic FC to provide insight into brain function. In this study, we investigated the reliability of dynamic FC summary measures derived using three commonly used estimation methods - sliding window (SW), tapered sliding window (TSW), and dynamic conditional correlations (DCC) methods. We applied each of these techniques to two publicly available rs-fMRI test-retest data sets - the Multi-Modal MRI Reproducibility Resource (Kirby Data) and the Human Connectome Project (HCP Data). The reliability of two categories of dynamic FC summary measures were assessed, specifically basic summary statistics of the dynamic correlations and summary measures derived from recurring whole-brain patterns of FC ("brain states"). The results provide evidence that dynamic correlations are reliably detected in both test-retest data sets, and the DCC method outperforms SW methods in terms of the reliability of summary statistics. However, across all estimation methods, reliability of the brain state-derived measures was low. Notably, the results also show that the DCC-derived dynamic correlation variances are significantly more reliable than those derived using the non-parametric estimation methods. This is important, as the fluctuations of dynamic FC (i.e., its variance) has a strong potential to provide summary measures that can be used to find meaningful individual differences in dynamic FC. We therefore conclude that utilizing the variance of the dynamic connectivity is an important component in any dynamic FC-derived summary measure.

Subject-specific regional measures of water diffusion are associated with impairment in chronic spinal cord injury.

PURPOSE: We aimed to identify non-invasive imaging parameters that can serve as biomarkers for the integrity of the spinal cord, which is paramount to neurological function. Diffusion tensor imaging (DTI) indices are sensitive to axonal and myelin damage, and have strong potential to serve as such biomarkers. However, averaging DTI indices over large regions of interest (ROIs), a common approach to analyzing the images of injured spinal cord, leads to loss of subject-specific information. We investigated if DTI-tractography-driven, subject-specific demarcation approach can yield measures that are more specific to impairment. METHODS: In 18 individuals with chronic spinal cord injury (SCI), subject-specific demarcation of the injury region was performed using DTI tractography, which yielded three regions relative to injury (RRI; regions superior to, at, and below injury epicenter). DTI indices averaged over each RRI were correlated with measures of residual motor and sensory function, obtained using the International Standard of Neurological Classification for Spinal Cord Injury (ISNCSCI). RESULTS: Total ISNCSCI score (ISNCSCI-tot; sum of ISNCSCI motor and sensory scores) was significantly (p < 0.05) correlated with fractional anisotropy and axial and radial diffusivities. ISNCSCI-tot showed strongest correlation with indices measured from the region inferior to the injury epicenter (IRRI), the degree of which exceeded that of those measured from the entire cervical cord-suggesting contribution from Wallerian degeneration. CONCLUSION: DTI tractography-driven, subject-specific injury demarcation approach provided measures that were more specific to impairment. Notably, DTI indices obtained from the IRRI region showed the highest specificity to impairment, demonstrating their strong potential as biomarkers for the SCI severity.

Presurgical brain mapping of the language network in patients with brain tumors using resting-state fMRI: Comparison with task fMRI.

PURPOSE: To compare language networks derived from resting-state fMRI (rs-fMRI) with task-fMRI in patients with brain tumors and investigate variables that affect rs-fMRI vs task-fMRI concordance. MATERIALS AND METHODS: Independent component analysis (ICA) of rs-fMRI was performed with 20, 30, 40, and 50 target components (ICA20 to ICA50) and language networks identified for patients presenting for presurgical fMRI mapping between 1/1/2009 and 7/1/2015. 49 patients were analyzed fulfilling criteria for presence of brain tumors, no prior brain surgery, and adequate task-fMRI performance. Rs-vs-task-fMRI concordance was measured using Dice coefficients across varying fMRI thresholds before and after noise removal. Multi-thresholded Dice coefficient volume under the surface (DiceVUS) and maximum Dice coefficient (MaxDice) were calculated. One-way Analysis of Variance (ANOVA) was performed to determine significance of DiceVUS and MaxDice between the four ICA order groups. Age, Sex, Handedness, Tumor Side, Tumor Size, WHO Grade, number of scrubbed volumes, image intensity root mean square (iRMS), and mean framewise displacement (FD) were used as predictors for VUS in a linear regression. RESULTS: Artificial elevation of rs-fMRI vs task-fMRI concordance is seen at low thresholds due to noise. Noise-removed group-mean DiceVUS and MaxDice improved as ICA order increased, however ANOVA demonstrated no statistically significant difference between the four groups. Linear regression demonstrated an association between iRMS and DiceVUS for ICA30-50, and iRMS and MaxDice for ICA50. CONCLUSION: Overall there is moderate group level rs-vs-task fMRI language network concordance, however substantial subject-level variability exists; iRMS may be used to determine reliability of rs-fMRI derived language networks.

New techniques for determining the longitudinal effects of local hemodynamics on the intima-media thickness in arteriovenous fistulae in an animal model.

Remodeling in the arteriovenous fistulas (AVFs) is believed to be a hemodynamic-driven process, which results in extreme changes in the diameter and intima-media thickening (IMT) of vessels over time. This study aims to describe the successful development of techniques that enabled correlation of changes in local and longitudinal wall shear stress (WSS) with the temporal variations of the diameter and IMT in the venous segment of AVFs. An AVF was created between the femoral artery and vein of a 50-kg pig. We have previously shown the successful use of CT-scan and ultrasound techniques for anatomical and flow measurements in AVFs, respectively. In this study, we developed new techniques involving markers (both in vivo and ex vivo), casting (ex vivo), and micro-MRI (ex vivo; 7 Tesla). A radiopaque marker (ROM) was sutured to the AVF at the day of surgery, which was visible in the CT-scan images, micro-MRI, and histology sections. Therefore, ROM served as a fixed local reference for both in vivo and ex vivo states of AVFs. Immediately after sacrificing the pig, a procedure was developed to create a cast from the AVF and thus, maintaining the in vivo state of the AVF during the histology process. Then, micro-MRI and histology techniques were conducted on the AVF to measure IMT in the vein. Along the ROM, the local changes in WSS levels for two cross-sections were tracked at 2D (D: days) and 28D post surgery. WSS levels reduced from 2D to 28D for both cross-sections. Also, the recirculation zones, which formed at 2D for both sections, became smaller in size at 28D. These hemodynamic changes were then mapped onto the corresponding IMT measurements from histology and micro-MRI. It was observed that the recirculation zones at 2D and 28D corresponded to the largest IMT in the two sections. In summary, the new methodologies allowed us to define a fixed local reference at all time points in the AVF, which enabled accurate tracking of local changes in hemodynamics (WSS), configuration (diameter), and structure (IMT) of the venous segment over time. This also empowered study of the interactions between these parameters, which could improve our understanding about the hemodynamic-driven remodeling in AVFs. From a clinical point of view, this information could be translated into local and early therapeutic interventions for dialysis patients.

Influence of temporal variation in wall shear stress on intima-media thickening in arteriovenous fistulae.

Arteriovenous fistula (AVF) failure is mainly due to venous stenosis characterized by significant amount of intima-media thickening (IMT), probably in the presence of negative (inward) remodeling. Our hypothesis is that the longitudinal changes in wall shear stress (WSS) within different configurations of AVF can influence remodeling factors (changes in luminal diameter (ΔDh ) and IMT) during its maturation process. Dh is an equivalent diameter for a noncircular conduit. A total of six AVFs with curved (C-AVF; n = 3) and straight (S-AVF; n = 3) configurations were created between the femoral artery and vein of three pigs, bilaterally. CT scans and ultrasounds were utilized to calculate local WSS at 2D (D: days), 7D, and 28D postsurgery. For each AVF, IMT was measured at four regions along the vein using morphometric analyses. At these regions, repeated measurements of WSS and luminal diameter of each AVF were obtained over time. The ΔD(h) between 7D and 28D was significantly larger for C-AVF than for S-AVF (2.27 ± 0.67 mm vs. 0.02 ± 0.55 mm; p < 0.05). Also, at 28D the amount of IMT in C-AVF (77.46 ± 7.10 units) was significantly greater (p < 0.05) when compared with S-AVF (53.71 ± 8.23 units). These structural changes were accompanied by significantly different gradients of WSS over time (τ') for C-AVF (-0.56 ± 0.60 dyne/cm(2)/day) in comparison with S-AVF (0.71 ± 0.39 dyne/cm(2)/day). Negative τ' for C-AVF corresponded to reduction in WSS level over time resulting in a physiological level of WSS at 28D (4.08 ± 5.08 dyne/cm(2)). In contrast, a positive τ' for S-AVF was associated with the increase in WSS levels over time causing high levels of WSS at 28D (36.68 ± 5.32 dyne/cm(2)). The decrease in WSS levels for the C-AVF over time was associated with outward remodeling of the venous wall (favorable to maturation). In contrast, for S-AVF, the increase in WSS levels over time was associated with inward remodeling and subsequently, venous stenosis. Thus, temporal gradients of WSS, which could be altered by the surgical configuration of AVF, may provide important information on the remodeling behavior of AVFs. Identification of an optimal AVF configuration, which results in a temporal decrease in WSS and an outward remodeling of the venous wall, may reduce AVF maturation failure.

Longitudinal assessment of hemodynamic endpoints in predicting arteriovenous fistula maturation.

Arteriovenous fistula (AVF) nonmaturation is currently a significant clinical problem; however, the mechanisms responsible for this have remained unanswered. Previous work by our group and others has suggested that anatomical configuration and the corresponding hemodynamic endpoints could have an important role in AVF remodeling. Thus, our goal was to assess the longitudinal (temporal) effect of wall shear stress (WSS) on remodeling process of AVFs with two different configurations. The hypothesis is that early assessment of hemodynamic endpoints such as temporal gradient of WSS will predict the maturation status of AVF at later time points. Two AVFs with curved (C-AVF) and straight (S-AVF) configurations were created between the femoral artery and vein of each pig. Three pigs were considered in this study and in total six AVFs (three C-AVF and three S-AVF) were created. The CT scan and ultrasound were utilized to numerically evaluate local WSS at 20 cross-sections along the venous segment of AVFs at 2D (D: days), 7D, and 28D postsurgery. These cross-sections were located at 1.5 mm increments from the anastomosis junction. Local WSS values at these cross-sections were correlated with their corresponding luminal area over time. The WSS in C-AVF decreased from 22.3 ± 4.8 dyn/cm(2) at 2D to 4.1 ± 5.1 dyn/cm(2) at 28D, while WSS increased in S-AVF from 13.0 ± 5.0 dyn/cm(2) at 2D to 36.7 ± 5.3 dyn/cm(2) at 28D. Corresponding to these changes in WSS levels, luminal area of C-AVF dilated (0.23 ± 0.14 cm(2) at 2D to 0.87 ± 0.14 cm(2) at 28D) with attendant increase in flow rate. However, S-AVF had minimal changes in area (0.26 ± 0.02 cm(2) at 2D to 0.27 ± 0.03 cm(2) at 28D) despite some increase in flow rate. Our results suggest that the temporal changes of WSS could have significant effects on AVF maturation. Reduction in WSS over time (regardless of initial values) may result in dilation (p < 0.05), while increase in WSS may be detrimental to maturation. Thus, creation of AVFs in a specific configuration which results in a decline in WSS over time may reduce AVF maturation failure.

Anatomic configuration affects the flow rate and diameter of porcine arteriovenous fistulae.

Although arteriovenous fistulae are currently the preferred form of vascular access, early failure is a significant problem. Since wall shear stress is thought to play an important role in the pathogenesis of early failure, and this stress varies markedly in different fistula configurations, we assessed the influence of configuration (curved or straight) on longitudinal changes of flow rate and lumen diameter in a porcine fistula model. Fistulae were created in eight pigs between the femoral artery and vein, with each animal having a curved and a straight configuration on opposite sides. Velocity measurements were obtained by ultrasound at the time of surgery and at intermediate time points up to 28 days. Quantification of both the configuration and the internal diameter of the fistulae was determined by CT scans. The overall rate of increased flow during each time interval (0 to 2 days, 2 to 7 days, and 7 to 28 days) was more pronounced with the curved fistulae. Moreover, the luminal diameter of curved fistulae had dilated more from the time of surgery to 28 days as compared to the straight fistulae. Thus, anatomical configuration of fistulae plays a major role in flow-mediated dilatation. Identifying the optimal configuration may result in increased diameter and consequently blood flow, and perhaps reduce the incidence of early failure.

Phase-locking of resting-state brain networks with the gastric basal electrical rhythm.

A network of myenteric interstitial cells of Cajal in the corpus of the stomach serves as its "pacemaker", continuously generating a ca 0.05 Hz electrical slow wave, which is transmitted to the brain chiefly by vagal afferents. A recent study combining resting-state functional MRI (rsfMRI) with concurrent surface electrogastrography (EGG), with cutaneous electrodes placed on the epigastrium, found 12 brain regions with activity that was significantly phase-locked with this gastric basal electrical rhythm. Therefore, we asked whether fluctuations in brain resting state networks (RSNs), estimated using a spatial independent component analysis (ICA) approach, might be synchronized with the stomach. In the present study, in order to determine whether any RSNs are phase-locked with the gastric rhythm, an individual participant underwent 22 scanning sessions; in each, two 15-minute runs of concurrent EGG and rsfMRI data were acquired. EGG data from three sessions had weak gastric signals and were excluded; the other 19 sessions yielded a total of 9.5 hours of data. The rsfMRI data were analyzed using group ICA; RSN time courses were estimated; for each run, the phase-locking value (PLV) was computed between each RSN and the gastric signal. To assess statistical significance, PLVs from all pairs of "mismatched" data (EGG and rsfMRI data acquired on different days) were used as surrogate data to generate a null distribution for each RSN. Of a total of 18 RSNs, three were found to be significantly phase-locked with the basal gastric rhythm, namely, a cerebellar network, a dorsal somatosensory-motor network, and a default mode network. Disruptions to the gut-brain axis, which sustains interoceptive feedback between the central nervous system and the viscera, are thought to be involved in various disorders; manifestation of the infra-slow rhythm of the stomach in brain rsfMRI data could be useful for studies in clinical populations.

Altered functional connectivity between sub-regions in the thalamus and cortex in schizophrenia patients measured by resting state BOLD fMRI at 7T.

The thalamus is a small brain structure that relays neuronal signals between subcortical and cortical regions. Abnormal thalamocortical connectivity in schizophrenia has been reported in previous studies using blood-oxygenation-level-dependent (BOLD) functional MRI (fMRI) performed at 3T. However, anatomically the thalamus is not a single entity, but is subdivided into multiple distinct nuclei with different connections to various cortical regions. We sought to determine the potential benefit of using the enhanced sensitivity of BOLD fMRI at ultra-high magnetic field (7T) in exploring thalamo-cortical connectivity in schizophrenia based on subregions in the thalamus. Seeds placed in thalamic subregions of 14 patients and 14 matched controls were used to calculate whole-brain functional connectivity. Our results demonstrate impaired thalamic connectivity to the prefrontal cortex and the cerebellum, but enhanced thalamic connectivity to the motor/sensory cortex in schizophrenia. This altered functional connectivity significantly correlated with disease duration in the patients. Remarkably, comparable effect sizes observed in previous 3T studies were detected in the current 7T study with a heterogeneous and much smaller cohort, providing evidence that ultra-high field fMRI may be a powerful tool for measuring functional connectivity abnormalities in schizophrenia. Further investigation with a larger cohort is merited to validate the current findings.

A magnetic resonance imaging-based articulatory and acoustic study of "retroflex" and "bunched" American English /r/.

Speakers of rhotic dialects of North American English show a range of different tongue configurations for /r/. These variants produce acoustic profiles that are indistinguishable for the first three formants [Delattre, P., and Freeman, D. C., (1968). "A dialect study of American English r's by x-ray motion picture," Linguistics 44, 28-69; Westbury, J. R. et al. (1998), "Differences among speakers in lingual articulation for American English /r/," Speech Commun. 26, 203-206]. It is puzzling why this should be so, given the very different vocal tract configurations involved. In this paper, two subjects whose productions of "retroflex" /r/ and "bunched" /r/ show similar patterns of F1-F3 but very different spacing between F4 and F5 are contrasted. Using finite element analysis and area functions based on magnetic resonance images of the vocal tract for sustained productions, the results of computer vocal tract models are compared to actual speech recordings. In particular, formant-cavity affiliations are explored using formant sensitivity functions and vocal tract simple-tube models. The difference in F4/F5 patterns between the subjects is confirmed for several additional subjects with retroflex and bunched vocal tract configurations. The results suggest that the F4/F5 differences between the variants can be largely explained by differences in whether the long cavity behind the palatal constriction acts as a half- or a quarter-wavelength resonator.

An image-processing toolset for diffusion tensor tractography.

Diffusion tensor imaging (DTI)-based fiber tractography holds great promise in delineating neuronal fiber tracts and, hence, providing connectivity maps of the neural networks in the human brain. An array of image-processing techniques has to be developed to turn DTI tractography into a practically useful tool. To this end, we have developed a suite of image-processing tools for fiber tractography with improved reliability. This article summarizes the main technical developments we have made to date, which include anisotropic smoothing, anisotropic interpolation, Bayesian fiber tracking and automatic fiber bundling. A primary focus of these techniques is the robustness to noise and partial volume averaging, the two major hurdles to reliable fiber tractography. Performance of these techniques has been comprehensively examined with simulated and in vivo DTI data, demonstrating improvements in the robustness and reliability of DTI tractography.

Advances in Spinal Functional Magnetic Resonance Imaging in the Healthy and Injured Spinal Cords.

PURPOSE OF REVIEW: This review provides an overview of the current spinal functional magnetic resonance imaging (fMRI) studies that investigate the healthy and injured spinal cords. RECENT FINDINGS: Spinal fMRI-derived outcome measures have previously been suggested to be sensitive to changes in neurological function in the spinal cord. A body of recent task-activated fMRI studies seems to confirm that detecting neural activity in the spinal cord using spinal fMRI may be feasible as well as reliable. Furthermore, a growing number of studies has shown that resting state fMRI in the spinal cord is also feasible, demonstrating that the investigation of changes in neural activity can also be performed in the absence of explicit tasks. SUMMARY: Current task-activated and resting state fMRI studies suggest that spinal fMRI has a strong potential to provide novel imaging biomarkers that can be used to investigate plastic changes in the injured spinal cord.

Reproducibility and variation of diffusion measures in the squirrel monkey brain, in vivo and ex vivo.

PURPOSE: Animal models are needed to better understand the relationship between diffusion MRI (dMRI) and the underlying tissue microstructure. One promising model for validation studies is the common squirrel monkey, Saimiri sciureus. This study aims to determine (1) the reproducibility of in vivo diffusion measures both within and between subjects; (2) the agreement between in vivo and ex vivo data acquired from the same specimen and (3) normal diffusion values and their variation across brain regions. METHODS: Data were acquired from three healthy squirrel monkeys, each imaged twice in vivo and once ex vivo. Reproducibility of fractional anisotropy (FA), mean diffusivity (MD), and principal eigenvector (PEV) was assessed, and normal values were determined both in vivo and ex vivo. RESULTS: The calculated coefficients of variation (CVs) for both intra-subject and inter-subject MD were below 10% (low variability) while FA had a wider range of CVs, 2-14% intra-subject (moderate variability), and 3-31% inter-subject (high variability). MD in ex vivo tissue was lower than in vivo (30%-50% decrease), while FA values increased in all regions (30-39% increase). The mode of angular differences between in vivo and ex vivo PEVs was 12 degrees. CONCLUSION: This study characterizes the diffusion properties of the squirrel monkey brain and serves as the groundwork for using the squirrel monkey, both in vivo and ex vivo, as a model for diffusion MRI studies.

Dynamic Functional Connectivity States Between the Dorsal and Ventral Sensorimotor Networks Revealed by Dynamic Conditional Correlation Analysis of Resting-State Functional Magnetic Resonance Imaging.

Functional connectivity in resting-state functional magnetic resonance imaging (rs-fMRI) has received substantial attention since the initial findings of Biswal et al. Traditional network correlation metrics assume that the functional connectivity in the brain remains stationary over time. However, recent studies have shown that robust temporal fluctuations of functional connectivity among as well as within functional networks exist, challenging this assumption. In this study, these dynamic correlation differences were investigated between the dorsal and ventral sensorimotor networks by applying the dynamic conditional correlation model to rs-fMRI data of 20 healthy subjects. k-Means clustering was used to determine an optimal number of discrete connectivity states (k = 10) of the sensorimotor system across all subjects. Our analysis confirms the existence of differences in dynamic correlation between the dorsal and ventral networks, with highest connectivity found within the ventral motor network.

A 3D high resolution ex vivo white matter atlas of the common squirrel monkey (Saimiri sciureus) based on diffusion tensor imaging.

Modern magnetic resonance imaging (MRI) brain atlases are high quality 3-D volumes with specific structures labeled in the volume. Atlases are essential in providing a common space for interpretation of results across studies, for anatomical education, and providing quantitative image-based navigation. Extensive work has been devoted to atlas construction for humans, macaque, and several non-primate species (e.g., rat). One notable gap in the literature is the common squirrel monkey - for which the primary published atlases date from the 1960's. The common squirrel monkey has been used extensively as surrogate for humans in biomedical studies, given its anatomical neuro-system similarities and practical considerations. This work describes the continued development of a multi-modal MRI atlas for the common squirrel monkey, for which a structural imaging space and gray matter parcels have been previously constructed. This study adds white matter tracts to the atlas. The new atlas includes 49 white matter (WM) tracts, defined using diffusion tensor imaging (DTI) in three animals and combines these data to define the anatomical locations of these tracks in a standardized coordinate system compatible with previous development. An anatomist reviewed the resulting tracts and the inter-animal reproducibility (i.e., the Dice index of each WM parcel across animals in common space) was assessed. The Dice indices range from 0.05 to 0.80 due to differences of local registration quality and the variation of WM tract position across individuals. However, the combined WM labels from the 3 animals represent the general locations of WM parcels, adding basic connectivity information to the atlas.

Intrinsic Visual-Motor Synchrony Correlates With Social Deficits in Autism.

BACKGROUND: Imitation, which is impaired in children with autism spectrum disorder (ASD) and critically depends on the integration of visual input with motor output, likely impacts both motor and social skill acquisition in children with ASD; however, it is unclear what brain mechanisms contribute to this impairment. Children with ASD also exhibit what appears to be an ASD-specific bias against using visual feedback during motor learning. Does the temporal congruity of intrinsic activity, or functional connectivity, between motor and visual brain regions contribute to ASD-associated deficits in imitation, motor, and social skills? METHODS: We acquired resting-state functional magnetic resonance imaging scans from 100 8- to 12-year-old children (50 ASD). Group independent component analysis was used to estimate functional connectivity between visual and motor systems. Brain-behavior relationships were assessed by regressing functional connectivity measures with social deficit severity, imitation, and gesture performance scores. RESULTS: We observed increased intrinsic asynchrony between visual and motor systems in children with ASD and replicated this finding in an independent sample from the Autism Brain Imaging Data Exchange. Moreover, children with more out-of-sync intrinsic visual-motor activity displayed more severe autistic traits, while children with greater intrinsic visual-motor synchrony were better imitators. CONCLUSIONS: Our twice replicated findings confirm that visual-motor functional connectivity is disrupted in ASD. Furthermore, the observed temporal incongruity between visual and motor systems, which may reflect diminished integration of visual consequences with motor output, was predictive of the severity of social deficits and may contribute to impaired social-communicative skill development in children with ASD.