{'fraction_non_alphanumeric': 0.03292804464819613, 'fraction_numerical': 0.013553916683276858, 'mean_word_length': 5.136497064579256, 'pattern_counts': {'":': 0, '<': 12, '<?xml version=': 0, '>': 3, 'https://': 1, 'lorem ipsum': 0, 'www.': 2, 'xml': 0}, 'pii_count': 0, 'substrings_counts': 0, 'word_list_counts': {'cursed_substrings.json': 27, 'profanity_word_list.json': 0, 'sexual_word_list.json': 0, 'zh_pornsignals.json': 0}}
Sex differences in structural brain asymmetry of children with autism spectrum disorders

Previous studies have confirmed the sex difference of gray matter asymmetry in typically developing controls and the abnormal gray matter asymmetry in autism spectrum disorders. However, whether and how sex differences of gray matter asymmetry exist in autism spectrum disorders remains studied. This paper analyzes the above issues and explores correlations between gray matter asymmetry and autistic symptoms. Data from 72 children (36 males and 36 females) with autism spectrum disorders and 72 typically developingcontrols (36 males and 36 females) between 8 and 14 years were included and obtained from the autism brain imaging data exchange repository (autism brain imaging data exchange I and autism brain imaging data exchange II). The voxel-based morphometry approach was used to assess gray matter asymmetry in T1-weighted brain data, and gray matter asymmetry was quantified as asymmetry index. A 2 × 2 analysis of covariance was used to identify asymmetry index differences among the four groups. Pearson correlation analysis was performed for asymmetry index values extracted from the clusters with significant differences between the four groups and autistic symptoms (social impairments, communication difficulties, and restricted, repetitive behaviors) measured by the revised autism diagnostic interview scale. Results showed that specific brain regions showed significant main effects for diagnosis in which autism spectrum disorders patients had more leftward asymmetry than typically developing-controls for the parahippocampal gyrus and the postcentral gyrus; specific brain regions showed significant main effects for sex in which females showed more rightward asymmetry for the middle temporal gyrus, inferior frontal gyrus, angular gyrus, and postcentral gyrus and minor rightward asymmetry for the superior frontal gyrus than males; significant diagnosis × sex interaction effects were identified in the angular gyrus and middle occipital gyrus. Pearson correlation analysis showed that males with autism spectrum disorders had a positive association between the asymmetry index value in the middle occipital gyrus and more significant verbal impairment measured by the revised autism diagnostic interview (r = 0.387, p = 0.026). The asymmetry index value in the parahippocampal gyrus was positively associated with more severe social impairment in females with autism spectrum disorders (r = 0.422, p = 0.020). We identified that the sex difference of gray matter asymmetry in children with autism spectrum disorders is qualitative rather than quantitative, which is relatively novel. Our findings provide the theoretical basis for conducting separate studies and using sexspecific diagnostic methods and treatments for males and females children with autism spectrum disorders.

Previous studies have confirmed the sex difference of gray matter asymmetry in typically developing controls and the abnormal gray matter asymmetry in autism spectrum disorders. However, whether and how sex differences of gray matter asymmetry exist in autism spectrum disorders remains studied. This paper analyzes the above issues and explores correlations between gray matter asymmetry and autistic symptoms. Data from 72 children (36 males and 36 females) with autism spectrum disorders and 72 typically developingcontrols (36 males and 36 females) between 8 and 14 years were included and obtained from the autism brain imaging data exchange repository (autism brain imaging data exchange I and autism brain imaging data exchange II). The voxel-based morphometry approach was used to assess gray matter asymmetry in T1-weighted brain data, and gray matter asymmetry was quantified as asymmetry index. A 2 × 2 analysis of covariance was used to identify asymmetry index differences among the four groups. Pearson correlation analysis was performed for asymmetry index values extracted from the clusters with significant differences between the four groups and autistic symptoms (social impairments, communication difficulties, and restricted, repetitive behaviors) measured by the revised autism diagnostic interview scale. Results showed that specific brain regions showed significant main effects for diagnosis in which autism spectrum disorders patients had more leftward asymmetry than typically developing-controls for the parahippocampal gyrus and the postcentral gyrus; specific brain regions showed significant main effects for sex in which females showed more rightward asymmetry for the middle temporal gyrus, inferior frontal gyrus, angular gyrus, and postcentral gyrus and minor rightward asymmetry for the superior frontal gyrus than males; significant diagnosis × sex interaction effects were identified in the angular gyrus and middle occipital gyrus. Pearson correlation analysis showed that males with autism spectrum disorders had a positive association between the asymmetry index value in the middle occipital gyrus and more significant verbal impairment measured by the revised autism diagnostic interview (r = 0.387, p = 0.026). The asymmetry index value in the parahippocampal gyrus was positively associated with more severe social impairment in females with autism spectrum disorders (r = 0.422, p = 0.020). We identified that the sex difference of gray matter asymmetry in children with autism spectrum disorders is qualitative rather than quantitative, which is relatively novel. Our findings provide the theoretical basis for conducting separate studies and using sex-

Introduction
Autism spectrum disorders (ASD) are a set of pervasive neurodevelopmental conditions characterized by social impairments, communication difficulties, and restricted, repetitive behaviors [1]. It is one of the most common developmental disorders, affecting about 1 in 160 people worldwide [2]. Despite its prevalence, studies have not been able to determine its exact etiology [3]. The imbalanced male-tofemale prevalence ratio is one of the most apparent characteristics of ASD [2], and the sex bias for males is even higher in neuroimaging studies [4,5]. There are also symptom differences between males and females with ASD [6]; males have lower levels of social impairments and externalizing problems [7] but more significant restricted behaviors compared to females with ASD [8]. Thus, sex is increasingly being recognized as the source of neuropathophysiological heterogeneity in ASD.
One of the most striking features of the human brain is its structural or functional lateralization or asymmetry. In most typically developing-controls (TDs), the left hemisphere specializes in motor control and language, while the right hemisphere is responsible for visual-spatial attention [9]. And in TDs, studies have found that numerous brain regions showed asymmetry differences in structure (especially gray matter) or function between males and females [10,11]. Alterations of lateralization or asymmetry have been found in various mental and neurocognitive disorders, including schizophrenia and dyslexia [12,13]. And atypical brain asymmetry has long been hypothesized for ASD, a supposition that is supported by multiple neuroimaging studies [14][15][16]. Several structure-based studies reported hemispheric asymmetry abnormalities in ASD, but discoveries of potential structural alterations for ASD have been inconsistent. For example, Groen et al. [17] reported increased hippocampal volume in ASD. On the other hand, Eilam-Stock et al. [18] found significantly decreased hippocampal volume in subjects with ASD compared with TDs. Therefore, further structure-based researches on the hemispheric brain asymmetry of ASD are necessary.
Due to the relatively high prevalence of ASD in males, almost all existing studies on the brain structure of ASD have focused on samples that are predominantly or exclusively male. However, there is growing evidence that males and females with ASD have different brain structural alteration patterns. For example, [19] reported that females might need higher detrimental loads, including brain structural alterations, before developing clinically relevant levels of autistic characteristics. In addition, ASD-associated neural anatomy in females may differ from that in males; Ecker et al. [20] detected that temporal lobe cortical thickness was decreased and increased in females and males with ASD, respectively. Unfortunately, potential sex differences in brain structure have often been ignored in ASD. This may be one of the main reasons for the inconsistent results of brain structure research in ASD. Thus, the inclusion of females with ASD is critical to help improve our understanding of brain structural alterations in ASD.
Based on the findings of the sex difference of gray matter (GM) asymmetry in TDs and the abnormality of GM asymmetry in ASD, quantitative or qualitative sex differences of the GM asymmetry alterations in ASD are probably to be expected [11,21]. However, if the effects on GM asymmetry in males and females with ASD are present in only one sex or the opposite form, previous research on GM asymmetry of ASD that did not consider the sex factor may be unreliable, as the above effects may be offset in these researches. This significantly reduced the possibility of detecting actual effects. Therefore, a sex-specific analysis is urgently needed to determine whether and how the sex differences of brain structural asymmetry might be altered in ASD to better understand the neuroanatomy of the situation.
Recently, Kurth et al. [22] established a fully automated voxel-based morphometry (VBM) method to specifically analyze brain GM asymmetry. A key advantage over the standard VBM workflow is establishing the voxel-wise hemispheric correspondence, which is ensured by using Diffeomorphic Anatomical Registration Lie (DARTEL) algebra to convert spatial normalization into a symmetric space. Several studies have used this approach to investigate the effects of sex, meditation, alcohol dependence, and schizophrenia on GM asymmetry [22][23][24].
Here, we used the advanced VBM method [25] to conduct a sex-specific study of GM asymmetry in ASD. GM asymmetry was quantified as asymmetry index (AI). We also explored correlations between AI values extracted from the clusters with significant differences between the four groups and autistic symptoms (social impairments, communication difficulties, and restricted, repetitive behaviors) measured by the revised autism diagnostic interview (ADI-R) scale.

Participants
The participants included 72 children (36 males and 36 females) with ASD and 72 TDs (36 males and 36 females) between 8 and 14.
We used Gpower (https://www.psychologie.hhu.de/arbeitsgruppen/allg emeine-psychologie-und-arbeitspsychologie/gpower.html) software to calculate the sample size. The power and the required effect size in this paper are 0.85 and 0.24, respectively. All participants were obtained from the autism brain imaging data exchange (ABIDE) repository (ABIDE I and ABIDE II). The dataset consists of 1060 individuals with ASD and 1166 TDs, aged 6-65, collected from 36 international sites. Data were collected from five sites (KKI, NYU, OHSU, YALE, GU) that included five or more females with ASD who met the criteria. All subjects with ASD had a clinician's Diagnostic and Statistical Manual of Mental Disorders (DSM)-IV-TR or DSM-V diagnosis. The majority of participants were assessed using Autism Diagnostic Observation Schedule (ADOS) modules 3 or 4, ADI-R, or both. The ABIDE database also provided information on medications and comorbidities. The female ASD group was formed first because most sites included a few females with ASD. To avoid biases due to the different sites, random subjects were selected to form age-, sex-and site-matched groups. All four groups were formed with the following criteria: (1) age range: 8-14 years, (2) full-scale intelligence quotient (FIQ) ≥70, and (3) right-handedness. Those with poor image quality were excluded. The Wechsler Intelligence Scale measured the FIQ for Children (KKI, GU, OHSU), the Differential Ability Scales (YALE), or the four subtests of the Wechsler Abbreviated Scale of Intelligence (NYU). Handedness was measured by the 22-item Edinburgh Handedness Inventory (KKI, NYU, OHSU) or parent report (YALE, GU). All data were provided by the site and published on the ABIDE website.

MRI data acquisition
All included data were acquired from 3.0-Tesla scanners, but the scanner types and image acquisition parameters varied across sites. The scan parameters and acquisition protocols are provided at http://fcon_1000.projects.nitrc.org/ind i/abide/.

Image analysis 2.3.1 MRI data processing
The high-resolution 3D T1-weighted images were preprocessed using the VBM8 toolbox (http://dbm.neuro.un i-jena.de/vbm8) package Statistical Parametric Mapping (SPM8, http://www.fil.ion.ucl.ac.uk/spm) software. Before preprocessing, visual checks were performed for all 3D T1weighted images. Then, all preprocessing steps were performed following the protocol by Kurth et al. [22]: (1) The 3D T1-weighted images were segmented into GM, white matter (WM), and cerebrospinal fluid; (2) the segmented images were flipped at the midline; (3) the original and flipped GM and WM versions of each participant were used to create a symmetric DARTEL anatomical template; (4) the original and flipped GM and WM tissue segments of all participants were warped into the created symmetric template to align the brain of each participant in a shared space; and (5) to limit further data analysis to the right hemisphere, we used the symmetric template space to create a right-hemispheric mask in MRIcron (http://people.cas.sc.edu/rorden/mricron/index.h tml). First, in the sagittal window of MRIcron, we selected one plane to the right in the midline line of GM. This plane was marked as a volume of interest using the drawing tool. Then, masking was repeated in the next plane until the right hemisphere was finished. Finally, the file was saved in NIfTI format.

Global GM asymmetry analysis
The AI value of GM volume was calculated with an automated procedure. This procedure used the 'calculate' script provided by [22]. First, type 'calculate' in MATLAB's command window to run the script and select 'Step 6', then the script asks for the original warped GM images and the hemispheric mask generated in the above steps [22]. This step generates the masked AI images of each participant. The formula in this script is [22]: where i1 is the original warped image, i2 is the flipped warped image, and i3 is the right-hemispheric mask image.
All AI maps were spatially smoothed with an 8-mm Gaussian smoothing kernel. All results were limited in the right hemisphere, and a positive AI indicated rightward asymmetry (larger right-hemispheric GM volume in a given cluster). Conversely, a negative AI indicated leftward asymmetry (larger left-hemispheric GM volume in a given cluster) [22,26].

Statistical analyses
A 2 × 2 analysis of covariance (ANCOVA) was used to assess the main effect of diagnosis, sex, and their interactions (Age, FIQ, and the site were covariates) for all AI maps SPM8. The Gaussian random field (GRF) theory correction was applied (voxel-level p < 0.01, cluster-level p < 0.05) to determine the significant regions [25].
The AI values were extracted using Data Preprocessing Assistant for rs-fMRI (DPARSF) software from the clusters with significant differences in the 2 × 2 ANCOVA. Then the simple-effect analysis was implemented with SPSS 20.0 software for extracted AI values of the interaction effect. Finally, Pearson correlation analysis was performed for all extracted AI values and ASD clinical symptoms (social impairments, communication difficulties, and restricted, repetitive behaviors) measured by ADI-R scale, p < 0.05 (uncorrected) was used to determine the statistically significant correlation.

Demographic and clinical details
Descriptive statistics are presented in Table 1. The four groups did not significantly differ in age (p > 0.05), but the ASD had lower FIQ than TDs (diagnosis effect) (F = 25.47, p < 0.001). There were no significant differences between females and males with ASD concerning medications, comorbidities, and scores of ADI-R and ADOS modules 3 or 4 (all p > 0.05).

The result of 2 × 2 ANCOVA for AI
All brain regions in results were reported using the Automated Anatomical Labeling template.

Diagnosis effect
Specific brain regions showed significant main effects for diagnosis. The ASD patients had more leftward asymmetry than TDs for the parahippocampal gyrus and postcentral gyrus (GRF correction, voxel-level p < 0.01, cluster-level p < 0.05; Fig. 1).

Sex effect
Relative to males, females showed significantly more rightward asymmetry in the middle temporal gyrus, inferior frontal gyrus, angular gyrus, and postcentral gyrus and minor rightward asymmetry in the superior frontal gyrus (GRF correction, voxel-level p < 0.01, cluster-level p < 0.05; Fig. 2).

Diagnosis-by-sex interaction
Significant diagnosis × sex interaction effects were identified in the angular gyrus and middle occipital gyrus (GRF correction, voxel-level p < 0.01, cluster-level p < 0.05; Fig. 3). Results of the simple-effect analysis were as follows: Females with ASD showed more, and males with ASD showed less rightward asymmetry relative to their respective TDs in the angular gyrus (all p < 0.05) and the middle occipital gyrus (p < 0.05 and p > 0.05 respectively). The detailed results of the simple-effect analysis are shown in Tables 2 and 3. All extracted AI values are listed in Table 4.

Brain-behavior relationships
We observed significant correlations between extracted AI values and clinical severity in the ASD group. Specifically, males with ASD showed a positive association between the AI value in the middle occipital gyrus and more significant verbal impairment measured by ADI-R (r = 0.387, p = 0.026; Fig. 4). In addition, the AI value in the parahippocampal gyrus was positively associated with more severe social impairment in females with ASD (r = 0.422, p = 0.020; Fig. 5 and Table 5).

Discussion
We mapped differences in GM asymmetry between ASD and TDs and between the male and female children with ASD. Previous studies on cortical asymmetry in ASD found variable differences or no differences [27]. Although Postema et al. [27] reported that total hemispheric results Table 1. Group characteristics for the included male and female participants, separately for the ASD and TD groups.

Males with ASD Females with ASD Males TDs Females TDs Diagnosis effect Sex effect Diagnosis × Sex n mean ± SD (rate) n mean ± SD (rate) n mean ± SD n mean ± SD F-value p-value F (t)-value p-value F-value p-value
Age (  showed generalized less leftward asymmetry of cortical thickness in ASD, the direction of GM asymmetrical alterations in specific brain regions remains to be established. Using the advanced VBM approach, we found more leftward asymmetry of GM volumes in the parahippocampal gyrus and postcentral gyrus in ASD. The parahippocampal gyrus is essential for ASD. Abnormalities in this brain region have been reported in previous studies on ASD. Mei et al. [28] reported decreased density in the parahippocampal gyrus in ASD. Yang et al. [29] found that the GM volume in the parahippocampal gyrus in adults with ASD was significantly increased relative to TDs. In addition, the AI value in the parahippocampal gyrus was positively associated with more severe social impairment in females with ASD (r = 0.422, p = 0.020). This result is consistent with previous studies that linked structural alterations in the parahippocampal gyrus to social deficits in ASD [30], indicating that structural alterations of the parahippocampal gyrus indeed serve a junctional role in ASD etiology.
Concerning the postcentral gyrus, volume abnormalities reported in this brain region in ASD have been inconsistent [31]. Mizuno et al. [31] reported decreased left postcentral gyrus volume in ASD , while Brieber et al. [32] described increased volume in this brain region in patients with ASD comparing with TDs. This may be due to the small sample sizes, the different proportions of females with ASD, and the  different image processing techniques among different studies. We used an advanced VBM method specifically designed to assess brain structural asymmetry, had a relatively large sample size, and considered the sex factor. Therefore, our findings are more reliable and can provide more substantial evidence to clarify the inconsistent findings reported earlier.
The postcentral gyrus participates in somatosensory functions. Abnormal sensory processing can lead to social impairments, communication difficulties, and restricted, repetitive behaviors, which are core characteristics of ASD [33]. There- fore, the abnormality of the postcentral gyrus may cause the core symptoms of ASD.
The sex difference of GM asymmetry was found in the angular gyrus in ASD. Males and females with ASD showed more minor and more rightward asymmetry in the angular gyrus, respectively, compared with their sex-matched TDs. A similar reduction in the angular gyrus was previously detected in [34] with the mixed-sex ASD sample (mainly males with ASD). Thus, less rightward asymmetry in the angular  gyrus seems to be unique to males with ASD. Interestingly, this brain area is involved in the theory of mind, which is essential to help the human brain reason about others and effectively communicate and navigate in the social world [35]. There are sex differences for empathy in ASD [36], and the theory of mind is necessary for empathy. Therefore, sex differences of GM asymmetry in the angular gyrus may be why males are less empathetic than females with ASD. The sex difference of GM asymmetry was also observed in the middle occipital gyrus in subjects with ASD. Like the angular gyrus, females with ASD showed more, and males with ASD showed less rightward asymmetry relative to their respective TDs in the middle occipital gyrus. This is the first work to find that males and females with ASD exhibited GM asymmetry alterations in the opposite direction relative to their sex-matched TDs. Moreover, males and females with ASD showed different correlativity between ASD symptom severity and AI value in several brain regions, including the middle occipital gyrus. This again suggests that the relationship between ASD and GM asymmetry partly varies by sex. The middle occipital gyrus mainly involves constructing visual and motorial perception, and aberrant perceptual processing might be associated with social dysfunction in ASD [37]. Considering the novelty of these results, they should be interpreted carefully, and they need to be replicated in future studies.
The results do not support the "Female Protective Effect" theory, which states that females need higher detrimental loads, including brain structural alterations, before developing clinically relevant levels of autistic characteristics [19]. Instead, we detected opposite sex-associated effects in ASD for both two brain regions. Similarly, [38] found that neither behavioral nor neural data supported the idea that females with ASD had more severe abnormalities. Thus, our research provides evidence that sex differences of GM asymmetry in ASD may be qualitative rather than quantitative. Based on these findings, we suggest that researchers should not conduct cross-sex neuroimaging studies for ASD because they may raise the risk of offsetting present results in only one sex or opposite directions in males and females.
There are limitations in the current research. First, it was based on multi-site datasets, which may introduce additional sources of variance associated with participant characteristics or scan parameters. Secondly, our sample size was still insufficient; further studies with larger sample sizes are needed to confirm sex-specific brain structure alterations in subjects with ASD. Finally, participants' clinical information was also limited, with ADOS and ADI-R scores available for only participants.  Using an advanced VBM approach, we identified diagnosis × sex effects of GM asymmetry in ASD. These effect patterns were qualitative rather than quantitative by showing that males and females with ASD exhibited GM asymmetry alterations in the opposite direction, which is relatively novel. Our findings highlight the necessity of separate studies for males and females with ASD and the potentials for sexspecific diagnostic methods and treatments.

Author contributions
CCL, MMN, PYF and HBX conceived and designed the experiments; CCL and MMN analyzed the data; CCL and HBX wrote the paper.

Ethics approval and consent to participate
The written informed consent was obtained from the site where they conducted the tests. The research was conducted with the Ethics Committee of the University of Zhongnan Hospital of Wuhan University (Ethical approval number: 2019072).