Handedness and its genetic influences are associated with structural asymmetries of the cerebral cortex in 31,864 individuals.

Roughly 10% of the human population is left-handed, and this rate is increased in some brain-related disorders. The neuroanatomical correlates of hand preference have remained equivocal. We resampled structural brain image data from 28,802 right-handers and 3,062 left-handers (UK Biobank population dataset) to a symmetrical surface template, and mapped asymmetries for each of 8,681 vertices across the cerebral cortex in each individual. Left-handers compared to right-handers showed average differences of surface area asymmetry within the fusiform cortex, the anterior insula, the anterior middle cingulate cortex, and the precentral cortex. Meta-analyzed functional imaging data implicated these regions in executive functions and language. Polygenic disposition to left-handedness was associated with two of these regional asymmetries, and 18 loci previously linked with left-handedness by genome-wide screening showed associations with one or more of these asymmetries. Implicated genes included six encoding microtubule-related proteins: TUBB, TUBA1B, TUBB3, TUBB4A, MAP2, and NME7-mutations in the latter can cause left to right reversal of the visceral organs. There were also two cortical regions where average thickness asymmetry was altered in left-handedness: on the postcentral gyrus and the inferior occipital cortex, functionally annotated with hand sensorimotor and visual roles. These cortical thickness asymmetries were not heritable. Heritable surface area asymmetries of language-related regions may link the etiologies of hand preference and language, whereas nonheritable asymmetries of sensorimotor cortex may manifest as consequences of hand preference.

Roughly 90% of the human population is right-handed and 10% left-handed, and this strong bias is broadly consistent across cultures, ethnicities, and history (1–5). Hand motor control is performed primarily by the contralateral brain hemisphere, such that right-handedness reflects left-hemisphere dominance for hand articulation, and vice versa (6). Behavioral precursors of hand preference are established in the developing human fetus (7, 8), probably through a genetically regulated program of asymmetrical brain development (9–12). Left-handedness occurs at increased frequencies in neurodevelopmental and psychiatric conditions, including intellectual disability (13), autism (14), and schizophrenia (15). This suggests overlapping genetic and developmental contributions to altered brain asymmetry and psychiatric conditions. However, it is important to stress that the large majority of left-handed people do not have these conditions.

Despite decades of research, the brain anatomical correlates of hand preference have remained uncertain. Some MRI studies have reported slightly altered laterality of the sensorimotor cortex around the central sulcus as well as temporal auditory cortex in left-handed people—while other studies found no effects at all (16–22). Most of these studies focused on hypothesis-driven regions of interest rather than mapping across the entire cerebral cortex, and they used different methods of analysis and measurement. The results were conflicting, not replicated or negative, possibly due to small sample sizes; the largest number of left-handers in any of these studies was 198 (16). Cortex-wide screening using atlas-defined regions has also not shown robust associations with hand preference (23, 24), even in as many as 608 left-handers versus 7,243 right-handers. One possibility is that brain anatomical correlates of handedness may be too focal to be captured by common atlas-defined parcellations. To this end, atlas-free mapping in a large population sample may provide new insights into how hand preference relates to cerebral cortical structural asymmetry.

The most statistically robust association of handedness with brain anatomy reported to date is altered average whole-hemispheric skew, or “torque,” in the horizontal and vertical planes found in 35,338 right-handed versus 3,712 left-handed adults from the UK Biobank (25). However, an association of whole-hemispheric torque with handedness does not identify specific brain regional asymmetries linked to this behavioral trait. Associations of left-handedness have also been reported with increased functional connectivity between left and right language networks in roughly 9,000 UK Biobank individuals (26) in an analysis of resting-state functional MRI data. The asymmetry of intrahemispheric functional connectivity during the resting state has also been associated with hand preference (27). The left hemisphere is dominant for language in more than 95% of the right-handed population but in only around 70% of the left-handed population, which suggests possible developmental and evolutionary relationships between these two functional asymmetries (28–32).

In the present study, we mapped cerebral cortical structural asymmetry with respect to hand preference using a large sample and an atlas-free approach: 28,802 adult right-handers and 3,062 left-handers from the UK Biobank measured for asymmetries of cortical surface area and thickness at each of 163,842 vertices in each hemisphere, before down sampling the asymmetry maps to 8,681 vertices to test associations with hand preference. Vertex-wise correspondence between the left and right hemispheres was achieved through resampling each individual’s cortical surface model to a symmetrical template created by interhemispheric coregistration (33, 34). This surface-based registration approach aligns cortical folding patterns across individuals and both hemispheres.

Using these data, we first aimed to identify specific clusters of vertices where cortical surface area or thickness asymmetries differed significantly at the group level between right and left-handers, for a uniquely well-powered and precision mapping analysis of the cortical correlates of hand preference. We then used meta-analyzed functional MRI (fMRI) data to annotate the cognitive and behavioral functions of the implicated cortical regions based on independent studies.

Handedness is a partly heritable trait with estimates of the heritability due to common genetic polymorphisms ranging from 1.2 to 5.9% across different large-scale studies and cohorts (35, 36) and around 25% in twin-based studies (37). Recent large-scale genome-wide association studies (GWAS) have made progress on the identification of genetic influences on handedness (26, 35, 36) and separately also on brain structural asymmetry (11). However, no study has previously investigated shared genetic influences on handedness and its specific cortical structural correlates. Using genome-wide genotype data for single-nucleotide polymorphisms (SNPs) in the UK Biobank individuals, we tested the SNP-based heritabilities of the cortical asymmetry measures at each regional cluster associated with hand preference, and also the genetic correlations between these regional cortical asymmetries. We then tested if polygenic disposition to left-handedness was associated with the specific cortical regional asymmetries that are associated with this behavioral trait, including the use of causal mediation analyses. In addition, we tested 39 individual SNPs in relation to the asymmetries of these cortical regions—the SNPs were previously implicated in left-handedness by a GWAS of 1,766,671 individuals (36). Together, these analyses would reveal numerous insights into the biology of gene–brain–behavioral links involving human hand preference.

Results

High-Resolution Mapping of Hand Preference in Relation to Cerebral Cortical Asymmetry in 31,864 Individuals.

For each of 28,802 right-handed and 3,062 left-handed adult individuals from the UK Biobank, we generated a cortical surface area asymmetry map and a cortical thickness asymmetry map, in which asymmetries were ultimately measured for each of 8,681 vertices across the cortical surface using an asymmetry index: AI = (Left − Right)/((Left + Right)/2) ( and see Materials and Methods).

For cortical surface area asymmetries, we found eight significant clusters (cluster-level corrected P < 0.05) where left-handers differed on average from right-handers: three on the anterior insula (most significant peak t-value = 4.98, P value = 6.55 × 10−7), one on the precentral gyrus (peak t = 4.35, P = 1.34 × 10−5), three in the anterior middle cingulate cortex (most significant peak t = 4.43, P = 9.32 × 10−6), and one on the inferior temporal/fusiform cortex (peak t = 4.74, P = 2.13 × 10−6; Fig. 1 and ). The most anterior of the insula clusters also extended onto the neighboring pars triangularis and pars orbitalis (Fig. 1 and ).

Fig. 1.

Handedness-associated cortical surface area asymmetries. (A) Group average differences of cortical surface area asymmetry between 28,802 right-handers and 3,062 left-handers shown for the whole cortex, without a threshold for significance. Red color indicates a lower asymmetry index in left-handers (weaker average leftward or stronger average rightward asymmetry in left-handers compared to right-handers). Blue indicates a higher asymmetry index in left-handers (stronger average leftward or weaker average rightward asymmetry in left-handers compared to right-handers). (B) Regional clusters with significant group average differences of surface area asymmetry between 28,802 right-handers and 3,062 left-handers. There were eight clusters that survived multiple testing correction (vertex-wise P < 0.001; cluster-wise P < 0.05). For all clusters, the effects were positive (weaker average leftward or stronger average rightward asymmetry in left-handers compared to right-handers). The clusters are shown on an inflated cortical surface model to visualize effects within folded regions. See also , which shows the same maps against Brodmann regional atlas boundaries, for reference.

For cortical thickness asymmetries, there were two clusters where left-handers differed on average from right-handers (cluster-level corrected P < 0.05): one on the postcentral gyrus (peak t = 4.05, P = 4.12 × 10−5) and one on the inferior occipital gyrus (peak t = 4.09, P = 4.39 × 10−5; Fig. 2 and ).

Fig. 2.

Handedness-associated cortical thickness asymmetries. (A) Group average differences of cortical thickness asymmetry between 28,802 right-handers and 3,062 left-handers shown for the whole cortex, without a threshold for significance. Red color indicates a lower asymmetry index in left-handers (weaker average leftward or stronger average rightward asymmetry in left-handers compared to right-handers). Blue indicates a higher asymmetry index in left-handers (stronger average leftward or weaker average rightward asymmetry in left-handers compared to right-handers). (B) Regional clusters with significant group average differences of cortical thickness asymmetry between 28,802 right-handers and 3,062 left-handers. There were two clusters that survived multiple testing correction (vertex-wise P < 0.001; cluster-wise P < 0.05). For both clusters, the effects were positive (weaker average leftward or stronger average rightward asymmetry in left-handers compared to right-handers). The clusters are shown on an inflated cortical surface model to visualize effects within folded regions. See also , which shows the same maps against Brodmann regional atlas boundaries, for reference.

Strikingly, for all surface areas and thickness asymmetries within these significant clusters, left-handers showed lower average asymmetry indexes than right-handers: in other words, if a cluster had an average rightward asymmetry in right-handers, it had a stronger average rightward asymmetry in left-handers, and if a cluster had an average leftward asymmetry in right-handers, it had a weaker average leftward asymmetry in left-handers (). These data indicate that left-handedness at the group level is associated with a relative shift of neural resources to the dominant right hemisphere for hand motor control, within all of the significant handedness-associated cortical clusters. The unilateral hemispheric group-average effects corresponding these clusters are in .

Functional Annotation of Cortical Regional Asymmetries Linked to Hand Preference.

We used meta-analyzed fMRI activation data from 14,371 studies compiled within the Neurosynth database (38) to annotate our handedness-associated structural asymmetry maps with cognitive and behavioral functions (see Materials and Methods).

The clusters with altered average surface area asymmetries in left-handedness are especially activated by tasks involving executive functions, language and reading, mood, and pain perception (Fig. 3 and ). The language-related annotations such as “word, phonological, orthographic, reading, and language” were likely contributed by the anterior insula and fusiform clusters in particular (Fig. 1 and ): clusters 2 and 8 on the anterior insula and cluster 1 on the fusiform cortex overlap with regions identified as being conjointly and left-asymmetrically activated in three sentence-level language task contrasts (39). Hemispheric language dominance is known to associate with hand preference, as atypical right-hemispheric language dominance is more frequent in left-handed than right-handed people (29, 40).

Fig. 3.

Functional annotation of regions that show altered asymmetry in left-handedness. Functional annotation based on meta-analyzed fMRI data for clusters showing altered surface area asymmetry (A) and cortical thickness asymmetry (B) in left-handedness. (Left) Coactivation maps for the handedness-associated clusters derived from the “decoder” function of Neurosynth (higher z-scores indicate greater coactivation). (Right) Cognitive terms associated with the coactivation maps for handedness-associated clusters. Only terms with correlations >0.2 between their meta-analytic activation maps and the coactivation maps for handedness-associated clusters are shown. See also .

We found that clusters with altered cortical thickness asymmetries in left-handedness are especially activated by tasks involving hand and finger movement as well as somatosensory and tactile tasks (Fig. 3 and ), likely contributed in particular by the postcentral gyrus cluster. The top functional terms were “movements, hand, execution, finger, action, motor, and motion” (Fig. 3 and ). These functional annotations based on independent fMRI data serve as a strong validation of the cortical thickness asymmetry associations with left-handedness that we observed (Fig. 3 and ). There were also numerous visual system–related functional annotations for the regions with altered cortical thickness asymmetries in left-handedness (), likely contributed especially by the inferior occipital cluster (Fig. 2 and ).

The term “motor” was annotated to the clusters showing altered surface area asymmetry and to those showing altered thickness asymmetry in left-handedness (), again supporting the validity of these anatomy–handedness associations.

Genetic Disposition to Left-Handedness Associates with Specific Cortical Asymmetries That Are Altered in Left-Handedness.

We made use of UK Biobank genome-wide genotype data to calculate the SNP-based heritabilities (41) of the cortical asymmetry measures, for each separate cluster that was associated with hand preference (Materials and Methods and ). These heritability values indicate the extent to which common genetic variation across the genome influences interindividual variation in each of the left-handedness–associated cortical regional asymmetries. The surface area asymmetries of the fusiform cluster, the three anterior insular clusters, and one anterior middle cingulate cluster were weakly but significantly heritable after false discovery rate (FDR) correction at 0.05, with pointwise heritability estimates ranging from 3.1 to 6.3% (Fig. 4 and ). The surface area asymmetry of the precentral gyrus cluster was not significantly heritable (pointwise heritability 0.0, P = 0.5). Likewise, the cortical thickness asymmetries associated with left-handedness were not significantly heritable, with pointwise estimates no greater than 1.4% (lowest P = 0.10) (Fig. 4 and ).

Fig. 4.

Genetic influences on hand preference impact some of the cortical asymmetries linked to hand preference. Five clusters of handedness-associated cortical surface area asymmetry showed significant heritabilities after FDR correction (i.e., three clusters in the anterior insular cortex [clusters 2, 4, and 8], one in the fusiform cortex [cluster 1], and one in the anterior middle cingulate cortex [cluster 6]). Other handedness-associated asymmetry clusters were not significantly heritable, including the pre- and postcentral clusters of the sensorimotor cortex (clusters 7 and 10). The associations of cortical asymmetries with polygenic disposition to left-handedness are indicated, as well as individual handedness-associated genes implicated in their variabilities.

Genetic correlation analysis indicated that the surface area asymmetries of the three anterior insula clusters were influenced by largely the same genetic variants across the genome, with significant pairwise genetic correlations from 0.85 to 1.0, whereas other pairwise genetic correlations between heritable regional asymmetries were weaker (see phenotypic and genetic correlations in ).

From an independent set of 272,673 right-handed and 33,704 left-handed individuals of the UK Biobank who were not overlapping and unrelated to those with brain image data, we derived effect sizes for the association of each SNP across the genome with hand preference (see Materials and Methods). We then applied these SNP-wise effects in mass combination to the 28,802 right-handed and 3,062 left-handed individuals that had genetic and brain image data for the present study, in order to quantify their polygenic dispositions to left-handedness, using the PRS-CS software (42). Polygenic disposition to left-handedness was significantly associated with hand preference in this sample as expected (logistic regression beta = 0.12, P = 6.63 × 10−8; see for more information), as well as cerebral cortical surface area asymmetry within the fusiform cluster (r = −0.02, P = 0.003), anterior insula cluster 2 (r = −0.02, P = 6.59 × 10−4) and anterior insula cluster 4 (r = −0.02, P = 9.01 × 10−5) after FDR correction at P < 0.05 (Fig. 4 and ). Specifically, higher polygenic disposition to left-handedness was associated with increased average rightward asymmetry in the fusiform cluster and decreased average leftward asymmetry in the anterior insula clusters (Fig. 4), consistent with the associations of these regional surface area asymmetries with the left-handedness trait itself. The polygenic disposition to left-handedness was not associated with any other regional cortical asymmetries linked to left-handedness, including the precentral surface area asymmetry cluster and the postcentral thickness asymmetry cluster (Fig. 4 and ).

Separately for the fusiform cluster and each anterior insula cluster that showed association with polygenic disposition to left-handedness, we applied mediation analyses according to two possible causal models: polygenic disposition > asymmetry > handedness and polygenic disposition > handedness > asymmetry (Materials and Methods). All of the models showed evidence of both direct and mediated effects ().

Individual Genomic Loci Associated with Hand Preference and Its Cerebral Cortical Correlates.

A previous GWAS based on 1,766,671 individuals found 41 individual genomic loci that were significantly associated with left-handedness after multiple testing correction over the whole genome (36). In the 28,802 right-handed and 3,062 left-handed individuals with both genetic and brain image data in the present study, we found that 18 of these genomic loci were associated with at least one of the cortical surface area asymmetry clusters that are linked to left-handedness after FDR correction at 0.05 (Fig. 4 and ). In particular, 17 out of these 18 SNPs were associated with surface area asymmetry of either the fusiform cluster or at least one of the anterior insula clusters (Fig. 4 and ).

The protein-coding genes located nearest to these SNPs in the genome include four that encode tubulin components of microtubules, TUBB, TUBA1B, TUBB3, and TUBB4A, as well as the microtubule-associated proteins MAP2 and NME7. These findings are in line with another recent study which implicated microtubule-related genes in brain structural asymmetry, in that case using atlas-based regional definitions across the whole brain but without considering handedness (11). Recessive mutations in NME7 cause situs inversus totalis, a rare condition in which the visceral organs (heart, liver, etc.) are reversed in their orientation on the left to right axis (43). Other genes implicated by the present study to affect both handedness and its specific cerebral cortical correlates include CNTN3, which mediates cell surface interactions during nervous system development (44), and ATXN2L, implicated in macrocephaly by de novo mutation (45).

Discussion

Cerebral Cortical Structural Asymmetry with Respect to Hand Preference.

In this study, we mapped cerebral cortical correlates of hand preference in 28,802 right-handed and 3,062 left-handed adults using an atlas-free, cortex-wide approach. This was by far the largest study of this question to date and brings clarity to a long-standing issue in human neuroscience, as previous studies had reported contradictory, equivocal, and/or negative results (see the introduction). Vertex-level mapping allowed us to detect and delineate relatively focal correlates of hand preference on the cortical surface, which may have been concealed in a previous consortium-scale study that used atlas-based regional parcellation—although the total sample size in that study was also smaller (7,243 right-handers and 608 left-handers), which would have affected the statistical power to detect subtle group-average effects (23, 46).

We found that, at the group level, left-handers had an average less leftward/more rightward shift of specific regional cortical surface area or thickness asymmetry than right-handers, distributed in several cortical regions: the fusiform, the anterior insula, the precentral, the postcentral, the anterior middle cingulate, and the inferior occipital cortex. For all of these regions, the average alterations were consistent with a relative shift of neural resources to the right hemisphere that controls the preferred hand in left-handers. Together, the implicated regions suggest a distributed neural network associated with controlling articulation of the preferred hand.

Functional annotation of the regions with altered average surface area asymmetries in left-handers, based on independent fMRI data from 14,371 meta-analyzed studies, pointed to their involvements in executive functions including working memory as well as language and reading, mood, and pain perception. Regarding language, previous studies have found an increase in atypical right-hemisphere language dominance among left-handed people (29, 40) (see the introduction). We found an association between hand preference and surface area asymmetry within the anterior insular cortex/pars triangularis corresponding to core regions of the language network that show left-lateralized activation during sentence-level language tasks (39). The anterior insula has been suggested to be especially important for the establishment of functional language lateralization during development through integrating across social, emotional, and attentional systems when acquiring language early in life (46). We found no significant associations of hand preference with asymmetries of superior temporal or middle temporal cortex, which are also core language network regions (39). Thus, hand preference apparently relates to cortical structural asymmetry of frontal and fusiform language regions in particular, rather than other language network regions that are closer to the auditory cortex of the temporal lobe.

fMRI-based annotation of the cortical regions with altered average thickness asymmetries in left-handers confirmed their involvements in hand sensorimotor functions with top annotations of “movements, hand, execution, finger,” and so on. This functional annotation based on independent data provides validation of the group-level cortical asymmetry alterations in left-handedness that we detected. As noted in the introduction, some previous studies targeted the central sulcus as a region of interest and reported altered cortical anatomy there in left-handers, but they were based on small samples and variable methods and measures with conflicting results, such that the findings remained uncertain. Here, we have implicated a specific region of the postcentral gyrus through cortex-wide, atlas-free mapping in a large sample and shown that cortical thickness asymmetry rather than surface area asymmetry is the feature most correlated with hand preference in this region. In contrast, we found that a specific region of precentral “premotor” cortex showed altered average surface area asymmetry in left-handers.

The regions that showed altered average cortical thickness asymmetry in left-handers also received numerous vision-related annotations based on fMRI data, such as “‘motor imagery, visual motion, eye movements, visual perception, eye fields,” and so on. These annotations were likely driven by the inferior occipital cluster. Left-handers have been reported to show more variable hemispheric involvement than right-handers in visual half-field tasks, which can index lateralized contributions to visual attentional tasks (47). The need for hand–eye coordination during complex manual tasks may drive this relationship. However, the inferior occipital cluster also lies within a language network region that shows left-lateralized activation during sentence-level tasks presented in different modalities (39).

Shared Genetic Influences on Hand Preference and Its Specific Cerebral Cortical Correlates.

We found that polygenic disposition to left-handedness was associated with specific handedness-linked asymmetries of surface area in the fusiform and anterior insular cortex. The directions of effect were as expected: higher polygenic disposition to left-handedness–associated with reduced average leftward or increased average rightward surface area asymmetry in these cortical regions—the same as for left-handedness itself. Notably, the regional cortical thickness asymmetries associated with hand preference (i.e., within the postcentral gyrus and inferior occipital cortex) were not significantly heritable and showed no associations with polygenic disposition to left-handedness. Neither did the surface area asymmetry of the precentral gyrus cluster. Therefore, some specific handedness-related cortical asymmetries share genetic influences with hand preference, while others apparently do not, at least to extents that could be detected in the current sample size.

In the adult cross-sectional data of the UK Biobank, it is not possible to establish cause–effect relations unambiguously. Mediation models were compatible with fusiform and anterior insula surface area asymmetries partly mediating the effect of polygenic disposition on hand preference, but also the other way around, i.e. hand preference partly mediating an effect of left-handedness polygenic disposition on these brain asymmetries. Both models are plausible, as brain anatomy may in principle be both causal to behavioral outcomes and shaped plastically by repeated behaviors, or indeed merely correlated because of shared underlying factors (48). However, the cortical asymmetries that share genetic influences with hand preference seem more likely to be involved in its etiology than those which do not, especially as some of the implicated genes may point to mechanisms of establishing the brain’s left to right axis in early development (see discussion of specific genes below). Furthermore, structural asymmetry of the anterior insula—the region with the strongest evidence for shared genetic influences on both its asymmetry and hand preference—is already detectable in the fetal brain (49).

In contrast, nonheritable correlates of hand preference in the adult cortex, including the sensorimotor cortex, may reflect plastic adaptation after the preference is already established developmentally. The lack of heritability of these asymmetries in the adult brain could then reflect their relatively distal relationships to genetically regulated mechanisms of left to right brain axis formation in the embryo. However, left-handedness has shown only limited associations with environmental and early life factors studied so far (50, 51) and is also only weakly to moderately heritable (35–37). This suggests that left-handedness arises largely from random (nonheritable) early developmental variation (52) occurring around a genetically regulated program that is biased toward a right-handed outcome. In other words, the genetically regulated program of asymmetrical brain development appears to give rise to right-handedness with roughly 90% probability and left-handedness with roughly 10% probability, while both interindividual genetic and environmental variation seem to have generally small impacts on these baseline probabilities. In this way, left-handedness might have some brain structural and functional correlates linked to its early development which are neither heritable nor affected by environmental variation.

This study pinpointed individual genomic loci which associate with both hand preference and some of its specific cortical structural correlates. The involvement of microtubule-related genes in human handedness and brain asymmetry suggests a possible mechanism whereby left to right asymmetry of the brain arises in the embryo from cellular chirality (an asymmetrical twisting of cellular shape specified by the cytoskeleton—internal cellular structure composed of protein filaments and microtubules) (11). Such mechanisms have been described with respect to other organs of other invertebrate and vertebrate species and can be organ intrinsic (i.e., arise independently from other developing organs or systems) (53–59). This seems especially pertinent to the brain, as left-hemisphere dominances for hand preference and language do not usually reverse in people with situs inversus totalis (i.e., reversal of the visceral organs on the left-right axis) when caused by genetic mutations which affect primary ciliary components (60–63). Primary cilia play an important role in the establishment of visceral asymmetry, in which the molecular chirality of the axoneme (microtubule-based central strand of the cilium) results in unidirectional ciliary motion within the embryonic node (64, 65).

In this context, it is especially notable that our data implicated NME7 in affecting both hand preference and handedness-associated cortical surface area asymmetries. Recessive mutations in this gene cause situs inversus totalis but without primary ciliary dyskinesia (i.e., ciliary function remains largely or wholly unaffected even while the visceral organs become reversed) (43). The functions of this gene are not well understood, but it is known to associate with γ-tubulin and the ciliary axoneme (43); thus, it has potential functions that involve the cilium but also the microtubule cytoskeleton more generally. NME7 provides a possible genetic link between brain and body asymmetry via as-yet-unidentified mechanisms which do not necessarily involve ciliary dysfunction and the nodal pathway of visceral laterality formation (60). All of the genes implicated here to affect both hand preference and its cerebral cortical correlates should be studied with gene-functional approaches, such as in cellular models and knockout mice, to understand their roles in affecting brain functional and structural laterality.

Limitations.

Vertex-wise mapping of asymmetry assumes that high-resolution correspondence can be achieved between contralateral vertices in the two hemispheres. For the most part the method appears to have worked well, insofar as it distinguished some regions of association with hand preference that make sense on a priori grounds, such as the sensorimotor cortex of the pre- and postcentral sulcus and language-related regions. However, we found it necessary to exclude a region surrounding the boundary with subcortical structures at the center of the brain (affecting entorhinal, parahippocampal, and cingulate regions; ), where extreme vertex-wise asymmetry values were present in many individuals (see Materials and Methods). This may indicate an artifact that involves extreme warping for this region to fit the symmetrical template.

The limitation of the UK Biobank cross-sectional data for making causal inferences from associations was already mentioned in the section Shared Genetic Influences on Hand Preference and Its Specific Cerebral Cortical Correlates. In addition, the UK Biobank is a volunteer cohort of middle-aged to older adults and therefore not fully representative of the population, with generally superior health and socioeconomic circumstances for this age range, and greater female participation (66). Higher risk for psychiatric disorders has been suggested to reduce participation in cohort studies (67), and we have previously found evidence for this in the UK Biobank (68).

In this study, we linked genetic, brain structural, and behavioral data and made use of independent, meta-analyzed fMRI data to understand the functions of the implicated sets of cortical regions. Future studies of hand preference in large datasets may seek to integrate more aspects of brain structure and function from the same individuals into the analyses, including task-based or resting-state fMRI and diffusion tensor imaging data on white matter structural connectivity.

The assessment of hand preference was by a simple questionnaire in which individuals could choose between four options: right-handed, left-handed, use both hands equally, or prefer not to answer. We found the “both hands” category to be relatively unreliable over repeat visits (Materials and Methods), which led us to use a binary left/right variable for the purposes of this study. Simple assessments such as this have been shown to capture an inherent dichotomy in hand preference that is also revealed by more quantitative, multi-item questionnaires (69). However, brain–behavior association results may have been different if using performance-based hand skill measures or semiquantitative multi-item handedness ratings, and ambidextrous people may form a somewhat distinct category with respect to behavioral and brain correlates (70, 71).

This study examined average group differences in cerebral cortical asymmetry between left- and right-handers. Future studies may consider the possibility that left-handedness associates with increased variability in terms of brain structure and function.

Summary.

This study used a symmetrical surface template to achieve atlas-free mapping of cerebral cortical asymmetry with respect to hand preference in an unprecedented sample size. Detailed and statistically reliable maps of the cortical correlates of human hand preference were produced, resolving a long-standing issue in human neuroscience. The implicated regions are especially involved in executive, language, motor, and visual functions. Across all the regions, the data uniformly indicated a less leftward/more rightward shift of neural resources in left-handers compared to right-handers, consistent with right-hemispheric control of the left hand. Hand preference was genetically associated with the asymmetries of frontal and fusiform language network regions, supporting developmental and evolutionary links between hand preference and language. Various other regional asymmetries associated with hand preference were not heritable and may reflect plastic changes in support of hand preference, including those of the primary sensorimotor cortex. Specific genes found to affect both hand preference and its cerebral cortical correlates included microtubule-related genes, one of which—NME7—provides a potential mechanistic link between brain and visceral asymmetries, which might not involve impaired ciliary function. This study provides an example of how linking across human genetics, brain structure, and behavior in large datasets can provide tangible new knowledge, as well as mechanistic and developmental hypotheses for future research on typical and atypical human brain function.

Materials and Methods

Participants.

This study was conducted under UK Biobank application 16066, with Clyde Francks as principal investigator. The UK Biobank is a general adult population cohort (72, 73). The UK Biobank received ethical approval from the National Research Ethics Service Committee North West - Haydock (reference 11/NW/0382), and all of their procedures were performed in accordance with the World Medical Association guidelines. Informed consent was obtained for all participants.

Imaging genetics dataset.

We used the UK Biobank neuroimaging data (74, 75) released in February 2020 together with the genotype data from the same participants. To achieve a sufficiently high degree of homogeneity for the genetic analyses, we restricted the sample to participants with “White British ancestry” defined by Bycroft et al. (“in.white.British.ancestry.subset”) (72). We excluded subjects with a mismatch of their self-reported and genetically inferred sex, with putative sex chromosome aneuploidies, or who were outliers based on heterozygosity (principal component corrected heterozygosity > 0.19) and genotype missingness (missing rate > 0.05) (72). We randomly excluded one subject from each pair with a kinship coefficient >0.0442 as defined by the UK Biobank (72).

Handedness was assessed by a touchscreen question with four choices: right-handed, left-handed, use both right and left hands equally, and prefer not to answer (UK Biobank field: 1707). We previously found that those indicating mixed handedness (roughly 2% of all individuals) had a high rate (41%) of changing their answers over repeat visits to the assessment centers, whereas those indicating either right-handedness or left-handedness at their first visit had consistency rates >97% on subsequent visits (31, 50). We therefore excluded mixed handers and focused only on those indicating either right- or left-handedness during their first visit.

These steps resulted in a final sample of 31,864 participants comprising 28,802 right-handers and 3,062 left-handers. The age range of these participants was from 45 to 81 y (mean 63.75), and 15,064 were male and 16,800 female.

Genome-wide screening dataset to obtain SNP-wise effect sizes for association with hand preference.

These additional UK Biobank individuals had genetic and handedness data but not brain image data. We applied the same genetic and handedness exclusion criteria as in the section Imaging genetics dataset with the extra requirement that none of these individuals should have a kinship coefficient >0.0442 with any of the individuals in the imaging genetic dataset. This yielded 272,673 right-handed and 33,704 left-handed participants not overlapping and unrelated to the set with brain imaging data.

MRI Processing.

Brain image data were derived from T1-weighted MRI scans (Siemens Skyra 3 Tesla MRI with 32-channel radio frequency receive head coil). We started from the Freesurfer 6.0 (76) “recon-all” cortical reconstructions generated by the UK Biobank imaging team (UK Biobank data field 20263, first imaging visit) but did not make use of image-derived phenotypes released by that team. Rather, the cortical thickness and surface area maps of both hemispheres of each individual were resampled into a standard symmetrical space created using interhemispheric coregistration (fsaverage_sym) (33), with 163,842 vertices per hemisphere (). This process required roughly 10 wk of processing on 12 cluster server nodes running in parallel and precisely aligned morphometry features between the two hemispheres to achieve vertex-wise correspondence. Separately for surface area and thickness maps, asymmetry for each left–right-matched pair of vertices in each participant was calculated according to the asymmetry index AI = (left − right)/((left + right)/2). The asymmetry maps for surface area and thickness of each individual were then down sampled to the left-hemisphere of template “fsaverage5” (10,242 vertices per hemisphere; ).

We noticed that many individuals had extreme values of asymmetry in a region wrapping around the boundary with subcortical structures near the center of the brain (), which likely arises from a registration artifact. For this region, we excluded any vertices from all individuals when they had surface area |AI|>1 in at least 1,000 (∼3%) individuals () (i.e., excluded from subsequent analyses of both surface area and cortical thickness asymmetries). This left 8,681 vertex-wise asymmetry measures per individual spanning the majority of the cerebral cortex.

Mapping Cortical Asymmetry Correlates of Hand Preference.

For each vertex, brain asymmetry differences between left- and right-handers were examined by two-sample t tests separately for surface area and cortical thickness while controlling for continuous covariates: age when attended assessment center (field 21003–2.0), quadratic age (age − mean_age)2, scanner position parameters (X, Y, and Z: fields 25756–2.0, 25757–2.0, and 25758–2.0), T1 signal-to-noise ratio (field 25734–2.0), and T1 contrast-to-noise ratio (field 25735–2.0) together with two categorical covariates: assessment center (field 54–2.0) and sex (field 31–0.0). Significant clusters were determined by the random field theory method at a cluster-corrected level of P < 0.05 (clusters forming with vertex-level P < 0.001) to minimize the likelihood of false positive results using the SurfStat toolbox (https://www.math.mcgill.ca/keith/surfstat/). Each cluster identified in this way has a peak vertex and additional contiguous vertices, with cluster-wise significance assessed against a random field with the same spatial correlation as the real map.

As a post hoc analysis, we also examined the unilateral effects corresponding to clusters with altered asymmetry in left-handers. Specifically, for each hemisphere and individual, we calculated the mean unilateral value across all vertices within each significant cluster and compared the differences between left- and right-handers by two-sample t tests while controlling for the same covariates mentioned above.

Functional Annotation of Cortical Regions Showing Altered Asymmetry in Left-Handers.

We made use of the Neurosynth database (38), a platform for large-scale automated synthesis of task-based fMRI data. This database defines brain-wide activation maps corresponding to specific cognitive or behavioral task terms using meta-analyzed functional activation maps from studies referring to those terms (38). At the time of accessing the database, it included 1,307 maps representing activation patterns corresponding to 1,307 cognitive or behavioral terms from 14,371 human neuroimaging studies. We created one bilateral mask in MNI152 standard space by labeling all clusters showing altered average surface area asymmetry in left-handers, and another bilateral mask in Montreal Neurological Institute (MNI) standard space that comprised all clusters showing altered average cortical thickness asymmetry in left-handers. This was achieved through applying the “mri_surf2vol” function in Freesurfer to the down-sampled, handedness-associated asymmetry cluster maps (defined in left hemisphere fsaverage5 surface space) and then identifying the mirrored voxels (based on the midline of left and right hemisphere) corresponding to these clusters in the right of the symmetrical standard space, to define masks with bilaterally labeled voxels.

These masks were then used separately as input to identify relevant cognitive and behavioral terms through the “decoder” function of the Neurosynth database. This tool first creates a coactivation map for each input mask based on all activation data in the database and then tests each term-specific activation map for spatial correlation with the coactivation map. Cognitive terms with correlations >0.2 are reported while excluding anatomical terms, nonspecific terms (e.g., “Tasks”), and one from each pair of virtually duplicated terms (such as “Words” and “Word”). This tool does not employ inferential statistical testing but rather ranks cognitive and behavioral terms according to the correlations of their meta-analyzed activation maps with the coactivation map that corresponds to a user-defined input mask.

Heritabilities and Genetic Correlations of Cortical Asymmetries Associated with Hand Preference.

From the imputed SNP genotype data released by the UK Biobank (March 2018), 9,516,135 autosomal variants with minor allele frequencies > 1%, INFO (imputation quality) score > 0.7 and Hardy–Weinberg equilibrium P > 1 × 10−7 were used to build a genetic relationship matrix using the Genome-wide Complex Trait Analysis (GCTA) software (77) (version 1.93.0beta). Specifically for analyses using GCTA, we further excluded one random participant from each pair having a kinship coefficient higher than 0.025 based on the calculated genetic relationship matrix (this analysis is particularly sensitive to higher levels of relatedness), resulting in 29,973 participants for this particular analysis (27,083 right-handed and 2,890 left-handed). Genome-based restricted maximum likelihood (GREML) analysis using GCTA (77) was used to estimate the SNP-based heritability for each cluster that showed altered asymmetry in left-handedness (using the mean asymmetry index across all vertices within each cluster), controlling for the abovementioned covariates [i.e., age, (age–mean_age)2, scanner position parameters, T1 signal-to-noise ratio, T1 contrast-to-noise ratio, assessment center, and sex] plus an additional binary covariate for genotyping array and the first ten genetic principal components capturing population genetic diversity (UK Biobank field ID: 22009–0.1∼22009–0.10) (72) and then applying FDR 0.05 across the 10 clusters. Bivariate GREML analysis (78) was used to estimate genetic correlations between pairs of significantly heritable cluster-wise asymmetries.

Polygenic Disposition to Left-Handedness.

In 272,673 right-handers and 33,704 left-handers who were not overlapping and unrelated to those with brain image data (see Genome-wide screening dataset to obtain SNP-wise effect sizes for association with hand preference), we again excluded genetic variants with minor allele frequencies < 1%, INFO score ≤ 0.7, and Hardy–Weinberg equilibrium P value ≤ 10−7. We then carried out an autosome-wide association scan for right- versus left-handedness under an additive genetic model with covariates for year of birth, sex, genotyping array, and the first 10 genetic principal components capturing population genetic diversity, using BGENIE (v1.2) (72). SNP-wise effects on left-handedness disposition from this sample were then applied in mass combination to the 28,802 right-handers and 3,062 left-handers with post-quality-control imaging and genetic data (see Imaging genetics dataset) to estimate their polygenic dispositions to left-handedness, using the PRS-CS software (42). This method uses a high-dimensional Bayesian regression framework to estimate posterior effect sizes of SNPs and has shown superior prediction statistics compared to approaches based on clumping by linkage disequilibrium and P value thresholding (79). In total, effect sizes of 1,103,636 SNPs spanning the autosomes were used based also on their presence in the 1,000 Genomes project European-descent dataset from where PRS-CS obtains its linkagedisequilibrium information (80). We used default parameters and the recommended global effect size shrinkage parameter phi = 0.01. The handedness polygenic scores in the imaging genetic dataset were z-standardized for subsequent analyses ().

We first tested whether polygenic disposition to left-handedness was associated with hand preference in the 28,802 right-handers and 3,062 left-handers using logistic regression with hand preference as the outcome variable and covariates of age, (age–mean_age)2, genotyping array, the first 10 genetic principal components (see Heritabilities and Genetic Correlations of Cortical Asymmetries Associated with Hand Preference), and sex.

We then extracted the mean asymmetry index over all vertices within each separate cluster that showed association of its asymmetry with left-handedness (eight clusters for surface area asymmetry, two for thickness asymmetry). Separately for each cluster’s mean asymmetry, we performed Pearson correlation analysis with polygenic disposition to left-handedness across individuals while controlling for the same covariates as in the SNP-based heritability analysis (see Heritabilities and Genetic Correlations of Cortical Asymmetries Associated with Hand Preference) and applying FDR correction at 0.05 across the 10 clusters tested.

For the asymmetry clusters that were significantly associated with polygenic disposition to left-handedness, we then applied the “mediation” toolbox in R https://www.rdocumentation.org/packages/mediation/versions/4.5.0 using the same covariates and bootstrapping 1,000 times in order to test causal mediation models framed with respect to either polygenic disposition > asymmetry > handedness or polygenic disposition > handedness > asymmetry separately for each cluster’s mean asymmetry.

Imaging Genetic Analysis of Individual Loci Associated with Left-Handedness.

A previous genome-wide association scan for left- versus right-handedness carried out in 1,766,671 individuals reported 41 individually significant genomic loci after correction for multiple testing across the whole genome (36). There were genotype data for 38 of these SNPs in the 28,802 right-handers and 3,062 left-handers of the imaging genetic sample in the present study, and we identified a proxy SNP for one more of these SNPs [rs11265393 as a proxy for rs66513715 with linkage disequilibrium r-squared of 1 in European-descent data according to LDproxy (81)]. We tested each of these 39 SNPs for association with each of the significantly heritable, cortical cluster-wise asymmetries (averaged over all vertices separately within each cluster) that showed alterations in left-handers using an additive model in BGENIE (v1.2) (72) with the same covariates as the heritability analysis (see Heritabilities and Genetic Correlations of Cortical Asymmetries Associated with Hand Preference). All 39 SNPs had imputation quality (INFO) scores > 0.7. We applied the FDR correction at 0.05 to control multiple testing over all of these SNP-cluster combinations.

Note that the previous GWAS of left- versus right-handedness based on 1,766,671 individuals (36) included the UK Biobank participants and therefore also included the participants with brain image data, although the brain image data were not used in that study. The sample overlap meant that we could not use those genome-wide, SNP-wise effect size estimates to calculate polygenic disposition to left-handedness in the current study, as the genome-wide screening and target datasets must not overlap for polygenic scoring analysis. In addition, the SNP-wise effect sizes have not been made publicly available from the GWAS of 1,766,671 individuals (36).