Introduction
A higher prevalence of mixed and left (non-right) handedness has been reported in preterm (<37 weeks’ gestation) and low birthweight (<2500 g) children compared with term born peers (Marlow, Roberts, & Cooke, Reference Marlow, Roberts and Cooke1989; Powls, Botting, Cooke, & Marlow, Reference Powls, Botting, Cooke and Marlow1996). Findings from a recent meta-analysis exploring handedness in preterm children found that, based on 18 case–control studies with a total of 1947 cases and 8170 controls, the rate of non-right handedness (NRH) in the preterm population was approximately 22%, compared with 12% in term children (Domellof, Johansson, & Ronnqvist, Reference Domellof, Johansson and Ronnqvist2011).
NRH has been associated with atypical development and poorer outcomes in preterm and low birthweight populations (Domellof et al., Reference Domellof, Johansson and Ronnqvist2011); however, past research has been limited and conflicting. Some studies have reported NRH to be associated with lower IQ, language delays, poorer motor coordination, and more difficult temperament (O’Callaghan, Burn, Mohay, Rogers, & Tudehope, Reference O’Callaghan, Burn, Mohay, Rogers and Tudehope1993; Ross, Lipper, & Auld, Reference Ross, Lipper and Auld1987), with stronger right-handed preference related to better cognitive development (Marlow, Hennessy, Bracewell, & Wolke, Reference Marlow, Hennessy, Bracewell and Wolke2007). Other studies, however, have reported similar cognitive, academic, and ball and balance skill performances between non-right and right-handed very low and extremely low birthweight children (Powls et al., Reference Powls, Botting, Cooke and Marlow1996; Saigal, Rosenbaum, Szatmari, & Hoult, Reference Saigal, Rosenbaum, Szatmari and Hoult1992), with little evidence that hand preference is related to the presence or severity of cerebral palsy or neurological classification (Ross et al., Reference Ross, Lipper and Auld1987). These inconsistent findings may reflect differences in measures of handedness and age of participants (4–6 years olds vs. 8–12 year olds), with studies involving older children finding little association. Hand preference may be more influenced by social, cultural, or educational effects in older compared with younger children (Medland, Perelle, De Monte, & Ehrman, Reference Medland, Perelle, De Monte and Ehrman2004). Furthermore, previous studies report on cohorts born several decades ago and may not be representative of more contemporary cohorts of very preterm (VPT) children in which survival is greater for those born at the limits of viability. This pattern of atypical development for NRH is not characteristically seen in the general population (Somers, Shields, Boks, Kahn, & Sommer, Reference Somers, Shields, Boks, Kahn and Sommer2015). However, the literature remains mixed even within the general population, with some studies reporting lower levels of achievement and cognitive performance for left-handed individuals (Johnston, Nicholls, Shah, & Shields, Reference Johnston, Nicholls, Shah and Shields2009; Nicholls, Chapman, Loetscher, & Grimshaw, Reference Nicholls, Chapman, Loetscher and Grimshaw2010; Resch et al., Reference Resch, Haffner, Parzer, Pfueller, Strehlow and Zerahn-Hartung1997).
Various theories have been proposed to explain the increased NRH prevalence and potential associated cognitive and motor deficits observed in preterm individuals, including genetic or hereditary factors. While the contribution of genetics and familial handedness to hand preference has mostly been explored in the general population (Annett, Reference Annett1985; McManus, Reference McManus1991; Medland et al., Reference Medland, Duffy, Wright, Geffen, Hay, Levy and Boomsma2009), these factors may influence the process of lateralization in VPT populations alike. One study maintained a higher incidence of NRH in preterm compared with term-born children even after accounting for familial handedness (Ross et al., Reference Ross, Lipper and Auld1987). However, further research is required to investigate the extent of the contribution of heritability to NRH in VPT children.
Currently, the prominent theory for increased NRH in preterm individuals is that NRH may reflect a disturbance of cerebral lateralization and underlying neurological abnormality (Bakan, Dibb, & Reed, Reference Bakan, Dibb and Reed1973; Coren & Porac, Reference Coren and Porac1980; Satz, Orsini, Saslow, & Henry, Reference Satz, Orsini, Saslow and Henry1985). However, previous findings have been inconsistent, with some studies illustrating an association between neurological abnormalities (abnormal neurological examination or cerebral ultrasound) and hand preference (Allin et al., Reference Allin, Rooney, Griffiths, Cuddy, Wyatt, Rifkin and Murray2006; Saigal et al., Reference Saigal, Rosenbaum, Szatmari and Hoult1992), and others failing to detect evidence of such a relationship (Marlow et al., Reference Marlow, Roberts and Cooke1989; Powls et al., Reference Powls, Botting, Cooke and Marlow1996).
It may be that the hemisphere involved in focal neurological damage is important to consider for NRH in preterm children, with one small study finding a trend for greater left-sided intraventricular hemorrhage in non-right handed preterm children based on cranial ultrasound (Powls et al., Reference Powls, Botting, Cooke and Marlow1996). Neonatal magnetic resonance imaging (MRI) is a sensitive diagnostic tool for detecting white matter punctate lesions, cerebellar injury, or subcortical insults in preterm children (Inder, Anderson, Spencer, Wells, & Volpe, Reference Inder, Anderson, Spencer, Wells and Volpe2003; Kidokoro et al., Reference Kidokoro, Anderson, Doyle, Woodward, Neil and Inder2014); however, the relationship between these pathologies and handedness has not yet been explored. In older children, MRI studies (Caldu et al., Reference Caldu, Narberhaus, Junque, Gimenez, Vendrell, Bargallo and Botet2006; Kesler et al., Reference Kesler, Vohr, Schneider, Katz, Makuch, Reiss and Ment2006; Lancefield et al., Reference Lancefield, Nosarti, Rifkin, Allin, Sham and Murray2006; Nosarti et al., Reference Nosarti, Rushe, Woodruff, Stewart, Rifkin and Murray2004), have found no association of handedness with fronto-occipital asymmetry (Lancefield et al., Reference Lancefield, Nosarti, Rifkin, Allin, Sham and Murray2006), temporal gyral width (Kesler et al., Reference Kesler, Vohr, Schneider, Katz, Makuch, Reiss and Ment2006), or corpus callosum size (Nosarti et al., Reference Nosarti, Rushe, Woodruff, Stewart, Rifkin and Murray2004) in VPT cohorts.
Given the possible association between neurological abnormality and NRH, usage of diffusion weighted imaging of implicated motor tracts, such as the corpus callosum (CC) and corticospinal tract (CST), may help clarify the relationship between handedness and cerebral morphology. Diffusion weighted imaging, and more specifically diffusion tensor imaging, enables examination of the microstructural properties of cerebral white matter tracts, providing several measures indicative of tract integrity, maturity, and organization. Specifically, fractional anisotropy (FA) is a measure of directional water diffusion and higher FA is indicative of greater white matter organization. The apparent diffusion coefficient, or mean diffusivity (MD), is an overall measure of diffusion, with axial diffusivity (AD) being a measure of diffusion parallel to white matter tracts, and radial diffusivity (RD) being an estimate of diffusion perpendicular to white matter tracts. Lower diffusivity is associated with greater maturation of white matter tracts and diffusivity is known to decrease with age (Huppi et al., Reference Huppi, Murphy, Maier, Zientara, Inder, Barnes and Volpe2001).
The corpus callosum (CC), a vital white matter structure that is regarded as a marker of functional laterality, has been associated with NRH in the normal population (Witelson, Reference Witelson1985). Some studies have shown that non-right handed individuals have a larger CC than their right handed counterparts (Tuncer, Hatipoglu, & Ozates, Reference Tuncer, Hatipoglu and Ozates2005; Witelson, Reference Witelson1985, Reference Witelson1989); however, this has not been consistently reported (Welcome et al., Reference Welcome, Chiarello, Towler, Halderman, Otto and Leonard2009; Westerhausen et al., Reference Westerhausen, Kreuder, Dos Santos Sequeira, Walter, Woerner, Wittling and Wittling2004). Furthermore, left-handedness has been reported to be associated with higher FA and lower MD in the CC (Westerhausen et al., Reference Westerhausen, Kreuder, Dos Santos Sequeira, Walter, Woerner, Wittling and Wittling2004, Reference Westerhausen, Walter, Kreuder, Wittling, Schweiger and Wittling2003). The posterior limb of the internal capsule (PLIC) is another region thought to be associated with handedness (Westerhausen, Huster, Kreuder, Wittling, & Schweiger, Reference Westerhausen, Huster, Kreuder, Wittling and Schweiger2007). Corticospinal tract (CST) fibers that arise from the primary motor and premotor areas and descend through the PLIC are known to be functionally involved in fine and gross movements (e.g., finger, hand, and arm). While much is known about handedness-related asymmetries in the cortical motor areas, such as the premotor and sensorimotor regions (Herve, Crivello, Perchey, Mazoyer, & Tzourio-Mazoyer, Reference Herve, Crivello, Perchey, Mazoyer and Tzourio-Mazoyer2006; Westerhausen et al., Reference Westerhausen, Huster, Kreuder, Wittling and Schweiger2007), findings regarding handedness-related asymmetries of CST connections, such as the PLIC, are controversial, requiring further investigation (Seizeur et al., Reference Seizeur, Magro, Prima, Wiest-Daessle, Maumet and Morandi2014; Westerhausen et al., Reference Westerhausen, Huster, Kreuder, Wittling and Schweiger2007).
Preterm birth is thought to alter the microstructure and morphology of the CC and PLIC. Research suggests that CC size, shape, FA, and diffusivity are reduced or altered in preterm populations compared with controls, particularly in the posterior subregions (Andrews et al., Reference Andrews, Ben-Shachar, Yeatman, Flom, Luna and Feldman2010; Caldu et al., Reference Caldu, Narberhaus, Junque, Gimenez, Vendrell, Bargallo and Botet2006; Constable et al., Reference Constable, Ment, Vohr, Kesler, Fulbright, Lacadie and Reiss2008; Lawrence et al., Reference Lawrence, Allen, Walshe, Allin, Murray, Rifkin and Nosarti2010; Nosarti et al., Reference Nosarti, Rushe, Woodruff, Stewart, Rifkin and Murray2004; Skiöld et al., Reference Skiöld, Horsch, Hallberg, Engström, Nagy, Mosskin and Ådén2010; Thompson et al., Reference Thompson, Inder, Faggian, Johnston, Warfield, Anderson and Egan2011). Similar findings have been reported within the CST, with VPT infants with moderate to severe brain abnormalities (intraventricular hemorrhage and ventriculomegaly on MRI) exhibiting impaired CST microstructural development as they reach term-equivalent age compared with VPT infants without abnormalities on MRI (Adams et al., Reference Adams, Chau, Poskitt, Grunau, Synnes and Miller2010). However, studies examining the association between CC and PLIC morphology and microstructure and NRH in VPT children are scarce, although one study reported no relationship between handedness and CC size (Nosarti et al., Reference Nosarti, Rushe, Woodruff, Stewart, Rifkin and Murray2004). We are not aware of any studies that have investigated the relationship between PLIC microstructural organization and NRH in VPT children.
The current study aimed to: (1) examine the association between handedness and neurodevelopmental outcomes at age 7 years and whether this relationship varies between VPT and term-born children; (2) explore the relationship between the hemisphere of neonatal brain injury and NRH in VPT children; and (3) evaluate the associations between neonatal CC and PLIC microstructural organization with NRH in VPT children.
It was hypothesized that stronger NRH would be associated with poorer neurodevelopmental outcomes. It was also expected that VPT children with left-sided neonatal cerebral injury would more likely exhibit stronger NRH at 7 years of age than VPT children with right-sided neonatal cerebral injury. Lastly, it was hypothesized that greater CC area, higher fractional anisotropy and lower diffusivity, particularly in the posterior CC subregions (i.e., isthmus and splenium), would be associated with stronger NRH in VPT children. Stronger rightward PLIC asymmetry (e.g., higher fractional anisotropy, lower diffusivity) was also expected to be associated with stronger NRH in VPT children.
Material and Methods
Participants
Participants were part of the Victorian Infant Brain Studies (VIBeS) cohort, a prospective, longitudinal study examining the development of children born VPT [<30 weeks’ gestational age (GA) and/or very low birthweight (<1250 g)]. Recruitment occurred between July 2001 and December 2003 from the Royal Women’s Hospital in Melbourne, Victoria, Australia. Two hundred twenty-seven eligible VPT infants without genetic or congenital abnormalities known to affect development were recruited; however, two infants subsequently died, and 1 was later excluded due to a late diagnosis of congenital infection known to affect developmental outcome, leaving 224 VPT infants. A term control group (37 to 42 weeks’ GA and/or birthweight ≥2500 g) consisting of 77 infants were recruited either at birth from the Royal Women’s Hospital or the community (n=46), or at 2 years of age from Maternal and Child Health Centers within metropolitan Melbourne (n=31). Neonatal brain MRI was performed on all VPT and term infants who were recruited during the neonatal period. Within the VPT sample, 111 infants had diffusion weighted imaging. This smaller sample was due to the late inclusion of diffusion sequences in the neonatal scanning protocol. Those with diffusion data were similar to those without in their neonatal and demographic characteristics. Of these, 79 VPT infants had useable diffusion data and handedness data at age 7 years (Figure 1). Follow-up assessments were conducted at 2, 5, and 7 years of age corrected for prematurity for all children (Reidy et al., Reference Reidy, Morgan, Thompson, Inder, Doyle and Anderson2013; Roberts, Lim, Doyle, & Anderson, Reference Roberts, Lim, Doyle and Anderson2011; Thompson et al., Reference Thompson, Wood, Doyle, Warfield, Lodygensky, Anderson and Inder2008; Treyvaud et al., Reference Treyvaud, Ure, Doyle, Lee, Rogers, Kidokoro and Anderson2013). At 7 years of age, 197 VPT children (88% of survivors) and 70 term children (91%) were seen for follow-up. Handedness data were available for 175 VPT and 69 term children at age 7 years (Figure 1). Reasons for missing handedness information in both groups were due to incomplete handedness inventories (n=17) or intellectual and/or physical disability that made inventory administration unreliable (n=6) (Figure 1). The original and follow-up studies were approved by the Human Research Ethics Committee of the Royal Women’s Hospital and the Royal Children’s Hospital in Melbourne, Australia. Written consent was obtained from families before data collection.

Fig. 1 Flowchart of participants in the current study. N.B. Data collection at the 2- and 5-year follow-up not included.
Procedure
During the perinatal period, a neonatal assessment and a term-equivalent MRI scan were performed. The 7-year follow-up included a neuropsychological assessment, and parent-report questionnaires relating to the child’s behavior and daily functioning.
Measures
Perinatal measures
Perinatal data were obtained from chart review and information pertaining to GA, sex, birth weight, multiple birth, length of hospital stay, bronchopulmonary dysplasia, cystic periventricular leukomalacia, sepsis, necrotizing enterocolitis, patent ductus arteriosus, and grade 3 or 4 intraventricular hemorrhage was recorded. Bronchopulmonary dysplasia was defined as the requirement for oxygen at or beyond 36 weeks’ GA. Necrotizing enterocolitis and sepsis comprised proven episodes, and intraventricular hemorrhage was recorded as the highest grade from serial ultrasound scans throughout the neonatal intensive care course (Shah et al., Reference Shah, Doyle, Anderson, Bear, Daley, Hunt and Inder2008).
MRI
At term-equivalent age, infants were scanned in a 1.5 Tesla General Electric MRI scanner (Signa LX Echospeed System; General Electric, Fairfield, Connecticut) at the Royal Children’s Hospital, without sedation. Infants underwent T 1 (0.8–1.6 mm coronal slices; flip angle 45°; repetition time 35 ms; echo time 9 ms; field of view 21×15 cm2; matrix 256×192), T 2 (1.7–3 mm coronal slices; repetition time 4000 ms; echo time 60/160 ms; field of view 22×16 cm2; matrix 256×192, interpolated 512×512) and line scan diffusion-weighted imaging (4- to 6-mm axial slices; two images at b=5 s/mm2; six non-collinear gradient directions at b=700 s/mm2).
Hemispheric injury assessment
Presence of neonatal brain injury on MRI was assessed and reviewed by a neonatal neurologist using a previously described scoring system for all VPT infants (Kidokoro, Neil, & Inder, Reference Kidokoro, Neil and Inder2013). To determine the hemisphere of injury, a single blinded rater (J.C.) reviewed reported cases of cerebral white matter and cerebellar unilateral injury based on the MRI scoring assessment (specifically the presence of focal unilateral cystic lesions and/or focal punctuate signal abnormality). Focal unilateral injury in either the left or right hemispheres was classified as left or right hemispheric injury accordingly.
CC segmentation and diffusion analysis
For those VPT infants with diffusion weighted imaging, corpora callosa were manually delineated on the mid-sagittal slice of the anterior commissure-posterior commissure (AC–PC) aligned T 1-weighted image using three-dimensional (3D) slicer software (www.slicer.org). CC tracing was performed by a trained operator who was blinded to all perinatal data and neurodevelopmental results. The method used to determine the cross-sectional area of the CC has been previously described (Thompson et al., Reference Thompson, Inder, Faggian, Johnston, Warfield, Anderson and Egan2011, Reference Thompson, Inder, Faggian, Warfield, Anderson, Doyle and Egan2012). Diffusion measures within the CC were obtained by overlaying the CC mask onto the diffusion map that had been co-registered with the AC-PC aligned T 1 image. MD, FA, AD, and RD were quantified within the whole CC as well as within the 6 CC subregions (genu, rostral body, anterior mid-body, posterior mid-body, isthmus, and splenium) (Thompson et al., Reference Thompson, Inder, Faggian, Johnston, Warfield, Anderson and Egan2011).
PLIC diffusion analysis
MD, FA, AD, and RD were calculated within manually selected regions of interest on an axial slice of the brain taken at the level through the basal ganglia and PLIC. Bilateral regions of interest included the middle third of the PLIC, as previously described (Cheong et al., Reference Cheong, Thompson, Wang, Hunt, Anderson, Inder and Doyle2009). PLIC asymmetry was examined by calculating an asymmetry quotient,
${{(R{\minus}L)} \over {(R{\plus}L){\times} 0.5}},$
for the four diffusion measures (MD, FA, AD, and RD). A negative asymmetry quotient indicates leftward asymmetry, while a positive asymmetry quotient indicates rightward asymmetry.
Neurodevelopmental measures
The neuropsychological assessment conducted at 7 years’ corrected age included measures of intellect, academic achievement, working memory, language, and gross and fine motor skills. All assessment results are reported using standardized scores based on the child’s corrected age to avoid bias in cognitive test scores (Wilson-Ching, Pascoe, Doyle, & Anderson, Reference Wilson-Ching, Pascoe, Doyle and Anderson2014). All assessments were performed by trained pediatric assessors, psychologists or pediatricians.
Handedness was assessed through a 12-item, self-report hand preference inventory (Briggs & Nebes, Reference Briggs and Nebes1975) based on Annett’s (Reference Annett1967) questionnaire. Scores on this measure range from −24 to +24, with higher scores representing more right-handedness. An overall handedness score of 8 or less indicated NRH (scores of −9 to 8 are considered mixed-handed; scores of −24 to −9 are considered left handed), whereas a score of 9 and above signified right-handedness. Familial left-handedness was assessed via questionnaire where the primary caregiver was asked to report if either parent was left-handed (Briggs & Nebes, Reference Briggs and Nebes1975). General intellectual functioning was assessed using the Full Scale IQ (FSIQ) score of the Wechsler Abbreviated Scale of Intelligence (Wechsler, Reference Wechsler1999). Language ability was measured using the Core Language score from the Clinical Evaluation of Language Fundamentals – Fourth edition Australian (Semel, Wiig, & Secord, Reference Semel, Wiig and Secord2006). Reading, spelling, and mathematical computational skills were measured through the Wide Range Achievement Test – Fourth edition (Wilkinson & Robertson, Reference Wilkinson and Robertson2006). Verbal working memory was assessed using the Backward Digit Recall subtest from the Working Memory Test Battery for Children (Pickering & Gathercole, Reference Pickering and Gathercole2001). Each of these tests has a normative mean of 100 and standard deviation (SD) of 15. Motor skills were assessed using the Movement Assessment Battery for Children – version 2 (Henderson, Sugden, & Barnett, Reference Henderson, Sugden and Barnett2007), where standard scores for manual dexterity (fine motor), aiming and catching, and balance (gross motor), were combined to form a total movement score (global motor) (normative mean=10; SD=3). Parents rated their child’s behavior using the Strengths and Difficulties questionnaire – Australian English version (Goodman, Reference Goodman1997). The total difficulty score was used as an estimate of behavioral problems in daily life, with a higher score reflecting greater behavioral difficulty (range: 0–40).
Social risk
Social risk at age 7 years was determined based on a composite measure assessing family structure, education of the primary caregiver, occupation and employment status of the primary income earner, language spoken at home, and maternal age when the child was born (Roberts et al., Reference Roberts, Howard, Spittle, Brown, Anderson and Doyle2008), where higher scores indicated greater social risk.
Data Analysis
Data were analyzed using Stata 13.0 (StataCorp, 2013). To address our first aim, the main and interaction effects of handedness score and birth group on neurodevelopmental outcomes at 7 years were examined using linear regression. Analyses were adjusted for corrected age at assessment, social risk, and parental left-handedness. Wilcoxon-Mann-Whitney tests were used to compare the distribution of handedness between VPT children with left- versus right-hemisphere injuries in the cerebral white matter and cerebellum (Aim 2). Lastly, linear regression models were used to examine the association between CC area and diffusion measures and PLIC diffusion asymmetry and handedness at 7 years (Aim 3). Analyses were adjusted for age at scan and parental left-handedness. Analyses involving CC area were also adjusted for intracranial volume. Regression models were fitted using Generalized Estimating Equations with an exchangeable correlations structure and are reported with robust standard errors to allow for correlations between twins/triplets in the study (Carlin, Gurrin, Sterne, Morley, & Dwyer, Reference Carlin, Gurrin, Sterne, Morley and Dwyer2005). When regression models could not be fitted using Generalized Estimating Equations, regression models were reported using robust (sandwich) estimates of standard errors to allow for correlations between twins/triplets. Due to the issue of multiple comparisons, results were interpreted based on the profile and magnitude of differences and associations, rather than simply relying on p-values.
Results
Sample Characteristics
Regarding perinatal and demographic characteristics, participants with handedness data at age 7 years were representative of the original VPT and term-born cohorts, except that rates of cerebral palsy at 2 years of age were lower in VPT participants compared with VPT non-participants (4% vs. 24%) (Table 1). Of VPT children with handedness data, children who did not have neonatal diffusion data had lower birth weights (mean difference −84; 95% CI [−149, −18]; p=.01) and poorer language performance at 7 years (mean difference −7.0; 95% CI [−12.1, −1.9]; p=.007) compared with VPT children with diffusion imaging. As expected, the VPT group differed from controls on perinatal medical variables and the proportion of singletons. VPT children were also more socially disadvantaged. Rates of parental left-handedness were similar between VPT and controls. There was a trend for stronger NRH in the VPT group compared with controls (regression coefficient (b) −2.97; 95% CI [−6.25, 0.31]; p=.08). Following adjustment for parental left-handedness this association disappeared (b −1.95; 95% CI [−5.67, 1.77]; p=.30).
Table 1 Perinatal and demographic characteristics of study participants

a For some variables, the sample size is less than the total sample due to missing data; actual sample reported in table.
NRH=non-right handedness; GA=gestational age; M=mean; SD=standard deviation.
NRH and Neurodevelopmental Outcomes at 7 Years of Age
As reported previously, VPT birth compared with birth at term was associated with poorer performance across all neurodevelopmental outcomes, even after adjusting for social risk and age at assessment (p<.05 for all outcomes; data not presented) (Omizzolo et al., Reference Omizzolo, Scratch, Stargatt, Kidokoro, Thompson, Lee and Anderson2014; Reidy et al., Reference Reidy, Morgan, Thompson, Inder, Doyle and Anderson2013).
Table 2 Association between handedness and neurodevelopmental outcomes by birth group at 7 years of age

ba=regression coefficient; a positive value indicates that a lower score is associated with stronger NRH. If the regression is statistically significant, the 95% CI will not contain 0.
b Adjusted for social risk at age 7 years, corrected age at assessment and parental left-handedness.
* p value for the interaction between birth group and handedness.
For reading and spelling, there was evidence of interactions between handedness and birth group, with stronger NRH associated with poorer performance in the VPT group, but slightly better performance in the term controls (Table 2). Poorer mathematics and working memory performances were associated with stronger NRH in the VPT group. There was a trend for poorer language outcomes and stronger NRH in the VPT group. There was little evidence of interactions between handedness and birth group in other neurodevelopmental outcomes.
Hemisphere of Neonatal Brain Injury and NRH at 7 Years of Age in Children Born VPT
Only a small number of children had severe unilateral brain injury on the neonatal MRI, and as such analyses were exploratory (see Table 3). Wilcoxon-Mann-Whitney tests found that right handed and non-right handed groups did not differ based on the hemisphere of unilateral cerebral white matter injury in VPT children (p>.05). However, there was a trend for a handedness effect for cerebellar hemorrhage (z=−1.94; p=0.05), with a higher handedness rank (greater right-handedness) for left-sided cerebellar hemorrhage.
Table 3 Frequency of neonatal MRI unilateral cerebral white matter and cerebellar injury based on NRH dichotomy in VPT children

a One right-handed child was found to have both right-sided cystic lesions and signal abnormality.
Neonatal CC Area and Diffusion and NRH at 7 Years of Age in Children Born VPT
Increased FA and decreased RD in the whole CC, as well as increased FA within the genu, were associated with stronger NRH (Table 4). Within the splenium, increased area, increased FA and decreased MD and RD were associated with stronger NRH.
Table 4 Neonatal area and diffusion measures within the CC and its subregions in VPT infants and the association with handedness at 7 years of age

ba=regression coefficient; a negative value indicates that higher score is associated with stronger NRH.
b Adjusted for age at scan, intracranial volume, and parental left handedness.
c Adjusted for age at scan and parental left handedness.
FA=fractional anisotropy; MD=mean diffusivity; AD=axial diffusivity; RD=radial diffusivity.
Neonatal PLIC Diffusion Asymmetry and NRH at 7 Years of Age in Children Born VPT
There was a trend for greater rightward FA asymmetry in the PLIC being associated with greater NRH in VPT children (b −0.37; 95% CI [−0.80, 0.06]; p=.09). There was little evidence for associations between MD, AD and RD asymmetry quotients and NRH in VPT children (p>.05; data not presented).
Discussion
The current study aimed to explore neurodevelopmental outcomes and neural mechanisms associated with handedness in VPT children. While we found weak evidence that VPT birth was associated with stronger NRH, this was in line with findings from a meta-analysis where 17 studies also reported weak evidence of more NRH in preterm children compared with term-born children (Domellof et al., Reference Domellof, Johansson and Ronnqvist2011). Detectable differences in NRH between VPT and term children may be related to the age at which handedness was assessed, as well as birth weight and gestational age characteristics of preterm participants. Absolute (right vs. left) hand preference may emerge with age, leading to greater representation of left-handers in older compared with younger cohorts, and stronger NRH may be expected with decreasing gestational age and/or birthweight (Domellof et al., Reference Domellof, Johansson and Ronnqvist2011). Of interest, larger-scale population studies have reported no evidence for a link between birth stress (e.g., preterm birth and low birth weight factors) and hand preference (Johnston et al., Reference Johnston, Nicholls, Shah and Shields2009; Nicholls, Johnston, & Shields, Reference Nicholls, Johnston and Shields2012; Van der Elst et al., Reference Van der Elst, Hurks, Wassenberg, Meijs, Van Boxtel and Jolles2011).
NRH and Neurodevelopmental Outcomes at 7 Years of Age
NRH was associated with poorer neurodevelopmental outcomes in VPT children at age 7 years but not for term born controls. These findings are similar to Marlow and colleagues (2007) who noted an association between hand preference and cognitive development at age 6 years in children born <25 weeks’ gestation. Considering that findings regarding the association between NRH and neuropsychological outcomes have largely been inconsistent in preterm children, and several studies were unable to show links between NRH and cognitive function (Luciana, Lindeke, Georgieff, Mills, & Nelson, Reference Luciana, Lindeke, Georgieff, Mills and Nelson1999; Powls et al., Reference Powls, Botting, Cooke and Marlow1996; Saigal et al., Reference Saigal, Rosenbaum, Szatmari and Hoult1992), our findings are important in providing further evidence that NRH is associated with poorer neuropsychological outcome in this population. The relationship between NRH and neuropsychological outcomes in VPT children in the current study was most pertinent for academic and working memory outcomes. In contrast to previous studies, we did not find motor abilities, in particular fine motor skills, were associated with handedness. It may be that the tasks of fine motor dexterity used were not sensitive to hand skill differences compared with other motor tasks (e.g., cutting, hammering) that could not be considered in the current study. Previous research has also suggested that the performance of the non-dominant hand should also be considered: Ross, Lipper, and Auld (Reference Ross, Lipper and Auld1992) found that non-right handed VPT children were less skilled at using their non-dominant hand compared with their right-handed peers. Regarding term-born children, our results are in contrast to a meta-analysis of typically developing children (Somers et al., Reference Somers, Shields, Boks, Kahn and Sommer2015), which showed a small advantage in verbal and spatial ability for right-handed children. It is likely that our study lacked sufficient power or the appropriate sample composition to detect a handedness effect on cognitive outcomes for our comparison group.
Potential Neural Mechanisms of NRH in VPT Children
Potential underlying neural mechanisms for NRH in VPT children were examined in the current study. Handedness was not related to the hemisphere of unilateral cerebral white matter injury, however there was a trend for stronger right-handedness in VPT children with left-sided cerebellar hemorrhage. While our results do not strongly support a pathological model of NRH, only a small number of VPT children had unilateral cerebral (n=6) and cerebellar (n=7) injuries, so the study was likely underpowered to detect relationships between lateralized brain injury and NRH. To our knowledge, this is the first study to specifically investigate lateralized cerebral white matter and cerebellar brain injury (detected via MRI) in relation to handedness in VPT children specifically. While further investigation using advanced imaging techniques and larger samples are needed, our study highlights the possibility that stronger NRH in VPT children might be associated with early brain injury.
The rate of severe brain abnormalities detected on cranial ultrasound is decreasing in contemporary preterm and low birth weight cohorts (Wilson-Costello et al., Reference Wilson-Costello, Friedman, Minich, Siner, Taylor, Schluchter and Hack2007). In our study, 3–4% of VPT children had severe brain abnormalities (grade 3 or 4 intraventricular hemorrhage or cystic periventricular leukomalacia), compared to a rate of 6–7% in VPT infants born in the late 1990s (Heuchan, Evans, Henderson Smart, & Simpson, Reference Heuchan, Evans, Henderson Smart and Simpson2002). However, despite this decline in severe brain pathology, the rate of later neurodevelopmental difficulties continues to be high. This may suggest that more diffuse brain injury or subtle changes in brain development, such as alterations in neonatal brain morphology or microstructure, could explain some of the poorer functional outcomes that emerge in this population. Considering this, we investigated whether the expression of NRH in VPT children was a consequence of a disruption to normal brain asymmetry development, by examining the microstructural integrity of major white matter tracts implicated in motor development, such as the CC and PLIC. In contrast to previous research (Caldu et al., Reference Caldu, Narberhaus, Junque, Gimenez, Vendrell, Bargallo and Botet2006; Nosarti et al., Reference Nosarti, Rushe, Woodruff, Stewart, Rifkin and Murray2004), we found CC morphology and microstructure to be associated with NRH in our VPT cohort. Specifically, we found that larger neonatal splenium area was associated with stronger NRH in VPT children. We further noted that increased FA in the whole CC, specifically in the genu and splenium, and decreased MD and RD in the splenium, were related to stronger NRH, which is similar to findings by Westerhausen and colleagues (Reference Westerhausen, Kreuder, Dos Santos Sequeira, Walter, Woerner, Wittling and Wittling2004). Our findings support early research by Witelson (Reference Witelson1989) that highlighted that non-right handed individuals have larger CC than their right-handed counterparts. It has been demonstrated that VPT birth interferes with the maturation and myelination processes of the CC, particularly for posterior subregions (Caldu et al., Reference Caldu, Narberhaus, Junque, Gimenez, Vendrell, Bargallo and Botet2006; Constable et al., Reference Constable, Ment, Vohr, Kesler, Fulbright, Lacadie and Reiss2008; Lawrence et al., Reference Lawrence, Allen, Walshe, Allin, Murray, Rifkin and Nosarti2010; Nosarti et al., Reference Nosarti, Rushe, Woodruff, Stewart, Rifkin and Murray2004; Skiöld et al., Reference Skiöld, Horsch, Hallberg, Engström, Nagy, Mosskin and Ådén2010; Thompson et al., Reference Thompson, Inder, Faggian, Johnston, Warfield, Anderson and Egan2011). During the late prenatal and early postnatal period, there is a naturally occurring process of pruning of callosal axons. This process is thought to result in pronounced functional lateralization and reduce the likelihood of left-handedness (Westerhausen et al., Reference Westerhausen, Kreuder, Dos Santos Sequeira, Walter, Woerner, Wittling and Wittling2004). VPT birth may disrupt this process, leading to a greater number, or density, of callosal axons being present during the neonatal period, increasing the prevalence of NRH in the VPT population. Therefore, our findings suggest that the primary effect of preterm birth on CC microstructure and morphology subsequently influences the development and expression of hand preference at age 7.
The current study is the first to explicitly examine neonatal PLIC microstructure in relation to handedness in VPT children. While there was a trend for NRH in VPT children to be associated with stronger rightward PLIC FA asymmetry, our overall PLIC diffusion results provided little evidence to support an association between PLIC microstructure and handedness in this population. Our findings are in line with Westerhausen and colleagues (2007), who also reported little association between CST asymmetries and handedness. While others have noted an association between handedness and CST morphology, these studies either included small samples (Seizeur et al., Reference Seizeur, Magro, Prima, Wiest-Daessle, Maumet and Morandi2014), or were based on other MRI techniques (i.e., volume based on apparent grey matter density) (Herve et al., Reference Herve, Leonard, Perron, Pike, Pitiot, Richer and Paus2009). Findings from our study therefore suggest that neonatal CST asymmetry characteristics are not early determinants of hand preference in VPT children. Handedness may instead reflect asymmetries of cortical motor areas rather than asymmetries of motor tracts or connections, such as the CST (Herve et al., Reference Herve, Crivello, Perchey, Mazoyer and Tzourio-Mazoyer2006).
While the current study did not have the power to categorize and examine weak (i.e., mixed handedness) versus strong handedness laterality (i.e., left and right handedness), some research has suggested the “strength” may be more important than the “direction” of handedness (Welcome et al., Reference Welcome, Chiarello, Towler, Halderman, Otto and Leonard2009). It has been suggested that mixed handed individuals exhibit greater inter-hemispheric connectivity and subsequently have larger CC compared with left- or right-handed individuals (Welcome et al., Reference Welcome, Chiarello, Towler, Halderman, Otto and Leonard2009; Witelson & Goldsmith, Reference Witelson and Goldsmith1991). Behavioral and cognitive discrepancies have been noted between mixed handed individuals and consistent-handed individuals (Welcome et al., Reference Welcome, Chiarello, Towler, Halderman, Otto and Leonard2009). It has also been suggested that mixed handedness may be associated with poorer behavior and intellectual ability in children, compared with left-handedness (Domellof et al., Reference Domellof, Johansson and Ronnqvist2011; Rodriguez et al., Reference Rodriguez, Kaakinen, Moilanen, Taanila, McGough, Loo and Järvelin2010; Rodriguez & Waldenström, Reference Rodriguez and Waldenström2008). Considering this, it is imperative that future studies examine and highlight the influence of atypical handedness on outcomes in VPT children.
The current study has some limitations. The study used a questionnaire measure of hand preference and did not include measures of lateral performance when evaluating handedness. It is important to distinguish hand preference from hand performance as they have been acknowledged to be different and separate outcomes (Porac & Coren, Reference Porac and Coren1981; Scharoun & Bryden, Reference Scharoun and Bryden2014). Furthermore, a relatively small sample size due to loss to follow-up and available MRI data may have resulted in a lack of power to identify differences between lateralized unilateral brain injury and hand preference, and associations between PLIC microstructure and NRH.
Findings from the current study have contributed to our understanding of handedness in VPT children and its relationship with neurodevelopmental outcomes. This study was unique in its specific examination of potential neural mechanisms underpinning the development of NRH in VPT children using advanced imaging techniques and a standardized measure of hand preference. Considering findings from the current study, it is important to examine handedness over time via longitudinal studies to determine whether expressions and correlates of NRH in VPT children persist or change dynamically with neurodevelopment. The data reported here suggest that handedness may be an important marker for some cognitive and academic difficulties in VPT children, and indicate that white matter connections may be differentially involved in hand preference. These findings have implications for VPT children particularly in the educational setting, where built-in support structures may be beneficial to assist non-right handed children strengthen and develop particular cognitive or developmental skills.
Acknowledgments
We acknowledge the input of the VIBeS research team, the Melbourne Children’s MRI Center, and all the families who participated in this study. This study was funded by Australia’s National Health & Medical Research Council [Project Grants (237117 & 491209), Early Career Award (1012236 to D.T.), research fellowship support (1081288 to P.A.; 1053787 to J.C.), Centre for Research Excellence in Newborn Medicine (1060733)]; and National Institutes of Health (HD058056). We also acknowledge the support provided by the Murdoch Childrens Research Institute, the Royal Children’s Hospital Foundation, and the Victorian Government’s Operational Infrastructure Support Program. This manuscript has never been published elsewhere and there are no conflicts of interest.