Automating the Human Connectome Project's Temporal ICA Pipeline.

Functional magnetic resonance imaging (fMRI) data are dominated by noise and artifacts, with only a small fraction of the variance relating to neural activity. Temporal independent component analysis (tICA) is a recently developed method that enables selective denoising of fMRI artifacts related to physiology such as respiration. However, an automated and easy to use pipeline for tICA has not previously been available; instead, two manual steps have been necessary: 1) setting the group spatial ICA dimensionality after MELODIC's Incremental Group-PCA (MIGP) and 2) labeling tICA components as artifacts versus signals. Moreover, guidance has been lacking as to how many subjects and timepoints are needed to adequately re-estimate the temporal ICA decomposition and what alternatives are available for smaller groups or even individual subjects. Here, we introduce a nine-step fully automated tICA pipeline which removes global artifacts from fMRI dense timeseries after sICA+FIX cleaning and MSMAll alignment driven by functionally relevant areal features. Additionally, we have developed an automated "reclean" Pipeline for improved spatial ICA (sICA) artifact removal. Two major automated components of the pipeline are 1) an automatic group spatial ICA (sICA) dimensionality selection for MIGP data enabled by fitting multiple Wishart distributions; 2) a hierarchical classifier to distinguish group tICA signal components from artifactual components, equipped with a combination of handcrafted features from domain expert knowledge and latent features obtained via self-supervised learning on spatial maps. We demonstrate that the dimensionality estimated for the MIGP data from HCP Young Adult 3T and 7T datasets is comparable to previous manual tICA estimates, and that the group sICA decomposition is highly reproducible. We also show that the tICA classifier achieved over 0.98 Precision-Recall Area Under Curve (PR-AUC) and that the correctly classified components account for over 95% of the tICA-represented variance on multiple held-out evaluation datasets including the HCP-Young Adult, HCP-Aging and HCP-Development datasets under various settings. Our automated tICA pipeline is now available as part of the HCP pipelines, providing a powerful and user-friendly tool for the neuroimaging community.

QuNex-An integrative platform for reproducible neuroimaging analytics.

Introduction: Neuroimaging technology has experienced explosive growth and transformed the study of neural mechanisms across health and disease. However, given the diversity of sophisticated tools for handling neuroimaging data, the field faces challenges in method integration, particularly across multiple modalities and species. Specifically, researchers often have to rely on siloed approaches which limit reproducibility, with idiosyncratic data organization and limited software interoperability. Methods: To address these challenges, we have developed Quantitative Neuroimaging Environment & Toolbox (QuNex), a platform for consistent end-to-end processing and analytics. QuNex provides several novel functionalities for neuroimaging analyses, including a "turnkey" command for the reproducible deployment of custom workflows, from onboarding raw data to generating analytic features. Results: The platform enables interoperable integration of multi-modal, community-developed neuroimaging software through an extension framework with a software development kit (SDK) for seamless integration of community tools. Critically, it supports high-throughput, parallel processing in high-performance compute environments, either locally or in the cloud. Notably, QuNex has successfully processed over 10,000 scans across neuroimaging consortia, including multiple clinical datasets. Moreover, QuNex enables integration of human and non-human workflows via a cohesive translational platform. Discussion: Collectively, this effort stands to significantly impact neuroimaging method integration across acquisition approaches, pipelines, datasets, computational environments, and species. Building on this platform will enable more rapid, scalable, and reproducible impact of neuroimaging technology across health and disease.

Empirical transmit field bias correction of T1w/T2w myelin maps.

T1-weighted divided by T2-weighted (T1w/T2w) myelin maps were initially developed for neuroanatomical analyses such as identifying cortical areas, but they are increasingly used in statistical comparisons across individuals and groups with other variables of interest. Existing T1w/T2w myelin maps contain radiofrequency transmit field (B1+) biases, which may be correlated with these variables of interest, leading to potentially spurious results. Here we propose two empirical methods for correcting these transmit field biases using either explicit measures of the transmit field or alternatively a 'pseudo-transmit' approach that is highly correlated with the transmit field at 3T. We find that the resulting corrected T1w/T2w myelin maps are both better neuroanatomical measures (e.g., for use in cross-species comparisons), and more appropriate for statistical comparisons of relative T1w/T2w differences across individuals and groups (e.g., sex, age, or body-mass-index) within a consistently acquired study at 3T. We recommend that investigators who use the T1w/T2w approach for mapping cortical myelin use these B1+ transmit field corrected myelin maps going forward.

Towards HCP-Style macaque connectomes: 24-Channel 3T multi-array coil, MRI sequences and preprocessing.

Macaque monkeys are an important animal model where invasive investigations can lead to a better understanding of the cortical organization of primates including humans. However, the tools and methods for noninvasive image acquisition (e.g. MRI RF coils and pulse sequence protocols) and image data preprocessing have lagged behind those developed for humans. To resolve the structural and functional characteristics of the smaller macaque brain, high spatial, temporal, and angular resolutions combined with high signal-to-noise ratio are required to ensure good image quality. To address these challenges, we developed a macaque 24-channel receive coil for 3-T MRI with parallel imaging capabilities. This coil enables adaptation of the Human Connectome Project (HCP) image acquisition protocols to the in-vivo macaque brain. In addition, we adapted HCP preprocessing methods to the macaque brain, including spatial minimal preprocessing of structural, functional MRI (fMRI), and diffusion MRI (dMRI). The coil provides the necessary high signal-to-noise ratio and high efficiency in data acquisition, allowing four- and five-fold accelerations for dMRI and fMRI. Automated FreeSurfer segmentation of cortex, reconstruction of cortical surface, removal of artefacts and nuisance signals in fMRI, and distortion correction of dMRI all performed well, and the overall quality of basic neurobiological measures was comparable with those for the HCP. Analyses of functional connectivity in fMRI revealed high sensitivity as compared with those from publicly shared datasets. Tractography-based connectivity estimates correlated with tracer connectivity similarly to that achieved using ex-vivo dMRI. The resulting HCP-style in vivo macaque MRI data show considerable promise for analyzing cortical architecture and functional and structural connectivity using advanced methods that have previously only been available in studies of the human brain.

Toblerone: Surface-Based Partial Volume Estimation.

Partial volume effects (PVE) present a source of confound for the analysis of functional imaging data. Correction for PVE requires estimates of the partial volumes (PVs) present in an image. These estimates are conventionally obtained via volumetric segmentation, but such an approach may not be accurate for complex structures such as the cortex. An alternative is to use surface-based segmentation, which is well-established within the literature. Toblerone is a new method for estimating PVs using such surfaces. It uses a purely geometric approach that considers the intersection between a surface and the voxels of an image. In contrast to existing surface-based techniques, Toblerone is not restricted to use with any particular structure or modality. Evaluation in a neuroimaging context has been performed on simulated surfaces, simulated T1-weighted MRI images and finally a Human Connectome Project test-retest dataset. A comparison has been made to two existing surface-based methods; in all analyses Toblerone's performance either matched or surpassed the comparator methods. Evaluation results also show that compared to an existing volumetric method (FSL FAST), a surface-based approach with Toblerone offers improved robustness to scanner noise and field non-uniformity, and better inter-session repeatability in brain volume. In contrast to volumetric methods, a surface-based approach negates the need to perform resampling which is advantageous at the resolutions typically used for neuroimaging.

Cerebral cortical folding, parcellation, and connectivity in humans, nonhuman primates, and mice.

Advances in neuroimaging and neuroanatomy have yielded major insights concerning fundamental principles of cortical organization and evolution, thus speaking to how well different species serve as models for human brain function in health and disease. Here, we focus on cortical folding, parcellation, and connectivity in mice, marmosets, macaques, and humans. Cortical folding patterns vary dramatically across species, and individual variability in cortical folding increases with cortical surface area. Such issues are best analyzed using surface-based approaches that respect the topology of the cortical sheet. Many aspects of cortical organization can be revealed using 1 type of information (modality) at a time, such as maps of cortical myelin content. However, accurate delineation of the entire mosaic of cortical areas requires a multimodal approach using information about function, architecture, connectivity, and topographic organization. Comparisons across the 4 aforementioned species reveal dramatic differences in the total number and arrangement of cortical areas, particularly between rodents and primates. Hemispheric variability and bilateral asymmetry are most pronounced in humans, which we evaluated using a high-quality multimodal parcellation of hundreds of individuals. Asymmetries include modest differences in areal size but not in areal identity. Analyses of cortical connectivity using anatomical tracers reveal highly distributed connectivity and a wide range of connection weights in monkeys and mice; indirect measures using functional MRI suggest a similar pattern in humans. Altogether, a multifaceted but integrated approach to exploring cortical organization in primate and nonprimate species provides complementary advantages and perspectives.

Ciftify: A framework for surface-based analysis of legacy MR acquisitions.

The preprocessing pipelines of the Human Connectome Project (HCP) were made publicly available for the neuroimaging community to apply the HCP analytic approach to data from non-HCP sources. The HCP analytic approach is surface-based for the cerebral cortex, uses the CIFTI "grayordinate" file format, provides greater statistical sensitivity than traditional volume-based analysis approaches, and allows for a more neuroanatomically-faithful representation of data. However, the HCP pipelines require the acquisition of specific images (namely T2w and field map) that historically have often not been acquired. Massive amounts of this 'legacy' data could benefit from the adoption of HCP-style methods. However, there is currently no published framework, to our knowledge, for adapting HCP preprocessing to "legacy" data. Here we present the ciftify project, a parsimonious analytic framework for adapting key modules from the HCP pipeline into existing structural workflows using FreeSurfer's recon_all structural and existing functional preprocessing workflows. Within this framework, any functional dataset with an accompanying (i.e. T1w) anatomical data can be analyzed in CIFTI format. To simplify usage for new data, the workflow has been bundled with fMRIPrep following the BIDS-app framework. Finally, we present the package and comment on future neuroinformatics advances that may accelerate the movement to a CIFTI-based grayordinate framework.

Classification of temporal ICA components for separating global noise from fMRI data: Reply to Power.

We respond to a critique of our temporal Independent Components Analysis (ICA) method for separating global noise from global signal in fMRI data that focuses on the signal versus noise classification of several components. While we agree with several of Power's comments, we provide evidence and analysis to rebut his major criticisms and to reassure readers that temporal ICA remains a powerful and promising denoising approach.

The Human Connectome Project 7 Tesla retinotopy dataset: Description and population receptive field analysis.

About a quarter of human cerebral cortex is dedicated mainly to visual processing. The large-scale spatial organization of visual cortex can be measured with functional magnetic resonance imaging (fMRI) while subjects view spatially modulated visual stimuli, also known as "retinotopic mapping." One of the datasets collected by the Human Connectome Project involved ultrahigh-field (7 Tesla) fMRI retinotopic mapping in 181 healthy young adults (1.6-mm resolution), yielding the largest freely available collection of retinotopy data. Here, we describe the experimental paradigm and the results of model-based analysis of the fMRI data. These results provide estimates of population receptive field position and size. Our analyses include both results from individual subjects as well as results obtained by averaging fMRI time series across subjects at each cortical and subcortical location and then fitting models. Both the group-average and individual-subject results reveal robust signals across much of the brain, including occipital, temporal, parietal, and frontal cortex as well as subcortical areas. The group-average results agree well with previously published parcellations of visual areas. In addition, split-half analyses show strong within-subject reliability, further demonstrating the high quality of the data. We make publicly available the analysis results for individual subjects and the group average, as well as associated stimuli and analysis code. These resources provide an opportunity for studying fine-scale individual variability in cortical and subcortical organization and the properties of high-resolution fMRI. In addition, they provide a set of observations that can be compared with other Human Connectome Project measures acquired in these same participants.

Extending the Human Connectome Project across ages: Imaging protocols for the Lifespan Development and Aging projects.

The Human Connectome Projects in Development (HCP-D) and Aging (HCP-A) are two large-scale brain imaging studies that will extend the recently completed HCP Young-Adult (HCP-YA) project to nearly the full lifespan, collecting structural, resting-state fMRI, task-fMRI, diffusion, and perfusion MRI in participants from 5 to 100+ years of age. HCP-D is enrolling 1300+ healthy children, adolescents, and young adults (ages 5-21), and HCP-A is enrolling 1200+ healthy adults (ages 36-100+), with each study collecting longitudinal data in a subset of individuals at particular age ranges. The imaging protocols of the HCP-D and HCP-A studies are very similar, differing primarily in the selection of different task-fMRI paradigms. We strove to harmonize the imaging protocol to the greatest extent feasible with the completed HCP-YA (1200+ participants, aged 22-35), but some imaging-related changes were motivated or necessitated by hardware changes, the need to reduce the total amount of scanning per participant, and/or the additional challenges of working with young and elderly populations. Here, we provide an overview of the common HCP-D/A imaging protocol including data and rationales for protocol decisions and changes relative to HCP-YA. The result will be a large, rich, multi-modal, and freely available set of consistently acquired data for use by the scientific community to investigate and define normative developmental and aging related changes in the healthy human brain.

Construction of a neonatal cortical surface atlas using Multimodal Surface Matching in the Developing Human Connectome Project.

We propose a method for constructing a spatio-temporal cortical surface atlas of neonatal brains aged between 36 and 44 weeks of post-menstrual age (PMA) at the time of scan. The data were acquired as part of the Developing Human Connectome Project (dHCP), and the constructed surface atlases are publicly available. The method is based on a spherical registration approach: Multimodal Surface Matching (MSM), using cortical folding for driving the alignment. Templates have been generated for the anatomical cortical surface and for the cortical feature maps: sulcal depth, curvature, thickness, T1w/T2w myelin maps and cortical regions. To achieve this, cortical surfaces from 270 infants were first projected onto the sphere. Templates were then generated in two stages: first, a reference space was initialised via affine alignment to a group average adult template. Following this, templates were iteratively refined through repeated alignment of individuals to the template space until the variability of the average feature sets converged. Finally, bias towards the adult reference was removed by applying the inverse of the average affine transformations on the template and de-drifting the template. We used temporal adaptive kernel regression to produce age-dependant atlases for 9 weeks (36-44 weeks PMA). The generated templates capture expected patterns of cortical development including an increase in gyrification as well as an increase in thickness and T1w/T2w myelination with increasing age.

The impact of traditional neuroimaging methods on the spatial localization of cortical areas.

Localizing human brain functions is a long-standing goal in systems neuroscience. Toward this goal, neuroimaging studies have traditionally used volume-based smoothing, registered data to volume-based standard spaces, and reported results relative to volume-based parcellations. A novel 360-area surface-based cortical parcellation was recently generated using multimodal data from the Human Connectome Project, and a volume-based version of this parcellation has frequently been requested for use with traditional volume-based analyses. However, given the major methodological differences between traditional volumetric and Human Connectome Project-style processing, the utility and interpretability of such an altered parcellation must first be established. By starting from automatically generated individual-subject parcellations and processing them with different methodological approaches, we show that traditional processing steps, especially volume-based smoothing and registration, substantially degrade cortical area localization compared with surface-based approaches. We also show that surface-based registration using features closely tied to cortical areas, rather than to folding patterns alone, improves the alignment of areas, and that the benefits of high-resolution acquisitions are largely unexploited by traditional volume-based methods. Quantitatively, we show that the most common version of the traditional approach has spatial localization that is only 35% as good as the best surface-based method as assessed using two objective measures (peak areal probabilities and "captured area fraction" for maximum probability maps). Finally, we demonstrate that substantial challenges exist when attempting to accurately represent volume-based group analysis results on the surface, which has important implications for the interpretability of studies, both past and future, that use these volume-based methods.

Using temporal ICA to selectively remove global noise while preserving global signal in functional MRI data.

Temporal fluctuations in functional Magnetic Resonance Imaging (fMRI) have been profitably used to study brain activity and connectivity for over two decades. Unfortunately, fMRI data also contain structured temporal "noise" from a variety of sources, including subject motion, subject physiology, and the MRI equipment. Recently, methods have been developed to automatically and selectively remove spatially specific structured noise from fMRI data using spatial Independent Components Analysis (ICA) and machine learning classifiers. Spatial ICA is particularly effective at removing spatially specific structured noise from high temporal and spatial resolution fMRI data of the type acquired by the Human Connectome Project and similar studies. However, spatial ICA is mathematically, by design, unable to separate spatially widespread "global" structured noise from fMRI data (e.g., blood flow modulations from subject respiration). No methods currently exist to selectively and completely remove global structured noise while retaining the global signal from neural activity. This has left the field in a quandary-to do or not to do global signal regression-given that both choices have substantial downsides. Here we show that temporal ICA can selectively segregate and remove global structured noise while retaining global neural signal in both task-based and resting state fMRI data. We compare the results before and after temporal ICA cleanup to those from global signal regression and show that temporal ICA cleanup removes the global positive biases caused by global physiological noise without inducing the network-specific negative biases of global signal regression. We believe that temporal ICA cleanup provides a "best of both worlds" solution to the global signal and global noise dilemma and that temporal ICA itself unlocks interesting neurobiological insights from fMRI data.

Dynamic patterns of cortical expansion during folding of the preterm human brain.

During the third trimester of human brain development, the cerebral cortex undergoes dramatic surface expansion and folding. Physical models suggest that relatively rapid growth of the cortical gray matter helps drive this folding, and structural data suggest that growth may vary in both space (by region on the cortical surface) and time. In this study, we propose a unique method to estimate local growth from sequential cortical reconstructions. Using anatomically constrained multimodal surface matching (aMSM), we obtain accurate, physically guided point correspondence between younger and older cortical reconstructions of the same individual. From each pair of surfaces, we calculate continuous, smooth maps of cortical expansion with unprecedented precision. By considering 30 preterm infants scanned two to four times during the period of rapid cortical expansion (28-38 wk postmenstrual age), we observe significant regional differences in growth across the cortical surface that are consistent with the emergence of new folds. Furthermore, these growth patterns shift over the course of development, with noninjured subjects following a highly consistent trajectory. This information provides a detailed picture of dynamic changes in cortical growth, connecting what is known about patterns of development at the microscopic (cellular) and macroscopic (folding) scales. Since our method provides specific growth maps for individual brains, we are also able to detect alterations due to injury. This fully automated surface analysis, based on tools freely available to the brain-mapping community, may also serve as a useful approach for future studies of abnormal growth due to genetic disorders, injury, or other environmental variables.

Neurite imaging reveals microstructural variations in human cerebral cortical gray matter.

We present distinct patterns of neurite distribution in the human cerebral cortex using diffusion magnetic resonance imaging (MRI). We analyzed both high-resolution structural (T1w and T2w images) and diffusion MRI data in 505 subjects from the Human Connectome Project. Neurite distributions were evaluated using the neurite orientation dispersion and density imaging (NODDI) model, optimized for gray matter, and mapped onto the cortical surface using a method weighted towards the cortical mid-thickness to reduce partial volume effects. The estimated neurite density was high in both somatosensory and motor areas, early visual and auditory areas, and middle temporal area (MT), showing a strikingly similar distribution to myelin maps estimated from the T1w/T2w ratio. The estimated neurite orientation dispersion was particularly high in early sensory areas, which are known for dense tangential fibers and are classified as granular cortex by classical anatomists. Spatial gradients of these cortical neurite properties revealed transitions that colocalize with some areal boundaries in a recent multi-modal parcellation of the human cerebral cortex, providing mutually supportive evidence. Our findings indicate that analyzing the cortical gray matter neurite morphology using diffusion MRI and NODDI provides valuable information regarding cortical microstructure that is related to but complementary to myeloarchitecture.

Multimodal surface matching with higher-order smoothness constraints.

In brain imaging, accurate alignment of cortical surfaces is fundamental to the statistical sensitivity and spatial localisation of group studies, and cortical surface-based alignment has generally been accepted to be superior to volume-based approaches at aligning cortical areas. However, human subjects have considerable variation in cortical folding, and in the location of functional areas relative to these folds. This makes alignment of cortical areas a challenging problem. The Multimodal Surface Matching (MSM) tool is a flexible, spherical registration approach that enables accurate registration of surfaces based on a variety of different features. Using MSM, we have previously shown that driving cross-subject surface alignment, using areal features, such as resting state-networks and myelin maps, improves group task fMRI statistics and map sharpness. However, the initial implementation of MSM's regularisation function did not penalize all forms of surface distortion evenly. In some cases, this allowed peak distortions to exceed neurobiologically plausible limits, unless regularisation strength was increased to a level which prevented the algorithm from fully maximizing surface alignment. Here we propose and implement a new regularisation penalty, derived from physically relevant equations of strain (deformation) energy, and demonstrate that its use leads to improved and more robust alignment of multimodal imaging data. In addition, since spherical warps incorporate projection distortions that are unavoidable when mapping from a convoluted cortical surface to the sphere, we also propose constraints that enforce smooth deformation of cortical anatomies. We test the impact of this approach for longitudinal modelling of cortical development for neonates (born between 31 and 43 weeks of post-menstrual age) and demonstrate that the proposed method increases the biological interpretability of the distortion fields and improves the statistical significance of population-based analysis relative to other spherical methods.

The relationship of cortical folding and brain arteriovenous malformations.

BACKGROUND: The pathogenesis of human intracranial arteriovenous malformations (AVMs) is not well understood; this study aims to quantitatively assess cortical folding in patients with these lesions. METHODS: Seven adult participants, 4 male and 3 female, with unruptured, surgically unresectable intracranial AVMs were prospectively enrolled in the study, with a mean age of 42.1 years and Spetzler-Martin grade range of II-IV. High-resolution brain MRI T1 and T2 sequences were obtained. After standard preprocessing, segmentation and registration techniques, three measures of cortical folding, the depth difference index (DDI), coordinate distance index (CDI) and gyrification index (GI)), were calculated for the affected and unaffected hemispheres of each subject as well as a healthy control subject set. RESULTS: Of the three metrics, CDI, DDI and GI, used for cortical folding assessment, none demonstrated significant differences between the participants and previously studied healthy adults. There was a significant negative correlation between the DDI ratio between affected and unaffected hemispheres and AVM volume (correlation coefficient r = -0.74, p = 0.04). CONCLUSION: This study is the first to quantitatively assess human brain cortical folding in the presence of intracranial AVMs and no significant differences between AVM-affected versus unaffected hemispheres were found in a small dataset. We suggest longitudinal, larger human MRI-based cortical folding studies to assess whether AVMs are congenital lesions of vascular development or de novo, dynamic lesions.

The Human Connectome Project's neuroimaging approach.

Noninvasive human neuroimaging has yielded many discoveries about the brain. Numerous methodological advances have also occurred, though inertia has slowed their adoption. This paper presents an integrated approach to data acquisition, analysis and sharing that builds upon recent advances, particularly from the Human Connectome Project (HCP). The 'HCP-style' paradigm has seven core tenets: (i) collect multimodal imaging data from many subjects; (ii) acquire data at high spatial and temporal resolution; (iii) preprocess data to minimize distortions, blurring and temporal artifacts; (iv) represent data using the natural geometry of cortical and subcortical structures; (v) accurately align corresponding brain areas across subjects and studies; (vi) analyze data using neurobiologically accurate brain parcellations; and (vii) share published data via user-friendly databases. We illustrate the HCP-style paradigm using existing HCP data sets and provide guidance for future research. Widespread adoption of this paradigm should accelerate progress in understanding the brain in health and disease.

A multi-modal parcellation of human cerebral cortex.

Understanding the amazingly complex human cerebral cortex requires a map (or parcellation) of its major subdivisions, known as cortical areas. Making an accurate areal map has been a century-old objective in neuroscience. Using multi-modal magnetic resonance images from the Human Connectome Project (HCP) and an objective semi-automated neuroanatomical approach, we delineated 180 areas per hemisphere bounded by sharp changes in cortical architecture, function, connectivity, and/or topography in a precisely aligned group average of 210 healthy young adults. We characterized 97 new areas and 83 areas previously reported using post-mortem microscopy or other specialized study-specific approaches. To enable automated delineation and identification of these areas in new HCP subjects and in future studies, we trained a machine-learning classifier to recognize the multi-modal 'fingerprint' of each cortical area. This classifier detected the presence of 96.6% of the cortical areas in new subjects, replicated the group parcellation, and could correctly locate areas in individuals with atypical parcellations. The freely available parcellation and classifier will enable substantially improved neuroanatomical precision for studies of the structural and functional organization of human cerebral cortex and its variation across individuals and in development, aging, and disease.

Analysis of cortical shape in children with simplex autism.

We used surface-based morphometry to test for differences in cortical shape between children with simplex autism (n = 34, mean age 11.4 years) and typical children (n = 32, mean age 11.3 years). This entailed testing for group differences in sulcal depth and in 3D coordinates after registering cortical midthickness surfaces to an atlas target using 2 independent registration methods. We identified bilateral differences in sulcal depth in restricted portions of the anterior-insula and frontal-operculum (aI/fO) and in the temporoparietal junction (TPJ). The aI/fO depth differences are associated with and likely to be caused by a shape difference in the inferior frontal gyrus in children with simplex autism. Comparisons of average midthickness surfaces of children with simplex autism and those of typical children suggest that the significant sulcal depth differences represent local peaks in a larger pattern of regional differences that are below statistical significance when using coordinate-based analysis methods. Cortical regions that are statistically significant before correction for multiple measures are peaks of more extended, albeit subtle regional differences that may guide hypothesis generation for studies using other imaging modalities.