Location via proxy:   [ UP ]  
[Report a bug]   [Manage cookies]                
Next Article in Journal
Improving Outcomes in People with Spinal Cord Injury: Encouraging Results from a Multidisciplinary Advanced Rehabilitation Pathway
Next Article in Special Issue
Social Brain Perspectives on the Social and Evolutionary Neuroscience of Human Language
Previous Article in Journal
Overview of the Current Knowledge and Conventional MRI Characteristics of Peri- and Para-Vascular Spaces
Previous Article in Special Issue
The Relationship between Language and Technology: How Screen Time Affects Language Development in Early Life—A Systematic Review
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Developmental Dyslexia: Insights from EEG-Based Findings and Molecular Signatures—A Pilot Study

by
Daniela Theodoridou
1,
Christos-Orestis Tsiantis
1,
Angeliki-Maria Vlaikou
2,3,
Vasiliki Chondrou
4,
Victoria Zakopoulou
5,
Pavlos Christodoulides
5,6,
Emmanouil D. Oikonomou
7,
Katerina D. Tzimourta
8,
Charilaos Kostoulas
9,
Alexandros T. Tzallas
7,
Konstantinos I. Tsamis
6,
Dimitrios Peschos
6,
Argyro Sgourou
4,
Michaela D. Filiou
2,3 and
Maria Syrrou
1,*
1
Laboratory of Biology, Faculty of Medicine, School of Health Sciences, University of Ioannina, 45110 Ioannina, Greece
2
Biomedical Research Institute, Foundation for Research and Technology-Hellas (FORTH), 45110 Ioannina, Greece
3
Laboratory of Biochemistry, Department of Biological Applications and Technology, School of Health Sciences, University of Ioannina, 45110 Ioannina, Greece
4
Laboratory of Biology, School of Science and Technology, Hellenic Open University, 26335 Patras, Greece
5
Department of Speech and Language Therapy, School of Health Sciences, University of Ioannina, 45110 Ioannina, Greece
6
Laboratory of Physiology, Faculty of Medicine, School of Health Sciences, University of Ioannina, 45110 Ioannina, Greece
7
Department of Informatics and Telecommunications, School of Informatics & Telecommunications, University of Ioannina, 47100 Arta, Greece
8
Department of Electrical and Computer Engineering, University of Western Macedonia, 50100 Kozani, Greece
9
Laboratory of Medical Genetics, Faculty of Medicine, School of Health Sciences, University of Ioannina, 45110 Ioannina, Greece
*
Author to whom correspondence should be addressed.
Brain Sci. 2024, 14(2), 139; https://doi.org/10.3390/brainsci14020139
Submission received: 23 November 2023 / Revised: 25 January 2024 / Accepted: 26 January 2024 / Published: 28 January 2024
(This article belongs to the Special Issue Neurodevelopmental Disorders and Early Language Acquisition)

Abstract

:
Developmental dyslexia (DD) is a learning disorder. Although risk genes have been identified, environmental factors, and particularly stress arising from constant difficulties, have been associated with the occurrence of DD by affecting brain plasticity and function, especially during critical neurodevelopmental stages. In this work, electroencephalogram (EEG) findings were coupled with the genetic and epigenetic molecular signatures of individuals with DD and matched controls. Specifically, we investigated the genetic and epigenetic correlates of key stress-associated genes (NR3C1, NR3C2, FKBP5, GILZ, SLC6A4) with psychological characteristics (depression, anxiety, and stress) often included in DD diagnostic criteria, as well as with brain EEG findings. We paired the observed brain rhythms with the expression levels of stress-related genes, investigated the epigenetic profile of the stress regulator glucocorticoid receptor (GR) and correlated such indices with demographic findings. This study presents a new interdisciplinary approach and findings that support the idea that stress, attributed to the demands of the school environment, may act as a contributing factor in the occurrence of the DD phenotype.

1. Introduction

Developmental dyslexia (DD) (Although Specific Learning Disorder is the generic term used in the DSM-5 to describe disorders characterized by difficulties in learning and academic skills, we adopt the term developmental dyslexia (DD) as it is an impairment that affects the development of decoding skills in reading, which is the core focus of this study, differentiating it from the other learning disorders of dyscalculia and disorders of written expression [1]) is a phenotypically and genetically heterogeneous, multifactorial, complex condition characterized by difficulties in reading, writing and spelling, which affects 5–20% of the general population [2] and 5–17% of school-age children [3]. Adults with DD present phonological processing impairments [4,5], working memory deficits [6], primary auditory processing abnormalities, impaired visual processing [7,8], attentional deficits [9] and irregularities in the processing of written information [10]. These difficulties, usually intertwined with a constant pattern of failure, provide strong conditions for emotional and psychological issues to arise (low self-esteem, anxiety, feelings of isolation) [11]. Although diagnosed in childhood, DD is a lifelong condition, possibly intensified by adverse experiences that occur during the lifetime, with putative additive and reinforcing effects on the phenotype [12,13,14,15].
Several genes and gene variants have been associated with DD predisposition, but, as in all complex conditions, the phenotypic and genotypic heterogeneity is the result of the interplay between genetic and non-genetic environmental factors, in the form of epigenetically encoded information (gene × environment interactions, G × E) [16,17,18,19,20,21,22,23]. According to recent hypotheses and research findings, among environmental factors, stress is widely documented as a contributing factor in the etiopathogenesis of complex disorders, due to its effects on brain plasticity and function, particularly during early developmental stages that are critical for neurodevelopment [24,25,26,27].
In response to a stressor, hypothalamic–pituitary–adrenal (HPA) axis-related genes and especially glucocorticoid signaling pathway genes mediate the stress response [28]. Glucocorticoid receptor (GR), encoded by the Nuclear Receptor Subfamily 3 Group C Member 1 (NR3C1) gene, and mineralocorticoid receptor (MR), encoded by the Nuclear Receptor Subfamily 3 Group C Member 2 (NR3C2) gene, mediate the physiological stress response and the organism’s homeostasis. NR3C1 expression levels are altered in anxiety and mood disorders and alterations have been already implicated in the DD phenotype in humans [29,30]. A methylation analysis of NR3C1 and its alternative, non-coding first exons revealed a strong association between stress and methylation levels and the occurrence of depressive episodes after prenatal adversity and adult adverse life experiences [31,32,33]. Glucocorticoid-induced leucine zipper (GILZ) is a GR-induced measure of the stress response that is transcriptionally activated by glucocorticoids and links the neuroendocrine stress response to immune-related disorders [34]. FK506 Binding Protein 5 (FKBP5) is a co-chaperone of GR, the levels of which were found altered in complex stress-associated phenotypes [35]. The functional single nucleotide variant (SNV) rs1360780 has been implicated in the occurrence of complex phenotypes and stress-associated psychopathology, with the (T) allele considered as a risk allele [26,36,37]. Solute Carrier Family 6 Member 4 (SLC6A4) encodes the serotonin transporter gene. The polymorphic 5-HTTLPR region located in its promoter, which also encompasses rs25531, has been linked to an increased risk of anxiety and depression-related disorders [38,39].
Glucocorticoids have been shown to regulate mitochondrial function [40] and GR has been reported to regulate mitochondrial DNA (mtDNA) gene expression [41]. Mitochondria have recently emerged as pathophysiology hubs and potential therapeutic targets for stress-related and anxiety disorders [42,43,44,45,46,47,48].
The current diagnostic approach for DD and specific learning disorders employs standardized clinical tests and tools to assess phonological skills, reading and spelling, while neuroimaging modalities are only used to exclude other possible diagnoses [49,50]. Diagnostic tools with high accuracy in internalizing disorders are also widely used in everyday clinical practice to test the association between reading problems and anxiety or depression [51,52]. Anxiety disorders have been repeatedly associated with alterations in the electroencephalography (EEG)-recorded neuronal network oscillations [53]. Research findings, on the other hand, focus on abnormalities observed on magnetic resonance imaging (MRI), EEG, magnetoencephalography (MEG) and biometric signals, such as heart rate, electrodermal activity (EDA) and eye tracking [54,55,56,57], which enhance our understanding of the pathophysiology and could be used for DD diagnosis and management [58,59,60,61].
EEG is a low cost, non-invasive method of recording the electrical activity on the scalp, which represents the activity of the cerebral cortex [62]. A set of distinct factors extracted from the spectral analysis of EEG signals that can provide bio-evidence of dyslexia was proposed [63]. Accordingly, in a recent report studying the functional brain connections in children with and without dyslexia using resting-state EEG data, brain network deficiencies were observed when performing various spelling, phonological and morphological awareness tasks [64]. Synchronized neuronal excitation produces explicit rhythms on EEG waves that have been related to distinct cognitive states and to brain disorders [65,66]. In studies focusing on EEG rhythms’ changes, increased frontal theta and beta rhythm activity was reported in the right parietal–occipital area in children with dyslexia when performing phonological tasks, while other studies have reported greater consistency in delta, theta and beta rhythms and lower consistency in alpha frequencies during relaxation [67,68,69,70].
In adults with dyslexia, research findings have underscored low brain activation in a variety of brain regions, such as reduced activation of the left inferior angular gyrus during phonological tasks [71,72], reduced activation in the frontal network of the left hemisphere associated with phonological working memory [73] and increased activation of the left inferior frontal gyrus during phonological tests [74].
The diagnostic coupling of the EEG method with machine learning techniques has given new insights into the neuronal network deficiencies underlying specific learning disorders, providing a wealth of information for their diagnosis and pathophysiology [64,75,76,77]. Despite the absence of a consensus on the methodology for implementation in the diagnosis [60,70,78], it has been shown that recently established EEG research protocols have diagnostic relevance in DD [56,79], contributing to the classification [80] and better understanding of brain function [81]. Such techniques were also applied in a study including children with and without dyslexia, which concluded that it is possible to identify underlying neurocognitive profiles for dyslexia and therefore to distinguish between the two populations through the information extracted from EEG [82]. In studies using algorithms to interpret EEG data, the classification accuracy reached 90% (study at rest) and 100% (study during testing), respectively [83,84].
Moreover, EEG data are also used by researchers (a) as part of multifactorial approaches to identify dyslexia [85], (b) as a means of recording emotions to analyze symptoms of both dyslexia and other disorders such as autism and ADHD [86], (c) as a means of comparing brain function between “competent” children with dyslexia (capable dyslexic) and severe cases of dyslexia (poor dyslexic) [87] and (d) as a means of recording changes in cerebral laterality after an intervention period in children with dyslexia [88].
The search for reliable biomarkers in medicine is a crucial issue for the development of effective personalized treatment strategies [89]. Although validating biomarkers for neurodevelopmental conditions carries difficulties, primarily due to the heterogeneity of these complex conditions, the use of the most widely accepted modalities to date—namely non-invasive neuroimaging, genetic and epigenetic data and the correlation of the findings—emerges as a promising strategy to advance the field [90]. EEG parameters have been combined with genetic and epigenetic findings of stress-related genes to characterize some neurological conditions. In a study, the FKBP5 rs1360780 and FKBP5 promoter methylation status was combined with EEG parameters in Alzheimer’s disease patients [91], while, in another study, EEG was used to explore the contribution of MR gene variants in habit learning under stress [92].
This approach necessitates a multidisciplinary teamwork, encompassing the co-evaluation of established diagnostic criteria complemented by neuropsychological assessment and advanced neuroimaging techniques, as well as evidence from molecular biomarkers. This multidisciplinary study correlates the already established diagnostic criteria for DD with EEG recordings and stress-associated molecular signatures.
The DD phenotype is not solely attributed to susceptibility genes, which have been identified but are not considered the major etiologic factor. The phenotype is the result of the interaction of environmental factors (risk or protective), along with a genetic predisposition, but the relative impact of each component is not easily traceable [13,93,94]. Stress appears often in the literature as an important factor in the context of DD, and the stress experienced by DD individuals may compromise the negative feedback mechanisms, thereby leading to HPA axis dysregulation [16,17,19,27].
As in all complex, heterogeneous phenotypes, the DD phenotype results from the interplay of genetic and environmental factors and the impact of each factor varies. The DD etiology might include a major- and a minor-impact component. According to our hypothesis, there are “stress-related DD phenotypes”, namely phenotypes where key HPA axis genes may interact in gene networks with DD susceptibility genes, and the major contributing and triggering factor is stress caused by the constant difficulties encountered and the struggles of DD individuals in everyday life and especially in the school environment. A gene network bioinformatic analysis of HPA axis genes and dyslexia susceptibility genes is performed and interactions between the two groups of genes are presented. This is a multidisciplinary pilot study that focuses on DD and stress and proposes an innovative procedure and protocol. In this study genetic, epigenetic and mitochondrial molecular findings are co-evaluated with EEG findings and DASS21 scores, in order to test the possible implications of the HPA axis in the coordination of a complex regulatory network of stress associated with DD.

2. Materials and Methods

2.1. Participants

The dyslexia group (DD) consisted of ten individuals (4 males, 6 females) with an average age of 21.3 years (mean + SD/21.3 + 2.1), who underwent a school-age intervention, and a control group of ten matched controls (2 males, 8 females) with an average age of 20.7 years (20.7 ± 1.2) (control). The sample participants were higher-education students, recruited from the Department of Speech & Language Therapy of the University of Ioannina, Greece (data published in [95]). The recruitment procedure regarding individuals with DD was based on their official diagnosis meeting the Diagnostic and Statistical Manual of Mental Disorders, Fifth Edition (DSM-5) [96]. All individuals who participated in this study provided written informed consent, in accordance with the Declaration of Helsinki (1975, revised 2013). They also consented to the analysis, processing and publication of the data. The study was approved by the Research Ethics and Conduct committee of the University of Ioannina, with issued protocol number 28622 (23 July 2019).

2.2. DASS21 Scale

DD and control group individuals were assessed with the Depression, Anxiety and Stress Scale—21 Items (DASS21) [96] for the following phenotypic traits: depression (dysphoria, hopelessness, devaluation of life, self-depreciation, lack of interest, anhedonia, inertia), anxiety (autonomic arousal, skeletal muscle effects, situational anxiety, subjective experience of anxious effect) and stress (levels of chronic non-specific arousal). The DASS21 was constructed not merely as another set of scales to measure conventionally defined emotional states, but to further advance the process of defining, understanding and measuring the ubiquitous and clinically significant emotional states usually described as depression, anxiety and stress [96,97]. The DASS21 has been standardized for the Greek population with an internal consistency reliability level of α = 0.965 [98].

2.3. Peripheral Blood Sample Collection, DNA and Total RNA Extraction

2.3.1. Peripheral Blood Sample Collection

Peripheral blood was used as it is the most widely studied and a relatively and easily accessible tissue [99]. All participants provided peripheral blood samples, after the provision of informed consent documents, in accordance with the principles of the Helsinki Declaration [100]. Per participant, 3 mL of peripheral blood was collected, and DNA and RNA were immediately extracted.

2.3.2. DNA Extraction

Genomic DNA was extracted from 1 mL whole blood with the NucleoSpin Blood Kit (Macherey-Nagel, Düren, Germany, #740951). DNA was further processed for the genotyping of SNVs of interest, namely rs1360780 (FKBP5), 5-HTTLPR and rs25531 (SCL6A4), as well as for the investigation of the methylation status of the NR3C1 gene and the evaluation of the mitochondrial DNA copy number (mtDNAcn).

2.3.3. Total RNA Extraction

Per participant, 2 mL samples of peripheral blood were layered on 1 mL of Lymphosep, Lymphocyte Separation Media (Biosera, Nuaille, Pays de la Loire, France, #LM-T1702), and peripheral blood mononuclear cells (PBMCs) were isolated following consecutive centrifugations (1800 rpm/20 min, 1600 rpm/10 min, 1400 rpm/10 min at 4 °C). The PBMC precipitate was resuspended in 500 μL of NucleoZol (Macherey-Nagel, Düren, Germany, #740404). The Nucleospin RNA set for Nucleozol (Macherey-Nagel, #740406) was used for total RNA extraction, following the manufacturer’s instructions. cDNA synthesis by reverse transcription was performed using 0.5 μg of total RNA, according to the standard protocol of the PrimeScript RT Reagent Kit, with the gDNA Eraser (Takara, Shiga, Japan, #RR047A).

2.4. Stress-Associated SNV Detection

Genotyping for the rs1360780 (T > C) SNV of FKBP5 was performed utilizing a TaqMan assay (Thermo Assay ID: C_8852038_10, #4351379, Applied Biosystems, Foster City, CA), based on the manufacturer’s instructions, and genotypes were identified with the StepOne Software v2.3 (Applied Biosystems, Waltham, MA, USA). Analysis of the polymorphic region 5-HTTLPR, characterized by a 44 bp insertion/deletion, and the SNV rs25531 (A > G), both located in the Solute Carrier Family 6 Member 4 (SLC6A4) gene, was carried out by a polymerase chain reaction–restriction fragment length polymorphism (PCR-RFLP) assay. PCR primer sequences were as follows: Forward: 5′-GGCGTTGCCGCTCTGAATGC-3′ and Reverse: 5′-GAGGGACTGAGCTGGACAACCAC-3′. Cycling conditions were as previously described [101]. The expected PCR amplicons were 529 bp and 486 bp, which were therefore named the Long (L) and Short (S) allele, respectively. PCR amplicons were digested with the MspI restriction enzyme (Takara, Otsu, Japan, #1150A) for the identification of the SNV rs25531, and, depending on the presence of the L or S allele, the digestion produced a discrete pattern of fragments, as previously described [101,102]. Genotyping was not possible for two control samples due to an insufficient amount of DNA.

2.5. Evaluation of mRNA Expression Levels of Stress-Related Genes

mRNA levels were quantified by real-time qPCR using the CFX Manager 3.1 software (Bio-Rad, Hercules, CA, USA). For each gene of interest (NR3C1, NR3C2, GILZ, FKBP5), primers were designed according to the corresponding NCBI GenBank sequences and selected using the NCBI primer design tool (Supplementary Table S1). Reactions were performed with the Universal SYBR Green Supermix (Kapa Biosystems, Wilmington, MA, USA, KK 4602), 400 nM of each set of primers and 20 ng of template cDNA. PCR cycling conditions were 94 °C for 3 min, followed by 35 cycles of 94 °C for 45 s 59 °C for 30 s, 72 °C for 1 min, and the final step at 72 °C for 10 min. All qPCR experiments were performed in duplicate. The Ct values of the GR (NR3C1), MR (NR3C2), GILZ (TSC22D3) and FKBP5 target genes were normalized to the Ct values of β-actin as a reference gene. Datasets were analyzed using the 2−ΔΔCt method and expressed as relative gene expression [103].

2.6. NR3C1 Gene Methylation Profile

2.6.1. Bioinformatic Analysis

The human NR3C1 gene is composed of eight coding exons and nine alternative, untranslated first exons. Seven of them (1D, 1J, 1E, 1B, 1F, 1C, 1H) lie lengthwise within NR3C1 CpG island 303, being located in the proximal promoter region. CpG island 303 was retrieved by the UCSC genome browser (genome.ucsc.edu) [104].

2.6.2. Methylation Analysis

The methylation profile analysis of two NR3C1 alternative non-coding first exons, which are suggested to act as alternative promoters [105], was performed with the Pyrosequencing CpG assay methodology. PCR primers as well as sequencing primers, specific to the two different regions, were designed with the PyroMark Assay Design Software version 2.0 (Qiagen, Venlo, The Netherlands) (Supplementary Table S2). Genomic DNA isolation was followed by the bisulfite conversion of all unmethylated cytosines to uraciles according to the manufacturer’s instructions (Qiagen, #59104). PCR amplification of bisulfite-converted DNA using a reverse biotinylated 5′ end primer was performed with the PyroMark PCR kit, according to the manufacturer’s instructions (Qiagen, #978703). Then, the methylation status of all CpG sites of each region was analyzed with PyroMark Q24 MDx technology (Qiagen, #970802). Methylated and unmethylated (bisulfite-converted) human DNA plus unmethylated (physiological) human DNA were used as internal controls in every round of the pyrosequencing reaction. The bisulfite conversion efficiency was checked using quality controls in every pyrosequencing run set. Variations in the methylation levels of various CpG sites were detected and presented in a sequence context. Two independent pyrosequencing reactions were performed for every DNA sequence tested.

2.7. Mitochondrial DNA Copy Number (mtDNAcn) Estimation in Blood

Blood mtDNAcn was evaluated by qPCR, using the CFX Connect Real-Time PCR Detection System (Bio-Rad). The mtDNAcn content in genomic DNA was quantified by the CFX Manager 3.1 software (Bio-Rad), using specific primers to target the nuclear DNA: GAPDH Forward 5′-GTGGTCTCCTCTGACTTCAACA-3′, GAPDH Reverse 5′-ACCACCCTGTTGCTGTAGCC-3′. In addition, specific primers were used to target the mtDNA [106]: L394 Forward 5′-CACCAGCCTAACCAGATTTC-3′, H475 Reverse 5′-GGGTTGTATTGATGAGATTAGT-3′. The PCR mixture contained 1× Universal SYBR Green Supermix (Kapa Biosystems, Cape Town, South Africa, #KK 4602), 400 nM primers and 10 ng of template genomic DNA. The PCR cycling program was the following: 95 °C for 10 min, followed by 40 cycles of 95 °C for 10 s, 60 °C for 30 s and 72 °C for 30 s. GAPDH was used as a reference housekeeping gene for the calculation of the relative expression values of mtDNAcn in each sample. All qPCR experiments were performed in duplicate. The cycle in which the product was detected was defined as Ct. The Ct values provided from the qPCR for the mtDNA were normalized to the Ct values of the reference gene (GAPDH), and the results were expressed by the 2−ΔΔCt method [103].

2.8. EEG Recordings

EEG signals for each channel and rhythm were recorded with a wearable EEG device, namely the Emotiv EPOC+, following the calculation of the relative band power (RBP). The power and energy spectral analysis of EEG oscillations is a well-established method for the study of cognitive processes [107,108,109,110]. The methodology for data collection and analysis was described in detail in our previous study, where the contribution of EEG recordings to the audio–visual recognition of words was explored in university students with DD [78]. For the present study, the same dataset was used, and the energy of EEG wavelets was correlated to the findings of the molecular analysis for participants with and without DD. The EEG channels that were explored followed the International 10–20 System (AF3, F3, F7, FC5, T7, P7, O1, AF4, F4, F8, FC6, T8, P8 and O2), with the exception of two electrodes that were isolated and rejected due to loss of connectivity (F8 and its corresponding electrode F7) (Supplementary Figure S1, used with permission from [78]).

2.8.1. EEG Data Acquisition

EEG recordings were conducted within a room designed to minimize sound and light, capturing data from each participant during the ongoing evaluation test. Before commencing, a skilled researcher briefed the participants on the experimental process. To familiarize them with the protocol and the equipment, each participant underwent a training session lasting 4–7 min. Throughout the evaluation and EEG recording, participants remained seated upright, in a relaxed state, with their eyes open. The duration of the EEG recording varied, spanning from 21 to 38 min (with an average of 28 min), contingent on individual test completion times. Recordings were halted if any participant experienced discomfort with the procedure or the device. A total of approximately 11.5 h of EEG data were amassed, comprising 5 h and 51 min from dyslexic participants and 5 h and 47 min from non-dyslexic participants. For data collection, the Emotiv EPOC was employed, a commercially available wearable EEG device. This device, known for its lifestyle applications, encompasses 14 sensors with corresponding felt pads positioned on the scalp following the International 10–20 System layout (including AF3, F3, F7, FC5, T7, P7, O1, AF4, F4, F8, FC6, T8, P8 and O2). Additionally, two electrodes were affixed to the mastoids to serve as reference channels. The EEG data were sampled at a frequency of 128 Hz, and electrode–scalp connections were established using a saline liquid solution applied to all sensor felt pads. The setup adhered to the guidelines of the EmotivPRO 3 software, with the connectivity quality routinely monitored at the beginning and throughout the recording process.

2.8.2. Preprocessing

The recordings were conducted using a montage based on the interconnected mastoids. Subsequent to each recording session, the EEG signals were exported in “.edf” format and subjected to processing utilizing the MATLAB platform along with the EEGLAB toolbox. To mitigate interference from 50 Hz power line noise oscillations, a Butterworth notch filter was applied to the EEG signals. Additionally, a high-pass FIR digital filter was implemented at 0.5 Hz to eliminate low-frequency oscillations. Subsequently, a set of five equiripple FIR filters was devised. These filters were designed to permit frequencies within specified ranges while attenuating frequencies outside of these ranges. The design of these five band-pass filters (within the range of 0.5–4 Hz, 4–8 Hz, 8–12 Hz, 13–30 Hz and 30–60 Hz) was tailored to correspond to the distinct EEG rhythms, thereby striving to extract spectral features within each specific frequency sub-band of interest. The EEG frequency bands are described in Table 1. Following this filtering process, each EEG recording was divided into non-overlapping epochs of 10 s, and spectral features were subsequently derived from these individual 10-s EEG segments.

2.8.3. Feature Extraction

From the preprocessed EEG recordings, the RBP was extracted as a feature. To derive the RBP, the power spectral density (PSD) of the signal within distinct frequency bands was acquired by employing the Welch method. This technique involves segmenting the signal into overlapping fragments and computing the squared magnitude of the discrete Fourier transform for each segment. By averaging these computed values, a comprehensive estimate of the PSD is generated. Subsequently, the relative proportion of PSD for each frequency band within each time epoch was computed, resulting in a feature matrix comprising five attributes for each row. To determine the relative ratio of PSD for a particular frequency band, the PSD of that band was first computed and subsequently divided by the PSD encompassing the entire frequency range of interest, specifically spanning from 0.5 Hz to 45 Hz.

2.9. Gene Network Analysis

GeneMANIA (http://genemania.org accessed on 21 July 2023) [111] is a tool that is used to generate hypotheses about gene functions, visualize gene function networks and analyze gene functions and interactions by integrating proteomic and genomic data from various sources. Given a query list, GeneMANIA lists genes that share properties or functional similarities with the original query. It also shows a functional relationship network, indicating relationships between the lists and curated genomic and proteomic data [111].

2.10. Statistical Analysis

Statistical analysis was performed using GraphPad Prism 8 (GraphPad 8.0.1 Software, San Diego, CA, USA). All data were examined for normality of distribution with the D’Agostino–Pearson and Shapiro–Wilk tests, and since all datasets did not follow a Gaussian distribution, they were then analyzed with non-parametric statistics. The DASS21 scores, mRNA levels and methylation status of DD and control individuals were compared with the non-parametric Mann–Whitney test, and molecular data were correlated with the results of the DASS21 questionnaire with Spearman’s correlation coefficient. Results are presented after correcting for and excluding outliers with the ROUT method. Data are expressed as mean ± SEM. Statistical significance was considered for p < 0.05 with 95% confidence intervals (CI). Correction for multiple statistical testing was performed using Hochberg’s correction by employing a free online tool available at https://multipletesting.com (accessed on 26 June 2023) [112].

3. Results

In the present study, all participants were assessed with the DASS21 questionnaire for depression, anxiety and stress status. Each participant provided a peripheral blood sample, which was processed for the quantification of the mRNA levels of stress-associated genes, the genotyping of stress-associated SNVs, the methylation status of NR3C1 and mtDNAcn evaluation. All participants underwent the EEG procedure, in order to correlate the molecular findings with EEG readouts.

3.1. Demographic Data, DASS21 Scores and SNVs

The psychological profiles of all participants were assessed using the DASS21 scale. Specifically, the scores for depression, anxiety and stress were measured for each participant, in order to evaluate the contributing effect and potential comorbidity of such phenotypes with DD occurrence. Descriptive statistics concerning the characteristics of DD and control human samples are presented in Table 2. Demographic information, including age (years), gender (males/females, M/F) and DASS21 scores and variants of interest are presented for individuals with DD, controls and the total sample. The control and DD groups presented non-statistically significant variations in depression and anxiety scales. Each participant was genotyped for rs1360780 of FKBP5, the 5-HTTLPR polymorphism and rs25531 of the SLC6A4 gene, albeit not significant. The allelic frequencies and statistical analysis are presented in Table 2.

3.2. Altered mRNA and Methylation Levels in NR3C1 Gene in DD

In order to study the HPA axis gene expression and methylation statuses, we analyzed the mRNA levels of the NR3C1, NR3C2, GILZ and FKBP5 genes, considered key players in the neuroendocrine stress response. As presented in Figure 1A, NR3C1 mRNA levels showed increased expression in the PBMCs of DD individuals compared to control samples (Mann–Whitney test, * p = 0.035, 95% CI: −0.003388 to −0.0001636), and the difference remained significant after applying Hochberg’s correction for multiple testing. To follow up on this finding, the methylation status across the NR3C1 5′ untranslated region (5′ UTR) was then examined (Figure 1B(i)). The 1D region (Figure 1B(ii)) showed a trend (T) for higher average methylation in DD compared to control individuals. Further statistical analysis revealed a trend towards increased methylation in CpGs 6 and 7, as well as a statistically significant increase in CpG 8 (Mann–Whitney test, * p = 0.004), in DD compared to control individuals. In the 1F region (Figure 1B(iii)), the average methylation percentage in DD samples was decreased. CpG 7 appeared hypomethylated, whereas CpGs 12 and 13 revealed a statistically significant decrease (Mann–Whitney test, * p12 = 0.036 and * p13 = 0.0304) in methylation percentages in DD compared to control samples. After applying the Hochberg correction, the differences that remained statistically significant were in CpG 8 (1D region) and CpG 13 (1F region). Data are presented in Supplementary Table S3.

3.3. Correlation of Stress-Related Molecular Changes and Behavioral Characteristics

To examine whether there was an association of the altered NR3C1 methylation profiles with the psychological characteristics of both the DD and control individuals, we correlated the methylation percentages of specific CpG sites in the regions 1D and 1F with the DASS21 scores. The methylation percentage of CpG 12 located in the 1F NR3C1 region was positively correlated with each individual’s stress scores only in DD (Spearman’s r = 0.6481 and * p = 0.0469) and not in control samples. Data are presented in Supplementary Table S4 and Supplementary Table S5, respectively. The statistical analysis of the mRNA levels of the stress-associated genes NR3C1, NR3C2, GILZ and FKBP5 did not reveal significant correlations (Supplementary Table S6).

3.4. Correlation of mRNA Levels with EEG Recordings

We correlated the individual EEG PSD from different brain areas with the mRNA levels of important stress mediators, using Spearman’s correlation (Table 3). The two groups studied showed associations in different brain areas. In DD individuals, the statistically significant associations were in the lower frequencies (delta and theta), whereas, in the control group, they were in the higher frequencies (alpha and beta) (Table 3). A schematic representation of the EEG channels that were correlated with stress-associated genes is provided in Supplementary Figure S2.

3.5. NR3C1 Changes Do Not Affect Mitochondrial DNA Copy Number (mtDNAcn)

NR3C1 has been found to play its biphasic roles by affecting mitochondrial function. Specifically, low doses of glucocorticoids functioning via GR in mitochondria have neuroprotective effects, while high doses exert neurotoxic effects [40,41]. In order to evaluate the GR impact on mitochondrial numbers, we assessed the mtDNAcn, a proxy for the assessment of mitochondrial numbers, in DD patients compared to the control group. We found no significant differences between the two groups (Mann–Whitney test, p = 0.3154, Supplementary Figure S3), suggesting that the NR3C1-observed changes did not impact mitochondrial abundance.

3.6. Gene Network Analysis

The gene network analysis included 19 target genes: genes associated with reading, speech and language skills (CMIP, CNTNAP2, CYP19A1, DCDC2, DIP2A, DYX1C1, GCFC2, KIAA0319, KIAA0319L, MRLP19, PCNT, PRMT2, S100B, ROBO1) [113,114,115] along with stress-associated genes (NR3C1, NR3C2, GILZ, FKBP5, SLC6A4). The resulting gene interaction network predicted 20 more genes that were linked to the 19 target genes (Figure 2). Table 4 lists the specific interactions of the stress-related genes with other genes in the network.

4. Discussion

DD is a complex neurodevelopmental condition. As in all multifactorial phenotypes, individual genetic variations and epigenetic differences should be considered in deciphering the underlying mechanisms. The etiological heterogeneity that characterizes these conditions is due to genetic (gene × gene interactions, G × G) and/or environmental factors (G × E) interactions but the relative impact of each factor is not always clear. The understanding of the basis of phenotype variation may aid classification and facilitate targeted interventions.
Stress is an important factor for the occurrence and progression of multifactorial phenotypes, and stress pathways are often included in relevant studies. The HPA axis includes important molecular components, such as the GR and MR receptors, which mediate the neuroendocrine response to stress and have been found altered in complex neurodevelopmental phenotypes, including DD, highlighting the existence of an inadequate stress response [29,93,116,117]. In phenotypes that are stress-related, a targeted intervention may prove effective [118,119,120]. A gene network analysis of stress- and dyslexia-associated genes revealed that the majority of the stress-associated genes studied showed some form of interaction with at least one of the known dyslexia-associated genes, supporting the notion that a number of DD susceptibility genes may respond to stress and the existence of a putative DD–stress phenotype [17,29,93] (Table 4).
This study focused on genetic and epigenetic parameters coupled with psychological assessments and brain activity. We investigated the expression levels of important HPA axis stress mediators (NR3C1, NR3C2, GILZ and FKBP5) and the methylation status of the NR3C1 promoter, coupled with DASS21 scores and EEG recordings in the two groups (DD–controls). The increased NR3C1 mRNA levels in the PBMCs of DD individuals compared to control samples (Mann–Whitney test, * p = 0.035) reinforces the epigenetically mediated effects of stress on the DD phenotype.
The DD and control groups showed no statistically significant differences in DASS21 scores, but alterations were observed in NR3C1 mRNA expression levels. This result was possibly due to the nature of the sample. The fact that the DD individuals in this study were students in higher-education institutions suggests that the longitudinal support received during their childhood was rather effective. Such an intervention may potentially contribute to the acquisition of mechanisms and strategies adequate for academic achievement, thus strengthening their stress management and positively affecting their general psycho-emotional and social development [121]. Therefore, stress effects are not reflected in the psychological evaluation and remain undetected with the current assessment techniques [122]. Exposure to adversities during early childhood could have a lifelong impact on HPA axis functioning and therefore the molecular approach probably permits a more sensitive analysis [123,124,125]. The alterations observed in the DD group vs. control group in NR3C1 mRNA expression levels could imply an inadequate stress response, mainly reflecting the impact of the memory of stressful DD experiences, characterized by struggles and failure in the school environment, which was not apparent in the DASS21.
The epigenetic analysis of the NR3C1 promoter region focused on the study of eight CpG sites within the 1D region and thirteen CpG sites within the 1F region and revealed alterations in methylation in CpG site 8 of 1D and CpG 13 of the 1F region. CpG site 8 of 1D showed higher, statistically significant, methylation levels, and CpG 13 of the 1F region showed lower, statistically significant, methylation levels in DD compared to the control group. Epigenetic modifications of the promoter region of the NR3C1 gene, including alternative non-translated first exons 1D and 1F, have been associated with GR expression in cases of stressful experiences and complex neuropsychiatric conditions [126,127,128,129,130]. Specifically, alterations in the methylation levels of CpGs across 1D and 1F are reported in the literature in cases of stressful life experiences, often related to different stressors and neuropsychiatric conditions [126,129,130,131,132,133,134,135]. The reported findings are not always reproducible because of the different conditions and pathologies studied and the variety of stressors, tissue samples and detection methods [127,131,136,137,138]. Although the reported results are sometimes contradictory, a number of studies are in agreement with this study’s results [128,132,139,140,141].
In complex heterogeneous conditions, such as DD, the variable phenotypes might be due to the incomplete penetrance and variable expressivity of the genes involved. Genomic regulation, expression and phenotypic variation are the result of the interplay of genetic and epigenetic interactions. Some of the key players in these procedures are inter- and intra-individual genomic variations, the genomic architecture, G × G interactions, the binding sites of different transcription factors and the recognition sites of epigenetic marks [142,143]. A variety of binding sites for transcription factors exist across the NR3C1 first exons [31,128,144,145,146]. Genetic variants located in binding sites may affect the methylation of different CpG sites and gene expression and contribute to phenotypic variation [147,148,149]. Moreover, the correlation between gene expression and CpG methylation can be positive or negative and depends on the CpG site location, even if it concerns neighboring CpG sites [150].
No differences concerning mtDNAcn were observed in DD individuals compared to controls in this study sample. This was probably due to the sample consisting of young adults who were able to manage the DD difficulties. According to unpublished data from our research group, different results were obtained in a study group with a different demographic composition.
Alterations in the EEG-recorded rhythms have been extensively described both in DD and stress [151,152]. In this context, unique DD patterns were identified in a sample of young adults, by analyzing EEG recordings [78]. The pairing between EEG energy and the mRNA levels of stress-related genes indicated a shift in the brain activity towards low-frequency rhythms in the sensory and association brain areas (occipital, parietal and temporal) in both hemispheres, correlated to the mRNA levels of stress-related genes (NR3C1, GILZ, FKBP5). This energy increase of lower frequencies and energy decrease of higher frequencies was observed both in the DD and control groups and aligned with previous EEG findings related to stress [153]. However, a systematic review of EEG alterations in DD revealed that this pattern was also typical for dyslexia [154], despite the broad spectrum of alterations described during specific reading tasks [85,155,156]. The observed overlap in EEG findings between different disorders underscores the common features and suggests that caution is needed in interpreting the study results [153]. The prevalence of slow waves in sensory and association brain areas could be attributed to cognitive processing driven by subcortical regions responsible for emotions and motivation [153,157]. Our results also indicate that the energy decrease of fast waves was more pronounced in the DD group than the control group. This reduced EEG energy for fast waves may be correlated with low attentional control and reduced mental focus in receiving and processing sensory information, as previously described [158], which could also be related to language processing difficulties. It is important, however, to note that significant correlations between EEG and mRNA levels for the target genes were not observed in the same brain area and frequency band combination between the control and DD groups. This discrepancy can be attributed to the differential network interactions that have been described in adults with DD [159]. Further EEG studies on brain connectivity in larger cohorts of well-characterized DD samples might illustrate a putative underlying network disruption, characteristic of a stress-related DD phenotype.
The DASS21 and EEG findings revealed a positive correlation of the methylation percentage in CpG 12 of the 1F region with the individual stress scores in the DD group. Moreover, changes in the expression patterns of the stress-related genes were also paired with a shift in brain activity towards low-frequency EEG rhythms in the affected sensory brain areas, which was possibly indicative of NR3C1-mediated HPA axis dysfunction in DD patients, which may act as a predisposing factor. The imbalance of the MR and GR populations’ distribution in typical stress centers, where they are both co-localized, could imply a putative stress-related DD phenotype.
Stress, and especially stress due to the constant difficulties and frustration associated with the school environment, may be an important factor in DD that can compromise negative feedback mechanisms, thereby leading to the persistent activation of the HPA axis [16,17]. In this study, the hypothesis that there is a stress-related DD phenotype was tested using a multidisciplinary approach composed of different modalities. In this context, the clinical diagnostic criteria of established methods, such as EEG and the DASS21 psychological assessment scale, were co-evaluated with genetic and epigenetic data, backed by bioinformatic gene network analysis data, in order to contribute to the literature an innovative model that may pave the way for the characterization of stress-specific molecular biomarkers, important for an individualized approach in DD diagnosis and carefully designed interventions. The putatively reversible nature of epigenetic alterations makes them a good candidate for stress-specific targeted interventions that aim to address environmental triggers. In the context of the hypothesis that there are DD phenotypes highly influenced by stress, this pilot study provides data and a demonstration of a new multidisciplinary approach that combines established clinical methods with EEG recordings and molecular findings of stress-associated genes.

5. Limitations

This was a pilot study with a rather small number of participants. Among the other limitations of the study, we highlight the heterogeneity of the sample population regarding the male:female ratio in the DD and control groups, and the fact that DD participants underwent an individualized intervention during their primary school years that targeted only the reading and writing difficulties and not stress management. Meanwhile, the reading, speech and language skills were not evaluated during the procedure; their academic performance was also not considered. It must be noted, however, that all participants were university students at the same institution, and they were allowed to take an oral or a written and oral exam, which may have alleviated some of the stress and anxiety. It must also be emphasized that this study was a pilot one, with the primary goal to present a new multidisciplinary transdiagnostic model and approach that includes stress in the occurrence of DD phenotypes, which hopefully will foster scientific interactions. Studies in substantially larger sample sizes, with a balanced sex distribution and carefully defined phenotypes (e.g., type of dyslexia and intervention), taking into account different types of educational systems and exams, as well as cohorts of people who have received not only speech and language interventions but also psychotherapy for stress, are required. Moreover, standardized genome-wide studies will be valuable in the understanding of the underlying mechanisms and the role and impact of stress in DD.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/brainsci14020139/s1, Figure S1: Regions of interest according to the electrode sites; Figure S2: Schematic presentation of EEG channels. Channels where EEG PSD is associated with mRNA levels of stress mediators are marked (orange—NR3C1, dark blue—GILZ, light blue—FKBP5). A) EEG channels that showed correlation in the DD group, B) EEG channels that showed correlation in the control group; Figure S3: mtDNAcn for the DD individuals compared to control (Mann–Whitney, p-value = 0.3154); Table S1: RT-qPCR primers for the genes identified; Table S2: Pyrosequencing primers for the CpG 303 in NR3C1 promoter region; Table S3: Statistical analysis of individual CG sites in 1D and 1F regions; Table S4: Spearman’s correlations for the DASS21 scale and methylation percentages for 1F region for DD individuals; Table S5: Spearman’s correlations for the DASS21 scale and methylation percentages for 1F region for control individuals; Table S6: Spearman’s correlations for the mRNA levels of stress-associated genes and DASS21 scores for DD and control individuals.

Author Contributions

Conceptualization, M.S.; methodology, D.T., V.Z., A.S., M.D.F. and M.S.; formal analysis, D.T., C.-O.T., A.-M.V., V.C., P.C., E.D.O., K.D.T., C.K. and A.T.T.; investigation, D.T., C.-O.T., A.-M.V., V.C., V.Z., P.C., K.D.T., C.K., A.T.T., K.I.T. and D.P.; resources, V.Z., K.I.T. and D.P.; data curation, V.Z., P.C., A.T.T., K.I.T., D.P., A.S., M.D.F. and M.S.; writing—original draft preparation, D.T., C.-O.T., A.-M.V. and V.C.; writing—review and editing, D.T., V.Z., P.C., A.T.T., K.I.T., A.S., M.D.F. and M.S.; visualization, D.T., C.-O.T., A.-M.V., P.C., E.D.O., K.D.T. and C.K.; supervision, V.Z., M.D.F. and M.S.; project administration, M.S.; funding acquisition, M.S. All authors have read and agreed to the published version of the manuscript.

Funding

We acknowledge the support of this work by the project “Establishment of ‘capacity building’ infrastructures in Biomedical Research (BIOMED-20)” (MIS 5047236), which is implemented under the Action “Reinforcement of the Research and Innovation Infrastructure”, funded by the Operational Program “Competitiveness, Entrepreneurship and Innovation” (NSRF 2014-2020) and co-financed by Greece and the European Union (European Regional Development Fund).

Institutional Review Board Statement

The study was conducted in accordance with the Declaration of Helsinki and approved by the Research Ethics and Conduct committee of the University of Ioannina (protocol code 28622, 23 July 2019).

Informed Consent Statement

Informed consent was obtained from all participants involved in the study.

Data Availability Statement

Raw data are available from the corresponding author upon request due to privacy/ethical restrictions.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Snowling, M.J.; Hulme, C.; Nation, K. Defining and understanding dyslexia: Past, present and future. Oxf. Rev. Educ. 2020, 46, 501–513. [Google Scholar] [CrossRef] [PubMed]
  2. Elliott, J.G.; Grigorenko, E.L. The Dyslexia Debate; Cambridge University Press: Cambridge, UK, 2014. [Google Scholar]
  3. D’Mello, A.M.; Gabrieli, J.D.E. Cognitive Neuroscience of Dyslexia. Lang. Speech Hear. Serv. Sch. 2018, 49, 798–809. [Google Scholar] [CrossRef] [PubMed]
  4. Wolf, M.; Bowers, P.G. Naming-speed processes and developmental reading disabilities: An introduction to the special issue on the double-deficit hypothesis. J. Learn. Disabil. 2000, 33, 322–324. [Google Scholar] [CrossRef]
  5. Paulesu, E.; Demonet, J.F.; Fazio, F.; McCrory, E.; Chanoine, V.; Brunswick, N.; Cappa, S.F.; Cossu, G.; Habib, M.; Frith, C.D.; et al. Dyslexia: Cultural diversity and biological unity. Science 2001, 291, 2165–2167. [Google Scholar] [CrossRef] [PubMed]
  6. Gathercole, S.E.; Baddeley, A.D. Working Memory and Language; Psychology Press: London, UK, 2014. [Google Scholar]
  7. Goswami, U. A temporal sampling framework for developmental dyslexia. Trends Cogn. Sci. 2011, 15, 3–10. [Google Scholar] [CrossRef] [PubMed]
  8. Nkomo, D.; Dube, B. The magnocellular theory of dyslexia: An overview, sub-categories and application in teaching. Eur. J. Spec. Educ. Res. 2019, 5. [Google Scholar] [CrossRef]
  9. Dubois, M.; Kyllingsbaek, S.; Prado, C.; Musca, S.C.; Peiffer, E.; Lassus-Sangosse, D.; Valdois, S. Fractionating the multi-character processing deficit in developmental dyslexia: Evidence from two case studies. Cortex 2010, 46, 717–738. [Google Scholar] [CrossRef]
  10. Snowling, M.J. Specific disorders and broader phenotypes: The case of dyslexia. Q. J. Exp. Psychol. 2008, 61, 142–156. [Google Scholar] [CrossRef]
  11. Livingston, E.M.; Siegel, L.S.; Ribary, U. Developmental dyslexia: Emotional impact and consequences. Aust. J. Learn. Difficulties 2018, 23, 107–135. [Google Scholar] [CrossRef]
  12. Menghini, D.; Finzi, A.; Benassi, M.; Bolzani, R.; Facoetti, A.; Giovagnoli, S.; Ruffino, M.; Vicari, S. Different underlying neurocognitive deficits in developmental dyslexia: A comparative study. Neuropsychologia 2010, 48, 863–872. [Google Scholar] [CrossRef]
  13. Peterson, R.L.; Pennington, B.F. Developmental dyslexia. Annu. Rev. Clin. Psychol. 2015, 11, 283–307. [Google Scholar] [CrossRef] [PubMed]
  14. Ramus, F.; Marshall, C.R.; Rosen, S.; van der Lely, H.K. Phonological deficits in specific language impairment and developmental dyslexia: Towards a multidimensional model. Brain 2013, 136, 630–645. [Google Scholar] [CrossRef] [PubMed]
  15. Cavalli, E.; Cole, P.; Leloup, G.; Poracchia-George, F.; Sprenger-Charolles, L.; El Ahmadi, A. Screening for Dyslexia in French-Speaking University Students: An Evaluation of the Detection Accuracy of the Alouette Test. J. Learn. Disabil. 2018, 51, 268–282. [Google Scholar] [CrossRef] [PubMed]
  16. Kershner, J.R. Dyslexia as an adaptation to cortico-limbic stress system reactivity. Neurobiol. Stress. 2020, 12, 100223. [Google Scholar] [CrossRef]
  17. Kershner, J.R. An Evolutionary Perspective of Dyslexia, Stress, and Brain Network Homeostasis. Front. Hum. Neurosci. 2020, 14, 575546. [Google Scholar] [CrossRef]
  18. Erbeli, F.; Rice, M.; Paracchini, S. Insights into Dyslexia Genetics Research from the Last Two Decades. Brain Sci. 2021, 12, 27. [Google Scholar] [CrossRef]
  19. Berretz, G.; Wolf, O.T.; Gunturkun, O.; Ocklenburg, S. Atypical lateralization in neurodevelopmental and psychiatric disorders: What is the role of stress? Cortex 2020, 125, 215–232. [Google Scholar] [CrossRef]
  20. Georgitsi, M.; Dermitzakis, I.; Soumelidou, E.; Bonti, E. The Polygenic Nature and Complex Genetic Architecture of Specific Learning Disorder. Brain Sci. 2021, 11, 631. [Google Scholar] [CrossRef]
  21. Skeide, M.A. Neurogenetic Insights into the Origins of Dyslexia and Dyscalculia. In The Cambridge Handbook of Dyslexia and Dyscalculia; Skeide, M.A., Ed.; Cambridge Handbooks in Psychology; Cambridge University Press: Cambridge, UK, 2022; pp. 155–161. [Google Scholar]
  22. Gialluisi, A.; Andlauer, T.F.M.; Mirza-Schreiber, N.; Moll, K.; Becker, J.; Hoffmann, P.; Ludwig, K.U.; Czamara, D.; Pourcain, B.S.; Honbolygo, F.; et al. Genome-wide association study reveals new insights into the heritability and genetic correlates of developmental dyslexia. Mol. Psychiatry 2021, 26, 3004–3017. [Google Scholar] [CrossRef]
  23. Hoeft, F.; Bouhali, F. Pre- and Postnatal Environmental Effects on Learning to Read and Mathematical Learning. In The Cambridge Handbook of Dyslexia and Dyscalculia; Skeide, M.A., Ed.; Cambridge Handbooks in Psychology; Cambridge University Press: Cambridge, UK, 2022; pp. 115–150. [Google Scholar]
  24. McEwen, B.S.; Gianaros, P.J. Stress- and allostasis-induced brain plasticity. Annu. Rev. Med. 2011, 62, 431–445. [Google Scholar] [CrossRef]
  25. Radley, J.; Morilak, D.; Viau, V.; Campeau, S. Chronic stress and brain plasticity: Mechanisms underlying adaptive and maladaptive changes and implications for stress-related CNS disorders. Neurosci. Biobehav. Rev. 2015, 58, 79–91. [Google Scholar] [CrossRef]
  26. Papadopoulou, Z.; Vlaikou, A.M.; Theodoridou, D.; Markopoulos, G.S.; Tsoni, K.; Agakidou, E.; Drosou-Agakidou, V.; Turck, C.W.; Filiou, M.D.; Syrrou, M. Stressful Newborn Memories: Pre-Conceptual, In Utero, and Postnatal Events. Front. Psychiatry 2019, 10, 220. [Google Scholar] [CrossRef]
  27. Mascheretti, S.; Bureau, A.; Battaglia, M.; Simone, D.; Quadrelli, E.; Croteau, J.; Cellino, M.R.; Giorda, R.; Beri, S.; Maziade, M.; et al. An assessment of gene-by-environment interactions in developmental dyslexia-related phenotypes. Genes. Brain Behav. 2013, 12, 47–55. [Google Scholar] [CrossRef]
  28. Sheerin, C.M.; Lind, M.J.; Bountress, K.E.; Marraccini, M.E.; Amstadter, A.B.; Bacanu, S.A.; Nugent, N.R. Meta-Analysis of Associations Between Hypothalamic-Pituitary-Adrenal Axis Genes and Risk of Posttraumatic Stress Disorder. J. Trauma. Stress. 2020, 33, 688–698. [Google Scholar] [CrossRef]
  29. Zakopoulou, V.; Vlaikou, A.M.; Darsinou, M.; Papadopoulou, Z.; Theodoridou, D.; Papageorgiou, K.; Alexiou, G.A.; Bougias, H.; Siafaka, V.; Zoccolotti, P.; et al. Linking Early Life Hypothalamic-Pituitary-Adrenal Axis Functioning, Brain Asymmetries, and Personality Traits in Dyslexia: An Informative Case Study. Front. Hum. Neurosci. 2019, 13, 327. [Google Scholar] [CrossRef]
  30. Bortoluzzi, A.; Blaya, C.; Rosa, E.D.d.; Paim, M.; Rosa, V.; Leistner-Segal, S.; Manfro, G.G. What can HPA axis-linked genes tell us about anxiety disorders in adolescents? Trends Psychiatry Psychother. 2015, 37, 232–237. [Google Scholar] [CrossRef] [PubMed]
  31. Oberlander, T.F.; Weinberg, J.; Papsdorf, M.; Grunau, R.; Misri, S.; Devlin, A.M. Prenatal exposure to maternal depression, neonatal methylation of human glucocorticoid receptor gene (NR3C1) and infant cortisol stress responses. Epigenetics 2008, 3, 97–106. [Google Scholar] [CrossRef] [PubMed]
  32. Smart, C.; Strathdee, G.; Watson, S.; Murgatroyd, C.; McAllister-Williams, R.H. Early life trauma, depression and the glucocorticoid receptor gene--an epigenetic perspective. Psychol. Med. 2015, 45, 3393–3410. [Google Scholar] [CrossRef] [PubMed]
  33. Park, C.; Rosenblat, J.D.; Brietzke, E.; Pan, Z.; Lee, Y.; Cao, B.; Zuckerman, H.; Kalantarova, A.; McIntyre, R.S. Stress, epigenetics and depression: A systematic review. Neurosci. Biobehav. Rev. 2019, 102, 139–152. [Google Scholar] [CrossRef]
  34. Thiagarajah, A.S.; Eades, L.E.; Thomas, P.R.; Guymer, E.K.; Morand, E.F.; Clarke, D.M.; Leech, M. GILZ: Glitzing up our understanding of the glucocorticoid receptor in psychopathology. Brain Res. 2014, 1574, 60–69. [Google Scholar] [CrossRef] [PubMed]
  35. Zannas, A.S.; Wiechmann, T.; Gassen, N.C.; Binder, E.B. Gene-Stress-Epigenetic Regulation of FKBP5: Clinical and Translational Implications. Neuropsychopharmacology 2016, 41, 261–274. [Google Scholar] [CrossRef]
  36. Matosin, N.; Halldorsdottir, T.; Binder, E.B. Understanding the Molecular Mechanisms Underpinning Gene by Environment Interactions in Psychiatric Disorders: The FKBP5 Model. Biol. Psychiatry 2018, 83, 821–830. [Google Scholar] [CrossRef]
  37. Zannas, A.S.; Binder, E.B. Gene-environment interactions at the FKBP5 locus: Sensitive periods, mechanisms and pleiotropism. Genes. Brain Behav. 2014, 13, 25–37. [Google Scholar] [CrossRef] [PubMed]
  38. Grabe, H.J.; Schwahn, C.; Appel, K.; Mahler, J.; Schulz, A.; Spitzer, C.; Barnow, S.; John, U.; Freyberger, H.J.; Rosskopf, D.; et al. Update on the 2005 paper: Moderation of mental and physical distress by polymorphisms in the 5-HT transporter gene by interacting with social stressors and chronic disease burden. Mol. Psychiatry 2011, 16, 354–356. [Google Scholar] [CrossRef]
  39. Palma-Gudiel, H.; Fananas, L. An integrative review of methylation at the serotonin transporter gene and its dialogue with environmental risk factors, psychopathology and 5-HTTLPR. Neurosci. Biobehav. Rev. 2017, 72, 190–209. [Google Scholar] [CrossRef]
  40. Du, J.; Wang, Y.; Hunter, R.; Wei, Y.; Blumenthal, R.; Falke, C.; Khairova, R.; Zhou, R.; Yuan, P.; Machado-Vieira, R.; et al. Dynamic regulation of mitochondrial function by glucocorticoids. Proc. Natl. Acad. Sci. USA 2009, 106, 3543–3548. [Google Scholar] [CrossRef] [PubMed]
  41. Hunter, R.G.; Seligsohn, M.; Rubin, T.G.; Griffiths, B.B.; Ozdemir, Y.; Pfaff, D.W.; Datson, N.A.; McEwen, B.S. Stress and corticosteroids regulate rat hippocampal mitochondrial DNA gene expression via the glucocorticoid receptor. Proc. Natl. Acad. Sci. USA 2016, 113, 9099–9104. [Google Scholar] [CrossRef]
  42. Picard, M.; McEwen, B.S. Psychological Stress and Mitochondria: A Systematic Review. Psychosom. Med. 2018, 80, 141–153. [Google Scholar] [CrossRef]
  43. Lopes, S.; Teplytska, L.; Vaz-Silva, J.; Dioli, C.; Trindade, R.; Morais, M.; Webhofer, C.; Maccarrone, G.; Almeida, O.F.X.; Turck, C.W.; et al. Tau Deletion Prevents Stress-Induced Dendritic Atrophy in Prefrontal Cortex: Role of Synaptic Mitochondria. Cereb. Cortex 2017, 27, 2580–2591. [Google Scholar] [CrossRef] [PubMed]
  44. Filiou, M.D.; Sandi, C. Anxiety and Brain Mitochondria: A Bidirectional Crosstalk. Trends Neurosci. 2019, 42, 573–588. [Google Scholar] [CrossRef]
  45. Nussbaumer, M.; Asara, J.M.; Teplytska, L.; Murphy, M.P.; Logan, A.; Turck, C.W.; Filiou, M.D. Selective Mitochondrial Targeting Exerts Anxiolytic Effects In Vivo. Neuropsychopharmacology 2016, 41, 1751–1758. [Google Scholar] [CrossRef] [PubMed]
  46. Filiou, M.D.; Zhang, Y.; Teplytska, L.; Reckow, S.; Gormanns, P.; Maccarrone, G.; Frank, E.; Kessler, M.S.; Hambsch, B.; Nussbaumer, M.; et al. Proteomics and metabolomics analysis of a trait anxiety mouse model reveals divergent mitochondrial pathways. Biol. Psychiatry 2011, 70, 1074–1082. [Google Scholar] [CrossRef] [PubMed]
  47. Christodoulides, P.; Zakopoulou, V. The contribution of EEG rhythms’ changes to the audiovisual recognition of words in university students with dyslexia. In Psychology Applications & Developments VII Advances in Psychology and Psychological Trends Series; Pracana, C., Wang, M., Eds.; inScience Press: Lisbon, Portugal, 2021; p. 181. [Google Scholar]
  48. Filiou, M.D.; Asara, J.M.; Nussbaumer, M.; Teplytska, L.; Landgraf, R.; Turck, C.W. Behavioral extremes of trait anxiety in mice are characterized by distinct metabolic profiles. J. Psychiatr. Res. 2014, 58, 115–122. [Google Scholar] [CrossRef] [PubMed]
  49. Petretto, D.R.; Masala, C. Dyslexia and specific learning disorders: New international diagnostic criteria. J. Child. Dev. Disord. 2017, 3, 4–19. [Google Scholar] [CrossRef]
  50. Hamilton, S.S.; Glascoe, F.P. Evaluation of children with reading difficulties. Am. Fam. Physician 2006, 74, 2079–2084. [Google Scholar]
  51. Francis, D.A.; Caruana, N.; Hudson, J.L.; McArthur, G.M. The association between poor reading and internalising problems: A systematic review and meta-analysis. Clin. Psychol. Rev. 2019, 67, 45–60. [Google Scholar] [CrossRef]
  52. Giovagnoli, S.; Mandolesi, L.; Magri, S.; Gualtieri, L.; Fabbri, D.; Tossani, E.; Benassi, M. Internalizing Symptoms in Developmental Dyslexia: A Comparison Between Primary and Secondary School. Front. Psychol. 2020, 11, 461. [Google Scholar] [CrossRef]
  53. Schoenberg, P.L. Linear and nonlinear EEG-based functional networks in anxiety disorders. Anxiety Disord. 2020, 1191, 35–59. [Google Scholar]
  54. Vajs, I.; Kovic, V.; Papic, T.; Savic, A.M.; Jankovic, M.M. Spatiotemporal Eye-Tracking Feature Set for Improved Recognition of Dyslexic Reading Patterns in Children. Sensors 2022, 22, 4900. [Google Scholar] [CrossRef]
  55. Jakovljevic, T.; Jankovic, M.M.; Savic, A.M.; Soldatovic, I.; Todorovic, P.; Jere Jakulin, T.; Papa, G.; Kovic, V. The Sensor Hub for Detecting the Developmental Characteristics in Reading in Children on a White vs. Colored Background/Colored Overlays. Sensors 2021, 21, 406. [Google Scholar] [CrossRef]
  56. Janković, M.M. Biomarker-based approaches for dyslexia screening: A review. In Proceedings of the 2022 IEEE Zooming Innovation in Consumer Technologies Conference (ZINC), Novi Sad, Serbia, 25–26 May 2022; pp. 28–33. [Google Scholar]
  57. Christoforou, C.; Fella, A.; Leppanen, P.H.T.; Georgiou, G.K.; Papadopoulos, T.C. Fixation-related potentials in naming speed: A combined EEG and eye-tracking study on children with dyslexia. Clin. Neurophysiol. 2021, 132, 2798–2807. [Google Scholar] [CrossRef]
  58. Evans, K.C.; Dougherty, D.D.; Pollack, M.H.; Rauch, S.L. Using neuroimaging to predict treatment response in mood and anxiety disorders. Ann. Clin. Psychiatry 2006, 18, 33–42. [Google Scholar] [CrossRef]
  59. Garcia Chimeno, Y.; Garcia Zapirain, B.; Saralegui Prieto, I.; Fernandez-Ruanova, B. Automatic classification of dyslexic children by applying machine learning to fMRI images. Biomed. Mater. Eng. 2014, 24, 2995–3002. [Google Scholar] [CrossRef]
  60. Ortiz, A.; Martinez-Murcia, F.J.; Luque, J.L.; Gimenez, A.; Morales-Ortega, R.; Ortega, J. Dyslexia Diagnosis by EEG Temporal and Spectral Descriptors: An Anomaly Detection Approach. Int. J. Neural Syst. 2020, 30, 2050029. [Google Scholar] [CrossRef]
  61. Thiede, A.; Glerean, E.; Kujala, T.; Parkkonen, L. Atypical MEG inter-subject correlation during listening to continuous natural speech in dyslexia. Neuroimage 2020, 216, 116799. [Google Scholar] [CrossRef] [PubMed]
  62. Xing, S.; McCardle, R.; Xie, S.S.Q. Reading the mind: The potential of electroencephalography in brain computer interfaces. In Proceedings of the 2012 19th International Conference on Mechatronics and Machine Vision in Practice (M2VIP), Auckland, New Zealand, 28–30 November 2012; pp. 275–280. [Google Scholar]
  63. Martínez-Murcia, F.; Ortiz, A.; Morales, R.; López, P.; Luque, J.; Castillo-Barnes, D.; Segovia, F.; Illan, I.; Ortega, J.; Ramírez, J.; et al. Periodogram Connectivity of EEG Signals for the Detection of Dyslexia. In Understanding the Brain Function and Emotions: 8th International Work-Conference on the Interplay Between Natural and Artificial Computation, IWINAC 2019, Almería, Spain, June 3–7, 2019, Proceedings; Springer International Publishing: Berlin/Heidelberg, Germany, 2019; pp. 350–359. [Google Scholar]
  64. Xue, H.; Wang, Z.; Tan, Y.; Yang, H.; Fu, W.; Xue, L.; Zhao, J. Resting-state EEG reveals global network deficiency in dyslexic children. Neuropsychologia 2020, 138, 107343. [Google Scholar] [CrossRef] [PubMed]
  65. Axmacher, N.; Mormann, F.; Fernandez, G.; Elger, C.E.; Fell, J. Memory formation by neuronal synchronization. Brain Res. Rev. 2006, 52, 170–182. [Google Scholar] [CrossRef] [PubMed]
  66. McCormick, D.A.; Contreras, D. On the cellular and network bases of epileptic seizures. Annu. Rev. Physiol. 2001, 63, 815–846. [Google Scholar] [CrossRef] [PubMed]
  67. Rippon, G.; Brunswick, N. Trait and state EEG indices of information processing in developmental dyslexia. Int. J. Psychophysiol. 2000, 36, 251–265. [Google Scholar] [CrossRef] [PubMed]
  68. Papagiannopoulou, E.A.; Lagopoulos, J. Resting State EEG Hemispheric Power Asymmetry in Children with Dyslexia. Front. Pediatr. 2016, 4, 11. [Google Scholar] [CrossRef]
  69. Kandel, S.; Lassus-Sangosse, D.; Grosjacques, G.; Perret, C. The impact of developmental dyslexia and dysgraphia on movement production during word writing. Cogn. Neuropsychol. 2017, 34, 219–251. [Google Scholar] [CrossRef]
  70. Perera, H.; Shiratuddin, M.F.; Wong, K.W. Review of EEG-based pattern classification frameworks for dyslexia. Brain Inform. 2018, 5, 4. [Google Scholar] [CrossRef] [PubMed]
  71. Peyrin, C.; Lallier, M.; Demonet, J.F.; Pernet, C.; Baciu, M.; Le Bas, J.F.; Valdois, S. Neural dissociation of phonological and visual attention span disorders in developmental dyslexia: FMRI evidence from two case reports. Brain Lang. 2012, 120, 381–394. [Google Scholar] [CrossRef] [PubMed]
  72. Steinbrink, C.; Groth, K.; Lachmann, T.; Riecker, A. Neural correlates of temporal auditory processing in developmental dyslexia during German vowel length discrimination: An fMRI study. Brain Lang. 2012, 121, 1–11. [Google Scholar] [CrossRef]
  73. Pecini, C.; Biagi, L.; Brizzolara, D.; Cipriani, P.; Di Lieto, M.C.; Guzzetta, A.; Tosetti, M.; Chilosi, A.M. How many functional brains in developmental dyslexia? When the history of language delay makes the difference. Cogn. Behav. Neurol. 2011, 24, 85–92. [Google Scholar] [CrossRef]
  74. Dufor, O.; Serniclaes, W.; Sprenger-Charolles, L.; Demonet, J.F. Left premotor cortex and allophonic speech perception in dyslexia: A PET study. Neuroimage 2009, 46, 241–248. [Google Scholar] [CrossRef] [PubMed]
  75. Nuwer, M.R. EEG topographic mapping and frequency analysis: Techniques and studies in clinical setting. Semin. Neurol. 1990, 10, 166–177. [Google Scholar] [CrossRef]
  76. Kershner, J.R. Neuroscience and education: Cerebral lateralization of networks and oscillations in dyslexia. Laterality 2020, 25, 109–125. [Google Scholar] [CrossRef]
  77. Di Liberto, G.M.; Peter, V.; Kalashnikova, M.; Goswami, U.; Burnham, D.; Lalor, E.C. Atypical cortical entrainment to speech in the right hemisphere underpins phonemic deficits in dyslexia. Neuroimage 2018, 175, 70–79. [Google Scholar] [CrossRef]
  78. Christodoulides, P.; Miltiadous, A.; Tzimourta, K.D.; Peschos, D.; Ntritsos, G.; Zakopoulou, V.; Giannakeas, N.; Astrakas, L.G.; Tsipouras, M.G.; Tsamis, K.I. Classification of EEG signals from young adults with dyslexia combining a Brain Computer Interface device and an Interactive Linguistic Software Tool. Biomed. Signal Process. Control. 2022, 76, 103646. [Google Scholar] [CrossRef]
  79. Eroglu, G.; Köprülü, M.; Karabacak, B. Developmental dyslexia biomarker detection with Quantitative electroencephalography (QEEG) data in children: Feasibility, acceptability, economic impact. Qeios 2022. [Google Scholar] [CrossRef]
  80. Kaisar, S. Developmental dyslexia detection using machine learning techniques: A survey. ICT Express 2020, 6, 181–184. [Google Scholar] [CrossRef]
  81. Ahire, N.; Awale, R.; Patnaik, S.; Wagh, A. A comprehensive review of machine learning approaches for dyslexia diagnosis. Multimed. Tools Appl. 2022, 82, 13557–13577. [Google Scholar] [CrossRef]
  82. Rezvani, Z.; Zare, M.; Žarić, G.; Bonte, M.; Tijms, J.; Molen, M.W.V.d.; González, G.F. Machine learning Classification of Dyslexic Children based on EEG Local Network Features. bioRxiv 2019. [Google Scholar] [CrossRef]
  83. Karim, I.; Qayoom, A.; Wahab, A.; Kamaruddin, N. Early Identification of Dyslexic Preschoolers Based on Neurophysiological Signals. In Proceedings of the 2013 International Conference on Advanced Computer Science Applications and Technologies, Kuching, Malaysia, 23–24 December 2013; pp. 362–366. [Google Scholar]
  84. Zainuddin, A.Z.A.; Lee, K.Y.; Mansor, W.; Mahmoodin, Z. Optimized KNN classify rule for EEG based differentiation between capable dyslexic and normal children. In Proceedings of the 2016 IEEE EMBS Conference on Biomedical Engineering and Sciences (IECBES), Kuala Lumpur, Malaysia, 4–8 December 2016; pp. 685–688. [Google Scholar]
  85. Jakovljevic, T.; Jankovic, M.M.; Savic, A.M.; Soldatovic, I.; Colic, G.; Jakulin, T.J.; Papa, G.; Kovic, V. The Relation between Physiological Parameters and Colour Modifications in Text Background and Overlay during Reading in Children with and without Dyslexia. Brain Sci. 2021, 11, 539. [Google Scholar] [CrossRef] [PubMed]
  86. Mammarella, I.C.; Cardillo, R.; Zoccante, L. Differences in visuospatial processing in individuals with nonverbal learning disability or autism spectrum disorder without intellectual disability. Neuropsychology 2019, 33, 123–134. [Google Scholar] [CrossRef] [PubMed]
  87. Mahmoodin, Z.; Mansor, W.; Lee, K.Y.; Mohamad, N.B. An analysis of EEG signal power spectrum density generated during writing in children with dyslexia. In Proceedings of the 2015 IEEE 11th International Colloquium on Signal Processing & Its Applications (CSPA), Kuala Lumpur, Malaysia, 6–8 March 2015; pp. 156–160. [Google Scholar]
  88. Penolazzi, B.; Spironelli, C.; Vio, C.; Angrilli, A. Brain plasticity in developmental dyslexia after phonological treatment: A beta EEG band study. Behav. Brain Res. 2010, 209, 179–182. [Google Scholar] [CrossRef]
  89. Slikker, W., Jr. Biomarkers and their impact on precision medicine. Exp. Biol. Med. 2018, 243, 211–212. [Google Scholar] [CrossRef] [PubMed]
  90. Abi-Dargham, A.; Moeller, S.J.; Ali, F.; DeLorenzo, C.; Domschke, K.; Horga, G.; Jutla, A.; Kotov, R.; Paulus, M.P.; Rubio, J.M.; et al. Candidate biomarkers in psychiatric disorders: State of the field. World Psychiatry 2023, 22, 236–262. [Google Scholar] [CrossRef]
  91. Kuznetsova, I.L.; Ponomareva, N.V.; Alemastseva, E.A.; Manakhov, A.D.; Andreeva, T.V.; Gusev, F.E.; Rogaev, E.I. The Interactive Effect of Genetic and Epigenetic Variations in FKBP5 and ApoE Genes on Anxiety and Brain EEG Parameters. Genes 2022, 13, 164. [Google Scholar] [CrossRef]
  92. Wirz, L.; Reuter, M.; Wacker, J.; Felten, A.; Schwabe, L. A Haplotype Associated with Enhanced Mineralocorticoid Receptor Expression Facilitates the Stress-Induced Shift from “Cognitive” to “Habit” Learning. eNeuro 2017, 4, 6. [Google Scholar] [CrossRef]
  93. Theodoridou, D.; Christodoulides, P.; Zakopoulou, V.; Syrrou, M. Developmental Dyslexia: Environment Matters. Brain Sci. 2021, 11, 7821. [Google Scholar] [CrossRef]
  94. Becker, N.; Vasconcelos, M.; Oliveira, V.; Santos, F.C.D.; Bizarro, L.; Almeida, R.M.M.; Salles, J.F.; Carvalho, M.R.S. Genetic and environmental risk factors for developmental dyslexia in children: Systematic review of the last decade. Dev. Neuropsychol. 2017, 42, 423–445. [Google Scholar] [CrossRef]
  95. Christodoulides, P.; Zakopoulou, V.; Tzimourta, K.D.; Tzallas, A.T.; Peschos, D. The contribution of EEG recordings to the audiovisual recognition of words in university students with dyslexia. Psychol. Appl. Conf. Trends 2021. [Google Scholar] [CrossRef]
  96. Lovibond, P.F.; Lovibond, S.H. The structure of negative emotional states: Comparison of the Depression Anxiety Stress Scales (DASS) with the Beck Depression and Anxiety Inventories. Behav. Res. Ther. 1995, 33, 335–343. [Google Scholar] [CrossRef] [PubMed]
  97. Lan, H.T.Q.; Long, N.T.; Hanh, N.V. Validation of Depression, Anxiety and Stress Scales (DASS-21): Immediate Psychological Responses of Students in the E-Learning Environment. Int. J. High. Educ. 2020, 9, 125–133. [Google Scholar] [CrossRef]
  98. Lyrakos, G.; Arvaniti, C.; Smyrnioti, M.; Kostopanagiotou, G. Translation and validation study of the depression anxiety stress scale in the Greek general population and in a psychiatric patient’s sample. Eur. Psychiatry 2011, 26, 1731. [Google Scholar] [CrossRef]
  99. Xie, H. Chapter 5—Personalized Epigenetics: Analysis and Interpretation of DNA Methylation Variation. In Personalized Epigenetics; Tollefsbol, T.O., Ed.; Academic Press: Boston, MA, USA, 2015; pp. 123–150. [Google Scholar]
  100. World Medical Association Declaration of Helsinki: Ethical principles for medical research involving human subjects. JAMA 2013, 310, 2191–2194. [CrossRef]
  101. Kosek, E.; Jensen, K.B.; Lonsdorf, T.B.; Schalling, M.; Ingvar, M. Genetic variation in the serotonin transporter gene (5-HTTLPR, rs25531) influences the analgesic response to the short acting opioid Remifentanil in humans. Mol. Pain. 2009, 5, 37. [Google Scholar] [CrossRef] [PubMed]
  102. Al-Eitan, L.N.; Jaradat, S.A.; Qin, W.; Wildenauer, D.M.; Wildenauer, D.D.; Hulse, G.K.; Tay, G.K. Characterization of serotonin transporter gene (SLC6A4) polymorphisms and its association with drug dependence in a Jordanian Arab population. Toxicol. Ind. Health 2014, 30, 598–610. [Google Scholar] [CrossRef]
  103. Livak, K.J.; Schmittgen, T.D. Analysis of relative gene expression data using real-time quantitative PCR and the 2(-Delta Delta C(T)) Method. Methods 2001, 25, 402–408. [Google Scholar] [CrossRef]
  104. Kent, W.J.; Sugnet, C.W.; Furey, T.S.; Roskin, K.M.; Pringle, T.H.; Zahler, A.M.; Haussler, D. The human genome browser at UCSC. Genome Res. 2002, 12, 996–1006. [Google Scholar] [CrossRef]
  105. Jiang, Z.; Qian, L.; Zou, H.; Jia, Y.; Ni, Y.; Yang, X.; Jiang, Z.; Zhao, R. Porcine glucocorticoid receptor (NR3C1) gene: Tissue-specificity of transcriptional strength and glucocorticoid responsiveness of alternative promoters. J. Steroid Biochem. Mol. Biol. 2014, 141, 87–93. [Google Scholar] [CrossRef]
  106. Chen, S.; Li, Z.; He, Y.; Zhang, F.; Li, H.; Liao, Y.; Wei, Z.; Wan, G.; Xiang, X.; Hu, M.; et al. Elevated mitochondrial DNA copy number in peripheral blood cells is associated with childhood autism. BMC Psychiatry 2015, 15, 50. [Google Scholar] [CrossRef]
  107. Turri, C.; Di Dona, G.; Santoni, A.; Zamfira, D.A.; Franchin, L.; Melcher, D.; Ronconi, L. Periodic and Aperiodic EEG Features as Potential Markers of Developmental Dyslexia. Biomedicines 2023, 11, 1607. [Google Scholar] [CrossRef]
  108. Chikhi, S.; Matton, N.; Blanchet, S. EEG power spectral measures of cognitive workload: A meta-analysis. Psychophysiology 2022, 59, e14009. [Google Scholar] [CrossRef] [PubMed]
  109. Basar, E.; Golbasi, B.T.; Tulay, E.; Aydin, S.; Basar-Eroglu, C. Best method for analysis of brain oscillations in healthy subjects and neuropsychiatric diseases. Int. J. Psychophysiol. 2016, 103, 22–42. [Google Scholar] [CrossRef] [PubMed]
  110. Katmah, R.; Al-Shargie, F.; Tariq, U.; Babiloni, F.; Al-Mughairbi, F.; Al-Nashash, H. A Review on Mental Stress Assessment Methods Using EEG Signals. Sensors 2021, 21, 5043. [Google Scholar] [CrossRef] [PubMed]
  111. Warde-Farley, D.; Donaldson, S.L.; Comes, O.; Zuberi, K.; Badrawi, R.; Chao, P.; Franz, M.; Grouios, C.; Kazi, F.; Lopes, C.T.; et al. The GeneMANIA prediction server: Biological network integration for gene prioritization and predicting gene function. Nucleic Acids Res. 2010, 38, W214–W220. [Google Scholar] [CrossRef] [PubMed]
  112. Menyhart, O.; Weltz, B.; Gyorffy, B. MultipleTesting.com: A tool for life science researchers for multiple hypothesis testing correction. PLoS ONE 2021, 16, e0245824. [Google Scholar] [CrossRef]
  113. Doust, C.; Gordon, S.D.; Garden, N.; Fisher, S.E.; Martin, N.G.; Bates, T.C.; Luciano, M. The Association of Dyslexia and Developmental Speech and Language Disorder Candidate Genes with Reading and Language Abilities in Adults. Twin Res. Hum. Genet. 2020, 23, 23–32. [Google Scholar] [CrossRef]
  114. Luciano, M.; Gow, A.J.; Pattie, A.; Bates, T.C.; Deary, I.J. The Influence of Dyslexia Candidate Genes on Reading Skill in Old Age. Behav. Genet. 2018, 48, 351–360. [Google Scholar] [CrossRef]
  115. Doust, C.; Fontanillas, P.; Eising, E.; Gordon, S.D.; Wang, Z.; Alagoz, G.; Molz, B.; andMe Research, T.; Quantitative Trait Working Group of the GenLang, C.; Pourcain, B.S.; et al. Discovery of 42 genome-wide significant loci associated with dyslexia. Nat. Genet. 2022, 54, 1621–1629. [Google Scholar] [CrossRef] [PubMed]
  116. Tsigos, C.; Chrousos, G.P. Hypothalamic-pituitary-adrenal axis, neuroendocrine factors and stress. J. Psychosom. Res. 2002, 53, 865–871. [Google Scholar] [CrossRef] [PubMed]
  117. Espin, L.; Garcia, I.; Del Pino Sanchez, M.; Roman, F.; Salvador, A. Effects of psychosocial stress on the hormonal and affective response in children with dyslexia. Trends Neurosci. Educ. 2019, 15, 1–9. [Google Scholar] [CrossRef]
  118. Kuracki, K.; Dluzniewska, A. Exam stress and the metacognitive strategies of reading in students with dyslexia: The role of motivational mechanisms and educational support. PLoS ONE 2023, 18, e0294255. [Google Scholar] [CrossRef] [PubMed]
  119. Francis, D.; Hudson, J.L.; Kohnen, S.; Mobach, L.; McArthur, G.M. The effect of an integrated reading and anxiety intervention for poor readers with anxiety. PeerJ 2021, 9, e10987. [Google Scholar] [CrossRef]
  120. Wilmot, A.; Hasking, P.; Leitao, S.; Hill, E.; Boyes, M. Understanding Mental Health in Developmental Dyslexia: A Scoping Review. Int. J. Environ. Res. Public. Health 2023, 20, 1653. [Google Scholar] [CrossRef]
  121. Zakopoulou, V.; Mavreas, V.; Christodoulides, P.; Lavidas, A.; Fili, E.; Georgiou, G.; Dimakopoulos, G.; Vergou, M. Specific learning difficulties: A retrospective study of their co morbidity and continuity as early indicators of mental disorders. Res. Dev. Disabil. 2014, 35, 3496–3507. [Google Scholar] [CrossRef]
  122. Allsopp, K.; Read, J.; Corcoran, R.; Kinderman, P. Heterogeneity in psychiatric diagnostic classification. Psychiatry Res. 2019, 279, 15–22. [Google Scholar] [CrossRef]
  123. Zankert, S.; Bellingrath, S.; Wust, S.; Kudielka, B.M. HPA axis responses to psychological challenge linking stress and disease: What do we know on sources of intra- and interindividual variability? Psychoneuroendocrinology 2019, 105, 86–97. [Google Scholar] [CrossRef] [PubMed]
  124. Hunter, A.L.; Minnis, H.; Wilson, P. Altered stress responses in children exposed to early adversity: A systematic review of salivary cortisol studies. Stress 2011, 14, 614–626. [Google Scholar] [CrossRef] [PubMed]
  125. Bunea, I.M.; Szentagotai-Tatar, A.; Miu, A.C. Early-life adversity and cortisol response to social stress: A meta-analysis. Transl. Psychiatry 2017, 7, 1274. [Google Scholar] [CrossRef] [PubMed]
  126. Roy, B.; Shelton, R.C.; Dwivedi, Y. DNA methylation and expression of stress related genes in PBMC of MDD patients with and without serious suicidal ideation. J. Psychiatr. Res. 2017, 89, 115–124. [Google Scholar] [CrossRef]
  127. Palma-Gudiel, H.; Cordova-Palomera, A.; Leza, J.C.; Fananas, L. Glucocorticoid receptor gene (NR3C1) methylation processes as mediators of early adversity in stress-related disorders causality: A critical review. Neurosci. Biobehav. Rev. 2015, 55, 520–535. [Google Scholar] [CrossRef]
  128. McGowan, P.O.; Sasaki, A.; D’Alessio, A.C.; Dymov, S.; Labonte, B.; Szyf, M.; Turecki, G.; Meaney, M.J. Epigenetic regulation of the glucocorticoid receptor in human brain associates with childhood abuse. Nat. Neurosci. 2009, 12, 342–348. [Google Scholar] [CrossRef]
  129. Wang, W.; Feng, J.; Ji, C.; Mu, X.; Ma, Q.; Fan, Y.; Chen, C.; Gao, C.; Ma, X.C.; Zhu, F. Increased methylation of glucocorticoid receptor gene promoter 1(F) in peripheral blood of patients with generalized anxiety disorder. J. Psychiatr. Res. 2017, 91, 18–25. [Google Scholar] [CrossRef]
  130. Labonte, B.; Yerko, V.; Gross, J.; Mechawar, N.; Meaney, M.J.; Szyf, M.; Turecki, G. Differential glucocorticoid receptor exon 1(B), 1(C), and 1(H) expression and methylation in suicide completers with a history of childhood abuse. Biol. Psychiatry 2012, 72, 41–48. [Google Scholar] [CrossRef]
  131. Turecki, G.; Meaney, M.J. Effects of the Social Environment and Stress on Glucocorticoid Receptor Gene Methylation: A Systematic Review. Biol. Psychiatry 2016, 79, 87–96. [Google Scholar] [CrossRef] [PubMed]
  132. Hompes, T.; Izzi, B.; Gellens, E.; Morreels, M.; Fieuws, S.; Pexsters, A.; Schops, G.; Dom, M.; Van Bree, R.; Freson, K.; et al. Investigating the influence of maternal cortisol and emotional state during pregnancy on the DNA methylation status of the glucocorticoid receptor gene (NR3C1) promoter region in cord blood. J. Psychiatr. Res. 2013, 47, 880–891. [Google Scholar] [CrossRef] [PubMed]
  133. Schur, R.R.; Boks, M.P.; Rutten, B.P.F.; Daskalakis, N.P.; de Nijs, L.; van Zuiden, M.; Kavelaars, A.; Heijnen, C.J.; Joels, M.; Kahn, R.S.; et al. Longitudinal changes in glucocorticoid receptor exon 1(F) methylation and psychopathology after military deployment. Transl. Psychiatry 2017, 7, e1181. [Google Scholar] [CrossRef]
  134. Romens, S.E.; McDonald, J.; Svaren, J.; Pollak, S.D. Associations between early life stress and gene methylation in children. Child. Dev. 2015, 86, 303–309. [Google Scholar] [CrossRef] [PubMed]
  135. van der Knaap, L.J.; Riese, H.; Hudziak, J.J.; Verbiest, M.M.; Verhulst, F.C.; Oldehinkel, A.J.; van Oort, F.V. Glucocorticoid receptor gene (NR3C1) methylation following stressful events between birth and adolescence. The TRAILS study. Transl. Psychiatry 2014, 4, e381. [Google Scholar] [CrossRef] [PubMed]
  136. Parade, S.H.; Ridout, K.K.; Seifer, R.; Armstrong, D.A.; Marsit, C.J.; McWilliams, M.A.; Tyrka, A.R. Methylation of the Glucocorticoid Receptor Gene Promoter in Preschoolers: Links With Internalizing Behavior Problems. Child. Dev. 2016, 87, 86–97. [Google Scholar] [CrossRef] [PubMed]
  137. Daskalakis, N.P.; Yehuda, R. Site-specific methylation changes in the glucocorticoid receptor exon 1F promoter in relation to life adversity: Systematic review of contributing factors. Front. Neurosci. 2014, 8, 369. [Google Scholar] [CrossRef] [PubMed]
  138. Dee, G.; Ryznar, R.; Dee, C. Epigenetic Changes Associated with Different Types of Stressors and Suicide. Cells 2023, 12, 1258. [Google Scholar] [CrossRef] [PubMed]
  139. Carmi, L.; Zohar, J.; Juven-Wetzler, A.; Desarnaud, F.; Makotkine, L.; Bierer, L.M.; Cohen, H.; Yehuda, R. Promoter methylation of the glucocorticoid receptor following trauma may be associated with subsequent development of PTSD. World J. Biol. Psychiatry 2023, 24, 578–586. [Google Scholar] [CrossRef] [PubMed]
  140. Tyrka, A.R.; Parade, S.H.; Welch, E.S.; Ridout, K.K.; Price, L.H.; Marsit, C.; Philip, N.S.; Carpenter, L.L. Methylation of the leukocyte glucocorticoid receptor gene promoter in adults: Associations with early adversity and depressive, anxiety and substance-use disorders. Transl. Psychiatry 2016, 6, e848. [Google Scholar] [CrossRef] [PubMed]
  141. Yehuda, R.; Flory, J.D.; Bierer, L.M.; Henn-Haase, C.; Lehrner, A.; Desarnaud, F.; Makotkine, I.; Daskalakis, N.P.; Marmar, C.R.; Meaney, M.J. Lower methylation of glucocorticoid receptor gene promoter 1F in peripheral blood of veterans with posttraumatic stress disorder. Biol. Psychiatry 2015, 77, 356–364. [Google Scholar] [CrossRef]
  142. Brenet, F.; Moh, M.; Funk, P.; Feierstein, E.; Viale, A.J.; Socci, N.D.; Scandura, J.M. DNA methylation of the first exon is tightly linked to transcriptional silencing. PLoS ONE 2011, 6, e14524. [Google Scholar] [CrossRef]
  143. Capp, J.P. Interplay between genetic, epigenetic, and gene expression variability: Considering complexity in evolvability. Evol. Appl. 2021, 14, 893–901. [Google Scholar] [CrossRef] [PubMed]
  144. Turner, J.D.; Muller, C.P. Structure of the glucocorticoid receptor (NR3C1) gene 5′ untranslated region: Identification, and tissue distribution of multiple new human exon 1. J. Mol. Endocrinol. 2005, 35, 283–292. [Google Scholar] [CrossRef] [PubMed]
  145. Tyrka, A.R.; Parade, S.H.; Eslinger, N.M.; Marsit, C.J.; Lesseur, C.; Armstrong, D.A.; Philip, N.S.; Josefson, B.; Seifer, R. Methylation of exons 1D, 1F, and 1H of the glucocorticoid receptor gene promoter and exposure to adversity in preschool-aged children. Dev. Psychopathol. 2015, 27, 577–585. [Google Scholar] [CrossRef] [PubMed]
  146. Giarraputo, J.; DeLoach, J.; Padbury, J.; Uzun, A.; Marsit, C.; Hawes, K.; Lester, B. Medical morbidities and DNA methylation of NR3C1 in preterm infants. Pediatr. Res. 2017, 81, 68–74. [Google Scholar] [CrossRef] [PubMed]
  147. Tseng, C.C.; Wong, M.C.; Liao, W.T.; Chen, C.J.; Lee, S.C.; Yen, J.H.; Chang, S.J. Genetic Variants in Transcription Factor Binding Sites in Humans: Triggered by Natural Selection and Triggers of Diseases. Int. J. Mol. Sci. 2021, 22, 4187. [Google Scholar] [CrossRef] [PubMed]
  148. Zhu, H.; Wang, G.; Qian, J. Transcription factors as readers and effectors of DNA methylation. Nat. Rev. Genet. 2016, 17, 551–565. [Google Scholar] [CrossRef] [PubMed]
  149. Hop, P.J.; Luijk, R.; Daxinger, L.; van Iterson, M.; Dekkers, K.F.; Jansen, R.; Consortium, B.; van Meurs, J.B.J.; t Hoen, P.A.C.; Ikram, M.A.; et al. Genome-wide identification of genes regulating DNA methylation using genetic anchors for causal inference. Genome Biol. 2020, 21, 220. [Google Scholar] [CrossRef]
  150. Spainhour, J.C.; Lim, H.S.; Yi, S.V.; Qiu, P. Correlation Patterns Between DNA Methylation and Gene Expression in The Cancer Genome Atlas. Cancer Inform. 2019, 18, 1176935119828776. [Google Scholar] [CrossRef]
  151. Sklar, B.; Hanley, J.; Simmons, W.W. An EEG experiment aimed toward identifying dyslexic children. Nature 1972, 240, 414–416. [Google Scholar] [CrossRef]
  152. Ulett, G.A.; Gleser, G. The Effect of Experimental Stress upon the Photically Activated EEG. Science 1952, 115, 678–682. [Google Scholar] [CrossRef]
  153. Newson, J.J.; Thiagarajan, T.C. EEG Frequency Bands in Psychiatric Disorders: A Review of Resting State Studies. Front. Hum. Neurosci. 2018, 12, 521. [Google Scholar] [CrossRef] [PubMed]
  154. Cainelli, E.; Vedovelli, L.; Carretti, B.; Bisiacchi, P. EEG correlates of developmental dyslexia: A systematic review. Ann. Dyslexia 2023, 73, 184–213. [Google Scholar] [CrossRef] [PubMed]
  155. Klimesch, W.; Doppelmayr, M.; Wimmer, H.; Gruber, W.; Rohm, D.; Schwaiger, J.; Hutzler, F. Alpha and beta band power changes in normal and dyslexic children. Clin. Neurophysiol. 2001, 112, 1186–1195. [Google Scholar] [CrossRef]
  156. João, A.; Benjamin, D.S.; Varghese, P.; Kanad, M.; Marina, K.; Annabel, M.; Fiona, G.; Angela, W.; Giovanni, M.D.L.; Denis, B.; et al. Atypical cortical encoding of speech identifies children with Dyslexia versus Developmental Language Disorder. bioRxiv 2022. [Google Scholar] [CrossRef]
  157. Knyazev, G.G.; Schutter, D.J.; van Honk, J. Anxious apprehension increases coupling of delta and beta oscillations. Int. J. Psychophysiol. 2006, 61, 283–287. [Google Scholar] [CrossRef]
  158. Ismail, L.E.; Karwowski, W. Applications of EEG indices for the quantification of human cognitive performance: A systematic review and bibliometric analysis. PLoS ONE 2020, 15, e0242857. [Google Scholar] [CrossRef]
  159. Turker, S.; Kuhnke, P.; Jiang, Z.; Hartwigsen, G. Disrupted network interactions serve as a neural marker of dyslexia. Commun. Biol. 2023, 6, 1114. [Google Scholar] [CrossRef]
Figure 1. mRNA levels of stress-associated genes of interest as well as methylation status of regions 1D and 1F in the NR3C1 gene promoter. (A) Examination of mRNA levels NR3C1 NR3C2, GILZ and FKBP5, DD samples vs. control samples. (B) (i) Schematic representation of the CpG island of NR3C1 gene with individual sites for 1D and 1F alternative untranslated first exons. (ii) Methylation percentages of the 1D CpG and (iii) 1F CpG sites. Asterisk (*) indicates a p-value < 0.05 and a double asterisk (**) indicates a p-value < 0.01. Abbreviations: GR, glucocorticoid receptor; MR, mineralocorticoid receptor; GILZ, glucocorticoid-induced leucine zipper.
Figure 1. mRNA levels of stress-associated genes of interest as well as methylation status of regions 1D and 1F in the NR3C1 gene promoter. (A) Examination of mRNA levels NR3C1 NR3C2, GILZ and FKBP5, DD samples vs. control samples. (B) (i) Schematic representation of the CpG island of NR3C1 gene with individual sites for 1D and 1F alternative untranslated first exons. (ii) Methylation percentages of the 1D CpG and (iii) 1F CpG sites. Asterisk (*) indicates a p-value < 0.05 and a double asterisk (**) indicates a p-value < 0.01. Abbreviations: GR, glucocorticoid receptor; MR, mineralocorticoid receptor; GILZ, glucocorticoid-induced leucine zipper.
Brainsci 14 00139 g001
Figure 2. Gene interaction network of dyslexia and stress-associated genes. Target genesare shown in the inner circle, while genes associated with the query genes are shown in the outer one. The left panel presents the different types of interactions with respective color coding, depicted with lines connecting genes in the network: physical interactions (pink), predicted (orange), co-expression (purple), genetic interaction (green) and shared protein domains (yellow). The right panel presents gene functions, represented with colored strips within the respective genes: hormone-mediated signaling pathway (red), cellular response to steroid hormone stimulus (green) and steroid hormone-mediated signaling pathway (blue), cilium assembly (orange) and response to steroid hormone (purple).
Figure 2. Gene interaction network of dyslexia and stress-associated genes. Target genesare shown in the inner circle, while genes associated with the query genes are shown in the outer one. The left panel presents the different types of interactions with respective color coding, depicted with lines connecting genes in the network: physical interactions (pink), predicted (orange), co-expression (purple), genetic interaction (green) and shared protein domains (yellow). The right panel presents gene functions, represented with colored strips within the respective genes: hormone-mediated signaling pathway (red), cellular response to steroid hormone stimulus (green) and steroid hormone-mediated signaling pathway (blue), cilium assembly (orange) and response to steroid hormone (purple).
Brainsci 14 00139 g002
Table 1. Frequency bands.
Table 1. Frequency bands.
Frequency RangeFrequency Band
0.5–4 HzDelta
4–8 HzTheta
8–12 HzAlpha
12–20 HzBeta1
20–30 HzBeta2
30–45 HzGamma
Table 2. Demographic characteristics regarding the DD and control individuals. DASS21 questionnaire scores and statistical values of the depression, anxiety and stress parameters examined (Mann–Whitney tests) and allelic frequencies of the genotyped variants of interest, namely FKBP5 rs1360780, SCL6A4 5-HTTLPR and SLC6A4 rs25531, for both the DD and control individuals. Abbreviations: DASS21, Depression, Anxiety and Stress Scale—21 Items; FKBP5, FK506 Binding Protein 5; SLC6A4, Solute Carrier Family 6 Member 4.
Table 2. Demographic characteristics regarding the DD and control individuals. DASS21 questionnaire scores and statistical values of the depression, anxiety and stress parameters examined (Mann–Whitney tests) and allelic frequencies of the genotyped variants of interest, namely FKBP5 rs1360780, SCL6A4 5-HTTLPR and SLC6A4 rs25531, for both the DD and control individuals. Abbreviations: DASS21, Depression, Anxiety and Stress Scale—21 Items; FKBP5, FK506 Binding Protein 5; SLC6A4, Solute Carrier Family 6 Member 4.
Demographic InformationDyslexia (n = 10)Control (n = 10)Total (n = 20)p-Value
Age (years), mean ± SD21.3 ± 2.120.7 ± 1.221 ± 1.70.6168
Sex (M/F)4/62/86/14-
DASS21 ScoresDyslexia (n = 10)Control (n = 10)Total (n = 20)p-Value
Depression, mean ± SD17 ± 10.713.8 ± 11.815.4 ± 11.10.4445
Anxiety, mean ± SD15.4 ± 12.48.8 ± 5.412.1 ± 9.90.2212
Stress, mean ± SD19.4 ± 9.619.6 ± 12.119.5 ± 10.60.8975
Allelic FrequenciesDyslexia (n = 10)Control (n = 8)Total (n = 18)
FKBP5
(rs1360780)
C0.70.6250.7
T0.30.3750.3
SLC6A4
(5-HTTLRP)
L0.40.6250.5
S0.60.3750.5
SLC6A4
(rs25531)
A111
G000
Table 3. Statistically significant correlations between EEG PSD and mRNA levels of stress-associated genes. Asterisk (*) indicates a p-value < 0.05 and a double asterisk (**) indicates a p-value < 0.01.
Table 3. Statistically significant correlations between EEG PSD and mRNA levels of stress-associated genes. Asterisk (*) indicates a p-value < 0.05 and a double asterisk (**) indicates a p-value < 0.01.
GeneBrain AreaRhythmCorrelationSpearman’s rp-Value
Control
NR3C1T7Beta 1Negative−0.783<0.05 *
T7Beta 2Negative−0.700<0.05 *
P7Alpha 2Negative-0.750<0.05 *
O2DeltaPositive0.821<0.05 *
GILZP7Alpha 2Negative−0.917<0.01 **
P8ThetaPositive0.786<0.05 *
O1Beta 2Negative−0.827<0.05 *
FKBP5P7Beta 2Negative−0.883<0.01 **
P8Alpha 1Negative−0.929<0.01 **
DD
NR3C1P8ThetaPositive0.883<0.01 **
P8Beta 1Negative−0.717<0.05 *
GILZO1ThetaPositive0.881<0.01 **
O1DeltaPositive0.786<0.05 *
FKBP5T8DeltaPositive0.717<0.05 *
Table 4. Interacting genes and networks of the studied stress-related genes according to the GeneMANIA gene interaction network. Genes associated with reading, speech and language skills are in bold.
Table 4. Interacting genes and networks of the studied stress-related genes according to the GeneMANIA gene interaction network. Genes associated with reading, speech and language skills are in bold.
GeneInteracting GeneNetworks
NR3C1NR3C2Physical Interactions, Shared Protein Domains
CNTNAP2Genetic Interactions
PCNTGenetic Interactions
CYP19A1Physical Interactions, Genetic Interactions
DIP2AGenetic Interactions
CMIPGenetic Interactions
FKBP5Physical Interactions
S100BGenetic Interactions
PAXBP1Genetic Interactions
AKAP9Co-Expression
EPB41L3Genetic Interactions
CNTN2Genetic Interactions
CPEGenetic Interactions
NR3C2NR3C1Physical Interactions, Shared Protein Domains
FKBP5Physical Interactions
SLC6A4Co-Expression
MIB1Genetic Interactions
CNTN2Co-Expression
CPECo-Expression
FKBP5NR3C1Physical Interactions
NR3C2Physical Interactions
ROBO1Co-Expression
GILZCo-Expression
TSC22D1Co-Expression, Genetic Interactions
CNTNAP2Genetic Interactions
PRMT2Genetic Interactions
P2RY1Genetic Interactions
SLC6A4NR3C2Co-Expression
PRMT2Co-Expression
CNTNAP2Genetic Interactions
MRFAP1Genetic Interactions
CNTN2Genetic Interactions
GILZFKBP5Co-Expression
AKR1B1Co-Expression
EPB41L3Co-Expression
ROBO1Genetic Interactions
TSC22D1Genetic Interactions, Shared Protein Domains
TSC22D2Physical Interactions, Shared Protein Domains
TSC22D4Physical Interactions, Shared Protein Domains
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Theodoridou, D.; Tsiantis, C.-O.; Vlaikou, A.-M.; Chondrou, V.; Zakopoulou, V.; Christodoulides, P.; Oikonomou, E.D.; Tzimourta, K.D.; Kostoulas, C.; Tzallas, A.T.; et al. Developmental Dyslexia: Insights from EEG-Based Findings and Molecular Signatures—A Pilot Study. Brain Sci. 2024, 14, 139. https://doi.org/10.3390/brainsci14020139

AMA Style

Theodoridou D, Tsiantis C-O, Vlaikou A-M, Chondrou V, Zakopoulou V, Christodoulides P, Oikonomou ED, Tzimourta KD, Kostoulas C, Tzallas AT, et al. Developmental Dyslexia: Insights from EEG-Based Findings and Molecular Signatures—A Pilot Study. Brain Sciences. 2024; 14(2):139. https://doi.org/10.3390/brainsci14020139

Chicago/Turabian Style

Theodoridou, Daniela, Christos-Orestis Tsiantis, Angeliki-Maria Vlaikou, Vasiliki Chondrou, Victoria Zakopoulou, Pavlos Christodoulides, Emmanouil D. Oikonomou, Katerina D. Tzimourta, Charilaos Kostoulas, Alexandros T. Tzallas, and et al. 2024. "Developmental Dyslexia: Insights from EEG-Based Findings and Molecular Signatures—A Pilot Study" Brain Sciences 14, no. 2: 139. https://doi.org/10.3390/brainsci14020139

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop