Cortical Morphometric Abnormality and Its Association with Working Memory in Children with Attention-Deficit/Hyperactivity Disorder
Article information
Abstract
Objective
Attention-deficit/hyperactivity disorder (ADHD) is a common neurodevelopmental disorder in children and adolescents. The present study investigated the cortical morphology features and their relationship with working memory (WM).
Methods
In the present study, a total of 36 medication naïve children with ADHD (aged from 8 to 15 years) and 36 age- and gender-matched healthy control (HC) children were included. The digit span test was used to evaluate WM. The magnetic resonance imaging (MRI) was used to examine the characteristics of cortical morphology. Firstly, we compared the cortical morphology features between two groups to identify the potential structural alterations of cortical volume, surface, thickness, and curvature in children with ADHD. Then, the correlation between the brain structural abnormalities and WM was further explored in children with ADHD.
Results
Compared with the HC children, the children with ADHD showed reduced cortical volumes in the left lateral superior temporal gyrus (STG) (p=6.67×10-6) and left anterior cingulate cortex (ACC) (p=3.88×10-6). In addition, the cortical volume of left lateral STG was positively correlated with WM (r=0.36, p=0.029).
Conclusion
Though preliminary, these findings suggest that the reduced cortical volumes of left lateral STG may contribute to the pathogenesis of ADHD and correlate with WM in children with ADHD.
INTRODUCTION
Attention-deficit/hyperactivity disorder (ADHD) is one of the most common neurodevelopmental disorders in children and adolescents. ADHD occurs in 7.2% of children [1], characterized by impaired levels of inattention, hyperactivity, and impulsivity which can affect their academic and social development [2]. The exact aetiology and pathophysiology of ADHD remains unclear. There is increasing evidence showing that the abnormalities in brain structure may play an important role in the pathophysiology of this disorder [3,4].
Structural magnetic resonance imaging (sMRI) has been widely used to explore the differences in brain structure between childhood ADHD patients and healthy controls. Smaller total brain volume [5,6] and structural abnormalities in fronto-striatal, fronto-temporo-parietal and fronto-cerebellar networks have been found in children with ADHD [7]. As far as brain cortical morphology’s concerned, ADHD has been categorized as a delay in cortical maturation [8,9]. For example, significantly smaller prefrontal cortical surface areas were found in children with ADHD [8,10]. Along with a decrease in surface area of over 7% in both hemispheres, the children with ADHD also displayed a decrease in cortical folding bilaterally [11]. Another sMRI study has contended that children with ADHD had global thinning of the cortex, particularly in the prefrontal and precentral regions [12]. ADHD children had reduced cortical thickness in the left superior frontal, left orbito-frontal and left dorsal anterior cingulate cortex (ACC) and decreased gray matter volume in the left orbito-frontal, left middle frontal/dorsolateral prefrontal, left middle temporal and left cerebellum in comparison to the control group [13].
Brain structural abnormalities might lead to impairment of cognitive functions, including working memory (WM), which has been widely studied and implicated in the neuropsychological model of ADHD [14]. In the existing literature, the relationship between brain structural abnormalities and WM were revealed for multiple psychiatric disorders. For instance, cortical thickness and surface area in the frontal, temporal, cingular and occipital brain regions were found to be significantly associated with cognitive function among first-episode psychosis patients, including WM [15]. For bipolar disorder, negative relationship between right temporal pole thickness and WM was discovered [16]. Moreover, a recent study indicated that WM was associated with larger volume and surface in frontotemporal regions across the psychosis spectrum [17]. Higher gray matter volume in several regions of frontal lobe, cerebellum and temporal lobe was associated with better WM in adult ADHD [18]. As for children with ADHD, they often experience WM difficulties that may be associated with a worse outcome. It is those very common behaviors associated with WM deficits that might be the main complaints that lead children with ADHD to seek treatment in hospitals [19]. Whilst a number of studies have examined brain structural abnormalities in children with ADHD, far fewer studies have further explored the association between the altered brain structures and deficit WM. The only one study that recruited 17 children with ADHD for cortical thickness analyses denoted a direct relationship between thinner cortical thickness of the left medial temporal cortex and deficient WM [20]. In addition to cortical thickness, more dimensional brain cortical morphology indicators are worth exploring simultaneously to uncover the alteration in brain morphology associated with ADHD, such as cortical volume, surface area, and cortical curvature. Cortical volume is the product of two independent components, cortical thickness and surface area. Decreased cortical volume might be associated with either a cortical thinning of the cortex or a decrease in the cortical surface area (or a combination of both) [11]. Cortical thickness is measured as the distance from the cortical white-gray matter boundary to the cortical surface (gray-cerebrospinal fluid boundary) [21,22]. The surface defining the boundary between gray matter and cerebrospinal fluid, or pial surface, was determined by deforming the gray/white surface outward [23,24]. Cortical curvature assesses differences in small areas (peak of gyri and sulci) [22]. Studies on cortical volume showed that brain regions associated with cognition function were more likely to show changes in ADHD [25]. Reduced cortical surface area and thickness indicated that there is a delay in the maturation of brain in ADHD [8,26]. Cortical curvature may be an incidental finding of uncertain clinical significance [22]. In all, the analyses of these four different structural indicators simultaneously would help us to understand the characteristics of brain structural abnormalities in ADHD more comprehensively [11].
The present study builds upon the limited research on cortical morphology in children with ADHD by involving 36 children with ADHD and 36 age- and gender- matched healthy controls to investigate the cortical morphology features (cortical volume, surface area, and cortical curvature) and their relationship with WM. We only included medication-naïve patients in this study to eliminate potential confounding effects of medication [5]. We hypothesized that: 1) children with ADHD would show regional cortical morphology abnormalities, and 2) cortical morphology abnormalities would be associated with WM within the ADHD group.
METHODS
Participants
All ADHD cases were recruited from child psychiatric clinics of our hospital. Diagnoses of ADHD were conducted according to the DSM-IV criteria based on a semi-structured interview by experienced child psychiatrists using the Kiddie Schedule for Affective Disorders and Schizophrenia-Present and Lifetime Version (K-SADS-PL) [27]. The inclusion criteria were as follows: 1) DSM-IV ADHD diagnosis, 2) age between 8 and 15 years, 3) full-scale intelligence quotient (FSIQ) ≥80 as measured by the Wechsler Intelligence Scale for Chinese Children-Revised (WISCC-R) [28], 4) medication-naïve, 5) right-handedness, as assessed by the Chinese Handedness Inventory [29]. The exclusion criteria for ADHD group were comorbidity with major neurological disorders, a diagnosis of schizophrenia, pervasive developmental disorder, tic disorder, emotional disorders, affective disorders, epilepsy, mental retardation, or other brain disorders, except for Learning Disorder (LD) and Disruptive Behavior Disorder (DBD). Forty-one ADHD children participated in the study. Five cases were excluded from analysis because of poor scan quality. Then, a total of 36 (28 males and 8 females) ADHD were enrolled for final analyses, including 19 ADHD inattentive type (ADHD-I), 1 ADHD hyperactive-impulsive type (ADHD-HI) and 16 ADHD combined type (ADHD-C) subjects. The mean age was 10.5 years (SD=1.8).
Healthy control (HC) subjects were all students from local schools, subjected to the Kiddie Schedule for Affective Disorders and Schizophrenia-Present and Lifetime Version (K-SADS-PL) by the same group of child psychiatrists that evaluated the ADHD group. None of the normal control subjects met the criteria for ADHD or any other psychiatric disorders. Finally, 36 age- and gender- matched individuals (28 males and 8 females) comprised the HC group. All normal control group were right-handed. The mean age was 10.5 years (SD=1.6).
This study and its procedures were approved by the Ethics Committee of Peking University Institute of Mental Health. The institutional review board (IRB) number is PKUH6 (2007)-(39). After procedures had been explained, written informed consent was obtained from all subjects and their parents. All procedures have been performed in accordance with the ethical standards laid down in the 1964 Declaration of Helsinki and its later amendments.
Psychiatric and neurocognitive assessments
ADHD symptoms were evaluated using ADHD Rating Scale-IV (ADHD-RS-IV), which consists of 18 items according to DSM-IV for ADHD [30]. Each symptom was scored on how often it occurred (i.e., if “never” presented a symptom, it would be coded as 0; if “occasionally,” 1;“often,” 2; and “always,” 3. We employed the Digit Span Tests from the WISCCR to measure WM, involving two subtests: Digit Span Forward (DSF) and Digit Span Backward (DSB). DSF was used to assess simple verbal memory and attention span [31]. DSB was applied for auditory working memory manipulation, testing the maximum number series an individual can aurally attend to, store and then repeat verbally in reverse order [31]. DSF and DSB should be treated as different constructs. DSF emerged as a signifificant predictor for one measure of attention [32], while DSB was considered to be a measure of WM because it demanded both storage and manipulation in order to retain and repeat the number in reverse order [20,32]. Because the digit span (DS) subtest of the Wechsler Intelligence Scale consists of DSF and DSB components that produce separate raw scores but are combined to produce a single scaled score [32]. Raw scores of DSB can be used to compare group differences in WM [33]. Therefore, the raw score of DSB and scaled score of DS were utilized in our study [20,34].
Structural MRI data acquisition and Image processing
All imaging data were acquired on a Siemens Trio 3T scanner (Siemens, Erlangen, Germany) at the State Key Laboratory of Cognitive Neuroscience and Learning, Beijing Normal University. High-resolution 3D T1-weighted images were obtained using a magnetization-prepared rapid gradient echo (MP-RAGE) sequence: repetition time (TR)=2,530 ms; echo time (TE)=3.39 ms; inversion time=1,100 ms; flip angle=7°; field of view (FOV)=256×256 mm; matrix=256×256; 128 slices, thickness=1.33 mm. MRI data were analyzed with atlas-based FreeSurfer software (http://surfer.nmr.mgh.harvard.edu, version 5.1.0). Freesurfer adopted a fully automated method to perform pre-processing steps including Talairach alignment, intensity normalization, and removal of skull and non-brain tissue with a hybrid watershed/surface deformation procedure, separation of the cerebellum and brainstem from the cerebrum and splitting of the left and right hemispheres. Each hemisphere was divided into 74 different cortical regions in compliance with the Destrieux’s atlas [35]. Images with motion-related artifacts were excluded. As for the quality control of the images, we followed the ENIGMA pipeline (http://enigma.ini.usc.edu/). All brain cortical regions with a volume, surface, thickness, or curvature >1.5 or <1.5 times the interquartile range were identified and visually inspected by overlaying their segmentations on the subjects’ anatomical images. Only subjects with images for which segmentation were judged to be accurate were included to further analysis. Cortical volume, cortical thickness, cortical surface area, and curvature were measured applying Freesurfer on all the subjects.
Statistical analysis
Statistical analyses were conducted in SPSS, version 19 (IBM Corp., Armonk, NY, USA) and R (https://www.r-project.org/). Demographics and clinical characteristics were compared using two-tailed t-tests for continuous variables or chi-square (χ2) tests in SPSS for categorical variables between patients and controls. Cortical volume, surface, thickness, and curvature were examined by the modified generalized linear model (GLM) in R, with age, sex, and FSIQ as covariates. Modified GLM, also known as glm2 in R, is a modified version of GLM in the stats package. It is identical to GLM except for minor modifications to change the default fitting method and more stable compared with GLM [36]. Packages of “haven” and “glm2” in the RStudio (https://www.rstudio.com/products/rstudio/download/) were used in our present study to analyze the data in the sav file. To adjust the person-to-person variations in head size, total intracranial volume (ICV) (summed volumes of total grey matter, total white matter and cerebrospinal fluid compartment derived from Freesurfer’s segmentation procedure) was set as a covariate. The level of significance was defined as p<0.05. Considering the number of analyzed cortical regions (n=74), Bonferroni corrections were performed with setting the adjusted significance at p<6.8×10-4 (0.05/74).
We further explored the correlation between abnormal brain structures in children with ADHD and their WM using Pearson partial correlation with age, sex and ICV as covariates. A statistical threshold of p=0.05 was employed for all correlational analyses. Taking into account of the number of analyzed cortical regions (n=2) and the number of WM indicators (n=2), the adjusted significance was set to p<0.0125 (0.05/4).
RESULTS
Demographic and clinical variables
The demographic and clinical characteristics of children with ADHD and healthy controls were shown in Table 1.
The age [(10.5±1.8) vs. (10.5±1.6), p=0.891] and gender (χ2<0.01, p=1.000) of participants in both children with ADHD and healthy controls were found no significant difference. HC subjects had a significantly higher FSIQ score than children with ADHD [(117.5±15.1) vs. (105.4±14.6), p=0.001]. Compared with HC subjects, significant group differences were observed on DSB and DS with age and sex as covariates. Children with ADHD performed significantly worse on DSB [(4.64±1.31) vs. (5.57±1.58), p=0.001] and DS [(10.92±2.66) vs. (13.20±3.08), p=3.62×10-4] than controls.
Cortical morphologic abnormality
Children with ADHD had a significant reduction of cortical volume in the left lateral aspect of the superior temporal gyrus (STG) (p=6.67×10-6) and left ACC (p=3.88×10-4) compared with controls with age, sex, FSIQ and ICV adjusted. In addition, children with ADHD exhibited a nominal reduction of cortical surface in the right ACC, middle-anterior part of ACC and sulcus (aMCC), supramarginal gyrus, planum polare of the STG, temporal pole, medial orbital sulcus (olfactory sulcus), superior temporal sulcus, and the left aMCC, inferior segment of the circular sulcus of the insula, parietooccipital sulcus. However, these regions did not survive Bonferroni corrections (Table 2).
There were no significant between-group differences in cortical surface, thickness or curvature after correction for multiple comparisons. Nominal reduction of cortical surface was found for children with ADHD in the bilateral cingulate cortices and the left parietal and temporal cortices before correction for multiple (Supplementary Table 1 in the online-only Data Supplement). In addition, children with ADHD manifested a nominal reduction of cortical thickness in the bilateral frontal and temporal cortices before correction for multiple (Supplementary Table 2 in the online-only Data Supplement). Children with ADHD showed a nominal decreased of cortical curvature in the right right orbital part of the inferior frontal gyrus before multiple corrections (Supplementary Table 3 in the online-only Data Supplement).
Correlations between abnormal cortical volumes and WM
The volume of left lateral aspect of the STG volume was positively correlated with DSB scores after controlling for age, sex and ICV (r=0.36, p=0.029). To be specific, the bigger the volume in children with ADHD was, the higher DSB scores were and the better WM was (Table 3). However, their correlation could not survive Bonferroni corrections. There were no correlations between the volume of ACC and DSB/DS scores.
DISCUSSION
As expected, the results of our present study provided strong evidence that children with ADHD showed cortical morphology abnormalities relative to healthy controls. Children with ADHD had smaller cortical volumes in the left lateral STG and left ACC. The decreased volume of the left lateral STG might be correlated with impaired WM in children with ADHD.
Cortical volume in the left lateral STG
The enrollment of temporal gyrus in ADHD has been illustrated in previous studies. One recent study recruiting 18 children with ADHD showed decreased volume in the left middle temporal gyrus in children with ADHD [13]. Boys with ADHD showed significantly lower gray volume in the right temporal regions [37]. Some other studies also argued that children and adolescents with ADHD had reduced volume in the middle temporal gyrus [38], the bilateral temporal polar and occipital cortices and right superior temporal sulcus [25], or the right middle and inferior temporal gyri [39]. Young adults with ADHD showed reduced volume in the right superior and middle temporal gyrus [4] and higher fractional anisotropy in the right STG [40]. Taken together with these preceding evidence, our present finding of the reduced volume in the left STG supported the importance of structural alteration of temporal gyrus in the etiology of ADHD. However, there was inconsistency on the exact location, which might be due to limited sample size or a gender-by-diagnosis interaction effect in the volume of brain abnormalities [41].
It’s worth noting that our findings of abnormalities in left STG in children with ADHD parallel with some fMRI findings. Children with ADHD displayed increased brain connectivity in the left STG after four weeks of equine-assisted activity and training that improved clinical symptoms of ADHD [42]. Our study group has found that decreases in negative functional connectivity with the putamen-ROIs was seen in STG in medication-naïve children with ADHD relative to controls [43]. An event-related fMRI revealed significantly reduced brain activation in left superior temporal lobes in patients compared with healthy subjects, and that the activation differences in superior temporal lobes correlated inversely with response variability in control subjects but not in patients [44]. Another fMRI study claimed that increased fractional amplitude of low-frequency fluctuations (fALFF) in the left STG was related to reductions in inattentive symptoms for children with ADHD treated with methylphenidate [45]. Besides, adolescents with ADHD also presented increased activation in the STG during successful response inhibition [46]. Interestingly, one study on quantifying patterns of brain activity between ADHD and their unaffected siblings also announced that STG is associated with risk for ADHD [47]. In summary, fMRI studies have expressed that abnormal functional relationships between STG and other brain regions may underlie the pathological basis of ADHD [43,48,49]. The results showed that the left lateral STG which we identified structural deficits in the children ADHD group also manifested functional deficits. Further studies are needed to explore the neural alterations in children with ADHD by combining and fusing both brain structural and functional imaging data to achieve more comprehensive interpretation.
Cortical volume in the left ACC
The reduced volume in the left ACC offered strong evidence that the left ACC is abnormal in children with ADHD. One study provided consistent evidence of the structural deficits observed in the left ACC for untreated children with ADHD [5]. A recent sMRI study including 18 children with ADHD and 18 healthy children also has identified that patients showed smaller volumes in the left ACC and the reductions of volume were correlated with patients’ inattention [50]. A meta-analysis of MRI studies, collecting 11 VBM studies from 2001 to 2011 in all, has found that untreated children with ADHD tended to show additional changes in the left ACC compared with treated children and healthy children [51]. Interestingly, the meta-analysis also contended that adults with persistent ADHD symptoms were still characterized by volume reductions in the left and right ACC. Meanwhile, a study by Matsuo et al. [52] that adult ADHD patients presented volume in the left ACC. Several studies stated that structural abnormalities in ADHD occur before the age of 10 and can persist into adulthood [51,53,54]. Therefore, these results resonate well with our finding reported here, which suggests that ADHD children show reduced volume in the left ACC. Inconsistent with our results, one previous study has argued that only the right ACC was significantly smaller in the treatment-naive children with ADHD-C compared with the treated children or healthy children [55]. However, further analyses from the same researchers have identified that the ADHD-C group had bilaterally smaller volumes of the ACC compared with ADHD-I or healthy children [56]. This discrepancy might be attributable to ADHD subtypes [56] or the use of medication [5] effects in the volume of brain abnormalities. In order to exclude the possible influence of ADHD subtypes on the results, we initially analyzed the difference in cortical morphology among ADHD-I subtypes, ADHD-C subtypes and healthy controls. The preliminary analyses among three groups indicated differences in cortical volume in the left lateral aspect of the STG , with age, sex, FSIQ and ICV adjusted, which is consistent with the results between children with ADHD and HC subjects. Post-hoc analysis indicated that both ADHD-I and ADHD-C were with smaller cortical volume in the left lateral STG than HC, however, no significant difference was found between these two subtypes (the significant region is shown in the Supplementary Table 4 in the online-only Data Supplement, and visualized in the Supplementary Figure 1 in the online-only Data Supplement).
Furthermore, several functional MRI studies have indicated that children with ADHD demonstrated marked reduced activation in the ACC extending to the supplementary motor area during response inhibition [57] or during attentional alerting [58]. Likewise, a systematic review of fMRI and electrophysiological (electroencephalography, EEG) studies has declared that children with ADHD showed an enhanced activation in the ACC [59]. A recent event-related potential source imaging study revealed reduced activation in the ACC children with ADHD during the stop-signal task [60]. These findings largely echo the results in the current study. In a nutshell, our findings confirm the abnormality of the left ACC in children with ADHD.
Cortical thickness, surface area and curvature
In this study, there were no differences in cortical surface, thickness, and curvature between children with ADHD and healthy controls. For the cortical surface, smaller cortical surface areas bilaterally [11] and prefrontal cortical surface areas [8,10] have been indicated in children with ADHD; however, negative results have also been reported as in our present study [22]. The inconsistency also existed for the studies of cortical thickness. some found that children with ADHD had decreased cortical thickness in the prefrontal, precentral and temporal regions [12,22,61]; wheras, some not [11]. As for cortical curvature, there were very few group differences [22]. Another study also spotted no differences in cortical gyrification index (a measure comparable to curvature) between children with ADHD and healthy controls [8]. Multiple factors might influence the consistency and reproducibility among different studies. Firstly, the small sample size might influence the statistical power and lead to false negative or false positive findings. Secondly, the status of comorbidities should also be considered [6]. In addition, although brain cortical development in children with ADHD lagged behind that of healthy children, children with ADHD might reach the peak cortical thickness several years later [26]. So, the effects from development are needed to be investigated precisely in the future. In all, we need subsequent research with a larger sample size to evaluate the influence of comorbidities and development on cortical morphology.
Cortical morphology abnormalities and WM
In the current study, cortical volume in the left lateral STG was positively correlated with DSB scores in children with ADHD, suggesting that the left lateral STG might be associated with WM. Children with ADHD are likely to exhibit obvious WM deficits according to a number of studies [62,63]. Fortunately, WM can be improved after working memory training with a commercially available and computerized working memory program (Cogmed) in children with ADHD [64]. The STG is significantly associated with inattention among children with ADHD in one sMRI study [65]. Meanwhile, WM deficits are strongly related to symptoms of inattention in children with ADHD [66]. Lesions in lateral temporal lobes have been proved being able to disrupt WM processes [67]. A recent study has claimed that human lateral temporal cortex is part of the neural system for verbal WM [68], which is in line with our findings to some extent. To this day, there have been many fMRI studies confirming that WM in ADHD is associated with brain activation changes and abnormal functional connectivity, such as less activation in the left ventrolateral prefrontal cortex (VLPFC), cerebellar, occipital regions and lower connectivity in the ventrolateral prefrontal cortex, the ACC, the superior parietal lobule, and the cerebellum [63], less activity in right frontal, temporal, and subcortical regions and left occipital and cerebellar regions in male ADHD subjects during performance of a verbal working memory task [69], and Default Mode Network (DMN) functional abnormalities [70-72]. Much research attention has been paid to the association between the prefrontal cortex and WM in ADHD [73-75], whereas this study should add to the evidence that ADHD deficits in WM are linked to several specific brain regions including left STG [76]. However, future research with more detailed indicators of WM might help better understand the correlations involved. Furthermore, the correlation between cortical volume in the left lateral STG and DSB scores was only indicated in children with ADHD, but not in healthy controls. This might be related to several factors such as stronger variance in children with ADHD and clinical characteristics that are specific to the patients group.
In conclusion, this study confirms previous findings of cortical morphology abnormalities in children with ADHD, demonstrating that children with ADHD have regional cortical morphology abnormalities, especially smaller cortical volumes in the left lateral STG and left ACC. Moreover, volumes in the left lateral STG is associated with verbal WM. This study suggests that reduced cortical volumes are associated with ADHD, and may be essential for understanding the pathophysiology of the disorder.
Limitations
Several limitations must be considered for the present study. First, because of the relatively small sample size of this study, it was difficult to analyze the effects of ADHD subtypes on the cortical morphology. Therefore, a larger sample size would be necessary to strengthen our results. Previous work has found brain structural differences between ADHD-C group and ADHD-I group, informing that different ADHD subtypes may be different phenotypes [56]. In our study, although no main diagnostic effect was found for cortical thickness, the preliminary analyses for subtypes indicated that ADHD-I group showed significantly increased cortical thickness in the left superior temporal sulcus, when compared with ADHD-C group (data not shown). The expanded sample size in the future might also help us to illustrate the potential effects of subtypes on cortical morphology abnormalities in children with ADHD. Second, a gender-by-diagnosis interaction effect might exist on the volume of brain abnormalities [41], such as ACC [77]. Although sex has been matched for two groups in this present study, a larger sample size may still be needed for betweengender comparisons to further illustrate the gender effects on cortical morphology. Third, we could not exclude the potential confounding contribution of LD and DBD to the results of the present study. Further investigation will be desired to evaluate the effects of comorbidities on cortical morphology.
Supplementary Materials
The online-only Data Supplement is available with this article at https://doi.org/10.30773/pi.2020.0333.
Acknowledgements
This work was supported by the National Science Foundation of China (81571340), the National Key Basic Research Program of China (973 program 2015CB856405), Beijing Municipal Science & Technology Commission (No.Z161100000516032), the National Natural Science Foundation of China (81641163, 81761148026, 81873802), Beijing Natural Science Foundation (7172245), and the Capital Foundation of Medical Developments (CFMD:2016-2-4113). We thank all the participants for their kind contributions to our research. Thank Nan Li, Feng-Mei Fan, and Zhi-Jiang Wang for their help in statistical analyses.
Notes
The authors have no potential conflicts of interest to disclose.
Author Contributions
Conceptualization: Qiu-Jin Qian, Lu Liu. Data curation: Fei-Fei Si. Formal analysis: Fei-Fei Si, Lu Liu. Funding acquisition: Qiu-Jin Qian, Lu Liu. Methodology: Fei-Fei Si, Lu Liu, Hai-Mei Li, Li Sun, Qing-Jiu Cao. Project administration: Hai-Mei Li, Li Sun, Qing-Jiu Cao. Resources: Yu-Feng Wang, Li Sun, Qing-Jiu Cao. Software: Fei-Fei Si. Supervision: Qiu-Jin Qian, Lu Liu, Hanna Lu. Validation: Fei-Fei Si, Lu Liu. Writing—original draft: Fei-Fei Si. Writing—review & editing: Qiu-Jin Qian, Lu Liu, Hanna Lu.