# Surface-Based Morphometry of Cortical Thickness and Surface Area Associated with Heschl's Gyri Duplications in 430 Healthy Volunteers

Damien Marie, Sophie Maingault, Fabrice Crivello, Bernard Mazoyer, Nathalie Tzourio-Mazoyer

► To cite this version:

Damien Marie, Sophie Maingault, Fabrice Crivello, Bernard Mazoyer, Nathalie Tzourio-Mazoyer. Surface-Based Morphometry of Cortical Thickness and Surface Area Associated with Heschl's Gyri Duplications in 430 Healthy Volunteers. Frontiers in Human Neuroscience, 2016, 10, pp.69. 10.3389/fnhum.2016.00069. hal-01382649

HAL Id: hal-01382649

https://hal.science/hal-01382649

Submitted on 17 Oct 2016

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_2.png" /></figure>

# Surface-Based Morphometry of Cortical Thickness and Surface Area Associated with Heschl's Gyri Duplications in 430 Healthy Volunteers

Damien Marie¹,²,³, Sophie Maingault¹,²,³, Fabrice Crivello¹,²,³, Bernard Mazoyer¹,²,³ and Nathalie Tzourio-Mazoyer¹,²,³*

¹Groupe d'Imagerie Neurofonctionnelle, Institut des Maladies Neurodégénératives UMR 5293, Université de Bordeaux, Bordeaux, France, ²Centre National de la Recherche Scientifique, Groupe d'Imagerie Neurofonctionnelle, Institut des Maladies Neurodégénératives UMR 5293, Bordeaux, France, ³Commissariat à l'Énergie Atomique et aux Énergies Alternatives, Groupe d'Imagerie Neurofonctionnelle, Institut des Maladies Neurodégénératives UMR 5293, Bordeaux, France

We applied Surface-Based Morphometry to assess the variations in cortical thickness (CT) and cortical surface area (CSA) in relation to the occurrence of Heschl's gyrus (HG) duplications in each hemisphere. 430 healthy brains that had previously been classified as having a single HG, Common Stem Duplication (CSD) or Complete Posterior Duplication (CPD) in each hemisphere were analyzed. To optimally align the HG area across the different groups of gyrification, we computed a specific surface-based template composed of 40 individuals with a symmetrical HG gyrification pattern (20 single HG, 10 CPD, 10 CSD). After normalizing the 430 participants' T1 images to this specific template, we separately compared the groups constituted of participants with a single HG, CPD, and CSD in each hemisphere. The occurrence of a duplication in either hemisphere was associated with an increase in CT posterior to the primary auditory cortex. This may be the neural support of expertise or great abilities in either speech or music processing domains that were related with duplications by previous studies. A decrease in CSA in the *planum temporale* was detected in cases with duplication in the left hemisphere. In the right hemisphere, a medial decrease in CSA and a lateral increase in CSA were present in HG when a CPD occurred together with an increase in CSA in the depth of the superior temporal sulcus (STS) in CSD compared to a single HG. These variations associated with duplication might be related to the functions that they process jointly within each hemisphere: temporal and speech processing in the left and spectral and music processing in the right.

**Keywords:** Heschl's gyrus, gyrification, anatomy, MRI, FreeSurfer, cortical thickness, cortical surface area, hemispheric specialization

## INTRODUCTION

Heschl's gyrus (HG), which hosts the Primary Auditory Cortex (PAC), exhibits very high anatomical variability in its size and gyrification across both individuals and hemispheres (Penhune et al., 1996; Leonard et al., 1998; Abdul-Kareem and Sluming, 2008; Marie et al., 2015). Three configurations of HG can be observed: a single HG, a Common Stem Duplication (CSD),

## OPEN ACCESS

**Edited by:**
Bogdan Draganski,
University of Lausanne, Switzerland

**Reviewed by:**
Stefan Elmer,
University of Zurich, Switzerland
Virginia Penhune,
Concordia University, Canada
Christian Lambert,
St George's University of London, UK

**\*Correspondence:**
Nathalie Tzourio-Mazoyer
nathalie.tzourio-mazoyer@u-bordeaux.fr;
ntzouriomazoyer@gmail.com

**Received:** 21 April 2015
**Accepted:** 11 February 2016
**Published:** 07 March 2016

**Citation:**
Marie D, Maingault S, Crivello F, Mazoyer B and Tzourio-Mazoyer N (2016) Surface-Based Morphometry of Cortical Thickness and Surface Area Associated with Heschl's Gyri Duplications in 430 Healthy Volunteers.
Front. Hum. Neurosci. 10:69.
doi: 10.3389/fnhum.2016.00069
corresponding to a partial split of the lateral part of the gyrus by the *sulcus intermedius* of Beck, and a Complete Posterior Duplication (CPD) that is a double HG (Abdul-Kareem and Sluming, 2008). Although, the PAC and HG do not exactly overlap at the microscopic level (Morosan et al., 2001). The PAC may seat in the medial two-thirds of a single HG and is limited to the first or anterior HG (aHG) in the case of duplication. MRI methods that provide estimates of myelination have demonstrated a highly myelinated area running over the medial two-thirds of HG in single gyrus and over the aHG in duplications (Dick et al., 2012). This myelinated area likely corresponds to the location of the PAC, primary areas being more heavily myelinated than the associative areas (Glasser and Van Essen, 2011). Using myelin-stained sections, Pfeifer reported one of the first demonstrations of a coupling between acoustic radiations and the gyrus in single HG, while only in aHG in case of duplication (Pfeifer, 1920). Post-mortem probabilistic maps of acoustic radiations further confirmed Pfeifer's pioneering observations (Rademacher et al., 2002).

There are numerous reports of associations between the anatomical characteristics of HG and long-term auditory learning for music and speech, or expertise in either domain (for reviews see Zatorre et al., 2012; Zatorre, 2013). These reports are within the frame of the spectro-temporal framework of acoustic processing (Zatorre and Belin, 2001; Hickok and Poeppel, 2007). For example, Warrier has demonstrated a direct anatomo-functional relationship in HG during sound processing, where the volume of activation during the temporal processing of sound is positively correlated with the volume of aHG on the left but not the right, and the opposite is observed during the spectral processing of sound (Warrier et al., 2009). Schneider combined MRI with MEG to demonstrate a left hemispheric specialization for (missing) fundamental pitch perception and a right hemispheric specialization for spectral pitch perception, consistent with this spectro-temporal framework of acoustic processing (Schneider et al., 2005). In fact, he demonstrated that listeners who completed the pitch change judgments on the fundamental pitch had larger gray matter (GM) volumes in the left HG than did listeners performing the pitch change judgment on the basis of the spectral envelope, who instead had an increase in the right HG GM volume.

Comparisons of the anatomical characteristics of HG with expertise also demonstrated lateralized results (summary in **Table 1**). In the left hemisphere, a larger volume of white matter (WM) and/or GM of aHG or HG have been associated with good performances in phonological awareness and literacy abilities in 10-year-old children (Sutherland et al., 2012), with expertise or proficiency in phonological processing of foreign speech sounds in adults (Golestani et al., 2007, 2011; Wong et al., 2008) and with early bilingualism (Ressel et al., 2012). Left HG duplications have been associated with phonological expertise (Golestani et al., 2011) but also with phonological disorders and dyslexia (Leonard et al., 1993; Altarelli et al., 2014). Moreover, a large second HG has been reported in dyslexic adults (Leonard et al., 2001).

On the right, a larger GM and/or WM volume of HG has been reported in professional musicians with absolute pitch compared to musicians without absolute pitch, and the proficiency of an

individual's absolute pitch ability has been shown to be positively correlated with GM volume (Wengenroth et al., 2013). A greater GM density was observed in HG of musicians compared to non-musicians (Bermudez et al., 2009), and a lateral expansion of HG was detected in 6-year-old children who completed 15 months of instrumental training compared to control children (Hyde et al., 2009). Applying surface-based morphometry (SBM) on anatomical images allowed several investigators to report variations in Cortical Thickness (CT) in relation to music expertise or learning in HG (notably reviewed in Zatorre, 2013), such as a CT increase in the right HG associated with musicianship (Bermudez et al., 2009) and pitch discrimination performances (Foster and Zatorre, 2010). Although, less direct, Liem who indicated that the N1 amplitude is an ERP component related to musical expertise, found a negative correlation between N1 amplitude and CT in HG (Liem et al., 2012). Finally, a high occurrence of right duplications had been notably observed in musicians having absolute pitch (Schneider et al., 2005).

The works cited above, which investigated the anatomo-functional relationships in HG in relation to performances during speech or musical sound processing, reported differences in the GM volume, CT, or the occurrence of duplications. Although we have recently shown that the presence of duplications is associated with a decrease in the surface area of aHG (Marie et al., 2015), the changes in CT or Cortical Surface Area (CSA) associated with the occurrence of duplications are still unknown. We applied SBM to the anatomical images of 430 participants who had previously been characterized according to their HG macroscopical anatomy (single HG, CSD, and CPD, Marie et al., 2015) to characterize the differences in CT and CSA in groups with varying HG duplication patterns.

# MATERIALS AND METHODS

## Participants

The characteristics of the 430 participants are fully described in our previous article (Marie et al., 2015). All were free of brain abnormalities, and the Basse-Normandie's ethical committee approved the study. All participants provided their informed written consent and received compensation for their participation in the framework of the BIL&GIN database (Mazoyer et al., 2015). Briefly, the population is balanced for handedness (self-report, 232 right-handers vs. 198 left-handers, Chi square = 2.69, *p* = 0.10), and gender (209 men vs. 221 women, Chi square = 0.33, *p* = 0.56). The mean age of the participants is 26 years (SD, 8). Each volunteer's skull perimeter was measured at the level of the eyebrows, passing at the top edge of the ears and along the occipital bump. The mean skull perimeter was 57 cm (SD, 2) and men had larger skull perimeters than women did (+3 cm, *t* = 17.1, *p* < 0.0001). Note, that none of these individual variables differed according to the left or right HG number (Marie et al., 2015).

## Image Acquisition

Anatomical images were acquired from 2007 to 2011 using the same 3T Philips Achieva scanner (Philips Medical Systems, Best, The Netherlands). High-resolution T1-weighted images were
<table><caption>TABLE 1 | Meta-analysis of reports on HG structure changes related to phonology, music, sound processing, or pathology.</caption><thead><tr><th></th><th>Population</th><th>Left HG</th><th>Sound temporal or phonological processing</th><th>Right HG</th><th>Sound tonal or musical processing</th></tr></thead><tbody><tr><td>Golestani et al., 2011</td><td>Adults</td><td>↑ Duplication<br/>↑ GM-WM totHG</td><td>Expertise phonetics</td><td>↑ GM-WM totHG</td><td></td></tr><tr><td>Sutherland et al., 2012</td><td>Children</td><td>↑ GM aHG</td><td>Phonological awareness, literacy abilities</td><td></td><td></td></tr><tr><td>Golestani et al., 2007</td><td>Adults</td><td>↑ WM aHG<br/>↑ Duplication</td><td>Foreign sound learning</td><td></td><td></td></tr><tr><td>Wong et al., 2008</td><td>Adults</td><td>↑ GM-WM aHG</td><td>Foreign sound learning</td><td></td><td></td></tr><tr><td>Ressel et al., 2012</td><td>Adults</td><td>↑ GM aHG</td><td>Early bilinguism</td><td></td><td></td></tr><tr><td>Schneider et al., 2005</td><td>Adults</td><td>↑ GM HG-PP</td><td>Pitch judgment on fundamental pitch</td><td>↑ Duplication<br/>↑ GM HG-PP</td><td>Pitch judgment on spectral envelope</td></tr><tr><td>Leonard et al., 1993</td><td>Adults</td><td>↑ Duplication</td><td>Dyslexia</td><td>↑ Duplication</td><td></td></tr><tr><td>Leonard et al., 2001</td><td>Adults</td><td>↑ CSA postHG</td><td>Dyslexia</td><td></td><td></td></tr><tr><td>Ma et al., 2015</td><td>Children</td><td></td><td>Dyslexia</td><td>↑ CT</td><td></td></tr><tr><td>Altarelli et al., 2014</td><td>Boys</td><td></td><td>Dyslexia</td><td>↑ Duplication</td><td></td></tr><tr><td>Warrier et al., 2009</td><td>Adults</td><td>↑ GM-WM aHG</td><td>Temporal processing</td><td>↑ GM-WM aHG</td><td>Spectral processing</td></tr><tr><td>Schneider et al., 2002</td><td>Adults</td><td>↑ GM aHG</td><td></td><td>↑ GM aHG<br/>↑ GM medial aHG</td><td>Musical expertise<br/>AMMA tonal test, N19m-P30m<br/>auditory ERP</td></tr><tr><td>Bermudez et al., 2009</td><td>Adults</td><td></td><td></td><td>↑ GM density<br/>↑ CT aHG</td><td>Musical expertise</td></tr><tr><td>Wengenroth et al., 2013</td><td>Adults</td><td></td><td></td><td>↑ GM HG</td><td>Musical expertise, AP</td></tr><tr><td>Hyde et al., 2009</td><td>Children</td><td></td><td></td><td>↑ GM-WM HG*</td><td>Instrumental training, melodic<br/>and rhythmic discrimination test</td></tr><tr><td>Foster and Zatorre, 2010</td><td>Adults</td><td></td><td></td><td>↑ GM and CT HS</td><td>Musical expertise</td></tr><tr><td>Liem et al., 2012</td><td>Adults</td><td></td><td></td><td>↓ CT HG and HS</td><td>N1 auditory ERP</td></tr></tbody></table>

aHG, anterior Heschl's gyrus; AP, Absolute Pitch; CT, Cortical Thickness; ERP, Event-Related Potentials; GM, Gray Matter; HG, Heschl's gyrus; HS, Heschl's sulcus; postHG, posterior Heschl's gyrus in case of duplication; PP, Planum Polare; STG, Superior Temporal Gyrus; totHG, total Heschl's gyrus including both convolutions in case of duplication; WM, White Matter; GM-WM, both GM and WM; *, the volume was measured using a measure of relative voxel expansion or contraction.

obtained using a 3D-FFE-TFE sequence (TR, 20 ms; TE, 4.6 ms; flip angle, 10°, inversion time, 800 ms; turbo field echo factor, 65; Sense factor, 2; field of view, 256 × 256 × 180 mm; and isotropic voxel, 1 mm³). For each participant, the line between the anterior and posterior commissures was identified on the midsagittal plane, and the T1-MRI volume was acquired after orienting the brain in the bi-commissural coordinate system.

## Image Analysis

### Definition of the HG Gyrification Patterns

A single expert (NT-M) identified HG and its gyrification type on individual images in the native space as fully described elsewhere (Marie et al., 2015). Briefly, the definition of the number of HG was optimized via the knife-cut method, which is based on the reconstruction of an oblique section plane passing through the Sylvian fissure (Kulynych et al., 1993; Tzourio-Mazoyer et al., 2010) computed with the aid of dedicated software (Voxeline, Diallo et al., 1998). An oblique section plane parallel to the Sylvian fissure was generated by uncovering the superior

temporal plane and passing through the point where HG is the largest. Axial, coronal, sagittal, and oblique slices are displayed simultaneously for identifying duplications.

The gyrification pattern of HG (single HG, common stem duplication or CSD, or complete posterior duplication, or CPD) was defined according to the anatomical criteria reviewed in Abdul-Kareem and Sluming (2008). In cases with a single HG, the anterior border is constituted by the temporal transverse sulcus, and the posterior border is Heschl's sulcus. In cases with a CSD, the gyrus is partially split by the *sulcus intermedius* of Beck (SI), which never reaches the internal border of the gyrus. We considered a CSD to be present if the SI length was at least one-third of the length of HG. There are two Heschl's sulci in cases with a CPD; the first Heschl's sulcus splits the gyrus medially into two parts and the second Heschl's sulcus becomes the posterior border of HG. As previously reported, we identified 284 single HG (66%), 98 CSDs (23%), and 48 CPDs (11%) in the left hemisphere, and 239 single HG (56%), 95 CSDs (22%), and 96 CPDs (22%) in the right hemisphere (Marie et al., 2015).
## Reproducibility of the Identification of the Gyrification Pattern

Reproducibility was measured by a second assessment of 40 randomly selected brains (80 hemispheres) by the same observer. The concordance between the two evaluations was measured with the chance corrected Kappa statistics (JMP pro 11). The Kappa statistics for the concordance of the two assessments of the number of gyri was 0.84 ± 0.06 (p < 0.0001). The Kappa statistics for the concordance of the two assessments of the HG pattern (single, CPD or CSD) was 0.82 ± 0.06 (p < 0.0001). According to Cohen's original paper, a Kappa value > 0.81 attests an excellent reproducibility (Cohen, 1960). Non-concordant classification occurred in eight cases (10% of the hemispheres). Four of the non-concordant associations corresponded to a confusion between CPD and CSD; in two cases, the non-concordance was between CSD and single, and in two others, it was between single and CPD in the right hemisphere.

## Surface-Based Morphometry

The cortical surfaces were reconstructed and the anatomical features were measured for each participant with the FreeSurfer 5.3.0 image analysis suite (http://surfer.nmr.mgh.harvard.edu/), (Dale et al., 1999; Fischl et al., 1999a,b, 2001, 2002, 2004; Fischl and Dale, 2000; Fischl, 2012). The FreeSurfer analysis stream includes intensity bias field removal optimized for the present 3T acquisitions and skull stripping. In the native anatomical space of the subject, FreeSurfer constructs models of the cortical surface that can be represented as a mesh of triangles and tessellated into over 160,000 vertices for each hemisphere. The vertex positions are adjusted such that the surface follows the T1 intensity gradient between the cortical WM and GM. This highly folded surface, called “white surface,” can be inflated to see buried surfaces within the sulci. A second surface, called the “pial surface,” is fit to the outer GM surface of the brain. The hemispheres are processed separately through a precise estimation of the inter-hemispheric fissure plane.

## Quality Control

The problems of estimating the surfaces around the thin temporal lobe strands that border the tip of the lateral ventricle have been identified with the FreeSurfer segmentation procedure (Dale et al., 1999). In this case, a reduction of the intensity of the voxels in the corresponding white matter leads to a misclassification of voxels and an inaccurate definition of the white and pial surfaces. To overcome this limitation, we visually checked (SM) each of the 430 pial and white surfaces on each axial, sagittal, and coronal section from each participant twice and identified surface reconstruction errors in 44 cases (10%). A correction procedure consisting of manually adding an average of nine control points (SD, 7) in the misclassified white matter was then applied in the 44 individuals. Then, the FreeSurfer procedure was launched a second time for these corrected individuals, resulting in more accurate surface reconstructions (Figure 1).

## Calculation of the Morphometric Variables in the Native Space of Each Participant

For each vertex, CT is defined as the distance between the white and pial surfaces and CSA is defined as the average area of

<figure><img src="image_3.png" /><figcaption>FIGURE 1 | An example of manual corrections applied in case of wrong estimation of pial and/or white cortical surfaces. In this individual, the FreeSurfer software did not correctly estimate both white and pial surfaces in the upper part of the temporal pole. Manually adding checkpoints to the white matter surface results in better surface estimation after a second extraction process (yellow line, white surface; red line, pial surface; green dots, checkpoints).</figcaption></figure>

the triangles to which the vertex belongs. The cortical surface curvature (CURV) parameter was also extracted and used to drive the non-linear surface-based inter-subject registration. The curvature is estimated as 1/r, where r is the radius of the inscribed circle to each point (or vertex) of the white surface.

## Definition of a Specific Template for Between-Subject Registration and Smoothing

A specific surface-based template for between-subject registration that included images from individuals with identical HG duplication patterns in the left and right hemispheres was computed to avoid a left alignment weighted toward the single HG scenario and a right alignment weighted toward the duplication scenario (HG anatomical pattern in Fsaverage, see Supplementary Figure 1). Forty individuals of the sample were included (26 men, 18 right-handers), corresponding to 20 bilateral single HG and 20 bilateral duplications (10 CSD and 10 CPD). The procedure for computing this specific surface-based template, called 40-average, is fully described and illustrated in Supplementary Material and Supplementary Figure 2.

Non-linear surface-based normalization of the 430 CT and CSA individual maps to this specific template was then applied. An additional step was performed by adding a Jacobian correction to the CSA maps to account for any stretching or compression in the registration that occurred during the nonlinear registration (Winkler et al., 2012). Note, that this correction is not needed for the CT maps because CT is measured along a vector normal to any stretching or compression.

After registration, the individual CT and CSA maps were surface-smoothed by a 10 mm FWHM Gaussian filter size, statistical power of both features benefiting from such optimal smoothing filter size as demonstrated by Liem et al. (2015). However, minimal detectable differences depending on smoothing filter size (Liem et al., 2015), individual CT maps of two groups (left CPD, and left single HG) were also surface-smoothed with a 5 mm FWHM Gaussian filter in a complementary analysis in order to evaluate the impact of a reduced filter size on structural differences detected between HG anatomical patterns.
## Curvature Variance Maps

To characterize the intrinsic anatomical variability of each HG configuration, we computed CURV variance maps in the 40-average space in both hemispheres. For each vertex, we calculated a CURV variance value using MATLAB 7.11.1 (R2010b), which was defined as the mean of the square minus the square of the mean for a given group (for example the calculation of curvature variance in the left hemisphere for the 284 subjects with single HG).

## Statistical Analysis

To compare the groups with varying HG duplication patterns, we ran two independent ANCOVA in each hemisphere on CT and CSA separately, with the HG pattern as the main factor. All analyses were performed with the QDEC module included in FreeSurfer, and the participant's skull perimeter was included as a confounding variable to account for differences in head size. The statistical threshold was set at *p* < 0.05 FDR corrected, and an additional cluster extension threshold was set at 20 mm².

To display the results on the supratemporal plane, a specific patch of the temporal lobe was computed for the mean image of each group and hemisphere by merging the corresponding labels of Desikan's parcellation provided by FreeSurfer (Desikan et al., 2006). The removal of the fusiform, entorhinal, and parahippocampal gyri was performed to clean up the patch (Figure 3).

Finally, the overlap of the significant variations in CT and CSA were computed by summing the binarized image of each anatomical feature (Figure 4).

# RESULTS

## Macroscopic Anatomical Variability of HG Configuration across Hemispheres: Analysis of the CURV Variance Maps

In cases with a single HG, the variability was low, except in the antero-lateral border of the left HG (Figure 2A). In cases with a CSD in the right hemisphere, a small cluster of high variability was observed in the *sulcus intermedius* of Beck (Figure 2B), similar to that on the left, which corresponded to a variation in the width of the mid-part of the second sulcus. A high variability of left CPD was observed and corresponded to the variation in the curvature pattern of the posterior position of the first Heschl's sulcus. On the right, a comparable pattern was present, but with a much lower amplitude.

## Morphometric Changes According to the Duplication Types
CT Differences

A larger CT was observed in the CSD group compared to the single HG group all along Heschl's sulcus and the second gyrus (Table 2, Figure 3).

A larger CT was also observed in the medial part of the second HG in both hemispheres in the CPD group compared to the single HG group, but it was more medially situated on the right.

Note, that there was no difference between CSD and CPD.

These increases in CT in case of duplication, whatever the side and whatever the type of comparison (either CSD or CPD minus

single), were located at the ventral limit of the second gyrus. In addition, in case of CSD, the increase in CT extended from the ventral part of the second gyrus toward the *sulcus intermedius* of Beck (Figure 4).

## CSA Differences

In the right hemisphere, the CSA was significantly smaller in the medial part of aHG in the CPD group compared with the single HG or CSD groups, whereas a larger CSA was observed in the lateral part of the second HG of the CPD group compared with the other two groups (Table 3, Figure 3). In the left hemisphere, a decreased CSA was detected on the medial part of the second gyrus of CPD as compared to the single HG group, corresponding to the location of the planum temporale (PT) in single HG cases.

In the middle part of Heschl's sulcus, smaller CSA values were observed in both hemispheres of the CSD group compared to the single HG group. Interestingly, the extent of the cluster was larger on the left and expanded in the PT.

In the medial part of the temporal transverse sulcus, a larger CSA was also detected in both hemispheres of the CSD group compared to the single HG group. Moreover, in the right hemisphere, a larger CSA was observed directly under the lateral part of HG in the depth of the middle part of the Superior Temporal Sulcus (STS) (Figure 3).

<figure><img src="image_4.png" /></figure>
<table><caption>TABLE 2 | Local variations of Cortical Surface Area (CSA) detected in “CPD-single HG”, “CSD-single HG,” and “CPD-CSD” contrasts.</caption><thead><tr><th>Anatomical localization</th><th>Cluster size (mm²)</th><th>pmax</th><th>x</th><th>y</th><th>z</th><th>Dupl. mean (SD) (mm²)</th><th>single HG mean (SD) (mm²)</th></tr></thead><tbody><tr><td colspan="8"><strong>L: CPD-SINGLE HG</strong></td></tr><tr><td>Medial 2nd HG</td><td>27</td><td>3.10<sup>-6</sup></td><td>-46</td><td>-30</td><td>7</td><td>26 (4)</td><td>30 (5)</td></tr><tr><td colspan="8"><strong>R: CPD-SINGLE HG</strong></td></tr><tr><td>Medial aHG</td><td>161</td><td>9.10<sup>-9</sup></td><td>40</td><td>-23</td><td>10</td><td>155 (28)</td><td>176 (31)</td></tr><tr><td>Lateral 2nd HG</td><td>176</td><td>10<sup>-6</sup></td><td>64</td><td>-18</td><td>6</td><td>192 (31)</td><td>171 (27)</td></tr><tr><td colspan="8"><strong>L: CSD-SINGLE HG</strong></td></tr><tr><td>HS-PT</td><td>187</td><td>2.10<sup>-6</sup></td><td>-58</td><td>-28</td><td>8</td><td>197 (33)</td><td>214 (36)</td></tr><tr><td>Medial TTS</td><td>52</td><td>3.10<sup>-5</sup></td><td>-30</td><td>-27</td><td>12</td><td>60 (7)</td><td>57 (6)</td></tr><tr><td colspan="8"><strong>R: CSD-SINGLE HG</strong></td></tr><tr><td>Middle STS</td><td>172</td><td>3.10<sup>-7</sup></td><td>53</td><td>-17</td><td>-3</td><td>193 (29)</td><td>177 (28)</td></tr><tr><td>Medial TTS</td><td>88</td><td>8.10<sup>-7</sup></td><td>32</td><td>-26</td><td>14</td><td>95 (8)</td><td>91 (8)</td></tr><tr><td>HS</td><td>43</td><td>3.10<sup>-6</sup></td><td>49</td><td>-24</td><td>7</td><td>39 (8)</td><td>45 (8)</td></tr><tr><td colspan="8"><strong>R: CPD-CSD</strong></td></tr><tr><td>Medial aHG</td><td>126</td><td>4.10<sup>-10</sup></td><td>38</td><td>-25</td><td>13</td><td>126 (19)</td><td>145 (23)</td></tr><tr><td>Lateral 2nd HG</td><td>156</td><td>4.10<sup>-8</sup></td><td>57</td><td>-25</td><td>10</td><td>188 (34)</td><td>159 (26)</td></tr></tbody></table>

Dupl., duplication; HG, Heschl's gyrus; aHG, anterior Heschl's gyrus; HS, Heschl's sulcus; PT, planum temporale; TTS, temporal transverse sulcus; SD, Standard Deviation; STS, superior temporal sulcus; *p* < 0.05 FDR corrected for multiple comparison; cluster extension threshold, 20 mm²; x, y, z, Talairach coordinates rounded to the whole number.

## Spatial Overlap between the CSA and CT Variations

The decreased CSA observed in the left and right Heschl's sulci overlapped with the area of the increased CT in both the CPD minus single and CSD minus single comparisons (Figure 4).

## DISCUSSION

Surface-Based Morphometry allowed us to investigate changes in macroscopical neuroanatomy without being limited to the intrinsic voxel resolution (Fischl and Dale, 2000). The high spatial normalization provided by the surface-based inter-subject alignment of curvature in Freesurfer enabled us to describe the changes in CT and CSA associated with the HG gyrification patterns, completing the results previously obtained by measuring the aHG surface area using manual delineation (Marie et al., 2015).

## Methodological Considerations

The identification of the HG gyrification patterns had an excellent reproducibility according to the Kappa measure. A first type of non-concordance was related to the distinction between a single gyrus and a CSD because the criterion for CSD was the minimal length of the sulcus intermedius of Beck (30%). A second type of non-concordance was the confusion between CPD and CSD, which was related to the determination of whether the sulcus splitting the gyrus reaches the internal limit of HG. This ambiguity is generally solved by applying the second criteria based on the shape of the duplication, which resembles an “M” on sagittal and coronal slices in cases with a CPD or a “heart” in cases with a CSD. The two cases of confusion between single and CPD were located in the right hemisphere, and corresponded to a difficulty in discriminating a second gyrus from the right PT.

The second methodological point concerns the use of a specific surface-based template in the present study. While we revealed that the use of such specific template improved

the inter-subject registration of HG anatomical pattern, in particular for the group of individuals presenting a complete HG duplication in the left hemisphere (see mean images presented in Supplementary Material and Supplementary Figure 1 for a detailed description of this improvement), we also showed, that the use of such specific template did not modulate the statistical differences detected between groups for both CT and CSA anatomical (see Supplementary Material and Supplementary Figure 3). The present absence of differences between results detected by statistical analyses obtained after a normalization performed with the specific 40-average template and the default Fsaverage template strongly suggests that the accurate alignment of HG anatomical pattern obtained with the specific 40-average surface-based template has no influence on the statistical maps resulting from group comparisons, mainly due to the lower spatial resolution of these statistical maps as compared to the resolution of the mean curvatures images. In addition, the results we obtained in the CPD minus single HG analysis with a smoothing filter size of 5 mm are similar to the results deriving from the analysis performed with a 10 mm smoothing filter, suggesting there is no impact of the size of the smoothing filter on CT variations related to HG structure changes (see Supplementary Results).

## Variability of Each HG Pattern

We conducted an analysis of variance on CURV, which clearly showed a high variance on the left, increasing from single to CSD to CPD. In fact, difference of variability between hemispheres might not be significant (due to the absence of any equivalence between the vertices of the two hemispheres, statistical comparison was not possible). In both hemispheres, a high variance was located in the most lateral part of HG for the single group and in the medial part for the CSD and CPD groups, and a low variance was observed in the internal part of HG. In the stereotactic space, Leonard et al. described the probability of the
<figure><img src="image_5.png" /></figure>

**FIGURE 3** | Significant variations of Cortical Thickness (CT) and Cortical Surface Area (CSA) in the “CPD–single HG”, “CSD–single HG” and “CPD–CSD” contrasts in the left and right hemispheres showing the spatial overlap of CT and CSA variations. (A). Significant variations in cortical thickness. (B). Significant variations in cortical surface area, the lateral view of the right temporal lobe is provided to illustrate the increase in CSA located in the mid-part of the STS in the “CSD–single HG” comparison. Significant variations are superimposed on the mean patch of each corresponding group (“CPD–single” and “CPD–CSD” contrasts: mean patch of 48 left CPD, mean patch of 96 right CPD; “CSD–single HG” contrast: mean patch of 98 left CSD and 95 right CSD; CSA: cortical surface area; CT, cortical thickness; CPD, common posterior duplication; CSD, common stem duplication; HG, Heschl’s gyrus; *p* < 0.05 FDR corrected for multiple comparison; cluster extension threshold, 20 mm²; hot color indicates a positive variation; cold color, indicates a negative variation; light gray, gyri; dark gray, sulci).

position of the *sulcus intermedius* of Beck and the second Heschl's sulcus in cases with duplications. They demonstrated differences in the position of the second sulcus across hemispheres, with a second sulcus that does not angle as sharply on the left as it does on the right (Figure 6 of Leonard et al., 1998), but they did not report any difference in the variability of sulci positions across hemispheres.

## Morphological Support of CT and CSA Variations

The present work reports both local differences directly related to morphological changes associated with duplication and more distant variations that may relate to differences in anatomical and/or functional organization associated with differences in gyrification. The location of CT increases, consistently located at the second gyrus in case of CPD, and along the sulcus intermediary of Beck in case of CSD, suggests that it corresponds to a morphological change from Heschl's sulcus to the second (complete or partial) gyrus. As it had been underlined by Van Essen (1997), CT might be larger in top of gyri than in the depth of sulci.

Concerning CSA the picture was different depending on hemispheres and on the type of duplication. The modifications we report could not be attributed solely to local morphological changes in relation with differences in gyrification for HG.

## Opposite Pattern of Modification between CT and CSA differences

Interestingly, when an overlap of CT and CSA modifications occurred, it always corresponded to an increase of CT in case of duplication associated with a local decrease of CSA, an opposite profile of structural modification of these two anatomical features. These results have to be linked to the study of Hogstrom on local correlations of CT and CSA through the whole cortex. By studying 322 healthy adults aged from 20 to 85 years old thanks to surface-based method, they observed that CSA and CT were negatively correlated in a per-vertex general linear model that took into account the age and gender of the participants (Hogstrom et al., 2013), corresponding, as in the present study, to opposite differences in CT and CSA variations.

## Duplications are Associated with a Larger CT in the Second HG: A Candidate Support for Increased Abilities Associated with Duplications in Healthy Individuals

Bilateral increases in the CT could represent a morphological support for the great abilities reported in individuals having duplications. As reviewed in the introduction, expertise in pitch processing based on the spectral envelope has been associated with a right duplication (Schneider et al., 2005), while the ability to rapidly learn foreign speech sounds (Golestani et al., 2007) and phonetic expertise (Golestani et al., 2011) were related to
### CPD - single

<figure><img src="image_6.png" /><figcaption>a. Cortical thickness</figcaption></figure>

Left hemisphere

Right hemisphere

<figure><img src="image_7.png" /><figcaption>b. Cortical surface area</figcaption></figure>

Left hemisphere

Right hemisphere

### CSD - single

<figure><img src="image_8.png" /><figcaption>a. Cortical thickness</figcaption></figure>

Left hemisphere

Right hemisphere

<figure><img src="image_9.png" /><figcaption>b. Cortical surface area</figcaption></figure>

Left hemisphere

Right hemisphere

### CPD - CSD

<figure><img src="image_10.png" /><figcaption>a. Cortical surface area</figcaption></figure>

Right hemisphere

**FIGURE 4 | Morphology and relationships between variations in CT and CSA.** Projection of the SBM results on the volumetric means of each group of individuals defined according to their gyration pattern. Results of “CPD-single” and “CPD-CSD” contrasts superposed on the mean images of 48 left CPD or 96 CPD for the right side. Results of “CSD-single HG” contrasts superimposed on the mean image of 98 CSD for the left hemisphere or 95 CSD for the right hemisphere. The yellow scale corresponds to increased values, and the blue scale to decreased values. From left to right a sagittal, a coronal, and an axial slice are presented passing through the significant clusters.

a left duplication. Musicianship was also associated to a larger cortical thickness all along the posterior supra-temporal plane surface posterior to HG, thereby including the second gyrus of duplication, when one is present (Bermudez et al., 2009). Finally, the CT of the right Heschl's sulcus has been positively associated with performances in a relative pitch detection task (Foster and Zatorre, 2010). The locations across studies are consistent and the larger CT we observed in the present study overlapping with the increased CT reported by Foster's study in cases with a CSD,

whereas it was in a slightly posterior position in cases with a CPD (Foster and Zatorre, 2010). According to the literature, increase in CT could relate to an increased number and/or size of cells within a column, their packing density, or the number of connections and extent of myelination (Rakic, 1995; Eickhoff et al., 2005).

On the other hand, the increased CT associated with duplications that we report here is similar to the CT of the Heschl's sulcus label of the Destrieux parcellation (Destrieux
<table><caption>TABLE 3 | Local variations of Cortical Thickness (CT) detected in “CPD-single HG”, “CSD-single HG” and “CPD-CSD” contrasts.</caption><thead><tr><th>Anatomical localization</th><th>Cluster size (mm²)</th><th>p<sub>max</sub></th><th>x</th><th>y</th><th>z</th><th>Dupl. mean (SD) (mm)</th><th>Single HG mean (SD) (mm)</th></tr></thead><tbody><tr><td colspan="8"><strong>L: CPD-SINGLE HG</strong></td></tr><tr><td>2<sup>nd</sup> HG</td><td>143</td><td>6.10<sup>-7</sup></td><td>-48</td><td>-28</td><td>6</td><td>2.70 (0.26)</td><td>2.43 (0.25)</td></tr><tr><td colspan="8"><strong>R: CPD-SINGLE HG</strong></td></tr><tr><td>2<sup>nd</sup> HG</td><td>36</td><td>3.10<sup>-7</sup></td><td>42</td><td>-31</td><td>12</td><td>2.71 (0.34)</td><td>2.44 (0.32)</td></tr><tr><td colspan="8"><strong>L: CSD-SINGLE HG</strong></td></tr><tr><td>HS</td><td>208</td><td>3.10<sup>-9</sup></td><td>-49</td><td>-24</td><td>5</td><td>2.72 (0.26)</td><td>2.52 (0.25)</td></tr><tr><td colspan="8"><strong>R: CSD-SINGLE HG</strong></td></tr><tr><td>Medial 2<sup>nd</sup> HG-HS</td><td>145</td><td>6.10<sup>-7</sup></td><td>46</td><td>-25</td><td>8</td><td>2.80 (0.24)</td><td>2.62 (0.26)</td></tr></tbody></table>

Dupl., duplication; HG, Heschl's gyrus; HS, Heschl's sulcus; SD, Standard Deviation; *p* < 0.05 FDR corrected for multiple comparison; cluster extension threshold, 20 mm²; x, y, z, Talairach coordinates rounded to the whole number.

et al., 2010), which has been shown to be anti-correlated with the amplitude of the N1 component of auditory-evoked potentials in the right hemisphere (Liem et al., 2012). As reviewed by Liem and colleagues, the N1 amplitude reflects the first step of auditory processing and is modulated by auditory experience, such as musical training (Pantev et al., 2001; Baumann et al., 2008). Liem et al. (2012) postulated that a refinement in the synaptic circuitry might have resulted in a thinner cortex, allowing neurons to be better synchronized and generate the electrophysiological signal. This hypothesis is in accord with Hyde's proposal that a thicker cortex might be less efficient in terms of neural processing (Hyde et al., 2007). Hyde et al. reported a thicker cortex in the right inferior frontal gyrus and the right auditory cortex in amusic patients compared to musically intact controls, and suggested that they were related to cortical malformations in the amusic brain. A recent study in dyslexics also showed an increase in the CT of the right superior temporal gyrus extending to HG and the PT, which resulted in a rightward CT asymmetry that was different from the healthy volunteers (Ma et al., 2015).

Thus, the increased CT observed in cases with HG duplications may or may not represent an advantage, depending on the population studied. A thicker cortex could be associated with cognitive impairments in cases with abnormal development, while, on the other hand, it would be associated with improved abilities in the context of normal development. This hypothesis needs further exploration, but it could partially explain the apparently contradictory associations of HG duplications with both expertise and dyslexia (Leonard et al., 1993, 1995, 1998, 2001; Altarelli et al., 2014).

## The Left CSD is Associated with a Smaller CSA in the PT, While, on the Right, it is Associated with a Larger CSA in the STS

In our previous work using manual delineation of the aHG (Marie et al., 2015), we have shown that left duplication is associated with decreased left aHG surface area. The present study also shows that left duplication is associated with a decreased left PT surface area. Such an association might have been previously overlooked because the second HG is usually included in the PT in cases with a complete duplication (as reviewed by Shapleske et al., 1999). Decrease of left aHG size

had been associated with a lower leftward lateralization of the temporal processing sounds (Warrier et al., 2009), and indirectly with lateralization during speech listening (Tzourio-Mazoyer et al., 2015), while a lower left PT surface area had also been associated with lower leftward asymmetry during speech listening (Tzourio et al., 1998; Josse et al., 2003, 2006). The question of whether these decreased surface areas and related anatomical asymmetries are detrimental to speech perception in healthy subjects is an open issue; the work conducted by Leonard and others suggests, that it might be the case in dyslexics that exhibit more HG duplications (Leonard et al., 1993, 1995, 1998, 2001; Altarelli et al., 2014).

Increased CSA was located in the middle part of the STS when a right CSD occurred. STS is rightward asymmetrical for gray matter volume (Watkins et al., 2001; Barrick et al., 2005) and sulcal depth (Van Essen, 2005; Leroy et al., 2015). Human-specific rightward STS asymmetry is organized early in development, as shown by Glasel et al. (2011). The rightward asymmetries of STS described in the abovementioned works are located ventrally to the right HG, as in the present study. From a functional point of view, the right HG, and STS are two major anatomical supports of the spectral processing and pitch analysis, and they constitute the neural support of the analysis of spectral envelope independently of the fine spectro-temporal structure. Warren suggested, that STS may have a generic role in the spectral analysis of sound (Warren et al., 2005) and specified that this role did not exclude the possibility that specific subregions of the STS are functionally specialized for higher order analysis of specific sound classes. In fact, the middle STS has been involved in syllables (Jäncke et al., 2002; Poeppel et al., 2004), speech (Binder et al., 2000), and voice processing (Belin et al., 2000). Interestingly, one study showed that the acoustic parameters of voices are extracted in the right middle STS for vocal tract parameters, and in the right HG for glottal fold parameters (Kreitewolf et al., 2014). Because these two distinct anatomical areas could share a close functional role, one might suggest they have associated anatomical variations. Regarding the diversity of auditory processing mediated by the middle part of the STS, further studies are needed to conclusively determine the functional correlates of this increase in CSA in the middle STS in cases with a CSD.

Compared with a single HG in both hemispheres, CSD was also characterized by an increase of medial CSA at the interface
between the temporal transverse sulcus and aHG. These increases in the CSA could relate to the growth of the wide common stem supporting the two merging gyri in cases with a CSD.

## CPD and CSD Exhibit Anatomical Differences in the CSA in the Right aHG

The differences across duplication patterns were limited to variations in CSA on the right. This result is only partially consistent with our previous report of an absence of significant differences in the reduced surface area of aHG between the CSD and CPD groups (Marie et al., 2015), that might be in relation to the increased sensitivity provided by the examination of CSA extracted by SBM as compared to manual delineation. aHG and the PAC primarily overlap in the medial part of aHG (Penhune et al., 1996; Rademacher et al., 2001), where a decreased CSA was observed in the CPD group. This decrease may correspond to a decrease in the number/spacing of the minicolumns of the PAC, according to the radial-unit hypothesis of cortical development initially proposed by Rakic (1995, 2009) and acknowledged as a potential interpretation of such variations by others (Elmer et al., 2013; Meyer et al., 2013). The radial-unit hypothesis postulates that the morphology of the CSA of a given brain region is principally driven by the number and spacing of the ontogenetic or radial cortical columns, rather than by the number and size of cells within a column, packing density, or number of connections (Rakic, 1995; Eickhoff et al., 2005), a functional cortical column of the brain corresponding to several radial columns (Rakic, 1988, 2009). Under this hypothesis, CSA variations presently observed could relate to changes in the number and spacing of minicolumns and, consequently, to differences in processing power. Further investigations are needed to test such hypothesis and evaluate the potential behavioral correlates. Actually there is no consensus regarding the underlying cellular mechanisms associated with variations in CSA, although numerous mechanisms such as neurogenesis or myelin changes have been proposed as potential candidates for CT/gray matter, or white matter variations, respectively, (Draganski and May, 2008; Zatorre et al., 2012). This lack of knowledge applies also to gray matter density evaluated by VBM that does not correlate with histological measurements of neuronal density (Eriksson et al., 2009) pointing toward a need for understanding of the microstructural changes that underlie these macrostructural variations detected by computational neuroanatomy.

## REFERENCES

* Abdul-Kareem, I. A., and Sluming, V. (2008). Heschl gyrus and its included primary auditory cortex: structural MRI studies in healthy and diseased subjects. *J. Magn. Reson. Imaging* 28, 287-299. doi: 10.1002/jmri.21445

* Altarelli, I., Leroy, F., Monzalvo, K., Fluss, J., Billard, C., Dehaene-Lambertz, G., et al. (2014). Planum temporale asymmetry in developmental dyslexia: revisiting an old question. *Hum. Brain Mapp.* 35, 5717-5735. doi: 10.1002/hbm.22579

* Barrick, T. R., Mackay, C. E., Prima, S., Maes, F., Vandermeulen, D., Crow, T. J., et al. (2005). Automatic analysis of cerebral asymmetry: an exploratory study

Laterally, a larger right CSA was observed in the second gyrus of the CPD group compared with both the CSD group and the single HG group. This location corresponds to the PT and the second HG of CPD, which is considered to belong to the PT when present (as reviewed by Shapleske et al., 1999). On the right side, the space devoted to the supratemporal plane, which includes HG and the PT, is smaller than that on the left side (Toga and Thompson, 2003), and CPD are more frequent (Marie et al., 2015). The present results indicate, that the increased surface area of the second HG might actually be at the cost of the right PT surface area when a CPD occurs.

## CONCLUSIONS

The present SBM study conducted in 430 participants explored the changes in CT and CSA associated with HG duplications. Variations, which significance and location were independent of the template used and of the smoothing of the images, were observed in HG and surrounding areas, i.e., the left PT and the right mid STS. Despite the lack of consensus on the neurophysiological basis of such morphological features changes measured with structural MRI, the present findings provide important information for neuroimaging research led on inter-individual anatomical variability of HG in relation to speech and/or music processing, their dysfunctions, as well as training and expertise in these domains.

## AUTHOR CONTRIBUTIONS

All authors listed, have made substantial, direct and intellectual contribution to the work, and approved it for publication.

## FUNDING

This work has been funded by our guardianships: the CNRS, CEA, and University of Bordeaux and DM was supported by a grant from the Aquitaine Regional Council (France).

## SUPPLEMENTARY MATERIAL

The Supplementary Material for this article can be found online at: http://journal.frontiersin.org/article/10.3389/fnhum.2016.00069

of the relationship between brain torque and planum temporale asymmetry. *Neuroimage* 24, 678-691. doi: 10.1016/j.neuroimage.2004.09.003

* Baumann, S., Meyer, M., and Jäncke, L. (2008). Enhancement of auditory-evoked potentials in musicians reflects an influence of expertise but not selective attention. *J. Cogn. Neurosci.* 20, 2238-2249. doi: 10.1162/jocn.2008.20157

* Belin, P., Zatorre, R. J., Lafaille, P., Ahad, P., and Pike, B. (2000). Voice-selective areas in human auditory cortex. *Nature* 403, 309-312. doi: 10.1038/35002078

* Bermudez, P., Lerch, J. P., Evans, A. C., and Zatorre, R. J. (2009). Neuroanatomical correlates of musicianship as revealed by cortical thickness and voxel-based morphometry. *Cereb. Cortex* 19, 1583-1596. doi: 10.1093/cercor/bhn196
Binder, J. R., Frost, J. A., Hammeke, T. A., Bellgowan, P. S. F., Springer, J. A., Kaufman, J. N., et al. (2000). Human temporal lobe activation by speech and nonspeech sounds. *Cereb. Cortex* 10, 512–528. doi: 10.1093/cercor/10.5.512

Cohen, I. A. (1960). The coefficient of agreement for nominal scales. *Educ. Psychol. Meas.* 20, 37–46. doi: 10.1177/001316446002000104

Dale, A. M., Fischl, B., and Sereno, M. I. (1999). Cortical surface-based analysis: I. Segmentation surface reconstruction. *NeuroImage* 9, 179–194. doi: 10.1006/nimg.1998.0395

Desikan, R. S., Ségonne, F., Fischl, B., Quinn, B. T., Dickerson, B. C., Blacker, D., et al. (2006). An automated labeling system for subdividing the human cerebral cortex on MRI scans into gyral based regions of interest. *Neuroimage* 31, 968–980. doi: 10.1016/j.neuroimage.2006.01.021

Destrieux, C., Fischl, B., Dale, A., and Halgren, E. (2010). Automatic parcellation of human cortical gyri and sulci using standard anatomical nomenclature. *Neuroimage* 53, 1–15. doi: 10.1016/j.neuroimage.2010.06.010

Diallo, B., Dolidon, F., Travere, J.-M., and Mazoyer, B. (1998). VoxLine: a software for 3D real-time visualization of biomedical images. *Comput. Med. Imaging Graph.* 22, 275–289. doi: 10.1016/S0895-6111(98)00040-8

Dick, F., Tierney, A. T., Lutti, A., Josephs, O., Sereno, M. I., and Weiskopf, N. (2012). *In vivo* functional and myeloarchitectonic mapping of human primary auditory areas. *J. Neurosci.* 32, 16095–16105. doi: 10.1523/JNEUROSCI.1712-12.2012

Draganski, B., and May, A. (2008). Training-induced structural changes in the adult human brain. *Behav. Brain Res.* 192, 137–142. doi: 10.1016/j.bbr.2008.02.015

Eickhoff, S., Walters, N. B., Schleicher, A., Kril, J., Egan, G. F., Zilles, K., et al. (2005). High-resolution MRI reflects myeloarchitecture and cytoarchitecture of human cerebral cortex. *Hum. Brain Mapp.* 24, 206–215. doi: 10.1002/hbm.20082

Elmer, S., Hänggi, J., Meyer, M., and Jäncke, L. (2013). Increased cortical surface area of the left planum temporale in musicians facilitates the categorization of phonetic and temporal speech sounds. *Cortex* 49, 2812–2821. doi: 10.1016/j.cortex.2013.03.007

Eriksson, S. H., Free, S. L., Thom, M., Symms, M. R., Martinian, L., Duncan, J. S., et al. (2009). Quantitative grey matter histological measures do not correlate with grey matter probability values from *in vivo* MRI in the temporal lobe. *J. Neurosci. Methods* 181, 111–118. doi: 10.1016/j.jneumeth.2009.05.001

Fischl, B. (2012). FreeSurfer. *Neuroimage* 62, 774–781. doi: 10.1016/j.neuroimage.2012.01.021

Fischl, B., and Dale, A. M. (2000). Measuring the thickness of the human cerebral cortex from magnetic resonance images. *Proc. Natl. Acad. Sci. U.S.A.* 97, 11050–11055. doi: 10.1073/pnas.200033797

Fischl, B., Liu, A., and Dale, A. M. (2001). Automated manifold surgery: constructing geometrically accurate and topologically correct models of the human cerebral cortex. *IEEE Trans. Med. Imaging* 20, 70–80. doi: 10.1109/42.906426

Fischl, B., Salat, D. H., Busa, E., Albert, M., Dieterich, M., Haselgrove, C., et al. (2002). Whole brain segmentation: automated labeling of neuroanatomical structures in the human brain. *Neuron* 33, 341–355. doi: 10.1016/S0896-6273(02)00569-X

Fischl, B., Sereno, M. I., and Dale, A. M. (1999a). Cortical surface-based analysis. II: inflation, flattening, and a surface-based coordinate system. *NeuroImage* 9, 195–207. doi: 10.1006/nimg.1998.0396

Fischl, B., Sereno, M. I., Tootell, R. B., and Dale, A. M. (1999b). High-resolution intersubject averaging and a coordinate system for the cortical surface. *Hum. Brain Mapp.* 8, 272–284.

Fischl, B., Van Der Kouwe, A., Destrieux, C., Halgren, E., Ségonne, F., Salat, D. H., et al. (2004). Automatically parcellating the human cerebral cortex. *Cereb. Cortex* 14, 11–22. doi: 10.1093/cercor/bhg087

Foster, N. E., and Zatorre, R. J. (2010). Cortical structure predicts success in performing musical transformation judgments. *Neuroimage* 53, 26–36. doi: 10.1016/j.neuroimage.2010.06.042

Glasel, H., Leroy, F., Dubois, J., Hertz-Pannier, L., Mangin, J. F., and Dehaene-Lambertz, G. (2011). A robust cerebral asymmetry in the infant brain:

the rightward superior temporal sulcus. *Neuroimage* 58, 716–723. doi: 10.1016/j.neuroimage.2011.06.016

Glasser, M. F., and Van Essen, D. C. (2011). Mapping human cortical areas *in vivo* based on myelin content as revealed by T1- and t2-weighted MRI. *J. Neurosci.* 31, 11597–11616. doi: 10.1523/JNEUROSCI.2180-11.2011

Golestani, N., Molko, N., Dehaene, S., LeBihan, D., and Pallier, C. (2007). Brain structure predicts the learning of foreign speech sounds. *Cereb. Cortex* 17, 575–582. doi: 10.1093/cercor/bhk001

Golestani, N., Price, C. J., and Scott, S. K. (2011). Born with an ear for dialects? Structural plasticity in the expert phonetician brain. *J. Neurosci.* 31, 4213–4220. doi: 10.1523/JNEUROSCI.3891-10.2011

Hickok, G., and Poeppel, D. (2007). The cortical organization of speech processing. *Nat. Rev. Neurosci.* 8, 393–402. doi: 10.1038/nrn2113

Hogstrom, L. J., Westlye, L. T., Walhovd, K. B., and Fjell, A. M. (2013). The structure of the cerebral cortex across adult life: age-related patterns of surface area, thickness, and gyration. *Cereb. Cortex* 23, 2521–2530. doi: 10.1093/cercor/bhs231

Hyde, K. L., Lerch, J., Norton, A., Forgeard, M., Winner, E., Evans, A. C., et al. (2009). Musical training shapes structural brain development. *J. Neurosci.* 29, 3019–3025. doi: 10.1523/JNEUROSCI.5118-08.2009

Hyde, K. L., Lerch, J. P., Zatorre, R. J., Griffiths, T. D., Evans, A. C., and Peretz, I. (2007). Cortical thickness in congenital amusia: when less is better than more. *J. Neurosci.* 27, 13028–13032. doi: 10.1523/JNEUROSCI.3039-07.2007

Jäncke, L., Wustenberg, T., Scheich, H., and Heinze, H. J. (2002). Phonetic perception and the temporal cortex. *Neuroimage* 15, 733–746. doi: 10.1006/nimg.2001.1027

Josse, G., Hervé, P. Y., Crivello, F., Mazoyer, B., and Tzourio-Mazoyer, N. (2006). Hemispheric specialization for language: brain volume matters. *Brain Res.* 1068, 184–193. doi: 10.1016/j.brainres.2005.11.037

Josse, G., Mazoyer, B., Crivello, F., and Tzourio-Mazoyer, N. (2003). Left planum temporale: an anatomical marker of left hemispheric specialization for language comprehension. *Cogn. Brain Res.* 18, 1–14. doi: 10.1016/j.cogbrainres.2003.08.007

Kreitewolf, J., Gaudrain, E., and von Kriegstein, K. (2014). A neural mechanism for recognizing speech spoken by different speakers. *Neuroimage* 91, 375–385. doi: 10.1016/j.neuroimage.2014.01.005

Kulynych, J. J., Vladar, K., Jones, D. W., and Weinberger, D. R. (1993). Three-dimensional surface rendering in MRI morphometry: a study of the planum temporale. *J. Comput. Assist. Tomogr.* 17, 529–535. doi: 10.1097/00004728-199307000-00003

Leonard, C. M., Eckert, M. A., Lombardino, L. J., Oakland, T., Kranzler, J., Mohr, C. M., et al. (2001). Anatomical risk factors for phonological dyslexia. *Cereb. Cortex* 11, 148–157. doi: 10.1093/cercor/11.2.148

Leonard, C. M., Martinez, P., Weintraub, B. D., and Hauser, P. (1995). Magnetic resonance imaging of cerebral anomalies in subjects with resistance to thyroid hormone. *Am. J. Med. Genet.* 60, 238–243. doi: 10.1002/ajmg.1320600314

Leonard, C. M., Puranik, C., Kuldau, J. M., and Lombardino, L. J. (1998). Normal variation in the frequency and location of human auditory cortex landmarks. Heschl's gyrus: where is it? *Cereb. Cortex* 8, 397–406. doi: 10.1093/cercor/8.5.397

Leonard, C. M., Voeller, K. K., Lombardino, L. J., Morris, M. K., Hynd, G. W., Alexander, A. W., et al. (1993). Anomalous cerebral structure in dyslexia revealed with magnetic resonance imaging. *Arch. Neurol.* 50, 461–469. doi: 10.1001/archneur.1993.00540050013008

Leroy, F., Cai, Q., Bogart, S. L., Dubois, J., Coulon, O., Monzalvo, K., et al. (2015). New human-specific brain landmark: the depth asymmetry of superior temporal sulcus. *Proc. Natl. Acad. Sci. U.S.A.* 112, 1208–1213. doi: 10.1073/pnas.1412389112

Liem, F., Mérillat, S., Bezzola, L., Hirsiger, S., Philipp, M., Madhyastha, T., et al. (2015). Reliability and statistical power analysis of cortical and subcortical FreeSurfer metrics in a large sample of healthy elderly. *Neuroimage* 108, 95–109. doi: 10.1016/j.neuroimage.2014.12.035

Liem, F., Zaehle, T., Burkhard, A., Jäncke, L., and Meyer, M. (2012). Cortical thickness of supratemporal plane predicts auditory N1 amplitude. *Neuroreport* 23, 1026–1030. doi: 10.1097/WNR.0b013e32835abc5c

Ma, Y., Koyama, M. S., Milham, M. P., Castellanos, F. X., Quinn, B. T., Pardoe, H., et al. (2015). Cortical thickness abnormalities associated with
dyslexia, independent of remediation status. *NeuroImage Clin.* 7, 177–186. doi: 10.1016/j.nicl.2014.11.005

Marie, D., Jobard, G., Crivello, F., Perchey, G., Petit, L., Mellet, E., et al. (2015). Descriptive anatomy of heschl's gyri in 430 healthy volunteers, including 198 left-handers. *Brain Struct. Funct.* 220, 729–743. doi: 10.1007/s00429-013-0680-x

Mazoyer, B., Mellet, E., Perchey, G., Zago, L., Crivello, F., Jobard, G., et al. (2015). BIL&GIN: a neuroimaging, cognitive, behavioral, and genetic database for the study of human brain lateralization. *NeuroImage* 124, 1225–1231. doi: 10.1016/j.neuroimage.2015.02.071

Meyer, M., Liem, F., Hirsiger, S., Jäncke, L., and Hänggi, J. (2013). Cortical surface area and cortical thickness demonstrate differential structural asymmetry in auditory-related areas of the human cortex. *Cereb. Cortex* 24, 2541–2552. doi: 10.1093/cercor/bht094

Morosan, P., Rademacher, J., Schleicher, A., Amunts, K., Schormann, T., and Zilles, K. (2001). Human primary auditory cortex: cytoarchitectonic subdivisions and mapping into a spatial reference system. *Neuroimage* 13, 684–701. doi: 10.1006/nimg.2000.0715

Pantev, C., Roberts, L. E., Schulz, M., Engelien, A., and Ross, B. (2001). Timbre-specific enhancement of auditory cortical representations in musicians. *Neuroreport* 12, 169–174. doi: 10.1097/00001756-200101220-00041

Penhune, V. B., Zatorre, R. J., MacDonald, J. D., and Evans, A. C. (1996). Interhemispheric anatomical differences in human primary auditory cortex: probabilistic mapping and volume measurement from magnetic resonance scans. *Cereb. Cortex* 6, 661–672. doi: 10.1093/cercor/6.5.661

Pfeifer, R. A. (1920). *Myelogenetisch-Anatomische Untersuchungen Über das Kortikale Ende der Hörleitung*. Leipzig: BG Teubner.

Poeppel, D., Guillemin, A., Thompson, J., Fritz, J., Bavelier, D., and Braun, A. R. (2004). Auditory lexical decision, categorical perception, and FM direction discrimination differentially engage left and right auditory cortex. *Neuropsychologia* 42, 183–200. doi: 10.1016/j.neuropsychologia.2003.07.010

Rademacher, J., Bürgel, U., and Zilles, K. (2002). Stereotaxic localization, intersubject variability, and interhemispheric differences of the human auditory thalamocortical system. *Neuroimage* 17, 142–160. doi: 10.1006/nimg.2002.1178

Rademacher, J., Morosan, P., Schormann, T., Schleicher, A., Werner, C., Freund, H. J., et al. (2001). Probabilistic mapping and volume measurement of human primary auditory cortex. *Neuroimage* 13, 669–683. doi: 10.1006/nimg.2000.0714

Rakic, P. (1988). Specification of cerebral cortical areas. *Science* 241, 170–176. doi: 10.1126/science.3291116

Rakic, P. (1995). A small step for the cell, a giant leap for mankind: a hypothesis of neocortical expansion during evolution. *Trends Neurosci.* 18, 383–388. doi: 10.1016/0166-2236(95)93934-P

Rakic, P. (2009). Evolution of the neocortex: a perspective from developmental biology. *Nat. Rev. Neurosci.* 10, 724–735. doi: 10.1038/nrn2719

Ressel, V., Pallier, C., Ventura-Campos, N., Díaz, B., Roessler, A., Ávila, C., et al. (2012). An effect of bilingualism on the auditory cortex. *J. Neurosci.* 32, 16597–16601. doi: 10.1523/JNEUROSCI.1996-12.2012

Schneider, P., Scherg, M., Dosch, H. G., Specht, H. J., Gutschalk, A. and Rupp, A. (2002). Morphology of Heschl's gyrus reflects enhanced activation in the auditory cortex of musicians. *Nat. Neurosci.* 5, 688–694. doi: 10.1038/nn871

Schneider, P., Sluming, V., Roberts, N., Scherg, M., Goebel, R., Specht, H. J., et al. (2005). Structural and functional asymmetry of lateral heschl's gyrus reflects pitch perception preference. *Nat. Neurosci.* 8, 1241–1247. doi: 10.1038/nn1530

Shapleske, J., Rossell, S. L., Woodruff, P. W., and David, A. S. (1999). The planum temporale: a systematic, quantitative review of its structural, functional and clinical significance. *Brain Res. Brain Res. Rev.* 29, 26–49. doi: 10.1016/S0165-0173(98)00047-2

Sutherland, M. E., Zatorre, R. J., Watkins, K. E., Hervé, P. Y., Leonard, G., Pike, B. G., et al. (2012). Anatomical correlates of dynamic auditory processing:

relationship to literacy during early adolescence. *Neuroimage* 60, 1287–1295. doi: 10.1016/j.neuroimage.2012.01.051

Toga, A. W., and Thompson, P. M. (2003). Mapping brain asymmetry. *Nat. Rev. Neurosci.* 4, 37–48. doi: 10.1038/nrn1009

Tzourio, N., Nkanga-Ngila, B., and Mazoyer, B. (1998). Left planum temporale surface correlates with functional dominance during story listening. *Neuroreport* 9, 829–833. doi: 10.1097/00001756-199803300-00012

Tzourio-Mazoyer, N., Marie, D., Zago, L., Jobard, G., Perchey, G., Leroux, G., et al. (2015). Heschl's gyration pattern is related to speech-listening hemispheric lateralization: FMRI investigation in 281 healthy volunteers. *Brain Struct. Funct.* 220, 1585–1599. doi: 10.1007/s00429-014-0746-4

Tzourio-Mazoyer, N., Simon, G., Crivello, F., Jobard, G., Zago, L., Perchey, G., et al. (2010). Effect of familial sinistrality on planum temporale surface and brain tissue asymmetries. *Cereb. Cortex* 20, 1476–1485. doi: 10.1093/cercor/bhp209

Van Essen, D. C. (1997). A tension-based theory of morphogenesis and compact wiring in the central nervous system. *Nature* 385, 313–318. doi: 10.1038/385313a0

Van Essen, D. C. (2005). A population-average, landmark- and surface-based (PALS) atlas of human cerebral cortex. *Neuroimage* 28, 635–662. doi: 10.1016/j.neuroimage.2005.06.058

Warren, J. D., Jennings, A. R., and Griffiths, T. D. (2005). Analysis of the spectral envelope of sounds by the human brain. *Neuroimage* 24, 1052–1057. doi: 10.1016/j.neuroimage.2004.10.031

Warrier, C., Wong, P., Penhune, V., Zatorre, R., Parrish, T., Abrams, D., et al. (2009). Relating structure to function: heschl's gyrus and acoustic processing. *J. Neurosci.* 29, 61–69. doi: 10.1523/JNEUROSCI.3489-08.2009

Watkins, K. E., Paus, T., Lerch, J. P., Zijdenbos, A., Collins, D. L., Neelin, P., et al. (2001). Structural asymmetries in the human brain: a voxel-based statistical analysis of 142 mri scans. *Cereb. Cortex* 11, 868–877. doi: 10.1093/cercor/11.9.868

Wengenroth, M., Blatow, M., Heinecke, A., Reinhardt, J., Stippich, C., Hofmann, E., et al. (2013). Increased volume and function of right auditory cortex as a marker for absolute pitch. *Cereb. Cortex* 24, 1127–1137. doi: 10.1093/cercor/bhs391

Winkler, A. M., Sabuncu, M. R., Yeo, B. T., Fischl, B., Greve, D. N., Kochunov, P., et al. (2012). Measuring and comparing brain cortical surface area and other areal quantities. *Neuroimage* 61, 1428–1443. doi: 10.1016/j.neuroimage.2012.03.026

Wong, P. C., Warrier, C. M., Penhune, V. B., Roy, A. K., Sadehh, A., Parrish, T. B., et al. (2008). Volume of left heschl's gyrus and linguistic pitch learning. *Cereb. Cortex* 18, 828–836. doi: 10.1093/cercor/bhm115

Zatorre, R. J. (2013). Predispositions and plasticity in music and speech learning: neural correlates and implications. *Science* 342, 585–589. doi: 10.1126/science.1238414

Zatorre, R. J., and Belin, P. (2001). Spectral and temporal processing in human auditory cortex. *Cereb. Cortex* 11, 946–953. doi: 10.1093/cercor/11.10.946

Zatorre, R. J., Fields, R. D., and Johansen-Berg, H. (2012). Plasticity in gray and white: neuroimaging changes in brain structure during learning. *Nat. Neurosci.* 15, 528–536. doi: 10.1038/nn.3045

**Conflict of Interest Statement:** 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.

Copyright © 2016 Marie, Maingault, Crivello, Mazoyer and Tzourio-Mazoyer. 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) or licensor 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.