Meta-analysis in human neuroimaging: computational modeling of large-scale databases.

Spatial normalization--applying standardized coordinates as anatomical addresses within a reference space--was introduced to human neuroimaging research nearly 30 years ago. Over these three decades, an impressive series of methodological advances have adopted, extended, and popularized this standard. Collectively, this work has generated a methodologically coherent literature of unprecedented rigor, size, and scope. Large-scale online databases have compiled these observations and their associated meta-data, stimulating the development of meta-analytic methods to exploit this expanding corpus. Coordinate-based meta-analytic methods have emerged and evolved in rigor and utility. Early methods computed cross-study consensus, in a manner roughly comparable to traditional (nonimaging) meta-analysis. Recent advances now compute coactivation-based connectivity, connectivity-based functional parcellation, and complex network models powered from data sets representing tens of thousands of subjects. Meta-analyses of human neuroimaging data in large-scale databases now stand at the forefront of computational neurobiology.

BrainMap VBM: An environment for structural meta-analysis.

The BrainMap database is a community resource that curates peer-reviewed, coordinate-based human neuroimaging literature. By pairing the results of neuroimaging studies with their relevant meta-data, BrainMap facilitates coordinate-based meta-analysis (CBMA) of the neuroimaging literature en masse or at the level of experimental paradigm, clinical disease, or anatomic location. Initially dedicated to the functional, task-activation literature, BrainMap is now expanding to include voxel-based morphometry (VBM) studies in a separate sector, titled: BrainMap VBM. VBM is a whole-brain, voxel-wise method that measures significant structural differences between or within groups which are reported as standardized, peak x-y-z coordinates. Here we describe BrainMap VBM, including the meta-data structure, current data volume, and automated reverse inference functions (region-to-disease profile) of this new community resource. CBMA offers a robust methodology for retaining true-positive and excluding false-positive findings across studies in the VBM literature. As with BrainMap's functional database, BrainMap VBM may be synthesized en masse or at the level of clinical disease or anatomic location. As a use-case scenario for BrainMap VBM, we illustrate a trans-diagnostic data-mining procedure wherein we explore the underlying network structure of 2,002 experiments representing over 53,000 subjects through independent components analysis (ICA). To reduce data-redundancy effects inherent to any database, we demonstrate two data-filtering approaches that proved helpful to ICA. Finally, we apply hierarchical clustering analysis (HCA) to measure network- and disease-specificity. This procedure distinguished psychiatric from neurological diseases. We invite the neuroscientific community to further exploit BrainMap VBM with other modeling approaches.

Implementation errors in the GingerALE Software: Description and recommendations.

Neuroscience imaging is a burgeoning, highly sophisticated field the growth of which has been fostered by grant-funded, freely distributed software libraries that perform voxel-wise analyses in anatomically standardized three-dimensional space on multi-subject, whole-brain, primary datasets. Despite the ongoing advances made using these non-commercial computational tools, the replicability of individual studies is an acknowledged limitation. Coordinate-based meta-analysis offers a practical solution to this limitation and, consequently, plays an important role in filtering and consolidating the enormous corpus of functional and structural neuroimaging results reported in the peer-reviewed literature. In both primary data and meta-analytic neuroimaging analyses, correction for multiple comparisons is a complex but critical step for ensuring statistical rigor. Reports of errors in multiple-comparison corrections in primary-data analyses have recently appeared. Here, we report two such errors in GingerALE, a widely used, US National Institutes of Health (NIH)-funded, freely distributed software package for coordinate-based meta-analysis. These errors have given rise to published reports with more liberal statistical inferences than were specified by the authors. The intent of this technical report is threefold. First, we inform authors who used GingerALE of these errors so that they can take appropriate actions including re-analyses and corrective publications. Second, we seek to exemplify and promote an open approach to error management. Third, we discuss the implications of these and similar errors in a scientific environment dependent on third-party software. Hum Brain Mapp 38:7-11, 2017. © 2016 Wiley Periodicals, Inc.

Multimodal Randomized Functional MR Imaging of the Effects of Methylene Blue in the Human Brain.

Purpose To investigate the sustained-attention and memory-enhancing neural correlates of the oral administration of methylene blue in the healthy human brain. Materials and Methods The institutional review board approved this prospective, HIPAA-compliant, randomized, double-blinded, placebo-controlled clinical trial, and all patients provided informed consent. Twenty-six subjects (age range, 22-62 years) were enrolled. Functional magnetic resonance (MR) imaging was performed with a psychomotor vigilance task (sustained attention) and delayed match-to-sample tasks (short-term memory) before and 1 hour after administration of low-dose methylene blue or a placebo. Cerebrovascular reactivity effects were also measured with the carbon dioxide challenge, in which a 2 × 2 repeated-measures analysis of variance was performed with a drug (methylene blue vs placebo) and time (before vs after administration of the drug) as factors to assess drug × time between group interactions. Multiple comparison correction was applied, with cluster-corrected P < .05 indicating a significant difference. Results Administration of methylene blue increased response in the bilateral insular cortex during a psychomotor vigilance task (Z = 2.9-3.4, P = .01-.008) and functional MR imaging response during a short-term memory task involving the prefrontal, parietal, and occipital cortex (Z = 2.9-4.2, P = .03-.0003). Methylene blue was also associated with a 7% increase in correct responses during memory retrieval (P = .01). Conclusion Low-dose methylene blue can increase functional MR imaging activity during sustained attention and short-term memory tasks and enhance memory retrieval. © RSNA, 2016 Online supplemental material is available for this article.

Intrinsic Resting-State Functional Connectivity in the Human Spinal Cord at 3.0 T.

PURPOSE: To apply resting-state functional magnetic resonance (MR) imaging to map functional connectivity of the human spinal cord. MATERIALS AND METHODS: Studies were performed in nine self-declared healthy volunteers with informed consent and institutional review board approval. Resting-state functional MR imaging was performed to map functional connectivity of the human cervical spinal cord from C1 to C4 at 1 × 1 × 3-mm resolution with a 3.0-T clinical MR imaging unit. Independent component analysis (ICA) was performed to derive resting-state functional MR imaging z-score maps rendered on two-dimensional and three-dimensional images. Seed-based analysis was performed for cross validation with ICA networks by using Pearson correlation. RESULTS: Reproducibility analysis of resting-state functional MR imaging maps from four repeated trials in a single participant yielded a mean z score of 6 ± 1 (P < .0001). The centroid coordinates across the four trials deviated by 2 in-plane voxels ± 2 mm (standard deviation) and up to one adjacent image section ± 3 mm. ICA of group resting-state functional MR imaging data revealed prominent functional connectivity patterns within the spinal cord gray matter. There were statistically significant (z score > 3, P < .001) bilateral, unilateral, and intersegmental correlations in the ventral horns, dorsal horns, and central spinal cord gray matter. Three-dimensional surface rendering provided visualization of these components along the length of the spinal cord. Seed-based analysis showed that many ICA components exhibited strong and significant (P < .05) correlations, corroborating the ICA results. Resting-state functional MR imaging connectivity networks are qualitatively consistent with known neuroanatomic and functional structures in the spinal cord. CONCLUSION: Resting-state functional MR imaging of the human cervical spinal cord with a 3.0-T clinical MR imaging unit and standard MR imaging protocols and hardware reveals prominent functional connectivity patterns within the spinal cord gray matter, consistent with known functional and anatomic layouts of the spinal cord.

Sulcal depth-position profile is a genetically mediated neuroscientific trait: description and characterization in the central sulcus.

Genetic and environmental influences on brain morphology were assessed in an extended-pedigree design by extracting depth-position profiles (DPP) of the central sulcus (CS). T1-weighted magnetic resonance images were used to measure CS length and depth in 467 human subjects from 35 extended families. Three primary forms of DPPs were observed. The most prevalent form, present in 70% of subjects, was bimodal, with peaks near hand and mouth regions. Trimodal and unimodal configurations accounted for 15 and 8%, respectively. Genetic control accounted for 56 and 66% of between-subject variance in average CS depth and length, respectively, and was not significantly influenced by environmental factors. Genetic control over CS depth ranged from 1 to 50% across the DPP. Areas of peak heritability occurred at locations corresponding to hand and mouth areas. Left and right analogous CS depth measurements were strongly pleiotropic. Shared genetic influence lessened as the distance between depth measurements was increased. We argue that DPPs are powerful phenotypes that should inform genetic influence of more complex brain regions and contribute to gene discovery efforts.

Concurrent TMS to the primary motor cortex augments slow motor learning.

Transcranial magnetic stimulation (TMS) has shown promise as a treatment tool, with one FDA approved use. While TMS alone is able to up- (or down-) regulate a targeted neural system, we argue that TMS applied as an adjuvant is more effective for repetitive physical, behavioral and cognitive therapies, that is, therapies which are designed to alter the network properties of neural systems through Hebbian learning. We tested this hypothesis in the context of a slow motor learning paradigm. Healthy right-handed individuals were assigned to receive 5 Hz TMS (TMS group) or sham TMS (sham group) to the right primary motor cortex (M1) as they performed daily motor practice of a digit sequence task with their non-dominant hand for 4 weeks. Resting cerebral blood flow (CBF) was measured by H2(15)O PET at baseline and after 4 weeks of practice. Sequence performance was measured daily as the number of correct sequences performed, and modeled using a hyperbolic function. Sequence performance increased significantly at 4 weeks relative to baseline in both groups. The TMS group had a significant additional improvement in performance, specifically, in the rate of skill acquisition. In both groups, an improvement in sequence timing and transfer of skills to non-trained motor domains was also found. Compared to the sham group, the TMS group demonstrated increases in resting CBF specifically in regions known to mediate skill learning namely, the M1, cingulate cortex, putamen, hippocampus, and cerebellum. These results indicate that TMS applied concomitantly augments behavioral effects of motor practice, with corresponding neural plasticity in motor sequence learning network. These findings are the first demonstration of the behavioral and neural enhancing effects of TMS on slow motor practice and have direct application in neurorehabilitation where TMS could be applied in conjunction with physical therapy.

Thalamic structural connectivity in medial temporal lobe epilepsy.

The thalamus has been implicated in various stages of medial temporal lobe epilepsy (MTLE) seizure evolution. The relative density and functional significance (in epileptogenesis) of thalamic projections to MTL subregions, however, remains to be determined. This study used structural and diffusion magnetic resonance imaging (MRI) to evaluate thalamic connection density with distinct MTL subregions in terms of location and volume. Nineteen MTLE patients with unilateral hippocampal sclerosis (HS; 12 right; 10 female) were compared to 19 age-matched controls. Five regions of interest (ROIs) per hemisphere were created in native space: thalamus, amygdala, entorhinal cortex, hippocampus, and parahippocampus. Separate probabilistic tractography analyses were performed between the thalamus and each ipsilateral MTL subregion (four per hemisphere). Individual connectivity profiles and regional volumes were assessed. The medial pulvinar consistently showed the highest connection density with the hippocampus in healthy controls and in MTLE patients. Decreased thalamic connected volume was observed for thalamohippocampal pathways in patients with MTLE, and indicates pathway-specific deafferentation. Regional hippocampal and thalamic atrophy was also observed, indicating gray and white matter loss in the thalamohippocampal pathway. Consistent localization of dense medial pulvinar (PuM) connectivity with the hippocampus suggests chronic PuM stimulation could modulate the MTLE seizure network. Decreased thalamic connected volume is a promising biomarker for epileptogenesis that merits longitudinal validation. A PowerPoint slide summarizing this article is available for download in the Supporting Information section here.

The effect of image quality on the assessment of child abuse photographs.

OBJECTIVES: Although child abuse pediatricians are frequently asked to evaluate risk of abuse based on photographs, the effect of photographic quality on this process is presently unknown. Photographs of abused children are often taken by professionals without photographic training, and quality varies widely. This article reports the first study of the effect of image quality on clinical assessment from photographs. METHODS: A total of 120 images depicting 60 cutaneous lesions were selected for the study. Paired images of single lesions varied in quality of focus, exposure, or framing. Seventy medical and nursing professionals were recruited from the Internet listservs focusing on child abuse. Subjects evaluated the images for quality (1-9 scale), opined if the image was "inadequate for interpretation," and answered a clinical question about the type of lesion displayed. Accuracy was defined as concordance between the subject and the live examiner's written documentation. Adequacy was defined as the proportion of subjects that did not indicate that the photograph was inadequate for interpretation. RESULTS: Mean accuracy among subjects was 64% and ranged from 35% to 84%. Accuracy was not predicted by subject profession, experience, or self-rated computer skill. Image quality and adequacy were independently associated with increased accuracy. CONCLUSIONS: Higher-quality images improved accuracy. An examiner's impression that an image is adequate did not guarantee an accurate interpretation. Reliance on photographs alone is not sufficiently accurate in the assessment of cutaneous trauma.

Electrophysiological and functional connectivity of the human supplementary motor area.

Neuro-imaging methods for detecting functional and structural inter-regional connectivity are in a rapid phase of development. While reports of regional connectivity patterns based on individual methods are becoming common, studies comparing the results of two or more connectivity-mapping methods remain rare. In this study, we applied transcranial magnetic stimulation during PET imaging (TMS/PET), a stimulation-based method, and meta-analytic connectivity modeling (MACM), a task-based method to map the connectivity patterns of the supplementary motor area (SMA). Further, we drew upon the behavioral domain meta-data of the BrainMap® database to characterize the behavioral domain specificity of two maps. Both MACM and TMS/PET detected multi-synaptic connectivity patterns, with the MACM-detected connections being more extensive. Both MACM and TMS/PET detected connections belonging to multiple behavioral domains, including action, cognition and perception. Finally, we show that the two connectivity-mapping methods are complementary in that, the MACM informed on the functional nature of SMA connections, while TMS/PET identified brain areas electrophysiologically connected with the SMA. Thus, we demonstrate that integrating multimodal database and imaging techniques can derive comprehensive connectivity maps of brain areas.

Abnormal resting state corticolimbic blood flow in depressed unmedicated patients with major depression: a (15)O-H(2)O PET study.

We investigated the differences in the resting state corticolimbic blood flow between 20 unmedicated depressed patients and 21 healthy comparisons. Resting state cerebral blood flow (CBF) was measured with H(2)(15)O PET. Anatomical MRI scans were performed on an Elscint 1.9 T Prestige system for PET-MRI coregistration. Significant changes in cerebral blood flow indicating neural activity were detected using an ROI-free image subtraction strategy. In addition, the resting blood flow in patients was correlated with the severity of depression as measured by HAM-D scores. Depressed patients showed decreases in blood flow in right anterior cingulate (Brodmann areas 24 and 32) and increased blood flow in left and right posterior cingulate (Brodmann areas 23, 29, 30), left parahippocampal gyrus (Brodmann area 36), and right caudate compared with healthy volunteers. The severity of depression was inversely correlated with the left middle and inferior frontal gyri (Brodmann areas 9 and 47) and right medial frontal gyrus (Brodmann area 10) and right anterior cingulate (Brodmann areas 24, 32) blood flow, and directly correlated with the right thalamus blood flow. These findings support previous reports of abnormalities in the resting state blood flow in the limbic-frontal structures in depressed patients compared to healthy volunteers.

A method to enhance the sensitivity of DTI analyses to group differences: a validation study with comparison to voxelwise analyses.

Studies of white matter (WM) abnormalities in psychiatric and neurological disorders often use the analysis package Tract-Based Spatial Statistics (TBSS). However, with small samples and/or subtle effects, a study using the standard TBSS approach can be underpowered. For such cases, a new method is presented that summarizes global differences between TBSS-derived fractional anisotropy (FA) images with a single paired t-statistic, estimating the degrees of freedom using spatial autocorrelation. The sensitivity of the method is demonstrated by using well-known aging effects on FA as a proxy for disease effects. Sixty healthy subjects were divided equally into younger- (YA), middle- (MA), and older-aged (OA) groups and significant global differences were demonstrated in the YA versus OA (all N ≥ 4, FA difference≈0.023), MA versus OA (all N≥4, FA difference≈0.017), and YA versus MA (FA difference≈0.005 at N=20) comparisons. In contrast, no significant difference could be detected in the YA versus MA comparison using voxelwise TBSS analysis with the full sample (N=20 per group). This method should facilitate localizing analyses in the direction of a proven group difference while providing clinically relevant information about pathophysiologic processes globally affecting WM.

Brain pathology recapitulates physiology: A network meta-analysis.

Network architecture is a brain-organizational motif present across spatial scales from cell assemblies to distributed systems. Structural pathology in some neurodegenerative disorders selectively afflicts a subset of functional networks, motivating the network degeneration hypothesis (NDH). Recent evidence suggests that structural pathology recapitulating physiology may be a general property of neuropsychiatric disorders. To test this possibility, we compared functional and structural network meta-analyses drawing upon the BrainMap database. The functional meta-analysis included results from >7,000 experiments of subjects performing >100 task paradigms; the structural meta-analysis included >2,000 experiments of patients with >40 brain disorders. Structure-function network concordance was high: 68% of networks matched (pFWE < 0.01), confirming the broader scope of NDH. This correspondence persisted across higher model orders. A positive linear association between disease and behavioral entropy (p = 0.0006;R2 = 0.53) suggests nodal stress as a common mechanism. Corroborating this interpretation with independent data, we show that metabolic 'cost' significantly differs along this transdiagnostic/multimodal gradient.

Comparison of the disparity between Talairach and MNI coordinates in functional neuroimaging data: validation of the Lancaster transform.

Spatial normalization of neuroimaging data is a standard step when assessing group effects. As a result of divergent analysis procedures due to different normalization algorithms or templates, not all published coordinates refer to the same neuroanatomical region. Specifically, the literature is populated with results in the form of MNI or Talairach coordinates, and their disparity can impede the comparison of results across different studies. This becomes particularly problematic in coordinate-based meta-analyses, wherein coordinate disparity should be corrected to reduce error and facilitate literature reviews. In this study, a quantitative comparison was performed on two corrections, the Brett transform (i.e., "mni2tal"), and the Lancaster transform (i.e., "icbm2tal"). Functional magnetic resonance imaging (fMRI) data acquired during a standard paired associates task indicated that the disparity between MNI and Talairach coordinates was better reduced via the Lancaster transform, as compared to the Brett transform. In addition, an activation likelihood estimation (ALE) meta-analysis of the paired associates literature revealed that a higher degree of concordance was obtained when using the Lancaster transform in the form of fewer, smaller, and more intense clusters. Based on these results, we recommend that the Lancaster transform be adopted as the community standard for reducing disparity between results reported as MNI or Talairach coordinates, and suggest that future spatial normalization strategies be designed to minimize this variability in the literature.

FASpecT/CT, A New SPECT/CT Acquisition With Higher Sensitivity and Efficiency in Radioiodine Thyroid Cancer Imaging.

PURPOSE: This article demonstrates the use of a new SPECT/CT acquisition protocol in patients with differentiated thyroid cancer (DTC). METHODS: SPECT/CT scans (FASpecT/CT) with fewer angle acquisitions were retrospectively reviewed in 30 DTC patients treated with radioiodine at University Hospital, San Antonio, Tex, from July 2017 to March 2019. This FASpecT/CT of 12 versus 60 to 64 sampled views for convention SPECT was made possible by iterative reconstruction. RESULTS: The FASpecT/CT protocol was judged to increase lesion detection in patients with low count rates. Furthermore, in patients with higher count rates, this technique reduced the acquisition time. FASpecT/CT patient images are shown as case examples in 4 of the 30 patients reviewed. CONCLUSIONS: This FASpecT/CT acquisition in radioiodine-treated DTC offers the potential of higher sensitivity for metastatic lymph node detection in low count rates and a significant decrease in imaging time in high count rates. These advantages make SPECT/CT imaging more acceptable for patients who have difficulty with longer imaging times, to include the pediatric population.

Lost in localization? The focus is meta-analysis.

The recent commentary by Derrfuss J, Mar RA. (2009). Lost in localization: the need for a universal coordinate database. Neuroimage, In Press proposed a universal coordinate database to archive functional neuroimaging results. In this response, we discuss our strategy in developing the BrainMap database, which was created as a mechanism to promote coordinate-based meta-analysis methods.

Protecting peer review: correspondence chronology and ethical analysis regarding Logothetis vs. Shmuel and Leopold.

Editors of scientific journals are ethically bound to provide a fair and impartial peer-review process and to protect the rights of contributing authors to publish research results. If, however, a dispute arises among investigators regarding data ownership and the right to publish, the ethical responsibilities of journal editors become more complex. The editors of Human Brain Mapping recently had the unusual experience of learning of an ongoing dispute regarding data-access rights pertaining to a manuscript already accepted for publication. Herein the editors describe the nature of the dispute, the steps taken to explore and resolve the conflict, and discuss the ethical principles that govern such circumstances. Drawing on this experience and with the goal of avoiding future controversies, the editors have formulated a Data Rights Policy and a Data Rights Procedure for Human Brain Mapping. Human Brain Mapping adopts this policy effective immediately and respectfully suggests that other journals consider adopting this or similar policies.

Amygdala hyperactivation in untreated depressed individuals.

The amygdala participates in the detection and control of affective states, and has been proposed to be a site of dysfunction in affective disorders. To assess amygdala processing in individuals with unipolar depression, we applied a functional MRI (fMRI) paradigm previously shown to be sensitive to amygdala function. Fourteen individuals with untreated DSM-IV major depression and 15 healthy subjects were studied using fMRI with a standardized emotion face recognition task. Voxel-level data sets were subjected to a multiple-regression analysis, and functionally defined regions of interest (ROI), including bilateral amygdala, were analyzed with MANOVA. Pearson correlation coefficients between amygdala activation and HAM-D score also were performed. While both depressed and healthy groups showed increased amygdala activity when viewing emotive faces compared to geometric shapes, patients with unipolar depression showed relatively more activity than healthy subjects, particularly on the left. Positive Pearson correlations between amygdala activation and HAM-D score were found for both left and right ROIs in the patient group. This study provides in vivo imaging evidence to support the hypothesis of abnormal amygdala functioning in depressed individuals.

The central sulcus: an observer-independent characterization of sulcal landmarks and depth asymmetry.

Studies of the central sulcus (CS) often use observer-dependent procedures to assess CS morphology and sulcal landmarks. Here, we applied a novel method combining automated sulcus reconstruction, surface parameterization, and an observer-independent depth measurement to study the CS. This facilitated the quantitative assessment of the spatial position and intersubject variability of several sulcal landmarks. Sulcal depth profiles also allowed us to develop an algorithm for the clear identification of several landmarks, including the pli de passage fronto-pariétal moyen (PPFM), first described by Broca. Using this algorithm, the PPFM was identified in the majority of sulci, but exhibited limited spatial variability. This appears to support Cunningham's theory that this landmark may be a developmental remnant, and may argue against its role as a guide to the more variable somatotopic hand area. Sulcal depth profiles were also utilized to assess the influence of sex, handedness, and age on CS morphology. These profiles revealed leftward depth asymmetry in the superior extent of the CS of male subjects and near the midpoint of the CS in female subjects. Age correlations were performed for these asymmetries, and a significant correlation was seen only in the male subgroup.

Perisylvian sulcal morphology and cerebral asymmetry patterns in adults who stutter.

Previous investigations of cerebral anatomy in persistent developmental stutterers have reported bilateral anomalies in the perisylvian region and atypical patterns of cerebral asymmetry. In this study, perisylvian sulcal patterns were analyzed to compare subjects with persistent developmental stuttering (PDS) and an age-, hand-, and gender-matched control group. This analysis was accomplished using software designed for 3-dimensional sulcal identification and extraction. Patterns of cerebral asymmetry were also investigated with standard planimetric measurements. PDS subjects showed a small but significant increase in both the number of sulci connecting with the second segment of the right Sylvian fissure and in the number of suprasylvian gyral banks (of sulci) along this segment. No differences were seen in the left perisylvian region for either sulcal number or gyral bank number. Measurements of asymmetry revealed typical patterns of cerebral asymmetry in both groups with no significant differences in frontal and occipital width asymmetry, frontal and occipital pole asymmetry, or planum temporale and Sylvian fissure asymmetries. The subtle difference in cortical folding of the right perisylvian region observed in PDS subjects may correlate with functional imaging studies that have reported increased right-hemisphere activity during stuttered speech.