- Open Access
Abnormal auditory ERP N100 in children with dyslexia: comparison with their control siblings
Behavioral and Brain Functionsvolume 5, Article number: 26 (2009)
Recent research has implicated deficits of the working memory (WM) and attention in dyslexia. The N100 component of event-related potentials (ERP) is thought to reflect attention and working memory operation. However, previous studies showed controversial results concerning the N100 in dyslexia. Variability in this issue may be the result of inappropriate match up of the control sample, which is usually based exclusively on age and gender.
In order to address this question the present study aimed at investigating the auditory N100 component elicited during a WM test in 38 dyslexic children in comparison to those of 19 unaffected sibling controls. Both groups met the criteria of the International Classification of Diseases (ICD-10). ERP were evoked by two stimuli, a low (500 Hz) and a high (3000 Hz) frequency tone indicating forward and reverse digit span respectively.
As compared to their sibling controls, dyslexic children exhibited significantly reduced N100 amplitudes induced by both reverse and forward digit span at Fp1, F3, Fp2, Fz, C4, Cz and F4 and at Fp1, F3, C5, C3, Fz, F4, C6, P4 and Fp2 leads respectively. Memory performance of the dyslexics group was not significantly lower than that of the controls. However, enhanced memory performance in the control group is associated with increased N100 amplitude induced by high frequency stimuli at the C5, C3, C6 and P4 leads and increased N100 amplitude induced by low frequency stimuli at the P4 lead.
The present findings are in support of the notion of weakened capture of auditory attention in dyslexia, allowing for a possible impairment in the dynamics that link attention with short memory, suggested by the anchoring-deficit hypothesis.
Dyslexia describes individuals having reading related difficulties despite normal intelligence, adequate training and satisfactory educational opportunities. Approximately 5–17% of the population is considered to have dyslexia [1, 2]. Recent research has implicated deficits of the working memory (WM) operation and attention within the patterns of cognitive characteristics of dyslexia [3–5]. However, the underlying mechanisms and the cognitive characteristics of dyslexia still remain unclear [6–8].
Event-related potentials (ERP) provide an objective index of cognition-related brain activity because they reflect aspects of information processing that are relatively inaccessible using self-report or traditional behavioral techniques. The auditory N100 component of ERPs has been conceptualized as the physiological correlate of both attentional and working memory operation [9–13].
Previous studies involving dyslexic subjects in comparison to controls have investigated hypotheses relating information processing deficiency with patterns of N100 component using ERP methodology. The results of these studies have been variable. Several studies have found that subjects with dyslexia have reduced N100 amplitude [14–16] and prolonged latency  as compared to controls. It has been argued that this reduced and slowed activity within these early stages of information processing contributed to the language deficits of dyslexic children, [18, 19]. On the other hand there are studies that failed to find a discrepancy between dyslexic children and controls  or even demonstrated opposite effects .
The reason for inconsistencies in the reported results may be the diversity in the study design and methodology. Variability is likely to arise at least partly from the fact that the normative data used possibly do not constitute a fair comparison for dyslexic subjects, because they are drawn from a control population selected to represent a not well-defined samples data except the matching of demographic data of age and sex. It is apparent, for example, that children of low socioeconomic status usually suffer more than do children of high socioeconomic status from problems that impair performance on cognitive tests. Therefore, a promising comparison might be an investigation between dyslexic subjects and their sibling controls, since both groups share similar family relations and socioeconomic status, thus eliminating most of the above confounding factors. At this point, it is worth noting that our team studying dyslexic children in comparison to their sibling controls found that dyslexic children manifested abnormal aspects of pre-attentive processing of information as they are reflected by P50 elicited during a working memory test .
In view of the above considerations, we hypothesized that the electrophysiological brain activity, as reflected by N100, in association with working memory operations, could be of significance in identifying possible pathophysiological mechanisms involved in dyslexia. Thus, the present study was designed to determine: (a) whether dyslexic children as compared to their sibling controls have similar or different features of N100 ERP component elicited during a working memory test; (b) whether dyslexic children have more difficulties than their normal siblings with regard to memory performance, and (c) the strength of association, if any, between memory performance and the amplitudes and latencies of the N100, and whether this association would be similar or different for dyslexic children and their sibling controls.
Fifty seven (57) children participated in the study. Thirty eight (26 boys and 12 girls) of them were outpatient cases who had been diagnosed as having dyslexia according to the 10th edition of the International Classification of Diseases (ICD-10) and the rest 19 children (7 boys and 12 girls) were sibling controls of the dyslexic group. The selection of the group of the sibling controls was also based on the criteria of the 10th edition of the (ICD-10). The mean ages and the standard deviations for the dyslexic children and for the controls were 11.47 ± 2.12 (range 7–16 years-old) and 12.21 ± 2.25 years (range 8–16 years-old), respectively (non-significant age difference, t-test p = 0.232). These differences remain non-significant in comparisons between groups separately for boys and girls, as well as between genders separately for the controls and dyslexics. In each case, the following assessments were performed: child psychiatric examination, psychological examination and educational evaluation. The Wechsler Intelligence Scale for Children – Third Edition (WISC-III)  was used to obtain the IQ of each child. The assessment of educational attainment included reading, comprehension, spelling and arithmetic ability. A standard procedure was followed using a special test for the Greek language .
Participants did not enter the study if they had (a) clinically notable neurological disease (including seizure disorder), (b) a history of head injury, (c) hearing difficulties and (d) attention deficit disorder and hyperkinetic syndrome. Prior to participation in the examination, parents were informed about the aims of the study, received a full description of the procedure, and provided written consent. Children were tested individually. The investigators explained to each child the procedure and the children also gave their consent. The study was approved by the local ethical committee.
Stimuli and procedure
The subjects were evaluated with the digit span Wechsler Auditory test [22, 25]. For each trial of the experiment, rest EEG signal was recorded for 500 msec. A single sound tone of either high (3000 Hz) or low frequency (500 Hz) was presented to the subjects through earphones, followed by the numbers which had to be memorized. If the frequency of the signal tone was low the subjects had to recall the numbers in the same order as that presented, else (high frequency tone) the subjects had to recall the numbers in the reverse order. The total task consisted of 52 repetitions for a period of about 45 min. An outline of the procedure is provided in Figure 1.
The children's EEG/ERP signals were recorded at 15 electrodes (Fp1, F3, C5, C3, Fp2, F4, C6, C4, O1, O2, P4, P3, Pz, Cz, Fz) according to the 10–20 international system, referred to both earlobes. Ag/AgCl electrodes were attached to the scalp with adhesive cream in order keep the electrode resistance below 5 kΩ. An electrode placed on the subject's forehead served as ground. The pass band of the amplifiers was set from 0.05 Hz to 35 Hz. During the recordings the subjects had their eyes closed in order to minimize eye movements and blinks. Eye movements were recorded through electro-oculogram (EOG) and recordings with EOG higher than 75 μV were rejected. Warning stimuli, as well as learning material (i.e. the numbers to recall), were presented binaurally via earphones at an intensity of 65 dB sound pressure level. The evoked biosignal was submitted to an analog-to-digital conversion, at a sampling rate of 1 kHz, and was averaged by a computerized system. An algorithm was used, which identified the N100 as the most negative peak in averaged lead curve, between 70 and 150 msec, after the warning stimulus. Peak amplitudes were measured relative to the mean amplitude of the 100-ms pre-stimulus baseline period and latency measurements were computed relative to stimulus onset (Figure 2).
As the N100 component is included in the array of early endogenous ERP components, which normally are modality specific, the ERP induced by the two modal stimuli were averaged separately . In particular, two varieties of N100 waveforms were obtained, one (low N100) evoked by the low frequency modality (26 trials) for each lead in all participants and another (high N100) evoked by the high frequency modality (26 trials) for each lead in all participants.
The behavioral performance refers to the number of correctly recalled digits. The total digits presented in each session were 298, 149 digits for the part of the session which engaged the low frequency stimuli and likewise 149 for the part which engaged the high frequency stimuli.
Preliminary analysis of the data showed that all the dependent variables did not deviate from the normal distribution. Furthermore, the amplitudes and the latencies were found to follow the multivariate normal distribution. This allowed for the application of parametric tests and specifically of multivariate parametric procedures.
The high and low frequency amplitudes and latencies of the N100 component, taken over the range of 70–150 msec, were separately subjected to multivariate analysis with age as covariate and group as the between-subjects factor, followed by multiple pairwise comparisons between the two groups with Bonferroni correction. Using age as a covariate controls for its effect on the ERP component. After finding the locations, at which there were statistically significant differences between the two groups, the leads were sorted in a descending order of the significance of the difference and stepdown procedures were applied. The purpose was to elucidate whether significant differences at more than one lead were due to high correlations between the lead amplitudes and/or latencies or were specific and independent for each lead. IQ scores were compared between the two groups with the t-test for independent samples. Finally, the association between the N100 amplitudes at the high and low frequency tone with memory performance for the two groups was analyzed using partial correlation, controlling for age. Differences between the correlation coefficients for the two groups were tested using the appropriate Fisher transformations. Statistical significance was set at the 0.05 level in all cases.
Amplitudes and latencies
No differences were observed for the mean latency values of the N100 component between the control and the dyslexics group, neither at the high nor at the low frequency tone. Conversely, as figure 3 shows, the mean absolute amplitudes of the N100 component for the control group was greater than for the dyslexics group for all leads, both at the high and at the low frequency tone. As post-hoc comparisons with Bonferroni correction showed, at the majority of the leads (7 out of 15 at high frequency, 9 out of 15 at low frequency), these differences achieved statistical significance. These results are shown in table 1 and table 2 respectively, where asterisks denote the leads, at which the mean absolute amplitude value for the control group is significantly greater than for the dyslexics group (one asterisk: p < 0.05, two asterisks: p < 0.01).
The dyslexics group can be divided into two subgroups: those children, whose siblings were also examined and dyslexics without normal siblings. Statistical analysis showed that these two subgroups were practically the same and differed in the same manner from the control group.
The N100 ERP component achieves its greater mean absolute amplitudes in frontal, central and parietal regions (Cz, C4, Fz, C6, P4, C3, Pz) for controls at both frequencies, in central, parietal regions (Cz, C4, P4, Pz, Fz, C3, C6) for dyslexics at high frequency and in central, parietal and occipital regions (Cz, Pz, P4, C4, P3, O2, O1) for dyslexics at low frequency stimulus type as it can be seen in Figure 3.
In both stimulus frequency types the lead with the maximal difference between the two groups was Fp2. Stepdown procedures proved that group differences found at the other leads could be attributed to the high covariance between the leads. This means that the leads at the N100 component act as an assembly, behaving in a coherent manner. Therefore dissimilar behaviour of the two groups is expected to be observed simultaneously over a large area of the scalp.
No differences were observed between the compared groups with regard to IQ scores (dyslexics: 90.74 ± 11.86, controls: 93.60 ± 10.11, p-value: 0.420). Similarly, with regard to the memory performance, the main finding of the study is that dyslexic children did not differ in comparison to controls for both stimulus types. Better performance was observed for all subjects in forward digit span (low frequency stimulus type) than reverse. In high frequency (reverse digit span) where the mental task is more demanding, difference in memory performance between controls and dyslexics is more intense (t-test, p = 0.084) but it does not achieve statistical significance. These can be illustrated in figure 4.
Correlations of N100 with memory
Table 3 shows the partial correlation coefficients for the two groups between N100 amplitudes at the high and low frequency tone with memory performance controlling for age.
Enhanced memory performance is related with increased N100 amplitudes (more negative value) for the control group and mainly in high frequency stimulus type. This relation is especially conspicuous for the amplitudes of the P4 electrode. Using the appropriate Fisher transformations, it was found that in high frequency stimulus type significant differences between the above correlation coefficients for the controls and dyslexics appear in electrodes C3, P4 and in low frequency stimulus type only in electrode P4.
The results showed that the dyslexic group demonstrated significantly reduced amplitudes of the N100 induced by the high frequency stimulus at Fp1, F3, Fp2, Fz, C4, Cz and F4 leads as well as significantly reduced amplitudes of the N100 induced by the low frequency stimulus at Fp1, F3, C5, C3, Fz, F4, C6, P4 and Fp2 leads. As far as memory performance is concerned, the main finding of the study is that dyslexic children did not fall short in a significant manner from their sibling controls,. However, enhanced memory performance in the control group is associated with increased N100 amplitude induced by high frequency stimuli at C5, C3, C6 and P4 leads and increased N100 amplitude induced by low frequency stimuli at P4 lead.
The patterns of the amplitudes of the N100 induced by both high and low frequency stimuli obtained in the present study suggest that children with dyslexia may demonstrate impairments in 'attentional operation' of information processing (as they are reflected by the N100 amplitudes elicited during a WM test) that involve or affect widespread brain areas. This suggestion seems to be in agreement with evidence elicited from recent research, demonstrating the implication of attentional alterations in the pathophysiology of dyslexia, as it is reflected by the N100 component [14, 16, 27].
In particular, Bonte and Blomert , investigating ERP correlates of implicit phonological processing during spoken word recognition in dyslexic and normally reading children (7–10 years), observed that the dyslexics manifested smaller N100 amplitudes at temporal electrodes than normal readers. These findings have been conceived of as indices of 'functional and/or anatomical anomalies in neural sources involved in speech processing in dyslexics'. Similarly, Pinkerton  found reduction of the N100 in a group of 14 boys with difficulties in reading, writing and spelling as compared to those of 18 matched controls (8–9 years). Interpreting these findings, the authors hypothesized that the reduced N100 in dyslexics might be understood as an index of impaired selective attention. A very recent study by Fosker and Thierry  investigated 12 dyslexic adults (mean age 20 ± 1 year) and 12 control adults (mean age 19 ± 1 year) using an oddball paradigm, reported that the dyslexic adults showed significantly smaller amplitudes of N100, concluding that these findings point to a deficit on the low-level perceptual processing capacity in association with attention. In corroboration to this idea seems reasonable to bring in mind that the N100 component can be elicited by 'pre-attentive' means .
Brunswick and Rippon  also found that normal controls exhibited significantly greater N100 amplitudes at the left temporal region during a dichotic listening test, while the dyslexic children exhibited approximately equivalent values of N100 amplitudes bilaterally. The authors regarded the divergence in N100 laterality as an indication of "inaccurate tuning of sensory information" leading to less reliable auditory representation.
Taken together – despite the use of divergent methodological approaches and eliciting procedures – these reports show that dyslexics demonstrate significantly smaller N100 amplitudes as compared to controls, hence indicating that dyslexia may be associated with abnormal processes mediating attentional operation. With regard to the locations that dyslexics demonstrate decreased N100 amplitudes, it is obvious that it concerns a widespread network within various brain regions.
As far as memory performance is concerned, the main finding of the study is that dyslexic children did not differ significantly in comparison to normal controls. This finding is in agreement with other related studies . However a very large body of evidence points to the presence of a deficit in phonological short term memory in developmental dyslexia. [30–32]. Sources for divergent findings may be conceived in terms of both the heterogeneity of the study samples and differences in the measurement of the procedures. However, this still leaves open the possibility that genetic factors, may contribute to the lack of significant difference in memory performance observed in the present study. Given that developmental dyslexia is genetically associated [33, 34] the possibility that the sibling controls of the current study were slightly dyslexic themselves cannot be completely ruled out. In this sense, the reports [35, 36] provide additional explanatory hints underscoring the importance of the maturational age trajectory.
Finally, the observed association between the amplitude of the N100 and memory performance in the control group, in contrast to the lack of such an association in dyslexic children could be explained taking into account the anchoring-deficit hypothesis, which proposes that the deficit of dyslexics resides in the dynamics that link attention with short memory through the implicit formation of stimulus-specific anchors .
Certain limitations of this investigation warrant consideration. Firstly, the main findings need to be replicated in independent samples including unrelated controls and it is to be determined, whether there is an association in a task-specific manner or across tasks. Secondly, beyond of the advantage which incurs of the use of unaffected sibling controls it should be taken into account that possible confounding parameters might be implicated through this fact, such as the common genetic factors with their possible consequences. Hence future studies controlling for genetic factors as well as trait and state parameters in combination with more experiments that combine the time resolution of event-related potentials with the spatial resolution of brain imaging techniques may lead to clearer definitions of the brain function in relation to the current findings.
The present findings provide evidence that dyslexic children may manifest impaired 'attentional operation' of information processing as reflected by N100 elicited during a working memory (WM) test. In addition these results acknowledge a possible impairment in the dynamics that link attention with short memory as suggested by the anchoring-deficit hypothesis.
Shaywitz S: Current concepts: Dyslexia. N Engl J Med. 1998, 338 (5): 307-312. 10.1056/NEJM199801293380507.
Shaywitz SE: Management of dyslexia, its rationale, and underlying neurobiology. Pediatr Clin North Am. 2007, 54 (3): 609-623. 10.1016/j.pcl.2007.02.013.
Montgomery JW: Working memory and comprehension in children with specific language impairment: what we know so far. J Commun Disord. 2003, 36: 221-231. 10.1016/S0021-9924(03)00021-2.
Toplak ME, Rucklidge JJ, Hetherington R, John SC, Tannock R: Time perception deficits in attention-deficit/hyperactivity disorder and comorbid reading difficulties in child and adolescent samples. J Child Psychol Psychiatry. 2003, 44 (6): 888-903. 10.1111/1469-7610.00173.
Webster RI, Shevell MI: Neurobiology of specific language impairment. J Child Neurol. 2004, 19: 471-481. 10.1177/08830738040190070101.
Fitch RH, Tallal P: Neural mechanisms of language-based learning impairments: insights from human populations and animal models. Behav Cogn Neurosci Rev. 2003, 2 (3): 155-178. 10.1177/1534582303258736.
Vellutino FR, Fletcher JM, Snowling MJ, Scanlon DM: Specific reading disability (dyslexia): what have we learned in the past four decades?. J Child Psychol Psychiatry. 2004, 45: 2-40. 10.1046/j.0021-9630.2003.00305.x.
Shaywitz SE, Shaywitz BA: Dyslexia (specific reading disability). Biol Psychiatry. 2005, 57 (11): 1301-1309. 10.1016/j.biopsych.2005.01.043.
Hillyard SA, Hink RF, Schwent VL, Picton TW: Electrical signs of selective attention in the human brain. Science. 1973, 182 (108): 177-180. 10.1126/science.182.4108.177.
Mulert C, Menzinger E, Leicht G, Pogarell O, Hegerl U: Evidence for a close relationship between conscious effort and anterior cingulate cortex activity. Int J Psychophysiol. 2005, 56 (1): 65-80. 10.1016/j.ijpsycho.2004.10.002.
Mulert C, Seifert C, Leicht G, Kirsch V, Ertl M, Karch S, Moosmann M, Lutz J, Möller HJ, Hegerl U, Pogarell O, Jäger L: Single-trial coupling of EEG and fMRI reveals the involvement of early anterior cingulate cortex activation in effortful decision making. Neuroimage. 2008, 42 (1): 158-68. 10.1016/j.neuroimage.2008.04.236.
Coull JT: Neural correlates of attention and arousal: insights from electrophysiology, functional neuroimaging and psychopharmacology. Prog Neurobiol. 1998, 55: 343-361. 10.1016/S0301-0082(98)00011-2.
Golob EJ, Starr A: Serial position effects in auditory event-related potentials during working memory retrieval. J Cogn Neurosci. 2004, 16: 40-52. 10.1162/089892904322755548.
Pinkerton F, Watson DR, McClelland RJ: A neurophysiological study of children with reading, writing and spelling difficulties. Dev Med Child Neurol. 1989, 31: 569-581.
Brunswick N, Rippon G: uditory event-related potentials, dichotic listening performance and handedness as indices of lateralization in dyslexic and normal readers. Int J Psychophysiol. 1994, 18: 265-275. 10.1016/0167-8760(94)90012-4.
Bonte ML, Blomert L: Developmental dyslexia: ERP correlates of anomalous phonological processing during spoken word recognition. Brain Res Cogn Brain Res. 2004, 21: 360-376. 10.1016/j.cogbrainres.2004.06.010.
Molfese DL: Predicting dyslexia at 8 years of age using neonatal brain responses. Brain Lang. 2000, 72: 238-245. 10.1006/brln.2000.2287.
Heim S, Keil A: Large-scale neural correlates of developmental dyslexia. Eur Child Adolesc Psychiatry. 2004, 13: 125-140. 10.1007/s00787-004-0361-7.
Neville HJ, Coffey SA, Holcomb PJ, Talla P: The neurobiology of sensory and language processing in language-impaired children. J Cogn Neurosci. 1993, 5: 235-253. 10.1162/jocn.1918.104.22.168.
Yingling CD, Galin D, Fein G, Peltzman D, Davenport L: Neurometrics does not detect 'pure' dyslexics. Electroencephalogr Clin Neurophysiol. 1986, 63: 426-430. 10.1016/0013-4694(86)90124-0.
Putter-Katz H, Kishon-Rabin L, Sachartov E, Shabtai EL, Sadeh M, Weiz R, Gadoth N, Pratt H: Cortical activity of children with dyslexia during natural speech processing: evidence of auditory processing deficiency. J Basic Clin Physiol Pharmacol. 2005, 16: 157-171.
Papageorgiou C, Anagnostopoulos D, Giannakakis GA, Sakelariou K, Tsiaparas N, Paraskevopoulou P, Nikita KS, Rabavilas A, Soldatos C: Preattentive deficits in developmental disorders of scholastic skills. Neuroreport. 2005, 16: 1829-1832. 10.1097/01.wnr.0000185963.96217.7f.
Watkins MW, Worrell FC: Diagnostic utility of the number of WISC-III subtests deviating from mean performance among students with learning disabilities. Psychol Schools. 2000, 37: 303-309. 10.1002/1520-6807(20002)37:4<311::AID-PITS1>3.0.CO;2-X.
Paraskevopoulos IN: 'Athina' test for diagnosis of learning disabilities (in Greek). Edited by: Helleneka Grammata. 1996, Athens: Faculty of Philosophy. University of Athens
Wechsler D: Manual for the Wechsler adult intelligence scale. 1955, New York: Psychological Corporation
Fabiani M, Gratton G, Coles M: Event-related potentials: methods, theory, and applications. 2000, New York: Cambridge University Press
Fosker T, Thierry G: Phonological oddballs in the focus of attention elicit a normal P3b in dyslexic adults. Brain Res Cogn Brain Res. 2005, 24: 467-475. 10.1016/j.cogbrainres.2005.02.019.
Kraus N, McGee T, Carrell TD, Sharma A: Neurophysiologic bases of speech discrimination. Ear Hear. 1995, 16 (1): 19-37. 10.1097/00003446-199502000-00003.
McNeil AM, Johnston RS: Word length, phonemic, and visual similarity effects in poor and normal readers. Mem Cognit. 2004, 32 (5): 687-95.
Ramus F, Szenkovits G: What phonological deficit?. J Exp Psychol (Colchester). 2008, 61 (1): 129-41. 10.1080/17470210701508822.
Kibby MY, Kroese JM, Morgan AE, Hiemenz JR, Cohen MJ, Hynd GW: The relationship between perisylvian morphology and verbal short-term memory functioning in children with neurodevelopmental disorders. Brain Lang. 2004, 89: 122-135. 10.1016/S0093-934X(03)00310-9.
Ahissar M: Dyslexia and the anchoring-deficit hypothesis. Trends Cogn Sci. 2007, 11: 458-465.
Shastry BS: Developmental dyslexia: an update. J Hum Genet. 2007, 52 (2): 104-109. 10.1007/s10038-006-0088-z.
Gibson CJ, Gruen JR: The human lexinome: genes of language and reading. J Commun Disord. 2008, 41 (5): 409-20. 10.1016/j.jcomdis.2008.03.003.
Dawes P, Bishop DV: Maturation of visual and auditory temporal processing in school-aged children. J Speech Lang Hear Res. 2008, 51 (4): 1002-1015. 10.1044/1092-4388(2008/073).
Brunner M, Baeumer C, Dockter S, Feldhusen F, Plinkert P, Proeschel U: Heidelberg Phoneme Discrimination Test (HLAD): normative data for children of the third grade and correlation with spelling ability. Folia Phoniatr Logop. 2008, 60 (3): 157-161. 10.1159/000121393.
This work has been funded by the project PENED 2003, Ministry of Development – General Secretariat of Research and Technology. The authors gratefully thank M. Kyprianou, Scientific Investigator, for his support on the statistical analysis of the experimental results. The authors would also thank K. Sakelariou, P. Paraskevopoulou, psychologists, for helpful comments and two anonymous reviewers for useful comments on the manuscript.
The authors declare that they have no competing interests.
CP participated in the design of the study, designed the experimental procedure and the clinical protocol, participated in data acquisition, participated in analysis and interpretation of data and drafted the manuscript. GAG participated in the design of the study, participated in acquisition, analysis and interpretation of data and drafted the manuscript. KSN conceived and coordinated this study, participated in the interpretation of data and corrected the final version of this manuscript. DA participated and conceived in the design of the study, participated in the interpretation of data. GNP participated and conceived in the design of the study, participated in the interpretation of data. AR: conceived and coordinated this study and corrected the final version of this manuscript. All authors have read and approved the final manuscript.