Takashi X. Fujisawa1,2*
Shota Nishitani1,2
Kai Makita1Akiko Yao1,2
Shinichiro Takiguchi3Shoko Hamamura2,3
Koji Shimada1,2,4
Hidehiko Okazawa1,2,4
Hideo Matsuzaki1,2,3Akemi Tomoda1,2,3*
- 1Research Center for Child Mental Development, University of Fukui, Fukui, Japan
- 2Division of Developmental Higher Brain Functions, United Graduate School of Child Development, Osaka University, Kanazawa University, Hamamatsu University School of Medicine, Chiba University, and University of Fukui, Osaka, Japan
- 3Department of Child and Adolescent Psychological Medicine, University of Fukui Hospital, Fukui, Japan
- 4Biomedical Imaging Research Center, University of Fukui, Fukui, Japan
The present study examined the relationship between DNA methylation differences and variations in brain structures involved in the development of attention-deficit hyperactivity disorder (ADHD). First, we used monozygotic (MZ) twins discordant (2 pairs of 4 individuals, 2 boys, mean age 12.5 years) for ADHD to identify candidate DNA methylation sites involved in the development of ADHD. Next, we tried to replicate these candidates in a case-control study (ADHD: N = 18, 15 boys, mean age 10.0 years; Controls: N = 62, 40 boys, mean age 13.9 years). Finally, we examined how methylation rates at those sites relate to the degree of local structural alterations where significant differences were observed between cases and controls. As a result, we identified 61 candidate DNA methylation sites involved in ADHD development in two pairs of discordant MZ twins, among which elevated methylation at a site in the sortilin-related Vps10p domain containing receptor 2 (SorCS2) gene was replicated in the case-control study. We also observed that the ADHD group had significantly reduced gray matter volume (GMV) in the precentral and posterior orbital gyri compared to the control group and that this volume reduction was positively associated with SorCS2 methylation. Furthermore, the reduced GMV regions in children with ADHD are involved in language processing and emotional control, while SorCS2 methylation is also negatively associated with emotional behavioral problems in children. These results indicate that SorCS2 methylation might mediate a reduced GMV in the precentral and posterior orbital gyri and therefore influence the pathology of children with ADHD.
Introduction
Attention-deficit hyperactivity disorder (ADHD) is one of the most common mental disorders in childhood, characterized by inattention, hyperactivity, and impulsivity, according to the 5th edition of the Diagnostic and Statistical Manual of Mental Disorders (DSM-5) (American Psychiatric Association, 2013) and often reaching into adulthood. The prevalence of ADHD in children worldwide was estimated to be 7.2% in a meta-analysis of 175 studies (Thomas et al., 2015). Patients with ADHD have difficulties in various cognitive domains, such as cognitive control, attention, timing, and working memory (Biederman et al., 1991; Rubia, 2018), as well as in other domains involved in emotional processing, such as motivation and timing processing, such as timing dissociation and delay-related impairments (Biederman et al., 1993; Sonuga-Barke et al., 2010; Rubia, 2018). According to a prospective follow-up study, approximately 50% of children with ADHD continue to have symptoms until adulthood, and if left untreated, they can be at higher risk of psychiatric problems such as depression, substance abuse, and social problems such as unemployment and criminal offenses (Biederman et al., 2006; Molina et al., 2009).
While the heritability of ADHD has been reported to be as high as 72–88% (Larsson et al., 2014), there has been obvious discordance in ADHD diagnosis between monozygotic (MZ) twin pairs, and often differences in severity within MZ concordant cases (Larsson et al., 2014), suggesting that epigenetic factors may be involved in the etiology. Epigenetic modifications regulate gene expression independently of changes in DNA sequence, primarily through DNA methylation and histone modifications (Henikoff and Matzke, 1997), and have been suggested to serve as a critical link between external environmental factors and long-lasting phenotypic changes (Malki et al., 2016). Epigenetic changes in the brain have been found to be involved in cognitive neurological processes, including psychiatric disorders, neurogenesis, and brain development. Candidate gene studies on ADHD’s DNA methylation profile based on peripheral samples such as blood or saliva have shown different methylation patterns of genes involved not only in dopaminergic, serotonergic, and neurotrophic systems including SLC6A4, DRD4, COMT, BDNF, and NGFR, but also neurotransmitter release or neurite outgrowth including ERC2 and CREB5 and associated with the symptoms and severity of ADHD (van Mil et al., 2014; Park et al., 2015; Xu et al., 2015; Dadds et al., 2016; Heinrich et al., 2017; Sengupta et al., 2017; Neumann et al., 2020).
Using disease discordant MZ twins for comparison in epidemiological epigenetic studies would be an ideal strategy, because the sex, age, perinatal environment, and other shared environmental factors that significantly influence the epigenome should be matched within MZ twins (Bell and Spector, 2011). Recent findings have revealed considerable epigenetic differences between MZ twins (Kaminsky et al., 2009), and such differences have been associated with phenotypic discordance between MZ twins, including psychiatric disorders (Kuratomi et al., 2008; Sugawara et al., 2011; Wong et al., 2014; Malki et al., 2016). Regarding ADHD, Chen et al. (2018) recently examined the relationship between brain structure and whole blood DNA methylation in 14 pairs of MZ-discordant cases, finding structural alterations in the striatum and cerebellum, as well as significant epigenetic differences in genes, such as γ-aminobutyric acid (GABA), dopamine and serotonin neurotransmitter systems, in these “discordant” brain structures (Chen et al., 2018). These findings support the role of DNA methylation in ADHD. However, given the high heterogeneity of ADHD, not only the study of DNA methylation associated significantly different brain regions in MZ discordant twins, but also examining the association between methylation array data analysis and structural alterations in the whole brain, including the cerebral cortex, can comprehensively elucidate the complex pathology of ADHD.
Here, we examined the relationship between DNA methylation differences on array data and variations in brain structures involved in the development of ADHD. Thus, our main hypothesis was that the DNA methylation sites nominated using identical MZ twins discordant for ADHD are ADHD-specific, and significantly associated with brain structures and symptoms observed in children with ADHD compared to typically developing children.
Materials and Methods
Participants
Two pairs of MZ twins discordant for ADHD were recruited from the Department of Child and Adolescent Psychological Medicine at the University of Fukui Hospital. The twins were 9-year-old males (pair 1) and 16-year-old females (pair 2). Eighteen children with ADHD (16 males and 2 females, mean age = 9.7 ± 1.6 years) were also recruited at the department for the case-control study. The diagnosis of ADHD was assessed by licensed child and adolescent psychiatrists according to DSM-5 criteria (American Psychiatric Association, 2013). Participants were also administered an assessment module of DSM-IV ADHD from the Schedule of Affective Disorders and Schizophrenia for School-Age Children, Epidemiologic version (K-SADS-E; Orvaschel and Puig-Antich, 1994). To further assess the core symptoms of ADHD (e.g., inattentive and hyperactive/impulsive symptoms), for the pairs of MZ twins discordant for ADHD parents were asked to complete the ADHD Rating Scale (ADHD-RS) (DuPaul et al., 1998) for all children with ADHD in the case-control study, the Swanson, Nolan, and Pelham Rating Scale (SNAP-IV) (Swanson, 1992). To exclude other psychiatric conditions (e.g., anxiety disorder), subjects were administered the Mini-International Neuropsychiatric Interview for Children and Adolescents (MINI-KID; Sheehan et al., 2010) by two licensed pediatric-psychological clinicians. Two existing Cohorts of 62 children (Cohort 1: n = 28, 21 males and 7 females, mean age = 14.9 ± 1.8 years; Cohort 2: n = 34, 19 males and 15 females, mean age = 13.1 ± 2.9 years) recruited from the local community in our previous study were used as controls for the case-control study (Takiguchi et al., 2015; Shimada et al., submitted1). All children had normal or corrected vision and normal hearing. All children, with the exception of the non-ADHD twins, were assessed using the Wechsler Intelligence Scale for Children-Fourth Edition (WISC-IV; Wechsler, 2003) or the Wechsler Adult Intelligence Scale-Third Edition (WAIS-III; Wechsler, 1997) and excluded if they had a full-scale intelligence quotient (FSIQ) <70. They were also excluded if they had any history of substance abuse, recent substance use, head trauma with loss of consciousness, significant fetal exposure to alcohol or drugs, perinatal or neonatal complications, neurological disorders, or medical conditions that might adversely affect growth and development. In the case-control study, behavioral and emotional problems were assessed using the Child Behavior Checklist (CBCL) in all children with ADHD and controls in Cohort 1 (Achenbach and Rescorla, 2001), and with the Strength and Difficulties Questionnaire (SDQ) in controls in Cohort 2 (Goodman, 1997).
Saliva Collection and DNA Extraction
Saliva samples were collected using the Oragene Discover OGR-500 kit (DNA Genotek Inc., Ottawa, ON, Canada). DNA was extracted using prepIT®•L2P reagent (DNA Genotek Inc.) and quantified using the Qubit™ dsDNA HS Assay Kit (Thermo Fisher Scientific Inc., Pittsburgh, PA, United States).
DNA Methylation Array
Genomic DNA (500 ng) was bisulfite-treated for cytosine-to-thymine conversion using the EZ DNA Methylation-Gold kit (Zymo Research, Irvine, CA, United States). The DNA was then whole-genome amplified, fragmented, and hybridized to the Human MethylationEPIC BeadChip (Illumina Inc., San Diego, CA, United States). BeadChips were scanned using iSCAN (Illumina Inc.), and the methylation level (β value) was calculated for each queried CpG locus using the GenomeStudio Methylation Module software, followed by the Psychiatric Genomics Consortium-Epigenome-Wide Association Studies quality control pipeline (Ratanatharathorn et al., 2017). Using CpGassoc (Barfield et al., 2012), samples with probe detection call rates <90% and those with an average intensity value of either <50% of the experiment-wide sample mean or <2,000 arbitrary units were excluded. Probes with detection P > 0.001 or those based on <3 beads were set to missing as were probes cross-hybridizing between autosomes and sex chromosomes (Teschendorff et al., 2013). CpG sites with missing data for >10% of samples within the dataset were excluded from the analysis. Probes containing single nucleotide polymorphisms (based on 1000 Genomes) within 10 base pairs of the target CpG were maintained in each dataset but flagged and tracked throughout the analysis pipeline. This decision was based on the growing recognition that sequence variants can influence DNA methylation patterns throughout the genome (Smith et al., 2014). Normalization of probe distribution and background differences between Type I and Type II probes was conducted using beta mixture quantile normalization (Teschendorff et al., 2013) after background correction. We did not remove the batch effect at this stage either for (1) MZ twins discordant for ADHD (Proband: N = 2, Non-proband: N = 2) and (2) ADHD cases (N = 18)–controls (N = 62) study since (1) those samples were scanned within the same chip and the row positional balance was identical [Pair 1: row 5 (Proband) vs. row 6 (Non-proband), pair 2: row 7 (Non-proband) vs. row 8 (Proband)], and (2) batches were completely confounded with case-control group status [Case chip ID (6 batches): 205111140162, 205111140170, 205111140171, 205134980172, 205134980191, and 205134980192; Control chip ID (8 batches): 203748260078, 203748260085, 203755070101, 203755080004, 203757350003, 203757350018, 203757350022, and 203757350023]. In such a case, it is not possible to remove technical signals when batches are confounded with variables of interest, even by employing tools such as ComBat (Johnson et al., 2007). As suggested by Nygaard et al. (2016) and Price and Robinson (2018), we decided to use chips and rows as additional covariates in our linear model instead of adjusting for batch effects in the initial processing to avoid P-value inflation. After quality control, 807,253 probes and 794,661 probes remained for (1) MZ twins discordant for ADHD and (2) the ADHD case-control study, respectively. We confirmed whether pair 1 and pair 2 were MZ twins using 59 “rs” probes on the EPIC chip using the R package ewastool (Heiss and Just, 2018), and found an identical genetic background (agreement: 0.9999891 and 0.9999893, respectively). As saliva contains a heterogeneous mixture of cell types of differing proportions in each sample, we used the EpiDISH method (Teschendorff et al., 2017) to estimate the proportion of epithelial cells derived from salivary DNA and entered it as a covariate in our statistical models.
Brain-Image Acquisition and Pre-processing
Image acquisition in the 52 participants in the case-control study (18 with ADHD, 34 controls in Cohort 2) was performed using a GE Signa PET/MR 3-Tesla scanner with an 8-channel head coil (GE Healthcare, Milwaukee, WI, United States). A T1-weighted anatomical dataset was obtained using a fast spoiled-gradient recalled imaging sequence (voxel size 1 × 1 × 1 mm, TE = 3.24 ms, TR = 8.46 ms, flip angle = 11°). Image acquisition for the other 28 controls in Cohort 1 in the case-control study was performed using a GE Discovery MR 750 3-Tesla scanner with a 32-channel head coil (GE Healthcare, Milwaukee, WI, United States). A T1-weighted anatomical dataset was obtained from each subject by a fast-spoiled gradient recalled imaging sequence (voxel size 1 × 1 × 1 mm, TE = 1.99 ms, TR = 6.38 ms, flip angle = 11°). VBM was performed as a global analytic approach using the Statistical Parametric Mapping version 12 software2 (Wellcome Department of Imaging Neuroscience, University College London, London, United Kingdom) implemented in MATLAB 2020b (Math Works Inc., Natick, MA, United States). T1-weighted images were segmented coarsely into gray matter (GM), white matter, cerebrospinal fluid, and skull/scalp compartments using tissue probability maps. The Diffeomorphic Anatomical Registration through Exponentiated Lie Algebra algorithm was applied to the segmented brain tissues to generate a study-specific template and to achieve an accurate inter-subject registration with improved realignment of smaller inner structures (Ashburner, 2007). The segmented GM images were spatially normalized, and written out with an isotropic voxel resolution of 1.5 mm. Any volume change induced by normalization was adjusted via a modulation algorithm. Spatially normalized GM images were smoothed by a Gaussian kernel of 6.2 mm full width at half maximum.
Statistical Analyses
First, to clarify epigenetic associations between proband and non-proband ADHD discordant MZ twins from methylation array data, multiple regression analysis was performed using CpGassoc (Barfield et al., 2012). In this analysis, DNA methylation at each CpG probe was entered as a dependent variable, and each group (proband or non-proband) entered as an independent variable. The proportion of epithelial cells was entered as a covariate, but we did not use age and sex as covariates because they were identical between groups.
Second, to confirm the reproducibility of the probes from the MZ twin discordant pair analysis in the case-control analysis, we examined the subset probes threshold set at P < 5.0E-05. DNA methylation at each CpG probe was entered as a dependent variable, and each group (case or control) as an independent variable. Age, sex, FSIQ, the proportion of epithelial cells, chip, and row for batch effect adjustments as explained previously, were entered as covariates, and results were threshold at false discovery rate (FDR) <0.05 by Benjamini-Hochberg.
Third, regional differences in gray matter volume (GMV) between groups were analyzed in SPM 12 using two-sample t-test models. Potential confounding effects of age, sex, FSIQ, scanner, and total GMV were modeled, and their attributed variances excluded from further analysis. Total GMV was calculated from the GM images obtained from pre-processing segmentation using the “Tissue Volumes” utility from the batching system in SPM12. A GM majority optimal threshold mask, created based on a study-specific sample, was applied to the analyses to eliminate voxels of non-GM for GMV-analyses (Ridgway et al., 2009). The resulting set of voxel values used for comparison generated a statistical parametric map of the t-statistic SPM{t} that was transformed to a unit normal distribution (SPM{Z}). The statistical threshold was set at P < 0.001 at the voxel level and P < 0.05, with a family wise error (FWE) correction for multiple comparisons. The anatomical localization of significant clusters was investigated using automated anatomical labeling and Brodmann area atlases implemented in the MRIcron software package (Rorden et al., 2007).
Finally, to further examine whether the ADHD-related GMV alterations were associated with DNA methylation, a correlation analysis for the residuals of each β value (methylation, GMV) controlled by control variables was performed. To this end, the adjusted eigenvariates, representing linearly transformed estimates of GMV, were extracted from the identified cluster. The significance level was set at P < 0.05. All statistical analyses were performed with R 3.6.3 (R Core Team, 2020), SPM 12, IBM SPSS Statistics for Windows version 26.0. (Armonk, NY: IBM Corp.).
Meta-Analytic Decoding of Regional Function Using NeuroSynth
The functional properties of structural regions with alterations between groups were decoded using a large-scale database-informed meta-analytic approach as implemented in NeuroSynth (Yarkoni et al., 2011). A meta-analytic map associated with the identified region coordinates was derived. Further, the terms (excluding terms for brain regions) ranked by the z-score were visualized using an online word cloud generator3.
Results
Clinical Status in Monozygotic Twin Discordant Pairs
The clinical status of the MZ twin discordant pairs is shown in Table 1. First, regarding IQ, both pairs of children with ADHD were in the 25–75 percentile range, and no significant defects in cognitive ability were observed. Next, regarding ADHD symptoms, the inattentiveness score was >90th percentile in both ADHD children, and <75th percentile in both control children, suggesting that the inattention symptoms were significantly stronger in children with ADHD. The hyperactive/impulsive score was <50th percentile in both ADHD and control children, suggesting that there were no significant hyperactivity symptoms.
Table 1. Clinical status between monozygotic twin discordant pairs.
Methylation Array Data Analysis of Monozygotic Twin Discordant Pairs
No CpG probes were detected above the genome-wide significance level under the number of EPIC array probes (P < 9.0E-08) (Tsai and Bell, 2015; Saffari et al., 2018) as a natural consequence of the extremely small sample size. We then threshold at P < 5.0E-05 (−log(P) = 4.3) by visual inspection of the Q-Q plot because the top probes over the threshold had residuals that steeply deviated from the expected line (Supplementary Figure 1). Sixty-one probes were above the threshold (Table 2 and Supplementary Figure 2), which were confirmed to be associated with ADHD in an independent case-control dataset.
Table 2. Top 61 differentially methylated CpG sites identified in ADHD-discordant monozygotic twin pairs, ranked by statistical significance and mean Δβ (calculated as DNA methylation level of control twin minus ADHD twin).
Demographic and Questionnaire Data in the Case-Control Study
The ADHD and control groups were matched for sex and handedness, but there was a significant difference in age between groups (T(78) = 5.88, P < 0.001). A two-sample t-test was used to compare the total FSIQ and CBCL total scores between groups. Compared to the control group, the ADHD group showed lower FSIQ (T(78) = 3.05, P = 0.003) and higher levels of ADHD-related emotional and behavioral problems (CBCL total, T(44) = −8.90, P < 0.001), although data for CBCL was only available from Cohort 1. In addition, in the control group of Cohort 2, SDQ total scores were not significantly different from the mean of standard Japanese children using one-sample t-test (SDQ total, T(33) = 0.47, P = 0.641), suggesting no notable emotional and behavioral problems. Multiple regression analysis was also performed to examine the effect of age on the differences in FSIQ and CBCL total score between groups, with each variable as the dependent variable; the results showed that the effect of group was significant (FSIQ: β = −0.38, T = −2.96, P = 0.004; CBCL total: β = 0.93, T = 5.95, P < 0.001) while that of age was not (FSIQ: β = −0.10, T = −0.77, P = 0.445; CBCL total: β = 0.16, T = 0.99, P = 0.326). Regarding ADHD symptoms, the inattention score of SNAP-IV in the ADHD group showed clinically mild symptoms on average, and the hyperactivity/impulsivity and opposition/defiance scores did not reach clinically significant levels. Since the ratio of epithelial cells in saliva samples affects the estimated methylation level (β value), we analyzed for the difference between groups finding no significant difference between the ADHD and control groups. We also estimated GMV, WMV, and total brain volume in the ADHD and control groups, respectively, but found no differences between groups in any of these parameters (Table 3).
Table 3. Demographic and clinical characteristics of ADHD and control groups.
Case-Control Subset Methylation Analysis
We extracted the top 61 probes from the initial analysis, among which 60 probes meeting the quality control criteria were included in the case-control dataset. Among these 60 probes, three probes (cg03700121, cg18430990, and cg19670431) in the coiled-coil domain containing 86 (CCDC86), transmembrane protein 240 (TMEM240), and sortilin-related Vps10p domain containing receptor 2 (SorCS2) genes were significantly associated with ADHD (FDRs < 0.05, Table 4). However, two of these probes (cg03700121 and cg18430990) were inconsistent regarding the direction of the effect, with increased methylation in twin studies and decreased methylation in case-control studies, while only the probe involved in SorCS2 showed consistent results in both studies. Hence, we used cg19670431 for subsequent epigenetic imaging analyses.
Table 4. Three CpG sites replicated in a sample of ADHD cases vs. controls.
Structural Brain-Image Data in the Case-Control Study
A whole-brain analysis with FWE correction at the cluster level was conducted to examine regional differences in GMV between the two groups (ADHD: n = 18, Controls: n = 62). Compared with the control group, the ADHD group showed reduced GMV in the left precentral gyrus (BA6; MNI coordinates, x = −44, y = 5, z = 57; cluster size = 437 voxels, P = 0.043, FWE corrected cluster level; Figure 1A) and the right posterior orbital gyrus (BA47; MNI coordinates, x = 24, y = 21, z = −21; cluster size = 555 voxels, P = 0.016, FWE corrected cluster level; Figure 1B).
Figure 1. Brain regions with significantly larger gray matter volume in the control group compared to the ADHD group. The statistical threshold for the contrasts was set at voxel-level P < 0.001 uncorrected for height and cluster-level P < 0.05 family wise error rate corrected for multiple comparisons. The color bar denotes the t-statistic range. (A) Left precentral gyrus (PrG; BA6; MNI coordinates, x = –44, y = 5, z = 57; cluster size = 437 voxels). (B) Right posterior orbital gyrus (POrG; BA47; MNI coordinates, x = 24, y = 21, z = –21; cluster size = 555 voxels). Reverse inference of functional properties in the left PrG (C) and right POrG (D) as decoded by NeuroSynth. The font size represents the rank according to the strength of the relationship between regions and terms.
Reverse inference on the functional properties related to the local regions where structural differences between groups were observed showed that most of the terms related to the left precentral gyrus are related to language functions such as semantic memory, working memory, theory of mind, and motion imagination (Figure 1C); most of the terms related to the right posterior orbital gyrus are related to–mainly negative–emotional information and their regulation (Figure 1D).
Relationship Between Sortilin-Related Vps10p Domain Containing Receptor 2 Methylation and Brain Structural Alterations
Sortilin-related Vps10p domain containing receptor 2 methylation was both positively correlated with GMV within a cluster in the precentral gyrus and the posterior orbital gyrus (Figures 2A,B). This result suggests that the more methylated SorCS2 is, the larger the GMV of the precentral and posterior orbital gyri. In addition, to verify the tissue specificity of the methylation pattern, we examined the brain-saliva correlation for the CpG probe (cg19670431) identified using a web tool based on human samples (Braun et al., 2019) and also confirmed a trend toward a positive correlation (ρ = 0.37, P = 0.09).
Figure 2. Correlation between SorCS2 methylation and gray matter volume in the left precentral gyrus (A) and the right posterior orbital gyrus (B). The red line shows the regression line. Adjusted residuals with respective covariates were used as values for methylation and local GMV of SorCS2. Note that the adjusted residuals for each covariate were used as estimates of SorCS2 methylation and local GMV.
Relationship Between Clinical Symptom Scores, Sortilin-Related Vps10p Domain Containing Receptor 2 Methylation and Brain Structural Alterations
We performed a correlation analysis investigate the association between the neurobiological basis of both SorCS2 methylation and local GMV alterations, and clinical symptoms (core symptoms based on SNAP-IV and emotional behavior problems based on CBCL tests) associated with ADHD. As a result, a significant negative correlation was confirmed between SorCS2 methylation and CBCL total score (r = −0.444, P = 0.002), while no significant correlation was found with SNAP total score (r = −0.227, P = 0.366); these scores also showed no correlation with both local GMV alterations. These results suggest that SorCS2 methylation may be involved in regulating emotional behavioral problems in children rather than ADHD-specific core symptoms.
Discussion
This study investigated the relationship between DNA methylation differences based on array data and brain structure involved variations in the development of ADHD. First, we investigated two pairs of MZ twins discordant for ADHD and identified 61 candidates for DNA methylation sites involved in the development of ADHD. Next, using these candidates in a case-control study we found that children with ADHD had elevated methylation in the SorCS2 gene body region. Finally, we observed that the ADHD group had significantly reduced GMV in the precentral gyrus and posterior orbital gyrus compared to the control group and that this volume reduction was positively associated with SorCS2 methylation. In addition, the reduced GMV regions in children with ADHD are involved in language processing and emotional control, and SorCS2 methylation is also negatively associated with emotional behavioral problems in children. These results indicate that SorCS2 methylation might mediate a reduced GMV in the precentral and posterior orbital gyri and therefore influence the pathology of children with ADHD.
We suggested that SorCS2 methylation is involved in ADHD through methylation array analysis of MZ twin discordant cases and case-control groups, while previous epigenome-wide studies found no evidence that SorCS2 methylation is involved in either children or adults with ADHD (van Dongen et al., 2019; Neumann et al., 2020; Rovira et al., 2020). Although SorCS2 is known to play a crucial role in neuronal viability and function (Glerup et al., 2016), human epidemiological studies have reported that single nucleotide polymorphisms in SorCS2 are associated with the risk of developing psychiatric disorders such as ADHD (Alemany et al., 2015), bipolar disorder (Baum et al., 2008; Ollila et al., 2009), and schizophrenia (Christoforou et al., 2011). Although recent human genome-wide association studies have suggested that a gene set related to dopamine signaling is involved not only in ADHD alone but also in the comorbidity of ADHD with obesity and narcolepsy (Mota et al., 2020; Takahashi et al., 2020); in animal studies, lack of SorCS2 reportedly induces ADHD-like behavior by altering the novelty response to psychostimulants and altering the dopaminergic firing pattern of the ventral tegmental area (Olsen et al., 2021). Other recent studies have revealed the significant roles of SorCS2 in brain derived neurotrophic factor (BDNF)-dependent plasticity and for social memory formation by N-methyl-D-aspartic (NMDA) receptor trafficking in hippocampal neurons (Glerup et al., 2016; Yang et al., 2021); in parallel, working memory deficits (Martinussen et al., 2005; Sowerby et al., 2011) and the involvement of BDNF or NMDA receptor signaling in ADHD have also been suggested (Bergman et al., 2011; Chang et al., 2014). Taken together, our results suggest that the SorCS2 gene methylation found in this study may affect the development and certain symptoms of ADHD by affecting dopaminergic, BDNF, and/or NMDA receptor signaling pathways.
Sortilin-related Vps10p domain containing receptor 2 methylation was positively associated with GMV in the precentral and posterior orbital gyri in the ADHD group, suggesting that unmethylated SorCS2 may lead to lower GMV. These results replicate prior results of surface area reduction in the left precentral and right orbital gyri found in another cohort of children with ADHD (Jung et al., 2019), as well as previous findings of cortical thickness reduction in the precentral gyrus and orbital gyri in a large-scale clinical sample of children (Hoogman et al., 2019). Big data analysis of structural magnetic resonance imaging of about 6,800 children found that the precentral gyrus surface area was one of the prominent local areas negatively associated with ADHD symptoms (Owens et al., 2021). Reverse inference showed that the precentral gyrus was associated with language function, and that children with ADHD have a higher risk of language problems (Sowerby et al., 2011; Hawkins et al., 2016; Korrel et al., 2017) which contributes to poor academic functioning (Sciberras et al., 2014). Hence, these findings suggest that SorCS2 gene methylation may induce language-related difficulties in children with ADHD via reduced GMV in the precentral gyrus. Next, reductions in GMV and cortical thickness in the right orbital gyrus of individuals with ADHD have often been reported (Vaidya, 2012; Hoogman et al., 2019; Jung et al., 2019). Although the function of the orbital gyri appears related to emotional information and calculations, emotional information may play an important role in decision-making and executive function (Bechara et al., 2000; Rolls and Grabenhorst, 2008). Numerous imaging studies using emotional and executive function tasks have reported reduced functional activation in the right orbital gyrus (Epstein et al., 2007; Rubia et al., 2009; Cubillo et al., 2012; Godinez et al., 2015). Taken together, these findings suggest that SorCS2 gene methylation may induce both reduced GMV in the orbital gyrus and emotional behavioral problems, although the direct association between the two could not be confirmed in this study.
Several limitations of the present study should be noted and taken into consideration for future studies. First, the sample size in this study was relatively small, including only two pairs of discordant MZ twins and 18 subjects in the ADHD group in the case-control study. In particular, the CpG sites identified in the array analysis of MZ twin discordant cases did not reach genome-wide significance level due to the small sample size; we only selected the top candidates based on their statistics. Although it is necessary to replicate the results with a larger sample size, it was particularly difficult to recruit MZ twins discordant by a single institution. Next, because we analyzed salivary DNA methylation, our data may not necessarily reflect the state of the brain due to the tissue specificity of methylation patterns (Smith et al., 2015). Regarding the tissue specificity of DNA methylation, although we tried to validate our results using a web tool that can investigate the correlation between methylation of brain, blood, saliva, and buccal cells collected from the same living human (Braun et al., 2019), another way to overcome the issue of tissue specificity is to directly examine the association between the methylation profile of SorCS2 and the precentral or posterior orbital gyri in the postmortem brains of children with ADHD. Finally, in the case-control study, the children’s age and cognitive abilities did not match and the batches for array analysis and the MR scanners for brain imaging were confounding between groups because we used two existing cohorts as control group. Although these factors were used as control variables in the statistical analysis, it made it difficult to distinguish whether the association between SorCS2 methylation and local GMV reduction was involved in the pathophysiology of ADHD or derived from demographic factors such as age and general cognitive abilities or research artifacts by batch and scanner effects. Although no significant difference was observed between the two groups in age-sensitive brain volume, there was a significant difference in emotional behavioral problems even after controlling for age, thus, future studies need to match age and cognitive abilities and exclude batch and scanner effects. Despite these limitations, this study sheds light on the some of the pathological mechanisms of ADHD in that it suggests DNA methylation candidates associated with brain structure specific to children with ADHD.
In conclusion, this study suggests that DNA methylation of the SorCS2 gene may induce language-related and emotional behavioral problems via brain structure alterations specific to children with ADHD. Some pharmacological or psychosocial interventions that enhance SorCS2 gene methylation may improve ADHD symptoms by interfering with the GMV reduction in the precentral and posterior orbital gyri. In the future, the elucidation of the molecular mechanism of local brain volume changes induced by SorCS2 methylation will be useful for understanding the pathophysiology of ADHD.
Data Availability Statement
The DNA methylation microarray data of monozygotic twins discordant samples have been deposited in the Gene Expression Omnibus database (GEO) with the primary accession code GSE186339 (https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi? acc=GSE186339). Other data that support the results of this study are included in the article/Supplementary Material, further inquiries can be directed to the corresponding author.
Ethics Statement
The study protocol was approved by the Research Ethics Committee of the University of Fukui, Japan, and all procedures were conducted in accordance with the Declaration of Helsinki and the Ethical Guidelines for Clinical Studies of the Ministry of Health, Labour, and Welfare of Japan. All parents provided written informed consent for participation in the study.
Author Contributions
TF, SN, and AT contributed to the conception and design of the study, performed the experiments and statistical data analysis, and wrote the first draft of the manuscript. KM, AY, ST, SH, KS, HO, and HM contributed substantially to performing the experiments and collecting the data, revised the manuscript critically for intellectual content, and approved the submitted version. All authors contributed to the article and approved the submitted version.
Funding
This work was supported by JSPS KAKENHI grant numbers JP19K21755, JP19H00617 (to AT), JP18K02480, JP20H04995 (to TF), and AMED under grant number JP20gk0110052 (to AT). The funders had no role in the study design, data collection and analysis, decision to publish, or preparation of the manuscript.
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s Note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fnins.2021.799761/full#supplementary-material
Footnotes
- ^ Shimada, K., Hamamura, S., Takiguchi, S., Makita, K., Yao, A., Fujisawa, T. X., et al. (submitted). Developmental differences in resting-state functional connectivity among youth in small-group family-like residential care.
- ^ http://www.fil.ion.ucl.ac.uk/spm/software/spm12/
- ^ https://www.wordclouds.com
References
Achenbach, T. M., and Rescorla, L. A. (2001). Manual for the ASEBA School-Age Forms & Profiles. Burlington, VT: University of Vermont, Research Center for Children, Youth, & Families.
Alemany, S., Ribasés, M., Vilor-Tejedor, N., Bustamante, M., Sánchez-Mora, C., Bosch, R., et al. (2015). New suggestive genetic loci and biological pathways for attention function in adult attention-deficit/hyperactivity disorder. Am. J. Med. Genet. B Neuropsychiatr. Genet. 168, 459–470. doi: 10.1002/ajmg.b.32341
American Psychiatric Association (2013). Diagnostic and Statistical Manual of Mental Disorders, Fifth Edition (DSM-5). Washington, DC: American Psychiatric Association.
Ashburner, J. (2007). A fast diffeomorphic image registration algorithm. Neuroimage 38, 95–113. doi: 10.1016/j.neuroimage.2007.07.007
Barfield, R. T., Kilaru, V., Smith, A. K., and Conneely, K. N. (2012). CpGassoc: an R function for analysis of DNA methylation microarray data. Bioinformatics 28, 1280–1281. doi: 10.1093/bioinformatics/bts124
Baum, A. E., Akula, N., Cabanero, M., Cardona, I., Corona, W., Klemens, B., et al. (2008). A genome-wide association study implicates diacylglycerol kinase eta (DGKH) and several other genes in the etiology of bipolar disorder. Mol. Psychiatry 13, 197–207. doi: 10.1038/sj.mp.4002012
Bechara, A., Damasio, H., and Damasio, A. R. (2000). Emotion, decision making and the orbitofrontal cortex. Cereb Cortex 10, 295–307. doi: 10.1093/cercor/10.3.295
Bell, J. T., and Spector, T. D. (2011). A twin approach to unraveling epigenetics. Trends Genet. 27, 116–125.
Bergman, O., Westberg, L., Lichtenstein, P., Eriksson, E., and Larsson, H. (2011). Study on the possible association of brain-derived neurotrophic factor polymorphism with the developmental course of symptoms of attention deficit and hyperactivity. Int. J. Neuropsychopharmacol. 14, 1367–1376. doi: 10.1017/S1461145711000502
Biederman, J., Faraone, S. V., Spencer, T., Wilens, T., Norman, D., Lapey, K. A., et al. (1993). Patterns of psychiatric comorbidity, cognition, and psychosocial functioning in adults with attention deficit hyperactivity disorder. Am. J. Psychiatry 150, 1792–1798. doi: 10.1176/ajp.150.12.1792
Biederman, J., Mick, E., Surman, C., Doyle, R., Hammerness, P., Harpold, T., et al. (2006). A randomized, placebo-controlled trial of OROS methylphenidate in adults with attention-deficit/hyperactivity disorder. Biol. Psychiatry 59, 829–835. doi: 10.1016/j.biopsych.2005.09.011
Biederman, J., Newcorn, J., and Sprich, S. (1991). Comorbidity of attention deficit hyperactivity disorder with conduct, depressive, anxiety, and other disorders. Am. J. Psychiatry 148, 564–577. doi: 10.1176/ajp.148.5.564
Braun, P. R., Han, S., Hing, B., Nagahama, Y., Gaul, L. N., Heinzman, J. T., et al. (2019). Genome-wide DNA methylation comparison between live human brain and peripheral tissues within individuals. Transl. Psychiatry 9:47. doi: 10.1038/s41398-019-0376-y
Chang, J. P., Lane, H. Y., and Tsai, G. E. (2014). Attention deficit hyperactivity disorder and N-methyl-D-aspartate (NMDA) dysregulation. Curr. Pharm. Des. 20, 5180–5185. doi: 10.2174/1381612819666140110115227
Chen, Y. C., Sudre, G., Sharp, W., Donovan, F., Chandrasekharappa, S. C., Hansen, N., et al. (2018). Neuroanatomic, epigenetic and genetic differences in monozygotic twins discordant for attention deficit hyperactivity disorder. Mol. Psychiatry. 23, 683–690. doi: 10.1038/mp.2017.45
Christoforou, A., McGhee, K. A., Morris, S. W., Thomson, P. A., Anderson, S., McLean, A., et al. (2011). Convergence of linkage, association and GWAS findings for a candidate region for bipolar disorder and schizophrenia on chromosome 4p. Mol. Psychiatry 16, 240–242. doi: 10.1038/mp.2010.25
Cubillo, A., Halari, R., Smith, A., Taylor, E., and Rubia, K. (2012). A review of fronto-striatal and fronto-cortical brain abnormalities in children and adults with Attention Deficit Hyperactivity Disorder (ADHD) and new evidence for dysfunction in adults with ADHD during motivation and attention. Cortex 48, 194–215. doi: 10.1016/j.cortex.2011.04.007
Dadds, M. R., Schollar-Root, O., Lenroot, R., Moul, C., and Hawes, D. J. (2016). Epigenetic regulation of the DRD4 gene and dimensions of attention-deficit/hyperactivity disorder in children. Eur. Child Adolesc. Psychiatry 25, 1081–1089.
DuPaul, G. J., Power, T. J., Anastopoulos, A. D., and Reid, R. (1998). ADHD Rating Scale-IV: Checklists, Norms, and Clinical Interpretation. New York, NY: Guilford Press.
Epstein, J. N., Casey, B. J., Tonev, S. T., Davidson, M. C., Reiss, A. L., Garrett, A., et al. (2007). ADHD- and medication-related brain activation effects in concordantly affected parent-child dyads with ADHD. J. Child Psychol. Psychiatry 48, 899–913. doi: 10.1111/j.1469-7610.2007.01761.x
Glerup, S., Bolcho, U., Mølgaard, S., Bøggild, S., Vaegter, C. B., Smith, A. H., et al. (2016). SorCS2 is required for BDNF-dependent plasticity in the hippocampus. Mol. Psychiatry 21, 1740–1751. doi: 10.1038/mp.2016.108
Godinez, D. A., Willcutt, E. G., Burgess, G. C., Depue, B. E., Andrews-Hanna, J. R., and Banich, M. T. (2015). Familial risk and ADHD-specific neural activity revealed by case-control, discordant twin pair design. Psychiatry Res. 233, 458–465. doi: 10.1016/j.pscychresns.2015.07.019
Goodman, R. (1997). The strengths and difficulties questionnaire: a research note. J. Child Psychol. Psychiatry 38, 581–586. doi: 10.1111/j.1469-7610.1997.tb01545.x
Hawkins, E., Gathercole, S., and Astle, D. The Calm Team Holmes, J. (2016). Language problems and ADHD symptoms: how specific are the links? Brain Sci. 6:50. doi: 10.3390/brainsci6040050
Heinrich, H., Grunitz, J., Stonawski, V., Frey, S., Wahl, S., Albrecht, B., et al. (2017). Attention, cognitive control and motivation in ADHD: Linking event-related brain potentials and DNA methylation patterns in boys at early school age. Sci. Rep. 7:3823. doi: 10.1038/s41598-017-03326-3
Heiss, J. A., and Just, A. C. (2018). Identifying mislabeled and contaminated DNA methylation microarray data: an extended quality control toolset with examples from GEO. Clin. Epigen. 10:73. doi: 10.1186/s13148-018-0504-1
Henikoff, S., and Matzke, M. A. (1997). Exploring and explaining epigenetic effects. Trends Genet. 13, 293–295.
Hoogman, M., Muetzel, R., Guimaraes, J. P., Shumskaya, E., Mennes, M., Zwiers, M. P., et al. (2019). Brain imaging of the cortex in ADHD: a coordinated analysis of large-scale clinical and population-based samples. Am. J. Psychiatry 176, 531–542. doi: 10.1176/appi.ajp.2019.18091033
Johnson, W. E., Li, C., and Rabinovic, A. (2007). Adjusting batch effects in microarray expression data using empirical Bayes methods. Biostatistics 8, 118–127. doi: 10.1093/biostatistics/kxj037
Jung, M., Mizuno, Y., Fujisawa, T. X., Takiguchi, S., Kong, J., Kosaka, H., et al. (2019). The effects of COMT polymorphism on cortical thickness and surface area abnormalities in children with ADHD. Cereb Cortex. 29, 3902–3911. doi: 10.1093/cercor/bhy269
Kaminsky, Z. A., Tang, T., Wang, S. C., Ptak, C., Oh, G. H., Wong, A. H., et al. (2009). DNA methylation profiles in monozygotic and dizygotic twins. Nat. Genet. 41, 240–245. doi: 10.1038/ng.286
Korrel, H., Mueller, K. L., Silk, T., Anderson, V., and Sciberras, E. (2017). Research Review: language problems in children with attention-deficit hyperactivity disorder – a systematic meta-analytic review. J. Child Psychol. Psychiatry 58, 640–654. doi: 10.1111/jcpp.12688
Kuratomi, G., Iwamoto, K., Bundo, M., Kusumi, I., Kato, N., Iwata, N., et al. (2008). Aberrant DNA methylation associated with bipolar disorder identified from discordant monozygotic twins. Mol. Psychiatry 13, 429–441. doi: 10.1038/sj.mp.4002001
Larsson, H., Chang, Z., D’Onofrio, B. M., and Lichtenstein, P. (2014). The heritability of clinically diagnosed attention deficit hyperactivity disorder across the lifespan. Psychol. Med. 44, 2223–2229. doi: 10.1017/S0033291713002493
Malki, K., Koritskaya, E., Harris, F., Bryson, K., Herbster, M., and Tosto, M. G. (2016). Epigenetic differences in monozygotic twins discordant for major depressive disorder. Transl. Psychiatry. 6:e839. doi: 10.1038/tp.2016.101
Martinussen, R., Hayden, J., Hogg-Johnson, S., and Tannock, R. (2005). A meta-analysis of working memory impairments in children with attention-deficit/hyperactivity disorder. J. Am. Acad. Child Adolesc. Psychiatry. 44, 377–384. doi: 10.1097/01.chi.0000153228.72591.73
Molina, B., Hinshaw, S. P., Swanson, J. M., Arnold, L. E., Vitiello, B., Jensen, P. S., et al. (2009). The MTA at 8 years: prospective follow-up of children treated for combined-type ADHD in a multisite study. J. Am. Acad. Child Adolesc. Psychiatry. 48, 484–500. doi: 10.1097/CHI.0b013e31819c23d0
Mota, N. R., Poelmans, G., Klein, M., Torrico, B., Fernàndez-Castillo, N., Cormand, B., et al. (2020). Cross-disorder genetic analyses implicate dopaminergic signaling as a biological link between Attention-Deficit/Hyperactivity Disorder and obesity measures. Neuropsychopharmacology 45, 1188–1195. doi: 10.1038/s41386-019-0592-4
Neumann, A., Walton, E., Alemany, S., Cecil, C., González, J. R., Jima, D. D., et al. (2020). Association between DNA methylation and ADHD symptoms from birth to school age: a prospective meta-analysis. Transl. Psychiatry 10:398. doi: 10.1038/s41398-020-01058-z
Nygaard, V., Rødland, E. A., and Hovig, E. (2016). Methods that remove batch effects while retaining group differences may lead to exaggerated confidence in downstream analyses. Biostatistics 17, 29–39. doi: 10.1093/biostatistics/kxv027
Ollila, H. M., Soronen, P., Silander, K., Palo, O. M., Kieseppä, T., Kaunisto, M. A., et al. (2009). Findings from bipolar disorder genome-wide association studies replicate in a Finnish bipolar family-cohort. Mol. Psychiatry 14, 351–353. doi: 10.1038/mp.2008.122
Olsen, D., Wellner, N., Kaas, M., de Jong, I., Sotty, F., Didriksen, M., et al. (2021). Altered dopaminergic firing pattern and novelty response underlie ADHD-like behavior of SorCS2-deficient mice. Transl. Psychiatry 11:74. doi: 10.1038/s41398-021-01199-9
Orvaschel, H., and Puig-Antich, J. (1994). Schedule for Affective Disorders and Schizophrenia for School-Age Children—Epidemiologic version (K-SADS-E), Fifth Edition. Fort Lauderdale, FL: Nova Southeastern University.
Owens, M. M., Allgaier, N., Hahn, S., Yuan, D., Albaugh, M., Adise, S., et al. (2021). Multimethod investigation of the neurobiological basis of ADHD symptomatology in children aged 9-10: baseline data from the ABCD study. Transl. Psychiatry 11:64. doi: 10.1038/s41398-020-01192-8
Park, S., Lee, J. M., Kim, J. W., Cho, D. Y., Yun, H. J., Han, D. H., et al. (2015). Associations between serotonin transporter gene (SLC6A4) methylation and clinical characteristics and cortical thickness in children with ADHD. Psychol. Med. 45, 3009–3017. doi: 10.1017/S003329171500094X
Price, E. M., and Robinson, W. P. (2018). Adjusting for batch effects in DNA methylation microarray data, a lesson learned. Front. Genet. 9:83. doi: 10.3389/fgene.2018.00083
Ratanatharathorn, A., Boks, M. P., Maihofer, A. X., Aiello, A. E., Amstadter, A. B., Ashley-Koch, A. E., et al. (2017). Epigenome-wide association of PTSD from heterogeneous cohorts with a common multi-site analysis pipeline. Am. J. Med. Genet. B Neuropsychiatr. Genet. 174, 619–630. doi: 10.1002/ajmg.b.32568
R Core Team. (2017). R: A Language and Environment for Statistical Computing. Vienna: R Foundation for Statistical Computing. Available online at: https://www.R-project.org/
Ridgway, G. R., Omar, R., Ourselin, S., Hill, D. L. G., Warren, J. D., and Fox, N. C. (2009). Issues with threshold masking in voxel-based morphometry of atrophied brains. Neuroimage 44, 99–111. doi: 10.1016/j.neuroimage.2008.08.045
Rolls, E. T., and Grabenhorst, F. (2008). The orbitofrontal cortex and beyond: from affect to decision-making. Prog. Neurobiol. 86, 216–244. doi: 10.1016/j.pneurobio.2008.09.001
Rorden, C., Karnath, H.-O., and Bonilha, L. (2007). Improving lesion-symptom mapping. J. Cogn. Neurosci. 19, 1081–1088. doi: 10.1162/jocn.2007.19.7.1081
Rovira, P., Sánchez-Mora, C., Pagerols, M., Richarte, V., Corrales, M., Fadeuilhe, C., et al. (2020). Epigenome-wide association study of attention-deficit/hyperactivity disorder in adults. Transl. Psychiatry 10:199. doi: 10.1038/s41398-020-0860-4
Rubia, K. (2018). Cognitive neuroscience of attention deficit hyperactivity disorder (ADHD) and its clinical translation. Front. Hum. Neurosci. 12:100. doi: 10.3389/fnhum.2018.00100
Rubia, K., Halari, R., Christakou, A., and Taylor, E. (2009). Impulsiveness as a timing disturbance: neurocognitive abnormalities in attention-deficit hyperactivity disorder during temporal processes and normalization with methylphenidate. Philos. Trans. R. Soc. Lond. B Biol. Sci. 364, 1919–1931. doi: 10.1098/rstb.2009.0014
Saffari, A., Silver, M. J., Zavattari, P., Moi, L., Columbano, A., Meaburn, E. L., et al. (2018). Estimation of a significance threshold for epigenome-wide association studies. Genet. Epidemiol. 42, 20–33. doi: 10.1002/gepi.22086
Sciberras, E., Mueller, K. L., Efron, D., Bisset, M., Anderson, V., Schilpzand, E. J., et al. (2014). Language problems in children with ADHD: a community-based study. Pediatrics 133, 793–800. doi: 10.1542/peds.2013-3355
Sengupta, S. M., Smith, A. K., Grizenko, N., and Joober, R. (2017). Locus-specific DNA methylation changes and phenotypic variability in children with attention-deficit hyperactivity disorder. Psychiatry Res. 256, 298–304. doi: 10.1016/j.psychres.2017.06.048
Sheehan, D. V., Sheehan, K. H., Shytle, R. D., Janavs, J., Bannon, Y., Rogers, J. E., et al. (2010). Reliability and validity of the mini international neuropsychiatric interview for children and adolescents (MINI-KID). J. Clin. Psychiatry 71, 313–326. doi: 10.4088/JCP.09m05305whi
Smith, A. K., Kilaru, V., Klengel, T., Mercer, K. B., Bradley, B., Conneely, K. N., et al. (2015). DNA extracted from saliva for methylation studies of psychiatric traits: evidence tissue specificity and relatedness to brain. Am. J. Med. Genet. B Neuropsychiatr. Genet. 168B, 36–44. doi: 10.1002/ajmg.b.32278
Smith, A. K., Kilaru, V., Kocak, M., Almli, L. M., Mercer, K. B., Ressler, K. J., et al. (2014). Methylation quantitative trait loci (meQTLs) are consistently detected across ancestry, developmental stage, and tissue type. BMC Genom. 15:145. doi: 10.1186/1471-2164-15-145
Sonuga-Barke, E., Bitsakou, P., and Thompson, M. (2010). Beyond the dual pathway model: evidence for the dissociation of timing, inhibitory, and delay-related impairments in attention-deficit/hyperactivity disorder. J. Am. Acad. Child Adolesc. Psychiatry 49, 345–355. doi: 10.1016/j.jaac.2009.12.018
Sowerby, P., Seal, S., and Tripp, G. (2011). Working memory deficits in ADHD: the contribution of age, learning/language difficulties, and task parameters. J. Atten. Disord. 15, 461–472. doi: 10.1177/1087054710370674
Sugawara, H., Iwamoto, K., Bundo, M., Ueda, J., Miyauchi, T., Komori, A., et al. (2011). Hypermethylation of serotonin transporter gene in bipolar disorder detected by epigenome analysis of discordant monozygotic twins. Transl. Psychiatry. 1:e24. doi: 10.1038/tp.2011.26
Swanson, J. M. (1992). School-Based Assessments and Interventions for ADD Students. Irvine, CA: KC Publishing.
Takahashi, N., Nishimura, T., Harada, T., Okumura, A., Choi, D., Iwabuchi, T., et al. (2020). Polygenic risk score analysis revealed shared genetic background in attention deficit hyperactivity disorder and narcolepsy. Transl. Psychiatry. 10:284. doi: 10.1038/s41398-020-00971-7
Takiguchi, S., Fujisawa, T. X., Mizushima, S., Saito, D. N., Okamoto, Y., Shimada, K., et al. (2015). Ventral striatum dysfunction in children and adolescents with reactive attachment disorder: functional (MRI) study. BJPsych. Open 1, 121–128. doi: 10.1192/bjpo.bp.115.001586
Teschendorff, A. E., Breeze, C. E., Zheng, S. C., and Beck, S. (2017). A comparison of reference-based algorithms for correcting cell-type heterogeneity in Epigenome-Wide Association Studies. BMC Bioinformatics 18:105. doi: 10.1186/s12859-017-1511-5
Teschendorff, A. E., Marabita, F., Lechner, M., Bartlett, T., Tegner, J., Gomez-Cabrero, D., et al. (2013). A beta-mixture quantile normalization method for correcting probe design bias in Illumina Infinium 450 k DNA methylation data. Bioinformatics 29, 189–196. doi: 10.1093/bioinformatics/bts680
Thomas, R., Sanders, S., Doust, J., Beller, E., and Glasziou, P. (2015). Prevalence of attention-deficit/hyperactivity disorder: a systematic review and meta-analysis. Pediatrics 135, e994–e1001. doi: 10.1542/peds.2014-3482
Tsai, P. C., and Bell, J. T. (2015). Power and sample size estimation for epigenome-wide association scans to detect differential DNA methylation. Int. J. Epidemiol. 44, 1429–1441. doi: 10.1093/ije/dyv041
Vaidya, C. J. (2012). Neurodevelopmental abnormalities in ADHD. Curr. Top. Behav. Neurosci. 9, 49–66. doi: 10.1007/7854_2011_138
van Dongen, J., Zilhão, N. R., Sugden, K., Bios Consortium, Hannon, E. J., Mill, J., et al. (2019). Epigenome-wide association study of attention-deficit/hyperactivity disorder symptoms in adults. Biol. Psychiatry 86, 599–607. doi: 10.1016/j.biopsych.2019.02.016
van Mil, N. H., Steegers-Theunissen, R. P., Bouwland-Both, M. I., Verbiest, M. M., Rijlaarsdam, J., Hofman, A., et al. (2014). DNA methylation profiles at birth and child ADHD symptoms. J. Psychiatr Res. 49, 51–59. doi: 10.1016/j.jpsychires.2013.10.017
Wechsler, D. (1997). Wechsler Adult Intelligence Scale, 3rd Edn. San Antonio, TX: Psychological Corporation.
Wechsler, D. (2003). Wechsler Intelligence Scale for Children, 4th Edn. San Antonio, TX: Psychological Corporation.
Wong, C. C., Meaburn, E. L., Ronald, A., Price, T. S., Jeffries, A. R., Schalkwyk, L. C., et al. (2014). Methylomic analysis of monozygotic twins discordant for autism spectrum disorder and related behavioural traits. Mol. Psychiatry 19, 495–503. doi: 10.1038/mp.2013.41
Xu, Y., Chen, X. T., Luo, M., Tang, Y., Zhang, G., Wu, D., et al. (2015). Multiple epigenetic factors predict the attention deficit/hyperactivity disorder among the Chinese Han children. J. Psychiatr. Res. 64, 40–50. doi: 10.1016/j.jpsychires.2015.03.006
Yang, J., Ma, Q., Dincheva, I., Giza, J., Jing, D., Marinic, T., et al. (2021). SorCS2 is required for social memory and trafficking of the NMDA receptor. Mol. Psychiatry 26, 927–940. doi: 10.1038/s41380-020-0650-7
Keywords: attention-deficit hyperactivity disorder (ADHD), DNA methylation, monozygotic twins, voxel-based morphometry (VBM), SorCS2
Citation: Fujisawa TX, Nishitani S, Makita K, Yao A, Takiguchi S, Hamamura S, Shimada K, Okazawa H, Matsuzaki H and Tomoda A (2022) Association of Epigenetic Differences Screened in a Few Cases of Monozygotic Twins Discordant for Attention-Deficit Hyperactivity Disorder With Brain Structures. Front. Neurosci. 15:799761. doi: 10.3389/fnins.2021.799761
Received: 22 October 2021; Accepted: 16 December 2021;
Published: 21 January 2022.
Edited by:
Noriyoshi Usui, Osaka University, JapanReviewed by:
Minyoung Jung, Korea Brain Research Institute, South KoreaEsther Walton, University of Bath, United Kingdom
Copyright © 2022 Fujisawa, Nishitani, Makita, Yao, Takiguchi, Hamamura, Shimada, Okazawa, Matsuzaki and Tomoda. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Takashi X. Fujisawa, tfuji@u-fukui.ac.jp; Akemi Tomoda, atomoda@u-fukui.ac.jp