Children With Dyslexia and Typical Readers: Sex-Based Choline Differences Revealed Using Proton Magnetic Resonance Spectroscopy Acquired Within Anterior Cingulate Cortex

Children with dyslexia exhibit slow and inaccurate reading, as well as problems in executive functions. Decreased signal activation in brain regions related to visual processing and executive functions has been observed with functional magnetic resonance imaging with reports of sex differences in brain patterns for visual processing regions. However, the underlying neurochemistry associated with deficits in executive functions for children with dyslexia has not been thoroughly evaluated. Reading ability and executive functions were assessed in fifty-three children [ages 8–12 years old, dyslexia (n = 24), and typical readers (n = 30)]. We employed short echo, single voxel, proton magnetic resonance spectroscopy to evaluate the perigenual anterior cingulate cortex (ACC). Pearson correlations were calculated between metabolite concentrations and measures of reading, processing speed, and executive function. Logistic regression models were used to determine the effects of brain metabolite concentrations, processing speed, and reading scores on dyslexia status. Differences by child’s sex were also examined. Compared to typical readers, higher global executive composite t-score is associated with greater odds for dyslexia (OR 1.14; 95% CI 1.05, 1.23); increased processing speed appears to be protective for dyslexia (OR 0.95; 95% 0.89–1.00). After adjustment for multiple comparisons, females with dyslexia showed strong and significant negative correlations between processing speed and myo-inositol (r = -0.55, p = 0.005) and choline (r = -0.54, p = 0.005) concentrations; effect modification by sex was confirmed in linear regression models (psex∗Cho = 0.0006) and (psex∗mI = 0.01). These associations were not observed for males or the group as a whole. These findings suggest that children with dyslexia share difficulty in one or more areas of executive function, specifically those related to response time. Also, metabolite changes in the ACC may be present in children with dyslexia, especially for females, and may hold value as possible markers for dyslexia.

Children with dyslexia exhibit slow and inaccurate reading, as well as problems in executive functions. Decreased signal activation in brain regions related to visual processing and executive functions has been observed with functional magnetic resonance imaging with reports of sex differences in brain patterns for visual processing regions. However, the underlying neurochemistry associated with deficits in executive functions for children with dyslexia has not been thoroughly evaluated. Reading ability and executive functions were assessed in fifty-three children [ages 8-12 years old, dyslexia (n = 24), and typical readers (n = 30)]. We employed short echo, single voxel, proton magnetic resonance spectroscopy to evaluate the perigenual anterior cingulate cortex (ACC). Pearson correlations were calculated between metabolite concentrations and measures of reading, processing speed, and executive function. Logistic regression models were used to determine the effects of brain metabolite concentrations, processing speed, and reading scores on dyslexia status. Differences by child's sex were also examined. Compared to typical readers, higher global executive composite t-score is associated with greater odds for dyslexia (OR 1.14; 95% CI 1.05, 1.23); increased processing speed appears to be protective for dyslexia (OR 0.95; 95% 0.89-1.00). After adjustment for multiple comparisons, females with dyslexia showed strong and significant negative correlations between processing speed and myo-inositol (r = −0.55, p = 0.005) and choline (r = −0.54, p = 0.005) concentrations; effect modification by sex was confirmed in linear regression models (p sex * Cho = 0.0006) and (p sex * mI = 0.01). These associations were not observed for males or the group as a

INTRODUCTION Reading Difficulties and Executive Functions
Reading difficulties (or dyslexia) are characterized by slow and inaccurate reading which continues into adulthood despite remedial intervention and exposure to the written language (International Dyslexia Association IDA, 2011). Specific challenges related to phonological and orthographical processing deficits, and more broadly the reading process, have been observed in individuals with dyslexia (Seidenberg and McClelland, 1989;Pugh et al., 2000). However, our accumulated data suggests that children with reading difficulties also demonstrate challenges in several higher order abilities, i.e., in executive functions (Horowitz-Kraus, 2014). More specifically, we found that children and adults with dyslexia demonstrate decreased error monitoring ability, as manifested using electroencephalographic (EEG)-event related potentials (ERP) amplitudes following a commission of an error (i.e., decreased error related negativity potential) compared to age matched typical readers Breznitz, 2008, 2013;Horowitz-Kraus, 2011, 2012. This alteration was not specific to reading materials but was also extended to tasks which do not contain verbal stimulation (i.e., the Wisconsin task) (Horowitz-Kraus, 2014). This ERP pattern, which is related to a mismatch between an actual and a desired response, is an evoked post-response and is related to a change in dopaminergic serge (Falkenstein et al., 1991(Falkenstein et al., , 2000. The error related negativity evoked from the anterior cingulate cortex (ACC) (Falkenstein et al., 2000) is part of the error detection system in the brain. The ACC is a critical part of the cingulo-opercular network which is related to a top-down monitoring process (Dosenbach et al., 2008). Support for the altered brain activation related to the ACC was found in our study showing a decreased functional connection within the cingulo-opercular network during reading in children with dyslexia and an increased functional connectivity following intervention that accompanied reading achievement (Horowitz-Kraus et al., 2015b). Since decreased ERPs are related to decreased neuronal firing, which in turn points at decreased signal activation on functional magnetic resonance imaging (fMRI), it is therefore not surprising that this network showed decreased functional connections with visual processing during a reading task in children with dyslexia, and altered functional connection during a Stroop task in this population (Levinson et al., 2018) emphasizing its critical role in monitoring during reading (Horowitz-Kraus et al., 2015a) as well as in the absence of a task (i.e., during rest) (Horowitz-Kraus et al., 2015b). Other ERP and fMRI studies revealed differences between individuals with dyslexia and typical readers, including manipulations aimed at improving reading through a specific triggering of executive functions, increased ERPs related to error monitoring (Horowitz-Kraus and Breznitz, 2014), and increased activation of the ACC  and of functional connectivity of the cingulo-opercular network (Horowitz-Kraus et al., 2015b). These EEG and fMRI findings led us to examine the specific differences in metabolite concentrations within the ACC mediating reading ability.
Are There Neurochemical Characteristics for Reading Difficulties? Lebel et al. (2016) observed using proton magnetic resonance spectroscopy (MRS) in typical preschool children that phonological processing ability was positively associated with glutamate (Glu), creatine (Cr), and myo-inositol (mI) concentrations in the pregenual anterior cingulate. As mI concentrations are thought to have a role as a marker of glial cells, the increased levels supported increased Glu neurotransmission and this metabolism provided a coupled relationship with phonological processing. However, most studies in the field of dyslexia focus on the role of phonological and orthographical routes (Shaywitz and Shaywitz, 2003a,b;Shaywitz et al., 2003). Pugh et al. (2014) characterized metabolite levels in children with dyslexia within the occipital cortex. This study employed English speaking, 6-10 year old typical readers and individuals with dyslexia and suggested a negative correlation between Cho and Glu concentrations with reading and phonological processing abilities (Pugh et al., 2014). The authors speculated that elevated Glu, a marker for hyperexcitability, may influence the coherence of neuronal networks involved in learning, which may also be the case for other pathologies. These researchers suggested that an unstable performance of children with reading difficulties is characterized with a moment to moment variance in performance and a lack of consistency (Pugh et al., 2014). This may provide a link to the monitoring challenges found in this population, which are not restricted to the reading domain. In a subsequent study, Del Tufo et al. (2018) found in children that cross-modal word matching mediates the relationship with increased Glu and increased Cho with poorer reading ability. Given that Cho levels represent membrane turnover as well as cellular and white matter density (Miller et al., 1996), elevated levels are in-turn related to increased connectivity or atypical myelination of the occipital cortex with additional regions in children with dyslexia (Pugh et al., 2014). Similarly, higher Cho levels in the angular gyrus are also associated with low reading scores in adults (Bruno et al., 2013). There it was suggested that Cho may be related to less efficient neuronal activity in this regions related to phonological processing and lower reading scores (Bruno et al., 2013). Bruno et al. (2013) also indicated that N-acetylaspartate (NAA) levels in the angular gyrus, a marker for intact neuronal ability, was related to a higher cognitive ability in adults (Jung et al., 2005). Del Tufo et al. (2018) also found that higher NAA predicted faster cross-modal matching reaction times in children. Increased NAA in the prefrontal cortex was also related to enhanced metabolic rates in these regions, probably due to differences in prefrontal maturation, as was observed in Asperger patients (Murphy et al., 2002). Higher NAA was previously related as a marker of viable neurons (Erickson et al., 2012) and there is conflicting evidence of positive and negative correlations (Mohamed et al., 2010) of NAA levels and cognitive control abilities. However, there is still a gap in knowledge as to the significance of metabolite concentrations in regions specifically related to executive functions in children with dyslexia.
Therefore, the goals of the current study were to: (1) determine whether metabolite concentrations in the ACC predict silent reading and sight word efficiency, processing speed or executive functions, (2) evaluate whether metabolite concentrations, reading ability, processing speed or executive functions vary in children with dyslexia compared to typical readers, and (3) determine if the associations vary by sex.

Participants
Participants were recruited from posted ads and through commercial advertisements. All participants gave their informed written assent and their parents provided informed written consents prior to inclusion in the study. This study was approved by the Cincinnati Children's Hospital Medical Center Institutional Review Board.

Baseline Reading Measures Used to Confirm Dyslexia Status
We confirmed the existence of reading difficulties using a battery of normative reading tests in English. Children with dyslexia were all diagnosed as having difficulty reading. Inclusion criteria for the dyslexia group were standard score of −1 and below or meeting the 25% or below cutoff in words reading, decoding and comprehension abilities. This battery included (a) words reading accuracy/orthography: the "Letter-Word" subset from the Woodcock and Johnson-III (WJ-III) battery (Woodcock and Johnson, 1989); (b) decoding: the "Word Attack" subset from the WJ-III (Woodcock and Johnson, 1989); and (c) reading comprehension subset from the WJ-III (Horowitz-Kraus et al., 2016). Participants in the typical readers group were age-matched students who volunteered for the study with fluent and accurate reading according to established normative levels for the WJ-III.

Reading Abilities
We further evaluated children's reading abilities using the Tests of Word Reading Efficiency (TOWRE) to assess the participant's ability to pronounce printed words (sight word efficiency) and the Tests of Silent Reading Efficiency and Comprehension (TORSEC) to assess the participant's fluency ability.

Executive Functions and Cognitive Abilities
To assess executive functions and cognitive abilities, we used the Behavioral Rating Inventory of Executive Functions (BRIEF) (Gioia et al., 2002) parents questionnaire as well as the speed of processing tests from the Wechsler Intelligence Scale for Children, Fifth Edition (WISC-V) (Wechsler, 2014). The BRIEF questionnaire includes questions covering the child's cognitive abilities (i.e., inhibition, organization, attention, monitoring, and emotional regulation). In the speed of processing subtests of the WISC-V, participants performed the symbol search and the coding sub-tests, both result in the Processing Speed Index (PSI) employed for the analyses.

Acquisition Methods
Brain magnetic resonance imaging (MRI) and spectroscopy (MRS) were acquired using a Philips Achieva MR scanner operating at 3 Tesla (3T) and equipped with a 32-channel head coil. A three-dimensional (3D), high-resolution, isotropic, T1weighted fast Fourier echo (FFE) anatomical imaging sequence was performed using 8.2 ms repetition time (TR), 3.7 ms echo time (TE), 1057 ms inversion time (TI), 8 degree flip angle, sensitivity encoding factor (SENSE) of 2, contiguous slices with a 1 mm thickness, and 1 × 1 mm voxel size. A single voxel, point resolved spectroscopy (PRESS) sequence was conducted using a 2000 ms TR, 30 ms TE, and 96 averages with water suppression along with an embedded unsuppressed water reference series of 16 averages. The 8 cubic centimeter single voxel was prescribed about the perigenual ACC within the medial frontal lobe localized from the 3D T1-FFE anatomical imaging sequence similar to the position described by Lebel et al. (2016) See Figure 1.
Participants were acclimated and desensitized to condition them for comfort inside the scanner (see Byars et al., 2002 for details). Head motions were controlled using elastic straps that were attached to either side of the head-coil apparatus. An MRI-compatible audio/visual system (Avotec, SS3150/ SS7100) was used for the presentation of a movie during the session.

MRS Data Analysis
The raw spectroscopy data were imported into LCModel (Provencher, 1993) commercial software for quantitative FIGURE 1 | (A) Representative location within the anterior cingulate cortex for the 8 cubic centimeter (2 cm per side) spectroscopic voxel positioned on T1 weighted imaging slices centered in the sagittal, coronal and axial plane orientations. (B) Representative single voxel, short echo, magnetic resonance spectrum with N-acetylasparate appearing on the x-axis with a peak at a chemical shift value of 2 parts per million (PPM), glutamate at 2.3-2.4 PPM, creatine at 3.0 PPM, choline at 3.2 PPM and myo-inositol at 3.5 PPM.
processing. Metabolite and water reference levels were determined in institutional units. The raw metabolite levels were adjusted for the tissue contributions from gray matter, white matter and cerebrospinal fluid (CSF) using FSL (Woolrich et al., 2009) adjusted to the T1 and T2 relaxation decay rate of the corrected water concentration and corrected for literature reported T1 and T2 relaxation decay rates of the primary metabolites including NAA, Cho, Cr, and mI (Wansapura et al., 1999;Traber et al., 2004;Tal et al., 2012). However, pure Glu as well as the combined glutamate and glutamine (GLX) concentrations were unadjusted for metabolite T1 and T2 relaxation decay (Gussew et al., 2012). Reported values were concentrations in units of millimolar (mM).

Statistical Analyses
First, two sample t-tests were used to examine differences in age and IQ measures, reading and executive function measures as well as metabolite concentrations among children with dyslexia and typical readers.
Second, Pearson (r) correlations were calculated between individual metabolite concentrations and reading scales, processing speed, and global executive function. Correlation differences by dyslexia status and sex were examined. Correction for multiple comparisons was conducted by False Discovery Rate (FDR) estimation.
Third, significant correlations (FDR cutoff of 0.05) were then evaluated by linear regression adjusting for child's sex and age at visit/test. Effect modification by sex was investigated by (1) including a two-way interaction term between metabolite concentration and child's sex in the model, and (2) performing sex-stratified linear regression if the two-way interaction term was significant (p < 0.05). Effect sizes (β) represent the change in outcome for an interquartile range (IQR) increase in metabolite concentration.
Fourth, to determine if metabolite concentrations, processing speed, and global executive function predict dyslexia status, logistic regression models were implemented. For logistic models of metabolites, all metabolites were included in the same model after the investigation of multi-collinearity (based on correlations among metabolites).
Lastly, to confirm reading difficulty among children with dyslexia, we evaluated the association between silent reading and sight word efficiency and dyslexia status using the TOWRE and TOSREC. While it is likely these measures would be associated with dyslexia, it should be noted that these specific measures were not used to define the individuals in our study (see section 2.2.1). Logistic models were adjusted for child's sex and age at visit/testing. Effect modification was evaluated by including a two-way interaction term between reading/behavior measure and child's sex in the model and by stratification if interaction term was significant. Analyses were conducted with SAS 9.4 (SAS Institute Inc., Cary, NC, United States) and SPSS13 (IBM, Armonk, NY, United States).
All participants were native English speakers, Caucasian, and with United States average socioeconomic status by income and education, as reported by the families using a socio-economic status questionnaire (Barratt, 2006) Participants were right-handed, displayed normal or correctedto-normal vision in both eyes, and had normal hearing. None had a history of neurological or mood disorders, and individuals with reported history of attention difficulties were excluded.

Behavioral Measures
Overall, children with dyslexia demonstrated significant lower scores in all reading measures taken in the study: these readers showed lower phonological and orthographical abilities. Executive functions and speed of processing in children with dyslexia were significantly lower than typical readers. See Table 1.

Neurochemical-Behavior Correlations
Pearson correlation coefficient (r) values between mI, NAA, Cr, Cho, Glu, GLX, reading scales, processing speed, and global executive function were evaluated. There were no significant correlations observed among the group as a whole; however, there were differences observed by dyslexia status and sex ( Table 2). Among children with dyslexia, there was a strong negative correlation between Cho and silent reading score (r = 0.51, p = 0.01); typical readers showed a strong negative correlation between mI and processing speed (r = −0.42, p = 0.02). After correction for multiple comparisons, no significant differences were found between group metabolite concentrations in children with dyslexia and typical readers.

Effect of Neurochemical Levels on Processing Speed
The overall effect of Cho and mI on processing speed was not significant; however there was indication of effect modification by child's sex. Figure 2 shows the effect of (a) Cho (p interaction = 0.0006) and (b) mI (p interaction = 0.01) on predicted mean scores for processing speed by child's sex. Higher levels of Cho (β −8.10; 95% CI −12.73, −3.45) and mI (β −5.22; 95%CI −8.22, −2.22) are associated with decreased processing speed among females but not males ( Table 3).

Effect of Global Executive Function, Processing Speed, Reading Measures, and Metabolites on Dyslexia Status
Children with a higher global executive composite t-score are more likely to have dyslexia (OR 1.14; 95% CI 1.05, 1.23) suggesting some difficulty in one or more areas of executive function is present among children with dyslexia compared to typical readers ( Table 4). We also observed a marginally significant (p = 0.049) association between increased processing speed and protective odds for dyslexia (OR 0.95; 95% 0.89-1.00) compared to typical readers. Metabolite concentrations did not predict dyslexia status. Metabolite concentrations were moderately to significantly correlated (r ranging from 0.22 to 0.59, r = 0.80 for GLX/Glu correlations) and given the strong correlation between GLX and Glu, a sensitivity analysis was conducted to determine if removal of one of these metabolites impacted the results; the findings remained insignificant. As expected, higher scores on silent reading and sight word efficiency are observed among typical readers as evident by the ORs less than 1 ( Table 4). There was no evidence of effect modification by child's sex (all interaction p-values > 0.20).  mI (myo-inositol), NAA (N-acetyl aspartate), Cr (creatine), Cho (Choline), Glu (glutamate), GLX (glutamate and glutamine), Silent reading (TOSREC), sight word efficiency (TOWRE-SWE), processing speed index (WISC), global executive composite (BRIEF). Bold text denotes significant correlation at p < 0.05. Correlations significant after correction for multiple comparisons are Cho/PSI (FDR = 0.05) and mI/PSI (FDR = 0.05). Cho and silent reading were correlated among children with dyslexia (r = −0.51, p = 0.01) and mI and processing speed were negatively correlated among typical readers (r = −0.42, p = 0.02); these were not significant after multiple comparisons.
FIGURE 2 | Effect of (A) choline and (B) myo-inositol on processing speed by sex. Graph represents the effect of (A) choline (Cho) and (B) myo-inositol (mI) on the predicted processing speed of participants by child's sex Sex-stratified regression models are adjusted for age at visit/testing. Blue and pink indicates the effect for males and females, respectively. Dashed lines represent 95% confidence intervals. Actual individual level data are overlaid for males (blue circles) and females (pink circles). Interaction p-value represents the p-value for the formal two-way interaction between metabolite level and child's sex in the linear regression model.

DISCUSSION
Over the past decade, multiple studies have implicated difficulties in executive functions for individuals with dyslexia. However, there is still a gap in knowledge as to the neurochemical characteristics related to reading difficulties in neural circuits associated with challenges in executive functions for this population. The current study aimed to characterize the metabolite concentrations related to reading ability in typical readers and in children with dyslexia, specifically focusing on the ACC. As hypothesized, lower reading ability was associated with high Cho levels, however, it was negatively associated with processing speed in females with dyslexia. The current study also demonstrated a negative relationship between processing speed and mI concentration. To our knowledge, this is among the first pediatric studies to look at MRS measures in the ACC and to report higher levels of Cho in the brain of females with dyslexia. Also, the study confirmed children, both boys and girls, with dyslexia demonstrated lower reading ability and executive functions scores compared to typical readers. The study did not observe relationships of the cognitive and reading measures with other metabolites, such as NAA, Glu, or GLX, nor did it observe metabolite changes in males. In this study, metabolite concentrations were determined with corrections for metabolite and water relaxation along with contributions from different tissue types contributing to the voxel. Other studies report relative metabolite levels to one another, usually to creatine, and fail to account for differences in tissue contributions, which may explain observed differences. The voxel placement within the pregenual ACC in the current study matches the location described by Lebel et al. (2016) however, it differs from that of other studies of individuals with dyslexia focused on the occipital lobe including lingual gyrus, calcarine sulcus, and cuneus.

High Cho and Reading Difficulty
Previous studies have pointed at high Cho levels related to reading levels in the occipital lobe (Pugh et al., 2014;Del Tufo et al., 2018) in the cerebellum (Laycock et al., 2008) and the left temporoparietal region (Bruno et al., 2013). The authors suggested that these high Cho levels in individuals with reading difficulties for reading-related regions reflect high membrane turnover, cellular density, and white matter density. These results are in line with evidence of impaired myelination in this population in white matter tracts passing the temporoparietal regions and the occipital lobe (Yeatman et al., 2012;Wandell and Yeatman, 2013). Our study is believed to be among the first to reveal metabolite changes in regions that were not traditionally included as part of the classical reading circuitry. However, in recent years there are additional models pointing at the critical role of executive functions-related neural circuits in reading, including the ACC (Horowitz- Horowitz-Kraus and Hutton, 2015). Our results provide additional evidence for the involvement of neural circuits related to executive functions in children with dyslexia extending the high Cho findings also to the ACC. The relationship to alterations in myelination (especially in the genu which passes though the ACC) should also be examined. As the ACC in a major hub in the cingulo-opercular network, (Dosenbach et al., 2008) it would be interesting to measure Cho levels in other sub-regions of the ACC (as outlined in) (de la Vega et al., 2016) and other brain regions that are related to the fronto-parietal network and the fast monitoring of cognitive processes also impaired in children with dyslexia.

Differences in mI Between Children With Dyslexia and Typical Readers
myo-inositol is a carbocyclic sugar molecule localized to glial cells. There are several known roles for mI in the brain. It is primarily considered a glial cell marker as it's concentration increases with glial cell based neoplasms and gliosis. However, it also functions as an osmolyte such that during periods of osmotic stress, as balance is maintained via mI transport across plasma membranes. However, in the context of the current study with Cho findings, there is support for the role of mI as a key precursor of phosphoinositides and phospholipids, cell membranes and myelin structures. The composition of these structures within the brain networks could directly influence reading abilities.

Sex-Specific Findings
We were able to explore sex-based differences as eight of twentyfour children with dyslexia and seventeen of thirty typical readers in our study were females. Sex differences are known in many disorders, including ADHD, and dyslexia, with greater male frequency with the diagnosis Shaywitz, 2003a,b, 2008). Multiple genetic, and non-inherited factors have been posited to explain this observation. The sex hormones, estrogen and progesterone, demonstrate cognitive and neuroprotective effects (Brann et al., 2007;Dumitriu et al., 2010). Arnett et al. (2017) concluded that the higher prevalence in males with reading difficulties can be explained by slower and more variable processing speed along with worse inhibitory control. Geschwind and Levitsky (1968) first reported structural differences and asymmetry in persons with dyslexia upon post-mortem brain examinations. Non-invasive neuroimaging investigations of brain structure have provided further insight into the neurobiological basis for sex-specific differences. Evans et al., found gray matter volume differences in adults and children, male and female with developmental dyslexia, specifically in the occipital lobe for girls (Evans et al., 2014).

Study's Limitations
The results of this study should be considered in light of the following limitations. First, the sample size of this study is relatively small. Second, the perigenual placement of the MRS voxel within the ACC limited the ability to relate to executive functions of the ACC assigned to more posterior aspects, such as cognitive control and conflict monitoring. Third, technical differences in the MRS acquisition and quantitation of metabolite levels, as previously discussed, can influence the results. Also, the relationship between BOLD and functional connectivity and our findings should be validated in a future study using a joint MRS-fMRI model combining the results of these two methodologies.

CONCLUSION
The current study's findings pinpoint metabolic differences related to the medial frontal lobe in females with dyslexia and typical readers. Further investigations are necessary to explore the metabolism and function of the ACC and how they influence reading abilities.