# Associations Between Individual Differences in Mathematical Competencies and Surface Anatomy of the Adult Brain

- Educational Neuroscience, Institute of Psychology, University of Graz, Graz, Austria

Previously conducted structural magnetic resonance imaging (MRI) studies on the neuroanatomical correlates of mathematical abilities and competencies have several methodological limitations. Besides small sample sizes, the majority of these studies have employed voxel-based morphometry (VBM)—a method that, although it is easy to implement, has some major drawbacks. Taking this into account, the current study is the first to investigate in a large sample of typically developed adults the associations between mathematical abilities and variations in brain surface structure by using surface-based morphometry (SBM). SBM is a method that also allows the investigation of brain morphometry by avoiding the pitfalls of VBM. Eighty-nine young adults were tested with a large battery of psychometric tests to measure mathematical competencies in four different areas: (1) simple arithmetic; (2) complex arithmetic; (3) higher-order mathematics; and (4) numerical intelligence. Also, we asked participants for their mathematics grades for their final school exams. Inside the MRI scanner, we collected high-resolution T1-weighted anatomical images from each subject. SBM analyses were performed with the computational anatomy toolbox (CAT12) and indices for cortical thickness, for cortical surface complexity, for gyrification, and sulcal depth were calculated. Further analyses revealed associations between: (1) the cortical surface complexity of the right superior temporal gyrus and numerical intelligence; (2) the depth of the right central sulcus and adults’ ability to solve complex arithmetic problems; and (3) the depth of the left parieto-occipital sulcus and adults’ higher-order mathematics competence. Interestingly, no relationships with previously reported brain regions were observed, thus, suggesting the importance of similar research to confirm the role of the brain regions found in this study.

## Introduction

Mathematical competencies are one of the key cognitive abilities in our modern societies, and they are crucial for our profession as well as social development (Parsons and Bynner, 2005). To date, an extensive number of neuroimaging studies have investigated the cognitive architecture of mathematical cognition. However, these studies have mainly focused on neurofunctional aspects of mathematical competencies, whereas literature on their neuroanatomical correlates is scarce. Especially, the surface characteristics of the gray matter (GM; i.e., cortical surface complexity, cortical thickness, gyrification, and sulcal depth) and their relationships with mathematical competencies have been widely ignored. The current study aimed to shed light on this aspect.

Numerous functional magnetic resonance imaging (fMRI) studies have demonstrated that numerical and mathematical cognition is associated with the activation of various frontal (e.g., ventral- and dorsolateral prefrontal cortex) and parietal brain regions (e.g., inferior and superior parietal cortex; for a review see Menon, 2014; for a meta-analysis see Arsalidou and Taylor, 2011). Also, it has been shown that the recruitment of these brain regions depends on the involved cognitive processes (e.g., quantity processing, arithmetic problem solving, etc.; Dehaene et al., 2003; Arsalidou and Taylor, 2011). For instance, while activation within the intraparietal sulcal (IPS) is thought to reflect quantity-based processes (Dehaene et al., 2003; Wilkey et al., 2017; Vogel et al., 2017a; or for a recent meta-analysis see Sokolowski et al., 2017), the controlled retrieval of a solution to a given arithmetic problem seems to rely on brain regions such as the inferior frontal gyrus (IFG), the supramarginal gyrus (SMG) and the angular gyrus (AG; Delazer et al., 2003, 2005; Ischebeck et al., 2006, 2009; Grabner et al., 2009; Klein et al., 2013, 2014; Menon, 2014; Peters and De Smedt, 2018; Heidekum et al., 2019). Furthermore, neuroimaging data indicate that activity in this network is also modulated by individual differences (e.g., Grabner et al., 2007; De Smedt and Gilmore, 2011; Berteletti et al., 2014). For instance, Grabner et al. (2007) compared two groups of young adults with different levels of mathematical competence who solved single-digit and multi-digit multiplication problems. The authors showed that individuals with higher mathematical competence more strongly activate the left AG while solving easy and more difficult multiplication problems. Grabner et al. (2007) concluded that mathematically more (compared to less) competent individuals rely more strongly on language-mediated processes (in particular fact retrieval) during arithmetic problem-solving.

In contrast to the large number of neurofunctional studies, only a few studies have used high-resolution structural MRI to examine how individual differences in brain morphology relate to variations in mathematical abilities and competencies (e.g., Isaacs et al., 2001; Aydin et al., 2007; Han et al., 2008, 2013; Rotzer et al., 2008; Rykhlevskaia et al., 2009; Lubin et al., 2013; Ranpura et al., 2013; Starke et al., 2013; Supekar et al., 2013; Cappelletti and Price, 2014; Evans et al., 2015; Price et al., 2016; Wilkey et al., 2018; Moreau et al., 2019). Many of those studies have focused on the comparison between dyscalculic individuals and controls (Han et al., 2008; Rotzer et al., 2008; Rykhlevskaia et al., 2009; Ranpura et al., 2013; Starke et al., 2013; Cappelletti and Price, 2014). For instance, it has been shown that dyscalculic adults have less GM volume in the right parietal cortex compared to controls (Cappelletti and Price, 2014) and that children with low mathematical skills relative to gender and age-matched controls show structural differences in the bilateral parietal lobes, right frontal lobe, and left occipital/parietal lobe (Han et al., 2008). These findings suggest an overlap of brain regions found in structural MRI studies and those reported in functional neuroimaging studies.

There are only a few studies (Supekar et al., 2013; Price et al., 2016; Wilkey et al., 2018) that investigated the neuroanatomical correlates of mathematics performance in typically developed individuals. However, most of these studies were conducted with children. For instance, Price et al. (2016) assessed the longitudinal and concurrent relations between GM volume and mathematical performance in a group of 50 children. GM volume maps derived from anatomical scans collected at the end of 1st grade and 2nd grade were related to performance on a standardized math test (i.e., Woodcock-Johnson III Tests of Achievement Calculation and Applied Problems). They found that left IPS GM volume at the end of 1st grade was positively related to math competence at the end of 2nd grade. Additionally, a positive association between GM volume in the same brain region at the end of 2nd grade and children’s concurrent math competence was observed. Wilkey et al. (2018) measured children’s performance in grade-level mathematics by a state-wide, school-based test of math achievement. They reported positive associations between differences in GM volume of the left and right hippocampal formations (including the hippocampus proper, entorhinal cortex, and subiculum) as well as the right IFG and children’s performance in grade-level mathematics.

Previous structural MRI studies on the neuroanatomical correlates of mathematical abilities and competencies have several methodological limitations. As highlighted above, the large majority of studies have involved children with mathematical learning disabilities (e.g., Han et al., 2013; Starke et al., 2013). Thus, they do not allow reliable inferences to typically developed populations because we cannot assume that neuroanatomical mechanisms that distinguish dyscalculic individuals from controls are the same that underlie individual differences in (typically developing) mathematical competence. Furthermore, we do not know whether the neuroanatomical mechanisms important during development are still playing a role in adults. Another methodological problem is the number of participants that were under investigation. Many findings are based on small sample sizes (Rotzer et al., 2008; Han et al., 2013), which is associated with various problems, such as low statistical power (Button et al., 2013) or an inflated false discovery rate (Colquhoun, 2014). Finally, most of these studies employed voxel-based morphometry (VBM) to investigate relationships between GM density and mathematical competencies. VBM is a technique that remains extremely popular because it is highly automated and therefore quick and easy to use. However, VBM also has been associated with important limitations in recent years. For example, recent studies have shown that differences in methodological choices like the registration algorithm (Peelle et al., 2012), or changes of user-specified parameters such as the smoothing kernel size (Henley et al., 2010) can lead to inconsistent VBM results. Furthermore, the interpretation of VBM findings is often difficult because the results may be driven by differences in cortical thickness, surface area (SA), cortical volume, and folding or any combination of these measures (Voets et al., 2008; Hutton et al., 2009).

A currently very promising alternative to VBM is surface-based morphometry (SBM), which also allows the investigation of brain morphometry. However, SBM uses brain surface meshes for spatial registration, which increases the accuracy of brain registration compared to mere volume-based registration (Desai et al., 2005). Further, it allows the calculation of additional measures of the neocortical surface structure, namely: cortical thickness (Dahnke et al., 2013), cortical surface complexity (quantification of the spatial frequency of gyrification and fissuration of the brain surface; Yotter et al., 2011), gyrification (Luders et al., 2006) and sulcal depth (Yotter et al., 2011). Due to these and further advantages, SBM has been frequently used in recent years, and a growing body of evidence suggests that cortical surface measures are particularly informative for individual differences in intellectual abilities (e.g., Geschwind and Rakic, 2013). For instance, longitudinal changes in cortical thickness have been related to variations in intelligence in children (Schnack et al., 2015) and a pronounced thickening of specific areas of the neocortex has been related to intellectual abilities in adulthood (Brans et al., 2010).

Only one study has used SBM in mathematical cognition so far. Moreau et al. (2019) investigated relationships between volumetric and surface characteristics of GM and dyslexia, dyscalculia and comorbid manifestation of both. They collected MRI data from four different groups of adults (i.e., dyslexic, dyscalculic, comorbid and control) and performed VBM as well as SBM analyses. By using Bayesian methods, Moreau et al. (2019) did not find any evidence for group differences in GM volume or any surface characteristics (i.e., cortical thickness, gyrification, sulcal depth, or cortical complexity). Therefore, the authors concluded that GM differences associated with these developmental disorders are not as reliable as previously suggested. However, these results were based on relatively small group samples (*N* = 12) and do not allow inferences to individual differences in typically developed populations.

Against this background, the current study aimed to investigate how structural variations in GM characteristics relate to individual differences in mathematical competencies in adults. Limitations of previous work were taken into account when designing the current study. First of all, to increase the power of the current study and the generalizability of its results, we collected data from a comparably large adult sample consisting of typically developed individuals. Thus, we attempted to identify neuroanatomical correlates of mathematical competencies in a population that has been underrepresented in previous studies. This will provide important insights into how the human brain mediates individual differences in mathematical cognition. Second, since previous studies showed that the neural correlates of mathematical competence depend on the mathematical demand (Price et al., 2016; Wilkey et al., 2018), we examined a broad range of mathematical competencies (e.g., arithmetic, higher-order mathematics) assessed using different instruments. Through the use of a large sample and different competence measures, we can investigate whether individual differences in mathematical competencies are related to a few brain regions of the classical model of mathematical cognition—as suggested by Dehaene et al. (2003)—or whether a more complex picture emerges. Finally, to avoid the above-described pitfalls of VBM we applied SBM to examine differences in brain surface structure (i.e., cortical thickness, cortical surface complexity, gyrification index, and sulcal depth).

## Materials and Methods

### Participants

Four-hundred and twenty-five German-speaking adults (female: *N* = 271; age: *M* = 23.13, *SD* = 5.58) participated in a larger behavioral test session, in which several psychometric measures were collected. Participants were recruited *via* social media, emails, and flyers. The majority of participants were enrolled as students at the University of Graz. All participants gave written informed consent before participation and received feedback regarding their intellectual abilities after testing as an incentive for taking part in the study.

From this pool, participants were recruited for two subsequent fMRI studies, in which, among other data, neurostructural images were collected. Both neuroimaging studies were carried out independently from each other. For the first neuroimaging study participants were randomly selected whereas in the second neuroimaging study participants were selected based on their arithmetic competencies (i.e., participants with lower and higher arithmetic competencies), which were measured in the behavioral test session. For both studies, all participants gave written informed consent before participation and were compensated with a minimum of € 20. The experimental procedure of both studies was approved by the ethics committee at the University of Graz, Austria.

For the current study, only participants from whom both data from the behavioral test session and neurostructural data were collected were included. In total 101 full data sets (Study 1: *N* = 46; Study 2: *N* = 55) were collected. However, due to partially missing data (*N* = 10) and other exclusion criteria (i.e., psychiatric disorders and left-handedness; *N* = 2) the final sample size comprises of 89 healthy young adults (Study 1: *N* = 36; Study 2: *N* = 53; age: *M* = 21.88; *SD* = 3.43; female: *N* = 57).

### Materials and Stimuli of the Behavioral Test Session

#### Berlin Intelligence Structure Test (BIS-T)

Participants’ numerical, verbal and figural intelligence was measured by the short version of the Berlin Intelligence Structure Test (Berliner Intelligenzstruktur-Test—BIS-T; Jäger et al., 1997). This test comprised 16 different tasks that either draw on the numerical (number of tests: 5), the verbal (number of tests: 6) or on the figural (number of tests: 5) component of intelligence. Additionally, each subtask can be assigned to four operational abilities, namely processing speed, memory, reasoning, and creativity. To assess numerical intelligence participants had to: (1) continue number series (reasoning); (2) cross out numbers in a matrix that were bigger by the factor of three in comparison to the preceding number (processing speed); (3) memorize pairs of digits (memory); (4) estimate the results of complex calculations (reasoning); and (5) to find different operands resulting in a given arithmetic solution (creativity). Tasks assessing verbal intelligence required participants to (6) fill in a missing letter resulting in a given word (processing speed); (7) list as many abilities a good salesman should not have (creativity); (8) complete analogies (reasoning); (9) decide whether a given statement represents a fact or an opinion (reasoning); (10) answer questions to a memorized text (memory); and (11) to cross out words that are meronyms of the preceding word (e.g., “year” preceded by the word “month”; processing speed). Finally, to assess figural intelligence participants had to (12) complete figural analogies (reasoning); (13) memorize marked buildings on a city map (memory); (14) design as many logos for a bike shop (creativity); (15) cross out “x” in an array of letters (processing speed); and (16) to complete figures (reasoning). For the subsequent analysis, we used sum scores for each intelligence component (potential maximum scores: numerical = 100; verbal = 128; figural = 181) that were based on raw scores (i.e., number of correct answers) of each subtask.

#### Arithmetic Fluency Task

Arithmetic competencies were assessed through a paper-pencil task developed by Vogel et al. (2017b) that is based on the French Kit Test (French et al., 1963). This test measures performance on simple and complex arithmetic problems. Participants were presented with 128 simple multiplications (i.e., consisting of two single-digit operands), followed by 64 simple additions (i.e., consisting of two single-digit operands), 128 simple subtractions (i.e., consisting of an operand <20 and a single-digit operand), 60 complex multiplications (i.e., consisting of a double-digit and a single-digit operand), 60 complex additions (i.e., consisting of three double-digit operands) and, finally, 60 complex subtractions (i.e., consisting of two double-digit operands). Participants were instructed to solve as many problems as possible within a given time. For solving simple arithmetic problems participants were given 90 s for each operation and 120 s for each operation type of the complex arithmetic problems. We then computed two scores indicating the numbers of correctly solved problems separately for each complexity level (i.e., simple vs. complex), where 320 was the potential maximum score for simple arithmetic problems and 180 for complex arithmetic problems.

#### Mathematics Test for Selection of Personnel (M-PA)

We assessed the subject’s performance in higher-order mathematics by using the short-version of the mathematics test for personnel selection (Mathematiktest für die Personalauswahl—M-PA; Jasper and Wagener, 2011). This test was developed to guarantee an optimal selection of suitable applicants based on their mathematical competencies. It is designed for teenagers and adults between 16 and 40 years who have at least a secondary school degree. The short version of the test comprises 31 complex mathematical problems presented as multiple-choice (MC) or open answer (OA) questions. Problems include fractions (3 OA), conversion of units (3 OA), exponentiation (7 OA), division with decimals (2 OA), algebra (1 MC), geometry (1 MC), roots (2 OA), and logarithm (7 OA). In total, participants had 15 min to solve the problems. For the subsequent analyses numbers of correctly solved problems were calculated (potential maximum score = 31).

#### Experimental Procedure

Behavioral data were collected in test sessions before the neuroimaging test sessions. Behavioral testing took place in laboratories of the Institute of Psychology at the University of Graz. One session took about 3 h, and participants were tested in small groups of a maximum of 12. After the greeting and general instruction participants were seated at a desk with a booklet, a computer screen and a keyboard on it. Desks were separated by partition walls, so that participants were able to work undisturbed. The booklet included a general instruction and all paper-pencil tests described in the “Materials and Methods” section. Additionally, the booklet included tests for measuring specific personality traits (i.e., the “NEO-Fünf-Faktoren Inventar—NEO-FFI” Borkenau and Ostendorf, 1993), the “German adaptation of the Abbreviated Math Anxiety Scale—AMAS-G” (Schillinger et al., 2018), the “Single-Item-Math Anxiety Scale—SIMA” (Núñez-Peña et al., 2014), a test measuring the attitude towards mathematics (Núñez-Peña et al., 2013, 2014), the “German Test Anxiety Inventory” (Prüfungsangstfragebogen, PAF; Hodapp et al., 2011), the “State-Trait Anxiety Inventory Trait scale—STAI-T” (Laux et al., 1981; Spielberger et al., 1983), and instructions for two computerized tasks (i.e., a single-digit multiplication task and an associative memory task), which participants also had to perform. The results of the latter tasks were not within the scope of the present study.

The order of tests was as following: BIS-T, M-PA, arithmetic fluency, computerized single-digit multiplication task, computerized associative memory task, NEO-FFI, PAF, attitudes towards mathematics, AMAS-G, SIMA, and STAI-T. Participants were instructed to work through the booklet page by page and to pause whenever they reached a page with a red stop sign. A trained experimenter took the time for all speeded tests (i.e., BIS-T: 45 min, M-PA: 15 min, arithmetic fluency: 10.5 min), and informed participants when they had to stop working on the respective test. At the end of the test booklet, participants were asked to fill out demographic questions according to their sex, age, field of study, and the mathematics grade in their final high school exam. In Austria, the grading system consists of a five-point scale, where 1 (“Sehr gut” – “very good”) is the best possible grade and 5 (“Nicht genügend” – “unsatisfactory”) is the lowest possible grade.

#### MRI Data Acquisition

Neuroimaging data were collected with a 3.0 T Siemens Skyra MRI scanner at the MRI Lab Graz. A 32-channel head coil and a Generalized Autocalibrating Partially Parallel Acquisitions (GRAPPA) sequence (TR = 1,950 ms, TE = 2.89 ms, 1 × 1 × 1 mm isotropic voxel resolution) was used to acquire the high-resolution T1-weighted anatomical images.

In addition to high-resolution T1-weighted anatomical images other neuroimaging data were collected in both studies: In study 1, functional and diffusion MRI (fMRI and dMRI) data were acquired. Participants had to perform three different tasks to measure the neural correlates of various cognitive processes. Results of the neurofunctional investigation are reported in Heidekum et al. (2019). In study 2 high-resolution T1-weighted anatomical images of participants were acquired within a training study that aimed to investigate the learning of arithmetic facts (i.e., multiplications). Participants underwent a 5-day multiplication fact training in which they intensively practiced a set of multiplications. Following the multiplication facts training, an MRI test session was conducted. In addition to the collection of anatomical data, functional imaging (i.e., task-based and resting-state activation data) as well as dMRI data were collected in this session. Since the behavioral and neurofunctional results of the training study are not within the scope of the present study they will be reported separately.

#### Analysis of Behavioral Data

Correlation coefficients were calculated to examine associations between all behavioral measurements. We calculated Spearman’s correlation coefficients for investigating associations with mathematics grade, whilst for the rest of the associations, Pearson’s correlation coefficients were calculated. *P*-values were Bonferroni corrected for 21 tested correlations (adjusted level of significance: *p*_{Bonf} < 0.00238). The normal distribution of each behavioral variable was verified by checking quantile-quantile plots and by calculating skewness and kurtosis. None of the variables significantly deviated from normality (i.e., *z*_{Skewness} > 1.96 or *z*_{Kurtosis} > 1.96).

#### Surface-Based Morphometry Analysis

SBM analyses were performed with the Computational Anatomy Toolbox (CAT12; Gaser and Dahnke, 2016), which is based on the Statistical Parametric Mapping (SPM, Welcome Department of Imaging Neuroscience, London, UK) software. In the first preprocessing step, T1-weighted anatomical images were normalized using “Diffeomorphic Anatomical Registration using Exponentiated Lie algebra” (DARTEL; Ashburner, 2007) and further segmented into GM, white matter (WM) and cerebrospinal fluid (CSF). In a next step, indices for cortical thickness (Dahnke et al., 2013), for cortical surface complexity (Yotter et al., 2011), for gyrification (Luders et al., 2006) and sulcal depth (Yotter et al., 2011) were calculated.

*Cortical Thickness*

Describes the distance between the inner surface (or boundary between GM and WM) and the outer surface (or boundary between GM and CSF; Dahnke et al., 2013). Cortical thickness, for instance, is an important biomarker for typical (O’Donnell et al., 2005) as well as atypical development (Thompson et al., 2004).

*Fractal Dimensionality Index (Cortical Surface Complexity)*

Characterizes the surface shape of the brain by quantifying the spatial frequency of gyrification and fissuration of the brain surface (Luders et al., 2004). More precisely, brain regions with a higher cortical surface complexity have a higher convolution level and therefore a larger surface area (Li et al., 2007; Yotter et al., 2011). Besides the investigation of age-related differences in brain structure (Madan and Kensinger, 2016), cortical complexity has also been successfully used to study differences in cognitive functions (King et al., 2010; Sandu et al., 2014).

*Gyrification Index*

Measures the regional curvature (i.e., convexity and concavity; Drury and Van Essen, 1997) of the brain and is defined as the ratio of the inner surface size to the outer surface size of an outer (usually convex) hull (Thompson et al., 1996; Luders et al., 2004, 2006). The gyrification index has typically been used to study group differences, such as between men and women (Luders et al., 2004) or between patients and controls (White et al., 2003; Shaw et al., 2012).

*Sulcal Depth*

Is based on the squared-transformed Euclidean distance between the central surface (average of the inner surface and the outer surface) and its convex hull (Lohmann, 1998; Yotter et al., 2011). Neuroanatomical studies using sulcal depth indices have shown, for instance, that it can account for differences in intelligence (Im et al., 2011; Yang et al., 2013).

Finally, full-width-at-half-maximum (FWHM) Gaussian kernels were used to smooth T1-weighted anatomical images. Following the matched-filter theorem, a 15.0 mm FWHM Gaussian kernel was applied for cortical thickness and a 20.0 mm one was used for folding measures (i.e., sulcal depth, gyrification, and cortical surface complexity).

In SPM we then performed second-level analyses to investigate whether structural variations in brain surface measures relate to individual differences in adults’ mathematical competencies. Multiple regression analyses for each mathematical competence measure (i.e., mathematics grade, simple arithmetic, complex arithmetic, higher-order mathematics, and numerical intelligence) were calculated separately. The multiple regression analyses were performed for each brain surface measurement (i.e., cortical thickness, cortical surface complexity, gyrification, and sulcal depth) and were controlled for the influence of verbal and figural intelligence as well as for age and sex (i.e., additional regressors of no interest). Statistical results of these whole-brain analyses are reported with family-wise error (FWE) corrected values at the cluster level (*p* < 0.05).

## Results

### Descriptive Analyses

In Table 1 means, standard deviations and ranges for all behavioral measurements are displayed.

### Correlation Analyses

Table 2 displays the correlations between all behavioral measurements. Almost all experimental measures that were thought to measure the same underlying component (i.e., mathematical cognition)—i.e., mathematics grade, simple arithmetic, complex arithmetic, higher-order mathematics and numerical intelligence—were significantly correlated with each other. The strongest association was found between performance on simple and complex arithmetic problems (*r* = 0.810, *p* < 0.001).

For the mathematics grade, we found negative correlations with performance on complex arithmetic problems (*r* = −0.358, *p* < 0.001) as well as knowledge in higher-order mathematics (*r* = −0.356, *p* < 0.001) and numerical intelligence (*r* = −0.346, *p* < 0.001). This means that individuals with a better mathematics grade (associated with a smaller number) had better scores on experimental tests measuring mathematical abilities.

Additionally, both simple and complex arithmetic was positively correlated with verbal and figural intelligence. However, there was no significant correlation of mathematics grade with verbal and figural intelligence and no significant correlation of higher-order mathematics with verbal and figural intelligence. Finally, as can be expected, all three intelligence components (i.e., numerical, verbal and figural intelligence) were positively correlated with each other.

### Structural Correlates of Mathematical Competencies

Multiple linear regression analyses were performed to investigate whether structural variations in brain surface measures (i.e., cortical thickness, cortical surface complexity, gyrification, and sulcal depth) relate to individual differences in adults’ mathematical abilities (i.e., mathematics grade, simple arithmetic, complex arithmetic, higher-order mathematics and numerical intelligence). The results are as followed:

Cortical thickness: no significant negative or positive associations were found.

Cortical surface complexity: analyses revealed that numerical intelligence was positively associated with the cortical surface complexity of the right superior temporal gyrus [MNI (x, y, z): 60, −9, −7; *p* < 0.05 FWE cluster corrected; *k* = 287; Figure 1A]. No significant negative associations were found.

**Figure 1. (A)** Significant positive association between numerical intelligence and cortical surface complexity of the right STG. **(B)** Significant positive association between performance on complex arithmetic problems and the depth of the right CS. **(C)** Significant positive association between higher-order mathematics competence and the depth of the left POS. Scatter plots showing partial correlations between **(A)** numerical intelligence and cortical surface complexity of the right STG, **(B)** performance on complex arithmetic problems and the depth of the right CS and **(C)** higher-order mathematics competence and the depth of the left POS. Colorbar: correlation coefficients (R-value). Abbreviations: STG, superior temporal gyrus; CSC, cortical surface complexity; CS, central sulcus; POS, parieto-occipital sulcus.

Gyrification: no significant negative or positive associations were found.

Sulcal depth: results showed a positive association between performance on complex arithmetic problems and the depth of the right central sulcus [MNI (x, y, z): 54, −9, 23; *p* < 0.05 FWE cluster corrected; *k* = 141; Figure 1B]. Additionally, multiple linear regression analysis revealed a positive relationship between adults’ competence in higher-order mathematics and the depth of the left parieto-occipital sulcus [MNI (x, y, z): −18, −70, 29; *p* < 0.05 FWE cluster corrected; *k* = 161; Figure 1C]. No significant negative associations were found.

## Discussion

The current study is the first to investigate the associations between variations in brain surface structure (i.e., cortical thickness, cortical surface complexity, gyrification index, and sulcal depth) and individual differences in a broad range of mathematical abilities (i.e., mathematics grade, performance on simple and complex arithmetic problems, higher-order mathematics and numerical intelligence) within a large sample of typically developed adults. Analyses revealed three brain regions that were associated with individual differences in mathematical abilities: (1) the cortical surface complexity of the right superior temporal gyrus was positively related to numerical intelligence; (2) the depth of the right central sulcus was positively associated with individual’s ability to solve complex arithmetic problems; and (3) the depth of the left parieto-occipital sulcus showed a positive relationship with the individual difference in performance on higher-order mathematics. Although SBM, as a tool to assess local brain morphology, does not allow any direct assumption of local brain function, this study supports the idea that these brain regions play a role in mathematical thinking. Surprisingly, we did not find any associations between adults’ mathematical abilities and the neuroanatomical differences of previously reported brain regions, such as the IPS (Li et al., 2013; Lubin et al., 2013; Price et al., 2016). Conversely, all three brain regions that were observed in the present study are not typically associated with mathematics-related cognitive abilities. Despite these unexpected findings, the following sections will be an attempt to link functional properties to the observed cortical regions. Since the present study can only speculate on the functional role of these regions, further research (applying SBM for investigating brain-behavior relationships in the field of mathematical cognition) is needed.

First, our analyses revealed a positive relationship between the cortical surface complexity of the right superior temporal gyrus [MNI (x, y, z): 60, −9, −7] and individual differences in numerical intelligence. Surface complexity was defined by the fractal dimensionality index, which quantifies the spatial frequency of gyrification and fissuration of the brain surface (Luders et al., 2004). Thus, the current finding indicates that in individuals with higher numerical intelligence the right superior temporal gyrus shows a higher convolution level and therefore a larger surface area. Cortical complexity has been successfully used to investigate age-related differences in brain structure (Madan and Kensinger, 2016) and, more importantly, to study differences in cognitive functions (King et al., 2010; Sandu et al., 2014). Unfortunately, the scarce neuroanatomical studies on the correlates of the right superior temporal gyrus are typically based on comparisons between various groups of patients and controls (e.g., schizophrenia: Honea et al., 2008; Nenadic et al., 2010; depression: Mak et al., 2009; Peng et al., 2011). Due to the methodological differences of these studies, they do not allow an interpretation of the current findings. Similarly, the results of neurofunctional studies are very heterogeneous. The right superior temporal gyrus has been implicated in several cognitive processes, such as spatial awareness (Karnath, 2001), emotion processing (Narumoto et al., 2001) or the detection of biological motion (Akiyama et al., 2006). Interestingly, previous work (Jung-Beeman et al., 2004) has also linked activity in the right superior temporal gyrus to creative problem solving, as it might be responsible for binding and accessing various types of available conceptual representations (Shen et al., 2017); cognitive processes that might also be important for problems in a numerical intelligence test. Besides these findings on domain-general cognitive functions, some studies have found activity in this region to be related to mathematical cognition (Fehr et al., 2008; Simos et al., 2008; Andres et al., 2011, 2012; Kucian et al., 2011; Gullick et al., 2012). Since the numerical intelligence test used in the current study required the correct application of arithmetic operations (i.e., continuing number series; crossing out numbers in a matrix that were bigger by the factor of three in comparison to the preceding number; finding different operands resulting in a given arithmetic solution), it is interesting that the majority of these studies (Fehr et al., 2008; Simos et al., 2008; Andres et al., 2011, 2012) has linked activation within this region to arithmetic problem-solving. For instance, in an fMRI block design, Andres et al. (2011) instructed a group of young males to multiply single Arabic digits (between 3 and 9) with, or subtract them from, a predefined digit. By contrasting multiplication problems with subtractions, the authors observed increased activity in right superior temporal areas (among other brain regions), which they related to the storing of semantic knowledge of arithmetic problems. Additionally, Fehr et al. (2008) investigated fMRI activation patterns during mental addition, subtraction, multiplication, and division in young adults. In contrast to Andres et al. (2011), greater right-hemispheric superior temporal activation was found by contrasting complex with simple division, which might have reflected the usage of procedural strategies for problem-solving. Accordingly, even though these findings are based on neurofunctional studies and they relate activity within the right superior temporal gyrus to different arithmetic strategies, they are in line with the result of the current study. Nevertheless, additional research is needed to confirm and to further specify the observed brain-behavior relationship.

Second, our results also revealed a positive correlation between the depth of the central sulcus [MNI (x, y, z): 54, −9, 23] and adults’ ability to solve complex arithmetic problems. Previously, sulcal depth indices (e.g., sulcal fundi, lines or pits) have been successfully implemented to study the neuroanatomical correlates of various cognitive mechanisms. For instance, variations in the properties of sulci can account for individual differences in intelligence (Im et al., 2011; Yang et al., 2013), and abnormalities in specific sulci, such as the IPS or the superior temporal sulcus, have been observed in different groups of atypically developed individuals (e.g., dyslexia: Steinbrink et al., 2008; Richlan et al., 2013; dyscalculia: Rotzer et al., 2008; Price et al., 2016). Morphological changes of the central sulcus, on the other hand, have been associated with neurodegenerative diseases, such as multiple sclerosis (Pagani et al., 2005; Ceccarelli et al., 2008). However, the central sulcus also acts as an extension of its adjacent areas, namely the pre- and postcentral gyrus. Whereas the precentral gyrus is involved in eye movements (Corbetta et al., 1998; Anderson et al., 2007), postcentral activations have been linked to grasping (Simon et al., 2002)—both processes are commonly involved in cognitive tasks that include visually presented stimuli. For that reason, it is not surprising that these brain areas have also been found to play a role in numerical and calculation tasks (for a meta-analysis see Arsalidou and Taylor, 2011). For instance, Kaufmann et al. (2008) showed that activation patterns within these regions were related to non-symbolic numerical as well as spatial processing. And conversely, Kesler et al. (2006) have linked pre- and postcentral activations to the use of arithmetic strategies such as subvocalization and finger counting. In the present study, we observed a relationship between performance on complex arithmetic problems and the depth of the central sulcus. As the central sulcus acts as an extension of the pre- and postcentral gyrus, we, therefore, conclude that it might affect its adjacent mathematic-related regions.

Finally, the results of our SBM analyses revealed a positive association between the depth of the parieto-occipital sulcus [MNI (x, y, z): −18, −70, 29] and adult’s higher-order mathematical competence. The parieto-occipital sulcus lies between the posterior parietal and occipital cortices. This brain region is a crucial node of the dorsal visual stream and is, therefore, involved in visuospatial processing that support spatial navigation and goal-directed actions (Caminiti et al., 1999; Hutchison et al., 2015; Tosoni et al., 2015; Richter et al., 2019). However, the parieto-occipital sulcus also separates the cuneus from the precuneus and lies in the direct vicinity of the IPS. Whereas the cuneus is associated with visual information processing (Vanni et al., 2001), the latter brain regions (i.e., precuneus and IPS) were found to contribute to higher neurocognitive functions, such as attention (Corbetta and Shulman, 2002; Rosen et al., 2015), working memory (Barton and Brewer, 2013; Bray et al., 2015) or memory retrieval (Hebscher et al., 2019; Heidekum et al., 2019). Moreover, the precuneus and the IPS are often found to be involved in mathematical cognition (Dehaene et al., 2003; Arsalidou and Taylor, 2011). In particular, previous studies in the field of mathematical cognition have related brain activity in the precuneus to mental calculation (i.e., addition, multiplication, and subtraction; for a recent meta-analysis see Arsalidou and Taylor, 2011) and IPS activation to the representation and manipulation of numerical magnitude (Dehaene et al., 2003; or for a recent meta-analysis see Sokolowski et al., 2017). For instance, greater activation in the precuneus is often observed when arithmetic problems are solved by procedural strategies, rather than by memory retrieval (Delazer et al., 2005; Fehr et al., 2008). And, IPS activity is typically found whenever participants explicitly or implicitly compare the magnitude of numerals or dot-arrays (Pinel et al., 2001; Ansari et al., 2005; Vogel et al., 2017a; Wilkey et al., 2017). Interestingly, the functional role of the IPS within mathematics has also been confirmed by neurostructural studies (Isaacs et al., 2001; Rotzer et al., 2008; Rykhlevskaia et al., 2009; Ranpura et al., 2013) showing reduced GM in children with dyscalculia. In general, these studies highlight the importance of parietal brain regions in mathematical thinking, which is in line with the result of the current study showing an association between the parieto-occipital sulcus and higher-order mathematics.

Taken together, the current study highlights the potential of SBM as an alternative to classical VBM for investigating brain-behavior relationships in the domain of mathematical cognition. Associations were found between: (1) the cortical surface complexity of the right superior temporal gyrus and numerical intelligence; (2) the depth of the right central sulcus and adults’ ability to solve complex arithmetic problems; and (3) the depth of the left parieto-occipital sulcus and adults’ higher-order mathematics competence. Although two out of three brain regions (i.e., central sulcus and parieto-occipital sulcus) lie in intermediate vicinity of cortical structures typically associated with numerical and mathematical cognition, none of these brain regions were observed in previous studies. This discrepancy could be due to the methodological approach of the current study: In contrast to previous work, the present study studied linear brain-behavior relationships in a large sample of typically developed adults, whereas previous studies are mainly based on group comparisons between children with and without mathematical learning disabilities. For that reason, our results suggest that in typically developed individuals different neuroanatomical structures might become important when it comes to the processing of numerical and mathematical information. Finally, the current study shows that the depth of sulci could be a good index to investigate brain-behavior relationships. Nevertheless, future studies implementing SBM in different groups of individuals are needed.

## Data Availability Statement

The datasets generated for this study are available on request to the corresponding author.

## Ethics Statement

The studies involving human participants were reviewed and approved by the ethics commitee of the University of Graz. The patients/participants provided their written informed consent to participate in this study.

## Author Contributions

AH: conceptualization, data curation, formal analysis, writing—original draft. SV: conceptualization, formal analysis, funding acquisition, supervision, writing—original draft, writing—review and editing. RG: conceptualization, funding acquisition, project administration, supervision, writing—original draft, writing—review and editing.

## Funding

This study was supported by a joint grant from the Research Foundation Flanders (Fonds Wetenschappelijk Onderzoek; FWO, G.0027.16) and the Austrian Science Fund (FWF, I 2425-G18).

## Conflict of Interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

## Acknowledgments

We thank Antonia Reuss, Dennis Wambacher and Thomas Zussner for assistance with data collection. Additionally, we thank Jochen Mosbacher for advice on the data analysis.

## References

Akiyama, T., Kato, M., Muramatsu, T., Saito, F., Nakachi, R., and Kashima, H. (2006). A deficit in discriminating gaze direction in a case with right superior temporal gyrus lesion. *Neuropsychologia* 44, 161–170. doi: 10.1016/j.neuropsychologia.2005.05.018

Anderson, E. J., Mannan, S. K., Husain, M., Rees, G., Sumner, P., Mort, D. J., et al. (2007). Involvement of prefrontal cortex in visual search. *Exp. Brain Res.* 180, 289–302. doi: 10.1007/s00221-007-0860-0

Andres, M., Michaux, N., and Pesenti, M. (2012). Common substrate for mental arithmetic and finger representation in the parietal cortex. *NeuroImage* 62, 1520–1528. doi: 10.1016/j.neuroimage.2012.05.047

Andres, M., Pelgrims, B., Michaux, N., Olivier, E., and Pesenti, M. (2011). Role of distinct parietal areas in arithmetic: an fMRI-guided TMS study. *NeuroImage* 54, 3048–3056. doi: 10.1016/j.neuroimage.2010.11.009

Ansari, D., Garcia, N., Lucas, E., Hamon, K., and Dhitalm, B. (2005). Neural correlates of symbolic number processing in children and adults. *Neuroreport* 16, 1769–1773. doi: 10.1097/01.wnr.0000183905.23396.f1

Arsalidou, M., and Taylor, M. J. (2011). Is 2+2=4? Meta-analyses of brain areas needed for numbers and calculations. *NeuroImage* 54, 2382–2393. doi: 10.1016/j.neuroimage.2010.10.009

Ashburner, J. (2007). A fast diffeomorphic image registration algorithm. *NeuroImage* 38, 95–113. doi: 10.1016/j.neuroimage.2007.07.007

Aydin, K., Ucar, A., Oguz, K. K., Okur, O. O., Agayev, A., Unal, Z., et al. (2007). Increased gray matter density in the parietal cortex of mathematicians: a voxel-based morphometry study. *Am. J. Neuroradiol.* 28, 1859–1864. doi: 10.3174/ajnr.A0696

Barton, B., and Brewer, A. A. (2013). Visual working memory in human cortex. *Psychology* 4, 655–662. doi: 10.4236/psych.2013.48093

Berteletti, I., Prado, J. Ô., and Booth, J. R. (2014). Children with mathematical learning disability fail in recruiting verbal and numerical brain regions when solving simple multiplication problems. *Cortex* 57, 143–155. doi: 10.1016/j.cortex.2014.04.001

Brans, R. G. H., Kahn, S., Schnack, H. G., Van Baal, G. C. M., Posthuma, D., Van Haren, N. E. M., et al. (2010). Brain plasticity and intellectual ability are influenced by shared genes. *J. Neurosci.* 30, 5519–5524. doi: 10.1523/JNEUROSCI.5841-09.2010

Bray, S., Almas, R., Arnold, A. E. G. F., Iaria, G., and MacQueen, G. (2015). Intraparietal sulcus activity and functional connectivity supporting spatial working memory manipulation. *Cereb. Cortex* 25, 1252–1264. doi: 10.1093/cercor/bht320

Button, K. S., Ioannidis, J. P. A., Mokrysz, C., Nosek, B. A., Flint, J., Robinson, E. S. J., et al. (2013). Power failure: why small sample size undermines the reliability of neuroscience. *Nat. Rev. Neurosci.* 14, 365–376. doi: 10.1038/nrn3475

Caminiti, R., Genovesio, A., Marconi, B., Mayer, A. B., Onorati, P., Ferraina, S., et al. (1999). Early coding of reaching: frontal and parietal association connections of parieto-occipital cortex. *Eur. J. Neurosci.* 11, 3339–3345. doi: 10.1046/j.1460-9568.1999.00801.x

Cappelletti, M., and Price, C. J. (2014). Residual number processing in dyscalculia. *Neuroimage Clin.* 4, 18–28. doi: 10.1016/j.nicl.2013.10.004

Ceccarelli, A., Rocca, M. A., Pagani, E., Colombo, B., Martinelli, V., Comi, G., et al. (2008). A voxel-based morphometry study of grey matter loss in MS patients with different clinical phenotypes. *NeuroImage* 42, 315–322. doi: 10.1016/j.neuroimage.2008.04.173

Colquhoun, D. (2014). An investigation of the false discovery rate and the misinterpretation of p-values. *R. Soc. Open Sci.* 1:140216. doi: 10.1098/rsos.140216

Corbetta, M., Akbudak, E., Conturo, T. E., Snyder, A. Z., Ollinger, J. M., Drury, H. A., et al. (1998). A common network of functional areas for attention and eye movements. *Neuron* 21, 761–773. doi: 10.1016/s0896-6273(00)80593-0

Corbetta, M., and Shulman, G. L. (2002). Control of goal-directed and stimulus-driven attention in the brain. *Nat. Rev. Neurosci.* 3, 201–215. doi: 10.1038/nrn755

Dahnke, R., Yotter, R. A., and Gaser, C. (2013). Cortical thickness and central surface estimation. *NeuroImage* 65, 336–348. doi: 10.1016/j.neuroimage.2012.09.050

De Smedt, B., and Gilmore, C. K. (2011). Defective number module or impaired access? Numerical magnitude processing in first graders with mathematical difficulties. *J. Exp. Child Psychol.* 108, 278–292. doi: 10.1016/j.jecp.2010.09.003

Dehaene, S., Piazza, M., Pinel, P., and Cohen, L. (2003). Three parietal circuits for number processing. *Cogn. Neuropsychol.* 20, 487–506. doi: 10.1080/02643290244000239

Delazer, M., Domahs, F., Bartha, L., Brenneis, C., Lochy, A., Trieb, T., et al. (2003). Learning complex arithmetic—an fMRI study. *Cogn. Brain Res.* 18, 76–88. doi: 10.1016/j.cogbrainres.2003.09.005

Delazer, M., Ischebeck, A., Domahs, F., Zamarian, L., Koppelstaetter, F., Siedentopf, C. M., et al. (2005). Learning by strategies and learning by drill—evidence from an fMRI study. *NeuroImage* 25, 838–849. doi: 10.1016/j.neuroimage.2004.12.009

Desai, R., Liebenthal, E., Possing, E. T., Waldron, E., and Binder, J. R. (2005). Volumetric vs. surface-based alignment for localization of auditory cortex activation. *NeuroImage* 26, 1019–1029. doi: 10.1016/j.neuroimage.2005.03.024

Drury, H. A., and Van Essen, D. C. (1997). Functional specializations in human cerebral cortex analyzed using the visible man surface-based atlas. *Hum. Brain Mapp.* 5, 233–237. doi: 10.1002/(sici)1097-0193(1997)5:4<233::aid-hbm5>3.0.co;2-4

Evans, T. M., Kochalka, J., Ngoon, T. J., Wu, S. S., Qin, S., Battista, C., et al. (2015). Brain structural integrity and intrinsic functional connectivity forecast 6 year longitudinal growth in children’s numerical abilities. *J. Neurosci.* 35, 11743–11750. doi: 10.1523/JNEUROSCI.0216-15.2015

Fehr, T., Code, C., and Herrmann, M. (2008). Auditory task presentation reveals predominantly right hemispheric fMRI activation patterns during mental calculation. *Neurosci. Lett.* 431, 39–44. doi: 10.1016/j.neulet.2007.11.016

French, J. W., Ekstrom, R. B., and Price, L. A. (1963). *Manual for Kit of Reference Tests for Cognitive Factors (Revised 1963).* Princeton, NJ: Educational Testing Service.

Gaser, C., and Dahnke, R. (2016). CAT – A computational anatomy toolbox for the analysis of structural MRI data. *OHBM Conference 2016*.

Geschwind, D. H., and Rakic, P. (2013). Cortical evolution: judge the brain by its cover. *Neuron* 80, 633–647. doi: 10.1016/j.neuron.2013.10.045

Grabner, R. H., Ansari, D., Reishofer, G., Stern, E., Ebner, F., and Neuper, C. (2007). Individual differences in mathematical competence predict parietal brain activation during mental calculation. *NeuroImage* 38, 346–356. doi: 10.1016/j.neuroimage.2007.07.041

Grabner, R. H., Ischebeck, A., Reishofer, G., Koschutnig, K., Delazer, M., Ebner, F., et al. (2009). Fact learning in complex arithmetic and figural-spatial tasks: the role of the angular gyrus and its relation to mathematical competence. *Hum. Brain Mapp.* 30, 2936–2952. doi: 10.1002/hbm.20720

Gullick, M. M., Wolford, G., and Temple, E. (2012). Understanding less than nothing: neural distance effects for negative numbers. *NeuroImage* 62, 542–554. doi: 10.1016/j.neuroimage.2012.04.058

Han, Z., Davis, N., Fuchs, L., Anderson, A. W., Gore, J. C., and Dawant, B. M. (2013). Relation between brain architecture and mathematical ability in children: a DBM study. *Magn. Reson. Imaging* 31, 1645–1656. doi: 10.1016/j.mri.2013.08.008

Han, Z., Fuchs, L., Davis, N., Cannistraci, C. J., Anderson, A. W., Gore, J. C., et al. (2008). “Analysis of anatomic variability in children with low mathematical skills,” in *Proceedings of the Medical Imaging 2008: Physiology, Function and Structure from Medical Images*, San Diego, CA, 69160S.

Hebscher, M., Meltzer, J. A., and Gilboa, A. (2019). A causal role for the precuneus in network-wide theta and gamma oscillatory activity during complex memory retrieval. *Elife* 8:e43114. doi: 10.7554/eLife.43114

Heidekum, A. E., Grabner, R. H., De Smedt, B., De Visscher, A., and Vogel, S. E. (2019). Interference during the retrieval of arithmetic and lexico-semantic knowledge modulates similar brain regions: evidence from functional magnetic resonance imaging (fMRI). *Cortex* 120, 375–393. doi: 10.1016/j.cortex.2019.06.007

Henley, S. M. D., Ridgway, G. R., Scahill, R. I., Klöppel, S., Tabrizi, S. J., Fox, N. C., et al. (2010). Pitfalls in the use of voxel-based morphometry as a biomarker: examples from Huntington disease. *AJNR Am. J. Neuroradiol.* 31, 711–719. doi: 10.3174/ajnr.A1939

Hodapp, V., Rohrmann, S., and Ringeisen, T. (2011). *PAF—Prüfungsangstfragebogen.* Göttingen: Hogrefe.

Honea, R. A., Meyer-Lindenberg, A., Hobbs, K. B., Pezawas, L., Mattay, V. S., Egan, M. F., et al. (2008). Is gray matter volume an intermediate phenotype for schizophrenia? A voxel-based morphometry study of patients with schizophrenia and their healthy siblings. *Biol. Psychiatry* 63, 465–474. doi: 10.1016/j.biopsych.2007.05.027

Hutchison, R. M., Culham, J. C., Flanagan, J. R., Everling, S., and Gallivan, J. P. (2015). Functional subdivisions of medial parieto-occipital cortex in humans and nonhuman primates using resting-state fMRI. *NeuroImage* 116, 10–29. doi: 10.1016/j.neuroimage.2015.04.068

Hutton, C., Draganski, B., Ashburner, J., and Weiskopf, N. (2009). A comparison between voxel-based cortical thickness and voxel-based morphometry in normal aging. *NeuroImage* 48, 371–380. doi: 10.1016/j.neuroimage.2009.06.043

Im, K., Choi, Y. Y., Yang, J. J., Lee, K. H., Kim, S. I., Grant, P. E., et al. (2011). The relationship between the presence of sulcal pits and intelligence in human brains. *NeuroImage* 55, 1490–1496. doi: 10.1016/j.neuroimage.2010.12.080

Isaacs, E. B., Edmonds, C. J., Lucas, A., and Gadian, D. G. (2001). Calculation difficulties in children of very low birthweight: a neural correlate. *Brain* 124, 1701–1707. doi: 10.1093/brain/124.9.1701

Ischebeck, A., Zamarian, L., Schocke, M., and Delazer, M. (2009). Flexible transfer of knowledge in mental arithmetic—an fMRI study. *NeuroImage* 44, 1103–1112. doi: 10.1016/j.neuroimage.2008.10.025

Ischebeck, A., Zamarian, L., Siedentopf, C., Koppelstätter, F., Benke, T., Felber, S., et al. (2006). How specifically do we learn? Imaging the learning of multiplication and subtraction. *NeuroImage* 30, 1365–1375. doi: 10.1016/j.neuroimage.2005.11.016

Jäger, A. O., Süß, H.-M., and Beauducel, A. (1997). *Berliner Intelligenzstruktur-Test [Berlin Intelligence Structure Test].* Göttingen: Hogrefe.

Jasper, F., and Wagener, D. (2011). *Mathematiktest für die Personalauswahl [Mathematics Test for Selection of Personnel].* Göttingen: Hogrefe.

Jung-Beeman, M., Bowden, E. M., Haberman, J., Frymiare, J. L., Arambel-Liu, S., Greenblatt, R., et al. (2004). Neural activity when people solve verbal problems with insight. *PLoS Biol.* 2:E97. doi: 10.1371/journal.pbio.0020097

Karnath, H. O. (2001). New insights into the functions of the superior temporal cortex. *Nat. Rev. Neurosci.* 2, 568–576. doi: 10.1038/35086057

Kaufmann, L., Vogel, S. E., Wood, G., Kremser, C., Schocke, M., Zimmerhackl, L. B., et al. (2008). A developmental fMRI study of nonsymbolic numerical and spatial processing. *Cortex* 44, 376–385. doi: 10.1016/j.cortex.2007.08.003

Kesler, S. R., Menon, V., and Reiss, A. L. (2006). Neurofunctional differences associated with arithmetic processing in Turner syndrome. *Cereb. Cortex* 16, 849–856. doi: 10.1093/cercor/bhj028

King, R. D., Brown, B., Hwang, M., Jeon, T., and George, A. T. (2010). Fractal dimension analysis of the cortical ribbon in mild Alzheimer’s disease. *NeuroImage* 53, 471–479. doi: 10.1016/j.neuroimage.2010.06.050

Klein, E., Moeller, K., Glauche, V., Weiller, C., and Willmes, K. (2013). Processing pathways in mental arithmetic—evidence from probabilistic fiber tracking. *PLoS One* 8:e55455. doi: 10.1371/journal.pone.0055455

Klein, E., Suchan, J., Moeller, K., Karnath, H.-O. O., Knops, A., Wood, G., et al. (2014). Considering structural connectivity in the triple code model of numerical cognition: differential connectivity for magnitude processing and arithmetic facts. *Brain Struct. Funct.* 221, 979–995. doi: 10.1007/s00429-014-0951-1

Kucian, K., Grond, U., Rotzer, S., Henzi, B., Schönmann, C., Plangger, F., et al. (2011). Mental number line training in children with developmental dyscalculia. *NeuroImage* 57, 782–795. doi: 10.1016/j.neuroimage.2011.01.070

Laux, L., Glanzmann, P., Schaffner, P., and Spielberger, C. (1981). *Das State-Trait-Angstinventar: STAI.* Weinheim: Beltz.

Li, Y., Hu, Y., Wang, Y., Weng, J., and Chen, F. (2013). Individual structural differences in left inferior parietal area are associated with schoolchildrens’ arithmetic scores. *Front. Hum. Neurosci.* 7:844. doi: 10.3389/fnhum.2013.00844

Li, X., Jiang, J., Zhu, W., Yu, C., Sui, M., Wang, Y., et al. (2007). Asymmetry of prefrontal cortical convolution complexity in males with attention-deficit/hyperactivity disorder using fractal information dimension. *Brain Dev.* 29, 649–655. doi: 10.1016/j.braindev.2007.04.008

Lohmann, G. (1998). Extracting line representations of sulcal and gyral patterns in MR images of the human brain. *IEEE Trans. Med. Imaging* 17, 1040–1048. doi: 10.1109/42.746714

Lubin, A., Rossi, S., Simon, G., Lanoë, C., Leroux, G., Poirel, N., et al. (2013). Numerical transcoding proficiency in 10-year-old schoolchildren is associated with gray matter inter-individual differences: a voxel-based morphometry study. *Front. Psychol.* 4:197. doi: 10.3389/fpsyg.2013.00197

Luders, E., Narr, K. L., Thompson, P. M., Rex, D. E., Jancke, L., Steinmetz, H., et al. (2004). Gender differences in cortical complexity. *Nat. Neurosci.* 7, 799–800. doi: 10.1038/nn1277

Luders, E., Thompson, P. M., Narr, K. L., Toga, A. W., Jancke, L., and Gaser, C. (2006). A curvature-based approach to estimate local gyrification on the cortical surface. *NeuroImage* 29, 1224–1230. doi: 10.1016/j.neuroimage.2005.08.049

Madan, C. R., and Kensinger, E. A. (2016). Cortical complexity as a measure of age-related brain atrophy. *NeuroImage* 134, 617–629. doi: 10.1016/j.neuroimage.2016.04.029

Mak, A. K. Y., Wong, M. M. C., Han, S. H., and Lee, T. M. C. (2009). Gray matter reduction associated with emotion regulation in female outpatients with major depressive disorder: a voxel-based morphometry study. *Prog. Neuropsychopharmacol. Biol. Psychiatry* 33, 1184–1190. doi: 10.1016/j.pnpbp.2009.06.025

Menon, V. (2014). “Arithmetic in the child and adult brain,” in *The Oxford Handbook of Mathematical Cognition*, eds R. Cohen Kadosh and A. Dowker (Oxford: Oxford University Press), 1–23.

Moreau, D., Wiebels, K., Wilson, A. J., and Waldie, K. E. (2019). Volumetric and surface characteristics of gray matter in adult dyslexia and dyscalculia. *Neuropsychologia* 127, 204–210. doi: 10.1016/j.neuropsychologia.2019.02.002

Narumoto, J., Okada, T., Sadato, N., Fukui, K., and Yonekura, Y. (2001). Attention to emotion modulates fMRI activity in human right superior temporal sulcus. *Cogn. Brain Res.* 12, 225–231. doi: 10.1016/s0926-6410(01)00053-2

Nenadic, I., Smesny, S., Schlösser, R. G. M., Sauer, H., and Gaser, C. (2010). Auditory hallucinations and brain structure in schizophrenia: voxel-based morphometric study. *Br. J. Psychiatry* 196, 412–413. doi: 10.1192/bjp.bp.109.070441

Núñez-Peña, M. I., Guilera, G., and Suárez-Pellicioni, M. (2014). The single-item math anxiety scale: an alternative way of measuring mathematical anxiety. *J. Psychoeduc. Assess.* 32, 306–317. doi: 10.1177/0734282913508528

Núñez-Peña, M. I., Suárez-Pellicioni, M., Guilera, G., and Mercadé-Carranza, C. (2013). A spanish version of the short mathematics anxiety rating scale (sMARS). *Learn. Individ. Differ.* 24, 204–210. doi: 10.1016/j.lindif.2012.12.009

O’Donnell, S., Noseworthy, M. D., Levine, B., and Dennis, M. (2005). Cortical thickness of the frontopolar area in typically developing children and adolescents. *NeuroImage* 24, 948–954. doi: 10.1016/j.neuroimage.2004.10.014

Pagani, E., Rocca, M. A., Gallo, A., Rovaris, M., Martinelli, V., Comi, G., et al. (2005). Regional brain atrophy evolves differently in patients with multiple sclerosis according to clinical phenotype. *Am. J. Neuroradiol.* 26, 341–346.

Parsons, S., and Bynner, J. (2005). *Does Numeracy Matter More?* London: National Research and Development Centre for Adult Literacy and Numeracy.

Peelle, J. E., Cusack, R., and Henson, R. N. A. (2012). Adjusting for global effects in voxel-based morphometry: gray matter decline in normal aging. *NeuroImage* 60, 1503–1516. doi: 10.1016/j.neuroimage.2011.12.086

Peng, J., Liu, J., Nie, B., Li, Y., Shan, B., Wang, G., et al. (2011). Cerebral and cerebellar gray matter reduction in first-episode patients with major depressive disorder: a voxel-based morphometry study. *Eur. J. Radiol.* 80, 395–399. doi: 10.1016/j.ejrad.2010.04.006

Peters, L., and De Smedt, B. (2018). Arithmetic in the developing brain: a review of brain imaging studies. *Dev. Cogn. Neurosci.* 30, 265–279. doi: 10.1016/j.dcn.2017.05.002

Pinel, P., Dehaene, S., Rivière, D., and LeBihan, D. (2001). Modulation of parietal activation by semantic distance in a number comparison task. *NeuroImage* 14, 1013–1026. doi: 10.1006/nimg.2001.0913

Price, G. R., Wilkey, E. D., Yeo, D. J., and Cutting, L. E. (2016). The relation between 1st grade grey matter volume and 2nd grade math competence. *NeuroImage* 124, 232–237. doi: 10.1016/j.neuroimage.2015.08.046

Ranpura, A., Isaacs, E., Edmonds, C., Rogers, M., Lanigan, J., Singhal, A., et al. (2013). Developmental trajectories of grey and white matter in dyscalculia. *Trends Neurosci. Educ.* 2, 56–64. doi: 10.1016/j.tine.2013.06.007

Richlan, F., Kronbichler, M., and Wimmer, H. (2013). Structural abnormalities in the dyslexic brain: a meta-analysis of voxel-based morphometry studies. *Hum. Brain Mapp.* 34, 3055–3065. doi: 10.1002/hbm.22127

Richter, M., Amunts, K., Mohlberg, H., Bludau, S., Eickhoff, S. B., Zilles, K., et al. (2019). Cytoarchitectonic segregation of human posterior intraparietal and adjacent parieto-occipital sulcus and its relation to visuomotor and cognitive functions. *Cereb. Cortex* 29, 1305–1327. doi: 10.1093/cercor/bhy245

Rosen, M. L., Stern, C. E., Michalka, S. W., Devaney, K. J., and Somers, D. C. (2015). Influences of long-term memory-guided attention and stimulus-guided attention on visuospatial representations within human intraparietal sulcus. *J. Neurosci.* 35, 11358–11363. doi: 10.1523/JNEUROSCI.1055-15.2015

Rotzer, S., Kucian, K., Martin, E., von Aster, M., Klaver, P., and Loenneker, T. (2008). Optimized voxel-based morphometry in children with developmental dyscalculia. *NeuroImage* 39, 417–422. doi: 10.1016/j.neuroimage.2007.08.045

Rykhlevskaia, E., Uddin, L. Q., Kondos, L., and Menon, V. (2009). Neuroanatomical correlates of developmental dyscalculia: combined evidence from morphometry and tractography. *Front. Hum. Neurosci.* 3:51. doi: 10.3389/neuro.09.051.2009

Sandu, A. L., Staff, R. T., McNeil, C. J., Mustafa, N., Ahearn, T., Whalley, L. J., et al. (2014). Structural brain complexity and cognitive decline in late life—a longitudinal study in the Aberdeen 1936 Birth Cohort. *NeuroImage* 100, 558–563. doi: 10.1016/j.neuroimage.2014.06.054

Schillinger, F. L., Vogel, S. E., Diedrich, J., and Grabner, R. H. (2018). Math anxiety, intelligence, and performance in mathematics: insights from the German adaptation of the Abbreviated Math Anxiety Scale (AMAS-G). *Learn. Individ. Differ.* 61, 109–119. doi: 10.1016/j.lindif.2017.11.014

Schnack, H. G., Van Haren, N. E. M., Brouwer, R. M., Evans, A., Durston, S., Boomsma, D. I., et al. (2015). Changes in thickness and surface area of the human cortex and their relationship with intelligence. *Cereb. Cortex* 25, 1608–1617. doi: 10.1093/cercor/bht357

Shaw, P., Malek, M., Watson, B., Sharp, W., Evans, A., and Greenstein, D. (2012). Development of cortical surface area and gyrification in attention-deficit/hyperactivity disorder. *Biol. Psychiatry* 72, 191–197. doi: 10.1016/j.biopsych.2012.01.031

Shen, W., Yuan, Y., Liu, C., and Luo, J. (2017). The roles of the temporal lobe in creative insight: an integrated review. *Think. Reason.* 23, 321–375. doi: 10.1080/13546783.2017.1308885

Simon, O., Mangin, J. F., Cohen, L., Le Bihan, D., and Dehaene, S. (2002). Topographical layout of hand, eye, calculation and language-related areas in the human parietal lobe. *Neuron* 33, 475–487. doi: 10.1016/s0896-6273(02)00575-5

Simos, P. G., Kanatsouli, K., Fletcher, J. M., Sarkari, S., Juranek, J., Cirino, P., et al. (2008). Aberrant spatiotemporal activation profiles associated with math difficulites in children: a magnetic source imaging study. *Neuropsychology* 22, 571–584. doi: 10.1037/0894-4105.22.5.571

Sokolowski, H. M., Fias, W., Mousa, A., and Ansari, D. (2017). Common and distinct brain regions in both parietal and frontal cortex support symbolic and nonsymbolic number processing in humans: a functional neuroimaging meta-analysis. *NeuroImage* 146, 376–394. doi: 10.1016/j.neuroimage.2016.10.028

Spielberger, C. D., Gorsuch, R. L., Lushene, P. R., Vagg, P. R., and Jacobs, A. G. (1983). *Manual for the State-Trait Anxiety Inventory.* Palo Alto, CA: Consulting Psychologists Press.

Starke, M., Kiechl-Kohlendorfer, U., Kucian, K., Pupp Peglow, U., Kremser, C., Schocke, M., et al. (2013). Brain structure, number magnitude processing and math proficiency in 6- to 7-year-old children born prematurely: a voxel-based morphometry study. *Neuroreport* 24, 419–424. doi: 10.1097/WNR.0b013e32836140ed

Steinbrink, C., Vogt, K., Kastrup, A., Müller, H. P., Juengling, F. D., Kassubek, J., et al. (2008). The contribution of white and gray matter differences to developmental dyslexia: insights from DTI and VBM at 3.0 T. *Neuropsychologia* 46, 3170–3178. doi: 10.1016/j.neuropsychologia.2008.07.015

Supekar, K., Swigart, A. G., Tenison, C., Jolles, D. D., Rosenberg-Lee, M., Fuchs, L., et al. (2013). Neural predictors of individual differences in response to math tutoring in primary-grade school children. *Proc. Natl. Acad. Sci. U S A* 110, 8230–8235. doi: 10.1073/pnas.1222154110

Thompson, P. M., Hayashi, K. M., Sowell, E. R., Gogtay, N., Giedd, J. N., Rapoport, J. L., et al. (2004). Mapping cortical change in Alzheimer’s disease, brain development and schizophrenia. *NeuroImage* 23, 2–18. doi: 10.1016/j.neuroimage.2004.07.071

Thompson, P. M., Schwartz, C., Lin, R. T., Khan, A. A., and Toga, A. W. (1996). Three-dimensional statistical analysis of sulcal variability in the human brain. *J. Neurosci.* 16, 4261–4274. doi: 10.1523/jneurosci.16-13-04261.1996

Tosoni, A., Pitzalis, S., Committeri, G., Fattori, P., Galletti, C., and Galati, G. (2015). Resting-state connectivity and functional specialization in human medial parieto-occipital cortex. *Brain Struct. Funct.* 220, 3307–3321. doi: 10.1007/s00429-014-0858-x

Vanni, S., Tanskanen, T., Seppä, M., Uutela, K., and Hari, R. (2001). Coinciding early activation of the human primary visual cortex and anteromedial cuneus. *Proc. Natl. Acad. Sci. U S A* 98, 2776–2780. doi: 10.1073/pnas.041600898

Voets, N. L., Hough, M. G., Douaud, G., Matthews, P. M., James, A., Winmill, L., et al. (2008). Evidence for abnormalities of cortical development in adolescent-onset schizophrenia. *NeuroImage* 43, 665–675. doi: 10.1016/j.neuroimage.2008.08.013

Vogel, S. E., Goffin, C., Bohnenberger, J., Koschutnig, K., Reishofer, G., Grabner, R. H., et al. (2017a). The left intraparietal sulcus adapts to symbolic number in both the visual and auditory modalities: evidence from fMRI. *NeuroImage* 153, 16–27. doi: 10.1016/j.neuroimage.2017.03.048

Vogel, S. E., Haigh, T., Sommerauer, G., Spindler, M., Brunner, C., Lyons, I. M., et al. (2017b). Processing the order of symbolic numbers: a reliable and unique predictor of arithmetic fluency. *J. Numer. Cogn.* 3, 288–308. doi: 10.5964/jnc.v3i2.55

White, T., Andreasen, N. C., Nopoulos, P., and Magnotta, V. (2003). Gyrification abnormalities in childhood- and adolescent-onset schizophrenia. *Biol. Psychiatry* 54, 418–426. doi: 10.1016/s0006-3223(03)00065-9

Wilkey, E. D., Barone, J. C., Mazzocco, M. M. M., Vogel, S. E., and Price, G. R. (2017). The effect of visual parameters on neural activation during nonsymbolic number comparison and its relation to math competency. *NeuroImage* 159, 430–442. doi: 10.1016/j.neuroimage.2017.08.023

Wilkey, E. D., Cutting, L. E., and Price, G. R. (2018). Neuroanatomical correlates of performance in a state-wide test of math achievement. *Dev. Sci.* 21:2. doi: 10.1111/desc.12545

Yang, J. J., Yoon, U., Yun, H. J., Im, K., Choi, Y. Y., Lee, K. H., et al. (2013). Prediction for human intelligence using morphometric characteristics of cortical surface: partial least square analysis. *Neuroscience* 246, 351–361. doi: 10.1016/j.neuroscience.2013.04.051

Keywords: right superior temporal gyrus, right central sulcus, left parieto-occipital sulcus, surface-based morphometry, mathematics

Citation: Heidekum AE, Vogel SE and Grabner RH (2020) Associations Between Individual Differences in Mathematical Competencies and Surface Anatomy of the Adult Brain. Front. Hum. Neurosci. 14:116. doi: 10.3389/fnhum.2020.00116

Received: 15 October 2019; Accepted: 13 March 2020;

Published: 27 March 2020.

Edited by:

Elise Klein, University of Tübingen, GermanyReviewed by:

Tanya M. Evans, University of Virginia, United StatesMario Pinto, Santa Lucia Foundation (IRCCS), Italy

Copyright © 2020 Heidekum, Vogel and Grabner. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Alexander E. Heidekum, alexander.heidekum@uni-graz.at