Gender differences in the structural connectome of the teenage brain revealed by generalized q-sampling MRI.

The question of whether there are biological differences between male and female brains is a fraught one, and political positions and prior expectations seem to have a strong influence on the interpretation of scientific data in this field. This question is relevant to issues of gender differences in the prevalence of psychiatric conditions, including autism, attention deficit hyperactivity disorder (ADHD), Tourette's syndrome, schizophrenia, dyslexia, depression, and eating disorders. Understanding how gender influences vulnerability to these conditions is significant. Diffusion magnetic resonance imaging (dMRI) provides a non-invasive method to investigate brain microstructure and the integrity of anatomical connectivity. Generalized q-sampling imaging (GQI) has been proposed to characterize complicated fiber patterns and distinguish fiber orientations, providing an opportunity for more accurate, higher-order descriptions through the water diffusion process. Therefore, we aimed to investigate differences in the brain's structural network between teenage males and females using GQI. This study included 59 (i.e., 33 males and 26 females) age- and education-matched subjects (age range: 13 to 14 years). The structural connectome was obtained by graph theoretical and network-based statistical (NBS) analyses. Our findings show that teenage male brains exhibit better intrahemispheric communication, and teenage female brains exhibit better interhemispheric communication. Our results also suggest that the network organization of teenage male brains is more local, more segregated, and more similar to small-world networks than teenage female brains. We conclude that the use of an MRI study with a GQI-based structural connectomic approach like ours presents novel insights into network-based systems of the brain and provides a new piece of the puzzle regarding gender differences.

Theoretical Bounds of Direct Binary Search Halftoning.

Direct binary search (DBS) produces the images of the best quality among half-toning algorithms. The reason is that it minimizes the total squared perceived error instead of using heuristic approaches. The search for the optimal solution involves two operations: (1) toggle and (2) swap. Both operations try to find the binary states for each pixel to minimize the total squared perceived error. This error energy minimization leads to a conjecture that the absolute value of the filtered error after DBS converges is bounded by half of the peak value of the autocorrelation filter. However, a proof of the bound's existence has not yet been found. In this paper, we present a proof that shows the bound existed as conjectured under the condition that at least one swap occurs after toggle converges. The theoretical analysis also indicates that a swap with a pixel further away from the center of the autocorrelation filter results in a tighter bound. Therefore, we propose a new DBS algorithm which considers toggle and swap separately, and the swap operations are considered in the order from the edge to the center of the filter. Experimental results show that the new algorithm is more efficient than the previous algorithm and can produce half-toned images of the same quality as the previous algorithm.

VBM Reveals Brain Volume Differences between Parkinson's Disease and Essential Tremor Patients.

Symptoms of essential tremor (ET) are similar to those of Parkinson's disease (PD) during their initial stages. Presently, there are few stable biomarkers available on a neuroanatomical level for distinguishing between these two diseases. However, few investigations have directly compared the changes in brain volume and assessed the compensatory effects of a change in the parts of the brain associated with PD and with ET. To determine the compensatory and/or degenerative anatomical changes in the brains of PD and ET patients, the present study tested, via two voxel-based morphometry (VBM) approaches (Basic vs. DARTEL VBM processing), the anatomical brain images of 10 PD and 10 ET patients, as well as of 13 age-matched normal controls, obtained through a 3T magnetic resonance scanner. These findings indicate that PD and ET caused specific patterns of brain volume alterations in the brains examined. In addition, our observations also revealed compensatory effects, or self-reorganization, occurring in the thalamus and the middle temporal gyrus in the PD and ET patients, due perhaps in part to the enhanced thalamocortical sensorimotor interaction and the head-eye position readjustment, respectively, in these PD and ET patients. Such a distinction may lend itself to use as a biomarker for differentiating between these two diseases.

Accelerating fourier volume rendering by polar coordinate data representation.

Volume rendering is an important tool to visualize three-dimensional data in biomedicine by projecting the data to a two-dimensional plane. The projection is done by ray casting and its complexity is proportional to the number of three-dimensional data points. To reduce complexity, Fourier volume rendering (FVR) uses slice projection theorem to facilitate the integration of voxels along the ray casting path. In this paper, we proposed a new method for FVR that stored and processed the frequency domain data in polar coordinate. By exploiting three aspects of data processing which is previously impossible in rectilinear coordinate, our new method is much faster than the previous methods. The first aspect is data regularity. When data are stored in polar coordinate, extracting a slice involves accessing data stored in adjacent memory location. This regularity makes memory access more efficient. The second aspect is to utilize the high data density near the origin in polar coordinate. We can obtain two benefits from this aspect. The first allows us to extract a slice by nearest-neighbor interpolation instead of more complex interpolation but without sacrificing image quality. The second allows us to trade off between image quality and memory storage. The third aspect is to recognize that converting from rectilinear coordinate to polar coordinate is a one-time process. Therefore, we can use a better interpolation kernel with larger support in coordinate conversion. In turn, most of the computation is shifted to the preprocessing stage and interactive rendering can be made very fast. In the experiments, we show that the speed in interactive visualization for our new method is independent of the size of the interpolation kernel, therefore, achieving comparable image quality at a faster rate than previous methods.

Three-dimensional elliptical reconstruction for stereoscopic magnetic resonance angiography.

Stereoscopic MRA acquires a pair of blood vessel projections at two different viewing angles. Previously, we have developed two algorithms to reconstruct 3-D blood vessels from stereoscopic MRA. The assumption we made was that blood vessels were tilting circular tubes and the shape of the vessel on every cross-section was an ellipse. Since an ellipse can be represented in either algebraic form or parametric form, our previous algorithms reconstructed the ellipses by representing them in these two forms. In this paper, we further improved the accuracy of our previous algorithms by an order through two enhancements. The first improvement we made was a better method to estimate the rotation angle of the major axis of an ellipse. Instead of using the center of two adjacent ellipses to estimate the rotation angle as in our previous method, the new method used the projection lengths of the two views to estimate the angle. The second improvement we made was the equation to describe the relationship between the major/minor axes and the projection lengths. In our experiments, the average estimation error for the parametric algorithm was improved from 0.471 pixels to 0.066 pixels. The average error for the algebraic algorithm was improved from 0.101 pixels to 0.014 pixels.

Temporally shifted hemodynamic response model helps to extract acupuncture-induced functional magnetic resonance imaging blood oxygenation-level dependent activities.

BACKGROUND: The onsets of needling sensation introduced by acupuncture stimulus can vary widely from subject to subject. This should be explicitly accounted for by the model blood oxygenation-level dependent (BOLD) time course used in general linear model (GLM) analysis to obtain more consistent across-subject group results. However, in standard GLM analysis, the model BOLD time course obtained by convolving a canonical hemodynamic response function with an experimental paradigm time course is assumed identical across subjects. Although some added-on properties to the model BOLD time course, such as temporal and dispersion derivatives, may be used to account for different BOLD response onsets, they can only account for the BOLD onset deviations to the extent of less than one repetition time (TR). METHODS: In this study, we explicitly manipulated the onsets of model BOLD time course by shifting it with -2, -1, or 1 TR and used these temporally shifted BOLD model to analyze the functional magnetic resonance imaging (fMRI) data obtained from three acupuncture fMRI experiments with GLM analysis. One involved acupuncture stimulus on left ST42 acupoint and the other two on left GB40 and left BL64 acupoints. RESULTS: The model BOLD time course with temporal shifts, in addition to temporal and dispersion derivatives, could result in better statistical power of the data analysis in terms of the average correlation coefficients between the used BOLD models and extracted BOLD responses from individual subject data and the T-values of the activation clusters in the grouped random effects. CONCLUSIONS: The GLM analysis with ordinary BOLD model failed to catch the large variability of the onsets of the BOLD responses associated with the acupuncture needling sensation. Shifts in time with more than a TR on model BOLD time course might be required to better extract the acupuncture stimulus-induced BOLD activities from individual fMRI data.

Carryover effects alter FMRI statistical analysis in an acupuncture study.

Carryover effects can contaminate ON/OFF BOLD contrasts designated in an fMRI experiment. Yet, the ON/OFF contrasts are essential to facilitate statistical analysis based on the significance of contrast levels. Here, we conducted an fMRI experiment with acupuncture stimulation applied on ST42 acupoint as well as with tactile stimulation on its skin surface. Experiment consisted of three two-block acupuncture and one two-block tactile fMRI runs. Each block started with 26-sec OFF period followed by either 26-sec needle manipulation in the acupuncture runs or by scratching skin surface with sand paper in the tactile. To test if carryover effects could alter the BOLD contrasts, we analyzed different portions of fMRI data using GLM method. Our results showed analyses on different portions of acupuncture fMRI data gave significantly different results. Statistical parametric maps of group random effects resulted from the analysis on the very first fMRI trial formed the broadest coverage of the active brain areas. BOLD model time course also best explained the adjusted raw time course at peak active voxel (coefficient of determination = 0.88). Analyses on other portions of fMRI data only selected subset of the active brain areas delineated by the analysis on the very first data trial and the BOLD model only mildly accounted for the adjusted raw time courses. In tactile runs, results were more consistent across analyses. Therefore, in fMRI experiments with strong carryover effects, a single-block experimental design with multiple repetitions, separated by long enough periods of time, should be more suitable to extract task BOLD effects.

Neuronal specificity of acupuncture response: a fMRI study with electroacupuncture.

Recently, neuronal correlates of acupuncture stimulation in human brain have been investigated by functional neuroimaging. The preliminary findings suggest that acupuncture at analgesic points involves the pain-related neuromatrix and may have acupoint-brain correlation. Although multiple models of control stimulations have been applied to address the specificity of the needling effect clinically, their impacts have not been evaluated by functional neuroimaging. With the advantage of objective parameter setting, electroacupuncture (EA) was used in this study to devise three distinct controls for real EA, i.e., mock EA (no stimulation), minimal EA (superficial and light stimulation), and sham EA (same stimulation as real EA) applied at nonmeridian points. Fifteen healthy volunteers received real EA at analgesic point Gallbladder 34 (Yanglinquan), sham EA, and one of either mock EA or minimal EA over the left leg in counter-balanced orders. Multisubject analysis showed that sham EA and real EA both activated the reported distributed pain neuromatrix. However, real EA elicited significantly higher activation than sham EA over the hypothalamus and primary somatosensory-motor cortex and deactivation over the rostral segment of anterior cingulate cortex. In the comparison of minimal EA versus mock EA, minimal EA elicited significantly higher activation over the medial occipital cortex. Single-subject analysis showed that superior temporal gyrus (encompassing the auditory cortex) and medial occipital cortex (encompassing the visual cortex) frequently respond to minimal EA, sham EA, or real EA. We concluded that the hypothalamus-limbic system was significantly modulated by EA at acupoints rather than at nonmeridian points, while visual and auditory cortical activation was not a specific effect of treatment-relevant acupoints and required further investigation of the underlying neurophysiological mechanisms.

Spoiled gradient-echo as an arterial spin tagging technique for quick evaluation of local perfusion.

PURPOSE: To introduce a simple gradient-echo arterial spin tagging (GREAST) technique available for most commercial magnetic resonance (MR) systems, for a quick evaluation of tissue perfusion. MATERIALS AND METHODS: The GREAST technique uses a combination of a short TR spoiled gradient-echo (SPGR) sequence with a selective presaturation radio frequency (RF) pulse that allows acquiring each tagged and control image within 10-20 seconds. The phantom and human studies were performed for evaluating the feasibility in measurement of local perfusion and the efficacy in alleviation of the asymmetric magnetization transfer (MT) and slice profile effects. RESULTS: Results show a good linear relationship between the signal attenuation caused by the presaturation pulse and flow rates in the phantom experiment and effective alleviation of the asymmetric MT and slice profile effects for various orientations of imaging slices. Human studies showed good perfusion results in brain imaging. Perfusion imaging on the liver and kidney were also conducted. The results could be significantly improved by effectively lessening motion-related artifacts. CONCLUSION: The GREAST technique is simple, easy to use, and applicable to examine local perfusion of the brain and other organs in the trunk.

Platform-independent image reconstruction for spiral magnetic resonance imaging.

The distinct feature of data acquisition for magnetic resonance imaging (MRI) is that the data are sampled in frequency domain instead of spatial domain. Therefore, the acquired data must be inverse Fourier transformed to generate images. To apply fast Fourier transform (FFT), the data are usually acquired on rectilinear grids. However, acquiring data on rectilinear grids is not very efficient in MRI. A spiral trajectory, which starts at the origin of the frequency domain and spins out to higher spatial frequency is more efficient and faster than the conventional method. Since the spiral trajectories do not sample on rectilinear grids, raw data must be re-interpolated onto rectilinear grids prior to inverse FFT. This re-gridding process is done using a reconstruction program. When the platforms to run the program grow, the efforts required to maintain the program become prohibitive. This problem can be solved through the platform-independent Java programming language. In this paper, we report on our attempt to implement the spiral MRI reconstruction program in Java. We show that the performance is not significantly affected and that it is practical to use a platform-independent reconstruction software.