17urn:lsid:arphahub.com:pub:8E638694-B4E0-570A-856A-746FF325BF6BResearch Ideas and OutcomesRIO2367-7163Pensoft Publishers10.3897/rio.3.e13391133917113Project ReportGradients of cortical hierarchy in AutismBethlehemRichard AIrb643@medschl.cam.ac.uk1FalkiewiczMarcel2FreybergJan3ParsonsOwen E1FarahibozorgSeyedeh-Rezvan14PretzschCharlotte3SoergelBjoern1MarguliesDaniel S2University of Cambridge, Cambridge, United KingdomUniversity of CambridgeCambridgeUnited KingdomMax Planck Institute for Human Cognitive and Brain Sciences, Leipzig, GermanyMax Planck Institute for Human Cognitive and Brain SciencesLeipzigGermanyKing's College London, London, United KingdomKing's College LondonLondonUnited KingdomMRC Cognition and Brain Sciences Unit, Cambridge, United KingdomMRC Cognition and Brain Sciences UnitCambridgeUnited Kingdom
Corresponding author: Richard AI Bethlehem (rb643@medschl.cam.ac.uk).
Academic editor:
2017020520173e1339114AE5B52-4404-5A32-8E54-BAB4AEA9251D57586024042017Richard AI Bethlehem, Marcel Falkiewicz, Jan Freyberg, Owen E Parsons, Seyedeh-Rezvan Farahibozorg, Charlotte Pretzsch, Bjoern Soergel, Daniel S MarguliesThis is an open access article distributed under the terms of the Creative Commons Attribution License 4.0 (CC-BY), which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.gradientsautismfunctional connectivityresting-state fMRIABIDEcortical hierarchyIntroduction & Aims
Autism Spectrum Disorder (ASD) is characterised by local and global disruptions of functional connectivity (FC) (Vissers et al. 2012). These disruptions are accompanied by a range of cognitive-behavioural symptoms, such as atypical socio-emotional or language processing. Several theories attempt to explain the nature of FC alterations and their impact on brain function: reduced long-distance (global) and increased short-distance (local) FC (Belmonte 2004), the 'underconnectivity hypothesis' (Just 2004) and reduced central coherence (Happé 2013). While these theories can explain local/regional anomalies in specific ASD subgroups, they do not generalize across the ASD spectrum. We propose to address these limitations with a different perspective on brain organization, based on the concept of connectivity gradients.
Connectivity gradients (Margulies et al. 2016) provide a low-dimensional approximation of whole-brain connectivity. Each gradient can be represented as a whole-brain spatial map that describes relationships between brain systems. Margulies et al. (2016) derived the gradients from healthy participants and found that the first ('principal') gradient is anchored in sensory regions on one side and transmodal regions on the other, including the default-mode network (DMN). This organization captures the major hierarchy of processing in the cortex, from sensory inputs to their abstract representations. We hypothesized that ASD might alter this hierarchy and the alteration would be reflected in the gradients. There are two ways the gradients could be affected: by changes of their scale or by the shifts of specific brain regions within the gradients. Here we explore the former option.
Hence, we set out to recreate the gradients from the orginal paper by Margulies et al. (2016) in a publically available and pre-processed autism imaging dataset (http://preprocessed-connectomes-project.org/abide/). We investigate potential differences in the scale of the gradient, its overlap with the DMN and its overlap with maps corresponding to the associated hierarchical cognitive processing as reported in figure 4 of Margulies et al. (2016).
Approach
To make the project practically feasible in the course of Brainhack, we used pre-processed male adult data (age range: 18-55) from the ABIDE dataset (Martino et al. 2013, http://preprocessed-connectomes-project.org/abide/). We chose to use the most fine-grained parcellation available (Craddock 400) that was pre-processed with the C-PAC pipeline and did not include global signal regression but did include bandpass filtering (Craddock et al. 2011). The processing pipelines for all analyses described here can be found on the Autism Research Centre github repository and, with more guidance, on the GitHub front end. In short, we selected only subjects without missing time-series to avoid dimension mismatches in the correlation matrix. We thus included 160 subjects: 91 with autism (age: 24.55±5.25) and 69 controls (26.15±7.69), matched for age (p = 0.14). We then ran diffusion embedding (Coifman and Lafon 2006) on the thresholded (>10%) matrices and backprojected the diffusion components for each individual subject using the pySTATIS package. We then saved the individual-level gradients in nifti format. We also created nifti maps of every 10th percentile of the primary gradient.
Our primary outcome measure was the scale of the gradient as estimated by the linear fit to the sorted gradient values. Individual fits were visually inspected to quantitatively assess the fit; examples can be found on the above mentioned GitHub website. As a secondary measure we computed a goodness of fit ratio for the gradient values inside and outside of brain masks obtained from NeuroSynth that accompanied the keywords listed in figure 4 of Margulies et al. 2016. For the binned principal gradient we also calculated the average z-score inside these masks in unthresholded meta-analysis maps from neurosynth for the same keywords.
Results
Gradient slopes show no distinguishable difference between the autism and neurotypical groups (Fig. 1), although the variability of the principal gradient seems somewhat larger in the autism group. The goodness of fit for the principal was compared against mask derived from thresholded NeuroSynth keywords as listed in Fig. 2. Both the autism as well as the neurotypical group show the highest goodness of fit with the default-mode network (DMN) as would be expected from Margulies et al. 2016. The overall order of the goodness of fit for the two groups did not differ. Analysis of the 10 percentile bins of the principal gradient also showed little divergence between the two groups (Fig. 3).
Limitations and future directions
In order to improve feasibility, the present study used one specific pre-processed parcellation template. This greatly improved processing time at the cost of reduced spatial specificity. The orginal paper by Margulies et al. (2016) used a voxel-based approach; it is very well possible that the parcellated resolution in the present study might not be high enough for an accurate goodness-of-fit measure. Future implementation of this pipeline will look into un-parcellated data. In addition, it is our intention to validate the current results against other parcellations of the same dataset. Furthermore, we will seek to explore other properties of the principal gradient and run more rigorous statistical tests. Finally, we will explore the second option for alterations of the gradients: shifts in the location of specific brain areas within the gradient space.
Autism is well known for its heterogeneity and it is possible that, by averaging over the entire cohort, potential subtle sub-group effects are lost. Future analysis will look into combining topographical gradient information with phenotypic information provided with ABIDE.
Acknowledgements
This work was completed during Brainhack Global - Cambridge 2017.
ReferencesBelmonteM. K.2004Autism and Abnormal Development of Brain Connectivity244292289231http://dx.doi.org/10.1523/jneurosci.3340-04.200410.1523/jneurosci.3340-04.2004CoifmanRonald R.LafonStéphane2006Diffusion maps211530http://dx.doi.org/10.1016/j.acha.2006.04.00610.1016/j.acha.2006.04.006CraddockR CameronJamesG AndrewHoltzheimerPaul EHuXiaoping PMaybergHelen S2011A whole brain fMRI atlas generated via spatially constrained spectral clustering.338191428http://dx.doi.org/10.1002/hbm.2133310.1002/hbm.21333HappéFrancesca2013The Weak Central Coherence Account of Autismhttp://dx.doi.org/10.1002/9780470939345.ch2410.1002/9780470939345.ch24JustM. A.2004Cortical activation and synchronization during sentence comprehension in high-functioning autism: evidence of underconnectivity127818111821http://dx.doi.org/10.1093/brain/awh19910.1093/brain/awh199MarguliesDaniel S.GhoshSatrajit S.GoulasAlexandrosFalkiewiczMarcelHuntenburgJulia M.LangsGeorgBezginGlebEickhoffSimon B.CastellanosF. XavierPetridesMichaelJefferiesElizabethSmallwoodJonathan2016Situating the default-mode network along a principal gradient of macroscale cortical organization113441257412579http://dx.doi.org/10.1073/pnas.160828211310.1073/pnas.1608282113MartinoA DiYanC-GLiQDenioECastellanosF XAlaertsKAndersonJ SAssafMBookheimerS YDaprettoMDeenBDelmonteSDinsteinIErtl-WagnerBFairD AGallagherLKennedyD PKeownC LKeysersCLainhartJ ELordCLunaBMenonVMinshewN JMonkC SMuellerSMüllerR-ANebelM BNiggJ TO'HearnKPelphreyK APeltierS JRudieJ DSunaertSThiouxMTyszkaJ MUddinL QVerhoevenJ SWenderothNWigginsJ LMostofskyS HMilhamM P2013The autism brain imaging data exchange: towards a large-scale evaluation of the intrinsic brain architecture in autism196659667http://dx.doi.org/10.1038/mp.2013.7810.1038/mp.2013.78VissersMarlies E.CohenMichael XGeurtsHilde M.2012Brain connectivity and high functioning autism: A promising path of research that needs refined models, methodological convergence, and stronger behavioral links361604625http://dx.doi.org/10.1016/j.neubiorev.2011.09.00310.1016/j.neubiorev.2011.09.003
Linear fit of gradient slopes for the top 10 gradients for each group (1 == neurotypical individuals; 2 == individuals with autism).
Goodness of fit for principal gradient. Ranked according to the median goodness of fit for both groups.
Average z-scores obtained from unthresholded reverse-inference neurosynth meta-analysis activation maps. We divided the principle gradient into percentiles and used this to create a mask. We then calculated the average z-score inside this mask on neurosynth maps.