Abstract
The canonical theory emphasizes fronto-parietal network (FPN) is key in mediating general fluid intelligence (gF). Meanwhile, recent studies show that multiple sensory regions in occipito-temporal border also play a key role in gF. However, the underlying mechanism is not yet clear. To investigate this issue, this study selects human MT complex (MT+), a region locates at the occipito-temporal border representing multiple sensory flows as a target brain area. Using ultra-high field magnetic resonance spectroscopy (MRS) to measure GABA/glutamate concentrations in MT+ combining resting-state fMRI functional connectivity (FC), behavioral examinations including MT+ perception suppression test and gF subtest in visuo-spatial component, we reveal that MT+ GABA and frontal-MT+ FC significantly correlate with the performance of visuo-spatial intelligence. Further, serial mediation model demonstrates that MT+ GABA predicting visuo-spatial gF fully mediated by reverberation effect between frontal and MT+ network. Our finding highlights that sensory cortex could integrate into complex cognition system as an intellectual hub.
eLife assessment
This study employed a comprehensive approach to examining how the MT+ region integrates into a complex cognition system in mediating human visuo-spatial intelligence. While the findings are useful, the experimental evidence is incomplete and the study design, hypothesis, analyses, writing, and presentation need to be improved. The work will be of interest to researchers in psychology, cognitive science, and neuroscience.
Introduction
General fluid intelligence (gF) is a current problem-solving ability, which shows high inter-individual differences in humans(Cattell & Raymond, 1963). At the beginning of the last century, Spearman(Spearman, 1904) proposed that some general or g factor contributes to our cognition. Visuo-spatial intelligence, usually tested by visual materials, is considered to have high g-loading(Colom et al., 2006; Deary et al., 2010; Jung & Haier, 2007). Recently, Duncan et al. show coding of gF in distributed regions including a new multiple-demand (MD) core, the multiple sensory regions at the temporo-occipital border(Assem et al., 2020; Duncan et al., 2020). However, more evidence is needed to uncover how the sensory cortex is a genuine cognitive core. The human extrastriate cortical region, middle temporal complex (MT+), locates at the temporo-occipital border(Dumoulin et al., 2000), and is a key region in the multiple representations of sensory flows (including optic, tactile, and auditory flows)(Bedny et al., 2010; Ricciardi et al., 2007); this ideally suits it to be a new MD core. To further investigate this issue, this work conducted multi-level examination including biochemical (glutamatergic - gabaergic in MT+), regional-systemic (brain connectivity with MT+ - based), and behavioral (visual motion function in MT+) levels to reveal if the function of MT+ contributes to the visuo-spatial component of gF.
Using a well-known visual motion paradigm (center-surround antagonism)(Liu et al., 2016; Tadin et al., 2003), Melnick et al found a strong link between suppression index (SI) of motion perception and the scores of the block design test (BDT, a subtest of the Wechsler Adult Intelligence Scale (WASI), which measures the visuo-spatial component (3D domain) of gF (Melnick et al., 2013). Based on these findings, MT+ is likely to play a key role in MD system of gF. Strong connectivity between the core regions suggests a mechanism for large-scale information exchange and integration in the brain(Barbey, 2018; Cole et al., 2012). This potential conjunctive coding may overlap with the inhibition and/or excitation mechanism of MT+. Taken together, we hypothesized that 3D visuo-spatial intelligence (the performance of BDT) might be predicted by the inhibitory and/or excitation mechanisms in MT+ and the integrative functions connecting MT+ with frontal cortex (Figure 1a). To verify the specificity of MT+, we used primary visual cortex (V1) - based GABA/Glu as control as it mediates the 2D rather than 3D visual domain(Born & Bradley, 2005).
![](https://prod--epp.elifesciences.org/iiif/2/97545%2Fv1%2Fcontent%2F585863v1_fig1.tif/full/max/0/default.jpg)
A priori hypothesis and experimental design. (a) Schematic of a priori hypothesis. The inhibition mechanism centered on MT+ GABA, including molecular level - the GABAergic inhibition in MT+ (green circle), brain connectivity of MT+ - frontal relating to MT+ GABA (blue circle), and perceptual suppression of visual motion (the behavioral presentation of MT+ function, red circle), contributes to the visuo-spatial component of gF (3D domain, yellow circle). (b) Schematic of experimental design. Session 1 (rectangle box of short line) was the functional MRI and MRS scanning at resting state. Session 2 (rectangle box of solid line) was another region of MRS acquisition. In the two sessions, the order of MRS scanning regions (MT+ and V1) was counter balanced across participants. There was a structure MRI scanning before each MRS data acquisition. Interval between the two sessions was used for behavioral measurement (rectangle box of dotted line): block design task (BDT) and psychophysical task - motion discrimination. Sold lines indicate the experiment sequence.
We employ ultra-high field (7T) magnetic resonance spectroscopy (MRS) technology to reliably resolve GABA and Glu concentrations(Ende, 2015; Liu et al., 2022; Song et al., 2021), and first demonstrate GABAergic inhibition mechanisms (but not excitatory Glu) in MT+ region involve in the 3D visuo-spatial ability. Further, analysis of functional brain connectivity at rest reveals that the network (between MT+ and frontal cortex) relating to MT+ GABA and perceptual suppression contribute the visuo-spatial intelligence. Our results provide direct evidence that inhibitory mechanisms centered on GABA levels in MT+ region (a sensory cortex) mediate multi-level visuo-spatial component (3D domain) of gF thus drawing a direct connection of biochemistry, brain connectivity, and behavioral levels.
Results
To determine whether the function of MT+ cortex contributes to visuo-spatial component (3D domain) of gF, we adopted the experimental design depicted in Figure 1b. Participants underwent two MRI sessions: the first encompassing resting-state fMRI and magnetic resonance spectra (MRS), and the second solely involving MRS. A 30-minute interval separated these sessions, during which participants performed motion discrimination tasks (using center-surround antagonism stimuli)(Tadin et al., 2003) and the block design test (BDT), which assesses the visuo-spatial ability (3D domain) of gF(Fangmeier et al., 2006). Both the volume-of-interests (VOIs) of MRS scanning in the left MT+ (targeted brain area) and the left primary visual cortex (V1, control brain area) had dimensions of 2×2×2 cm3, and the MRS scanning sequences were randomized across the two sessions. The MT+ MRS VOIs were demarcated using an anatomical landmark(Dumoulin et al., 2000). For 14 subjects, we also utilized fMRI to functionally pinpoint the MT+ to validate the placement of the VOI (Figure 2a, b). The V1 MRS VOIs were anatomically defined (Methods). Here, MRS data after extensive quality control (31/36 in MT+, and 28/36 in V1) were taken for further analysis (Methods).
![](https://prod--epp.elifesciences.org/iiif/2/97545%2Fv1%2Fcontent%2F585863v1_fig2.tif/full/max/0/default.jpg)
MT+ localizer scans and MT+ MRS VOI placement. (a) Single task block designs. First: a cross fixation on the center of the screen (10s). Second: a moving grating (2°) toward left last 10s. Third: the grating keeps static for 10s. Fourth: the grating moves toward right last 10s. The localizer scans consist of 8 blocks. (b) MT+ location and MRS VOI placement. The upper template is the horizontal view. The lower templates from left to right are coronal and sagittal views. The warm color indicates the overlap of fMRI activation of MT+ across 14 subjects, the cold color bar indicates the overlap of MRS VOIs across all subjects.
GABA and Glu concentrations in MT+ and V1 and their relation to SI and BDT
An example of MRS from voxel located in MT+ is shown in Figure 3a. LCModel fittings for GABA spectra from all subjects in MT+ (n = 31) and V1 (n = 28) are illustrated in Figure 3b (color scale presents the BDT scores). We discerned a significant association between the inter-subjects’ BDT scores and the GABA levels in MT+ voxels, but not in V1 voxels. Quantitative analysis displayed that BDT significantly correlates with GABA concentrations in MT+ voxels (r = 0.39, P = 0.03, n = 31, Figure 3c). After using partial correlation to control for the effect of age, the relationship remains significant (rpartial = 0.43 P = 0.02, 1 participant excluded due to the age greater than mean + 2.5SD) In contrast, there was no obvious correlation between BDT and GABA levels in V1 voxels (figure supplement 1a). We show that SI significantly correlates with GABA levels in MT+ voxels (r = 0.44, P = 0.01, n = 31, Figure 3d). In contrast, no significant correlation between SI and GABA concentrations in V1 voxels was observed (figure supplement 1b). This finding supports the hypothesis that motion perception is associated with neural activity in MT+ area, but not in V1(Schallmo et al., 2018). LCModel fittings for Glu spectra from all subjects in MT+ (n = 31) and V1 (n = 28) voxels are presented in figure supplement 2a.
![](https://prod--epp.elifesciences.org/iiif/2/97545%2Fv1%2Fcontent%2F585863v1_fig3.tif/full/max/0/default.jpg)
MRS spectra and the relationships between GABA levels and SI / BDT. (a) Example spectrum from the MT+ voxel of one participant. The first line is the LCModel fitting result of all metabolites, and the following lines show the Glu and GABA spectra fitting with LCModel, and then the baseline. (b) Individual participants fitted GABA MRS spectra from the MT+ (top) and V1 (bottom) voxels from baseline measurement. The colors of the GABA spectra represent the individual differences of BDT. The color bar represents the scores of BDT. (c) and (d) Pearson’s correlations showing significant positive correlations between MT+ GABA and BDT scores (c), between MT+ GABA and SI (d). The ribbon between dotted lines represents the 95% confidence interval, and the black regression line represents the Pearson’s correlation coefficient (r). GABA and Glu concentrations (Conc.) are absolute, with units of mmol per kg wet weight (Methods).
Unlike in the case of GABA, no significant correlations between BDT and Glu levels were found in both MT+ and V1 voxels (figure supplement 2b, c). While, as expected(Song et al., 2021), we observed significant positive correlations between GABA and Glu concentrations in both MT+ (r = 0.62, P < 0.001, n = 31) and V1 voxels (r = 0.56, P = 0.002, n = 28) (figure supplement 3a, b). Additionally, a significant correlation between SI and BDT was discerned (r = 0.51, P = 0.002, n = 34, figure supplement 4a, rpartial = 0.67, P < 0.001, 1 participant excluded due to the age greater than mean + 2.5SD), corroborating previous conclusions(Melnick et al., 2013). Two outliers evident in Figure 3d were excluded, with consistent results depicted in figure supplement 4b.
MT - frontal FC relates to SI and BDT
We next took the left MT+ as the seed region and separately measured interregional FCs between the seed region and each voxel in the frontal regions (a priori search space). These measurements were correlated with performance in 3D visuospatial ability (BDT) to identify FCs with significant correlations. Results from connectivity-BDT analysis are summarized in Table 1 and shown in Figure 4a. We found that brain regions with FC strength to the seed region (left MT+) significantly correlated with BDT scores were situated within the canonical cognitive cores of MD system or FPN (Brodmann areas (BAs) 6, 9, 10, 46, 47)(Assem et al., 2020; Deary et al., 2010; Duncan et al., 2020; Duncan et al., 2000; Gray et al., 2003; Jung & Haier, 2007). Across the whole brain search, the similar FCs (between MT+ and frontal cognitive cores) still showed significant correlations with BDT scores (Table supplement 1) (also shown in figure supplement 5a). Additionally, we identified certain parietal regions (BAs 7, 39, 40) with significant correlations between their connectivity to the left MT+ and the BDT scores (Table supplement 1) (also shown in figure supplement 5a). These significant connections between MT+ and FPN/MD system suggest that left MT+ is involved in the efficient information integration network, serving as a maker of intellectual hub.
![](https://prod--epp.elifesciences.org/iiif/2/97545%2Fv1%2Fcontent%2F585863v1_fig4.tif/full/max/0/default.jpg)
Significant FCs from connectivity-behavior analyses in a priori search space. The seed region is the left MT+. The significant FCs are obtained from a priori space (frontal cortex). (a) The significant FCs obtained from connectivity-BDT analysis. Single voxel threshold P < 0.005, adjacent size ≥ 23 (AlphaSim correcting, Methods). (b) The significant FCs obtained from connectivity-SI analysis. Single voxel threshold P < 0.005, adjacent size ≥ 22 (AlphaSim correcting, Methods). Positive correlations are shown in warm colors, while, negative correlations are shown in cold colors.
![](https://prod--epp.elifesciences.org/iiif/2/97545%2Fv1%2Fcontent%2F585863v1_tbl1.tif/full/max/0/default.jpg)
FC of voxels showing significant correlation with BDT scores across subjects in frontal cortex.
Then, we correlated MT+ - based global FCs with SI. Though spatial suppression during motion perception (quantified by SI) is considered to be the function of area MT+(Melnick et al., 2013), the top-down modulation from the frontal cortex can increase surround suppression(Liu et al., 2016). Our connectivity-SI analysis in the frontal regions (a priori search space) displayed 3 brain regions in which FCs strength significantly correlated with SI: right BA4/6, left BA6, and right BA46 (summarized in Table supplement 2, and shown in Figure 4b). Across the whole brain search, we identified total 7 brain regions in which FCs strength significantly correlated with SI, and 3 of these were in the frontal cortex. This is consistent with the results obtained by connectivity-SI analysis in a priori search space (frontal cortex) (Table supplement 3 and figure supplement 5b).
Local MT+ GABA acts on SI and BDT via global MT-frontal connectivity
To determine whether local neurotransmitter levels (such as GABA and Glu) in the MT+ region mediate the broader 3D visuo-spatial ability of BDT, which as a component of gF, is linked to the frontal cortex(Fangmeier et al., 2006), we correlated the significant FCs of MT – frontal in Figure 4a (also shown in Table 1) with the GABA and Glu levels in MT+ region. The results revealed that only two FCs significantly correlated with inhibitory GABA levels in MT+: 1) the FC of left MT+ - right BA 46 (significantly negative correlation, r = - 0.56, P = 0.02, n = 29, FDR correction, Figure 5a left); 2) the FC of left MT+ - right BA 6 (significantly positive correlation, r = 0.69, P = 0.002, n = 29, FDR correction, Figure 5b left) (also shown in Table 2). There were no significant correlations between these FCs and the excitatory Glu levels in MT+ (Table 2). Across the whole brain search, we obtained the same two MT+ - frontal FCs significantly correlating with both MT+ GABA levels and BDT (Table supplement 4), this is consistent with the results in a priori search space (frontal cortex) (Table 2). We then correlated the significant FCs in Figure 5b (also in Table supplement 2) with GABA and Glu concentrations in MT+ and found that almost all the correlations are significant except one (between the FC of left MT+ - right BA46 and the Glu levels in MT+) (Table supplement 5). Among the three FCs, the clusters of two FCs have substantial voxel overlap with the FCs we found by the connectivity-BDT analysis (Figure 5a, b). Across the whole brain search, there were total 7 brain regions in which FCs strength were significantly correlated with SI, all the 7 FCs significantly correlated the MT+ GABA levels, while, no FC had significant correlation with the MT+ Glu levels (Table supplement 6). Taken together, our results displayed that the overlap FCs from the analyses of connectivity - behavior (BDT and SI)-GABA are the MT+ - BA 46 and MT+ - BA 6 (Figure 5a, b). These results suggest that the FCs of MT+ - frontal regions (BA 46 and BA 6) coupling with local MT+ GABA underlie the neural basis for both the simple motion perception (quantified by SI) and the complex 3D visuo-spatial ability (quantified by BDT).
![](https://prod--epp.elifesciences.org/iiif/2/97545%2Fv1%2Fcontent%2F585863v1_fig5.tif/full/max/0/default.jpg)
Local MT+ GABA acts on SI and BDT via global MT-frontal connectivity. (a) Significant negative correlation between the FC of left MT+ - right DLPFC (BA46) and MT+ GABA (FDR correction). (b) Significant positive correlation between the FC of left MT+ - right (pre) motor cortex (BA4/6) and MT+ GABA (FDR correction). In (a) and (b), left: the significant FCs obtained from connectivity-BDT analysis; right: the significant FCs obtained from connectivity-SI analysis. (c) Significant pathways: MT+ GABA→ FC (left MT+ - right BA46, negative correlation) → SI (negative correlation) → BDT (positive correlation). This pathway can explain 34% of the variance in BDT. (d) Significant pathways: MT+ GABA → FC (left MT+ - right BA6, positive correlation) → BDT (positive correlation).The bolded lines represent the hypothesized mediation effect. The dotted lines represent alternative pathways. *: P < 0.05; **: P < 0.01; ***: P < 0.001.
![](https://prod--epp.elifesciences.org/iiif/2/97545%2Fv1%2Fcontent%2F585863v1_tbl2.tif/full/max/0/default.jpg)
Correlations between FC in Table 1 and GABA/Glu concentrations in MT+.
In order to fully investigate the potential roles of multivariable contributing to BDT scores, serial mediation analyses(Hayes, 2013) were applied to the MR and behavioral data. Following our hypothesis, the independent variable (X) is MT+ GABA, the dependent variable (Y) is BDT scores, the covariate is the age, and the mediators are FC (M1) and SI (M2). We used the overlap clusters from the analyses of connectivity-BDT-GABA and connectivity-SI-GABA to compute the FC of MT+ - BA46, 1 participant was excluded due to his age greater than mean + 2.5SD. The serial mediation model is shown in Figure 5c. GABA levels in MT+ significantly negatively correlated with the FC of MT+ - BA46 (β = - 0.32, P < 0.001), which in turn significantly negatively correlated with SI (β = - 0.19, P = 0.035), and consequently, significantly positively correlated with BDT (β = 38.5, P = 0.009). Critically, bootstrapped analyses revealed that our hypothesized indirect effect (i.e., MT+ GABA → FC of MT+ - BA46 → SI → BDT) was significant (β = 2.28, SE = 1.54, 95% CI = [0.03, 5.94]). The model accounted for 34% of the variance in BDT. However, when considering the MT-BA6 FC as the mediator M1, the serial model does not show a significant indirect effect. Consequently, we explored a mediation model, which revealed that the MT-BA6 FC totally mediates the relationship between GABA and BDS. (Figure 5d). For sensitivity purposes, we tested the alternative models, in which the order of the mediators was reversed. The pathway that MT+ GABA was predicted to be associated with SI, followed by the FC of MT+ - BA 46, and then BDT, did not yield the chained mediation effects on BDT (figure supplement 6).
To summarize (shown in Figure 6), the results from the serial mediation analyses are consistent with our hypothesis. That is, higher GABAergic inhibition in MT+ relates to stronger negative FC between MT+ and BA46, leading to enhanced ability for surround suppression (filtering out irrelevant information(Tadin, 2015), ultimately resulting in more efficient visual 3D processing of gF (the higher BDT scores).
![](https://prod--epp.elifesciences.org/iiif/2/97545%2Fv1%2Fcontent%2F585863v1_fig6.tif/full/max/0/default.jpg)
Sketch depicting the multi-level inhibitory mechanisms centered on MT+ GABA contributing to visuo-spatial intelligence. Inhibitory GABA in MT+ (a sensory cortex, shown in green circle), coupling with the conjunctive coding between MT+ and BA46 (cognitive control core, shown in purple circle), and mediated by motion surround suppression (shown in red circle), contributes to visuo-spatial intelligence (BDT, 3D domain, shown in red and white building blocks). In this sketch, the two-colored parallel arrows show the negative FC between MT+ and BA46, the colored arrows below the green circle display the inhibition mechanisms centered on MT+ GABA (2), filtered the irrelevant information in (1) and focused on the efficient visuo-spatial information (3). Black long arrows display the direction of information flow: from input information (1) to visuo-spatial intelligence (4).
Discussion
Here, we provide evidence that MT+ inhibitory mechanisms mediate processing in the visuo-spatial component (3D domain) of gF on multiple levels, that is, from molecular over brain connectivity to behavior. First, this study found that higher MT+ inhibitory GABA levels (but not excitatory Glu) relate to FC between MT+ and BA 46 that contribute to both SI and BDT. Our serial mediation analyses indicate that the inhibitory mechanisms linking to MT+, including GABA levels in MT+ (but not Glu), FCs of MT+ - BA46 coupling with MT+ inhibitory GABA (but not excitatory Glu), and behavior (SI indexing perceptual suppression in MT+) predict the inter-subject variance in the 3D gF task (BDT) (Figure 5c). Second, we demonstrate discrete GABAergic inhibition mechanisms in MT+ that mediate the strong FCs between MT+ - frontal regions (BA46 and BA6): significant negative correlation with the FC of MT+ - BA 46 (Figure 5a), whereas significant positive correlation with the FC of MT+ - BA 6 (Figure 5b). This indicates that different frontal regions, DLPFC (BA 46) and premotor cortex (BA6), contribute uniquely to gF through MT+ - based inhibitory mechanisms.
The goal of our research is to reveal that the inhibitory (not excitatory) mechanism in MT+ contributes to multi-level processing in 3D visuo-spatial ability (BDT). Monkey electrophysiological experiments revealed that selective attention gates the visual cortex, including area MT, effectively suppressing the irrelevant information(Everling et al., 2002; Treue & Maunsell, 1996). These findings aligns with the “neural efficiency” hypothesis of intelligence(Haier et al., 1988), which puts forward the human brain’s ability to suppress the repetition of information. Neural suppression is associated with the balance between excitation and inhibition (EIB), usually represented by covariation between Glutamate and GABA(Ozeki et al., 2009). Here, this study exploited the high spectral resolution afforded by ultrahigh field (7T) MRS to reliably resolve GABA measurement, adequately discriminate the glutamate and glutamine signals, and resolve the high accuracy Glu measurement(Ende, 2015).
This work implemented the MRS scanning in MT+ (3D visual domain) and V1 (2D visual domain) regions and found that MT+ inhibitory GABA (but not excitatory Glu) significantly correlated with BDT, i.e., the higher GABA levels in MT+ (rather than excitatory Glu) relate to higher visual 3D processing (BDT) (Figure 3c). We searched the global MT+ - based FCs with the connectivity-BDT analyses (in priori search space and whole brain search to valid), and then, correlated these significant FCs with the GABA and Glu concentrations in MT+, and found two FCs (MT+ - BA46 and MT+ - BA6) significantly correlating with MT+ inhibitory GABA (whereas no FC significantly correlated with MT+ excitatory Glu). Accordingly, our results emphasize the importance of MT+ inhibitory GABA (but not excitatory Glu) in processing the 3D visual-spatial intelligence (BDT). Crucially, a previous study has reported a strong link between motion spatial suppression (SI) and the performance of BDT(Melnick et al., 2013). Our recent human study(Song et al., 2021) and other study’s animal experiments(Ozeki et al., 2009; Sato et al., 2016) demonstrated that the conjoint action of inhibition (GABA) and excitation (Glu) underlies visual spatial suppression.
In this work, our novel data show the chained mediation effects from local MT+GABA to more global BDT: MT+ GABA → FC (MT+ and BA46) → SI → BDT. Thereby, our data indicate that inhibitory mechanisms in MT+, from the biochemistry over FC to the behavioral level, can predict the inter-subject variance in the 3D gF task (BDT) (Figure 5c).
Another interesting finding reveal that GABAergic inhibition in MT+ coupling with distinct brain connectivity mechanisms between BA 46-MT+ and BA6-MT+. A previous human fMRI experiment found that the positive and negative correlations between BDT and the activation of frontal regions appeared at different reasoning phases (validation or integration phases during reasoning)(Fangmeier et al., 2006). On the one hand, a monkey electrophysiological experiment reported the delayed modulation from PFC (especially in DLPFC (BA 46))) to area MT during a visual motion task(Zaksas & Pasternak, 2006). Computational models converged with empirical data of awake monkey experiments present slow temporal modulation from PFC to MT/MST(Donner et al., 2009; Siegel et al., 2015; Wang, 2002; Wimmer et al., 2015). On the other hand, human MEG studies(Donner et al., 2009; Wilming et al., 2020) reported that the gamma-band activity in the visual cortex (including area MT) exhibited high coherence with the activity in (pre-) motor regions (BA 4/6). These works suggest that the relation of long-range FC and local inhibition mechanism (MT+ GABA) might indicate that inhibition in MT+ contributes to efficient long-range integration and coordination in distant brain areas like the prefrontal and premotor cortex.
How does MT+ assemble into the MD system as an intellectual hub rather than a simple input module? The results in Figure 5a, c showed that the overlap brain regions from the analyses of connectivity-BDT-GABA/connectivity-SI-GABA are the MT+ - BA 46. This overlap couples with local visual suppression (SI) and consequently plays an important role in intelligence (BDT). The direction discrimination task in this work (the visual motion paradigm of center-surround antagonism) was previously considered a mainly local function of MT+(Melnick et al., 2013; Tadin, 2015; Tadin et al., 2003). However, our results with connectivity-SI analyses revealed that both local (FC within BA 18) and global brain connectivity (FC between MT+ and frontal regions) contribute to SI (Table supplement 3). In human psychophysical experiments(Melnick et al., 2013; Tadin et al., 2003) the brief stimulus duration (∼100 ms) in motion discrimination precludes most top-down attentional effects(Wang, 2002; Zaksas & Pasternak, 2006), while, attention, which predicted the performance of the motion discrimination task, was sustained throughout the stimulus intervals(Siegel et al., 2015). Furthermore, animal experiments have revealed that the local circuits in the visual cortex combing with the top-down modulation and intracortical horizontal connection mediate the visual-spatial suppression(Angelucci et al., 2002; Keller et al., 2020; Li et al., 2019; Zhang et al., 2014). Our results (shown in Figure 5a, b, right) present the intrinsic binding of local GABAergic inhibition in MT+, which suppressing redundancy of visual motion processing (SI), and the activity of brain connectivity between MT+ and frontal regions. These individually inherent traits may contribute to the individual difference in 3D visuo-spatial ability (Figure 5a, b, left). A candidate divisive normalization model(Carandini & Heeger, 2011; Reynolds & Heeger, 2009) can explain how such reverberation affects the process of suppressing the irrelevant information, from perception to intelligence(Melnick et al., 2013; Tadin, 2015). We summarize a framework (Figure 6) to indicate and visualize our findings.
Together, this study offers a comprehensive insight into how the information exchange and integration between the sensory cortex (MT+) and cognition core of BA 46, coupling with the MT+ GABA, can predict the performance of 3D visuo-spatial ability (BDT). Our results provide direct evidence that a sensory cortex (MT+), at biochemistry, connectivity, and behavior levels, can assemble into complex cognition as an intellectual hub.
Materials and Methods
Subjects
Thirty-six healthy subjects (18 female, mean age: 23.6 years±2.1, range: 20 to 29 years) participated in this study, they were recruited from Zhejiang University. All subjects had normal or corrected-to-normal vision. In addition, they reported no psychotropic medication use, no illicit drug use within the past month, no alcohol use within 3 days prior to scanning, and right-handed. This experiment was approved by the Ethics Review Committee of Zhejiang University and conducted in accordance with the Helsinki Declaration. All participants signed informed consent forms prior to the start of the study and were compensated for their time. All subjects participated in the motion spatial suppression psychophysical, resting-state fMRI and MRS (MT+ and V1 regions, in random sequence) experiments, but only part of the MRS data (31/36 in MT+ region and 28/36 in V1 region) survived quality control (see the part of MRS date processing). The sample size is determined by the statistic requirement (30 sample for Person correlation statistical analysis).
Motion surrounding suppression measurement
All stimuli were generated using Matlab (MathWorks, Natick, MA) with Psychophysics Toolbox(Brainard, 1997), and were shown on a linearized monitor (1920×1080 resolution, 100-Hz refresh rate, Cambridge Research System, Kent, UK). The viewing distance was 72 cm from the screen, with the head stabilized by a chinrest. Stimuli were drawn against a gray (56 cd per m-2) background.
A schematic of the stimuli and trial sequences is shown in our recent study(Song et al., 2021). The stimulus was a vertical drifting sinusoidal grating (contrast, 50%; spatial frequency, 1 cycle/°; speed, 4°/s) of either small (diameter of 2°) or large (diameter of 10°) size. The edge of the grating was blurred with a raised cosine function (width, 0.3°). A cross was presented in the center of the screen at the beginning of each trial for 500ms, and participants were instructed to fixate at the cross and to keep fixating at the cross throughout the trial. In each trial, a grating of either large or small size was randomly presented at the center of the screen. The grating drifted either leftward or rightward, and participants were asked to judge the perceived moving direction by a key press. Response time was not limited. The grating was ramped on and off with a Gaussian temporal envelope, and the grating duration was defined as 1 SD of the Gaussian function. The duration was adaptively adjusted in each trial, and duration thresholds were estimated by a staircase procedure. Thresholds for large and small gratings were obtained from a 160-trial block that contained four interleaved 3-down/1-up staircases. For each participant, we computed the correct rate for different stimulus durations separately for each stimulus size. These values were then fitted to a cumulative Gaussian function, and the duration threshold corresponding to the 75% correct point on the psychometric function was estimated for each stimulus size.
Stimulus demonstration and practice trials were presented before the first run. Auditory feedback was provided for each wrong response. To quantify the spatial suppression strength, we calculated the spatial suppression index (SI), defined as the difference of log10 thresholds for large versus small stimuli(Schallmo et al., 2018; Tadin et al., 2003):
Block design task measurement
The block design task was administered in accordance with the WAIS-IV manual(Wechsler, 2008). Specifically, participants were asked to rebuild the figural pattern within a specified time limit using a set of red and white blocks. The time limits were set as 30 s to 120 s according to different levels of difficulty. The patters were presented in ascending order of difficulty, and the test stopped if two consecutive patterns were not constructed in the allotted time. The score was determined by the accomplishment of the pattern and the time taken. A time bonus was awarded for rapid performance in the last six patterns. The score ranges between 0 and 66 points, with higher scores indicating better perceptual reasoning.
MR experimental procedure
MR experiments were performed in a 7T whole body MR system (Siemens Healthcare, Erlangen, Germany) with a Nova Medical 32 channel array head coil. Sessions included resting-state functional MRI, fMRI localizer scan, structural image scanning, and MRS scan. Resting-state scans were acquired with 1.5-mm isotropic resolution (transverse orientation, TR/TE = 2000/20.6 ms, 160 volumes, slice number = 90, flip angle = 70°, eyes closed). Structural images were acquired using a MP2RAGE sequence (TR/ TI1/ TI2 = 5000/901/3200 ms) with 0.7-mm isotropic resolution. MRS data were collected within two regions (MT+ and V1) for each subject, and we divided them into two sessions to avoid discomfort caused by long scanning. The order of MRS VOIs (MT+ and V1) in the two sessions was counterbalanced across participants. Interval between two sessions was used for block design and motion discrimination tasks. One session included fMRI localizer scan, structural image scanning, and MRS scan for the MT+ region; the other session included structural image scan, and MRS scan for the V1 region. Spectroscopy data were acquired using a 1H-MRS single-voxel short-TE STEAM (Stimulated Echo Acquisition Mode) sequence(Frahm et al., 1989) (TE/TM/TR = 6/32/7100ms) with 4096 sampling points, 4-kHz bandwidth, 16 averages, 8 repetitions, 20×20×20 mm3 VOI size, and VAPOR (variable power and optimized relaxation delays) water suppression(Tkáč et al., 1999). Prior to acquisition, first- and second-order shims were adjusted using FASTMAP (fast, automatic shimming technique by mapping along projections(Gruetter, 1993). Two non-suppressed water spectra were also acquired: one for phase and eddy current correction (only RF pulse, 4 averages) and another for metabolite quantification (VAPOR none, 4 averages). Voxels were positioned based on anatomical landmarks using a structural image scan collected in the same session, while avoiding contamination by CSF, bone, and fat. The MT+ VOIs were placed in the ventrolateral occipital lobe, which was based on anatomical landmarks(Dumoulin et al., 2000; Schallmo et al., 2018). We did not distinguish between the middle temporal (MT) and medial superior temporal (MST) areas in these MT+ VOIs(Huk et al., 2002). For 14 subjects, we also functionally identified MT+ as a check on the placement of the VOI. A protocol was used with a drifting grating (15% contrast) alternated with a static grating across blocks (10 s block duration, 160 TRs total). Using fMRI BOLD signals, these localizer data were processed online to identify the MT+ voxels in the lateral occipital cortex, which responded more strongly to moving vs. static gratings. In addition, we only used the left MT+ as the target region to scan, which was motivated by studies showing that left MT+ was more effective at causing perceptual effects(Tadin et al., 2011). For V1 region, the VOI was positioned on each subject’s calcarine sulcus on the left side(Tadin et al., 2011) based on anatomical landmarks(Boucard et al., 2007; Dumoulin et al., 2000)
MRS data processing
Spectroscopy data were preprocessed and quantified using magnetic resonance signal processing and analysis, https://www.cmrr.umn.edu/ downloads/mrspa/), which runs under MATLAB and invokes the interface of the LCModel (Version 6.3-1L) (Chen et al., 2019). First, we used the non-suppressed water spectra to perform eddy current correction and frequency/phase correction. Second, we checked the quality of each FID (16 averages) visually and removed those with obviously poor quality. Third, the absolute concentrations of each metabolite were quantitatively estimated via the Water-Scaling method. For partial-volume correction, the tissue water content was computed as follows(Ernst et al., 1993):
where fgm, fwm, and fcsf were the GM/WM/CSF volume fraction in MRS VOI and we used FAST (fMRI’s automated segmentation tool, part of the FSL toolbox)(Zhang et al., 2001) to segment the three tissue compartments from the T1-weighted structural brain images. For water T2 correction, we set water T2 as 47ms (Marjańska et al., 2012). Our concentrations were mM per kg wet weight. Furthermore, LCModel analysis was performed on all spectra within the chemical shift range of 0.2 to 4.0 ppm.
Poor spectral quality was established by a Cramer-Rao Lower Bound (CRLB) of more than 20%(Cavassila et al., 2001), and some data were excluded from further analysis. The details were described in our recently paper(Song et al., 2021).
Rs-fMRI data processing and analysis
Resting-state functional image was analyzed in the Data Processing and Analysis for Brain Imaging DPABI toolbox(Yan et al., 2016) based on SPM 12 (http://www.fil.ion.ucl.ac.uk/spm/). The preprocessing steps included discard of the first five volumes, slice timing, realignment to the 90th slice, coregistration of each subject’s T1-weighted anatomical and functional images, segmentation of the anatomical images into six types of tissues using DARTEL, linear detrend, regressing nuisance variables (including realignment Friston 24-parameter, global signal, white matter and CSF signal) (Friston et al., 1996), normalization to the standard Montreal Neurological Institute (MNI) space with the voxel size of 1.5× 1.5× 1.5 mm3 using DARTEL, spatial smoothing with a Gaussian kernel of 3 mm full-width-half-maximum (FWHM), and band-pass filtering with Standard frequency band (SFB, 0.01–0.1 Hz). Spherical ROI with a radius of 6mm was placed in left MT. The coordinate for left MT (−46, −72, −4, in MNI space) was obtained by our localizer fMRI experiment. We calculated the seed-to-voxel whole brain FC map for each subject. All the FC values were Fisher-Z-transformed.
We did a similar connectivity-behavior analysis to a previous study (Song et al., 2008). First, we computed the Pearson’s correlation coefficient between BDT scores and the FC values across subjects in a voxel-based way. Then, to evaluate the significance, we transformed the r-value into t-value
Statistical Analysis
PROCESS version 3.4, a toolbox in SPSS, was used to examine the mediation model. There are some prerequisites for mediation analysis: the independent variable should be a significant predictor of the mediator, and the mediator should be a significant predictor of the dependent variable.
SPSS 20 (IBM, USA) was used to conduct all the remaining statistical analysis in the study. We evaluated the correlation of variables (GABA, Glu, SI, BDT) using Pearson’s correlation analysis. Differences or correlations were considered statistically significant if P < 0.05. Significances with multiple comparisons were tested with false discovery rate (FDR) correction. The effect of age on intelligence was controlled for by using partial correlation in the correlation analysis and was taken as a covariate in the serial mediation model analysis.
Data availability
All data needed to evaluate the conclusions in the paper are present in the paper and/or the Supplementary information. Source data are provided with this paper and have been archived at GitHub. They could be downloaded with reasonable request.
Code availability
This paper does not report original code.
Acknowledgements
The authors thank Prof. Dost Ongur and Fei Du for guidance on the MRS data processing, and thank Zhejiang University 7T Brain Imaging Research Center. This work was supported by STI 2030—Major Projects (2021ZD0200401 to X.M.S., 2022ZD0206000 to R.B.), the National Natural Science Foundation of China Grants (61876222, 32000761, 82222032), Humanities and Social Sciences Ministry of Education (20YJC880095, 18YJA190001), the Key R&D Program of Zhejiang (2022C03096 to X.M.S.), the European Union’s Horizon 2020 Framework Program for Research and Innovation under the Specific Grant Agreement No. 785907 (Human Brain Project SGA2 to G.N.), and the MOE Frontier Science Center for Brain Science & Brain-Machine Integration, Zhejiang University.
Competing interests
The authors declare no competing interests.
![](https://prod--epp.elifesciences.org/iiif/2/97545%2Fv1%2Fcontent%2F585863v1_figs1.tif/full/max/0/default.jpg)
Relationships between GABA concentration in V1 and BDT/SI. (a) There is no significant relationship between BDT and GABA concentrations in V1 region. (b) There is also no obvious correlation between SI and GABA concentrations in V1 region. GABA concentration (Conc.) is absolute, with units of mmol per kg wet weight (Methods).
![](https://prod--epp.elifesciences.org/iiif/2/97545%2Fv1%2Fcontent%2F585863v1_figs2.tif/full/max/0/default.jpg)
Individual Glu MRS spectra from MT+ / V1 regions and relationships between BDT and Glu concentrations in MT+ / V1 regions. (a) Individual participant fitted Glutamate MRS spectra from the MT+ (top, n = 31) and V1 (bottom, n = 28) voxels from baseline measurement. The colors of the Glutamate spectra represent the individual differences of BDT. The color bar represents the scores of BDT. (b) There is no significant relationship between BDT and Glu concentrations in MT+ region. c There is also no significant correlation between BDT and Glu concentrations in V1 region. Glutamate concentration (Conc.) is absolute, with units of mmol per kg wet weight (Methods).
![](https://prod--epp.elifesciences.org/iiif/2/97545%2Fv1%2Fcontent%2F585863v1_figs3.tif/full/max/0/default.jpg)
Correlations between GABA and Glu concentrations in MT+ and V1 regions. (a) There is significant correlation between GABA and Glu concentrations in MT+ region. (b) The levels of GABA also significantly correlate with Glu in V1 region. GABA and Glu concentrations (Conc.) are absolute, with units of mmol per kg wet weight (Methods).
![](https://prod--epp.elifesciences.org/iiif/2/97545%2Fv1%2Fcontent%2F585863v1_figs4.tif/full/max/0/default.jpg)
Two linear correlations. (a) Significant positive correlation between suppression index and BDT scores (took out two outlines). (b) BDT also significantly correlates with GABA concentration in MT+ region (without two outliers, having the similar result shown in Figure 2d).
![](https://prod--epp.elifesciences.org/iiif/2/97545%2Fv1%2Fcontent%2F585863v1_figs5.tif/full/max/0/default.jpg)
Significant FCs searched from connectivity-behavior (BDT/SI) analyses in the whole brain. The seed region is the left MT+. The significant FCs are obtained from the entire brain search, single voxel threshold P < 0.01, adjacent size ≥ 37 voxels (AlphaSim correcting, Methods). Positive correlations are shown in warm colors, while, negative correlations are shown in cold colors. (a) the significant FCs obtained from connectivity-BDT analysis. (b) the significant FCs obtained from connectivity-SI analysis.
![](https://prod--epp.elifesciences.org/iiif/2/97545%2Fv1%2Fcontent%2F585863v1_figs6.tif/full/max/0/default.jpg)
The alternative serial mediation models from local MT+ GABA to global performance of BDT. The pathway that MT+ GABA was predicted to be associated with SI, followed by the FC of MT+ - BA 46, and then BDT, did not yield the chained mediation effects on BDT. *: P < 0.05; **: < 0.01; ***: P < 0.001.
![](https://prod--epp.elifesciences.org/iiif/2/97545%2Fv1%2Fcontent%2F585863v1_tbls1.tif/full/max/0/default.jpg)
FC of voxels showing significant correlation with BDT scores across subjects in whole brain.
![](https://prod--epp.elifesciences.org/iiif/2/97545%2Fv1%2Fcontent%2F585863v1_tbls2.tif/full/max/0/default.jpg)
FCs of voxels showing significant correlation with SI across subjects in frontal cortex.
![](https://prod--epp.elifesciences.org/iiif/2/97545%2Fv1%2Fcontent%2F585863v1_tbls3.tif/full/max/0/default.jpg)
FCs of voxels showing significant correlation with SI across subjects in whole Brain.
![](https://prod--epp.elifesciences.org/iiif/2/97545%2Fv1%2Fcontent%2F585863v1_tbls4.tif/full/max/0/default.jpg)
Correlations between FC in Table supplement 1 and GABA/Glu concentrations in MT+
![](https://prod--epp.elifesciences.org/iiif/2/97545%2Fv1%2Fcontent%2F585863v1_tbls5.tif/full/max/0/default.jpg)
Correlations between FCs in Table supplement 2 and GABA/Glu concentrations in MT+.
![](https://prod--epp.elifesciences.org/iiif/2/97545%2Fv1%2Fcontent%2F585863v1_tbls6.tif/full/max/0/default.jpg)
Correlations between FCs in Table supplement 3 and GABA/Glu concentrations in MT+.
References
- 1.Circuits for local and global signal integration in primary visual cortexJournal of Neuroscience 22:8633–8646
- 2.A domain-general cognitive core defined in multimodally parcellated human cortexCerebral Cortex 30:4361–4380
- 3.Network neuroscience theory of human intelligenceTrends in cognitive sciences 22:8–20
- 4.Sensitive period for a multimodal response in human visual motion area MT/MSTCurrent Biology 20:1900–1906
- 5.Structure and function of visual area MTAnnu Rev Neurosci 28:157–189https://doi.org/10.1146/annurev.neuro.26.041002.131052
- 6.Occipital proton magnetic resonance spectroscopy (1H-MRS) reveals normal metabolite concentrations in retinal visual field defectsPLoS One 2
- 7.The Psychophysics ToolboxSpat Vis 10:433–436
- 8.Normalization as a canonical neural computationNat Rev Neurosci 13:51–62https://doi.org/10.1038/nrn3136
- 9.Theory of fluid and crystallized intelligence: A critical experimentJournal of Educational Psychology 54:1–22
- 10.Cramér–Rao bounds: an evaluation tool for quantitationNMR in Biomedicine: An International Journal Devoted to the Development and Application of Magnetic Resonance In Vivo 14:278–283
- 11.Regional GABA concentrations modulate inter-network resting-state functional connectivityCerebral Cortex 29:1607–1618
- 12.Global connectivity of prefrontal cortex predicts cognitive control and intelligenceJ Neurosci 32:8988–8999https://doi.org/10.1523/JNEUROSCI.0536-12.2012
- 13.Distributed brain sites for the g-factor of intelligenceNeuroimage 31:1359–1365https://doi.org/10.1016/j.neuroimage.2006.01.006
- 14.The neuroscience of human intelligence differencesNature reviews neuroscience 11:201–211
- 15.Buildup of choice-predictive activity in human motor cortex during perceptual decision makingCurr Biol 19:1581–1585https://doi.org/10.1016/j.cub.2009.07.066
- 16.A new anatomical landmark for reliable identification of human area V5/MT: a quantitative analysis of sulcal patterningCereb Cortex 10:454–463https://doi.org/10.1093/cercor/10.5.454
- 17.Integrated Intelligence from Distributed Brain ActivityTrends Cogn Sci 24:838–852https://doi.org/10.1016/j.tics.2020.06.012
- 18.A neural basis for general intelligenceScience 289:457–460https://doi.org/10.1126/science.289.5478.457
- 19.Proton Magnetic Resonance Spectroscopy: Relevance of Glutamate and GABA to NeuropsychologyNeuropsychol Rev 25:315–325https://doi.org/10.1007/s11065-015-9295-8
- 20.Absolute quantitation of water and metabolites in the human brain. I. Compartments and water. Journal of magnetic resonance, Series B 102:1–8
- 21.Filtering of neural signals by focused attention in the monkey prefrontal cortexNature neuroscience 5:671–676
- 22.fMRI evidence for a three-stage model of deductive reasoningJournal of cognitive neuroscience 18:320–334
- 23.Localized high- resolution proton NMR spectroscopy using stimulated echoes: initial applications to human brain in vivoMagnetic resonance in medicine 9:79–93
- 24.Movement-related effects in fMRI time-seriesMagnetic resonance in medicine 35:346–355
- 25.Neural mechanisms of general fluid intelligenceNat Neurosci 6:316–322https://doi.org/10.1038/nn1014
- 26.Automatic, localized in vivo adjustment of all first-and second-order shim coilsMagnetic resonance in medicine 29:804–811
- 27.Cortical glucose metabolic rate correlates of abstract reasoning and attention studied with positron emission tomographyIntelligence 12:199–217
- 28.Introduction to Mediation, Moderation, and Conditional Process Analysis: A Regression-Based Approach
- 29.Retinotopy and functional subdivision of human areas MT and MSTJournal of Neuroscience 22:7195–7205
- 30.The Parieto-Frontal Integration Theory (P-FIT) of intelligence: converging neuroimaging evidenceBehav Brain Sci 30:135–154https://doi.org/10.1017/S0140525X07001185
- 31.Feedback generates a second receptive field in neurons of the visual cortexNature 582:545–549
- 32.Subdomains within orientation columns of primary visual cortexScience advances 5
- 33.From Molecular to Behavior: Higher Order Occipital Cortex in Major Depressive DisorderCereb Cortex 32:2129–2139https://doi.org/10.1093/cercor/bhab343
- 34.A neural basis for the spatial suppression of visual motion perceptionElife 5
- 35.Localized 1H NMR spectroscopy in different regions of human brain in vivo at 7 T: T2 relaxation times and concentrations of cerebral metabolitesNMR in Biomedicine 25:332–339
- 36.A strong interactive link between sensory discriminations and intelligenceCurr Biol 23:1013–1017https://doi.org/10.1016/j.cub.2013.04.053
- 37.Inhibitory stabilization of the cortical network underlies visual surround suppressionNeuron 62:578–592
- 38.The normalization model of attentionNeuron 61:168–185
- 39.The effect of visual experience on the development of functional architecture in hMT+Cerebral Cortex 17:2933–2939
- 40.An excitatory basis for divisive normalization in visual cortexNature neuroscience 19:568–570
- 41.Suppression and facilitation of human neural responsesElife 7https://doi.org/10.7554/eLife.30334
- 42.Cortical information flow during flexible sensorimotor decisionsScience 348:1352–1355https://doi.org/10.1126/science.aab0551
- 43.Brain spontaneous functional connectivity and intelligenceNeuroimage 41:1168–1176https://doi.org/10.1016/j.neuroimage.2008.02.036
- 44.Reduction of higher-order occipital GABA and impaired visual perception in acute major depressive disorderMol Psychiatry 26:6747–6755https://doi.org/10.1038/s41380-021-01090-5
- 45.General Intelligence” Objectively Determined and MeasuredAmerican Journal of Psychology 15:201–293
- 46.Suppressive mechanisms in visual motion processing: From perception to intelligenceVision Res 115:58–70https://doi.org/10.1016/j.visres.2015.08.005
- 47.Perceptual consequences of centre-surround antagonism in visual motion processingNature 424:312–315https://doi.org/10.1038/nature01800
- 48.Improved motion perception and impaired spatial suppression following disruption of cortical area MT/V5Journal of Neuroscience 31:1279–1283
- 49.In vivo 1H NMR spectroscopy of rat brain at 1 ms echo timeMagnetic Resonance in Medicine: An Official Journal of the International Society for Magnetic Resonance in Medicine 41:649–656
- 50.Attentional modulation of visual motion processing in cortical areas MT and MSTNature 382:539–541
- 51.Probabilistic decision making by slow reverberation in cortical circuitsNeuron 36:955–968https://doi.org/10.1016/s0896-6273(02)01092-9
- 52.Wechsler, D. (2008). Wechsler Memory Scale–Fourth Edition (WMS-IV) technical and interpretive manual.Wechsler Memory Scale–Fourth Edition (WMS-IV) technical and interpretive manual
- 53.Large-scale dynamics of perceptual decision information across human cortexNat Commun 11https://doi.org/10.1038/s41467-020-18826-6
- 54.Sensory integration dynamics in a hierarchical network explains choice probabilities in cortical area MTNat Commun 6https://doi.org/10.1038/ncomms7177
- 55.DPABI: data processing & analysis for (resting-state) brain imagingNeuroinformatics 14:339–351
- 56.Directional signals in the prefrontal cortex and in area MT during a working memory for visual motion taskJ Neurosci 26:11726–11742https://doi.org/10.1523/JNEUROSCI.3420-06.2006
- 57.Long-range and local circuits for top-down modulation of visual cortex processingScience 345:660–665
- 58.Segmentation of brain MR images through a hidden Markov random field model and the expectation-maximization algorithmIEEE transactions on medical imaging 20:45–57
Article and author information
Version history
- Sent for peer review:
- Preprint posted:
- Reviewed Preprint version 1:
- Reviewed Preprint version 2:
- Reviewed Preprint version 3:
Copyright
© 2024, Gao et al.
This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.
Metrics
- views
- 260
- downloads
- 18
- citations
- 0
Views, downloads and citations are aggregated across all versions of this paper published by eLife.