Regional Grey Matter Structure Differences between Transsexuals and Healthy Controls—A Voxel Based Morphometry Study

Abstract

Gender identity disorder (GID) refers to transsexual individuals who feel that their assigned biological gender is incongruent with their gender identity and this cannot be explained by any physical intersex condition. There is growing scientific interest in the last decades in studying the neuroanatomy and brain functions of transsexual individuals to better understand both the neuroanatomical features of transsexualism and the background of gender identity. So far, results are inconclusive but in general, transsexualism has been associated with a distinct neuroanatomical pattern. Studies mainly focused on male to female (MTF) transsexuals and there is scarcity of data acquired on female to male (FTM) transsexuals. Thus, our aim was to analyze structural MRI data with voxel based morphometry (VBM) obtained from both FTM and MTF transsexuals (n = 17) and compare them to the data of 18 age matched healthy control subjects (both males and females). We found differences in the regional grey matter (GM) structure of transsexual compared with control subjects, independent from their biological gender, in the cerebellum, the left angular gyrus and in the left inferior parietal lobule. Additionally, our findings showed that in several brain areas, regarding their GM volume, transsexual subjects did not differ significantly from controls sharing their gender identity but were different from those sharing their biological gender (areas in the left and right precentral gyri, the left postcentral gyrus, the left posterior cingulate, precuneus and calcarinus, the right cuneus, the right fusiform, lingual, middle and inferior occipital, and inferior temporal gyri). These results support the notion that structural brain differences exist between transsexual and healthy control subjects and that majority of these structural differences are dependent on the biological gender.

About the Authors:

Lajos Simon

Affiliation: Semmelweis University, Department of Psychiatry and Psychotherapy, Budapest, Hungary

Lajos R. Kozák

Affiliation: Semmelweis University, Magnetic Resonance Imaging Research Center, Budapest, Hungary

Viktória Simon

* E-mail: [email protected]

Affiliation: Semmelweis University, Department of Psychiatry and Psychotherapy, Budapest, Hungary

Pál Czobor

Affiliations Semmelweis University, Department of Psychiatry and Psychotherapy, Budapest, Hungary, Nathan Kline Institute for Psychiatric Research, Orangeburg, New York, United States

Zsolt Unoka

Affiliation: Semmelweis University, Department of Psychiatry and Psychotherapy, Budapest, Hungary

Ádám Szabó

Affiliation: Semmelweis University, Magnetic Resonance Imaging Research Center, Budapest, Hungary

Gábor Csukly

Affiliation: Semmelweis University, Department of Psychiatry and Psychotherapy, Budapest, Hungary

Introduction

Transsexualism is a heterogenous condition both in its manifestation and etiology. There are numerous chromosomal abnormalities or well-defined biological causes that can lie behind the incongruence between an individual's biological gender (i.e. that based on the sex chromosomes and/or the manifestation of the sexual organs) and gender identity (i.e., someone's sense and perception of being male or female). However, there is a group of individuals who do not show such known genetic or somatic abnormality and yet experience strong incongruency between their assigned biological gender and their gender identity. In the Diagnostic and Statistic Manual of Mental Disorders 4th edition - Text Revision [1] this is referred as Gender Identity Disorder (GID). GID is characterized by a long-standing and strong feeling of being a member of another gender, long-standing distress about the assigned gender and feeling of incongruity with the assigned gender-roles causing significant clinical discomfort and impairment in both the individual's social and professional life and in other life areas. Based on DSM-IV-TR, GID cannot be diagnosed if transgender experiences are associated with any physical intersex conditions.

The background of transsexualism has been the topic of debate for decades. Recently, mainly early developmental disturbances have been suggested by the pertinent literature [2]. According to a recent review about the sexual differentiation of the human brain, transsexualism might be the result of the fact that the development of the sexual organs in the fetal life occurs well before the sexual differentiation of the brain. Thus, if something disturbs the sexual differentiation of the brain, the fetus already has sexual organs according to his/her assigned sex, while his/her brain might develop differently [2]. These authors suggest that the disturbance of the testosterone surge that masculinize the fetal brain might be at the background of GID in certain cases. Furthermore, they emphasize that there is no compelling evidence that postnatal environmental factors play a crucial role in sexual orientation and gender identity [2]. The theories about the neurobiological background of GID were partly based on earlier neuroanatomical findings of the same group, Swaab and collegues. This group described differences in the brain structure of Male to Female (MTF) and Female to Male (FTM) transgender subjects and controls (post mortem) in regions of the brain that showed sex differences regarding their size [3]. Specifically, they found that the size and the number of neurons in the bed nucleus of striata terminalis (BSTc) and the third interstitial nucleus of the anterior hypothalamus (INAH-3) of MTF transgender patients were typical for the size and neuron numbers found generally in females [3]–[5].

In the last two decades, however, structural imaging studies reported controversial results about structural brain differences in transgender and control subjects. Four structural MRI studies focused on white matter differences in transgender patients. Emory et al. [6] found no difference in the whole corpus callosum or the splenium region between transsexuals and controls. More recently, Yokota et al. [7] concluded that the pattern of corpus callosum shape in transsexuals was closer to that in subjects with the same gender identity than to that in subjects with the same biological sex. Rametti et al. conducted two separate studies in untreated FTM and MTF patients, applying diffusion tensor imaging (DTI) [8], [9]. They found that in case of FTM transsexuals certain fasciculi involved in higher cognitive functions the white matter microstructure pattern was closer to the pattern of subjects who shared their gender identity (males) than to those who shared their biological gender (females) [8]. In case of MTF transsexuals, DTI findings showed that the white matter microstructure pattern was intermediate between male and female controls. Based on the direction of the differences the authors suggested that certain fasciculi did not complete the masculinization process during brain development [9]. Rametti et al. (2012) [10] further investigated FTM transsexuals in order to examine the effect of testosterone treatment on brain structure. They found that testosterone treatment changed white matter microstructure in FTMs. Specifically, the fractional anisotropy (FA) values in the superior longitudinal fascicule and the right corticospinal tract increased after 7 months of testosterone treatment, compared with the FA values before treatment. Furthermore, this effect was predicted by the pretreatment, free testosterone levels [10].

Four structural MRI studies investigated the structure of grey matter (GM) and yielded in mixed results. One study, using voxel based morphometry (VBM) concluded that the right putamen was “feminized” in MTF transsexuals, another study applying VBM and MR volumetry concluded that their data did not support the hypotheses according to which MTF transsexuals show feminized brain structure in certain areas, however reported on particular features of the brain structure of nonhomosexual MTF transsexuals including areas that are involved in body perception [11], [12]. The third study, measuring GM thickness reported thicker GM in MTF transsexuals compared with control males in several brain areas both in the left and the right hemisphere, along both the medial and lateral surface of the cortex [13]. This study concluded that their results support the notion that gender identity and brain anatomy are associated and that a “shift” exists between MTF transsexuals and gender congruent males with regard to brain structure [13]. The most recent study in this field investigated cortical thickness and volumetric differences of certain subcortical regions in both MTF and FTM transsexuals by MR volumetry. This study found that FTMs showed subcortical grey matter masculinization (right putamen), while MTFs showed feminization with regard to cortical thickness, as well as right hemisphere localized differences compared with male controls [14].

Except for the most recent study of Zubiaurre-Elorza et al (2012) [14], all of the aforementioned three studies obtained their data only on MTF transsexuals. Thus, there is scarcity of imaging data on FTM transsexuals. Our objective was to compare the regional cortical structure of both untreated MTF and FTM transgender subjects with that of male and female control subjects, applying VBM.

Our hypothesis was that the regional structural parameters of the brain of transsexual subjects will be different from that of control subjects with the same biological gender.

Methods

This investigation was carried out in accordance with the latest version of the Declaration of Helsinki. The study was approved by the local ethics committee (Semmelweis University Regional and Institutional Committee of Scientific and Research Ethics) and all included subjects provided written informed consent.

Subjects

All consecutively arriving patients of the Transgender Special Outpatient Service of the Psychiatry and Psychotherapy Department of the Semmelweis University (Budapest, Hungary) diagnosed with GID, based on the DSM-IV TR diagnostic criteria [1] were approached to enter a neuroimaging study, between 2007 September and 2009 March. Both MTF and FTM patients were eligible for the study, but only those with homosexual orientation. The rationale for this choice was based on the Blanchard typology [15] which considers two fundamentally different types of transsexualism: homosexual and nonhomosexual. Homosexual transsexual individuals are sexually attracted to the same biological gender, while nonhomosexual transsexual individuals are attracted to either the opposite gender or show no sexual orientation/attraction at all. According to Blanchard, homosexual transsexuals are usually younger at initial presentation of gender identity disorder and show more pronounced and frequent childhood femininity, as well as different anthropometric data [16], [17]. One might argue that mixing individuals from both transsexual groups in one study targeting the neurobiological background of transsexualism might bias the results by introducing heterogeneity in the sample. Thus, in our study, only homosexual transsexual individuals were included preventing our findings from the aforementioned bias. Evaluation of the subjects' sexual orientation was based on self-report.

Further exclusion criteria were: previous cross-hormonal treatment, any known chromosomal or hormonal disorder in the background of transgender identity, any neurological disorder in the anamnesis.

Healthy volunteers were recruited to serve as controls from among medical students, colleagues and friends of the research team who were free from any symptoms of GID or any psychiatric disorders. The presence of symptoms of GID was evaluated based on a free clinical interview asking simple questions targeting the symptoms of GID listed in DSM-IV-TR, while the presence of psychiatric symptoms was assessed by SCL-90 [18]. Control subjects were selected to represent a population matched in age and gender identity to the patient group.

Only data from the structural imaging findings are presented in this paper, results of the functional imaging findings will be reported in upcoming publications.

Diagnosis of GID

All GID subjects underwent a detailed diagnostic interview with an expert psychiatrist in the field and also filled out a test battery assessing transgender identity disorder symptoms and associated behaviors and psychiatric comorbidity in order to confirm the diagnosis and exclude the presence of other mental disorder behind the symptoms of GID. Sexual orientation of the patients was assessed by self-report. During the clinical interview basic demographic data, family history, psychiatric history and psychiatric status were also collected.

2.3 MRI image acquisition

High resolution anatomical images of all participants were collected at the MR Research Center (Semmelweis University, Budapest, Hungary) on a 3 Tesla Philips Achieva whole body clinical MRI scanner (Philips, Best, The Netherlands) equipped with an 8-channel SENSE head-coil. Whole brain anatomical images were obtained using a T1 weighted three dimensional spoiled gradient echo (T1W 3D TFE) sequence as provided by the vendor, but fine-tuned to provide the best possible separation between white matter and GM; 180 contiguous slices were acquired from each subject with the following imaging parameters: TR = 9.7 ms; TE = 4.6 ms; flip angle  = 8°; FOV of 240 mm×240 mm; voxel size of 1.0×1.0×1.0 mm.

Data processing

We performed voxel based morphometry analysis [19] on the imaging data to investigate differences in local GM volume between our subject groups. Data preprocessing and analysis were performed within the SPM8 software framework (http://www.fil.ion.ucl.ac.uk/spm/) using the VBM8 toolbox (http://dbm.neuro.uni-jena.de/vbm/). We applied the default processing parameters of the VBM8 toolbox using the default preprocessing pipeline consisting of the following steps: (a) segmentation of the different tissue classes, (b) linear (i.e., affine) and nonlinear (i.e., Dartel) registration [20] of the subjects' brains to the Montreal Neurological Institute (MNI) atlas space, [21] with 1.5×1.5×1.5 mm resolution using the ICBM152 and DARTEL templates supplied with the SPM8 toolbox and (c) modulation of the GM tissue segments by the nonlinear normalization parameters to account for individual brain size differences. The segmentation procedure was refined by accounting for partial volume effects [22], applying adaptive maximum a posteriori estimations, and denoising using a hidden Markov random field model [23], and also by using a spatially adaptive non-local means filters [24]. Next, the whole data set was evaluated for outliers by checking sample homogeneity using covariances between image pairs; none of the images were excluded from further analysis based on this check. Finally, the realigned and normalized GM segments were smoothed with an 8 mm FWHM Gaussian kernel.

Statistical analysis

Voxel-wise image intensities (henceforth GM intensities) of the smoothed warped modulated GM compartments representing regional GM volume were compared using a 2x2 ANCOVA model specified in SPM8, with Biological gender, and the presence of GID as main effects and age as a covariate of no interest. Upon whole brain model estimation, F contrasts were calculated to assess the sites of main effects and interactions (GID, Biological gender and GID × Biological gender). Statistical maps were considered significant at the level of p<0.001 uncorrected, due to the limited statistical power we present unthresholded F-statistical maps along with thresholded maps, to show that the significant clusters are indeed markedly different from the background [25]. To limit the possibility that local image artifacts or noise affect the outcome of statistical comparisons a cluster size threshold of 30 voxels was set; this way only clusters with a minimum volume of 101.25 cubic millimeters were considered relevant. Anatomical labeling of significant clusters was performed using the xjView toolbox (http://www.alivelearn.net/xjview).

Post hoc analyses were performed on significant clusters for the Biological gender, GID, and GID × Biological gender interaction contrasts. To this end, GM intensity values corresponding to the statistically significant clusters were extracted from the individual brains, and averaged within subject using the MarsBaR SPM toolbox (http://marsbar.sourceforge.net/) [26] and processed further in Matlab as follows: in case of the Biological gender and GID contrasts post hoc F and T statistics were calculated; for these comparisons differences were considered significant at a level of p<0.05. In case of clusters of the GID × Biological gender contrast a 2×2 ANOVA was performed, followed by all possible pairwise comparisons of the control females, control males, FTM transsexuals and MTF transsexuals, using unpaired two-sample Student's tests between the groups. The significance levels for these post hoc Student's tests were corrected for multiple comparisons using the Bonferroni method, i.e., only results with p<0.05/6 = 0.008 (as six comparisons were performed for a given cluster), were considered as significant.

Results

Basic demographic data

Seventeen patients with the diagnosis of GID (10 MTF patients and 7 FTM patients), based on the DSM-IV TR [1] diagnostic criteria, before cross-hormonal treatment, and 18 age matched control subjects (11 females and 7 males) were included in the study. Mean age of subjects with GID was: 28.5 years (SD:7.69) for males and 24.8 years (SD:6.45) for females, while mean age of controls was 27.1 years (SD: 5.54) for males and 23.9 years (SD:3.42) for females.

3.2 Results of whole brain VBM

The whole brain analysis showed 3 significant clusters for the GID effect, i.e. where the regional GM volume was found to be significantly different between patients and controls (Table 1 and Figure 1). The Biological gender main effect, accounting for gender differences but not for GID status was significant in 6 clusters (Table 2). Significant GID × Biological gender interactions were found in 4 clusters (Table 3, and Figure 2).

[Figure omitted. See PDF.]

Figure 1. Main effect of GID.

The figure represents clusters with significant GM volume difference, depending on GID status. The results are based on a group-wise 2×2 ANCOVA model, estimated upon the whole brain, where GID status and Biological gender were the main factors and age served as a covariate of no interest. Differences were considered significant at p<0.001, uncorrected with a cluster size threshold of 101.25 cubic millimeters (left panels). In the right panels the unthresholded SPMF maps are shown. Color coding of clusters in the left panels and maps in the right panels is based on F-values, and is similar across clusters; the p<0.001 threshold is marked with a dotted line on the color bar. GID: Gender Identity Disorder; GM: grey matter

https://doi.org/10.1371/journal.pone.0083947.g001

[Figure omitted. See PDF.]

Figure 2. The sites of GID × Biological gender interaction.

Clusters with significant GID × Biological gender interaction are shown, based on a 2×2 ANCOVA model where age served as a covariate of no interest. Interactions were considered significant upon whole brain model estimation at p<0.001 uncorrected with a cluster size threshold of 101.25 cubic millimeters (left panels). In the upper rows, those clusters are visible in which GM volume was larger in subjects with male gender identity, while the lower rows depict clusters in which GM volume was larger in subjects with female gender identity. In the right panels, the unthresholded SPMF maps are shown. Color coding of clusters in the left panels and maps in the right panels is based on F-values, and is similar across clusters; the p<0.001 threshold is marked with a dotted line on the color bar. GID: Gender Identity Disorder; GEN: Biological gender; GM: Grey matter

https://doi.org/10.1371/journal.pone.0083947.g002

[Figure omitted. See PDF.]

Table 1. Significant effect of GID on regional cortical structurea.

https://doi.org/10.1371/journal.pone.0083947.t001

[Figure omitted. See PDF.]

Table 2. Significant effects of biological gender on regional cortical structurea.

https://doi.org/10.1371/journal.pone.0083947.t002

[Figure omitted. See PDF.]

Table 3. Regions of interest with significant GID × Biological gender interactiona.

https://doi.org/10.1371/journal.pone.0083947.t003

Results of the post hoc analyses

Results of the post hoc analyses performed on the clusters for the GID effect indicated that local GM volume in healthy controls was higher in all significant regions compared with GID patients (229 voxels affecting the left cerebellum, 76 voxels affecting the right cerebellum, and a 32 voxel cluster including the left angular gyrus and the left inferior parietal lobule; Table 1).

In significant clusters of the Biological gender main effect, higher local GM volume was found in males compared to females in 3 of the clusters (239 voxels affecting the Right and Left Posterior Cingulate and Precuneus; 51 voxels affecting the Left Parahippocampal Gyrus, Left Cerebellum Anterior Lobe, Left Culmen; 33 voxels affecting the Right Inferior Occipital Gyrus, Right Lingual Gyrus), while the opposite direction of effect was found in the other clusters (44 voxels affecting the Right Superior Temporal Gyrus; 35 voxels affecting the Right Brainstem and Pons; 31 voxels affecting the Left Superior Temporal Gyrus, and Left Insula) (Table 2).

The post hoc analyses of the significant clusters with GID × Biological gender interaction showed that transgender patients had significantly different regional GM volume from that of controls sharing their biological gender and did not differ significantly from that of controls sharing their gender identity in the examined clusters (Figures 3).

[Figure omitted. See PDF.]

Figure 3. Post hoc analysis for GID × Biological gender contrast.

Results of the post hoc analysis of the GID × Biological gender ANCOVA contrasts are shown. For this purpose voxel-wise GM intensity values (representing local GM volume), corresponding to the statistically significant clusters were extracted from the individual brains, then averaged within subject and submitted to a 2×2 ANOVA followed by all possible pair-wise comparisons of the control females, control males, FTM transsexuals and MTF transsexuals, using unpaired two-sample Student's tests between the groups. Each panel demonstrates the post hoc analysis results for an individual cluster (using the same labeling scheme as in Figure 2 previously). White bars show the local grey matter volume for biological females and black bars show the local grey matter volume for biological males. Dashed lines indicate significance level at p<0.001, while solid lines indicate significance level at p<0.0083 (Bonferroni corrected). GID: Gender Identity Disorder; FTM: female to male transsexual subject; MTF: male to female transsexual subject

https://doi.org/10.1371/journal.pone.0083947.g003

Specifically, in a 304 voxel sized cluster affecting the left pre- and postcentral gyri (including the somatosensory cortex and the primer motor cortex); a 32 voxel sized cluster affecting the left posterior cingulate, calcarine gyrus, and the precuneus, showed lower GM volume in MTF transgender patients and female controls compared to FTM transgender patients and male controls. The opposite direction of differences could be observed in a 123 voxel sized cluster in the right occipital lobe involving the middle and inferior occipital, the fusiform, and the lingual gyri, in a 42 voxel sized cluster affecting the right inferior temporal gyrus, where regional GM volume proved to be higher in MTF transgender patients and female controls compared with FTM transgender patients and male controls (please see Table 3 and Figure 3 for details).

Discussion

Early studies investigating the brain structure of transgender patients were based on post mortem neuroanatomical methodology [3]–[5], but recently in vivo imaging studies gained importance. While the findings are diverse, the general trend of the results points towards a distributed pattern of neuroanatomical/structural differences subtending multiple brain areas of transgender subjects compared with controls from the same biological gender [6]–[14]. It has to be noted that direct comparison of the results of these studies is hampered by the different methodological approaches used.

We chose VBM as our method of investigation, and to our knowledge the present study is unique with respect to examining both FTM and MTF transsexuals with this approach. Hence, we were able to fully explore GID × Biological gender interactions, and we identified 4 clusters where local GM volume of transgender patients did not differ significantly from controls sharing their gender identity but were different from those sharing their biological gender. Specifically, in the area of the right middle and inferior occipital gyri, the fusiform, and the lingual gyri, the right inferior temporal gyrus, regional GM volume was higher in those subjects who had female gender identity (MTF transsexuals and female controls), while in the left pre-and postcentral gyri, the left posterior cingulate, calcarine gyrus and precuneus, regional GM volume was higher in those subjects with male gender identity (FTM transsexuals and male controls).

Previous studies that used VBM approach in a fashion somewhat similar to our study were limited to investigating the MTF type of GID. For example, Luders et al. [11] compared MTF transsexuals before cross-hormonal treatment (n = 24) to male (n = 30) and female controls (n = 30). Their VBM study showed that only the right putamen had significantly larger volume in MTF transsexuals compared with male control subjects, and was within the range found in females. Luders et al applied a relatively conservative approach when restricted their statistical maps to clusters ≥123 voxels, which, together with the fact that they examined MTF transgender subjects only, makes it difficult to directly compare our results. Nevertheless, there were only 2 brain areas that showed significant GID × Biological gender interaction and 1 brain area that showed significant GID effect at the cluster size limit of ≥123 voxels in our study: a cluster involving the left pre- and postcentral gyri (304 voxels) and a cluster (123 voxels) involving the right middle and inferior occipital gyri, the fusiform and lingual gyri and a cluster (229 voxels) involving the left cerebellum, respectively.

In another study also limited to MTF transsexuals Savic and Arver, [12] reported no “feminization” of any brain region with regard to structure. Nonetheless, certain brain areas (clusters ≥100 voxels) showed characteristic structural features in the transsexual group compared with both male and female control groups. Specifically, they found reduced thalamus and putamen volumes and increased GM volumes in the insular and inferior frontal cortex and in the right temporo-parietal junction (angular gyrus and superior temporal gyrus) in the transsexual group compared with both control groups. In our study, however only the angular gyrus (but in the left hemisphere) was affected among these areas, showing lower regional GM concentration in both FTM and MTF transgender subjects compared to controls, independent of their biological gender. When comparing the results reported by Savic and Arver to either our study or to other imaging studies in the literature of transsexualism, it has to be taken into consideration that their reported results were obtained from a solely nonhomosexual transsexual group of patients. The lack of real overlap between our and Savic and Arvers' findings, despite the very similar methodology used, might at least in part be explained by the difference of the sexual orientation of the two samples.

In a recent paper, the GM thickness of MTF transsexuals (n = 24) and age matched male controls (n = 24) was compared [13] and significant between group differences were reported in several areas including the right pre- and postcentral gyri, the left postcentral- and paracentral gyri, the right precuneus, the right inferior temporal, lingual and fusiform gyri [13]. In all of the reported regions, MTF transsexuals showed thicker GM compared to male controls and no areas were reported where controls showed thicker GM compared to MTF transsexuals [13]. Although they examined only biologically male subjects, while our study examined both male and female subjects, there is partial overlap between the affected brain areas found in our study and that of Luders et al [13]. Specifically, in both studies areas in the right inferior temporal, lingual, and fusiform gyri were affected and showed thicker/more concentrated GM in MTF transsexuals compared with male control subjects. Zubiaurre-Elorza et al's findings [14] are also consistent with these latter findings from both our and Luders et al's [13] study. These authors investigated both MTF (n = 18) and FTM transsexuals (n = 24) with MR volumetry, and compared the results with male (n = 29) and female control (n = 23) subjects [14]. They found that FTM transsexuals showed subcortical grey matter masculinization in the right putamen, while MTF transsexuals showed feminized cortical thickness data. Specifically, MTFs did not differ from female controls with regard to cortical thickness but differed from male controls in the regions of the lateral and medial orbito-frontal regions, the insula, and the medial occipital cortex of the right hemisphere [14]. Nonetheless, a caveat is needed when comparing cortical volumetry results with VBM because the latter approach can estimate either regional GM concentration or regional GM volume depending on the processing parameters; moreover due to normalization and modulation steps VBM represents both the volume and the shape of the local cortical sheet.

We found a robust difference between transsexuals and controls affecting both the anterior and posterior lobe of the left cerebellum (229 voxels), an effect that is independent from the biological gender. Such a difference in the cerebellar region has not been documented in the literature.

The regions found affected in our study are mainly involved in neural networks playing role in body perception, including memory retrieval, self-awareness, visual processing, body and face recognition and sensorimotor functions [27]–[30]. Findings from the studies of Savic and Arver [12] and Luders et al. [13] also implicated that these functional networks are affected. We must note, however, that it is still unclear how the observed structural differences might translate into functional differences. Furthermore, it is also unclear how the observed differences could be interpreted in terms of a causal chain: as parts of the neurobiological background or perhaps as the consequences of transsexualism [12].

Methodological strengths and limitations

Our reported results have to be interpreted in light of certain limitations. Most importantly, given the limited availability of the patient population, and our strict exclusion criteria, the number of GID patients who finally underwent imaging was 17, which might be considered low, despite that this sample size is comparable to similar imaging studies in the literature [8]–[14]. The limited sample size may indeed limit statistical power, however it helps to pinpoint clear trends in the data that can be validated by further investigations.

We used VBM for investigating in vivo, regional structural brain differences between subjects with GID and controls. VBM is considered to be an objective and reliable measure for comparing local structural differences voxelwise, by pinpointing changes in the local composition of the cortex after discounting macroscopic differences in brain shape by means of normalization [19]. Nevertheless, VBM is hampered by certain methodological shortcomings that make it difficult to compare results across studies. First, the intensity data captured and analyzed with the VBM method cannot easily be translated into the “size” of a particular brain area, since it is rather a mixture of local volume and shape; moreover the results of segmentation heavily depend on data quality [31]. Second, the choice of processing parameters (e.g. whether to use standard or optimized VBM, whether to use modulation, whether to clean up GM boundaries, etc.), and the choice of statistical model can heavily influence results and interpretation, e.g. VBM on unmodulated GM compartments represents local grey matter concentration, while VBM on GM compartment modulated by the Jacobian normalization transformation represents the absolute volume [19], [31]. Interpreting VBM results may also be problematic, as it is often not easy to explain connection between grey matter volume and brain function.

Although there are no widely adopted standards in the VBM literature, we followed the guidelines recommended by Ridgway et al [25]. Due to the limited sample size, an exploratory whole brain analysis was performed first at p<0.001 uncorrected significance level, followed by post hoc-analyses at more stringent statistical criteria. We adopted this strategy since, on one hand, presenting uncorrected results does not necessarily limit the validity of the findings in case of VBM, as shown by Ashburner and Friston [19], who concluded that the validity of statistical tests based on uncorrected statistics was not severely compromised. On the other hand, in order to show that the significant clusters are indeed markedly different from the background, we presented unthresholded F-statistical maps, as well [25]. The post-hoc analyses were performed at the cluster level, hence being less susceptible to the nonstationarity of the smoothness of VBM data [19]. If a difference was found to be significant after the within-cluster averaging and the Bonferroni-corrected post-hoc comparisons, then it was supposed to be stronger than just a trend-effect.

The necessity and the proper way of cluster size thresholding is also an unsolved issue in the VBM literature. Based on the fact that the spatially variant smoothness distribution of VBM data limits the validity of such thresholding approaches [32], Ashburner and Friston stated that voxel-based extent statistic should not be used in VBM [19]. Nevertheless, most authors tend to utilize cluster size thresholding albeit with very different, and most often arbitrary parameters [33], [34]. The reason behind the popularity of such thresholding is that it can decrease the chance of Type I errors, even at lower than usual statistical thresholds, while it can also limit the graininess of the results. Based on our relatively low sample size, we decided to use cluster size thresholding. We set the threshold to 30 voxels, which is slightly more than 100 cubic millimeters. This choice is arbitrary, but it is in the range of previous papers [33], [34] and provides an easy to interpret size scale for the reader.

According to the pertinent literature, age is an important factor that might play a role in the inconsistent results of neuroimaging studies with regard to sex dimorphism [35], [36], i.e. age related changes might confound results. In order to minimize the age associated bias of our findings we included age matched controls in our sample which was originally from a young cohort with a narrow age range (17–38 years).

FTM and MTF transgender patients were both included in this study, which is rare in the literature. Including both MTF and FTM subjects in our analyses allowed us to differentiate between findings that are connected to the “patient” status (GID effect) from those findings that are connected to a GID × Biological gender interaction (i.e. that would suggest sex dimorphic changes). This statistical analytic approach is also unique in the pertinent literature.

Conclusion

Our findings support the notion that structural differences exist between subjects with GID and controls from the same biological gender. We found that transsexual subjects did not differ significantly from controls sharing their gender identity but were different from those sharing their biological gender in their regional GM volume of several brain areas, including the left and right precentral gyri, the left postcentral gyrus (including the somatosensory cortex and the primary motor cortex), the left posterior cingulate, precueneus and calcarinus, the right cuneus, the right fusiform, lingual, middle and inferior occipital, and inferior temporal gyri. Additionaly, we also found areas in the cerebellum and in the left angular gyrus and left inferior parietal lobule that showed significant structural difference between transgender subjects and controls, independent from their biological gender.

There is only small number of studies in the field of structural imaging of transgender subjects and although the applied methods are different, sample sizes are modest and results are therefore yet inconclusive in details, significant structural differences were found between transgender patients and controls repeatedly and in some cases, in overlapping brain areas. These initial results, including the results of our study, need to be further replicated and refined in future studies on larger samples, as well as followed by functional imaging studies that might clarify how these structural differences impact the process of the disturbed evolution of gender identity and/or how disturbed gender identity affect brain structure and functions.

Author Contributions

Conceived and designed the experiments: LS GCS ZSU PC LRK. Performed the experiments: LS GCS LRK ASZ. Analyzed the data: LRK ASZ VS. Contributed reagents/materials/analysis tools: LS LRK PC GCS ZSU. Wrote the paper: VS LRK PC GCS LS ZSU.

Citation: Simon L, Kozák LR, Simon V, Czobor P, Unoka Z, Szabó Á, et al. (2013) Regional Grey Matter Structure Differences between Transsexuals and Healthy Controls—A Voxel Based Morphometry Study. PLoS ONE 8(12): e83947. https://doi.org/10.1371/journal.pone.0083947

References

1. American Psychiatric Association (2000) Diagnostic and Statistical Manual of Mental Disorders, Fourth Edition, Text Revision. APA, Washington, DC.

2. Bao AM, Swaab DF (2011) Sexual differentiation of the human brain: Relation to gender identity, sexual orientation and neuropsychiatric disorders. Front Neuroendocrinol 32: 214–226.

3. Zhou JN, Hofman MA, Gooren LJ, Swaab DF (1995) A sex difference in the human brain and its relation to transsexuality. Nature 378: 68–70.

4. Garcia-Falgueras A, Swaab DF (2008) A sex difference in the hypothalamic uncinate nucleus: relationship to gender identity. Brain 131: 3132–3146.

5. Kruijver FP, Zhou JN, Pool CW, Hofman MA, Gooren LJ, et al. (2000) Male-to-female transsexuals have female neuron numbers in a limbic nucleus. J Clin Endocrinol Metab 85: 2034–2041.

6. Emory LE, Williams DH, Cole CM, Amparo EG, Meyer WJ (1991) Anatomic variation of the corpus callosum in persons with gender dysphoria. Arch Sex Behav 20: 409e17.

7. Yokota Y, Kawamura Y, Kameya Y (2005) Callosal shapes at the midsagittal plane: MRI differences of normal males, normal females, and GID. In: Proceedings of the 2005 IEEE, Engineering in Medicine and biology 27th Annual Conference; p.3055e8..

8. Rametti G, Carrillo B, Gomez-Gil E, Junque C, Segovia S, et al. (2011a) White matter microstructure in female to male transsexuals before cross-sex hormonal treatment. A diffusion tensor imaging study. J Psychiatr Res 45(2): 199–204.

9. Rametti G, Carrillo B, Gómez-Gil E, Junque C, Zubiaurre-Elorza L, et al. (2011b) The microstructure of white matter in male to female transsexuals before cross-sex hormonal treatment. A DTI study. J Psychiatr Res 45(7): 949–54.

10. Rametti G, Carillo B, Gómez-Gil E, Junque C, Zubiaurre-Elorza L, et al. (2012) Effects of androgenization on the white matter microstructure of female-to-male transsexuals. A diffusion tensor imaging study. Psychoneuroendocrinol 37: 1261–1269.

11. Luders E, Sanchez FJ, Gaser C, Toga AW, Narr KL, et al. (2009) Regional gray matter variation in male-to-female transsexualism. Neuroimage 46(4): 904–907.

12. Savic I, Arver S (2011) Sex Dimorphism of the Brain in Male-to-Female Transsexuals. Cerebral Cortex 21: 2525–2533.

13. Luders E, Sánchez FJ, Tosun D, Shattuck DW, Gaser C, et al. (2012) Increased Cortical Thickness in Male-to-Female Transsexualism. J Behav Brain Sci 2: 357–362.

14. Zubiaurre-Elorza L, Junque C, Gómez-Gil E, Segovia S, Carillo B, et al.. (2012) Cortical Thickness in Untreated Transsexuals. Cereb Cortex doi:10.1093/cercor/bhs267

15. Blanchard R, Clemmensen LH, Steiner BW (1987) Heterosexual and homosexual gender dysphoria. Arch Sex Behav 16: 139–152.

16. Blanchard R (1989) The classification and labeling of nonhomosexual gender dysphorias. Arch Sex Behav 18: 315–334.

17. Blanchard R, Dickey R, Jones CL (1995) Comparison of height and weight in homosexual versus nonhomosexual male gender dysphorics. Arch Sex Behav 24: 543–554.

18. Derogatis LR, Rickels K, Rock AF (1976) The SCL-90 and the MMPI: a step in the validation of a new self-report scale. Br J Psychiatry 128: 280–289.

19. Ashburner J, Friston KJ (2000) Voxel-based morphometry-The methods. Neuroimage 11: 805–821.

20. Ashburner J (2007) A fast diffeomorphic image registration algorithm. Neuroimage 38(1): 95–113.

21. Evans AC, Collins DL, Mills SR, Brown ED, Kelly RL, et al.. (1993) 3D statistical neuroanatomical models from 305 MRI volumes. Proc. IEEE Nucl. Sci. Symp. Med. Imaging Conf. 1813–1817.

22. Tohka J, Zijdenbos A, Evans A (2004) Fast and robust parameter estimation for statistical partial volume models in brain MRI. Neuroimage 23: 84–97.

23. Rajapakse JC, Giedd JN, Rapoport JL (1997) Statistical Approach to Segmentation of Single-Channel Cerebral MR Images. IEEE Trans. Med. Imag 16: 176–186.

24. Manjon JV, Coupé P, Martí-Bonmatí L, Collins DL, Robles M (2010) Adaptive Non-Local Means Denoising of MR Images With Spatially Varying Noise Levels. J Magn Reson Imaging 31: 192–203.

25. Ridgway GR, Henley SM, Rohrer JD, Scahill RI, Warren JD, et al. (2008) Ten simple rules for reporting voxel-based morphometry studies. NeuroImage 40: 1429–1435.

26. Brett M, Anton J-L, Valabregue R, Poline J-B (2002) Region of interest analysis using an SPM toolbox. 8th International Conference on Functional Mapping of the Human Brain. 2002 June 2-6, Sendai, Japan

27. Gallagher HL, Frith CD (2003) Functional imaging of “theory of mind”. Trends Cogn Sci 7: 77–83.

28. Hodzic A, Kaas A, Muckli L, Stirn A, Singer W (2009) Distinct cortical networks for the detection and identification of human body. Neuroimage 45: 1264–1271.

29. Northoff G, Heinzel A, de Greck M, Bermpohl F, Dobrowolny H, et al. (2006) Self-referential processing in our brain-a meta-analysis of imaging studies on the self. Neuroimage 31: 440–457.

30. Rekkas PV, Constable RT (2005) Evidence that autobiographic memory retrieval does not become independent of the hippocampus: an fMRI study contrasting very recent with remote events. J Cogn Neurosci 17: 1950–1961.

31. Mechelli A, Price CJ, Friston KJ, Ashburner J (2005) Voxel-Based Morphometry of the Human Brain: Methods and Applications. Curr Med Imaging Rev 1 2pp 105–113.

32. Woermann FG, Free SL, Koepp MJ, Ashburner J, Duncan JD (1999) Voxel-by-voxel comparison of automatically segmented cerebral grey matter-A rater-independent comparison of structural MRI in patients with epilepsy. NeuroImage 10: 373–384.

33. Henley SM, Ridgway GR, Scahill RI, Klöppel S, Tabrizi SJ, et al. (2010) Pitfalls in the use of voxel-based morphometry as a biomarker: examples from huntington disease. Am J Neuroradiol 31: 711–719.

34. Ferreira LK, Diniz BS, Forlenza OV, Busatto GF, Zanetti MV (2011) Neurostructural predictors of Alzheimer's disease: a meta-analysis of VBM studies. Neurobiol Aging 32: 1733–1741.

35. Chen X, Sachdev PS, Wen W, Ansteyc KJ (2007) Sex differences in regional GM in healthy individuals aged 44-48 years: A voxel-based morphometric study. Neuroimage 36: 691–699.

36. Raz N, Gunning-Dixon F, Head D, Rodrigue KM, Williamson A, et al. (2004) Aging, sexual dimorphism, and hemispheric asymmetry of the cerebral cortex: replicability of regional differences in volume. Neurobiol Aging 25: 377–396.